發(fā)動(dòng)機(jī)原理大作業(yè)

1、發(fā)動(dòng)機(jī)原理大作業(yè)110511 11051136題目：1、編制一個(gè)單軸渦輪噴氣發(fā)動(dòng)機(jī)設(shè)計(jì)點(diǎn)性能計(jì)算程序，要求：輸入給定的發(fā)動(dòng)機(jī)狀態(tài)參數(shù)（最好以文件的形式輸入），可以正確的得到發(fā)動(dòng)機(jī)各截面的總溫、總壓、質(zhì)量流量，特征截面的流通面積，發(fā)動(dòng)機(jī)推力、單位推力、耗油率。2、在增壓比為1.5-60的范圍內(nèi)，計(jì)算單位推力與耗油率隨增壓比的變化關(guān)系。解：程序：#include<stdio.h>#include<math.h>double e=2.71828;int main()float Ma,H,PIc,YITAc,T4,YITAb,SIGMAb,YITAt,SIGMAc,YITAm,

2、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)pri

3、ntf("can't 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,&am

4、p;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

5、;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)

6、*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-Der

7、ta)*(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=

8、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-q

9、m3*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 油氣比 速度系數(shù) Ma n");fprintf(fq,"一n");fprintf(fq," 大氣環(huán)境 %7.2f %6.2f n",p0

10、,T0);fprintf(fq,"一n");fprintf(fq," 滯止參數(shù) %7.2f %6.2f n",p1,T1);fprintf(fq,"一n");fprintf(fq,"壓氣機(jī)進(jìn)口%7.2f %6.2f %5.2f %8.6f 0 %8.6f%4.2fn",p2,T2,qm3,A2,Lmd2,Ma2);fprintf(fq,"一n");fprintf(fq,"燃燒室進(jìn)口%7.0f %6.2f %5.2f %8.6f 0 %8.6f%4.2fn",p3,T3,qm3

11、a,A3,Lmd3,Ma3);fprintf(fq,"一n");fprintf(fq," 渦導(dǎo)進(jìn)口 %7.0f %6.2f %8.5f %8.6f n",p4,T4,qm4,f);fprintf(fq,"一n");fprintf(fq," 渦導(dǎo)出口 %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

12、.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,&q

13、uot;一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,"

14、排氣速度 %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&qu

15、ot;);if(!fp1)printf("can't open filen");return -1;while(!feof(fp1)fscanf(fp1,"%f",&ci+);fclose(fp1);for(i=0;i<118;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=c13

16、3;Derta=c134;qm3=c135;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=SIG

17、MAi*p1;for(i=0;i<118;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+fn

18、i)+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*Hu-cpg*T4);FILE *fq1=fopen("out2.txt","w");while(!feof(fq1)fprintf(