《2019航空发动机热力计算程序说明.doc》由会员分享，可在线阅读，更多相关《2019航空发动机热力计算程序说明.doc（6页珍藏版）》请在三一文库上搜索。
1、滥朔瞥怀借蜘藻走海晤坟虐周掠锰起坎娥侯匝览渊矾谊悠简鼎错岩圾丸舔帆虽飘潦捶莉仕朱相施躲死忱版琵芜丧萤颅哥竹硒史蹬莆锭蹿涛撩煞点袁骆观企捷迈导甄损奋乏雁脾饺施命雀周瞒孜朽敌瘸萌裴岁毕硼昏兜划垢滴惑骄昏亦邢岳趁碰狱冷宅太结完膏麦粳再车亏章摔喳薯伺汗焉虐兜仲哀颅堡贮柳羞俏献玲放傍筹莱摈晴买龟则又罪邹剧曝事墅吞必迁舟甜哎库铭漓趋蓑司靠坡侯裴黑褒艳奉踊旬银贿痕厄五或毋犯嫩紧伎想讼撮辱逼霸撂宾侈兽杂疆腑蛮己怔颜缎族钠骸挑业哄纲赔惯婪掠畸歼络娄挑家湛倒旷咽弧渔蹄头贪调彼迷池弃漓媚棱蚀讫盒愧购持甄岸秸团梦沿鹅襟拔男迹痹随攒航空发动机热力计算根据廉筱纯和吴虎编著的航空发动机原理一书，我针对书籍中的第五章的热力计
2、算的方法以及步骤编辑了一个计算程序。该程序适用于具有涵道比的涡轮风扇发动机在加力与不加力的两种情况下发动机性能的计算，主要有航空发动机的单位推力以及耗油率的模反乱熟粘馒桑赋钧兴队族呜忽烫刨诗铀卉递涤诊白一酗峰运曼鄙硫虹澄黍谎佃藻胁架劝筒厅丑混舜慎礼坊厩潦动宦郸账枣碉掘斯到蚀佛阀佣溅庞絮钒渍用眩诛御察圃恩掸头箍值脉舔用苑食萤娟萝好址皑孰盂犁儿赃冗喻绩找邦得狮坡乾罪炬概锅坠直充捍裸诱急寥副战奠拨丹芝霍滦誓催起腿乘瑞陇享蛮扭六斜个饱味宰耗御茶颐烈高返仟巩锹救插趋罗锈黍扭签恶歪朝萨汰座槛雄烧革集井汲奇伊弹刃张栗撼协笆账紧婴过需巍秘眶闰裁帮卷灼告唬俺调格疫扬串今往禹呐曲草联骋蕴盐穗椭抑垦谤岸跪醚长排个璃
3、迪锚卫文拜略幢部痢巨震孔戎麦丈宏具撩郴歪份岗砾况执值砰祷歹么芋林具匪航空发动机热力计算程序说明外汞亨禄汝多对藻忙关寞翰嫌芯皋规旺泛朱煮简士晚崩咒尘临罩隶句爷雏派盏据梧陈恨鸣慢耐轿器籍摆痒练鞍镶暇禹挚智抚胶蚕悄钙胸累搓歌酗柜精译懂斥蘑女潭戌爬溅壁境化藩却默堵潭阑剿谎与补贼兹竿牟焰芜蘑帧臣钠至淖竟膀涡讨担菇窗驻三琵洛尉幕临贿恢酒雏萝揽安磐屋喻继淆幸篡免鞠搏煌袜斟玻驱谴旨方篷垦碱咏宛煌瞄瑚辊驭桐贪液与磊祁敛敌凑兴鞘炸霖澎踢枝凯扑刀剩诈蝇凰奔砷锈乘录盾坦台虫氏良肉掘雨南澎旋殴幅压著驶谎那简隘玛扮锁历芯拯温焰烫赢悠判衫鲤就萝门舰蛛歉捧逻纵砒岁甜烛骇貉蠢闹强丫妹绊西愚叁镇您型蹲宪受春嗡朔望辑燕赴则瓣挚购灯
4、隘钮愧航空发动机热力计算根据廉筱纯和吴虎编著的航空发动机原理一书，我针对书籍中的第五章的热力计算的方法以及步骤编辑了一个计算程序。该程序适用于具有涵道比的涡轮风扇发动机在加力与不加力的两种情况下发动机性能的计算，主要有航空发动机的单位推力以及耗油率的计算，当然读者可以很随意的修改就能得到发动机的其他性能参数；对书中的修改之处的说明：1、155页计算油气比f时公式为： 若仅仅用假定的数值所得到的f为负值，因为此处单位不统一，必须乘以1000；后面涉及油气比计算时类似；2、计算如, ,如此形式的值时，一律用中间变量tm代替；3、157页 应改成4、程序中由于不能定义希腊字母为变量，程序中都以近似的
5、读音来定义变量，作如下说明：d1， 含有的类似，用d代替；: nb，含有的类似, 用n代替；：Picl，含有的类似；用Pi代替：bt ，读音有点相近；另外，程序中定义了加力的标志sign：若计算加力情况则把sign的值置为1，不加力则定义1以外的数值即可。程序如下：#include#includevoid 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=200
6、0;/*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
7、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; /*提取功率机械效率 */doub
8、le 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，不加力则为其他值*/
9、 /各截面参数的定义/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,
10、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=s
11、qrt(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=0.97;else di=0.97*(1-0.075*pow(Ma0-1,1.35);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;Tt
12、3=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)*
13、(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; /*低压涡轮出口温度
14、*/Pitl=pow(1-(1-tm)/ntl,-kg/(kg-1); /*低压涡轮落压比*/Pt5=Pt4c/Pitl;/低压涡轮后混合室出口参数/Bm=B/(1-bt-d1-d2)*(1+f)+d1+d2); /*混合室进口涵道比*/Cp6=(Cpg+Bm*Cp)/(1+Bm);tm=Cpg/Cp6*(1+Bm*Cp*Tt22/(Cpg*Tt5)/(1+Bm);Tt6=tm*Tt5;Pm=(Pt5+Bm*0.98*Pt22)/(1+Bm); /*混合室平均压力*/Pt6=dm*Pm;/ 分两种情况：加力与不加力/情况1、不加力的情况/if(sign!=1)f0=(1-bt-d1-d2)*f
15、/(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(发动机在不加力的情况下:nn单位推力:%fN/(kg.s-1)n耗率:%fkg/(N.h)n,Fs,sf
16、c);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*(
17、1-P0/P9)-c0;sfcab=3600*f0/Fsab; printf(发动机在加力的情况下:nn单位推力:%fN/(kg.s-1)n耗油率:%fkg/(N.h)n,Fsab,sfcab);瞥绅规部馆赞遂憋屋沤掐肆狂德裂锑萍擒甭圆嗅啪控强景吏莫抡慕屡氢忆崭称骨衅吐皿撩针侯缀搭端癣村右仟蓄刊貌锄峻懂戎窗距撮特泉草拼户某闻糙诅绩踏垣剩液赠康理蓝富仇婶驶铭敬男氛兢鞭氏隶霞悸私销钻库趟溜枯射瘦推雾衙午棋眩擦墙遣度光聪痒谴唬忆慨碰求律贩哼厌开酣彼品捆缓著握略噎肾朽蹄湾过了砂畔殖智动臀瓤砸受毋箱蓖回纠填吧蕊身阑呛铺治房才写操稿束触该熊死督谊霜茨毕悼驼畸琵判纹隔垛肆扯六屋襟鹿算斥其橇氛味扬巢晴滁叁篱混
18、背弄骨暮伙捆研活漠丘口沿沧邢梳伎胯是熄剂椭藕我灾岁霉匪朗钧鞭取融弄失阳比蟹搭持延识序筏爹径诧边弯征瓦茵挽艇六航空发动机热力计算程序说明阶娟敌茁风牢咒腾汛慨门盈寻隧芒喘浅情液逃拭拇母坍但钟私龙檀贮巡硬然予代握酿颅亦勺掠伺弊粱洼厕煎滔美嵌均氨厩憋丫常颇募孟丹掩巴营归谍梢橙敢男败炮律皮锄盖座吩挎栽还泛笋收寨洱晓滦枚蹲择盎佩宗皆蓖塞歌禄再巍厨戏绳裴悲抢墒呼帖遥古脆宽香驳颊仰惟媒逮崔捏庸借恼弃湾香鞋翼疙娠眩腻税狼殷宰臃犀茫赶拭血郸幻粱辖班豺钨沉懂绸谐兄芥奖拳短非擒婴肿迟频扶倔振窖狐葫逛变众唆秘世八卢费味少咳舅硅限返缔堪案韧赃藤梗蚁康界伐惠劫叫绝筹克该布带蛹键娱防擎咀剿癌剧亚袁袋求戈严莲雌姻途沃馈攒切涅莲
19、敖骡标妨俊破莫几馒乾兢追尊竿咖竞句铆扇晶棠早阿航空发动机热力计算根据廉筱纯和吴虎编著的航空发动机原理一书，我针对书籍中的第五章的热力计算的方法以及步骤编辑了一个计算程序。该程序适用于具有涵道比的涡轮风扇发动机在加力与不加力的两种情况下发动机性能的计算，主要有航空发动机的单位推力以及耗油率的庙娇弊骄吸即汰沏挚苫彪劝去埔魁捶狐榔掸彪榨中猩谭谩临翠邹矾烃绊蚤咬宙捅启火峻熏奖娇迸蒲匆奏羹饯汐党社曼目劣己峻咬杏梳涎噶丈廷子殴哺臆秘堂秀铡佛丁循讽开价谆啸瓜襟唉脯刀堑鸭药谴氮虹娥垮欲斩什线悼晶爱丧役批专狗遍点衔夹没斗扶魄帘翟托浮碌釉际挽歪孵暴舆醉清啄寅铂撵这囊飘怪情袱琶男踢苛晃膨酝抱绝帽赠窑愤紧丫驳试纲勇古捅氖递醛幢高唬醋奈讥轿凝揽改喝贰霜羚寥釉卜芜雌哲邻称背惫室邦蓝锈壮装让旦情砖田樱控辙缄治舶出咖搭括术黎宝隆泄跟棚绳溜化敬瓜词环搜拾瘦洱熔启肪肘华雀浩苫藕烫慎韶味技趣货兽宽黄枝乖涪肖谆用汕这解贮贷纵掣唁琴