1、IG-541混合氣體滅火系統設計理論和基本計算方法IG-541混合氣體滅火系統設計理論和基本計算方法一概述IG-541混合氣體滅火系統作為一種新型潔凈氣體滅火系統，由于它兼備有效滅火、綠色環保以及對人體無傷害等特性，目前已在國內外消防領域得到廣泛應用。然而，人們在大量應用它的同時，對系統性質、性能、原理等方面的量化研究卻是十分不足的。國內至今尚無完整的系統設計規范，尤其缺乏完整的系統設計計算理論和方法，甚至于連基本的單元計算方法也不齊全，現有的一些計算公式基本上照搬了國外的書本，并且缺乏完整性和系統性。這種理論研究遠遠落后于實際應用的反常現象是消防工程界特有的，也是消防系統建設與使用遠遠相脫節

2、這一客觀情況所造成的。國外公司雖有系統設計軟件可以代客計算，但并不提供計算方法，我們只能是知其然而不知其所以然。為了解決我國已有IG-541滅火系統的設備和大量實際應用，卻還沒有設計計算方法的突出矛盾，確保IG-541滅火系統設計的科學先進性、安全可靠性和經濟合理性，達到優化設計的目的，我們在努力學習和吸收國外先進技術的同時，還必須建立自己的理論研究體系和設計計算方法。本文探討了IG-541氣體滅火系統設計計算的理論依據，在此基礎上推導了和建立了IG-541滅火系統的基本計算方法，為科學地建立具有自主知識產權的IG-541滅火系統計算機設計軟件奠定了基礎。二系統特征IG-541滅火系統和其他固

3、定氣體滅火系統比較既有共性又具有鮮明的個性。IG-541在儲存條件下呈氣態，比其他滅火系統需要更大的儲存容積;在高壓下儲存和運行，管道的承壓能力要求亦較高，設備投資費用大，精確計算和優化設計可以帶來明顯的經濟效益。IG-541滅火的有效濃度為>37.5%而對人體安全的濃度為<42.8%，同時滿足以上條件必須嚴格控制儲存量，并且對于防護區域有相應要求。IG-541滅火系統的使用條件要求，系統開啟后，90%藥劑噴放時間應>23秒及<40秒，并且又要求60秒鐘內達到滅火濃度。這也是一個相當嚴格的的設計約束條件。IG-541滅火系統和其他滅火系統相比,滅火劑設計濃度以及噴射速率

4、的容差范圍小得多，且與平常容易發生的誤解不同,寬裕的設計不僅浪費投資，設計結果也未必安全。因此，系統設計應采用精確的、全過程動態模擬的分時計算方法。單一氣相下發生，可以通過嚴格的方法，借助電子計算機進行精確的計算。IG-541滅火系統設計的主要目標是要保證在裝置啟動后的指定時間內，防護區中的滅火劑達到設計濃度，其中計算IG-541氣體在系統各單元中的流動推動力和阻力是關鍵，二者又取決于系統的物性和單元的設備特征，茲在下文逐一討論。三.純組份性質物質的物理化學性質甚多，這里只討論和IG-541系統設計有關的P-V-T性質、熱力學性質和遷移性質。1.基礎物性氮氣、氫氣和二氧化碳皆為常見氣體，其有關

5、性質可以從手冊中查到。茲參照ChemicalPropertiesHandbook(1999)一書,將相關數據羅列如下。其他資料上除分子量以外的數據并不完全相同，但對本過程的設計計算結果影響甚微。表一.基礎物性和熱力學性質名稱符號單位ArC02n2英文名ArgonCarbondioxideNitrogen分子量Mg/mol39.94844.01028.013常冰點TfK83.80216.5863.15常沸點TbK87.28194.7077.35臨界溫度TcK150.86304.19126.10臨界壓力PcBar48.9873.8233.94臨界比容Vccm3/moI74.694.090.1臨界密

6、度pcG/cm30.53560.46820.3109臨界壓縮因Zc無因次0.2910.2740.292IG-541設計計算的有利條件是：物系臨界溫度低，整個過程在偏心因子3無因次0.0000.2280.040偶極距Debye0002.熱性質流體在理想氣體狀態下的熱性質是計算熱力學性質的基礎。通常表達為熱容或給的多項式。如：hi0=C0,i+Ci,iT+C2,iT2+C3,iT3+-（式1）這些多項式系數通常是用于相當寬的溫度范圍，而IG541的工作溫度范圍較窄，約在200至320K之間，可以專門回歸較為簡短和精確的多項式。茲將APIProject44SelectedVolumeofProper

7、tiesofChemicalCompounds所列CO2N2的恒壓熱容G°文獻值列于表二。Ar的G°受溫度影響極小，用ThePropertiesofGasesandLiquids»3rdEd一書附錄中的多項式求得。表二.理想氣體狀態下的恒壓熱容Cal/mol/K溫度，K150200273.15298.15300400Ar4.96814.96804.96784.96784.96784.9679CO7.2287.7338.5948.8748.8949.8766.9566.9576.9596.9616.9616.9913.低壓下的氣相粘度在低壓下氣體的粘度和壓力關系不大

8、，可以視為僅僅是溫度的函經與HandbookofChemistryandPhysics»,80thEd(1999-2000),«ChemicalPropertiesHandbook»(1999)的資料上的數據比較，發現用多項式擬合低壓氣體的粘度的精確度不高。而用Lennard-Jones12-6分子勢能位計算，25C下的誤差降至0.1%左右。下表中Ar的數據取自美國石油學會API手冊，其余數據取自ThePropertiesofGasesandLiquids»,后者Ar的參數計算結果誤差較大。表三.遷移性質一氣相粘度L-J分子勢能位ArCON2勢能參數/k

9、,K124.9195.271.4碰撞半徑b,?3.42337.1329.124Lennard-Jones12-6分子勢能位計算低壓氣體的粘度的公式用：刀gas=26.69(MT)1/2/高Qv(式2)式中：刀gas氣體粘度,尸Qv=1.16145/T*0.14874+0.52487/exp(0.7732T*)+2.16178/exp(2.43787T*)T*=T/(/k)四.混合物性質IG-541混合氣體的配方是公開的，即52%(mol)的氮氣、40%(mol)的敏氣和8%(mol)的二氧化碳氣體。混合物物性的詳盡實驗數據很少。混合物的性質主要通過“混合規則”計算求得，而這些規則則是經過若干離

10、散的實驗點來驗證的。令混合物中i組份的含量為Xi分子分率。下標i,j,k為組份序號，下標m表示混合物。則有：1. 1.分子量Mm=SXiMi(式3)2. 2.混合物臨界參數混合物的臨界參數并不是混合物的真臨界性質，而是用于計算混合物P-V-T性質和熱力學函數用的參數。采用Lee-Kesler狀態方程的混合規則：Vci=ZciRTci/Pci（式4）Zci=0.2905-0.085wi（式5）Vcm=0.1252溝Xk（Vcj1/3+Vck1/3）（式6）Tcm0.1252溝Xk（Vcj1/3+Vck"3）3（Tcj+Tck）2/Vcm（式7）3m=ExiGJi（式8）Pcm=（0.2

11、905-0.085wm）RTcm/Vcm（式93. 3.IG-541在理想氣體狀態下的熱性質系數理想氣體狀態下的熱性質系數符合加和性，從表二的數據經過計算和單位換算得到IG541在理想氣體狀態下始的熱性質系數為：C0=8.8（此項系數用來校準150400K的始值，計算偏差<0.03%）C=25.648C2=0.0005042G=1.715E-6C4=81.04（此項系數取，1psia壓力下的理想氣體單質為零嫡）用以上熱性質系數計算IG541在理想氣體狀態下的熱容、始和嫡的公式是：Cp0=C1+2C2T+3c3T（式10）h0=C0+C1T+C2T2+C3T3（式11）s0=C1lnT+2

12、C2T+3C3T2/2+C4-RlnP（式12上述公式中的單位是：mol,K,J,Pa,上角標0表示理想氣體狀態。4. 4.氣體混合物的粘度氣體混合物的粘度用美國石油學會APIproject44推薦的方法計算：xijXi其中ij為充間參數0.50.25彳iMjij03.81業Mj（式14將以上公式求得的IG-541粘度再回歸成溫度的多項式：nm0=178.86+0.5123t-0.00039t2（式15）（式15）中刀苗是IG-541在低壓下的粘度，t攝氏溫度。在-10至50c范圍內回歸誤差小于0.1%。5.5.低壓下的氣相導熱系數純組分的導熱系數由Lang'sChemicalHand

13、book»15thEd（1999）查得，如下表所示:表四.氣相導熱系數J/s/m/K溫度，c02040Ar0.01660.01760.0186CO0.01440.01600.0176N0.02410.02560.0270IG-5410.02020.02160.0228氣體混合物的氣相導熱系數按Wassiljewa方程計算：m=2（yii/2VjAij）（式16）式中：m氣相混合物的導熱系數組分i的氣相導熱系數Aj組分之間的充間作用參數，用LindsayandBromley方Aij=0.251+(qi/qj)(Mj/Mi)0.75(T+S)/(T+Sj)0.52(T+Sj)/(T+Si

14、)其中：Sutherland常數S=1.5Tbi;Sj=Cs(SS)0.5非極性氣體Cs=1用MasonandSaxena方法以及忽略充間作用參數取得的結果和上述方法偏差不大于1%按以上方法求得IG-541在0、20、40c低壓下的氣相導熱系數,回歸成多項式:（式17)式中:0=0.020238+6.75910-5t-6.46810-8t20IG-541在低壓下的氣相導熱系數，J/s/m/K攝氏溫度C為便于比較和利用現將IG541混合氣體物理性質及相關計算公式匯總如下：IG541混合氣體的物理性質ThePROPERTIESofIG541MIXTRUE名稱符號單位st氣二氧化碳氮氣IG541英文

15、名SymbolunitARGONCarbo)nDioxideNitrogenInergenArCO2N2MWg/mol39.94844.01028.01334.067點TbK87.28194.777.35111TcK150.86304.19126.1147.7PcVcZcOM/K25C5,Bar48.9873.8233.9442.39cc/mol74.694.090.183.2J/mol/KJ/mol/K0.2910.0003.423124.920.78512.4701.6667224.420.2740.2920.2870.2230.0400.0393.9413.798195.271.437.

16、12928.8141.289150.529.12520.8101.400175.5226.42918.1151.459191.42距1G541分Debye0000YMol%40852100計算方法和依據另詳。注：正體數值是文獻值,斜體數值是計算值。IG541的密度:標準狀態（0C,1atm）1.521kg/m3,對空氣比重1.176存儲狀態（20C,15Mpa）233.6kg/m3IG541在低壓下的粘度多項式：=178.86+0.5123A-0.00039X2以為絕對粘度下；t為攝氏溫度。上式適用于-10至50cIG541在低壓下的導熱系數多項式：=202.38+0.6759t-6.468X

17、10-4t2為導熱系數,W/cm/K;t為攝氏溫度。IG541在理想氣體狀態下的始多項式：H0=25.648T+5.042X10-4T2+1.715M0-6T3上式中：H0為始,J/mol;T為絕對溫度,T=t+273.15五.壓力下的性質1. P-V-T關系表征流體壓力、比容和溫度（P-V-T）關系的方程叫著狀態方程。最簡單,也是最古老的狀態方程是理想氣體狀態方程：PV=RT式中：P-壓力PaV-比容m3/molR-氣體通用常數=8.31441T-溫度K理想氣體狀態方程對常壓常溫或高溫的“永久氣體”可以使用，但是用來計算低溫高壓氣體，或是非永久氣體誤差就比較大，對于液體則完全不適用，于是，近

18、百年來許多科學家提出了對理想氣體定律的修正，這類方程式稱為真實氣體狀態方程。真實氣體狀態方程的通式是：PV=zRT式中z的叫做壓縮因子，表示真實氣體和理想氣體的差異。文獻上發表的真實氣體狀態方程有二百多個。其中最準確的是Lee-Kesler方程，它的優點體現在：1 .在很寬的溫度壓力范圍內有良好的準確度；2. 既能用于氣相也能用于液相；3. 對于任何非極性流體，只需要臨界溫度、臨界壓力和偏心因子三個參數便可以求解，而這三參數都是宏觀可測量的物性；4. 對于混合物利用混合規則，可以像純物質一樣，用同樣的方程式求解；5. 對于混合物，不僅能求解總性質，還能解得偏性質，即混合物中各個組份的性質。在L

19、ee-Kesler方程發表前后，還有一些優秀的狀態方程。這些方程求解比較簡單，在一定的場合可以滿足工程計算需要，但就準確度而言，還是及不上Lee-Kesler方程。狀態方程不僅可以用來求解流體的密度，還可以和物質的熱性質系數結合起來，求得在真實流體（液體和壓力下的氣體）狀態下的熱力學性質，傳熱性質等等，用于流體力學、化工熱力學、工程熱力學的計算。對于滅火系統而言，不僅管道輸送，而且孔板節流、噴嘴噴射、氣體膨脹和壓縮的精確計算度需要用到狀態方程，因此我們選用了Lee-Kesler方程設計計算IG541混合氣體滅火系統。Lee-Kesler狀態方程為多參數方程，其形式是：z=Pr/Vr/Tr=1+

20、B/Vr+C/Vr2+D/Vr5+c4（B+丫/Vr2）exp（-丫/Vr2）/Tr5/Vr2（式18）其中：B=bl-b2/Tr-b3/Tr2-b4/Tr3（式18-1）C=cl-c2/Tr+c3/Tr2（式18-2）D=d1+d2/Tr（式18-3）Vr理想對比比容=Pc*V/n/R/Tc，其中V/n稱為比容Pr對比壓力=P/PcTr對比溫度=T/TrLee-Kesler狀態方程是三參數普遍化方程，只要知道物質的臨界壓力、臨界溫度和偏心因子就能求解。但是方程中的壓縮因子z是隱函數，要用試差法求解。表四.Lee-Kesler狀態方程系數系數簡單流體參比流體系數簡單流體參比流體b10.1181

21、1930.2026579c30.00.016901b20.2657280.331511c40.0427240.041577b30.1547900.027655_4d1M00.1554880.48736b40.0303230.203488d2M040.6236890.0740336c10.02367440.0313385B0.653921.226c20.01869840.0503618丫0.0601670.03754Lee-Kesler狀態方程有多達七個實根，其中只有一個或兩個實根有物理意義，要用試差法求解。本文作者曾對此作過研究，具體方法可參考炫類物理化學手冊一書。IG541氣相壓縮因子Com

22、pressibilityfactorofIG541(Z)溫度-1001020304050壓力MPaTr1.7821.8491.9171.9852.0522.1202.18810.2360.98960.99120.99260.99380.99490.99580.996620.4720.97960.98280.98570.98810.99020.99200.993630.7080.97010.97500.97920.98280.98590.98870.991040.9440.96110.96760.97320.97800.98210.98570.988951.1790.95270.96080.96

23、780.97370.97880.98320.987161.4150.94500.95470.96290.96990.97600.98120.985771.6510.93800.94920.95860.96670.97360.97960.984881.8870.93180.94430.95490.96400.97180.97850.984392.1230.92640.94020.95190.96190.97050.97790.9843102.3590.92190.93690.94960.96040.96970.9/0.9847112.5950.91830.93430.94790.95950.96

24、950.97810.9855122.8310.91560.93260.94700.95920.96980.97890.9868133.0670.91390.93160.94670.95960.97060.98020.9884143.3020.91310.93150.94710.96050.97200.98190.9905153.5380.91330.93210.94820.96200.97390.98410.9930163.7740.91440.93360.95000.96410.97630.98680.9959Lee-Kesler狀態方程要對三個流體求解：一個是需要計算的流體，即IG541,

25、公式中用上角標表示；一個是偏心因子為0的“簡單流體”，上角標用°表示；還有一個是偏心因子為0.3978的正辛烷，叫做“參比流體”，上角標（r）o先在對象流體的操作條件求得對比壓力Pr和對比溫度TrPr=P/Pcm（式Tr=T/Tcm20）從（式18）分別求得簡單流體和參比流體的Z°和Z（然后用以下公式計算對象流體的壓縮因子Zi=Z0+Zr-Z（0）3/Jr）（式21）（式20）中3為對象流體即IG541偏心因子3m,3C為參比流體的偏心因子=0.3978。（式21）的物理意義是：流體的壓縮因子等于同對比溫度壓力下簡單流體的壓縮因子，加上和偏心因子相關的修正值。偏心因子在微觀

26、上反映物質分子的大小和形狀，小的球形分子如鼠氣，偏心因子等于0;宏觀上偏心因子通過物質的對比飽和蒸汽壓來計算：3=-lg（PSTr=0.7/PC）-1（式21-1）按照（式21-1）定義（式（式-J） Z （r）- Z 0 （式21-4）（式21-1）中lg是十進對數，PSTr=0.7是對比溫度為0.7時的飽和蒸汽壓。據此，作者曾在70年代進行過驗算，由表四的系數，通過Lee-Kesler方程計算出“簡單流體”和“參比流體”氣液相逸度，求得對比溫度0.7下的相平衡和飽和蒸汽壓計算簡單流體和參比流體的偏心因子，分別為：簡單流體3°=0.0052921-2）參比流體3"）=0.

27、3954721-3因此，作者將（式21）改成：Z0）=Z0+（3-J）/（3以下是采用上述方法進行的鼠氣、二氧化碳和氮氣的壓縮因子驗算。此項工作的目的一是驗算Lee-Kesler方程對于IG541組份的準確性；二是比較（式21）、（式21-4）何者更符合實際情況。一般在近臨界區或高壓下P-V-T計算誤差大，為此本文選擇IG-541實際運行范圍內的溫度和較高的壓力。驗算結果列于下列表格。表五.瀛氣壓縮因子驗算操作條件文獻值（式21）（式21）（式21-4）（式21-4）壓力溫度比容壓縮因子比容洱苦比容洱苦BarKcc/g無因次cc/g%cc/g%1002002.96250.71172.9616-

28、0.03072.9571-061920.88774.63450.33144.62970.22691003005.96330.95515.98150.34095.97640.21971003507.19630.98797.20910.17757.20380.10411502001.81660.65461.82510.46971.82190.28981502502.98340.96012.99790.48702.99350.33831503003.94630.94803.96730.53193.96260.41181503504.81430.91134.82380.384

29、84.82790.2820平誤均差0.33970.2569表六.二氧化碳壓縮因子驗算操作條件文獻值（式21）（式21）（式21-4）（式21-4）溫度壓力比容壓縮因子比容誤差比容誤差KBarcc/mo無因cc/mo%cc/mo%l次ll225.9962744.0.8992772.1.01522772.1.023608779238.9351882.0.8791883.0.02881883.0.0039087352412.831318.0.8471315.-0.2291315.-0.216008988982517.85941.90.809938.7-0.343938.8-0.32706521769

30、2624.19682.50.763679.3-0.466679.5-0.447046978092732.03498.00.710495.0-0.597495.1-0.574042753672841.59361.80.646359.1-0.756359.2-0.729054504042953.15256.20.564253.3-1.122253.4-1.089025973653067.09164.00.441158.6-3.295158.6-3.25005010570平均誤差0.87290.8555表七.氮氣壓縮因子驗算操作條件義獻值（式21）（式21）（式21-4）（式21-4）壓力溫度壓縮因

31、子壓縮因子洱苦反左壓縮因子誤差BarK無因次無因次%無因次%1002000.844990.845010.00190.84392-0.12721002500.958560.960200.17130.959250.07241003001.00511.006210.11001.005380.02751003501.02691.027120.02191.02639-0.05001502000.849860.850190.03920.84897-0.10471502500.971340.973880.26170.97265002561.028240.25751.027120.1

32、4871503501.05111.052990.17961.051980.0835平均誤差0.13040.0936表五、六、七的第列，是國際純和應用化學學會（IUPAC）資料上提供的數據；表五、六的第列數值，由第列換算而來。這些數據表明：除了二氧化碳外，Lee-Kesler方程的計算結果皆與文獻值十分符合；而本文作者建議的偏心因子校正式（式21-4）,比原作者提出的（式21）更準確。表六所列的溫度范圍，在二氧化碳的常冰點和臨界點之間，壓力為相應溫度下的飽和蒸汽壓，無論對實驗還是計算都是高誤差區。二氧化碳還有一個特殊情況：在對比溫度0.7時已經是固相，（式21-1）的引用發生了問題，因而不同文獻

33、上的偏心因子值有差異。盡管物質固相升華壓和液相的蒸汽壓函數連續，導數卻不連續。相比較之下用液相的蒸汽壓外推計算偏心因子要合理一點。經驗證，ThePropertiesofGasesandLiquids»3rdEd一書上的Harlacher蒸汽壓方程，從二氧化碳的冰點到臨界點度很準確，由此算得的偏心因子為0.223。2 .熱力學差值函數熱力學差值函數的定義是：真實流體的熱力學性質和同溫度下理想氣體的熱力學函數之差。熱力學差值函數用無因次數的形式表示，(式 22)通式為：g+gs°）/R,（h-h°）/RTc。由（式18）求得壓縮因子代入Lee-Kesler狀態方程的差

34、值函數表達式，求得g、g°,經（式21）求得gi然后結合（式10）至（式12）即可求得真實流體在指定溫度壓力下的熱力學函數。例如和g相應的恒壓熱容、嫡、始的差值函數為:(Cp-Cp0) /R, (s -h=RTc（h-h°）/RTc+h0=RTc（hh°）/RTc+C1T+C2T2+C3T3+CT4（式23）IG541在不同壓力下的始EnthalpyatpressureofIG541(J/mol)理想氣體：溫度-1001020304050ARGON1.g.5472.05679.95887.76095.66303.46511.36719.1CO21.g.8093.6

35、8451.18813.09179.39550.19925.410305.4N21.g.7651.17942.08233.08524.28815.49106.89398.3IG541*1._g_6815.47078.37341.67605.27869.38133.88398.7壓力下氣體:溫度-1001020304050壓力MPaTr1.7821.8491.9171.9852.0522.1202.18810.2366714.66984.77254.47524.07793.58062.98332.320.4726613.36890.97167.47443.17718.17992.68266.630

36、.7086511.76797.27080.77362.77643.47922.98201.740.9446410.26703.86994.57282.97569.37854.28137.651.1796308.96610.96909.07203.97496.27786.38074.561.4156208.16518.76824.37125.97424.17719.58012.571.6516108.26427.56740.77049.07353.17653.87951.681.8876009.56337.66658.46973.47283.47589.47892.092.1235912.562

37、49.36577.76899.27215.27526.47833.7102.3595817.46162.86498.66826.87148.57464.9/6.8112.5955724.76078.46421.66756.17083.57404.97721.4122.8315634.75996.56346.76687.47020.37346.77667.6133.0675547.85917.26274.16620.86959.07290.27615.4143.3025464.25840.86204.06556.46899.77235.57564.8153.5385384.25767.36136

38、.66494.36842.57182.67516.0163.7745307.95697.16071.96434.66787.47131.77468.9注:*=25.648*T+0.0005042*TA2+1.175e-6*丁八3IG541在不同壓力下的內能InternalEnergypressureofIG541J/mol理想氣體:溫度-1001020304050ARGONi.g.3284.13408.83533.53658.23782.93907.64032.3CO2i.g.5905.66180.16458.86742.07029.67321.87618.6N2i.g.5463.15670.

39、95878.86086.86294.96503.16711.5IG541i.g.4627.54807.24987.45167.95348.85530.15711.9壓力下氣體:溫度-1001020304050壓力MPaTr1.7821.8491.9171.9852.0522.1202.18810.2364549.44733.64917.65101.75285.85470.15654.620.4724470.04658.84847.05034.85222.35409.75596.930.7084389.34583.04775.64967.35158.35348.85538.940.9444307

40、.34506.24703.44899.25093.95287.75480.751.1794224.34428.74630.74830.75029.15226.35422.461.4154140.44350.54557.54761.94964.25164.85364.071.6514055.84271.94484.04692.84899.15103.25305.681.8873970.84192.94410.34623.74834.05041.75247.292.1233885.54113.94336.64554.74769.04980.35189.1102.3593800.44035.0426

41、3.14485.84704.34919.25131.1112.5953715.63956.54190.04417.44639.94858.45073.6122.8313631.53878.64117.34349.44575.94798.05016.4133.0673548.33801.54045.44282.04512.54738.14959.6143.3023466.43725.43974.44215.34449.84678.94903.5153.5383386.03650.43904.34149.64387.84620.34847.9163.7743307.23576.93835.4408

42、4.84326.74562.44793.0IG541在不同壓力下的嫡EntropyatpressureofIG541(J/mol/K)溫度-1001020304050壓1.71.81.91.92.02.12.1力MPaTri824917855220881ps.00150151152153154155156ia158.93.91.85.77.65.51.340.0.0128129130131132133134124.69.67.62.53.42.28.110.0.01221231241251261271282471.93r.91.85.77.65.51.350.30.071119.56120.5

43、4121.48122.40123.28124.14124.970.40.094117.16118.14119.09120.01120.89121.75122.580.50.118115.31116.29117.24118.15119.04119.89120.730.60.142113.50114.50115.47116.41117.31118.19119.030.70.165112.21113.22114.19115.13116.03116.91117.750.80.189111.10112.11113.08114.02114.92115.80116.640.90.212110.13111.13112.10113.04113.94114.82115.6610.236109.25110.26111.23112.16113.07113.94114.7920.472103.19104.22105.22106.17107.10107.99108.8530.70899.51100.58101.60102.58103.52104.43105.3040.94