Листинг 6. Расчет внутренних КПД цикла Ренкина
для различных начальных температур
Program CIKL_RENK; Uses Crt;
Const Rw=0.46151; ym1=0.21; ym2=0.27315; K1=1e3; Tk0=273.0;
p2=0.004; tvs=30; {Temper. pri p2} t12=600;{pred. temp.} kh0=3.39e-8;
p1=13.0;{davL.nas.} t1=400;{temp.nas.} kp0T=0.85; kpH=0.87;
bp0=3.237e-4; bp1=2.5e-4; bp2=–1.1354e-3; bp3=-4.381e-4;
cp0=5.6084e-6; cp1=-2.5994e-6; cp2=–1.2604e-8;
kp0=2.12787e3; kp1=1.48285e3; kp2=3.79026e2; kp3=4.6174e1;
kp4=1.081161e4; as0=–7.821541; as1=8.286568e1;
as2=10.28003;as3=–11.48776; bs0=2.20732; bs1=–2.117e-1;
bs2= –2.16605e-3; bs3=1.619692e-4; bs4=4.89980e-5; bs5=3.691725e-6;
c10= –3.15399e3; c11=2.913765e4; c12=1.224973e5; c13=2.984568e5;
c14=–3.632168e5; c15=1.785296e5; c20=6.0102277e3; c21=–4.7493e4;
c22=2.388416e5; c23=-5.704046e5; c24=6.772865e5; c25=–3.264862e5;
d10=2.71288e-3; d11=–2.51341e-2; d12=1.590277e-1; d13=–5.625152e-1;
d14=1.16296; d15=–1.299799; d16=0.6110896; d20=0.9997; d21=–0.029;
d22=–0.2; d23=–10.0; d24=-4.4e12; e10=–11.54816; e11=96.15764;
e12=–341.8428; e13=7.197764e2; e14=-7.973969e2; e15=3.640519e2;
e20=2.960815e1; e21=–1.327532e2; e22=1.68014e2; e23=6.151844e2;
e24=–2.409461e3; e25=3.125479e3; e26=–1.470736e3; an0=9.771e-4;
an1=1.774e-4; an2=2.52e-5; an3=2.96e-6; bv0=3.225e-6; bv1=1.3436e-6;
bv2=1.684e-8; bv3=1.432e-7; cv0=3.7e-8; cv1=3.558e-8; cv2=-4.05e-13;
eh0=49.4; eh1=4.225e2; eh2=4.767; eh3=3.333e-2; fh0=–9.25; fh1=1.67;
fh2=7.36e-3; fh3=–8e-3; gh0=–7.3e-2; gh1=7.9e-2; gh2=6.8e-4;
Var Tk1, Tk2, y, y2, y3, y6, y8, y9, y14, y15, y21, y31, ys, ys2, ps, Ly,
Ls, Sbs, Ap0, Ap1, Ap2, Ap3, Ap4, Ap5, hp0, sp0, vpp, hpp, spp, ts, Tks,
v21, v22, h21, h22, s21, s22, Kh, Gh, Fh, Eh, tn1, pn1, tn2, tn3, tn6, tn12,
Ltn, hw,x2, h2x, H0, x2d, h2d, h3, h3d, LtH, LdH, kpt, kpi, kpd, dti, did:real;
Procedure TDSPP(t,p:real); {Peregr par}
Begin y:=Tk1/K1; y2:=sqr(y); y21:=sqr(y–ym1); Ly:=ln(y);
y3:=exp(3*Ly); y31:=exp(3*ln(y-ym1)); y6:=exp(6*Ly);
y8:=exp(8*Ly);y9:=exp(9*Ly); y14:=exp(14*Ly);y15:=exp(15*Ly);
Ap0:=bp0+bp1*y+bp2/y2+bp3/y21; Ap1:=cp0+cp1/y8+cp2/y14;
Ap2:=bp0+3*bp2/y2+3*bp3/y21+2*bp3*ym1*y31;
Ap3:=cp0+9*cp1/y8+15*cp2/y14; Ap4:=–bp1+2*bp2/y3+2*bp3/y31;
Ap5:=8*cp1/y9+14*cp2/y15; vpp:=Ap0+Ap1*p+Rw*y/p;
hp0:=kp0+kp1*y+kp2*y2+kp3*Ly; hpp:=hp0+(Ap2+Ap3*p/2)*K1;
sp0:=(kp1*Ly+2*kp2*y-kp3/y+kp4)/K1;
spp:=sp0+(Ap4+Ap5*p/2)*p–Rw*ln(K1*p); end;
Procedure TDSPH(ts:real);{Nas par}