Air-Strate Otto Çevrim Modeli Anlamak

Otto döngüsü, dört yüzüncü kıvılcım ateşleme motorunun termodinamik idealizasyonudur. Sürekli hacimde, sıkıştırılma, yanma ve egzoz işlemlerinin temel enerji dönüşüm mekanizmalarının izole edilmesi ile yapılır.Bu basitleştirme, doğrudan kapalı-formasyonda, basınçlar, çalışma sıvısı ve ısı verimliliğin hesaplanmasında, sürekli olarak gerekli ısının değiştirilmesi ile ilgili olarak, sürekli olarak gerekli bir adım olarak değişir.

döngüsü oluşturan dört işlem şunlardır:

  • [FONT:0)Isentropic kompresyonu (1+2)).: piston ısı transfer olmadan şarjı sıkıştırır, sıcaklık ve baskı sağlar.
  • [FONT:0)Constantvol-ume ısı ek (2+3)[[[Dönemli: Enerjinin anlık salıverilmesi (enerjinin açıklanması) en üst ölü merkezde sıcaklık ve basınçlarını zirve değerlerine yükseltmektedir.
  • [FONT:0)Isentropic genişleme (3+4)): Yüksek basınçlı gazlar, pistonu aşağı itir, net çalışmayı krankshaft'a taşır.
  • [0]Constantvol-ume ısı reddedilmesi (4+1)[[Dön 1: 1)[Dönemli içerik, ilk duruma geri ısı kaybeder, döngüyü kapat.

Model sürekli belirli bir ısı oranını kullanır ve gaz değişimi göz ardı eder, aynı zamanda verimlilikte uygun bir üst sınır sağlar; gerçek motorlar her zaman MATLAB'de esneklik, ısının nasıl etkilendiğini araştırmak için esnek, etkileşimli bir araç sunar.

Bir Yapılı MATLAB Çevresini Belirle

Yeni bir senaryo veya işlev dosyası yaratarak başlayın. Temiz bir çalışma alanı, silinme ve gelecekteki değişiklikleri basitleştirir. sadece [[Döntücü bir başlık kullanın, amaç, giriş parametreleri ve çıktıları. Pre-allocate dizileri için dinamik bir yeniden yüklemeden kaçınmak için.

% 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

MATLAB'in inşa edilmiş-inurFLT:4) veya [[DörtÜyetim: Bu, daha sonra arsa veya ihracat sonuçlarınızı geliştirirken, bir yapı dizisi oluşturmak için kullanılabilirlik sağlar.

Eleştirel motor parametrelerini tanımlamak

Simülasyonun doğruluğu giriş verilerine bağlıdır. En azından, sıkıştırma oranını, basınç ve sıcaklık oranını, belirli ısı oranını ve enerji bir benzin motoru için birim başına eklenmiştir.

  • [FONT=0]Compression oranı r[DÜDÜDÜDÜDÜDÜDÜDÜDÜDÜDÜDÜSÜDÜDÜDÜDÜDÜDÜDÜDÜDÜDÜDÜDÜDÜDÜDÜDÜDÜDÜSÜSÜSÜDÜSÜDÜ: 0,0)Komünlik oranı, 14'ün yüksek orandaki oranlar bazen gelişmiş yakıtlarla araştırılıyor.
  • [0]Kaptan sıcaklığı T1[[Dönemli) veya 350 K, eğer iskan ve oturma gazlarından hava ısıtma için muhasebe.
  • [[Düzücükler için baskı P1[[Dönetici: 101325 Pa (1 atm) doğal olarak aspirasyonlu olarak; 200 kPa veya daha fazla turboşarj motor için.
  • [FONT:0)Mutfak ısı oranı ⁇ [DÜT:1): 1. Kuru hava için 1.4; CO2 ve H2O gibi tri-atomik moleküllerin varlığı nedeniyle yakıt-hava karışımları ve yanma ürünleri için 1.35.
  • [FONT:0]Specific Heat Ek Q in[DÜT:1): Yaklaşık 2.5-3 MJ /kg, daha düşük ısıtma değeri ve stoichiometrik hava yakıt oranına göre, tipik bir benzin için, daha düşük ısıtma değeri 44 MJ/kg ile karıştırılır, ancak hava ile karıştırılır (yaklaşık 14.7:1 hava yakıt yakıt oranı), karışım enerji yoğunluğu yaklaşık 2.9 MJ / kg karışıma göre.

Bunları senaryonun üst kısmındaki MATLAB değişkenlerini tanımlamak, böylece hassas çalışmalar için kolayca değiştirilebilirler.

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)

Ayrıca hava kütlesini de dahil edebilirsiniz (28.97 g/mol) daha sonra mol bazlı hesaplamalara geçiş yaparsanız gaz sabit açık kaçınmaları birim hatalarına dikkat edin.

Dört Köşe ABD'yi Hesaplamak

Parametreler ayarlandığında, hava standart ilişkileri doğrudan devlet puanlarını verir. Analiz birim başınadır, bu yüzden belirli hacimlerle çalışırız.Bağış hacmini birliğe ayarlarız (v2 = 1) göreceli ölçeklendirmeyi basitleştirir; o zaman v1 = r çünkü kompresyon oranı v1/v2 = r. Bu normalleştirme uygun çünkü tüm hacimler tespit hacmine göre ifade edilir.

Devlet 1 (Intake, Bottom Dead Centre)

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

State 2 (Isentropic Essay to Top Dead Centre)

Sürekli ⁇ ile bir izentropic süreci için:

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

Bu ilişkiler ideal gaz yasasından ve politropik ilişkiden kaynaklanmaktadır:0)P v ⁇ = sabit) kompresyondaki sıcaklık artışı önemlidir: r=9 ve ⁇ =1.35, T2 ⁇ 300 × 9 .(0.35) ⁇ 690 K.

Devlet 3 (Constant-Volume Heat Addition)

İlk olarak sabit hacimdeki özel ısıyı hesaplamak:

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

Devlet 3'teki sıcaklık tipik benzinli yanma için 2500 K'yi geçebilir. Peak basıncı P3 genellikle 60-80 bar'e ulaşır ve bu da motor blok ve kafa için yapısal gereksinimleri belirler.

State 4 (Isentropic Dahili Cilte Genişleme)

TDC'den BDC'ye genişleme oranı bir r, ancak belirli hacimler açısından v3 = v2 = 1 ve v4 = v1 = r, bu yüzden hacim oranı 1/r.

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

Bu dört nokta, sıcaklık, baskı ve belirli hacim açısından döngüyü tamamen tanımlar. Daha sonra arsa ve hesaplama için bunları diziler. yararlı bir doğrulama, P4'ün P1'in biraz üzerinde olduğunu kontrol etmek; farklar egzoz darbesini arttırır.

Net Çalışmayı hesaplamak, Termal Verimliliği ve Etkili Baskı

Net özel çalışma, ısı eklenmiş ve ısı reddedilen farkdır. Hava standart modeli için, işlemde reddedilen ısı 4+1:

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

Termal verimlilik:

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

[FONT:0)Mean etkili baskı (MEP)), boş işi yerinden edilmeye karşı döngüye kıyasla karşılaştıran bir performans ölçümdür. Hava standart döngüsü için, özel yerinden edilme hacmi v1 – v2 = r – 1.

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

MEP motor büyüklüğünden bağımsızdır ve motor tasarımı kalitesinin uygun bir göstergesidir. Doğal olarak aspirasyonlu benzin motorları için tipik değerler 8-12 bar (0.8-1.2 MPa). Çıktı MEP, daha tam bir resim sağlamak için verimlilik sağlar. Ayrıca okuma kabiliyeti için bara dönüştürebilirsiniz:

Geliştirilmiş doğruluk için değişken ısıları genişletin

Sürekli ⁇ varsayımı en yüksek verime sahip olabilir çünkü gerçek egzoz gazları yüksek sıcaklıklarda belirli ısı oranlarına sahiptir.Bu nedenle, sürekli ⁇ 'yi bir sıcaklık bağımlısı ile değiştirin. sağlam bir yöntem NASA polinom katlarını kullanmak içindir.) Bu katta standart olarak ısıtılabilir. ) [O2, N2, CO2, H2O).

Değişkenleri uygulamak, iteratif bir çözüm gerektirir çünkü isentropic exponentin kendisi bilinmeyen sıcaklığa bağlıdır. MATLAB'sULLD::20 veya basit bir döngü. Örneğin, isentropic kompresyonu sırasında, ilişki olur:

% 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);

[FONTA polinom kaynağı[[DÜDÜT:1) Bu katları yüklemek için senaryonuzu kopyalayın ve bunları devlet hesaplamaları sırasında arayınız.Kaynak belirli ısıları, tahmini verimliliği 2-5 oranında azaltacaksınız, bu da gerçek motor verileriyle yakından uyumlu.

Tam polinom entegrasyon olmadan hızlı bir yaklaşım için, lineer olarak sıcaklık ile değişen belirli bir ısı oranı kullanın: ⁇ (T) = 1.4 – 0.00005* (T – 300). Bu basit korelasyon trendi yakalar ve bir süre içinde birlikte uygulanabilir.

Bir modüler Fonksiyonlar-Driven Analizi

Bir monolithic senaryodan ziyade, bir MATLAB işlevinde döngü mantığını ortaya koyar. Bu, parametre süpürücüler için yeniden kullanılabilir hale getirir ve kod duplikasyonunu azaltır.

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

Değişken özel ısılar için, bir sıvı bayrağı ekleyin:

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

Yerdeki işlevle, bir sürücü senaryosu sıkıştırma oranlarıyla süpürebilir:

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--');

Bu modüler yaklaşım, 12'nin ötesinde artan sıkıştırma oranını hemen gösteriyor. arsa önemli bir ticarete işaret ediyor: verimlilik artıyor ama bir azalma oranına rağmen, MEP genellikle tipik ısı girişleri için r=10 civarında zirveler.

Tanık P-V ve T-S Diagrams

Görselleştirme, modeli doğrulamaya ve sonuçları iletişim kurmaya yardımcı olur. MATLAB'in arsa işlevleri yayın kaliteli diyagramlar üretmektedir.

Basınç-Volume (P-V) Diagram

Depres ve genişleme eğrileri takip eder:0)P v ⁇ = sabit[Döntgen: 1) Düz eğriler için bir dizi orta noktası:

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;

Kapalı alan net çalışmayı temsil eder. Sayısal olarak sayısal olarak kullanarak döngüyü entegre edebilirsiniz.(örneğin, üst eğriyi (kompresyon + genişleme) entegre etmek ve daha düşük eğriyi (sıratıcı sıralama hattı) çıkarmak için.Rezerv:0MATLAB belgeleri sayısal entegrasyon için).

SıcaklıkEntropy (T-S) Diagram

Sürekli belirli ısılarla ideal bir gaz için, entropi 1 eyalete göre değişir:

% 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);

Uygun hat stilleri ile dört işlem yapın. değişken için ⁇ , ısıtma eğrisinin sayısal entegrasyonunu kullanarak ısının eklendiğini ortaya koyar. T-S diagramı, işlem 1-2 ve 3-4 (vertical hatları) ve sürekli-volume ısıtma/rejeksiyonu (gerekli yollar) için, ısının altında alan ısının reddedildiği gibi; farkları net çalışmadır.

T-S arsasını geliştirmek için, sürekli baskı hatları (isobars) ideal gaz ilişkisinden kaynaklanmaktadır. Bu, döngüyü daha geniş bir termodinamik çerçevede bağlamaya yardımcı olur.

Hassasiyet Araştırmaları

Temel model işe alındığında, nüfuzlarını anlamak için girdi parametrelerini sistematik olarak değiştirir. Ortak çalışmalar şunları içerir:

  • [FONT=0]Compression oranı süpürücü[Dönder: Grafik verimliliği ve MEP ile R. Verimliliği düşük r ve sonra plasler hızla yükselir; MEP genellikle baskı artışı ve ısı eki arasında üst üste zirveler alır. Aynı rakamda hem ölçümler göstermek için alt kat kullanın.
  • Ek varyasyon: Değişen Q in modelleri farklı yakıt enerjileri veya hava yakıt oranları.Q in üst baskı ve iş çıkarır, ancak maddi sınırları aşabilir; P3 < 150 bar (15 MPa) tipik motor gücünü yansıtacak bir kısıtlama içerir.
  • [FONT:0)Mutfak ısı oranı etkisi[[DÜT:1): Aşağı ⁇ , düzeltmeyi göstermek için verimlilik azaltır. Bu, özellikle doğal gaz veya hidrojen gibi alternatif yakıtları analiz ederken önemlidir.

İki parametre süpürücü için, bir matrixte nested döngüler ve mağaza sonuçları kullanın. Sonra birleşik etkiyi görselleştirmek için bir kontras arsa veya renkli harita oluşturun:

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 (%)');

Bu tür hassasiyet haritaları, motor tasarımında ticaretten ayrılan analiz için paha biçilmezdir.Onlar, verimlilikin yüksek olduğu işletme bölgelerinin yüksek olduğunu tespit etmeye yardımcı olurlar.

Sonuçlarınızı Geçerlileştirmek

Mühendis kararları için model kullanmadan önce, bilinen ölçmelere karşı kontrol edin. Hava standart verimlilik formülü ⁇ = 1 - 1 / ^ ( ⁇ 1) Sürekli ⁇ davası için doğrudan geçerlilik sağlar. Ek olarak, hesaplanan üst basınç P3, yayınlanan motor haritaları veya ders kitaplarının verileriyle karşılaştırıldığında, Stone,ENFLT:0)Introduction to Internal Combustion Engines).

Başka bir doğrulama tekniği, [[0) enerji dengesi) kontrol etmektir: net çalışma artı reddedilen ısı, yüzen hassaslık içinde giriş ısısına eşit olmalıdır. Ayrıca, sıkıştırma ve genişleme süreçlerinin tatmin edici olduğunu doğrulamayı gerektirir.

Gerçekçi Etkileri için Modeli uzatın

Temel çerçeve gerçek dünya fenomenlerini içerecek şekilde genişletilebilir:

  • [FONT=0]Finite yanma süresi[[Dönetici:0))[değiştir | kaynağı değiştir] = 1 – s-a * (( ⁇ 0) / MATLAB'in [3 ) [Döneticileri ile birlikte enerji denklemini entegre etmek için kullanılır.
  • [FONT:0)Heat transfer kayıpları): InTOR a Woschni-tip korelasyon için uygun ısı geçişi silindir duvarlarına kadar. Bu, crank-kırık çözülme sıcaklık ve baskı gerektirir.
  • [FONT:0)Zizlenmiş gaz recirculation (EGR)[Dönetici: Belirli ısı oranını ve gaz kompozisyonunu dillendirmek için ayarlamanız için ayarlayabilirsiniz.TheurFLT:2).MATLAB Simscape Engine Silindir bloğu) bu tür çalışmalar için hazır bir bileşen sunar, ancak ayrıca belirli hava ısılarını ve gaz kompozisyonunu lineer olarak birleştirerek basitleştirilmiş bir dilselleştirme modelini de uygulayabilirsiniz.

Hızlı bir uzatma için, Wiebe işlevi ile başlayın. Örnek yaklaşımlar [[0)MATLAB Dosya Değişimi) ile birlikte bulunabilir. Bu basit ek, verimlilik ve zirve basıncı üzerindeki kıvılcım zamanlamasının etkisini ele alır.

Analizleri büyük bir iş akışına entegre etmek

Bir mühendislik projesinde, döngü analizi araç düzeyinde simülasyonlar veya rapor nesli ile beslenmelidir. MATLAB'in 03.30'u kullanın: devlet masalarını ve performans ölçümleri Excel'de daha fazla işlem için metrikleri dışa aktarmak.Çalışanlar işlevle kodunuzu bir PDF olarak birleştirin.

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

Ayrıca Python veya Simulink'ten MATLAB'ı diğer alt sistemlerle termodinamik modelini çiftleştirmek için MATLAB MATLAB isteyebilirsiniz.Bir hibrit araç çalışması için, Otto döngüsü verimliliği haritası Simulink'te bir gözlüğü tablo olarak ithal edilebilir.

Sonuç Sonuç Sonuç Sonuç Sonuç Sonuç Sonuç Sonuç

Bir hava standart Otto döngüsü modeli MATLAB'de inşa etmek, tasarım seçimlerinin nasıl etkilediğine dair güçlü bir platform sağlar.Sisteme göre motor termodinamikleri sistematik olarak tanımlamakta ve tasarım mantığı olmadan performans ölçümleri hesaplamak ve hesaplamak, net çalışma, termal verimlilik ve etkili bir baskı anlamına gelir, motor gelişimini etkileyen basit bir yaklaşıma anında bilgi verirsiniz.