2024年5月25日发(作者：宰父融雪)
#include
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;
double nmp=0.98;
double Cp=1.005;
double k=1.4;
double Cpg=1.244;
double kg=1.3;
double Hu=42900;
double d1=0.05;
double d2=0.05;
double bt=0.01;
int sign=1;
/* 尾喷管总压恢复系数 */
/* 高压轴机械效率 */
/* 低压轴机械效率 */
/* 相对功率提取系数 ,单位为 kj/kg*/
/* 提取功率机械效率 */
/* 空气的定压比热容，单位为
/* 空气的比热比 */
/* 燃气的定压比热容，单位为 j/（kg 。 K）*/
/* 燃气的比热比 */
/* 燃油的低热值，单位为 kj/kg */
/* 高压涡轮的相对冷却空气量 */
/* 低压涡轮的相对冷却空气量 */
/* 飞机相对引气量 */
/* 加力标记，若加力则为 1，不加力则为其他值 */
j/（kg 。 K）*/
// 各截面参数的定义 //
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(H<11)
{
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));
/* 进口总压 */
2024年5月25日发(作者：宰父融雪)
#include
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;
double nmp=0.98;
double Cp=1.005;
double k=1.4;
double Cpg=1.244;
double kg=1.3;
double Hu=42900;
double d1=0.05;
double d2=0.05;
double bt=0.01;
int sign=1;
/* 尾喷管总压恢复系数 */
/* 高压轴机械效率 */
/* 低压轴机械效率 */
/* 相对功率提取系数 ,单位为 kj/kg*/
/* 提取功率机械效率 */
/* 空气的定压比热容，单位为
/* 空气的比热比 */
/* 燃气的定压比热容，单位为 j/（kg 。 K）*/
/* 燃气的比热比 */
/* 燃油的低热值，单位为 kj/kg */
/* 高压涡轮的相对冷却空气量 */
/* 低压涡轮的相对冷却空气量 */
/* 飞机相对引气量 */
/* 加力标记，若加力则为 1，不加力则为其他值 */
j/（kg 。 K）*/
// 各截面参数的定义 //
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(H<11)
{
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));
/* 进口总压 */