כיצד לבצע ניתוח התרמודינמיקה של מנוע אוטו-מחזור באמצעות Matlab

הבנת דגם Air-Standard אוטו-מחזור

מחזור אוטו הוא אידיאליזציה תרמודינמית של מנוע הפחתת הניצוץ 4 שכבות.זה מחליף את הצריכה המורכבת, דחיסה, דחיסה, תהליכים ממצה עם ארבעה שלבים בלתי הפיכים המבודדים את מנגנוני ההמרה האנרגיה החיוניים. במודל סטנדרטי אוויר, נוזל העבודה הוא אוויר (טיפול כמו גז אידיאלי עם חום קבוע), כפילציה הוא קרוב כמו תוספת מיידית חום, ורמת ביצועים קבועים, מאפשר פשטה, לחץ קבוע, הוא קבוע, הוא קבוע של פעילות גופנית, הוא קבוע, הוא קבוע, הוא קבוע, הוא קבוע, הוא קבוע, הוא קבוע, הוא קבוע, החליקציה, הוא תהליך ניקוי, הוא קבוע של פעילות גופנית, הוא קבוע, הוא קבוע של לחץ, הוא קבוע, הוא קבוע של לחץ, הוא קבוע, הוא קבוע, החלמה, החלמה, הוא קבוע, הוא קבוע של לחץ, החליקציה, הוא קבוע של פעילות גופנית, הוא קבוע, הוא קבוע, הוא קבוע, החליקציה, החליקציה, החליקציה, החלמה, החליקציה, הוא קבוע של פעילות גופנית קבועה של פעילות גופנית, הוא קבוע, החליקציה, הוא קבוע של פעילות גופנית, הוא קבוע, הוא קבוע של פעילות גופנית, הוא קבוע, הוא קבוע, הוא קבוע, החליקציה, הוא קבוע,

ארבעת התהליכים שמרכיבים את המעגל הם:

• (ב) דחיסה של Isentropic (1 מ)1.10.10.10.10.10: הפיסטון מחסחסם את המטען ללא העברה חום, גידול טמפרטורה ולחץ.
• (FLT:0) תוספת חום של מים (2 מ) 30303303FLT: שחרור מיידי של אנרגיה (ייצוג של בעירה) במרכז המתים העליון מעלה את הטמפרטורה והלחץ לערכי השיא שלהם.
• (FLT:0) התרחבות איסנטרופית (3=1-4FLT:1): גזים בלחץ גבוה דוחפים את הפיסטון למטה, ומספקים עבודה נטו ל- crankshaft.
• (ב) דחיית חום (למעלה מ-1): תכולת הצלינדר מאבדת חום חזרה למצב הראשוני, סגירת הלולאה.

מכיוון שהמודל משתמש יחס חום קבוע ומתעלמו מחילופי גז, הוא מספק מדד מהיר שניתן לחדד באופן שיטתי.מיישם את המחזור הזה ב MATLAB נותן לך כלי גמיש ואינטראקטיבי לבדיקת יחס דחיסה, תוספת חום ונכסים נוזליים המשפיעים על הביצועים.

בסביבה של MATLAB

החל על ידי יצירת קובץ חדש או פונקציה.מרחב עבודה נקי עושה פיזור ושינוי עתידי פשוט. השתמש ראש תיאורי המעדכנת את המטרה, פרמטרים קלט, ופלטים.מערכים קדם-כולים למשתנים ממשלתיים כאשר גורף מעל פרמטרים כדי למנוע קידוד דינמי מעל פני מעלה, אימוץ שמות משתנים משמעותיים כגון LT:0 ולא רק LTF:1 (למרות סטנדרטיים) הוא יכול להיות חסימות סטנדרטיות (R)

% otto_cycle_analysis.m
% Air‑standard Otto cycle analysis
% Inputs:
% r – compression ratio (dimensionless)
% T1, P1 – intake state (K, Pa)
% gamma – specific heat ratio (cp/cv)
% Q_in – specific heat addition (J/kg)
% Outputs:
% state_table – array of T, P, v at four corners
% W_net, eta – specific net work (J/kg) and thermal efficiency

(ב) שימוש ב-MTLAB's Built-inFLT:4 או ב-FLT:5 ⁇ כדי לאחסן את ארבע המדינות משפר את יכולת הקריאה כאשר אתה מבסס או את תוצאות היצוא מאוחר יותר, למשל, ליצור מבנה מערך של מערך 6FLT 6 עם שדות עם שדות (FLT 7,FLT:8, FLT:9, לעומת 1 עד 4, זה הופך את התכונות האינטואיטיביות של המדינה:

Defining the Critical Engine Parameters

הדיוק של הסימולציה תלוי בנתונים קלט.מינימום, ציין את יחס הדחיסה, לחץ הצריכה וטמפרטורה, יחס החום הספציפי, ואת האנרגיה שנוספה לכל מסה יחידה.

• יחס דיכוי (FLT:0) יחס היחס לצמצום 1: 8-12 עבור מנועי הניצוץ המודרניים; יחס גבוה יותר משפר את יעילות אך מגביר את הסיכון להורדת לחץ, במנועי מחקר, יחס גבוה ככל ש-14 נחקרים לעיתים עם דלקים מתקדמים.
• (FLT:0) צריכת הטמפרטורה T1FLT:1: 300 K (עמ') או 350 K אם יש צורך בצריכת אוויר מגזים כפולים וגורשים.
• (ב) לחץ P1FLT: 101325 פא (1 בשעה) באופן טבעי מחלחל; עד 200 kPa או יותר עבור מנועים מכווצים.
• יחס חום ⁇ FLT:1 [: 1.4 עבור אוויר יבש; 1.3-1.35 עבור תערובת אוויר דלק ומוצרי הבעירה בשל נוכחות של מולקולות 3 atomic כמו CO2 ו- H2O.
• (FLT:0) תוספת חום של Q inFearLT:1; בערך 2.5-3 MJ / ק"ג עבור בנזין בהתבסס על ערך חימום נמוך יותר ו יחס דלק אוויר סטוכמטרי. עבור בנזין טיפוסי, הערך הנמוך יותר הוא בערך 44 מ"ג / ק"ג, אך מעורבב עם אוויר (התוצאה 14.7:1 דלק אוויר), התערובת הופכת בערך 2.9J / ק"ג.

Define אלה כמו משתנים MATLAB בחלק העליון של התסריט כך שניתן לשנות בקלות למחקרי רגישות.

r = 9.0; % compression ratio
T1 = 300; % K
P1 = 101325; % Pa
gamma = 1.35; % approximate for fuel‑air mixture
Q_in = 2500e3; % J/kg
R_air = 287; % J/(kg·K)

אתה יכול גם לכלול את מסה של אוויר (28.97 גרם / מול) אם אתה מאוחר יותר לעבור חישובים המבוססים על שומנים.שמירה על הגז קבוע באופן מפורש נמנע שגיאות יחידה.

ארבע המדינות קורנות

עם הפרמטרים שנקבעו, יחסי האוויר סטנדרטיים נותנים את המדינה נקודות ישירות.הניתוח הוא לכל מסה יחידה, כך שאנו עובדים עם נפח מסוים.קביעת נפח ההנקה לאחדות (v2= 1) מפשט את הסקאלה היחסית; ואז v1= r כי יחס הדחיסה v1/v2= r. זה נורמליזציה נוח כי כל הנפח מובע ביחס לנפח הנקה.

מדינת 1 (כניסה, מרכז מתים)

v1 = r; % relative to clearance volume
% T1 and P1 are inputs

המדינה 2 (Isentropic Compression to Top Dead Center)

עבור תהליך אנטרופי עם קבוע ⁇ :

v2 = 1;
T2 = T1 * r^(gamma-1);
P2 = P1 * r^gamma;

יחסים אלה נובעים מחוק הגז האידיאלי והיחס הפוליטרופי:0 (p v ⁇ = קבועFLT:1 ). הטמפרטורה עולה על דחיסה היא משמעותית: עבור r=9 ו ⁇ 1.3=5, T2 ⁇ 300 × 9(0.35) 690 K.

המדינה 3 (המשך לתוספת חום)

ראשית, לשים את החום הספציפי בנפח קבוע:

cv = R_air / (gamma-1);
T3 = T2 + Q_in / cv;
P3 = P2 * (T3 / T2);
v3 = v2; % no volume change

הטמפרטורה במצב 3 יכולה לעלות על 2500 K עבור דלק טיפוסית צריכת הדלק P3 לעתים קרובות מגיע 60-80 בר, אשר קובע דרישות מבניות עבור בלוק המנוע וראש.

המדינה 4 (התרחבות איסנטרופית לכרך המקורי)

יחס ההתרחבות מ- TDC ל- BDC הוא r, אך מבחינת נפח מסוים v3= v2= 1 ו- v4= v1= r, כך יחס נפח הוא 1/r.

v4 = v1;
T4 = T3 * (1/r)^(gamma-1);
P4 = P3 * (1/r)^gamma;

ארבעת הנקודות הללו מגדירות במלואן את המחזור במונחים של טמפרטורה, לחץ ונפח ספציפי.חנות אותם במערךים עבור מזחלות ו חישוב מאוחר יותר. אימות שימושי הוא לבדוק כי P4 הוא מעט מעל P1; ההבדל מניע את הפיצוץ הממצה.

חישוב העבודה נטו, יעילות תרמית, ומשמעות לחץ יעיל

העבודה הספציפית הנקייה היא ההבדל בין חום מוסף חום נדחה.עבור המודל הסטנדרטי, החום נדחה בתהליך 4.1 הוא:

Q_out = cv * (T4 - T1);
W_net = Q_in - Q_out;

יעילות תרמית היא:

eta = W_net / Q_in; % from energy balance
eta_direct = 1 - 1/r^(gamma-1); % direct formula
% Both should match within numerical round-off

(FLT:0) Mean יעיל לחץ (MEP)FearLT:1 הוא מדד ביצועים המשווה את העבודה נטו לכל מחזור לנפח העקירה.עבור מחזור סטנדרטי האוויר, נפח העקירה הספציפי הוא v1 - v2= r - 1. לכן: 1.

mep = W_net / (v1 - v2); % Pa

MEP הוא עצמאי בגודל המנוע והוא אינדיקטור נוח של איכות עיצוב המנוע.ערכים אופייניים עבור מנועי בנזין באופן טבעי נשפים הם 8-12 בר (0.8-1.2 MPa) מחוץ ל- MEP לצד יעילות לספק תמונה מלאה יותר.You עשוי גם להמיר לבריות: FLT:19 .

כולל טמפרטורות ספציפיות ליעילות משופרת

(ההנחה הקבועה של ⁇ ) יכולה להגביר את היעילות מכיוון שלגזים ממצה אמיתיים יש יחסי חום נמוכים יותר בטמפרטורות גבוהות.לחשוב על כך, להחליף את ה ⁇ הקבועה עם מתאם תלות בטמפרטורה. שיטה חזקה היא להשתמש בתקני פולינומיים פולינומיים של נאס"א עבור מוצרים סטנדרטיים:0cLTFLT:2,2FLT3 של מינים עיקריים (O2, CO2, C2, C2N2, C2, C2, C2N)

יישום ⁇ משתנה דורש פתרון הרהרטיבי כי ה- isentropic exponent עצמו תלוי הטמפרטורה הלא ידועה. השתמש MATLAB'sFLT:20 או לולאה פשוטה.

% Use an anonymous function to find T2 such that
% integral of cp dT from T1 to T2 = R * ln(v1/v2)
% This replaces the constant‑gamma formula.
% For air, cp(T) can be approximated by:
% cp = a0 + a1*T + a2*T^2 + a3*T^3 + a4*T^4
% Using NASA 7-coefficient polynomials.
fun = @(T2) integral(@(T) polyval(cp_coeffs, T), T1, T2) - R_air * log(r);
T2 = fzero(fun, T1 * r^0.4);

ה-FLT:0 NASA פולינומאלי משאב FIRLT:1 מציע מזהמים להורדה עבור מינים אוויריים ובעירה. Modify את התסריט שלך כדי לטעון את האפקטיביות הללו ולקרוא להם במהלך חישובים המדינה.You תבחן כי חום ספציפי משתנה להפחית את היעילות הצפויה על ידי 2-5 נקודות ביחס דחיסות טיפוסיות, אשר תואמות הדוק עם מנוע נתונים אמיתי.

עבור היערכות מהירה ללא שילוב פולינומי מלא, השתמש יחס חום מסוים ממוצע משתנה ליניארי עם טמפרטורה: ⁇ (T) = 1.4 -0.00005*(T - 300) זה פשוט לוכד את המגמה ויכול להיות מיושם בלולאה תוך זמן עד ההתכנסות.

בניית ניתוח פונקציונלי מודולרי -Driven Analysis

במקום תסריט מונוליטי, מטביע את לוגיקה מחזורית בתפקוד MATLAB.זה הופך אותו לזמין עבור פרמטר גורפים ולהפחית את שכפול הקוד. a פונקציה השלד מוצג מתחת:

function [states, performance] = ottoCycle(r, T1, P1, Q_in, gamma, R_air)
 % states: structure with T, P, v for states 1-4
 % performance: structure with W_net, eta, mep
 % (calculation code from previous section goes here)
end

עבור חום ספציפי משתנה, להוסיף דגל נוזלי:

function [states, performance] = ottoCycle(r, T1, P1, Q_in, fluid_type)
 % fluid_type: 'constant_gamma' or 'variable_cp'
 % (logic branches based on flag)
end

עם הפונקציה במקום, תסריט נהג יכול לטאטא באמצעות יחסי דחיסה:

r_values = 6:1:14;
for i = 1:length(r_values)
 [~, perf] = ottoCycle(r_values(i), 300, 101325, 2500e3, 'constant_gamma');
 eta(i) = perf.eta;
 mep(i) = perf.mep;
end
plot(r_values, eta*100, 'b-', r_values, mep/1e5, 'r--');

גישה מודולרית זו מציגה מיד את ההחזרים ההפחתה של יחסי הדחיסה הגדלים מעבר ל-12.המזימה מדגישה את משחק הסחר המרכזי: יעילות ממשיכה לעלות אך בקצב ההפחתה, בעוד ש-MeP לעתים קרובות עולה סביב r=10 עבור קלטות חום טיפוסיות.

יצירת אבחון P-V ו- T-S Diagrams

Visualization מסייע לאמת את המודל ולתקשר תוצאות.תפקודי העלילה של MATLAB מייצרים דיאגרמות באיכות פרסום.

לחץ-Volume (P-V) Diagram

דחיסה והתרחבות מתעמימות על ההרחבה (FLT:0)P v ⁇ =קבוע =קבוע = ReFLT:1.

n = 50;
v_comp = linspace(v2, v1, n); % reverse direction for plotting
P_comp = P2 * (v_comp ./ v2).^(-gamma);
v_exp = linspace(v3, v4, n);
P_exp = P3 * (v_exp ./ v3).^(-gamma);
figure;
hold on;
plot(v_comp, P_comp, 'b', 'LineWidth', 1.5);
plot(v_exp, P_exp, 'b', 'LineWidth', 1.5);
plot([v2 v3], [P2 P3], 'r--', 'LineWidth', 1.2);
plot([v4 v1], [P4 P1], 'r--', 'LineWidth', 1.2);
xlabel('Specific volume (relative)'); ylabel('Pressure (Pa)');
title('Air‑Standard Otto Cycle P–V Diagram');
legend('Compression', 'Expansion', 'Heat addition', 'Heat rejection');
grid on;

האזור הסגור מייצג את העבודה נטו.You יכול לשלב את הלולאה באמצעות ההרחבה (FLT):26 כדי לאמת את W net הנחשק.לדוגמה, לשלב את העקומה העליונה (דחיסה + הרחבה) ולהכפיל את העקומה התחתונה (עמוד דחיית חום) לתיעוד FLT:0MATLAB על אינטגרציה מספרית:1 לפרטים.

טמפרטורה-Entropy (T-S) Diagram

עבור גז אידיאלי עם חום ספציפי קבוע, שינוי אנטרופי ביחס למדינה 1 הוא:

% s = cp * log(T / T1) - R * log(P / P1)
cp = gamma * R_air / (gamma - 1);
s1 = 0;
s2 = cp * log(T2/T1) - R_air * log(P2/P1);
s3 = cp * log(T3/T1) - R_air * log(P3/P1);
s4 = cp * log(T4/T1) - R_air * log(P4/P1);

העלילה של ארבעת התהליכים עם סגנונות קו מתאימים. עבור ⁇ משתנה, להשתמש שילוב מספרי של dQ / T כדי compute entropy. ⁇ T-S מגלה את האופי האנטרופי של תהליכים 1-2 ו 3-4 (קווים מלוכלכים) ואת קבוע חימום כרכים קבוע / כרית (נתיבים מעוקלים).

כדי לשפר את העלילה של T-S, קווי לחץ מתמיד (isobars) מיחס הגז האידיאלי.זה עוזר להקשר בין המחזור במסגרת תרמודינמית רחבה יותר.

ביצוע Sרגישות מחקרים

ברגע שמודל הבסיס עובד, באופן שיטתי לשנות את הפרמטרים קלט כדי להבין את השפעתם.מחקרים משותפים כוללים:

• (FLT:0) יחס דיכוי גורף 1:1: יעילות פשטות ו- MEP לעומת r. Efficiency עולה במהירות על שער נמוך ולאחר מכן מישורים; meP לעתים קרובות מגיע לשיא של שער מתון בגלל ההחלפה בין לחץ על עלייה תוספת חום. השתמש subplots כדי להראות שני מדדים על אותו הדמות.
• <חזק>Heat תוספת וריאציות: שינוי Q in מודלים שונים של אנרגיית דלק או יחסי דלק אוויר. הגדלת Q in מעלה לחץ גבוה ועבודה אך עשוי לעלות על גבולות חומריים; כולל מעצורים ש- P3 < 150 בר (15 MPa) כדי לשקף כוח מנוע טיפוסי.
• (FLT:0) יחס חום משמעותי אפקט יחס חום 1FLT: ⁇ התחתונה מפחיתה את היעילות. overlay קבוע סיבולת ואפקטים משתנים כדי להראות את התיקון.זה חשוב במיוחד כאשר מנתחים דלקים חלופיים כמו גז טבעי או מימן.

עבור גורפת דו-מטר, השתמש בלולאות מאונן ותוצאות בחנות במריצה. ואז ליצור מקהלת או מפת צבע כדי לדמיין את ההשפעה המשולבת:

r_vec = 6:1:14;
Q_vec = 2e6:0.2e6:3e6;
for i = 1:length(r_vec)
 for j = 1:length(Q_vec)
 [~, perf] = ottoCycle(r_vec(i), 300, 101325, Q_vec(j), 1.35, 287);
 eta_mat(i,j) = perf.eta;
 end
end
contourf(r_vec, Q_vec/1e6, eta_mat', 20);
xlabel('Compression ratio'); ylabel('Q_in (MJ/kg)');
colorbar; title('Thermal Efficiency (%)');

מפות רגישות כאלה הן בלתי הולמות לניתוח של סחרחורת בעיצוב מנוע.הם מסייעים לזהות אזורים תפעוליים שבהם יעילות גבוהה בעוד הלחץ העליון נשאר מקובל.

אימות התוצאות שלך

לפני השימוש במודל להחלטות הנדסיות, לבדוק נגד מדדים ידועים.נוסחת יעילות סטנדרטית אוויר ⁇ = 1- 1 /r( ⁇ 1) מספק אימות ישיר של המקרה הקבוע של ⁇ , בנוסף, משווה לחץ שיא מחושב P3 נגד נתונים ממפות מנוע או דוגמאות ספרים שפורסמו (למשל, סטון, FLT:0 Introduction to Internal Combustion EnginesFLTs) אשר לעתים קרובות עולה מ-S-Firdexitation (למשל, לדוגמה, לדוגמה, CLT) חומרים מורכבים של מנועים (למשל, כולל חומרים מורכבים של מנועים) ו-S (למשל, LT) LT) סימולציות (D.

טכניקה נוספת של אימות היא לבדוק את ה-FLT:0 מאזן האנרגיה של ההרחבה 1: העבודה נטו בתוספת חום נדחה צריך שווה את חום הקלט בתוך דיוק צף נקודות, גם לאמת כי תהליכי הדחיסה וההתרחבות מספקים את היחס האנטרופי (FLT:2T v( ⁇ 1) = קבועFLT:3 כדי לזהות שגיאות קידוד.

פיתוח המודל לאפקטים ריאליים

ניתן להרחיב את המסגרת הבסיסית כדי לכלול תופעות בעולם האמיתי:

• (ב) [ה]: [ה] [ה]] [ה]] [ה] [ה]] [ה]]][ה]]][ה]]][ה]]]][ה]]] ל[[המאה ה-20]], [ה[[המאה ה-20]]], [[המאה ה-20]], [[המאה ה-20]]]], [[המאה ה-20]]]]]]]], [[המאה ה-20]], [[המאה ה-20]].
• (FLT:0) העברת הפסדי FLT:1: שילוב מתאם מסוג Woschni עבור העברה חום אחידה לקירות גלילי.זה דורש טמפרטורה ולחץ מטושטש.הת הקורלציה מעניקה את הפחתת החום כתפקוד של עווית, מהירות פיסטון, טמפרטורה ולחץ.
• (ב) [15] ,0) , 000 ⁇ גז (EGR)FLT:1: התאמת יחס חום ספציפי ורכב גז כדי לשקף דילולציה.ה-FLT:2MATLABscape Engine Cylinder חוסם FLT 3: מספק מרכיב מוכן עבור מחקרים כאלה, אך ניתן גם ליישם מודל מלוטש פשוט על ידי שילוב של חום ספציפי של גזים ספציפיים.

עבור הרחבה מהירה, להתחיל עם הפונקציה Wiebe.דוגמה, ניתן למצוא את התוספת הפשוטה של תזמון ניצוץ על יעילות ולחץ שיא.

שיפור הניתוח לתוך זרימת עבודה גדולה יותר

בפרויקט הנדסי, ניתוח מחזורי חייב להאכיל סימולציות ברמת הרכב או בדור הדו"ח. השתמש ב- MATLAB'sFLT:30 לייצא שולחנות המדינה ומדבירי ביצועים ל- CSV לעיבוד נוסף ב- Excel.שלב מזימה ל- PDF באמצעות הפונקציה FLT:31.

results = table(r_values', eta', mep', ...
 'VariableNames', {'CompressionRatio','Efficiency','MEP_bar'});
writetable(results, 'otto_sweep.csv');

ניתן גם לקרוא MATLAB מ Python או Simulink כדי ליצור את המודל התרמודדי עם תת-מערכות אחרות.עבור מחקר רכב היברידי, מפת יעילות מחזור אוטו ניתן לייבא כשולחן מבט ב- Simulink.

מסקנה

בניית מודל מחזור סטנדרטי של אוטו MATLAB מספקת פלטפורמה חזקה אך פשוטה לחקור thermoדינמית מנוע ניצוץ. על ידי הגדרת פרמטרים, נקודות מצב מחשוב, ו חישוב מדדי ביצועים כגון עבודה נטו, יעילות תרמית, ומשמעות לחץ יעיל, אתה מקבל תובנה מיידית כיצד בחירות עיצוב משפיעות על יישום.