基于MATLAB的干氣密封彈簧比壓計算分析

趙彥東， 孫 鐵， 潘 瑩, 周長茂

(遼寧石油化工大學(xué) 機械工程學(xué)院，遼寧 撫順 113001)

基于MATLAB的干氣密封彈簧比壓計算分析

趙彥東， 孫 鐵， 潘 瑩, 周長茂

(遼寧石油化工大學(xué) 機械工程學(xué)院，遼寧 撫順 113001)

根據(jù)Muijderman的螺旋槽窄槽理論，在已知螺旋槽進出口壓力的條件下，對干氣密封的彈簧比壓做了計算推導(dǎo)。借助MATLAB軟件，求解槽區(qū)氣膜壓力方程，得出了關(guān)于半徑和壓力的二維數(shù)組并擬合方程。結(jié)合經(jīng)典算例，計算了槽根壓力、開啟力，與文獻比較并分析相對誤差，驗證了推理計算的可靠性。求得工況下的彈簧比壓,為彈簧比壓的設(shè)計以及運行時的維護提供依據(jù)。

彈簧比壓； MATLAB軟件； 干氣密封

干氣密封是一種非接觸式機械密封，其中以螺旋槽密封最為典型[1-2]。在正常工況下，相同結(jié)構(gòu)的干氣密封，通入不同壓力的氮氣，對應(yīng)不同的彈簧壓力。如果壓力過低，達不到密封效果；壓力過高可能對彈簧造成破壞。因此，在已知通入氮氣壓力的情況下，求出彈簧比壓，明確彈簧是否超負(fù)荷運轉(zhuǎn),對提高密封性能，延長使用壽命[3]，以及在日常設(shè)計和操作過程中有重要的指導(dǎo)意義。

目前主要研究如何減少泄漏量，增大氣膜剛度、端面開啟力，提高密封性能。彭旭東等[4-5]引入仿生學(xué)的觀點，對螺旋槽進行設(shè)計；胡文績等[6]采用CFD的方法，對端面流場進行分析；宋鵬云等[7-8]比較泵入式干氣密封和泵出式干氣密封的開啟力，并對密封端面在非平行工況下的壓力分布進行計算。目前，對干氣密封彈簧比壓的計算，公開發(fā)表的論文還不多。

本文根據(jù)Muijderman的螺旋槽窄槽理論[9]，計算干氣密封彈簧比壓，并對實例進行計算，驗證計算方法的可行性,為干氣密封彈簧比壓設(shè)計計算提供依據(jù)。

1 計算模型

1.1 物理模型

圖1為干氣密封正常運行情況下的簡圖。假定動環(huán)和靜環(huán)的彈性模量足夠大，能夠保證端面間隙平行[10]，不考慮壓差和溫差作用使密封平行端面破壞[8]。圖1中h0為氣膜厚度，R0為螺旋槽外徑，Ri為螺旋槽內(nèi)徑，Rg為螺旋槽槽根半徑，Rb為螺旋槽平衡半徑，單位均為μm。

圖1 干氣密封平行端面

Fig.1Drygassealparallelendface

1.2 數(shù)學(xué)模型

滿足上述假設(shè)，密封端面的開啟力為：

式中，p1(r)為螺旋槽干氣密封端面壩區(qū)壓力分布；p2(r)為螺旋槽干氣密封端面開槽區(qū)壓力分布。

根據(jù)Muijderman的螺旋槽窄槽理論[9]，給出了螺旋槽干氣密封端面平行時，氣膜壓力分布的表達式。

槽區(qū)的氣膜壓力為：

壩區(qū)的氣膜壓力為：

式中，η為氣體的黏度；m為泄漏量；T為氣體的開氏溫度；Rc為氣體常數(shù)；h0為壩區(qū)氣膜厚度;ω為密封環(huán)旋轉(zhuǎn)角速度;h1為槽區(qū)氣膜厚度，h1=h0+hg，其中hg為螺旋槽深;g1(α，H，γ)，g2(α，H，γ) 為螺旋槽系數(shù)，分別為：

式中，H為槽區(qū)、壩區(qū)膜厚比;α為螺旋線螺旋角度;γ為螺旋槽臺寬與槽寬的比值。

在正常工況下，閉合力等于彈簧的彈力與內(nèi)外壓強差力的和?？汕蟮瞄]合力為：

式中，Rb為平衡半徑;psp為彈簧比壓。

2 計算結(jié)果與討論

根據(jù)實體模型以及查閱文獻，得知方程的邊界條件以及實際數(shù)值。

(1) 邊界條件：當(dāng)r=Ri時，p=pi; 當(dāng)r=R0時，p=p0；根據(jù)質(zhì)量守恒定律，流過槽區(qū)與壩區(qū)的氣體質(zhì)量相等，也就是說式(2)和式(3)中的m值相等。

(2) 查閱文獻：由文獻[11]可知：Ri=5.842×104μm，R0=7.778×104μm，Rg=6. 900×104μm，Rb=6.130×104μm，α=15°，γ=1，hg= 5 μm，pi= 1.013×105Pa，p0=4. 585×106Pa，T=303.15 K，ω=1 087. 08 rad/s，溫度303.15 K時空氣黏度η=1.86×10-5Pa·s，氣體常數(shù)Rc= 286. 706 J/(kg·K)。經(jīng)過實驗測得，當(dāng)氣膜厚度大于2.644 8 μm時，端面處于非平行運行，確定此時臨界狀態(tài)的彈簧比壓有很重要的研究意義。

2.1 壩區(qū)方程求解并討論

由式(3)積分得：

代入數(shù)據(jù)求解密封端面質(zhì)量流量為m=4.209 44×10-4kg/s。與文獻[8]比較，得出相對誤差為δ1=0.000 017%。誤差十分小，說明計算準(zhǔn)確。

將壩區(qū)方程轉(zhuǎn)化為常微分方程初值問題。應(yīng)用MATLAB軟件中的4階Runge-Kutta公式，對式(7)進行求解。設(shè)定自變量為r，因變量為p，自變量取值范圍[Ri,Rg],輸入初值p=pi, 即pi=101 300 Pa。計算結(jié)果如圖2所示，求得pg1=4 838 139 Pa。

Fig.2Sealdampressure-radiuscurve

計算得槽根氣膜壓力與文獻[11-13]對比，結(jié)果如表1所示。從表1中可以看出，本文計算的膜厚介于2.030和3.050之間，其壓力介于Gabriel、宋鵬云在相同膜厚時所計算的壓力之間，超出陳侃計算的壓力范圍。由于陳侃應(yīng)用Labview軟件，本文應(yīng)用MATLAB軟件，所以存在誤差。由表1還可以看出，應(yīng)用MATLAB計算相對精確。

表1 文中求解槽根壓力與文獻值比較Table 1 Comparison of the root pressure in the paper and the literature value

用常規(guī)數(shù)學(xué)的方法，直接對式(3)積分，得：

2.2 槽區(qū)方程求解并討論

求解式(10)時，應(yīng)用常規(guī)的數(shù)值計算方法很難得出方程的解。此處應(yīng)用MATLAB軟件，引用4階Runge-Kutta算法，求解。設(shè)定自變量為

r，因變量為p，自變量取值范圍[Rg,R0],輸入初值p=pg1, 即pg1=4 838 139 Pa。由于pg1是近似值，不能作為已知條件，以免增大誤差。因此不斷改變初值，直到當(dāng)r=7.778×104μm時，p0=4 585 200 Pa為止。此時輸入的初值為pg2=5 181 652 Pa。比較用式(8)和式(10)計算出的槽根壓力，相對誤差δ3=6.94%，誤差相對較大。所以本文取兩計算結(jié)果的平均值作為初值，即pg3=5 009 895.5 Pa。計算出槽區(qū)壓力分布，如圖3所示。

圖3 密封槽區(qū)壓-徑曲線

Fig.3Sealgroovepressure-radiuscurve

觀察圖3可知，由于計算精度高，所以許多點組成的縱坐標(biāo)p與橫坐標(biāo)r呈線性分布。因此，應(yīng)用工具箱cftool，對求得的點集進行曲線擬合，求得近似的線性關(guān)系[14]。擬合的情況見表2。

表2 方程擬合的適合度Table 2 Goodness of fit of equation

方差為0.000 351 9，相關(guān)系數(shù)為1，說明曲線擬合良好。

求得方程為：

p2(r)=p1×r5+p2×r4+p3×r3+

式(11)系數(shù)見表3。

表3 式(11)系數(shù)Table 3 Equation (11) coefficient

2.3 彈簧比壓求解及討論

應(yīng)用MATLAB積分命令，求解式(1)得：

與文獻[6]的開啟力比較，相對誤差δ4=

0.086%。

最終求得在膜厚為2.644 8μm時的彈簧比壓為：

文獻[11]中給出的端面比壓psp=1.00×104Pa，是安裝時的彈簧比壓。由于氣膜厚度在2～3μm，所以彈簧發(fā)生形變，導(dǎo)致比壓變化。本文計算在膜厚為2.644 8μm時，對應(yīng)的彈簧比壓psp=1.548×104Pa。與安裝時的彈簧比壓在同一個數(shù)量級。由式(13)可以看出分母的數(shù)量級為10-3。假設(shè)分子的計算值與實際值差1，除以分母后實際值與計算值的差值達到1 000。所以，本文計算較為準(zhǔn)確。

3 結(jié)論

應(yīng)用MATLAB軟件，建立了壩區(qū)和槽區(qū)氣膜壓力常微分方程的一種求解方法。畫出關(guān)于半徑和壓強的二維圖，從圖中可以看出壩區(qū)壓強隨著半徑的增加呈拋物線狀增加；槽區(qū)壓強隨著半徑增加呈線性減少，直至減少到通入氣體的壓強。也就是說，適當(dāng)?shù)卦黾油ㄈ霘怏w的壓強，可以提高槽根處的壓強，提高密封效果。

基于MATLAB軟件，對經(jīng)典實例進行計算分析，得出端面開啟力和槽根壓力。槽區(qū)方程和壩區(qū)方程計算出的槽根壓力誤差為6.94%。主要是由于槽、壩區(qū)方程基于螺旋槽不同結(jié)構(gòu)推導(dǎo)，導(dǎo)致誤差增大。為了減小誤差，槽根壓強應(yīng)取兩個方程計算出的平均值。

正常工況下，已知進出口壓強、螺旋槽結(jié)構(gòu)等數(shù)據(jù)，可以求出在設(shè)備運行時的彈簧比壓。與設(shè)計的彈簧比壓進行比較，可推斷是否此時處于彈簧的最佳運行狀態(tài)，是否此時已經(jīng)超負(fù)荷運轉(zhuǎn)。對彈簧的保護、密封性能的提高有一定的指導(dǎo)意義。

[1] 陳秀琴,朱維兵,王和順.干氣密封技術(shù)研究現(xiàn)狀及發(fā)展趨勢[J].液壓與氣動,2008(2): 52-55．

ChenXiuqin,ZhuWeibing,WangHeshun.Researchstatusanddevelopingtrendofdrygassealtechnique[J].ChineseHydraulics&Pneumatics,2008(2): 52-55．

[2] 郝木明.機械密封技術(shù)及應(yīng)用[M].北京:中國石化出版社, 2010:5-7.

[3] 孫見君,顧伯勤,魏龍.彈簧比壓對機械密封性能影響的分形分析[J].潤滑與密封,2006,(6):67-70.

SunJianjun,GuBoqin,WeiLong.Analysisofinfluenceofspringpressureonperformanceofmechanicalsealsbyfractaltheory[J].LubricationEngineering, 2006(6):67-70.

[4] 彭旭東,呼延晨龍,白少先,等.基于鳥翼輪廓的干式氣體密封仿生型槽設(shè)計[J].摩擦學(xué)學(xué)報,2012,32(6):563-569.

PengXudong，HuyanChenlong，BaiShaoxian，etal.Designofabiomorphicgroovedrygassealbasedonbirdwingoutlines[J]．Tribology,2012,32(6):563-569.

[5] 江錦波,彭旭東,白少先,等.仿鳥翼微列螺旋槽干氣密封性能分析與選型[J]. 摩擦學(xué)學(xué)報, 2015, 35(3):274-281.

JiangJinbo，PengXudong，BaiShaoxian，etal.Performanceanalysisandselectionofabionicbirdwingmulti-arrayspiralgroovedrygasseal[J]．Tribology，2015，35(3): 274-281．

[6] 胡文績,陳秀琴,余向東.T型槽干氣密封端面流場的數(shù)值模擬[J].潤滑與密封,2008, 33(11):20-23.

HuWenji，ChenXiuqin，YuXiangdong.NumericalsimulationoffaceflowfieldonT-shapegroovedrygasseal[J].LubricationEngineering, 2008, 33(11):20-23.

[7] 宋鵬云,丁志浩.螺旋槽泵出型干氣密封端面氣膜壓力近似解析計算[J]. 潤滑與密封,2011,36(4):1-3.

SongPengyun,DingZhihao.Anapproximateanalyticalmethodofthegasfilmpressureoftheoutwardpumpingspiralgroovedrygasseals[J].LubricationEngineering, 2011,36(4):1-3.

[8] 宋鵬云,李英,馬愛琳.螺旋槽干氣密封端面非平行間隙壓力分布的近似解析計算[J]. 排灌機械工程學(xué)報,2015,33(12):1050-1055.

SongPengyun,LiYing,MaAilin.Anapproximateanalyticalmethodofgasfilmpressureinnon-parallelgapofspiralgroovedrygasseal[J]．JournalofDrainageandIrrigationMachineryEngineering，2015,33(12):1050-1055.

[9]MuijdermanEA.Spiralgroovebearings[M].NewYork:Springer-Verlag,1966.

[10] 許靜.螺旋槽干氣密封高壓端面變形機理的理論與實驗研究[D].杭州:浙江工業(yè)大學(xué),2014．

[11]GabrielRP.Fundamentalsofspiralgroovenon-con-tactingfaceseals[J].LubricationEngineering,1994,50(3): 215-224．

[12] 陳侃，黃澤沛，姚進，等.不等溫假設(shè)下的螺旋槽干氣密封端面氣膜壓力分布近似解析計算[J].排灌機械工程學(xué)報，2013，31(7):605-610．

ChenKan，HuangZepei，YaoJin，etal．Approxi-matelyestimatingpressuredistributioningaslubricatedspiralgroovefacesealundernon-isothermalconditions[J].JournalofDrainageandIrrigationMachineryEngineering，2013，31(7):605- 610．

[13] 宋鵬云.螺旋槽干氣密封端面氣膜壓力計算方法討論[J].潤滑與密封，2009，34(7):7-9．

SongPengyun．Discussionaboutthecalculationmethodsofthegasfilmpressureofthespiralgroovedrygasseals[J].LubricationEngineering，2009，34 (7):7-9．

[14] 孫杰, 金珊. 電容式位移傳感器在氣閥測試中的應(yīng)用[J]. 石油化工高等學(xué)校學(xué)報, 2005, 18(1):55-57.

Sunjie,JinShan.Applicationofcapcity-typeshifttransduceratgasvalvetesting[J].JournalofPetrochemicalUniversities, 2005, 18(1):55-57.

Calculation Method of Dry Gas Seal Spring Pressure Based on MATLAB

Zhao Yandong， Sun Tie， Pan Ying， Zhou Changmao

(SchoolofMechanicalEngineering，LiaoningShihuaUniversity，F(xiàn)ushunLiaoning113001，China)

According to the theory of the spiral groove of the Muijderman,the MATLAB software is used to calculate the spring pressure of the dry gas seal.In this process,the gas film pressure equation is solved and a two-dimensional array about the radius and pressure is obtained. MATLAB curve fitting tool was used to fit the equation of radius and pressure.Combined with the classic example, the root pressure and the opening force of the spiral groove are calculated, and the relative error is compared with the literature, which verifies the reliability of inference calculation.The spring pressure of the working condition is obtained, which provides the basis for the design of the spring pressure and the maintenance of the spring.

Spring pressure; MATLAB software; Dry gas seal

2017-02-23

2017-03-07

中國石油化工股份有限公司科技攻關(guān)項目(311084)。

趙彥東(1992-)，男，碩士研究生，從事新型高效石化裝備設(shè)計方面研究；E-mail:2578686997@qq.com。

孫鐵(1964-)，男，碩士，教授，從事新型高效石化裝備設(shè)計方面研究；E-mail:15141366088@139.com。

1006-396X(2017)06-0092-05

投稿網(wǎng)址：http://journal.lnpu.edu.cn

TE94

A

10.3969/j.issn.1006-396X.2017.06.017

(編輯 王亞新)