%%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');