发动机原理大作业

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

1、发动机原理大作业11051111051136题目：1、编制一个单轴涡轮喷气发动机设计点性能计算程序，要求：输入给定的发动机状态参数(最好以文件的形式输入)，可以正确的得到发动机各截面的总温、总压、质量流量，特征截面的流通面积，发动机推力、单位推力、耗油率。2、在增压比为1.5-60的范围内，计算单位推力与耗油率随增压比的变化关系。解：程序：#include#includedoublee=2.71828;intmain(){floatMa,H,Plc,YITAc,T4,YITAb,SIGMAb,YITAt,SIGMAc,YITAm,k,cp,kgzcp

2、g,R,Rg,Hu,Derta,qm3,Lmd2,Lmd3,qLmd2,A2,qLmd3,A3,Lmd4a,qLmd4a,A4a,Lmd9,qLmd9,A9,fl;floatT0/p0/V0/Tl/pl/SIGMAi/T2/p2/p3/T3/qm3a/Lc,tp4/qm4/T4a/p4a/qm4a/qm5/T5/Plt/p5/p9/T9/P9/Ma9/t9/c9/V9/qm9/F/Fs/sfczMa2=0.55/Ma3=0.3/Ma4a=l/ma9=l;FILE*fp=fopen("in.txt","r");if(!fp){printf("can'topenfile、n”);retu

3、rnJ;}while(!feof(fp)){fscanf(fp,H%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%化&Ma,&H,&Plc,&YITAc,&T4,&YITAb,&SIGMAb,&YITAt,&SIGMAc,&YITAm,&k,&cp,&kg,&cpg,&R,&Rg,&Hu,&Deg&qm3);}fclose(fp);if(H<=ll){T0=288.15-6.5*H;p0=101325*pow(l-H/44.308,5.2553);}if(H>11){T0=216.7;p0=0.227*pow(e,(ll-H)/6.338);}V0=Ma*s

4、qrt(k*R*1000*T0);Tl=T0*(l+(k-l)/2*Ma*Ma);pl=p0*pow(l+(k-l)/2*Ma*Ma,k/(k-l));if(Ma<=l)SIGMAi=0.97;if(Ma>l)SIGMAi=0.97*(l-0.075*pow(Ma-l,1.35));T2=T1;p2=SIGMAi*pl;Lmd2=sqrt(((k+l)/2*Ma2*Ma2)/(l+(k-l)/2*Ma2*Ma2));qLmd2=pow((k+l)/2,l/(k-l))*Lmd2*pow((l-(k-l)/(k+l)*Lmd2*Lmd2),l/(k-l));A2=qm3*sqrt(T2

5、)/qLmd2/p2/0.04042;p3=Plc*p2;T3=T2*(l+(pow(Plc,(k-l)/k)-l)/YITAc);Lmd3=sqrt(((k+l)/2*Ma3*Ma3)/(l+(k-l)/2*Ma3*Ma3));qLmd3=pow((k+l)/2,l/(l<-l))*Lmd3*pow((l-(k-l)/(k+l)*Lmd3*Lmd3),l/(k-l));qm3a=qm3*(l-Derta);A3=qm3a*sqrt(T3)/qLmd3/p3/0.04042;Lc=cp*(T3-T2);p4=SIGMAb*p3;T4=1500;f=(Cpg*T4-cp*T3)/(YI

6、TAb*Hu-cpg*T4);qm4=qm3*(l-Derta)*(l+f);T4a=T4*(((l-Derta)*(l+f)+cp*Derta*T3/cpg/T4)/((l-Derta)*(l+f)+Derta));fl=(cpg*T4-cp*T3)/(YITAb*Hu-cpg*T4)*qm3a/qm3;p4a=p4;qm4a=qm3*((l-Derta)*(l+f)+Derta);Lmd4a=sqrt(((k+l)/2*Ma4a*Ma4a)/(l+(k-l)/2*Ma4a*Ma4a));qLmd4a=pow((k+l)/2,l/(k-l))*Lmd4a*pow((l-(k-l)/

7、(k+l)*Lmd4a*Lmd4a),l/(k-l))；A4a=qm4a*sqrt(T4a)/qLmd4a/p4a/0.03938;Lc=cp*(T3-T2);qm5=qm4a;T5=T4a*(l-cp*(T3-T2)/(((l-Derta)*(l+f)+Derta)*YITAm*cpg*T4a));Plt=pow(l-(l-T5/T4a)/YITAt,-kg/(kg-l));p5=p4a/Plt;p9=SIGMAc*p5;T9=T5;P9=p0;Ma9=