PMV-VB计算代码

以下程序由VB编写：ta=Val(ta0)“输入室内温度”fi=Val(fi0)“输入湿度”va=Val(va0)“输入室内空气流速”tr=Val(tr0)“输入室内平均辐射温度”icl=0.155*Val(icl0)“输入服装热阻”m=Val(m0)“输入人体代谢强度”w=Val(w0)“输入做功强度”pa=Exp(16.6536-4030.183/(ta+235))*fi*10“求水蒸气分压力”Ificl0.078Thenfcl=1+1.29*iclElsefcl=1.05+0.645*iclEndIf“求穿衣服与裸体时的人体表面积比”hcf=12.1*va^0.5mw=m-wtra=tr+273taa=ta+273tcla=taa+(35.5-ta)/(3.5*(6.45*icl*0.15+0.1))xn=tcla/100p1=icl*fclp2=p1*3.96p3=p1*100p4=p1*taap5=308.7-0.028*mw+p2*(tra/100)^4“------------------------------以下为迭代求服装表面温度tcl---------------------------”Doxf=xnxf=(xf+xn)/2hcn=2.38*Abs(100*xf-taa)^0.25IfhcfhcnThenhc=hcfElsehc=hcnEndIf“求对流换热系数”xn=(p5+p4*hc-p2*xf^4)/(100+p3*hc)“迭代公式”eps=0.0001“迭代收敛条件”LoopUntilAbs(xn-xf)epstcl=100*xn–273“服装表面温度摄氏度”ts=0.303*Exp(-0.036*m)+0.028hl1=3.05*0.001*(5733-6.99*mw-pa)Ifmw58.15Thenhl2=0.42*(mw-58.15)Elsehl2=0EndIfhl3=1.7*0.00001*m*(5867-pa)hl4=0.0014*m*(34-ta)hl5=3.96*fcl*(xn^4-(tra/100)^4)hl6=fcl*hc*(tcl-ta)pmv0=ts*(mw-hl1-hl2-hl3-hl4-hl5-hl6)“求PMV”PMV=Format(pmv0,0.0)“设置PMV的显示格式”PPD=100-95*Exp(-0.03353*PMV^4-0.2179*PMV^2)“求PPD”