涡扇发动机程序.docVIP

#includestdio.h #includemath.h void main() { //假设飞行条件// double Ma0=0.8,H=12; //发动机工作的一些参数// double B=0.57,Picl=3.6,Pich=9;/* B为涵道比，Picl为风扇的增压比，Pich为高压压气机增压比*/ double Pi=24.5,Tt4=1800,Ttab=2000;/*Pi为总增压比，Tt4为燃烧室出口温度，Ttab为加力燃烧室出口温度*/ double R=287.06,Rg=308.66; double Lcl,Lch,f,tm; /*风扇处每千克空气消耗的功*/ double Wc,W4,W4a,W4c; /*各截面流量*/ //预计的部件效率或损失系数// double di=0.97; /* 进气道总压恢复系数*/ double ncl=0.868; /* 风扇绝热效率*/ double nch=0.878; /*高压压气机绝热效率 */ double nb=0.98; /*主燃烧效率*/ double db=0.97; /* 主燃烧室总压恢复系数*/ double nth=0.89; /* 高压涡轮效率*/ double ntl=0.91; /* 低压涡轮效率*/ double dm=0.97; /*混合室总压恢复系数 */ double nbab=0.97; /*加力燃烧效率 */ double dbab=0.96; /* 加力燃烧室总压恢复系数*/ double de=0.98; /*尾喷管总压恢复系数 */ double nmh=0.98; /*高压轴机械效率 */ double nml=0.98; /* 低压轴机械效率*/ double Ct0=3; /* 相对功率提取系数,单位为kj/kg*/ double nmp=0.98; /*提取功率机械效率 */ double Cp=1.005; /* 空气的定压比热容，单位为 j/(kg。K)*/ double k=1.4; /* 空气的比热比*/ double Cpg=1.244; /* 燃气的定压比热容，单位为 j/(kg。K)*/ double kg=1.3; /*燃气的比热比 */ double Hu=42900; /*燃油的低热值，单位为kj/kg */ double d1=0.05; /* 高压涡轮的相对冷却空气量*/ double d2=0.05; /* 低压涡轮的相对冷却空气量*/ double bt=0.01; /* 飞机相对引气量*/ int sign=1; /*加力标记，若加力则为1，不加力则为其他值*/ //各截面参数的定义// double Pt0,P0,Tt0,T0,a0,c0; /*进口截面*/ double Pt2,P2,Tt2,T2; double Pt3,P3,Tt3,T3; double Pt4,P4,T4; double Pt4a,P4a,Tt4a; /*高压涡轮前截面*/ double Pt4c,P4c,Tt4c; /*低压涡轮前截面*/ double Pt5,P5,Tt5,T5; double Pt6,P6,Tt6,T6; double Pt7,P7,Tt7,T7; double Pt8,P8,Tt8,T8; double Pt9,P9,Tt9,T9,Ma9,a9,c9; double Pt22,P22,Tt22,T22; double Pt45,P45,Tt45,T45; //涡轮参数定义// double Pith,Pitl;/*高低压涡轮落压比*/ //涡轮后混合室参数定义// double Bm,Cp6,Pm; //加力燃烧室参数// double fab,f0; /*加力燃烧室油气比fab,总油气比f0*/ //发动机性能参数定义// double Fs,Fsab,sfc,sfcab; //0-0截面压力和温度// if(H11) { T0=(288.15-6.5*H); P0=101325*pow((1-H/44.308),5.2553); } else { T0=216.7; P0=22700*exp((H-11)/6.338); } 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)