蒸汽特性表.docx
《蒸汽特性表.docx》由会员分享，可在线阅读，更多相关《蒸汽特性表.docx（33页珍藏版）》请在冰豆网上搜索。
蒸汽特性表
'thisprogramusedtoexcelforsat/wat/steregion
'*************************************************
'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）>steamwithBTN
'p_ste=thp2（t,h）;p_ste=tsp2（t,s）>steamwithBTN
't_ste=pht（p,h）;t_ste=pst（p,s）>steamwithIFC
'p_ste=thp（t,h）;p_ste=tsp（t,s）>steamwithIFC
'hht=hht（p0,t0,p1）>steamwithIFC
'hhh=hhh（p0,h0,p1）>steamwithIFC
'itai=itai（p0,t0,p1,t1）>steamwithIFC
'hsp=hspt（h,s,1）;hst=hspt（h,s,2）>steamwithIFC
'*************************************************
'createdby张传明2004.07.30
'***function1sat_r***
Functionps（t）
Dimot,opk,o1,o11,p11,p12,p1,p2AsDouble
Dimk（1To9）AsDouble
k
（1）=-7.691234564
k
（2）=-26.08023696
k（3）=-168.1706546
k（4）=64.23285504
k（5）=-118.9646225
k（6）=4.16711732
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
Ifp12>=0Then
opk=Exp（o1*（p12））
Else
opk=1/Exp（o1*Abs（p12））
EndIf
ps=opk*22.12
EndFunction
'***function2sat_r***
Functionts1（p）
Dimy,y1AsDouble
DimiAsInteger
Dimb（5）AsDouble
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）
Fori=1To5
y1=y1*y+b（5-i）
Nexti
ts1=1#/y1*1000#-273.15
EndFunction
Functionts（p）
Dimt1,t2,t3,p1,p2,p3,a,b,c,dAsDouble
DimiAsInteger
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-p1*p1）
b=（d*（t2-t1）-（t3-t1））/（d*（p2-p1）-（p3-p1））
a=（（t2-t1）-b*（p2-p1））/（p2*p2-p1*p1）
c=t1-a*p1*p1-b*p1
t4=a*p*p+b*p+c
'IfAbs（t1-t4）/t4>0.0001Then
't1=t4
'GoToagain_k
'EndIf
ts=t4
EndFunction
'***function3wat_r***
Functionwat（t,p,para）
DimiAsInteger
Dimot1,ot2,ot7,ot10,ot11,ot17,ot18,ot19,ot20,ot21AsDouble
Dimxt,xx,xx1,y,y1,z,oh,oh0,os,ovAsDouble
Dimot,op,oh1,os1,ov1AsDouble
Consta0=-8.4820150677145E-06,a1=-1.6551695754548E-06
Dima（0To22）,aa（1To12）AsDouble
a（0）=6824.687714
a
（1）=-542.2063673
a
（2）=-20966.66205
a（3）=39412.86787
a（4）=-67332.77739
a（5）=99023.81028
a（6）=-109391.1774
a（7）=85908.41667
a（8）=-45111.68742
a（9）=14181.38926
a（10）=-2017.271113
a（11）=7.982692717
a（12）=-0.02616571843
a（13）=0.00152241179
.022********
a（15）=242.1647003
a（16）=1.269716088E-10
a（17）=2.074838328E-07
a（18）=2.17402035E-08
a（19）=1.105710498E-09
a（20）=12.93441934
a（21）=0.00001308119072
a（22）=6.047626338E-14
aa
（1）=0.8438375405
aa
（2）=0.0005362162162
aa（3）=1.72
.0734*******
aa（5）=0.0497585887
aa（6）=0.65371543
aa（7）=0.00000115
aa（8）=0.000015108
aa（9）=0.14188
aa（10）=7.002753165
aa（11）=0.0002995284926
aa（12）=0.204
ot=（t+273.15）/647.3
op=p/22.12
SelectCasepara
Case1
ot1=Log（ot）
ot2=ot*ot
ot7=Exp（7*ot1）
ot11=Exp（11*ot1）
ot18=ot11*ot7
ot19=ot18*ot
ot20=ot18*ot2
y=1#-（aa
（1）*ot+aa
（2）/ot7）*ot
y1=-2#*aa
（1）*ot+6#*aa
（2）/ot7
z=y+Sqr（aa（3）*y*y+2#*（-aa（4）*ot+aa（5）*op））
oh=a0+a（0）*ot
xt=1/ot
Fori=1To10
xt=xt*ot
oh=oh-（i-2）*a（i）*xt
Nexti
oh0=z*（17*（z/29-y/12）+5*ot*y1/12）
oh=oh+a（11）*（oh0+（aa（4）-（aa（3）-1）*y*y1）*ot）*Exp（（-5/17）*Log（z））
xx=aa（6）-ot
xx1=xx*xx*xx*xx
xx=xx1*xx1*xx
oh0=a（12）-a（14）*ot2+a（15）*（9*ot+aa（6））*xx
oh=oh+（oh0+a（16）*（20*ot19+aa（7））/（aa（7）+ot19）/（aa（7）+ot19））*op
oh=oh-（12*ot11+aa（8））*（a（17）+（a（18）+a（19）*op）*op）*op/（aa（8）+ot11）/（aa（8）+ot11）
oh=oh+a（20）*（17*aa（9）*ot18+19*ot20）*（Exp（（-3）*Log（aa（10）+op））+aa（11）*op）
oh1=oh+（a（21）*aa（12）+21*a（22）*op/ot20）*op*op*op
wat=oh1*70.1204
Case2
ot1=Log（ot）
ot7=Exp（7*ot1）
ot10=Exp（10*ot1）
ot11=ot10*ot
ot17=ot7*ot10
ot18=ot17*ot
ot19=ot18*ot
ot21=ot19*ot*ot
y=1-（aa
（1）*ot+aa
（2）/ot7）*ot
y1=-2*aa
（1）*ot+6*aa
（2）/ot7
z=y+Sqr（aa（3）*y*y+2#*（-aa（4）*ot+aa（5）*op））
os=-a1+a（0）*ot1
xt=1/ot
Fori=2To10
xt=xt*ot
os=os-（i-1）*a（i）*xt
Nexti
os=os+a（11）*（（（5/12）*z-（aa（3）-1）*y）*y1+aa（4））*Exp（（-5/17）*Log（z））
xx=aa（6）-ot
xx1=xx*xx*xx*xx
xx=xx1*xx1*xx
os=os+（-a（13）-2*a（14）*ot+10*a（15）*xx+19*a（16）*ot18/（aa（7）+ot19）/（aa（7）+ot19））*op
os=os-11*ot10*（a（17）+（a（18）+a（19）*op）*op）*op/（aa（8）+ot11）/（aa（8）+ot11）
os=os+a（20）*（18*aa（9）*ot17+20*ot19）*（Exp（-3*Log（aa（10）+op））+aa（11）*op）
os1=os+（a（21）+20*a（22）*op/ot21）*op*op*op
wat=os1*0.1083275143
Case3
ot1=Log（ot）
ot7=Exp（7*ot1）
ot11=Exp（11*ot1）
ot18=ot7*ot11
ot19=ot18*ot
ot20=ot19*ot
y=1#-（aa
（1）*ot+aa
（2）/ot7）*ot
z=y+Sqr（aa（3）*y*y+2#*（-aa（4）*ot+aa（5）*op））
ov=a（11）*aa（5）*Exp（（-5/17）*Log（z））
xx=（aa（6）-ot）*（aa（6）-ot）
xx=xx*xx*xx*xx*xx
ov=ov+a（12）+（a（13）+a（14）*ot）*ot+a（15）*xx+a（16）/（aa（7）+ot19）
ov=ov-（a（17）+（2#*a（18）+3#*a（19）*op）*op）/（aa（8）+ot11）
ov=ov-a（20）*（ot18*aa（9）+ot20）*（-3#*Exp（-4*Log（aa（10）+op））+aa（11））
ov1=ov+（3#*a（21）*（aa（12）-ot）+4#*a（22）*op/ot20）*op*op
wat=ov1*0.00317
CaseElse
wat=0
EndSelect
EndFunction
'***function4usedforste_r***
Functionolp（ot）
Constl0=15.74373327,l1=-34.17061978,l2=19.31380707
olp=l0+l1*ot+l2*ot*ot
EndFunction
'***function5usedforste_r***
Functionolp1（ot）
Constl1=-34.17061978,l2=19.31380707
olp1=l1+2*l2*ot
EndFunction
'***function6ste_r***
Functionste（t,p,para）
Dimi,jAsInteger
Dimx1,op0,op1,op2,oh,oh0,oh1,oh3,oh4,oh5,ot1,oppAsDouble
Dimos,os0,os1,os3,os4,os5AsDouble
Dimov,ov1,ov0,ov3AsDouble
Dimot,op,oh2,os2,ov2AsDouble
Consta0=-8.4820150677145E-06,a1=-1.6551695754548E-06
Consti1=4.260321148,bb0=0.7633333333
Dimn（1To8）AsInteger
n
（1）=2
n
（2）=3
n（3）=2
n（4）=2
n（5）=3
n（6）=2
n（7）=2
n（8）=2
Diml（6To8）AsInteger
l（6）=1
l（7）=1
l（8）=2
Dimz（1To8,1To3）AsInteger
z（1,1）=13
z（1,2）=3
z（1,3）=0
z（2,1）=18
z（2,2）=2
z（2,3）=1
z（3,1）=18
z（3,2）=10
z（3,3）=1
z（4,1）=25
z（4,2）=14
z（4,3）=0
z（5,1）=32
z（5,2）=28
z（5,3）=24
z（6,1）=12
z（6,2）=11
z（6,3）=0
z（7,1）=24
z（7,2）=18
z（7,3）=0
z（8,1）=24
z（8,2）=14
z（8,3）=0
Dimx（6To8,1To2）AsInteger
x（6,1）=14
x（6,2）=0
x（7,1）=19
x（7,2）=0
x（8,1）=54
x（8,2）=27
Dimb（0To9,0To6）AsDouble
b（0,0）=16.83599274
b（0,1）=28.56067796
b（0,2）=-54.38923329
b（0,3）=0.4330662834
b（0,4）=-0.6547711697
.0856*******
b（0,6）=0
b（1,0）=0
b（1,1）=0.06670375918
b（1,2）=1.388983801
b（1,3）=0
b（1,4）=0
b（1,5）=0
b（1,6）=0
b（2,0）=0
.0839*******
b（2,2）=0.02614670893
b（2,3）=-0.03373439453
b（2,4）=0
b（2,5）=0
b（2,6）=0
b（3,0）=0
b（3,1）=0.4520918904
b（3,2）=0.1069036614
b（3,3）=0
b（3,4）=0
b（3,5）=0
b（3,6）=0
b（4,0）=0
b（4,1）=-0.5975336707
b（4,2）=-0.08847535804
b（4,3）=0
b（4,4）=0
b（4,5）=0
b（4,6）=0
b（5,0）=0
b（5,1）=0.5958051609
b（5,2）=-0.5159303373
b（5,3）=0.2075021122
b（5,4）=0
b（5,5）=0
b（5,6）=0
b（6,0）=0
b（6,1）=0.1190610271
b（6,2）=-0.09867174132
b（6,3）=0
b（6,4）=0
b（6,5）=0
b（6,6）=0
b（7,0）=0
b（7,1）=0.1683998803
.0580*******
b（7,3）=0
b（7,4）=0
b（7,5）=0
b（7,6）=0
b（8,0）=0
b（8,1）=0.006552390126
b（8,2）=0.0005710218649
b（8,3）=0
b（8,4）=0
b（8,5）=0
b（8,6）=0
b（9,0）=193.6587558
b（9,1）=-1388.522425
b（9,2）=4126.607219
b（9,3）=-6508.211677
b（9,4）=5745.984054
b（9,5）=-2693.088365
b（9,6）=523.5718623
Dimbb（6To8,1To2）AsDouble
bb（6,1）=0.4006073948
bb（6,2）=0
bb（7,1）=0.08636081627
bb（7,2）=0
.853*******
bb（8,2）=0.3460208861
ot=（t+273.15）/647.3
op=p/22.12
SelectCasepara
Case1
x1=bb0*（1#-ot）
op0=Log（op）
op1=olp（ot）
op2=olp1（ot）
oh=a0+b（0,0）*ot
ot1=1/ot
Fori=1To5
ot1=ot1*ot
oh=oh-（i-2）*b（0,i）*ot1
Nexti
opp=1#
oh1=0#
Fori=1To5
oh0=0#
Forj=1Ton（i）
oh0=oh0+（1#+z（i,j）*bb0*ot）*b（i,j）*Exp（z（i,j）*x1）
Nextj
opp=opp*op
oh1=oh1+opp*oh0
Nexti
oh0=0#
Fori=6To8
oh3=0#
Forj=1Tol（i）
oh3=oh3+x（i,j）*bb（i,j）*Exp（x（i,j）*x1）
Nextj
oh4=0#
Forj=1Tol（i）
oh4=oh4+bb（i,j）*Exp（x（i,j）*x1）
Nextj
oh3=bb0*ot*oh3/（Exp（（2-i）*op0）+oh4）
oh5=0#
Forj=1Ton（i）
oh5=oh5+b（i,j）*Exp（z（i,j）*x1）*（（1#+z（i,j）*bb0*ot）-oh3）
Nextj
oh1=oh1+oh5/（Exp（（2-i）*op0）+oh4）
Nexti
oh=oh-oh1
oh1=0#
Fori=0To6
oh1=oh1+（1#+ot*（10#*op2/op1+
升级会员 