蒸汽特性表.doc

this program used to excel for sat/wat/ste region ************************************************* p=ps(t) sature t=ts(p) sature hw=wat(t,p,1);sw=wat(t,p,2);vw=wat(t,p,3) water hs=ste(t,p,1);ss=ste(t,p,2);vs=ste(t,p,3) steam t_ste=pht2(p,h);t_ste=pst2(p,s) steam with BTN p_ste=thp2(t,h);p_ste=tsp2(t,s) steam with BTN t_ste=pht(p,h);t_ste=pst(p,s) steam with IFC p_ste=thp(t,h);p_ste=tsp(t,s) steam with IFC hht=hht(p0,t0,p1) steam with IFC hhh=hhh(p0,h0,p1) steam with IFC itai=itai(p0,t0,p1,t1) steam with IFC hsp=hspt(h,s,1);hst=hspt(h,s,2) steam with IFC ************************************************* created by 张传明 2004.07.30 *** function 1 sat_r *** Function ps(t) Dim ot, opk, o1, o11, p11, p12, p1, p2 As Double Dim k(1 To 9) As Double k(1) = -7.691234564 k(2) = -26 k(3) = -168.1706546 k(4) = 64 k(5) = -118.9646225 k(6) = 4 k(7) = 20.9750676 k(8) = 1000000000# k(9) = 6# ot = (t + 273.15) / 647.3 o1 = 1# - ot o11 = o1 * o1 p11 = k(1) + k(2) * o1 + o11 * (k(3) + k(4) * o1 + k(5) * o11) p12 = ot * (1 + k(6) * o1 + k(7) * o11) p1 = p11 / p12 p2 = 1 / (k(8) * o11 + k(9)) p12 = p1 - p2 If p12 = 0 Then opk = Exp(o1 * (p12)) Else opk = 1 / Exp(o1 * Abs(p12)) End If ps = opk * 22.12 End Function *** function 2 sat_r *** Function ts1(p) Dim y, y1 As Double Dim i As Integer Dim b(5) As Double b(0) = 2.20732 b(1) = -0.2117187 b(2) = -0.002166605 b(3) = 0.0001619692 b(4) = 0.000048996 b(5) = 0.000003691725 y = Log(p) y1 = b(5) For i = 1 To 5 y1 = y1 * y + b(5 - i) Next i ts1 = 1# / y1 * 1000# - 273.15 End Function Function ts(p) Dim t1, t2, t3, p1, p2, p3, a, b, c, d As Double Dim i As Integer again_k: t1 = ts1(p) p1 = ps(t1) t2 = ts1(p1) p2 = ps(t2) t3 = ts1(p2) p3 = ps(t3) d = (p3 * p3 - p1 * p1) / (p2 * p2 -