当前位置：文档之家› 发动机原理大作业单轴涡轮喷气发动机设计点性能计算程序

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

发动机原理大作业

110511 11051136

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

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

解：程序：

#include

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");