航空发动机原理大作业.doc

《航空发动机原理大作业.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、航空发动机原理大作业某涡轮喷气发动机参数设计班级：0207102姓名:吴吉昌学号：0207102132010年12月20日作业题目：在海平面、静止状态、标准大气条件、最大工作状态时，对有关涡喷发动机的F，SFC的要求如下表所示，它们均采用收敛喷管，Vcol为压气机出口处的相对引气量，Vr为涡轮中的相对回气量。试选择有关参数，计算并画出Fs,SFC及Qma随兀k*(或T3*)的变化曲线，并确定满足性能要求的工作过程参数。发动机ABCDEFVcol0.030.02500.0300Vr0.020.0200.0200要FdN>=求SFC

2、Kg/dN.hr<=260012507002900265022000.951.01.20.961.21.1本报告选定C组参数进行设计讨论MATLAB源程序：function[Fs,SFC,Qma,f]=LIANXI2(Pk*)T0=288.2;P0=1.0133e5;C0=0;F=7000;C=340.3;Hn=42900000;di=1;nm=0.99;Vcol=0;Vr=0;sb=0.97;dt=0.98;db=0.96;nt*=0.93;nk*=0.83;k=1.4;R=287.3;k’=1.33;R’=288;Cp=k*

3、R/(k-1);Cp’=k’*R’/(k’-1);XH=(1250:10:1550);YH=[3137.53172.03206.83241.63276.83311.83346.93382.13417.53453.03488.73524.23560.03595.93632.03668.13704.23740.53777.03813.53850.23886.83923.53960.83997.44034.54071.34108.24245.24282.54219.7];Xh2=(350:10:850);Yh2=[350.5360.6

4、370.7380.8390.9401.4411.1421.3431.5441.7451.9462.1472.3482.6492.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.2855.2866.2877.2];Xh3=(1250:10:1550);Y

5、h3=[1336.61348.41360.21372.11383.91395.81407.71419.61431.61443.51455.51467.41479.41491.31503.31515.31527.31539.21551.31563.21575.31587.21599.11611.01623.01635.01647.31659.51671.81684.11696.5];T3*=(1250:5:1550);M0=C0/C;T0*=T0*(1+(k-1)*M0^2/2);P1*=di*P0*(1+(k-1)*M0^2/2

6、)^(k/(k-1));T1*=T0_;VTk*=T1**(Pk*^((k-1)/k)-1)/nk*;T2*=T1*+VTk*;P2*=P1**Pk*;P3*=P2**db;h3a*=interp1(Xh3,Yh3,T3*)*1000;h2a*=interp1(Xh2,Yh2,T2*)*1000;H3*=interp1(XH,YH,T3*)*1000;f=(h3a*-h2a*)./(sb.*Hn-H3*+h2a*);VTt*=Cp*VTk*./(Cp**(1-Vcol+f)*nm);T4*=T3*-VTt*;Pt*=(1-VTt

7、*./(T3**nt*)).^(-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)).^0.5;elseP9=P0;C9=(2*Cp’*T4**(1-P0./P9*).^((k’-1)/k’)).^0.5;endT9=T4*-C9.^2/2/Cp’;Fs=(1-Vcol+f+Vr).*(C9-C0./(1-Vcol+f+Vr)+R’*T9/C9*(1-P0./P9));SFC=3600*f*(1-

8、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