Projeto Conceitual de Aeronaves de Transporte

O código MATLAB que foi escrito para gerar os gráficos das Figs. I.5 e I.6 é dado a
seguir:
%*********************************************************************
% Estimativa do consumo de combustivel em cruzeiro do Fokker 100
%*********************************************************************
clc
% constantes
kg2lb = 2.20462262; % fator de conversao coversao de kg para libra
ft2m = 0.3048;% fator de conversao de p´s para metros
cfe = 0.0030; % coef de fricao equivalente (Dado do Roskam)
% Atmosfera
rho35 = 0.37959; % Densidade do ar em 35000 pes (kg/m3)
vsom35= 296.5388; % Velocidade do som em 35000 pes (m/s)
rho31 = 0.44165; % Densidade do ar em 31000 pes (kg/m3)
vsom31= 301.86; % Velocidade do som em 31000 pes (m/s)
rho30 = 0.45831; % Densidade do ar em 30000 pes (kg/m3)
vsom30= 303.177; % Velocidade do som em 30000 pes (m/s)
% Dados da aeronave
arw= 8.43; % Alongamento da asa
sw=93.5; % Area da asa
mtow=43090; % Masssa maxima de decolagem
ftakeoff=0.995;
fclimb=0.98;
rangekm=1730; % (km)
ct=0.69; % consumo especifico do motor em cruzeiro (1/h)
%
mtowlb=mtow*kg2lb;
rangem=rangekm*1000;
swft2=sw/(ft2m*ft2m);
areamolhadaft2=0.0199+0.7531*log10(mtowlb);
areamolhadaft2=10^areamolhadaft2;
cd0=cfe*areamolhadaft2/swft2;
oswa = 0.65;
ncont=1;
while oswa < 0.85
e=oswa;
oswplo(ncont)=oswa;
cl_long=sqrt(pi*e*arw*cd0/3);
cd =4*cd0/3;
ld=cl_long/cd;
% Velocidade de cruzeiro de longo alcance no inicio do cruz
masscruzi=ftakeoff*fclimb*mtow; % (kg)
vcruz35i=sqrt(masscruzi*9.81/(0.50*rho35*cl_long*sw));
vcruz31i=sqrt(masscruzi*9.81/(0.50*rho31*cl_long*sw));
vcruz30i=sqrt(masscruzi*9.81/(0.50*rho30*cl_long*sw));
% fracao de massa no cruzeiro
t1=sqrt(2/(rho30*sw));
t2=(sqrt(cl_long))/cd;
t3=3600/ct;
sqqwf=-rangem/(2*t1*t3*t2)+sqrt(masscruzi*9.81);
massfin=sqqwf^2/9.81;
fmassa=massfin/masscruzi;
mcomb=masscruzi-massfin;
mcomb30plo(ncont)=mcomb;
mach30plo(ncont)=vcruz30i/vsom30;
teralt=rangem*ct/(3600*vcruz30i*ld);
fmassalt=1/exp(teralt);
%