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

航空发动机热力计算 根据廉筱纯和吴虎编著的《航空发动机原理》一书，我针对书籍中的第五章的热力计算的方法以及步骤编辑了一个计算程序。该程序适用于具有涵道比的涡轮风扇发动机在加力与不加力的两种情况下发动机性能的计算，主要有航空发动机的单位推力以及耗油率的计算，当然读者可以很随意的修改就能得到发动机的其他性能参数； 对书中的修改之处的说明： 1、155页计算油气比f时公式为： QUOTE 若仅仅用假定的数值所得到的f为负值，因为此处单位不统一， QUOTE 必须乘以1000；后面涉及油气比计算时类似； 2、计算如 QUOTE , QUOTE , QUOTE ,如此形式的值时，一律用中间变量tm代替； 3、157页 QUOTE 应改成 4、程序中由于不能定义希腊字母为变量，程序中都以近似的读音来定义变量，作如下说明： = 1 \* GB3 ① QUOTE ：d1， 含有 QUOTE 的类似，用d代替 QUOTE ； = 2 \* GB3 ② QUOTE : nb，含有的类似, 用n代替； = 3 \* GB3 ③ QUOTE ：Picl，含有 QUOTE 的类似；用Pi代替 QUOTE = 4 \* GB3 ④ QUOTE ：bt ，读音有点相近； 另外，程序中定义了加力的标志sign：若计算加力情况则把sign的值置为1，不加力则定义1以外的数值即可。 程序如下： #includestdio.h #includemath.h void main() { //假设飞行条件// double Ma0=1.6,H=11; //发动机工作的一些参数// double B=0.4,Picl=3.8,Pich=4.474;/* B为涵道比，Picl为风扇的增压比，Pich为高压压气机增压比*/ double Pi=17,Tt4=1800,Ttab=2000;/*Pi为总增压比，Tt4为燃烧室出口温度，Ttab为加力燃烧室出口温度*/ double R=287.06,Rg=287.4; 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=1005; /* 空气的定压比热容，单位为 j/(kg。K)*/ double k=1.4; /* 空气的比热比*/ double Cpg=1244; /* 燃气的定压比热容，单位为 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=2; /*加力标记，若加力则为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; /*高压涡轮前截面