張偉政， 李水平， 丁雪興， 席喜林

(蘭州理工大學 石油化工學院，甘肅 蘭州 730050)

靜壓干氣密封微尺度氣膜力學特性

張偉政， 李水平， 丁雪興， 席喜林

(蘭州理工大學 石油化工學院，甘肅 蘭州 730050)

根據(jù)靜壓干氣密封的特點，以氣體雷諾方程為依據(jù)，通過伽遼金法得到潤滑氣膜壓力分布的變分方程，由氣膜的邊界條件，得到用有限元法求解穩(wěn)態(tài)下氣膜雷諾方程。并利用多項式擬合得到不同工況下的壓力擬合曲線，計算得到開啟力F。由不同膜厚下開啟力的大小擬合出開啟力與膜厚的函數(shù)關(guān)系式，及對膜厚求導得到氣膜剛度與氣膜厚度的關(guān)系式。結(jié)果表明，在靜環(huán)密封端面徑向上壓力，由節(jié)流孔處向內(nèi)徑和外徑處呈現(xiàn)拋物線型遞減趨勢，在節(jié)流孔處出現(xiàn)壓力最高峰；隨著氣膜厚度的增加，開啟力逐漸減小，氣膜剛度也逐漸減小。

數(shù)值模擬； 氣膜剛度； 穩(wěn)態(tài); 伽遼金法； 微尺度

目前在高轉(zhuǎn)速設(shè)備密封中，干氣密封是發(fā)展比較迅速的新型密封類型。為了使干氣密封的應用更加廣泛，適用于更多工程應用，提出靜壓式氣膜密封，用于低轉(zhuǎn)速設(shè)備中[1-2]。在低轉(zhuǎn)速下，不能形成足夠的氣膜開啟力推開密封端面，進而高速動壓效應理論已經(jīng)不能適用，這就需要建立一套適用低轉(zhuǎn)速運行下的干氣密封潤滑理論來指導工程方面的應用[3-4]。

干氣密封主要依靠密封軸的轉(zhuǎn)速，在高速轉(zhuǎn)動下形成動壓效應，從而使動靜環(huán)分開。然而在低轉(zhuǎn)速設(shè)備中，如反應釜，不能形成穩(wěn)定工作下所需要的氣膜剛度。靜壓式干氣密封適用于低轉(zhuǎn)速設(shè)備軸端密封，其穩(wěn)定運行的條件之一是有穩(wěn)定的氣膜剛度?？紤]到靜壓干氣密封結(jié)構(gòu)，靜環(huán)端面的開槽深度及在端面的位置、靜環(huán)端面節(jié)流孔的特性參數(shù)、密封間隙等都成為影響穩(wěn)定運行的因素。供氣壓力、出口壓力、節(jié)流孔直徑及個數(shù)等節(jié)流孔特性參數(shù)是靜壓干氣密封能否穩(wěn)定運行的關(guān)鍵，也是提高密封使用壽命的關(guān)鍵。

為使計算結(jié)果更符合實際工況，將密封端面的節(jié)流孔特性考慮進去，這增加了求解方程的難度。而解決這一難題的方法主要是通過有限元法[5-6]或者差分法[7-8]來求解加入氣源后的雷諾方程。杜建軍等[9]利用有限元法，計算出氣膜在軸承內(nèi)部的壓力分布方程，并分析出氣膜壓力在軸承內(nèi)部的變化規(guī)律。金朝旭等[10]根據(jù)氣體潤滑理論，使用有限元法求解靜壓下干氣密封端面的雷諾方程，討論分析了不同轉(zhuǎn)速和不同靜壓結(jié)構(gòu)下，氣膜剛度等因素對密封性能的影響規(guī)律。余建平[11]利用有限差分法來求解氣膜密封的雷諾方程，研究分析了不同結(jié)構(gòu)參數(shù)和不同的操作參數(shù)對靜壓干密封端面性能的影響。

綜上所述，帶節(jié)流孔修正后的二階偏微分雷諾方程的求解難度大，有限差分法的收斂差，計算效率低。有限元法在氣體軸承的使用，結(jié)合精確計算軸承內(nèi)氣體雷諾方程，采用伽遼金法將雷諾方程進行變分化簡。利用有限元法設(shè)置壓力分布的插值函數(shù)，得到離散化的雷諾方程，再通過MATLAB編程，得到外加壓式靜壓干氣密封端面壓力分布的計算程序，根據(jù)不同工況下運行情況，分析氣膜剛度的穩(wěn)定及其變化規(guī)律。通過上述理論和方程計算，為靜壓干氣密封在穩(wěn)定下運行提供參考，也為生產(chǎn)出良好的產(chǎn)品，生產(chǎn)技術(shù)的改進提供理論依據(jù)。

1 工作原理

在低轉(zhuǎn)速運行下的干氣密封的主要基礎(chǔ)理論是靜壓氣體潤滑軸承理論[12]，通過外界提供的氣體壓力，在密封端面間形成穩(wěn)定可靠的且使動靜環(huán)分開的力，從而實現(xiàn)端面間的氣膜密封。因氣源供氣方式的不同，將靜壓干氣密封分為自加壓式和外加壓式兩類。靜壓干氣密封的外加壓式典型結(jié)構(gòu)如圖1所示[13]。

圖1 外加壓式干氣密封

Fig.1Externalpressurizedaerostaticdrygasseal

靜環(huán)端面的結(jié)構(gòu)如圖2所示。在靜環(huán)的端面開有節(jié)流孔，背部用作氣源氣體流通的通道，通過外界供氣，在端面形成穩(wěn)定壓力。因為有獨立的供氣系統(tǒng)，端面間的壓力比較穩(wěn)定，更加可靠。所以，選擇對外加壓式干氣密封的結(jié)構(gòu)進行理論分析。

圖2 靜環(huán)端面

Fig.2Staticringsealface

2 理論模型

2.1 力學模型

對外加壓式靜壓干氣密封進行受力分析(見圖3)，當靜壓干氣密封系統(tǒng)處于穩(wěn)態(tài)運行時，端面間壓力的變化引起氣膜厚度的變化、開啟力的變化?；跉饽ず穸扰c開啟力之間關(guān)系，建立靜壓干氣密封系統(tǒng)在穩(wěn)態(tài)工況條件下的力學模型。

圖3 靜壓干氣密封受力示意圖

Fig.3Theforcediagramofaerostaticdrygasseal

2.2 數(shù)學方程

2.2.1 穩(wěn)態(tài)下的雷諾方程 利用靜壓氣體軸承潤滑理論和流體力學等相關(guān)理論[14]，考慮節(jié)流孔氣體流量的影響，推導出適用于外加壓式靜壓干氣密封在穩(wěn)定運行下的雷諾方程：

式中，η為氣體動力黏度;h為氣膜厚度;p為氣膜壓力;r(θ)為密封端面的極坐標;δi為Kronecker數(shù)，在無孔處為0，在有孔處為1。Q為節(jié)流孔處的流量，Q=24pv[12],v為節(jié)流孔氣體流速。

邊界條件為:p|r=r1=pa，p|r=r2=pn。r1、r2分別為密封環(huán)端面內(nèi)徑和外徑。

2.2.2 變分方程的建立 利用伽遼金方法[15]將上述所得到的雷諾方程進行化簡，得到無量綱的伽遼金方程：

式中，ε為方程余量，δf為氣膜無量綱壓方的變分，Ω為解域。

將式(4)代入到伽遼金方程得：

利用分步積分和高斯公式來簡化解域為Ω的變分方程，可得：

式(6)即為有限元求解的出發(fā)點。

再取變換式：

在式(7)的變換之下，圖4(b) 的扇形將變換為圖4(c)所示的矩形，相應的網(wǎng)格亦成為矩形，進一步劃分為三角形的有限單元。

圖4 計算區(qū)域

Fig.4Calculationarea

2.2.4 方程的離散 通過對Ω解域進行劃分，形成M個三角形單元后，變分方程(6)對整個解域的積分就可以簡化為M個單元，面積為Δe的積分之和。即為：

解域內(nèi)任意點的氣膜壓力可以由解域內(nèi)劃分節(jié)點的壓力來表示。通過線性插值函數(shù)，得到在每個單元上的剛度矩陣，再組合起來寫成矩陣形式，其形式為：

最后，通過超松弛迭代法進行求解，即可求出壓力分布。

2.2.5 端面開啟力和氣膜剛度特性的計算

(1) 開啟力

密封處于穩(wěn)定運行時，開啟力等于閉合力：

根據(jù)有限元法算得端面壓力分布，通過二次多項式擬合得到壓力函數(shù)p(r)，進而得到開啟力，即

(2)氣膜剛度

根據(jù)考慮的方向不同分別有軸向氣膜剛度和角向氣膜剛度。本文以軸向氣膜剛度為例分析，即：

3 分析與討論

3.1 算例

圖5 靜環(huán)端面壓力分布

Fig.5Pressuredistributionofstaticringsealface

從圖5中可知，在靜環(huán)端面徑向上，壓力由節(jié)流孔處向內(nèi)徑、外徑呈拋物線型遞減趨勢；在周向上，由節(jié)流孔處逐漸遞減至兩節(jié)流孔中間。在節(jié)流孔處壓力達到最大值為0.352 182 MPa。在兩節(jié)流孔中間內(nèi)徑位置處取最小值為0.181 328 MPa。

3.2 不同工況下開啟力分布

在不同工況計算得到端面壓力分布函數(shù)，通過開啟力公式計算得到圖6曲線。從圖6中可以看出，不同工況下隨著氣膜厚度的增加，節(jié)流孔徑的增加，開啟力逐漸減小。分析得出，氣體流量一定時，氣膜厚度增大，引起開啟力減?。粴饽ず穸纫欢ǖ那闆r下，由節(jié)流孔流量計算公式可知，隨著節(jié)流孔徑的增大，氣體流量增加，使端面開啟力增加。因此，隨氣膜厚度逐步增加，開啟力呈現(xiàn)遞減的趨勢。

圖6 不同工況下開啟力分布

Fig.6Openforcedistributionofdifferentconditions

3.3 不同工況下氣膜剛度分布

不同工況下氣膜剛度分布如圖7所示。從圖7中可以看出，靜壓干氣密封處于穩(wěn)態(tài)運行時，隨著氣膜厚度的增加，氣膜的剛度呈現(xiàn)遞減趨勢。隨著節(jié)流孔的增加，氣膜剛度在增加。節(jié)流孔徑較小時，隨著氣膜厚度的增加，氣膜剛度遞減的速度較快。為使靜壓干氣密封在運行時起到更好的密封效果，應使穩(wěn)定運行下的氣膜厚度達到最小值?？紤]到較小的節(jié)流孔不易加工，應使節(jié)流孔直徑增大一些，得到足夠的氣膜剛度。

圖7 不同工況下氣膜剛度分布

Fig.7Filmstiffnessdistributionofdifferentconditions

4 結(jié)論

(1) 針對靜壓干氣密封系統(tǒng)，建立了系統(tǒng)的分析模型，給出了建立變分方程的方法和過程，進一步給出了基于有限元數(shù)值方法求解氣膜壓力分布的方法和過程，通過編寫程序算出了端面氣膜壓力分布圖。

(2) 在靜環(huán)端面徑向上，由節(jié)流孔出來的氣體壓力向內(nèi)徑和外徑逐漸遞減，節(jié)流孔處的壓力為最大值。在2 μm氣膜厚度的工況下，壓力最大值為0.352 182 MPa，最小值為0.181 328 MPa。

(3) 徑向壓力由節(jié)流孔向內(nèi)徑和外徑方向呈近似拋物線狀下降，壓力最大值出現(xiàn)在節(jié)流孔處，在氣膜厚度為2 μm時。

(4) 根據(jù)壓力圖的分布趨勢，利用二次多項式擬合，得到不同工況下不同膜厚端面壓力擬合函數(shù)關(guān)系式與函數(shù)圖，分析不同氣膜厚度對開啟力與氣膜剛度的影響。

[1] 劉飛.靜壓干氣密封研究[D].成都:西華大學,2010.

[2] Stolarski T A, Xue Y. Performance study of a back-depression mechanical dry gas seal [J]. Journal of Engineering Tribology, 1998, 212(4): 279-290.

[3] Mohamid F, Yong T, Marc B. Prediction of the stability of air thrust bearings by numerical, analytical and experimental methods [J]. Wear, 1996, 198(1-2):1-6．

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

Chen Xiuqin,Zhu Weibing,Wang Heshun. Research status and developing trend of dry gas seal technique.[J]. Chinese Hydraulics & Pneumatics, 2008(2): 52-56.

[5] 丁雪興,王文鼎,金海俊,等. 干氣密封摩擦副啟停階段摩擦特性的仿真研究[J]. 石油化工高等學校學報, 2017, 30(2): 91-96.

Ding Xuexing,Wang Wending,Jin Haijun,et al.A numerical simulation study of frictional characteristics of seal faces on the starting and stopping process of a dry gas seal[J]. Journal of Petrochemical Universities, 2017, 30(2): 91-96.

[6] 周圣人,朱維兵,王和順.節(jié)流孔特性對靜壓干氣密封性能的影響[J].潤滑與密封,2010,35(11):69-73.

Zhou Shengren,Zhu Weibing,Wang Heshun. Influence of orifice characteristics on static pressure dry gas seal performance[J]. Lubrication Engineering,2010,35(11):69-73.

[7] Zheng Shufei, Jiang Shuyun. Improved finite difference method for pressure distribution of aerostatic bearing[J]. Journal of Southeast University (English Edition),2009,25(2): 501-505.

[8] 閆飛,陳德林,丁雪興,等.靜壓干氣密封的開啟力與氣膜剛度的計算[J].甘肅科學學報,2014,26(4): 119-122.

Yan Fei,Chen Delin,Ding Xuexing et al.Calculation of opening force and gas film stiffness of aerostatic gas seal[J]. Journal of Gansu Sciences,2014,26(4):119-122.

[9] 杜建軍,劉暾,姚英學.狹縫節(jié)流氣體靜壓軸頸-止推軸承靜態(tài)特性分析[J].摩擦學學報,2002,22(1):66-70.

Du Jianjun,Liu Dun,Yao Yingxue.Analysis of the static performance of externally pressurized gas journal-thrust bearing with slot restrictors[J].Tribology,2002,22(1):66-70.

[10] 金朝旭,李雙喜,蔡紀寧,等. 可調(diào)控型氣膜潤滑密封靜壓結(jié)構(gòu)參數(shù)優(yōu)化[J].化工學報, 2015, 66(4):1425-1433.

Jing Zhaoxu, Li Shuangxi, Cai Jining, et al. Hydrostatic structural parameters optimization of regulatable gas lubricated seal[J].CIESC Journal, 2015, 66(4):1425-1433.

[11] 余建平. 靜壓氣體潤滑機械密封特性研究[D].昆明: 昆明理工大學, 2008．

[12] 劉墩, 劉育華, 陳世杰. 靜壓氣體潤滑[M]. 哈爾濱: 哈爾濱工業(yè)大學出版社, 1990．

[13] Fujiwara S，F(xiàn)use T. Advanced aerostatic dry gas seal[C]// 16th International Conference on Fluid Sealing,BHR Group 2000 Fluid Sealing.Brugge, Belgium：[s.n.],2000: 483-499.

[14] 楊沛然. 流體潤滑數(shù)值分析[M]. 北京: 國防工業(yè)出版社, 1998.

[15] 老大中. 變分法基礎(chǔ)[M]. 2版.北京:國防工業(yè)出版社, 2007.

[16] 楊環(huán).靜壓干氣密封流場數(shù)值模擬及其參數(shù)優(yōu)化[D].蘭州:蘭州理工大學,2012

Mechanical Properties of Aerostatic Dry Gas Seal Microscale Gas Film

Zhang Weizheng， Li Shuiping， Ding Xuexing， Xi Xilin

(CollegeofPetrochemicalEngineering,LanzhouUniversityofTechnology，LanzhouGansu730050,China)

According to the characteristics of aerostatic dry gas seals and the gas Reynolds equation, Galerkin method was used to derive the variational equations of pressure distribution of the gas film. Based on the gas film boundary conditions, the derivation process using the finite element method of the steady-state Reynolds equations was given. The pressure fitting curves of different thicknesses of the gas films were obtained by using polynomial fitting. In addition, the opening force was carried out. The fitting formula of the opening force and the film thinckness was obtained according to the opening force of the different film thicknesses.The relationship of the film stiffness and the film thickness was obtained by the derivation of the film thickness.The pressure distribution of the stationary seal ring end face was obtained. A parabolic decreasing trend is showed from the orifice to the inner and outer diameters in the radial direction of the stationary ring seal face. Pressure peak occurs at the orifice. With the increase of the gas film thickness, opening force and gas film stiffness decrease.

Numerical simulation; Gas flim stiffness; Steady state; Galerkin method; Microscale

2017-04-11

2017-04-21

國家自然科學基金項目(51565029,51165020);甘肅省自然科學基金項目(145RJZA083)。

張偉政(1978-)，男，博士，副教授，從事流體動密封與閥門相關(guān)技術(shù)研究；E-mail：zhangweiz@163.com。

丁雪興(1964-)，男，博士，教授，從事流體動密封方面研究；E-mail：xuexingding@163.com。

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

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

TE964； TQ 051

A

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

(編輯 王亞新)