O Programa Para HP Prime

Programa Para HP Prime

EXPORT Turbina_a_Gás()

BEGIN

LOCAL T1,T1a,T2,T2a,T2real,Txc,Txd;

LOCAL T3,T3a,T4,T4a,T4real,Tff,Tfq;

LOCAL ηcp,ηcp1,ηt,ηt1,ηg,ηg1,ηi,ηi1;

LOCAL ηcarnot,ηcarnot1;

LOCAL MENU,menu1,menu2;

LOCAL Wcp,Wt,Wele,Wg,Wtr,Wgr;

LOCAL Mar,Cp,Mc,TX,c,PCI;

LOCAL ηcarnot,ηcc1,ηcc,h;

REPEAT

CHOOSE(MENU,"Turbina a Gás","Compressor","Turbina","Gerador Ideal","Potencia Eletrica","Camara de Combustão","Potencia Corrigida da Turbina","Potencia Eletrica Corrigida","Rendimento da Turbina a Gás (Instalação)","Consumo Especifico","Rendimento de Carnot","Sair");

IF MENU==1 THEN

REPEAT

CHOOSE(menu1,"Compressor","Temperatura ideal","Rendimento","Potencia","Voltar");

IF menu1==1 THEN

INPUT({Txc,T1a,T2a},"Temperatura ideal",{"P2/P1","Temperatura de entrada T1 (ºC)","Temperatura de saida T2 (ºC)"});

IF Txc==0 THEN

T2:=T2a+273;

T1:=T1a+273;

Txc:=(T2/T1)^(7/2);

MSGBOX("P2/P1=(T2/T1)^(7/2)");

MSGBOX("P2/P1="+Txc);

END;

IF T1a==0 THEN

T2:=T2a+273;

T1:=T2/Txc^(2/7);

T1a:=T1-273;

MSGBOX("T1=T2/(P2/P1)^(2/7)");

MSGBOX("T1 (K)="+T1);

MSGBOX("T1 (ºC)="+T1a);

END;

IF T2a==0 THEN

T1:=T1a+273;

T2:=Txc^(2/7)*T1;

T2a:=T2-273;

MSGBOX("T2=(P2/P1)^(2/7)*T1");

MSGBOX("T2 (K)="+T2);

MSGBOX("T2 (ºC)="+T2a);

END;

END;

IF menu1==2 THEN

INPUT({ηcp1,T1a,T2a,T2real},"Rendimento",{"Rendimento do Compressor (%)","Temperatura de entrada T1 (ºC)","Temperatura de saida T2 (ºC)","Temperatura Real (T2real)(ºC)"});

IF ηcp1==0 THEN

ηcp:=(T2a-T1a)/(T2real-T1a);

ηcp1:=ηcp*100;

MSGBOX("ηcp=((T2-T1)/(T2real-T1))*100");

MSGBOX("Rendimento do Compressor (%)="+ηcp1);

END;

IF T1a==0 THEN

ηcp:=ηcp1/100;

T1a:=((T2real*ηcp)-T2a)/(ηcp-1);

MSGBOX("T1=((T2real*ηcp)-T2)/(1-ηcp)");

MSGBOX("T1 (ºC)="+T1a);

END;

IF T2a==0 THEN

ηcp:=ηcp1/100;

T2a:=T1a+(T2real-T1a)*ηcp;

MSGBOX("T2=T1+(T2real-T1)*ηcp");

MSGBOX("Temperatura teorica (T2)(ºC)="+T2a);

END;

IF T2real==0 THEN

ηcp:=ηcp1/100;

T2real:=T1a+(T2a-T1a)/ηcp;

MSGBOX("T2real:=(T2-T1+T1*ηcp)/ηcp");

MSGBOX("Temperatura Real (T2)(ºC)="+T2real);

END;

END;

IF menu1==3 THEN

INPUT({Txc,T1a,Mar,Wcp,ηcp1},"Potencia teorica do Compressor",{"Taxa de compressão (P2/P1)","Temperatura de entrada T1 (ºC)","Massa de Ar (Kg/h)","Potencia do Compressor (kW)","Rendimento do Compressor (%)"});

IF Wcp==0 THEN

T1:=T1a+273;

Cp:=0.24;

ηcp:=ηcp1/100;

Wcp:=Mar*Cp*T1*(Txc^(2/7)-1)*(1)/(860*ηcp);

MSGBOX("Wcp=Mar*Cp*T1*(Txc^(2/7)-1)*(1)/(860*ηcp)");

MSGBOX("Wcp (kW)="+Wcp);

END;

IF Mar==0 THEN

T1:=T1a+273;

Cp:=0.24;

ηcp:=ηcp1/100;

Mar:=(Wcp*ηcp*860)/(Cp*T1*(Txc^(2/7)-1));

MSGBOX("Mar=(Wcp*ηcp*860)/(Cp*T1*(Txc^(2/7)-1))");

MSGBOX("Mar (Kg/h)="+Mar);

END;

IF T1a==0 THEN

Cp:=0.24;

ηcp:=ηcp1/100;

T1:=(Wcp*ηcp*860)/(Cp*Mar*(Txc^(2/7)-1));