水蒸汽计算说明

1、97水蒸汽计算说明1在计算中时常返回一个数组，其数组的格式如下，相应的字母，在函数命名上会体现。使用时，只需要将相应的dll文件复制到mathlab下的work里面即可。因为可以同时计算定压热容、定容热容、单位体积比容、声速、热力学能等参数。如果你计算中需要使用这些参数，给我说一声，我再给出相应的函数。字母类型单位t温度p压强MPax干度，-1表示非饱和状态，（01）表示正常数据。h焓Kj/kgs熵Kj/Cp定压热容Kj/(kg.)Cv定容热容Kj/(kg.)函数列表函数条件返回值错误判断pT（p）->t饱和区参数不对返回-273.15tP(t)->p饱和区参数不对返回0pxH(p

2、,x)->H饱和区if（x<0）x=0和if（x>1）x=1，如果p不在饱和区，则返回0pxS(p,x)->S饱和区if（x<0）x=0和if（x>1）x=1，如果p不在饱和区，则返回0txH(t,x)->H饱和区if（x<0）x=0和if（x>1）x=1，如果t不在饱和区，则返回0txS(t,x)->S饱和区if（x<0）x=0和if（x>1）x=1，如果t不在饱和区，则返回0AphT（p，h）->t所有区参数不对返回-273.15ApsT(p，s) ->t所有区参数不对返回-273.15AthP（t，h）-

3、>p所有区参数不对返回0AtsP（t，s）->p所有区参数不对返回0AtpH（t，p）->h除饱和区的其它区参数不对返回0AtpS（t，p）->s除饱和区的其它区参数不对返回0AtpCp(t,p)->Cp除饱和区的其它区参数不对返回0AtpCv(t,p)->Cv除饱和区的其它区参数不对返回0phA（p，h）->（t,p,x,h,s）所有区参数不对，返回（0,0,-1,0,0）,不在饱和区x=-1psA（p，s）->（t,p,x,h,s）所有区参数不对，返回（0,0,-1,0,0）,不在饱和区x=-1thA（t，h）->（t,p,x,h,s）

4、所有区参数不对，返回（0,0,-1,0,0）,不在饱和区x=-1tsA（t，s）->（t,p,x,h,s）所有区参数不对，返回（0,0,-1,0,0）,不在饱和区x=-1tpA（t，p）->（t,p,x,h,s）除饱和区的其它区参数不对，返回（0,0,-1,0,0）ftpH(t,p)->h水区，f表示fluid参数不对返回0，请链接ftpS(t,p)->S水区，f表示fluid参数不对返回0，请链接gtpH(t,p) ->h汽区，f表示gas参数不对返回0，请链接gtpS(t,p) ->S汽区，f表示gas参数不对返回0，请链接tab=5 4 3 2;0.5

5、 0.4 0.3 0.2x=0.25Lagr（tab,x）->y知道对应关系，根据拉格朗日内插出y值详细请链接ExpEf(t0,p0,t1,p1)->efef是膨胀效率详细请链接SmokeH(t,RO2,NO2,H2O,Vair,K)->hy hy是烟气焓详细请链接dif(level,different,matrix)->有效平均值level：基准值different：偏差matrix：数据矩阵详细请链接2 考虑到现场温度和压强，比如快到饱和区的水的参数温度t，在存有误差的情况下，可能会跳到饱和温度上面去，这时用Ctp计算时，会计算气态的参数，所以提供了2个函数，ftp

6、H和ftpS专门计算水，gtpH和gtpS专门计算气态的参数和临界以后的参数，计算方法是根据原有的趋势，计算相应的值，比如ftpH，存有一个h与（t，p）的函数关系，当t和p在处于水状态下，是完全对的，但是当t和p偏离水状态下，其原有的h与（t，p）的函数关系，在（t，p）偏离不是很大的情况下，依然有很好的准确性，所以会根据原有的函数关系，计算出相应的值。注意的是，这两类主要用在温度低于临界点373.946以下的处于水和气边界下，存有汽化潜热的条件下，如果当超出这个温度范围，即在临界温度（t0,p0）即（373.946，22.064MPa）以下，当仪器测得的t0多了一点点误差，会因计算处气态的

7、参数，此时存在汽化潜热，会引入特别大误差，当处于（373.946，22.064MPa），因为不存在汽化潜热，所以，计算后的结果，跟实际的误差很小。即主要用在除氧器和凝汽器这两个存有临界区的条件下参数计算上。3水蒸气或水状态下的计算，不进行边界检测，主要考虑现场的误差，注意尽量此函数用于计算临界点（t0，p0）以下的情况。当然也可以计算临界点以上。其中ftpH和ftpS覆盖了区域1，gtpH和gtpS覆盖了区域3，2，5，计算时，会先判断是否是在区域5内，如果在就安5号区域的计算方法计算结果。然后判断是否在区域3内，如果在，就安3号区域的方法计算结果，否则就安2号区域的方法机型计算。4ftpH(

8、t,p):由温度和压强计算水的焓值，当判断温度和压强不是水状态时，会安旧有的关系，计算出相应的焓值，主要是考虑到现场测量仪器存有误差。5ftpS(t,p): 由温度和压强计算水的熵值，当判断温度和压强不是水状态时，会安旧有的关系，计算出相应的熵值，主要是考虑到现场测量仪器存有误差。6gtpH(t,p): 由温度和压强计算汽的焓值，当判断温度和压强不是汽状态时，会安旧有的关系，计算出相应的焓值，主要是考虑到现场测量仪器存有误差。7gtpS(t,p): 由温度和压强计算汽的熵值，当判断温度和压强不是汽状态时，会安旧有的关系，计算出相应的熵值，主要是考虑到现场测量仪器存有误差。9赠送一个函数：这个是

9、拉格朗日插值方法，即线性插值，每一段看成是直线。免去我们还要在excel做捏合的麻烦。比如：主气流量t/h214120781411.48928.96643.714主汽压力MPa25.425.3220.6713.98.99Mathlab语言计算：tab=2141 2078 1411.48 928.96 643.714; 25.4 25.32 20.67 13.9 8.99lagr(tab， 1500)即计算出对应于主气流量看为1500t/h时的线性插值后的压强值，注意l 此函数的tab第一行数组，在这里是主气流量，必须是从大到小的顺序，否则会报错不予计算。l 当对应的流量，比如在这里lagr（t

10、ab，2200），此时的流量已经超过了tab里的上下限，如果超出上限会根据第一段的斜率，来计算2200时的值，如果超出下限，会根据最后一段的斜率进行计算。l 因为我们在用Excel捏合计算时，有些基本点不能经过，误差或太大时，用此函数比较合适。l 如果需要三次样条曲线捏合和Hermit插值的，再联系我。10ExpEf(t0,p0,t1,p1)，计算从参照点(t0,p0)膨胀到（t1,p1）点的膨胀效率，在计算高压缸和中压缸内效率的时候有用，因为低压缸排气可能出于湿蒸汽区，不能用（t1，p1） 来唯一确定一参照点，故不能用此计算。需要手动计算。在参数不对时，会返回-1，在使用中，必须注意，(t0

11、，p0)在(t1,p1)的膨胀前，否则计算结果任然安等效线来计算。如图：11SmokeH(t,RO2,NO2,H2O,Vair,K)由温度，以及RO2的含量(标况下m3/kg)，NO2的含量(标况下m3/kg)，H2O的含量(标况下m3/kg)，及空气的含量Vair(标况下m3/kg)，以及过量空气系数，计算出锅炉烟气的焓值Kj/kg。即结果=RO2*（RO2单位体积焓值）+ NO2*（NO2单位体积焓值）+H2O*（H2O单位体积焓值）+（k-1）*(Vair*空气标况单位体积焓值)，如下图，即输入蓝色部分，得到红色部分结果。12dif(level,different,valmatrix)该函数计算一组数据的平均值，该值放在valmatrix参数中，该参数可以是矩阵，比如566,567,568;pms1,pms2,12.5，其中pms1和pms2是变量。Level是基准值，different是偏离值，当level-valmatrx的绝对值超过different，那么会忽略改点，如果没有一个合适点，那么就用level替换掉。举例