蒸汽特性表.docx
- 文档编号:5352395
- 上传时间:2023-05-08
- 格式:DOCX
- 页数:33
- 大小:20.46KB
蒸汽特性表.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+
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 蒸汽 特性