航空發(fā)動(dòng)機(jī)原理大作業(yè)

1、航空發(fā)動(dòng)機(jī)原理大作業(yè)某渦輪噴氣發(fā)動(dòng)機(jī)參數(shù)設(shè)計(jì)班級(jí)：0207102姓名：吳吉昌學(xué)號(hào)：0207102132010年 12月 20 日作業(yè)題目：在海平面、靜止?fàn)顟B(tài)、標(biāo)準(zhǔn)大氣條件、最大工作狀態(tài)時(shí)，對有關(guān)渦噴發(fā)動(dòng)機(jī)的F, SFC勺要求如下表所示，它們均采用收斂噴管，Vcol為壓氣機(jī)出口處的相對引氣量，Vr為渦輪中的相對回氣量。試選擇有關(guān)參數(shù)， 計(jì)算并畫出Fs,SFC及Qm隨兀k*（或Ts*）的變化曲線，并確定滿足性能要 求的工作過程參數(shù)。發(fā)動(dòng)機(jī)ABCDEFVcol0.030.02500.0300Vr0.020.0200.0200要 FdN=26001250700290026502200求 SFC Kg/

2、dN.hrv二0.951.01.20.961.21.1本報(bào)告選定C組參數(shù)進(jìn)行設(shè)計(jì)討論MATLA源程序：fun ctio n Fs,SFC,Qma,f=LIANXI2(Pk*)T0=288.2;P0=1.0133e5;C0=0;F=7000;C=340.3;H n=42900000;di=1; nm=0.99;Vcol=0;Vr=0;sb=0.97;dt=0.98;db=0.96; nt*=0.93 ;n k*=0.83;k=1.4;R=287.3;k =1.33;R =288;Cp=k*R/(k- 1);Cp =k *R /(k -1);XH=(1250:10:1550);YH=3137.53

3、172.03206.83241.63276.83311.83346.93382.13417.53453.03488.7 3524.2 3560.0 3595.9 3632.0 3668.1 3704.2 3740.5 3777.0 3813.53850.2 3886.8 3923.5 3960.8 3997.4 4034.5 4071.3 4108.2 4245.2 4282.54219.7;Xh2=(350:10:850);Yh2=350.5 360.6 370.7 380.8 390.9 401.4 411.1 421.3 431.5 441.7 451.9462.1472.3482.64

4、92.9503.2513.5523.8534.1544.5554.9565.3575.7586.2596.6607.1617.6628.2638.8649.3660.0670.6681.3692.0702.7713.4724.2735.0745.9756.7767.6778.5789.4800.3811.6822.2833.2844.2 855.2 866.2 877.2;Xh3=(1250:10:1550);Yh3=1336.6 1348.4 1360.2 1372.1 1383.9 1395.8 1407.7 1419.6 1431.6 1443.5 1455.5 1467.4 1479.

5、4 1491.3 1503.3 1515.3 1527.3 1539.2 1551.31563.2 1575.3 1587.2 1599.1 1611.0 1623.0 1635.0 1647.3 1659.5 1671.8 1684.1 1696.5;T3*=(1250:5:1550);M0=C0/C;T0*=T0*(1+(k-1)*M0A2/2);P1*=di*P0*(1+(k-1)*M0A2/2)A(k/(k-1);T1*=T0_;VTk*=T1*(Pk*A(k-1)/k)-1)/nk*;T2*=T1*+VTk*;P2*=pi*pk*；P3*=P2*db；h3a*=i nterp1(Xh

6、3,Yh3,T3*)*1000;h2a*=i nterp1(Xh2,Yh2,T2*)*1000;H3*=i nterp1(XH,YH,T3*)*1000;f=(h3a*-h2a*)./(sb.*H n-H3*+h2a*);VTt*=Cp*VTk*./(Cp*(1-Vcol+f)* nm);T4*=T3*-VTt*;Pt*=(1-VTt*./(T3* nt*).A(-k*/(k*-1);P4*=P3*./Pt*;P9*=P4*dt;a=P9*/P0;if(a=1.85)P9=P9*./1.85;C9=(2*k *R *T4 * ./(k +1).A0.5;elseP9=P0;C9=(2*Cp *

7、T4 *(1-P0./P9*).A(k - 1)/k ).人0.5;endT9=T4*- C9.A2/2/Cp ；-P0./P9);Fs=(1-Vcol+f+Vr).*(C9-C0./(1-Vcol+f+Vr)+R *T9/C9*(1SFC=3600*f*(1-Vcol)./Fs;Qma=F./Fs;a1,b1,c1,d1=LIANXI2(9);a2,b2,c2,d2=LIANXI2(10);a3,b3,c3,d3=LIANXI2(11);T3*=(1250:5:1550);figure(1)plot(T3*,a1,g*,T3*,a2,ko,T3*,a3,bx);xlabel(T3*);yla

8、bel(Fs);title( 渦噴發(fā)動(dòng)機(jī)性能參數(shù) Fs在Pk* 一定時(shí)隨T3*的變化情況);figure(2)plot(T3*,b1,g*,T3*,b2,ko,T3*,b3,bx);xlabel(T3*);ylabel(SFC);title( 渦噴發(fā)動(dòng)機(jī)性能參數(shù) SFC在Pk* 一定時(shí)隨T3*的變化情況);figure(3)plot(T3*,c1,g*,T3*,c2,ko,T3*,c3,bx);xlabel(T3*);ylabel(Qma);title( 渦噴發(fā)動(dòng)機(jī)性能參數(shù) Qma在 Pk* 一定時(shí)隨T3*的變化情況);figure(4)plot(T3*,d1,g*,T3*,d2,ko,T3*

9、,d3,bx);xlabel(T3*);ylabel(f);title( 渦噴發(fā)動(dòng)機(jī)性能參數(shù)f在Pk* 一定時(shí)隨T3*的變化情況);運(yùn)行結(jié)果如下渦噴發(fā)動(dòng)機(jī)性能券數(shù)刊在定時(shí)隨T3的變化情況940 11111920x e900x-些*680 r翼-$孳 *060 r&40 r-S2D異f900 -陽曲-7601111111250 13DO15G01550T3_渦噴岌動(dòng)機(jī)性能蔘數(shù)EFC在Pk_走時(shí)隨T3的變化情況0 105D 1丄co0.0950 09n.oes0.115 |111110.11o-rD.081Ll11112501300135014001 陌 D16001550T3% a嘰-%嘰渦噴發(fā)動(dòng)機(jī)性能參數(shù)Qma在Pk_定時(shí)隨73的變化情況9.49.298.88.6e 8.48.287.87.67.41250130013501400145015001550T3_渦噴發(fā)動(dòng)機(jī)性能參數(shù)HPk_定時(shí)隨T3的變化情況0.030.028 ?0.0260.0240.0220.020.0180.0161250130013501400145015001550T316附圖：飛行條件下 兀k*、T3*對Fs、SFC勺影響0.08500600700800電於牛頓-