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

发动机原理大作业110511 11051136题目：1、编制一个单轴涡轮喷气发动机设计点性能计算程序，要求：输入给定的发动机状态参数（最好以文件的形式输入），可以正确的得到发动机各截面的总温、总压、质量流量，特征截面的流通面积，发动机推力、单位推力、耗油率。2、在增压比为1.5-60的范围内，计算单位推力与耗油率随增压比的变化关系。解：程序：#include#includedouble 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,p9,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(cant open filen);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(h11)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(ma1) 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面积m2 油气比 速度系数 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.2fn,p2,t2,qm3,a2,lmd2,ma2);fprintf(fq,一n);fprintf(fq,燃烧室进口%7.0f %6.2f %5.2f %8.6f 0 %8.6f%4.2fn,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.2fn,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.2fn,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 pi118,fsn118,sfcn118,c136,p3n118,fn118,t3n118,p4n118,fn118,qm4n118,t4an118,p4an118,qm4an118,qm5n118,t5n118,pitn118,p5n118,t9n118,p9n118,ma9n118,t9n118,c9n118,v9n118,qm9n118;int i=0;file *fp1=fopen(in1.txt,r);if(!fp1)printf(cant open filen);return -1;while(!feof(fp1)fscanf(fp1,%f,&ci+);fclose(fp1);for(i=0;i118;i+) pii=ci;ma=c118;h=c119;yitac=c120;t4=c121;yitab=c122;sigmab=c123;yitat=c124;sigmac=c125;yitam=c126;k=c127;cp=c128;kg=c129;cpg=c130;r=c131;rg=c132;hu=c133;derta=c134;qm3=c135;if(h11)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(ma1) sigmai=0.97*(1-0.075*pow(ma-1,1.35);t2=t1;p2=sigmai*p1;for(i=0;i118;i+)p3ni=pii*p2;t3ni=t2*(1+(pow(pii,(k-1)/k)-1)/yitac);p4ni=sigmab*p3ni;t4=1500;fni=(cpg*t4-cp*t3ni)/(yitab*hu-cpg*t4);qm3a=qm3*(1-derta);qm4ni=qm3*(1-derta)*(1+fni);t4ani=t4*(1-derta)*(1+fni)+cp*derta*t3ni/cpg/t4)/(1-derta)*(1+fni)+derta);p4ani=p4ni;qm4ani=qm3*(1-derta)*(1+fni)+derta);qm5ni=qm4ani;t5ni=t4ani*(1-cp*(t3ni-t2)/(1-derta)*(1+fni)+derta)*yitam*cpg*t4ani);pitni=pow(1-(1-t5ni/t4ani)/yitat,-kg/(kg-1);p5ni=p4ani/pitni;p9ni=sigmac*p5ni;t9ni=t5ni;p9=p0;ma9ni=sqrt(2/(kg-1)*(pow(p9ni/p9,(kg-1)/kg)-1);t9ni=t9ni/(1+(kg-1)/2*ma9ni*ma9ni);c9ni=sqrt(kg*rg*1000*t9ni);v9ni=c9ni*ma9ni;qm9ni=qm5ni;fni=qm9ni*v9ni-qm3*v0;fsni=fni/qm3;sfcni=3600*(cpg*t4-cp*t3ni)*qm3*(1-derta)/fni/(yitab*