(整理)航空发动机热力计算程序说明.docx

: .
航空发动机热力计算
根据廉筱纯和吴虎编著的《航空发动机原理》一书，我针对书籍中的第五章的热力计算p6,Pm;//加力燃烧室参数//
doublefab,fO;/*加力燃烧室油气比fab,总油气比fO*///发动机性能参数定义//
doubleFs,Fsab,sfc,sfcab;//O-O截面压力和温度//
if(H<11){
TO=(-*H);
PO=1O1325*pow((1-H/),);}else
{
TO=;
P0=22700*exp((H-11)/);}
a0=sqrt(k*R*T0);c0=Ma0*a0;
Pt0=P0*pow(1+(k-1)*Ma0*Ma0/2,k/(k-1));/*进口总压*/Tt0=T0*(1+(k-1)*Ma0*Ma0/2);/*进口总温*/
//进气道出口总温和总压//if(Ma0<=1)
{
di=;}
elsedi=*(1-*pow(Ma0-1,));Pt2=di*Pt0;
Tt2=Tt0;//风扇出口参数//
Pt22=Pt2*Picl;Tt22=Tt2*(1+(pow(Picl,(k-1)/k)-1)/ncl);
Lcl=Cp*(Tt22-Tt2);//高压压气机出口总温和总压//
Pt3=Pt22*Pich;Tt3=Tt22*(1+(pow(Pich,(k-1)/k)-1)/nch);
Lch=Cp*(Tt3-Tt22);//主燃烧室出口参数//
f=(Cpg*Tt4-Cp*Tt3)/(nb*Hu*1000-Cpg*Tt4);Pt4=Pt3*db;
Tt4=1800;//高压涡轮出口参数//
tm=((1-bt-d1-d2)*(1+f)+Cp*d1*Tt3/(Cpg*Tt4))/((1-bt-d1-d2)*(1+f)+d1);Tt4a=Tt4*tm;/*4a代表内外涵气流在高压涡轮前混合之后的截面，tm为中间值，此处为Tt4a/Tt4*/
Pt4a=Pt4;tm=1-Cp*(Tt3-Tt22)/(((1-bt-d1-d2)*(1+f)+d1)*nmh*Cpg*Tt4a);
Tt45=tm*Tt4a;Pith=pow(1-(1-tm)/nth,-kg/(kg-1));/*高压涡轮落压比*/
Pt45=Pt4a/Pith;//低压涡轮出口参数//
tm=((1-bt-d1-d2)*(1+f)+d1+d2*Cp*Tt3/(Cpg*Tt45))/((1-bt-d1-d2)*(1+f)+d1+d2);Tt4c=tm*Tt45;
Pt4c=Pt45;tm=1-(Cp*(Tt22-Tt2)+Ct0/nmp)*(1+B)/(nml*Cpg*Tt4c*(1-bt-d1-d2)*(1+f)+d1+d2);
Tt5=tm*Tt4c;/*低压涡轮出口温度*/Pitl=pow(1-(1-tm)/ntl,-kg/(kg-1));/*低压涡轮落压比*/
Pt5=Pt4c/Pitl;//低压涡轮后混合室出口参数//
Bm=B/((1-bt-d