涡扇发动机程序.docxVIP

#in cludestdio.h # includemath.h void main() { 〃假设飞行条件〃 double Ma0=0.8rH=12; 〃发动机工作的一些参数〃 double B=0?57,Picl=36Pich=9;/“ B为涵道比,Piel为风扇的増压比.Pich为高压压气机增压比*/ double Pi=24.5Jt4=1800zTtab=2000;/*Pi为总增压比.Tt4为燃烧室出口温度.Ttab为加力燃烧室出I I温度*/ double R=287.06,Rg=308.66; double Lcl,Lch,f,tm; 厂风扇处每干克空气消耗的功*7 double Wc,W4,W4a,W4c; 厂乞截面流量食/ 〃预计的部件效率或损失系数〃 double di=0.97; 厂进气道总压恢复系数*/ double ncl=0.868; /?风扇绝热效率*/ double nch=0.87 厂高压压气机绝热效率?/ double nb=0.98; 厂主燃烧效率拿/ double db=0.97; /*主燃烧室总压恢复系数次/ double nth=0.89; /*高压涡轮效率?/ double ntl=0.91; 厂低压涡轮效率?/ double dm=0.97; 厂混合室总压恢复系数V double nbab=0.97; 严加力燃烧效率?/ double dbab=0.96; /*加力燃烧室总压恢圮系数袁/ double de=0?98; 厂尾喷管总压恢复系数7 double nmh=0.98; 厂高压轴机械效率7 double nml=0.98; /*低压轴机械效率?/ double CtO=3; 厂相对功率提取系数,单位为kj/kg*/ double nmp=0.98; 厂提取功率机械效率T 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 dl=0.05; /*高压涡轮的相对冷却空气甘F/ double d2=0.05; 厂低压涡轮的相对冷却空气戢?/ double bt=O.Ol; /?飞机相对引气量*7 int sign=l; 厂加力标记，若加力则为1，不加力则为其他值?/ 〃各栽面参数的定义〃 double PtOfPOJtOzTO,aO,cO;厂进口截面?/ double Pt2fP2Jt2J2; double Pt3rP3Jt3J3; double Pt4,P4,T4; double Pt4a,P4a；K4a;广高压涡轮前裁面專/ double Pt4c/P4c/Tt4c;厂低压涡轮前截面拿/ double Pt5rP5Jt5,T5; double Pt6rP6Jt6J6; double Pt7rP7Jt7zT7; double PtP8,Tt8zT8; double Pt9,P9；Tt9,T9,Ma9?a9?c9; double Pl22,P22,Tt22,T22; double Pt45,P45,Tt45J45; 〃涡轮参数定义〃 double PigPitl;厂高低压涡轮落压比?/ 〃涡轮后混合室参数定义〃 double BmzCp6zPm; 〃加力燃饶室参数〃 double fabJO; 厂加力燃烧室油气比fab,总油气比fO4/ 〃发动机性能参数定义〃 double Fs,Fsab,sfc?sfcab; //0-0截面压力和温度〃 if(Hll) { T0=(288?15?6.5?H); P0=101325*pow((l-HM308),5.2553); } else { T0=216.7; P0=22700*exp((H-ll)/6.338); } aO=sqrt(ksR#TO); cO=MaO*aO; /?进口总压*7PtO 二 P0Sow(l+(k?l)?Ma(TMa02k/(k?：L)); printf(H当前高度:\n\n 温度:%fk\n 压强:％fk\n 声速:%fm/s\n 速度:%fm/s\n 进口总ffi:%fPa\n /?进口总压*7 //进气道出U1总船和总压// if(MaO=l) { di=0.97; } elsedi=0.97*(l-0.075*pow(Ma0-14.35)); Pt2=diMPtO; Tt2=TtO; printf(进气道出口总温和总压:\n\n温度:％fk\n压力:%fPa\n*\Tt2zPt2); //风扇岀口参数〃 P