plik


%%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 turbina
Projekty EC projekt19
projekt wiazara ec algorytm plyta warstwowa
BIGCC projekt EC
Projekty EC projekt
projekt wiazara ec algorytm
Projekty EC projekt12 04
Projekty EC projekt26
Projekt pracy aparat ortodontyczny ruchomy
Projekt mgif
projekt z budownictwa energooszczednego nr 3

więcej podobnych podstron