tbc, Tbk, Ns1, Ns2, Ns3, Ns4, Ns5, qz1, qz2, qz3, qz4, qzz, qx1, qx2, qc,
qxz, Gc, cv, cp, cn, Gw, Gp, qp, rp, gwp, gwk, t2, ps1, ps2, pp2, TT1,
Vbc1, xs, Vbc2, Vbc3, Vbc4, Vbc5, rp1, rp4, rp5, rp5k, hs1, hs2, qxzs,
ys, ts, Nss, rp0: real;
Function Pns(ts:real):real; Begin
ys:=ts/kj; pns:=exp(a0s+a1s/ys+a2s*ys+a3s*ln(ts));
end;
Function hns1(ts:real):real; Begin
ys:=ts/kj; hns1:=c0+ys*(c1+ys*(c2+ys*(c3+ys*(c4+c5*ys))));
end;
Function hns2(ts:real):real; Begin
ys:=ts/kj; hns2:=f0+ys*(f1+ys*(f2+ys*(f3+ys*(f4+f5*ys))));
end;
Procedure Krp(p,t:real; var rp:real);
Begin
p:=p/kd; if p<300 then rp0:=kr0+p*(kr1+kr2*p)
else rp0:=kr3+kr4*p; rp:=1.03*(1+t/273)*rp0; end;
Begin ClrScr;
en:=(n–1)/n; en1:=(n1–1)/n1;
сv:=Ru/(mm*(k–1)); cp:=cv*k;
xs:=exp((1/z)*ln(pka/p1a)); kps:=exp(en*ln(xs));
tbc:=ts1+dtw+dtb; Tbk:=tbc+273; Tk1:=t1+273; TTk1:=Tbk/Tk1;
p1k:=p1a*xs; ps1:=pns(Tk1); ps2:=pns(Tbk);
hs1:=hns1(Tk1); hs2:=hns2(Tk1); rp:=hs2–hs1;
pp2:=fib*ps1*p1k/p1a; p2k:=p1k*xs; p3k:=p2k*xs; p4k:=p3k*xs;
p5k:=p4k*xs; T1k:=Tk1*kps; T2k:=Tbk*kps; T3k:=T2k; T4k:=T2k;
T5k:=T2k; t5c:=T5k–273;
Krp(p1a,t1,rp1); Krp(p4k,tbc,rp4);
Krp(p5k,tbc,rp5); Krp(p5k,t5c,rp5k);
Vbc1:=V1*(p1a–fib*ps1)/p1a;
Vbc2:=Vbc1*TT1*p1a/(p1k–ps2);
Vbc3:=Vbc1*TT1*p1a/(p2k–ps2);
Vbc4:=Vbc1*(TT1)*p1a/p3k;
Vbc5:=Vbc1*TT1*(rp4/rp1)*p1a/p4k;
Lv1:=1–e01*(exp(ln(xs)/n1)–1);
Lv2:=1–e02*(exp(ln(xs)/n2)–1);
Lv3:=1–e03*(exp(ln(xs)/n3)–1);
Lv4:=1–e04*(exp(ln(xs)/n4)–1);
Lv5:=1–e05*((rp4/rp5k)*xs–1);
Vh1:=V1/(Lv1*2*n0);
Vh2:=Vbc2/(Lv2*2*n0);