GasTurbine

T1=input('Ingrese la temperatura de entrada en la primera etapa de compresión en kelvin: ');
P1=input('Ingrese la presión de entrada en la primera etapa de compresión en kPa: ');
T3=input('Ingrese la temperatura de entrada en la segunda etapa de compresión en kelvin: ');
P2=input('Ingrese la presión de entrada en la segunda etapa de compresión en kPa: ');
rc=input('Ingrese la relación decompresiones extremas: ');
T6=input('Ingrese la temperatura de entrada en la primera etapa de la turbina en kelvin: ');
T8=input('Ingrese la temperatura de entrada en la segunda etapa de la turbina en kelvin: ');
P7=input('Ingrese la presión de entrada en la segunda etapa de expansión en kPa: ');
re=input('Ingrese la relación de expansiones extremas: ');
AC=input('Ingrese el flujo volumetricodel compresor m^3/s: ');
clc
P3=P2;
P8=P7;


%Solicitud de eficiencia de la turbina de gas para cálculo de ciclo Brayton

efi=menu('¿EL ciclo de Turbina de Gas es Ideal en cuanto a la turbina y el compresor?','No','Si');
switch efi
case 1

%Eficiencia para la turbina de gas de dos etapas
muc1=input('Ingrese la eficiencia del compresor de la primera etapa: ');muc1=muc1/100;
mut1=input('Ingrese la eficiencia de la turbina de la primera etapa: ');
mut1=mut1/100;
muc2=input('Ingrese la eficiencia del compresor de la segunda etapa: ');
muc2=muc2/100;
mut2=input('Ingrese la eficiencia de la turbina de la segunda etapa: ');
mut2=mut2/100;
h1=interp1(T,h,T1);
Pr1=interp1(T,Pr,T1);Pr2=(P2/P1)*Pr1;
T2=interp1(Pr,T,Pr2);
h2=interp1(Pr,h,Pr2);
h2a=h1+((h2-h1)/muc1);
T2a=interp1(h,T,h2a);
h3=interp1(T,h,T3);
Pr3=interp1(T,Pr,T3);
Pr4=((P1*rc)/P3)*Pr3;
T4=interp1(Pr,T,Pr4);
h4=interp1(Pr,h,Pr4);
h4a=h3+((h4-h3)/muc2);
T4a=interp1(h,T,h4a);
h6=interp1(T,h,T6);Pr6=interp1(T,Pr,T6);
Pr7=Pr6*(P7/(re*P1));
T7=interp1(Pr,T,Pr7);
h7=interp1(Pr,h,Pr7);
h7a=h6-((h6-h7)*mut1);
T7a=interp1(h,T,h7a);
h8=interp1(T,h,T8);
Pr8=interp1(T,Pr,T8);
Pr9=Pr8*(P1/P8);
T9=interp1(Pr,T,Pr9);
h9=interp1(Pr,h,Pr9);
h9a=h8-((h8-h9)*mut2);
T9a=interp1(h,T,h9a);%Análisis final para el cálculo de potencia del ciclo Brayton

Reg=menu('¿La turbina tiene regenerador? En que caso de ser cierta su respuesta indique la eficiencia','Si','No');
switch Reg
case 1
mur=input('Eficiencia del regenerador: ');
mur=mur/100;
h5=h4a+mur*(h9a-h4a);
T5=interp1(h,T,h5);
h10a=h4a+h9a-h5;T10a=interp1(h,T,h10a);
Wc=(h2a-h1)+(h4a-h3);
Wt=(h6-h7a)+(h8-h9a);
qe=(h6-h5)+(h8-h7a);
qs=h9a-h1;
mu=((Wt-Wc)/qe)*100;
rw=(Wc/Wt)*100;
m=((AC*P1)/((8314/28.97)*(T1)))*1000;
wciclo=Wt-Wc;
Wciclo=m*(Wt-Wc);
clc
disp('Solución de Análisis aire-estandar')disp('__________________________________')
disp('Parámetro Análisis aire-estandar')
disp('__________________________________')
fprintf(' T2a %6.3f k \n',T2a)
fprintf(' T4a %6.3f k \n',T4a)
fprintf(' T7a %6.3f k \n',T7a)
fprintf(' T9a %6.3f k \n',T9a)
fprintf(' T5 %6.3f k \n',T5)
fprintf(' T10a %6.3f k\n',T10a)
fprintf(' Wc %6.3f kJ/kg \n',wciclo)
fprintf(' rw %6.3f %%\n',rw)
fprintf(' m %6.3f kg/s\n',m)
fprintf('Eficiencia %6.3f %%\n',mu)
fprintf(' Potencia %4.3f kW\n',Wciclo)
T=[T1,T2,T3,T4,T5,T6,T7,T8,T9,T1];
S=[2,2,1,1,3.5,5,5,6,6,2];
Ta=[T1,T2a,T3,T4a,T5,T6,T7a,T8,T9a,T10a,T1];...