该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
H=6;
[T0,P0]=ambient(H);
k=1.4;
kg=1.33;
const=-kg/(kg-1);
EFFb=0.98;
EFFm=0.98;
deltaH=43000;%KJ/kg
ccloss=0.03;
exloss=0.215;
P5=P0/(1-exloss);
ERfpt=2;
ERfptacl=1;
ERmap=1.5;
ERcal=1.1;
RPM=50500;
BETA=1;
while(abs(ERfptacl-ERfpt)/ERfpt>0.01)
BETA=BETA+0.01;
[M1,PR,EFF]=compr_map(BETA,RPM);
NDM1=M1;
T1=T0;
P1=P0;
P2=P1*PR;
T2=T1*(1+(PR^((k-1)/k)-1)/EFF);
PWRc=Air_H(T2)-Air_H(T1);
Cp2=Air_Cp(T2);
f=0.019;
while (abs(ERmap-ERcal)/ERmap>0.01)
f=f+0.0001;
T3=T2;
t=1100;
while (abs(T3-t)/t>0.002)