%% Dane
P1 = 120;
P2 = 5;
etha_w = 0.8;
Qc = 230;
u = 0.5;
etha_k = 0.75;
Qw = 20;
etha_mech = 0.98;
etha_gen = 0.99;
%% Czynniki robocze
i1 = XSteam('h_pT',P1, 530);
s1 = XSteam('s_pT', P1, 530);
i2t = XSteam('h_ps', P2, s1);
i2 = i1-etha_w*(i1-i2t);
s2 = XSteam('s_ph', P2, i2);
i3 = XSteam('hL_p', P2);
s4 = XSteam('s_ph', P2, i1);
D = Qc/(1.25*(i2-i3));
Ds = 0.25*D;
G = (Ds*(i2-i1))/(i3-i1);
DR = Ds-G;
Dk = DR+D;
Pg = D*(i1-i2)*etha_mech*etha_gen;
%% Sprawność termodynamiczna
etha_td = (Pg+u*Qc)/(Dk*(i1-i3));
B = (Dk*(i1-i3))/(etha_k*Qw);
%% Wykresy
for s=1:10
i1x(s) = XSteam('h_ps', P1, s);
i2x(s) = XSteam('h_ps', P2, s);
end
s=1:10;
figure
plot(s,i1x,s,i2x);
grid on
ylim ([1500 4000]);
xlim ([4 9]);
legend ('12 MPa', '0.5 MPa');
xlabel('Entropia');
ylabel('Entalpia');
hold on
sx=[s1 s1 s1 s2 s2 s4 s4 s1];
ix=[i2t i1 i1 i2 i2 i1 i1 i1];
text(7, 3500, 'redukcja');
plot(sx, ix, '-r*', 'LineWidth',2);
hold off