发动机原理大作业单轴涡轮喷气发动机设计点性能计算程序

发动机原理大作业
110511 11051136
题目：1、编制一个单轴涡轮喷气发动机设计点性能计算程序，要求：输入给定的发动机状态参数（最好以文件的形式输入），可以正确的得到发动机各截面的总温、总压、质量流量，特征截面的流通面积，发动机推力、单位推力、耗油率。

2、在增压比为1.5-60的范围内，计算单位推力与耗油率随增压比的变化关系。

解：程序：
#include<stdio.h>
#include<math.h>
double e=2.71828;
int main()
{float
Ma,H,PIc,YITAc,T4,YITAb,SIGMAb,YITAt,SIGMAc,YITAm,k,cp,kg,cpg,R,Rg,Hu,Derta,qm3, Lmd2,Lmd3,qLmd2,A2,qLmd3,A3,Lmd4a,qLmd4a,A4a,Lmd9,qLmd9,A9,f1;
float
T0,p0,V0,T1,p1,SIGMAi,T2,p2,p3,T3,qm3a,Lc,f,p4,qm4,T4a,p4a,qm4a,qm5,T5,PIt,p5,p 9,T9,P9,Ma9,t9,c9,V9,qm9,F,Fs,sfc,Ma2=0.55,Ma3=0.3,Ma4a=1,ma9=1;
FILE *fp=fopen("in.txt","r");
if(!fp)
{printf("can't open file\n");return -1;}
while(!feof(fp)){fscanf(fp,"%f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f % f %f",
&Ma,&H,&PIc,&YITAc,&T4,&YITAb,&SIGMAb,&YITAt,&SIGMAc,&YITAm,&k,&cp, &kg,&cpg,&R,&Rg,&Hu,&Derta,&qm3);}
fclose(fp);
if(H<=11)
{T0=288.15-6.5*H;
p0=101325*pow(1-H/44.308,5.2553);}
if(H>11)
{T0=216.7;
p0=0.227*pow(e,(11-H)/6.338);}
V0=Ma*sqrt(k*R*1000*T0);
T1=T0*(1+(k-1)/2*Ma*Ma);
p1=p0*pow(1+(k-1)/2*Ma*Ma,k/(k-1));
if(Ma<=1) SIGMAi=0.97;
if(Ma>1) SIGMAi=0.97*(1-0.075*pow(Ma-1,1.35));
T2=T1;
p2=SIGMAi*p1;
Lmd2=sqrt(((k+1)/2*Ma2*Ma2)/(1+(k-1)/2*Ma2*Ma2));
qLmd2=pow((k+1)/2,1/(k-1))*Lmd2*pow((1-(k-1)/(k+1)*Lmd2*Lmd2),1/(k-1));
A2=qm3*sqrt(T2)/qLmd2/p2/0.04042;
p3=PIc*p2;
T3=T2*(1+(pow(PIc,(k-1)/k)-1)/YITAc);
Lmd3=sqrt(((k+1)/2*Ma3*Ma3)/(1+(k-1)/2*Ma3*Ma3));
qLmd3=pow((k+1)/2,1/(k-1))*Lmd3*pow((1-(k-1)/(k+1)*Lmd3*Lmd3),1/(k-1));
qm3a=qm3*(1-Derta);
A3=qm3a*sqrt(T3)/qLmd3/p3/0.04042;
Lc=cp*(T3-T2);
p4=SIGMAb*p3;
T4=1500;
f=(cpg*T4-cp*T3)/(YITAb*Hu-cpg*T4);
qm4=qm3*(1-Derta)*(1+f);
T4a=T4*(((1-Derta)*(1+f)+cp*Derta*T3/cpg/T4)/((1-Derta)*(1+f)+Derta));
f1=(cpg*T4-cp*T3)/(YITAb*Hu-cpg*T4)*qm3a/qm3;
p4a=p4;
qm4a=qm3*((1-Derta)*(1+f)+Derta);
Lmd4a=sqrt(((k+1)/2*Ma4a*Ma4a)/(1+(k-1)/2*Ma4a*Ma4a));
qLmd4a=pow((k+1)/2,1/(k-1))*Lmd4a*pow((1-(k-1)/(k+1)*Lmd4a*Lmd4a),1/(k-1 ));
A4a=qm4a*sqrt(T4a)/qLmd4a/p4a/0.03938;
Lc=cp*(T3-T2);
qm5=qm4a;
T5=T4a*(1-cp*(T3-T2)/(((1-Derta)*(1+f)+Derta)*YITAm*cpg*T4a));
PIt=pow(1-(1-T5/T4a)/YITAt,-kg/(kg-1));
p5=p4a/PIt;
p9=SIGMAc*p5;
T9=T5;
P9=p0;
Ma9=sqrt(2/(kg-1)*(pow(p9/P9,(kg-1)/kg)-1));
Lmd9=sqrt(((k+1)/2*ma9*ma9)/(1+(k-1)/2*ma9*ma9));
qLmd9=pow((k+1)/2,1/(k-1))*Lmd9*pow((1-(k-1)/(k+1)*Lmd9*Lmd9),1/(k-1));
t9=T9/(1+(kg-1)/2*Ma9*Ma9);
c9=sqrt(kg*Rg*1000*t9);
V9=c9*Ma9;
qm9=qm5;
A9=qm9*sqrt(T9)/qLmd9/p9/0.03938;
F=qm9*V9-qm3*V0;
Fs=F/qm3;
sfc=3600*(cpg*T4-cp*T3)*qm3*(1-Derta)/F/(YITAb*Hu-cpg*T4);
FILE *fq=fopen("out1.txt","w");
while(!feof(fq))
{fprintf(fq,"┌一────┬────┬────┬────┬────┬────┬────┬──┐\n");
fprintf(fq,"│截面│压力pa │温度K │流量kg/s│面积m^2
│油气比│速度系数│Ma │\n");
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│大气环境│%7.2f│%6.2f │──│──│──│││\n",p0,T0);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│滞止参数│%7.2f│%6.2f │──│──│──│││\n",p1,T1);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│压气机进口│%7.2f│%6.2f │%5.2f │%8.6f│0 │%8.6f│%4.2f│\n",p2,T2,qm3,A2,Lmd2,Ma2);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│燃烧室进口│%7.0f │%6.2f │%5.2f │%8.6f│0 │%8.6f│%4.2f│\n",p3,T3,qm3a,A3,Lmd3,Ma3);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│涡导进口│%7.0f │%6.2f │%8.5f│──│%8.6f│││\n",p4,T4,qm4,f);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│涡导出口│%7.0f │%6.2f │%8.5f│%8.6f│%8.6f│%8.6f │%4.2f│\n",p4a,T4a,qm4a,A4a,f1,Lmd4a,Ma4a);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│涡轮出口│%7.1f│%6.2f │%8.5f│──│%8.6f│││\n",p5,T5,qm5,f1);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│尾喷管出口│%7.1f│%6.2f │%8.5f│%8.6f│%8.6f│%8.6f │%4.2f│\n",p9,T9,qm9,A9,f1,Lmd9,ma9);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│飞行速度│%7.3f│││││││\n",V0);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│飞行Ma │%4.1f │││││││\n",Ma);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│飞行高度│%2.0f │││
││││\n",H);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│推力│%8.2f│││││││\n",F);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│单位推力│%7.3f│││││││\n",Fs);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│耗油率│%8.6f│││││││\n",sfc);
fprintf(fq,"├一────┼────┼────┼────┼────┼────┼────┼──┤\n");
fprintf(fq,"│排气速度│%8.3f│││││││\n",V9);
fprintf(fq,"└一────┴────┴────┴────┴────┴────┴────┴──┘\n");
break;}
fclose(fq);
float
PI[118],Fsn[118],sfcn[118],c[136],p3n[118],Fn[118],T3n[118],p4n[118],fn[118],qm4n [118],T4an[118],p4an[118],qm4an[118],qm5n[118],T5n[118],PItn[118],p5n[118],T9n [118],p9n[118],Ma9n[118],t9n[118],c9n[118],V9n[118],qm9n[118];
int i=0;
FILE *fp1=fopen("in1.txt","r");
if(!fp1)
{printf("can't open file\n");return -1;}
while(!feof(fp1)){fscanf(fp1,"%f",&c[i++]);}
fclose(fp1);
for(i=0;i<118;i++) PI[i]=c[i];
Ma=c[118];
H=c[119];
YITAc=c[120];
T4=c[121];
YITAb=c[122];
SIGMAb=c[123];
YITAt=c[124];
SIGMAc=c[125];
YITAm=c[126];
k=c[127];
cp=c[128];
kg=c[129];
cpg=c[130];
R=c[131];
Rg=c[132];
Hu=c[133];
Derta=c[134];
qm3=c[135];
if(H<=11)
{T0=288.15-6.5*H;
p0=101325*pow(1-H/44.308,5.2553);}
if(H>11)
{T0=216.7;
p0=0.227*pow(e,(11-H)/6.338);}
V0=Ma*sqrt(k*R*1000*T0);
T1=T0*(1+(k-1)/2*Ma*Ma);
p1=p0*pow(1+(k-1)/2*Ma*Ma,k/(k-1));
if(Ma<=1) SIGMAi=0.97;
if(Ma>1) SIGMAi=0.97*(1-0.075*pow(Ma-1,1.35));
T2=T1;
p2=SIGMAi*p1;
for(i=0;i<118;i++)
{p3n[i]=PI[i]*p2;
T3n[i]=T2*(1+(pow(PI[i],(k-1)/k)-1)/YITAc);
p4n[i]=SIGMAb*p3n[i];
T4=1500;
fn[i]=(cpg*T4-cp*T3n[i])/(YITAb*Hu-cpg*T4);
qm3a=qm3*(1-Derta);
qm4n[i]=qm3*(1-Derta)*(1+fn[i]);
T4an[i]=T4*(((1-Derta)*(1+fn[i])+cp*Derta*T3n[i]/cpg/T4)/((1-Derta)*(1+fn[i])+D erta));
p4an[i]=p4n[i];
qm4an[i]=qm3*((1-Derta)*(1+fn[i])+Derta);
qm5n[i]=qm4an[i];
T5n[i]=T4an[i]*(1-cp*(T3n[i]-T2)/(((1-Derta)*(1+fn[i])+Derta)*YITAm*cpg*T4an[i] ));
PItn[i]=pow(1-(1-T5n[i]/T4an[i])/YITAt,-kg/(kg-1));
p5n[i]=p4an[i]/PItn[i];
p9n[i]=SIGMAc*p5n[i];
T9n[i]=T5n[i];
P9=p0;
Ma9n[i]=sqrt(2/(kg-1)*(pow(p9n[i]/P9,(kg-1)/kg)-1));
t9n[i]=T9n[i]/(1+(kg-1)/2*Ma9n[i]*Ma9n[i]);
c9n[i]=sqrt(kg*Rg*1000*t9n[i]);
V9n[i]=c9n[i]*Ma9n[i];
qm9n[i]=qm5n[i];
Fn[i]=qm9n[i]*V9n[i]-qm3*V0;
Fsn[i]=Fn[i]/qm3;
sfcn[i]=3600*(cpg*T4-cp*T3n[i])*qm3*(1-Derta)/Fn[i]/(YITAb*Hu-cpg*T4);} FILE *fq1=fopen("out2.txt","w");
while(!feof(fq1))
{fprintf(fq1,"┌一───┬一─────┬一─────┐\n");
fprintf(fq1,"│增压比│单位推力N/kg│耗油率kg/N.h│\n");
for(i=0;i<118;i++)
{fprintf(fq1,"├一───┼一─────┼一─────┤\n");
fprintf(fq1,"│%4.1f │%8.3f │%8.6f │\n",PI[i],Fsn[i],sfcn[i]);}
fprintf(fq1,"└一───┴──────┴──────┘\n");
break;}
fclose(fq1);
return 0;
}
输入文件：
输出文件：。