航空发动机热力计算程序说明.doc

、不加力的情况// if(sign!=1) { f0=(1-bt-d1-d2)*f/(1+B); /* 总油气比*/ Pt7=0.98*Pt6; Tt7=Tt6; Tt9=Tt7; Pt9=de*Pt7; P9=P0; Ma9=sqrt(2/(kg-1)*(pow(Pt9/P9,(kg-1)/kg)-1)); T9=Tt9/(1+(kg-1)*Ma9*Ma9/2); a9=sqrt(kg*Rg*T9); c9=Ma9*a9; Fs=(1+f0-bt/(1+B))*(c9+R*T9/c9*(1-P0/P9))-c0; sfc=3600*f*(1-bt-d1-d2)/(Fs*(1+B)); printf(" 发动机在不加力的情况下:\n\ n 单位推力:%fN/(kg.s-1)\ n 耗率:%fkg/(N.h)\n",Fs,sfc); } else { Tt7=2000; fab=(1+f*(1-bt-d1-d2)/(1+B-bt))*((Cpg*Tt7-Cp6*Tt6)/(nbab*Hu*1000-Cpg*Tt7)); f0=((1-bt-d1-d2)*f+(1+B+bt)*fab)/(1+B); Pt7=0.96*Pt6; Tt9=Tt7; Pt9=de*Pt7; P9=P0; Ma9=sqrt(2/(kg-1)*(pow(Pt9/P9,(kg-1)/kg)-1)); T9=Tt9/(1+(kg-1)*Ma9*Ma9/2); a9=sqrt(kg*Rg*T9); c9=Ma9*a9; Fsab=(1+f0-bt/(1+B))*(c9+R*T9/c9*(1-P0/P9))-c0; sfcab=3600*f0/Fsab; printf(" 发动机在加力的情况下:\n\n 单位推力:%fN/(kg.s-1)\n 耗油率:%fkg/(N.h)\n",Fsab,sfcab); }}