变工况课程设计

1、目录一、课程设计目的2二、课程设计内容2三、课程设计中使用的计算方法23.1 逐级热力核算方法23.2 整机热力核算3四、课程设计中的约定6五、N300-16.7/537/537汽轮机介绍75.1汽轮机本体和系统基本组成75.2 汽轮机N300主要结构参数和设计工况下的热力参数8六、计算106.1 初压不变P为60%流量下N300-16.7/537/537的变工况计算106.2 应用程序14七、分析21八、总结21一、课程设计目的（1） 通过课程设计，系统地总结、巩固、加深在汽轮机原理课程中已学到知识，进一步了解汽轮机的工作原理。（2） 在尽可能考虑制造、安装和运行的要求下，进行某一机组的变工

2、况热力计算，掌握汽轮机热力计算的原理、方法和步骤。（3） 通过课程设计对电站汽轮机建立整体的、量化的概念，掌握查阅的使用各种设计资料、标准、手册等参考文献的技巧。（4） 培养综合应用书本知识、自主学习、独立工作的能力，培养与其他人相互协作的工作作风。二、课程设计内容 以某种型号的汽轮机为对象，在已知结构参数和非设计工况新蒸汽参数和流量的条件下，进行通流部分热力校核计算，求出该变工况下级的内功率、相对内效率等全部特征参数，并与设计工况作对比分析。主要工作如下：（1） 非设计工况下通流部分各级热力过程参数计算。（2） 轴端汽封漏汽量校核计算。（3） 与设计工况性能和特征参数作比较分析。三、课程设计

3、中使用的计算方法3.1 逐级热力核算方法汽机变工况下的核算方法很多，当新工况偏离设计工况不远时，可采用近似估计方法。当汽轮机的新工况偏离设计工况较远时，或者在特殊工况下，就需要进行逐级核算求取级的各项参数。汽轮机逐级热力核算一般使用两种方法：一、从已知新汽参数开始，逐级向低压段进行核算，称为顺算法；二、从估算的终参数开始，逐级向高压段进行核算。就具体级而言，两种方法可概述如下：a) 由级前向级后核算由已知的级前初参数开始，根据流量比求得压力比，然后求出喷嘴后压力，求出相应的喷嘴理想比焓降，计算出喷嘴出口气流理想速度。然后进行动叶核算，方法与上述类同，求出新工况下动叶后压力，求出动叶理想比焓降，

4、计算出动叶出口气流速度。上述方法适用于在喷嘴和动叶后的压力，继续使用该方法进行计算，但需要根据已知条件进行假定压力的校核。一般来说，顺算法在超临界流动中很少。b) 由级后向级前推算从已知的级后压力开始，假设新工况下排气比焓、级后各项损失，求出动叶后参数，并对前面假设的余速损失进行校核及修正；继而假定动叶进口相对速度，求出动叶前参数、喷嘴后参数，再对动叶进口相对速度进行校核并修正；最后确定喷嘴前蒸汽状态点，并对最初估计的各项损失进行校核及修正；待各项假设假设都校核通过，即可求出新工况下级的反动度、级效率和内效率。判断动叶和喷嘴内是否出现超音速流动，当级内出现超音速流动时，还需确定临界状态点，并计

5、算出口气流偏转角。3.2 整机热力核算 多级汽轮机整机的热力核算可以从已知新汽参数开始，是顺算法，由高压级逐渐向低压级进行核算；也可同时使用逆算法，由已知的假定的排汽点逐渐向初参数靠拢。当计算的初比焓值（对喷嘴调节剂汽轮机，一般为第一压力级前比焓值）与已知的初比焓值不符，需对末级排汽比焓进行修正；当计算的初比焓高于已知出比焓时，说明假设的末级排汽比焓值过高，反之，说明末级排汽比焓假设的过低；在初比焓相差不大时，可采用平移热力过程曲线的方法是比焓值重合，若计算得出比焓值与已知的相差较大，则需重新假设末级比焓进行计算，直至满足精度要求。 在具体问题中究竟采用哪种方式进行核算，要根据所给定的新工况的

6、条件及要求的精确程度来决定。当新工况与设计工况下的各级均未出现超音速流动时，采用流量比压力比关系式，先确定各级喷嘴前后压力，再根据各级反动度与比焓降的变化规律确定各级喷嘴后压力，进行核算。当工况变动不大时，仅调节各级与一、二级的热力过程曲线有较大变化，可仅对末级进行详细核算，中间级通过流量与压力关系式确定级前压力，然后逐级平移热力过程线即可。这就是近似估算法。当工况变动很大时或级由亚音速流动转变为超音速流动，必须进行逐级核算顺算法的计算步骤：在确定满足使用顺算法的条件以后，可以根据下列计算步骤进行热力核算：（1） 求取喷嘴新工况的压力比，流量比及临界流量；临界流量 流量比 压力比 （2） 求取

7、新工况下喷嘴中的各项参数：喷嘴后压力 根据水蒸气性质求出喷嘴理想出口比焓喷嘴理想比焓降 喷嘴理想出口速度 喷嘴损失 喷嘴实际比焓降 求出喷嘴出口比焓，温度喷嘴实际出口 轮周速度 动叶进口相对速度 与速度对应的比焓 动叶进汽角 冲角 撞击损失 动叶进口滞止焓 根据水蒸气性质求得动叶进口滞止参数（3） 求取动叶的压力比，流量比，及临界流量； 临界流量 流量比 压力比 动叶后压力 （4）求新工况下的动叶中各项参数： 根据水蒸气热力性质求得 动叶理想比焓降 动叶出口理想速度 速度系数查表得到 动叶损失 动叶出口焓 根据水蒸气性质动叶出口温度 动叶实际比焓降 动叶出口相对速度 动叶出口绝对速度 动叶出汽

8、角： 画出速度三角形，求出速度三角形中各项参数 （5）求级的其他各项损失、内效率、内功率： 无限长叶片轮周比焓降 叶高损失 （将扇形损失一并考虑时， ） ：喷嘴高度，单位取,见附表3-7。 轮周比焓降： 级的理想可用能量： （为下一级对本级的余速利用系数） 扇形损失 ：动叶高度 ：动叶平均直径，以上参见附表3-7 动叶摩擦损失 （为圆周速度） 湿气损失 （为极平均蒸汽干度） 隔板漏气损失 ：隔板汽封平均直径：隔板汽封间隙：隔板汽封齿数：喷嘴出口面积，以上参见附表3-7叶頂漏气损失 ：与叶顶轴向间隙和围周带边厚度有关的经验系数：与叶顶轴向间隙及速度比有关的经验系数:与反动度及径高比有关的经验系数

9、：叶顶轴向间隙与叶高之比值，以上参数见附表3-7，余数损失 级有效焓降 相对内效率 级的内效率 本级完成计算四、课程设计中的约定（1） 各级回热抽汽量正比例于主汽流量；（2） 门杆漏气和调门开启重叠度不计；（3） 压力损失的参考值：高、中压连通管的压损为初压的5%，再热器的压损为高压缸排气压力的12%，中、低压连通管的压损为中压缸排气压力的2%。准确的数据可参考具体型号机组的资料。（4） 余速利用系数的参考值为：调节级后的第一压力级、前面有抽气口的压力级利用上一级余速的系数为0.4，其它压力级为0.8。（5） 第一次计算，用弗留格尔公式确定调节级后压力；（6） 假定汽机排气压力为设计工况下的值

10、，用平移设计工况热力过程线方法初步确定排气点。（7） 待通流部分计算结束后，校核计算轴端汽封的漏气量，轴端汽封的漏气量正比例于汽封前的压力。五、N300-16.7/537/537汽轮机介绍5.1汽轮机本体和系统基本组成 本机是按美国西屋公司的技术制造的300MW亚临界，中间再热式、高中压合缸、双缸双排汽、单轴凝气式汽轮机，如图4-2所示。它由高中压积木块BB0243与低压缸积木块BB074组合而成。为了进一步提高机组的经济性，对原引进技术做了改进设计，而且低压缸末级叶片采用905mm的长叶片。机组型号为N300-16.7/537/537，工厂产品号为D156. 主要技术参数：额定功率300MW

11、；主汽门前额定参数16.7MPa、537C，再热器汽门前温度537C；工作转速3000rpm；额定被压5.39KPa；回热级数3高、4低、1除氧；额定工况蒸汽流量910.2t/h、热耗7937KJ/(KWh)。 本机通流部分共35级叶片，其中高压缸1+11级，中压缸9级低压缸2*7级。动叶片除低压缸末三级为扭转叶片外，其余均为等截面叶片。高、中、低压缸隔板静叶均为扭转叶片。有六个高压调压调节阀，每阀控制21个喷嘴，调节阀全开时部分进气度为0.95455。汽轮机在额定参数下5阀全开时可以带额定负荷。汽轮机在6个调节阀全开、新汽参数16.7MPa、537C（超压5%），时运行，这一工况定义为最大负

12、荷工况。通流部分结构参数如表4-5所示。最大工况和额定工况下的热力参数如表4-6所示。 主汽门、调节阀、进气管的压损的4%，再热器及管道为10%，中联们及管道为2.5%，重低压连通管为2%。2号高加抽汽分别来自左右侧低压缸排汽，除氧器抽器来自中亚缸排气，第5、6级抽汽分别来自左、右侧低压缸的第3、5级前，因此低压缸为不对称抽汽。 静叶与动叶汽封：高压静叶3齿，平齿，中压静叶5齿，2高3低，动叶5齿，平齿，间隙都是0.75mm；低压静叶3齿，1高2低，动叶1齿，平齿，间隙都是1.1mm。 气缸抽汽汽封结构数据见表4-45.2 汽轮机N300主要结构参数和设计工况下的热力参数表3-5：N300-1

13、6.7/537/537汽轮机通流部分的结构参数缸别级号喷嘴动叶喷嘴数目平均直径出口高度出口面积叶片数目平均直径出口高度出口面积/mmmmcm2/mmmmcm2高压缸0.27671261061.322.9201.69.3490801064.327.94311.201.251572844.268.4456.24.249092854.268.82459.853.251692864.269.2472.69.253392874.269.56483.924.251994884.270.0489.81.253594894.470.52502.305.256696904.470.9515.90.2539969

14、15.071.96525.226.257198948.076.2583.47.2551100959.177.66596.957.2578100970.879.6625.86.2488106983.482.66635.398.2544104991.481.0641.80.24891081002.482.42646.029.2550106101484.5686.41.24671121026.787.50696.2510.2526110103990.3744.54.24481161051.793.36755.1511.2536112106495.7810.94.24591181077.099.368

15、26.71中压缸1.35141001078.9101.41207.7.35461021089.7103.81259.52.3566981100.5106.11308.3.35731021111.3108.51353.23.3589981122.0110.91402.4.3607761134.8113.21456.74.3672941145.6115.61527.3.3672741159.8119.31595.85.3676961173.0124.11680.8.3601801190.7131.31768.66.34691161204.9137.11800.0.3482901220.2142.0

16、1895.47.34741181233.4146.81976.1.3485921251.1154.02109.78.34971201265.3159.92245.7.35117921283.0167.12368.29.3522901302.1176.12567.2.3487781327.4188.62742.5低压缸1.22301961915.286.361158.7.25352021921.090.681380.12.25252061929.099.061511.2.28251421937.4107.11827.43.29242261953.0122.72187.8.29701961971.

17、9141.62578.44.29942201993.9163.63051.8.32201262018.9188.63800.95.31931382051.2220.94465.3.33101242096.3265.85656.56.2450802250.7421.97119.5.33501202299.8474.6114197.3203662604.0845.821128.4839962642.9905.536380表3-6：N300/-16.7/537/537汽轮机最大工况和设计工况下的热力参数缸别级号最大工况额定工况级前压力MPa级前温度内功率kW级流量t/h级前压力MPa级前温度内功率K

18、W高压缸016.853716930905.116.03534.318555113.1492.76087885.611.6483.85365212.1479.76365885.610.7471.15608311.1466.26805885.69.87457.95992410.1451.87191885.69.01443.8632559.18436.67456885.68.17428.9655168.28420.97815885.67.83413.5685577.41404.58219885.66.61397.6719486.59387.57163823.85.89381.0634595.8837

19、1.57521754.35.26365.26653105.21354.77849754.34.66348.46930114.58337.38180754.34.10331.87203中压缸13.47533.68723754.33.12533.6781823.10515.38987754.32.78515.3805232.75496.39514754.32.47496.4883142.41476.19869754.32.16476.4883152.10455.110245754.31.88455.5915861.81433.210093721.21.63433.7906271.55410.510

20、758721.21.39411.1965181.30386.211469721.21.17386.91027691.08360.211927721.20.970361.110666低压缸10.852330.39545329.30.769331.4861020.599286.210227329.30.541287.2919730.412242.79501311.80.372243.7858540.257193.610241311.80.232194.5936250.147140.09152301.90.133140.8829760.081794.614691285.50.073691.81322

21、270.024365.612782273.50.022863.411308六、计算6.1 初压不变P为60%流量下N300-16.7/537/537的变工况计算 1）喷嘴设计工矿核心工况下的压力比、，流量比 2）求取新工况下喷嘴中的各项参数： 喷嘴后压力 根据水蒸气求出喷嘴理想出口比焓 喷嘴中滞止理想比焓 喷嘴出口速度 喷嘴损失 喷嘴实际比焓降 出口比焓 喷嘴实际出口速度 轮周速度 ： 动叶进口相对速度 与速度对应的比焓 =7.06 动叶进汽角 冲角 撞击损失 动叶进口滞止焓 根据水蒸气性质求得动叶进口滞止参数 3） 求取动叶的压力比，流量比，及临界流量设计工况下，动叶进口状态点比焓，比熵由程

22、序得：滞止压力临界流量 流量比 压力比 动叶后压力 4） 求新工况下的动叶中各项参数： 根据水蒸气热力性质求得 由程序得 所以 动叶理想比焓降 动叶出口理想速度 速度系数查表得到 动叶损失 动叶出口焓 根据水蒸气性质动叶出口温度 动叶实际比焓降 动叶出口相对速度 动叶出口绝对速度 动叶出汽角 5） 求级的其他各项损失、内效率、内功率： 无限长叶片轮周比焓降 叶高损失（将扇形损失一并考虑时， ）：喷嘴高度，单位取mm,见附表3-7。轮周比焓降：级的理想可用能量：（为下一级对本级的余速利用系数）扇形损失：动叶高度：动叶平均直径，以上参见附表3-7动叶摩擦损失（为圆周速度）湿气损失（为极平均蒸汽干度

23、）隔板漏气损失：隔板汽封平均直径：隔板汽封间隙：隔板汽封齿数：喷嘴出口面积，以上参见附表3-7叶頂漏气损失 ：与叶顶轴向间隙和围周带边厚度有关的经验系数：与叶顶轴向间隙及速度比有关的经验系数:与反动度及径高比有关的经验系数：叶顶轴向间隙与叶高之比值，以上参数见附表3-7，余数损失级有效焓降相对内效率级的内功率6.2 应用程序第一程序function vsteam,hsteam,ssteam=ptsteam(p,t)function vsteam,hsteam,ssteam=ptsteam(p,t)%ptsteam.m%caculating v,h,s from pressure and tem

24、perature knownR=0.461526;JO=0,1,-5,-4,-3,-2,-1,2,3;nO=-0.96927686500217e1;0.10086655968018e2;-0.56087911283020e-2;0.71452738081455e-1;-0.40710498223928;0.14240819171444e1;-0.43839511319450e1;-0.28408632460772;0.21268463753307e-1;pr=p;tr=540/t;rO=log(pr);rOt=0;for i=1:9 rO=rO+nO(i)*trJO(i); rOt=rOt+n

25、O(i)*JO(i)*tr(JO(i)-1);end rOp=pr(-1);I=1;1;1;1;1;2;2;2;2;2;3;3;3;3;3;4;4;4;5;6;6;6;7;7;7;8;8;9;10;10;10;16;16;18;20;20;20;21;22;23;24;24;24;J=0;1;2;3;6;1;2;4;7;36;0;1;3;6;35;1;2;3;7;3;16;35;0;11;25;8;36;13;4;10;14;29;50;57;20;35;48;21;53;39;26;40;58;n=-0.17731742473213e-2;-0.17834862292358e-1;-0.45

26、996013696365e-1;-0.57581259083432e-1;-0.50325278727930e-1;-0.33032641670203e-4;-0.18948987516315e-3;-0.39392777243355e-2;-0.43797295650573e-1;-0.26674547914087e-4;0.20481737692309e-7;0.43870667284435e-6;-0.32277677238570e-4;-0.15033924542148e-2;-0.40668253562649e-1;-0.78847309559367e-9;0.12790717852

27、285e-7;0.48225372718507e-6;0.22922076337661e-5;-0.16714766451061e-10;-0.21171472321355e-2;-0.23895741934104e2;-0.59059564324270e-17;-0.12621808899101e-5;-0.38946842435739e-1;0.11256211360459-10;-0.82311340897998e1;0.19809712802088e-7;0.10406965210174e-18;-0.10234747095929e-12;-0.10018179379511e-8;-0

28、.80882908646985e-10;0.10693031879409;-0.33662250574171;0.89185845355421e-24;0.30629316876232e-12;-0.42002467698208e-5;-0.59056029685639e-25;0.37826947613457e-5;-0.12768608934681e-14;0.73087610595061e-28;0.55414715350778e-16;-0.94369707241210e-6;rt=0;rtt=0;rtp=0;for i=1:43 rt=rt+n(i)* prI(i)*(tr-0.5)

29、J(i); rtt=rtt+n(i)* prI(i)* J(i)*(tr-0.5)(J(i)-1); rtp=rtp+n(i)*I(i)* pr(I(i)-1)*(tr-0.5)J(i);endvsteam=pr*(r0p+rtp)/p*R*t/1000;hsteam=tr*(r0t+rtt)*R*t;ssteam=(tr*(r0t+rtt)-(r0+rt)*R;第二程序：function vwater,hwater,swater=ptwater(p,t)function vwater,hwater,swater=ptwater(p,t)%ptwater,m%caculating v,h,s

30、from pressure and temperature knownR=0.461526;I=0,0,0,0,0,0,0,0,1,1,1,1,1,1,2,2,2,2,2,3,3,3,4,4,4,5,8,8,21,23,29,30,31,32;J=-2,-1,0,1,2,3,4,5,-9,-7,-1,0,1,3,-3,0,1,3,17,-4,0,6,-5,-2,10,-8,-11,-6,-29,-31,-38,-39,-40,-41;n=0.14632971213167;-0.84548187169114;-0.37563603672040e1;0.33855169168385e1;-0.95

31、791963387872;0.15772038513228;-0.16616417199501e-1;0.81214629983568e-3;0.28319080123804e-3;-0.60706301565874e-3;-0.18990068218419e-1;-0.32529748770505e-1;-0.21841717175414e-1;-0.52838357969930e-4;-0.47184321073267e-3;-0.30001780793026e-3;0.47661393906987e-4;-0.44141845330846e-5;-0.72694996297594e-15

32、;-0.31679644845054e-4;-0.28270797985312e-5;-0.85205128120103e-9;-0.22425281908000e-5;-0.65171222895601e-6;-0.14341729937924e-12;-0.40516996860117e-6;-0.12734301741641e-8;-0.17424871230634e-9;-0.68762131295531e-18;0.144783078288521e-19;0.26335781662795e-22;-0.11947622640071e-22;0.18228094581404e-23;-

33、0.93537087292458e-25;% caculating specific volumepr=p/16.53;tr=1386/t;r=0;rt=0;rp=0;for i=1:34 r=r+n(i)*(7.1-pr)I(i)*(tr-1.222)J(i); rt=rt+n(i)*(7.1-pr)I(i)*J(i)*(tr-1.222)(J(i)-1); rp=rp-n(i)*I(i)*(7.1-pr)(I(i)-1)*(tr-1.222)J(i);vwater=pr*rp/p*R*t/1000;hwater=tr*rt*R*t;swater=(tr*rt-r)*R;第三程序：funct

34、ion ts=tsaturation(p)n=0.11670521452767e4;-0.72421316703206e6;-0.17073846940092e2;0.12020824702470e5-0.32325550322333e7;0.14915108613530e2;-0.48232657361591e4;0.40511340542057e6;-0.23855557567849;0.65017534844798e3;p=p0.25;E=p2+n(3)*p+n(6);F=n(1)*p2+n(4)*p+n(7);G=n(2)*p2+n(5)*p+n(8);D=2*G/(-F-(F2-4*

35、E*G)5);ts=(n(10)+D-(n(10)+D)2-4*(n(9)+n(10)*D)0.5)/2第四程序：function ps=psaturation(t)n=0.11670521452767e4;-0.72421316703203e6;-0.17073846940092e2; 0.12020824702470e5;-0.32325550322333e7;0.14915108613530e2; -0.48232657361591e4;0.40511340542057e6;-0.23855557567849; 0.65017534844798e3;t=t+n(9)/(t-n(10);A

36、=t2+n(1)*t+n(2);B=n(3)*t2+n(4)*t+n(5);C=n(6)*t2+n(7)*t+n(8);ps=(2*C/(-B+(B2-4*A*C)0.5)4;第五程序：functiont,v,s,x =ph(p,h)% ph.m% caculating t,v,s from pressure and enthalpy knownts=tsaturation(p);vwater,hwater,swater=ptwater(p,ts);vsteam,hsteam,ssteam=ptsteam(p,ts);if h<=hwater'supercooled water&

37、#39;I=0,0,0,0,0,0,1,1,1,1,1,1,1,2,2,3,3,4,5,6;J=0,1,2,6,22,32,0,1,2,3,4,10,32,10,32,10,32,32,32,32;n=-0.23872489924521e3;0.40421188637945e3;0.11349746881718e3;-0.58457616048039e1;-0.15285482413140e-3;-0.10866707695377e-5; -0.13391744872601e2;0.43211039183559e2;-0.54010067170506e2; 0.30535892203916e2

38、;-0.65964749423638e1;0.93965400878363e-2;0.11573647505340e-6;-0.25858641282073e-4;-0.40644363084799e-8;0.66456186191635e-7;0.80670734103027e-10;-0.93477771213947e-12; 0.58265442020601e-14;-0.15020185953503e-16;hr=h/2500;t=0;for i=1:20 t=t+n(i)*pI(i)*(hr+1)J(i);endv,h,s=ptwater(p,t);x=t-ts;elseif h&g

39、t;=hsteam 'superheated steam' hr=h/2000;N=0.90584278514723e3;-0.67955786399241;0.12809002730136e-3;0.26526571908428e4;0.45257578905948e1;if p<4 '2a sector' I=0;0;0;0;0;0;1;1;1;1;1;1;1;1;1;2;2;2;2;2;2;2;2;3;3;4;4;4;5;5;5;6;6;7; J=0;1;2;3;7;20;0;1;2;3;7;9;11;18;44;0;2;7;36;38;40;42;

40、44;24;44;12;32;44;32;36;42;34;44;28;n=0.10898952318288e4;0.84951654495535e3;-0.10781748091826e3; 0.33153654801263e2;-0.74232016790248e1;0.11765048724356e2; 0.18445749355790e1;-0.41792700549624e1;0.62478196935812e1; -0.17344563108114e2;-0.20058176862096e3;0.27196065473796e3; -0.45511318285818e3;0.309

41、1968804755e4;0.25226640357872e6; -0.61707422868339e-2;-0.31078046629583;0.11670873077107e2; 0.12812798404046e9;-0.98554909623276e9;0.28224546973002e10; -0.35948971410703e10;0.17227349913197e10;-0.13551334240775e5; 0.12848734664650e8;0.13865724283226e1;0.23598832556514e6; -0.13105236545054e8;0.739998

42、35474766e4;-0.55196697030060e6; 0.37154085996233e7;0.19127729239660e5;-0.41535164835634e6; -0.62459855192507e2;t=0;for i=1:34t= t+n(i)*pI(i)*(hr-2.1)J(i);endelseif p<N(1)+ N(2)*h+N(3)*h2 '2b sector' I=0;0;0;0;0;0;0;0;1;1;1;1;1;1;1;1;2;2;2;2;3;3;3;3;4;4;4;4;4;4;5;5;5;6;7;7;9;9;J=0;1;2;12;1

43、8;24;28;40;0;2;6;12;18;24;28;40;2;8;18;40;1;2;12;24;2;12;18;24;28;40;18;24;40;28;2;28;1;40;n=0.14895041079516e4;0.74307798314034e3;-0.97708318797837e2; 0.24742464705674e1;-0.63281320016026;0.11385952129658e1; -0.47811863648625;0.85208123431544e-2;0.93747147377932; 0.33593118604916e1;0.33809355601454

44、e1;0.16844539671904; 0.73875745236695;-0.47128737436186;0.15020273139707; -0.21764114219750e-2;-0.21810755324761e-1;-0.10829784403677; -0.46333324635812e-1;0.71280351959551e-4;0.11032831789999e-3; 0.18955248387902e-3;0.30891541160537e-2;0.13555504554949e-2; 0.28640237477456e-6;-0.10779857357512e-4;-

45、0.76462712454814e-4; 0.14052392818316e-4;-0.31083814331434e-4;-0.10302738212103e-5; 0.28217281635040e-6;0.12704902271945e-5;0.73803353468292e-7; -0.11030139238909e-7;-0.81456365207833e-13;-0.25180545682962e-10; -0.17565233969407e-17;0.86934156344163e-14;t=0;for i=1:38t=t+n(i)*(p-2)I(i)*(hr-2.6)J(i);

46、endelse'2c sector'I=-7;-7;-6;-6;-5;-5;-2;-2;-1;-1;0;0;1;1;2;6;6;6;6;6;6;6;6;J=0;4;0;2;0;2;0;1;0;2;0;1;4;8;4;0;1;4;10;12;16;20;22;n=-0.32368398555242e13;0.73263350902181e13;0.35825089945447e12; -0.58340131851590e12;-0.10783068217470e11;0.20825544563171e11; 0.61074783564516e6;0.85977722535580e

47、6;-0.25745723604170e5; 0.31081088422714e5;0.12082315865936e4;0.48219755109255e3; 0.37966001272486e1;-0.10842984880077e2;-0.45364172676660e-1; 0.14559115658698e-12;0.11261597407230e-11;0.17804982240686e-10; 0.12324579690832e-6;-0.11606921130984e-5;0.2846367088554e-4;-0.59270038474176e-3;0.12918582991

48、878e-2;t=0;for i=1:23 t=t+n(i)*(p+25)I(i)*(hr-1.8)J(i); end end v,h,s=ptsteam(p,t); x=t-ts;else 'wet steam' x=(h-hwater)/(hsteam-hwater) t=ts; v=vwater+x*(vsteam-vwater); s=swater+ x*(ssteam-swater);end第六程序：functiont,v,h,x =ps(p,s)% ps.m% caculating t,v,s from pressure and enthalpy knownts=t

49、saturation(p);vwater,hwater,swater=ptwater(p,ts);vsteam,hsteam,ssteam=ptsteam(p,ts);if s<=swater'supercooled water' I=0;0;0;0;0;0;1;1;1;1;1;1;2;2;2;2;2;3;3;4; J=0;1;2;3;11;31;0;1;2;3;12;31;0;1;2;9;31;10;32;32; n=0.17478268058307e3;0.34806930892873e2;0.65292584978455e1; 0.33039981775489;-0

50、.19281382923196e-6;-0.24909197244573e-22; -0.26107636489332;0.22592965981586;-0.64256463395226e-1; 0.78876289270526e-2;0.35672110607366e-9;0.17332496994895e-23; 0.56608900654837e-3;-0.32635483139717e-3;0.44778286690632e-4; -0.51322156908507e-9;-0.42522657042207e-25;0.26400441360689e-12; 0.7812460045972