哈尔滨工业大学课程设计说明书（论文）
Ymh=0.99 %低压轴机械效率 Yml=0.99
%功率提取轴机械效率 Ymp=0.99 %空气定熵指数 k=1.4
%空气定压比热容 cp=1.005 %燃气定熵指数 kg=1.3
%燃气定压比热容 cpg=1.244 %气体常数 R=287 %燃油低热值 Hf=42900
%冷却高压涡轮空气量系数 Z=0.15 %飞机引气系数 beit=0.01
%0-0截面的温度和压力 T0=288.15-6.5*H
p0=1.0133*(1-H/44.308)^5.2553*10^5 a0=(k*R*T0)^0.5 c0=a0*Ma0
pt0=p0*(1+0.2*Ma0^2)^3.5 Tt0=T0*(1+0.2*Ma0^2) %进气道出口总温和总压 Xi=Ximax*(1-0.075*(Ma0-1)^1.35) pt2=pt0*Xi Tt2=Tt0 %风扇出口参数 pt22=pt2.*PAIcl
Tt22=Tt2.*(1+(PAIcl.^(0.4/k)-1)./Ycl) Lcl=cp.*(Tt22-Tt2)
%高压压气机出口总温和总压
哈尔滨工业大学课程设计说明书（论文）
pt3=pt22.*PAIch
Tt3=Tt22.*(1+(PAIch.^(0.4/k)-1)./Ych) Lch=cp.*(Tt3-Tt22) %燃烧室出口参数
f=(cpg*Tt4-cp.*Tt3)./(Yb*Hf-cpg*Tt4) pt4=pt3.*Xb %高压涡轮
tm=((1-beit-Z).*(1+f)+cp.*Z.*Tt3./(cpg*Tt4))./((1-beit-Z).*(1+f)+Z) Tt4a=Tt4.*tm pt4a=pt4
t4a=1-Lch./(((1-beit-Z).*(1+f)+Z).*Ymh.*cpg.*Tt4a) Tt45=Tt4a.*t4a
PAIth=(1-(1-t4a)./Yth).^(-(kg/(kg-1))) pt45=pt4a./PAIth %低压涡轮参数计算
t45=1-((Lcl+3.0/Ymp).*(1+B))./(((1-beit-Z).*(1+f)+Z).*Yml.*cpg.*Tt45) Tt5=Tt45.*t45
PAItl=(1-(1-t45)./Ytl).^(-(kg/(kg-1))) pt5=pt45./PAItl %混合室出口参数 Bm=B./((1-beit-Z).*(1+f)+Z) cp6=(cpg+Bm.*cp)./(1+Bm)
t5=(cpg./cp6).*(1+Bm.*cp.*Tt22./(cpg.*Tt5))./(1+Bm) Tt6=Tt5.*t5
pt6=Xm*(pt5+Bm.*0.98.*pt22)./(1+Bm) %加力燃烧室参数计算 Tt72=Tt6
fab=(1+f.*(1-beit-Z)./(1+B-beit)).*(cpg*Tt71-cp6.*Tt6)./(Yab*Hf-cpg*Tt71) f01=((1-beit-Z).*f+(1+B-beit).*fab)./(1+B) f02=((1-beit-Z).*f)./(1+B) pt71=pt6.*Xab1 pt72=pt6.*Xab2 %尾喷管出口参数 pt91=pt71.*Xc pt92=pt72.*Xc Tt91=Tt71 Tt92=Tt72
哈尔滨工业大学课程设计说明书（论文）
p9=p0
Ma91=((2/(kg-1)).*((pt91./p0).^((kg-1)/kg)-1)).^0.5 Ma92=((2/(kg-1)).*((pt92./p0).^((kg-1)/kg)-1)).^0.5 T91=Tt91.*(1+((kg-1)/2).*Ma91.^2).^(-1) T92=Tt92.*(1+((kg-1)/2).*Ma92.^2).^(-1) a91=(kg.*R.*T91).^0.5 a92=(kg.*R.*T92).^0.5 c91=a91.*Ma91 c92=a92.*Ma92 %发动机单位性能参数 Fsab=(1+f01-beit/(1+B)).*c91-c0 Fs=(1+f02-beit/(1+B)).*c92-c0 sfcab=3600.*f01./Fsab sfc=3600.*f02./Fs plot(PAIcl,sfcab) hold on plot(PAIcl,sfc)