變工況課程設(shè)計

1、目錄一、課程設(shè)計目的2二、課程設(shè)計內(nèi)容2三、課程設(shè)計中使用的計算方法23.1 逐級熱力核算方法23.2 整機熱力核算3四、課程設(shè)計中的約定6五、N300-16.7/537/537汽輪機介紹75.1汽輪機本體和系統(tǒng)基本組成75.2 汽輪機N300主要結(jié)構(gòu)參數(shù)和設(shè)計工況下的熱力參數(shù)8六、計算106.1 初壓不變P為60%流量下N300-16.7/537/537的變工況計算106.2 應(yīng)用程序14七、分析21八、總結(jié)21一、課程設(shè)計目的（1） 通過課程設(shè)計，系統(tǒng)地總結(jié)、鞏固、加深在汽輪機原理課程中已學(xué)到知識，進一步了解汽輪機的工作原理。（2） 在盡可能考慮制造、安裝和運行的要求下，進行某一機組的變工

2、況熱力計算，掌握汽輪機熱力計算的原理、方法和步驟。（3） 通過課程設(shè)計對電站汽輪機建立整體的、量化的概念，掌握查閱的使用各種設(shè)計資料、標(biāo)準(zhǔn)、手冊等參考文獻的技巧。（4） 培養(yǎng)綜合應(yīng)用書本知識、自主學(xué)習(xí)、獨立工作的能力，培養(yǎng)與其他人相互協(xié)作的工作作風(fēng)。二、課程設(shè)計內(nèi)容 以某種型號的汽輪機為對象，在已知結(jié)構(gòu)參數(shù)和非設(shè)計工況新蒸汽參數(shù)和流量的條件下，進行通流部分熱力校核計算，求出該變工況下級的內(nèi)功率、相對內(nèi)效率等全部特征參數(shù)，并與設(shè)計工況作對比分析。主要工作如下：（1） 非設(shè)計工況下通流部分各級熱力過程參數(shù)計算。（2） 軸端汽封漏汽量校核計算。（3） 與設(shè)計工況性能和特征參數(shù)作比較分析。三、課程設(shè)計

3、中使用的計算方法3.1 逐級熱力核算方法汽機變工況下的核算方法很多，當(dāng)新工況偏離設(shè)計工況不遠(yuǎn)時，可采用近似估計方法。當(dāng)汽輪機的新工況偏離設(shè)計工況較遠(yuǎn)時，或者在特殊工況下，就需要進行逐級核算求取級的各項參數(shù)。汽輪機逐級熱力核算一般使用兩種方法：一、從已知新汽參數(shù)開始，逐級向低壓段進行核算，稱為順?biāo)惴?；二、從估算的終參數(shù)開始，逐級向高壓段進行核算。就具體級而言，兩種方法可概述如下：a) 由級前向級后核算由已知的級前初參數(shù)開始，根據(jù)流量比求得壓力比，然后求出噴嘴后壓力，求出相應(yīng)的噴嘴理想比焓降，計算出噴嘴出口氣流理想速度。然后進行動葉核算，方法與上述類同，求出新工況下動葉后壓力，求出動葉理想比焓降，

4、計算出動葉出口氣流速度。上述方法適用于在噴嘴和動葉后的壓力，繼續(xù)使用該方法進行計算，但需要根據(jù)已知條件進行假定壓力的校核。一般來說，順?biāo)惴ㄔ诔R界流動中很少。b) 由級后向級前推算從已知的級后壓力開始，假設(shè)新工況下排氣比焓、級后各項損失，求出動葉后參數(shù)，并對前面假設(shè)的余速損失進行校核及修正；繼而假定動葉進口相對速度，求出動葉前參數(shù)、噴嘴后參數(shù)，再對動葉進口相對速度進行校核并修正；最后確定噴嘴前蒸汽狀態(tài)點，并對最初估計的各項損失進行校核及修正；待各項假設(shè)假設(shè)都校核通過，即可求出新工況下級的反動度、級效率和內(nèi)效率。判斷動葉和噴嘴內(nèi)是否出現(xiàn)超音速流動，當(dāng)級內(nèi)出現(xiàn)超音速流動時，還需確定臨界狀態(tài)點，并計

5、算出口氣流偏轉(zhuǎn)角。3.2 整機熱力核算 多級汽輪機整機的熱力核算可以從已知新汽參數(shù)開始，是順?biāo)惴?，由高壓級逐漸向低壓級進行核算；也可同時使用逆算法，由已知的假定的排汽點逐漸向初參數(shù)靠攏。當(dāng)計算的初比焓值（對噴嘴調(diào)節(jié)劑汽輪機，一般為第一壓力級前比焓值）與已知的初比焓值不符，需對末級排汽比焓進行修正；當(dāng)計算的初比焓高于已知出比焓時，說明假設(shè)的末級排汽比焓值過高，反之，說明末級排汽比焓假設(shè)的過低；在初比焓相差不大時，可采用平移熱力過程曲線的方法是比焓值重合，若計算得出比焓值與已知的相差較大，則需重新假設(shè)末級比焓進行計算，直至滿足精度要求。 在具體問題中究竟采用哪種方式進行核算，要根據(jù)所給定的新工況的

6、條件及要求的精確程度來決定。當(dāng)新工況與設(shè)計工況下的各級均未出現(xiàn)超音速流動時，采用流量比壓力比關(guān)系式，先確定各級噴嘴前后壓力，再根據(jù)各級反動度與比焓降的變化規(guī)律確定各級噴嘴后壓力，進行核算。當(dāng)工況變動不大時，僅調(diào)節(jié)各級與一、二級的熱力過程曲線有較大變化，可僅對末級進行詳細(xì)核算，中間級通過流量與壓力關(guān)系式確定級前壓力，然后逐級平移熱力過程線即可。這就是近似估算法。當(dāng)工況變動很大時或級由亞音速流動轉(zhuǎn)變?yōu)槌羲倭鲃?，必須進行逐級核算順?biāo)惴ǖ挠嬎悴襟E：在確定滿足使用順?biāo)惴ǖ臈l件以后，可以根據(jù)下列計算步驟進行熱力核算：（1） 求取噴嘴新工況的壓力比，流量比及臨界流量；臨界流量 流量比 壓力比 （2） 求取

7、新工況下噴嘴中的各項參數(shù)：噴嘴后壓力 根據(jù)水蒸氣性質(zhì)求出噴嘴理想出口比焓噴嘴理想比焓降 噴嘴理想出口速度 噴嘴損失 噴嘴實際比焓降 求出噴嘴出口比焓，溫度噴嘴實際出口 輪周速度 動葉進口相對速度 與速度對應(yīng)的比焓 動葉進汽角 沖角 撞擊損失 動葉進口滯止焓 根據(jù)水蒸氣性質(zhì)求得動葉進口滯止參數(shù)（3） 求取動葉的壓力比，流量比，及臨界流量； 臨界流量 流量比 壓力比 動葉后壓力 （4）求新工況下的動葉中各項參數(shù)： 根據(jù)水蒸氣熱力性質(zhì)求得 動葉理想比焓降 動葉出口理想速度 速度系數(shù)查表得到 動葉損失 動葉出口焓 根據(jù)水蒸氣性質(zhì)動葉出口溫度 動葉實際比焓降 動葉出口相對速度 動葉出口絕對速度 動葉出汽

8、角： 畫出速度三角形，求出速度三角形中各項參數(shù) （5）求級的其他各項損失、內(nèi)效率、內(nèi)功率： 無限長葉片輪周比焓降 葉高損失 （將扇形損失一并考慮時， ） ：噴嘴高度，單位取,見附表3-7。 輪周比焓降： 級的理想可用能量： （為下一級對本級的余速利用系數(shù)） 扇形損失 ：動葉高度 ：動葉平均直徑，以上參見附表3-7 動葉摩擦損失 （為圓周速度） 濕氣損失 （為極平均蒸汽干度） 隔板漏氣損失 ：隔板汽封平均直徑：隔板汽封間隙：隔板汽封齒數(shù)：噴嘴出口面積，以上參見附表3-7葉頂漏氣損失 ：與葉頂軸向間隙和圍周帶邊厚度有關(guān)的經(jīng)驗系數(shù)：與葉頂軸向間隙及速度比有關(guān)的經(jīng)驗系數(shù):與反動度及徑高比有關(guān)的經(jīng)驗系數(shù)

9、：葉頂軸向間隙與葉高之比值，以上參數(shù)見附表3-7，余數(shù)損失 級有效焓降 相對內(nèi)效率 級的內(nèi)效率 本級完成計算四、課程設(shè)計中的約定（1） 各級回?zé)岢槠空壤谥髌髁浚唬?） 門桿漏氣和調(diào)門開啟重疊度不計；（3） 壓力損失的參考值：高、中壓連通管的壓損為初壓的5%，再熱器的壓損為高壓缸排氣壓力的12%，中、低壓連通管的壓損為中壓缸排氣壓力的2%。準(zhǔn)確的數(shù)據(jù)可參考具體型號機組的資料。（4） 余速利用系數(shù)的參考值為：調(diào)節(jié)級后的第一壓力級、前面有抽氣口的壓力級利用上一級余速的系數(shù)為0.4，其它壓力級為0.8。（5） 第一次計算，用弗留格爾公式確定調(diào)節(jié)級后壓力；（6） 假定汽機排氣壓力為設(shè)計工況下的值

10、，用平移設(shè)計工況熱力過程線方法初步確定排氣點。（7） 待通流部分計算結(jié)束后，校核計算軸端汽封的漏氣量，軸端汽封的漏氣量正比例于汽封前的壓力。五、N300-16.7/537/537汽輪機介紹5.1汽輪機本體和系統(tǒng)基本組成 本機是按美國西屋公司的技術(shù)制造的300MW亞臨界，中間再熱式、高中壓合缸、雙缸雙排汽、單軸凝氣式汽輪機，如圖4-2所示。它由高中壓積木塊BB0243與低壓缸積木塊BB074組合而成。為了進一步提高機組的經(jīng)濟性，對原引進技術(shù)做了改進設(shè)計，而且低壓缸末級葉片采用905mm的長葉片。機組型號為N300-16.7/537/537，工廠產(chǎn)品號為D156. 主要技術(shù)參數(shù)：額定功率300MW

11、；主汽門前額定參數(shù)16.7MPa、537C，再熱器汽門前溫度537C；工作轉(zhuǎn)速3000rpm；額定被壓5.39KPa；回?zé)峒墧?shù)3高、4低、1除氧；額定工況蒸汽流量910.2t/h、熱耗7937KJ/(KWh)。 本機通流部分共35級葉片，其中高壓缸1+11級，中壓缸9級低壓缸2*7級。動葉片除低壓缸末三級為扭轉(zhuǎn)葉片外，其余均為等截面葉片。高、中、低壓缸隔板靜葉均為扭轉(zhuǎn)葉片。有六個高壓調(diào)壓調(diào)節(jié)閥，每閥控制21個噴嘴，調(diào)節(jié)閥全開時部分進氣度為0.95455。汽輪機在額定參數(shù)下5閥全開時可以帶額定負(fù)荷。汽輪機在6個調(diào)節(jié)閥全開、新汽參數(shù)16.7MPa、537C（超壓5%），時運行，這一工況定義為最大負(fù)

12、荷工況。通流部分結(jié)構(gòu)參數(shù)如表4-5所示。最大工況和額定工況下的熱力參數(shù)如表4-6所示。 主汽門、調(diào)節(jié)閥、進氣管的壓損的4%，再熱器及管道為10%，中聯(lián)們及管道為2.5%，重低壓連通管為2%。2號高加抽汽分別來自左右側(cè)低壓缸排汽，除氧器抽器來自中亞缸排氣，第5、6級抽汽分別來自左、右側(cè)低壓缸的第3、5級前，因此低壓缸為不對稱抽汽。 靜葉與動葉汽封：高壓靜葉3齒，平齒，中壓靜葉5齒，2高3低，動葉5齒，平齒，間隙都是0.75mm；低壓靜葉3齒，1高2低，動葉1齒，平齒，間隙都是1.1mm。 氣缸抽汽汽封結(jié)構(gòu)數(shù)據(jù)見表4-45.2 汽輪機N300主要結(jié)構(gòu)參數(shù)和設(shè)計工況下的熱力參數(shù)表3-5：N300-1

13、6.7/537/537汽輪機通流部分的結(jié)構(gòu)參數(shù)缸別級號噴嘴動葉噴嘴數(shù)目平均直徑出口高度出口面積葉片數(shù)目平均直徑出口高度出口面積/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汽輪機最大工況和設(shè)計工況下的熱力參數(shù)缸別級號最大工況額定工況級前壓力MPa級前溫度內(nèi)功率kW級流量t/h級前壓力MPa級前溫度內(nèi)功率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）噴嘴設(shè)計工礦核心工況下的壓力比、，流量比 2）求取新工況下噴嘴中的各項參數(shù)： 噴嘴后壓力 根據(jù)水蒸氣求出噴嘴理想出口比焓 噴嘴中滯止理想比焓 噴嘴出口速度 噴嘴損失 噴嘴實際比焓降 出口比焓 噴嘴實際出口速度 輪周速度 ： 動葉進口相對速度 與速度對應(yīng)的比焓 =7.06 動葉進汽角 沖角 撞擊損失 動葉進口滯止焓 根據(jù)水蒸氣性質(zhì)求得動葉進口滯止參數(shù) 3） 求取動葉的壓力比，流量比，及臨界流量設(shè)計工況下，動葉進口狀態(tài)點比焓，比熵由程

22、序得：滯止壓力臨界流量 流量比 壓力比 動葉后壓力 4） 求新工況下的動葉中各項參數(shù)： 根據(jù)水蒸氣熱力性質(zhì)求得 由程序得 所以 動葉理想比焓降 動葉出口理想速度 速度系數(shù)查表得到 動葉損失 動葉出口焓 根據(jù)水蒸氣性質(zhì)動葉出口溫度 動葉實際比焓降 動葉出口相對速度 動葉出口絕對速度 動葉出汽角 5） 求級的其他各項損失、內(nèi)效率、內(nèi)功率： 無限長葉片輪周比焓降 葉高損失（將扇形損失一并考慮時， ）：噴嘴高度，單位取mm,見附表3-7。輪周比焓降：級的理想可用能量：（為下一級對本級的余速利用系數(shù)）扇形損失：動葉高度：動葉平均直徑，以上參見附表3-7動葉摩擦損失（為圓周速度）濕氣損失（為極平均蒸汽干度

23、）隔板漏氣損失：隔板汽封平均直徑：隔板汽封間隙：隔板汽封齒數(shù)：噴嘴出口面積，以上參見附表3-7葉頂漏氣損失 ：與葉頂軸向間隙和圍周帶邊厚度有關(guān)的經(jīng)驗系數(shù)：與葉頂軸向間隙及速度比有關(guān)的經(jīng)驗系數(shù):與反動度及徑高比有關(guān)的經(jīng)驗系數(shù)：葉頂軸向間隙與葉高之比值，以上參數(shù)見附表3-7，余數(shù)損失級有效焓降相對內(nèi)效率級的內(nèi)功率6.2 應(yīng)用程序第一程序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