航空发动机原理大作业

南京航空航天大学020810324赵祥成第三章作业第八题:在海平面,静止状态,标准大气条件,最大工作状态时,涡喷发动机,收敛喷管,选择相关参数,计算并画出Fs,SFC及qma随*k(或*3T)的变化曲线.解:选择D组参数:.02.0,03.0Rcol基本计算步骤:(1)进气道出口参数:)211(200*0MkTT*0*11200*0*0*1*)211(PPMkPPTTikk(2)压气机出口参数:*1**1*1*2/)1()(kkkkppkTCTTCW**1*2**1*2*1**1*1*2*/)1(kkkkkkkPPTTTTTTT(3)燃烧室出口气流参数:bPP*2*3(4)单位质量空气的燃油量(油气比):*2*3*2*3arbaahHHhhf(5)涡轮出口气流参数:**3*41**3****3*4'**'')1()1(tkkttttmcolpkptPPTTTTTfCTCT(6)尾喷管出口气流参数:*4*9*4*9TTPPe根据0*9/PP判断燃气在尾喷管中的工作状态,若燃气在尾喷管中为亚临界状态,则燃气在尾喷管中完全膨胀091*40*4'9))(1(2''PPPPTCCkkep若为临界或超临界状态1'*49*4'''9'')21/(12kkekPPTRkkC'29*492pCCTT(7)结论:*999*9'9099'09)1()()1(3600))1(1)(1(TfqAPKqFfSFCPPCTRfCCfFRcolmmascolRcolRcols采用matlab编程(一):在*3T为定值时,Fs,SFC及qma随*k的变化H=0;%地面试车,海平面,静止,标准大气条件,进口参数T0=288.2;C0=0;Ce=340.3;P0=1.0133e+5;ro=1.255;miu=1.789e-5;niucol=0.03;niuR=0.02;M0=0;sigemai=1;%收敛喷管,M2,长度较短k1=1.4;R1=287.3;Cp1=1005;k2=1.33;R2=288;Cp2=1160.7;T0z=T0*(1+(k1-1)/2*M0^2);%进气道出口参数T1z=T0z;P0z=P0*(1+(k1-1)/2*M0^2)^(k1/(k1-1));P1z=P0z*sigemai;yitakz=0.82;%压气机出口参数Pikz=2:0.5:121;%对压气机压比赋值Wk=Cp1*T1z*(Pikz.^((k1-1)/k1)-1)/yitakz;detaTkz=T1z*(Pikz.^((k1-1)/k1)-1)/yitakz;T2z=T1z+detaTkz;P2z=P1z*Pikz;sigmab=0.92;%燃烧室出口气流参数P3z=P2z*sigmab;h2az=Cp1*T2z;%单位质量空气的燃油量(油气比)T3z=1500;%对T3z赋值h3az=Cp1*T3z;kesib=0.98;Hr=4.29e+7;H3z=4034500;%对H3z赋值f=(h3az-h2az)./(kesib*Hr-H3z+h2az);yitam=0.99;%涡轮出口气流参数yitatz=0.9;detaTtz=Cp1*detaTkz./(Cp2*(1-niucol+f)*yitam);T4z=T3z-detaTtz;Pitz=(1-detaTtz./(T3z*yitatz)).^(-k2/(k2-1));P4z=P3z./Pitz;sigemae=0.97;%尾喷管出口气流参数P9z=P4z*sigemae;T9z=T4z;ifP9z/P0=1.85%判断喷管出口状态(亚临界,临界,超临界),确定各参数数值C9=sqrt(2*k2*R2*T4z/(k2+1));T9=T4z-C9.^2/(2*Cp2);P9=sigemae*P4z/((k2+1)/2)^(k2/(k2-1));M9=1;lamda9=1;q_lamda9=1;elseC9=sqrt(2*Cp2*T4z.*(1-(P0./(sigemae*P4z)).^((k2-1)/k2)));T9=T4z-C9.^2/(2*Cp2);P9=P0;M9=C9./(sqrt(k2*R2*T9));lamda9=sqrt((M9).^2*(k2+1)/(2+(k2-1)*(M9).^2));q_lamda9=((k2+1)/2)^(1/(k2-1))*lamda9.*(1-(k2-1)*(lamda9).^2/(k2+1))^(1/(k2-1));endFs=(1-niucol+f+niuR).*(C9+R2*T9.*(1-P0./P9)./C9);%结论SFC=3600*f.*(1-niucol)./Fs;A9=0.02;%对A9赋值Km2=0.03965;qma=Km2*P9z.*A9*q_lamda9./((1-niucol+f+niuR).*sqrt(T9z));subplot(3,1,1)%输出图形plot(Pikz,Fs);xlabel('Pikz');ylabel('Fs(dN*s/kg)');title('在T3z为定值(1500K)时,Fs,SFC及qma随Pikz的变化曲线');subplot(3,1,2)plot(Pikz,SFC);xlabel('Pikz');ylabel('SFC(kg/(dN*h))');subplot(3,1,3)plot(Pikz,qma);xlabel('Pikz');ylabel('qma(kg/s)');02040608010012014005001000PikzFs(dN*s/kg)在T3z为定值(1500K)时,Fs,SFC及qma随Pikz的变化曲线0204060801001201400.050.10.150.20.25PikzSFC(kg/(dN*h))020406080100120140246810Pikzqma(kg/s)(二):在*k为定值时,Fs,SFC及qma随*3T的变化对(一)的代码做如下修改:Pikz=25;%对压气机压比赋值h2az=Cp1*T2z;%单位质量空气的燃油量(油气比)T3z=1000:1:2000;%对3z赋值h3az=Cp1*T3z;kesib=0.98;Hr=4.29e+7;H3z=2311200:3656.9:5968100;%对H3z赋值f=(h3az-h2az)./(kesib*Hr-H3z+h2az);subplot(3,1,1)%输出图形plot(T3z,Fs);xlabel('T3z(K)');ylabel('Fs(dN*s/kg)');title('在Pikz为定值(25)时,Fs,SFC及qma随T3z的变化曲线');subplot(3,1,2)plot(T3z,SFC);xlabel('T3z(K)');ylabel('SFC(kg/(dN*h))');subplot(3,1,3)plot(T3z,qma);xlabel('T3z(K)');ylabel('qma(kg/s)');10001100120013001400150016001700180019002000050010001500T3z(K)Fs(dN*s/kg)在Pikz为定值(25)时,Fs,SFC及qma随T3z的变化曲线100011001200130014001500160017001800190020000.060.080.10.12T3z(K)SFC(kg/(dN*h))10001100120013001400150016001700180019002000051015T3z(K)qma(kg/s)(三):*k与*3T都为变量时,Fs,SFC,F的变化以及取点采用hold语句:先对*k取5,10,15,20,25,对*3T取1200K~1800K,步长取为1,得出(Fs,SFC)五条图线,再对*3T取1200K,1400K,1600K,1800K,对*k取为5~25,步长取为0.02,得出(Fs,SFC)四条图线,将这九条图线用hold语句绘制在一幅坐标图上:500600700800900100011000.060.070.080.090.10.110.120.130.140.15同理得出(F,SFC)图线,取点如下: