%%Entalpie
i0 = XSteam('h_pT',100, 500);
s0 = XSteam('s_pT', 100, 500);
iz = XSteam('hL_p', 10);
is = XSteam ('hL_p', 2);
iut = XSteam ('h_ps', 10, s0);
ippt = XSteam ('h_ps', 2, s0);
etha_wt = 0.85;
iu = i0 - etha_wt*(i0 - iut);
ipp = i0 - etha_wt*(i0 - ippt);
alfa = (iz - is)/(iu - is);
%%Sprawność termodynamiczna (układ z podgrzewem regeneracyjnym)
etha_or = (alfa*(i0-iu)+(1-alfa)*(i0-ipp))/(i0-iz);
%%Sprawność termodynamiczna (układ bez podgrzewu)
etha = (i0 - ipp)/(i0 - is);
%%Moc na turbinie
m = 100;
PT = m*(alfa*(i0 - iu)+(1-alfa)*(i0-ipp));
%%Moc na wymienniku ciepłowniczym
Qwc = (1-alfa)*m*(ipp-is);
%%Sprawność elektrociepłowni
Qpal = m*(i0-iz);
etha_EC = (PT + Qwc)/Qpal;
%%Wykresy i-s dla różnych ciśnień
for s=1:10
i0(s) = XSteam('h_ps',100, s);
iu(s) = XSteam('h_ps', 10, s);
ipp(s) = XSteam('h_ps', 2, s);
end
s=1:10;
figure
% plot(s,i0,'--b*',s,iu,'-r^',s,ipp,'-.g');
plot(s,i0,s,iu,s,ipp);
grid on
legend ('10MPa', '1MPa', '0.2MPa');
xlabel('Entropia');
ylabel('Entalpia');
Wyszukiwarka
Podobne podstrony:
Projekty EC turbinaProjekty EC projekt19projekt wiazara ec algorytm plyta warstwowaBIGCC projekt ECProjekty EC projektprojekt wiazara ec algorytmProjekty EC projekt12 04Projekty EC projekt26Projekt pracy aparat ortodontyczny ruchomyProjekt mgifprojekt z budownictwa energooszczednego nr 3więcej podobnych podstron