張樹強， 王 良， 陳 杰， 趙偉剛

(西安航天動力研究所，西安 710100)

氣體端面密封以其不可替代的優(yōu)勢在高速旋轉的設備中被廣泛采用，穩(wěn)健的密封氣膜是其正常工作的重要保證。理想情況下密封動、靜環(huán)端面均垂直于轉軸，密封氣膜穩(wěn)定，然而在實際工作中，由于制造誤差、裝配誤差等、熱變形、力變形以及轉軸振動等原因，密封氣膜是隨時間動態(tài)變化的。如若補償環(huán)(一般為靜環(huán))不能對非補償環(huán)(一般為動環(huán))的振動提供良好的動態(tài)跟蹤響應(稱為追隨性)，密封端面間的氣膜將會發(fā)生失穩(wěn)，出現(xiàn)密封泄漏量增大或者端面接觸摩擦的現(xiàn)象，最終導致密封失效。

諸多學者對氣體端面密封的動態(tài)特性進行了深入研究。國內(nèi)，張樹強等[1]基于攝動法研究了動靜壓混合式氣體端面密封追隨性；宋鵬云等[2]研究了實際氣體效應對T形槽干氣密封動態(tài)剛度和阻尼系數(shù)的影響規(guī)律；劉向鋒等[3]等基于半解析法研究了極端工況下干氣密封動態(tài)特性；丁雪興等[4]研究了熱耗散變形下干氣密封系統(tǒng)的軸向振動穩(wěn)定性。國外，Green等[5]研究了不同轉軸速度和密封環(huán)錐角下氣體端面密封的瞬態(tài)響應；Miller等[6-7]先后采用直接數(shù)值求解法和半解析法研究了螺旋槽氣體端面密封的動態(tài)特性；Ruan[8]基于攝動法研究了螺旋槽氣體端面密封的追隨性，并且采用直接數(shù)值計算法驗證了求解結果；Blasiak等[9]采用直接數(shù)值求解法對比研究了不同槽型氣體端面密封動態(tài)特性。眾多研究主要集中在密封操作和結構參數(shù)對其動態(tài)特性的影響，雖然考慮了彈簧和補償環(huán)用輔助密封圈(通常為O形圈)剛度和阻尼的影響，但均在給定值下進行的，并未進行系列化的優(yōu)化研究。鑒于此，本文研究了彈簧和補償環(huán)用O形圈剛度和阻尼對目前常見的三種氣體端面密封(動壓式密封、靜壓式密封以及動靜壓混合式密封)動態(tài)追隨性的影響。

目前，氣體端面密封動態(tài)特性有兩種研究方法：一種是直接數(shù)值求解法(Direct Numerical Simulation Method,DNSM)[10-11]；另一種是線性化的數(shù)值求解方法，包括攝動法(Small Perturbation Method,SPM)[12]、步進法(Step Jump Method,SJM)[13]以及直接數(shù)頻響應法(Direct Numerical Frequency Response Method,DNFRM)[14]。直接數(shù)值求解法雖然能夠完全獲得密封動態(tài)特性，但是數(shù)值求解耗時，且不利于動態(tài)特性影響因素參數(shù)化的設計研究，實踐證明線性化的數(shù)值求解方法對于密封的動態(tài)特性研究也是準確和實用，因而本文基于攝動法進行彈簧和O形圈剛度和阻尼對氣體端面密封追隨性的影響研究。

1 密封結構及工作原理

氣體端面密封結構簡圖如圖1所示，動環(huán)固連于轉軸，靜環(huán)通過O形圈和彈簧安裝在密封殼體中，O形圈提供輔助密封，彈簧提供閉合力。密封工作時可以在動環(huán)和靜環(huán)端面間形成穩(wěn)定的氣膜，密封氣膜具有剛度和阻尼特性，在阻止密封腔介質(壓力為po)泄漏的同時，保證密封端面非接觸運轉，避免摩擦磨損，且具有較低的泄漏率。

1-動環(huán);2-靜環(huán);3-螺旋槽;4-均壓槽;5-節(jié)流孔;6-O形圈;7-彈簧;8-彈簧座;9-轉軸

圖2所示為三種不同結構的氣體端面密封，圖2(a)為動壓式氣體端面密封，在動環(huán)端面上加工有螺旋槽，依靠流體動壓效應實現(xiàn)非接觸運轉；圖2(b)為靜壓式氣體端面密封，在靜環(huán)端面上加工有節(jié)流孔和均壓槽，二者相互連通，節(jié)流孔的作用是將外部阻封氣(氣源壓力為ps，經(jīng)過節(jié)流后壓力變?yōu)閜d)導入密封端面，提供流體靜壓承載能力，節(jié)流孔實際上產(chǎn)生流動阻力增加氣膜的剛度和穩(wěn)定性，均壓槽的作用是將通過節(jié)流孔輸入的阻封氣在周向進行均布；圖2(c)為動靜壓混合式氣體端面密封，靜環(huán)端面結構與靜壓式密封完全一致，在動環(huán)端面上加工有“魚骨”形螺旋槽，動靜壓混合式密封同時具備動壓式密封和靜壓式密封的特點。本文研究的參數(shù)如表1所示，包括密封結構參數(shù)、操作參數(shù)以及動環(huán)微擾幅值等。

(a) 動壓式密封

(b) 靜壓式密封

(c) 動靜壓混合式密封

圖2 密封端面結構

Fig.2 Ring geometry of gas face seal

2 追隨性研究數(shù)學模型

為研究氣體端面密封的追隨性，將密封裝置簡化為彈簧-阻尼-質量系統(tǒng)(如圖3所示)：密封氣膜具有剛度(k)和阻尼(c)特性；O形圈具有阻尼(cs)特性，彈簧阻尼和O形圈阻尼相比可以忽略不計；彈簧和O形圈共同具有剛度(ks)特性，本文稱為綜合剛度，即兩者剛度之和；密封環(huán)具有質量(m)特性。基于攝動法求解控制氣膜壓力分布的動態(tài)非線性Reynolds方程，得到氣膜微擾壓力，進而求解表征密封氣膜動態(tài)特性的剛度和阻尼系數(shù)。在動環(huán)微擾下建立靜環(huán)動力學方程，求解方程研究密封的動態(tài)追隨性。由于本文目的是研究彈簧和O形圈剛度和阻尼對氣體端面密封追隨性的影響規(guī)律，因而基于攝動法進行氣體端面密封追隨性的詳細數(shù)值求解過程這里不再敘述。

2.1 氣膜潤滑方程

控制密封端面間流體流動，且能計及小孔節(jié)流動靜壓混合作用，等溫常黏度的可壓縮二維柱坐標非穩(wěn)態(tài)Reynolds方程為[15-16]

表1 密封分析參數(shù)

圖3 氣體端面密封簡化模型

(1)

式(1)具有如下的邊界條件

p(ri,θ,t)=pi

p(ro,θ,t)=po

p(r,0,t)=p(r,2π,t)

(2)

如圖4所示，密封在擾動頻率υ下有軸向和兩個角向的微擾運動z(t)、α(t)、β(t)，式(1)中的瞬態(tài)氣膜厚度為

h=ho+z(t)+rsinθα(t)-rcosθβ(t)+〈hg〉+〈hd〉

(3)

式中：ho為穩(wěn)態(tài)氣膜厚度；hg為動壓槽深度；hd為均壓槽深度；〈〉為只在槽區(qū)計入氣膜厚度公式。

引入的微擾運動為簡諧運動

z(t)=|Δz|eiυt,α(t)=|Δα|eiυt,β(t)=|Δβ|eiυt

(4)

微擾運動可在穩(wěn)態(tài)氣膜壓力p0的基礎上引起微擾壓力Δp，利用一階泰勒級數(shù)展開，瞬態(tài)氣膜壓力(含節(jié)流孔后壓力)表示為

(5)

圖4 密封動力學模型

2.2 氣膜動態(tài)剛度和阻尼定義

微擾時氣膜對靜環(huán)的軸向力Fz和角向偏轉力矩Mx、My的增量由微擾壓力Δp引起，表達為

(6)

將Δp代入式(6)可得

(7)

由于氣膜被視為具有剛度和阻尼特性，所以氣膜對靜環(huán)作用力的增量又有

(8)

式中：kij，cij(i，j=z,x,y)為氣膜剛度和阻尼系數(shù)，負號代表氣膜阻止微擾運動。

對比式(7)和式(8)，可得到由擾動壓力所表達的氣膜剛度和阻尼系數(shù)

(9)

氣體端面密封氣膜18個動態(tài)性能參數(shù)剛度和阻尼系數(shù)具有如下特征：

(1) 由兩個角向(x和y)微擾引起的角向剛度和阻尼滿足kxx=kyy、kxy=-kyx、cxx=cyy和cxy=-cyx的關系。

(2) 軸向和角向耦合剛度和阻尼系數(shù)(kzx、kzy、kxz、kyz、czx、czy、cxz、cyz)均為零。

(3) 氣膜動態(tài)剛度和阻尼系數(shù)受擾動頻率υ影響。

因此，氣膜動態(tài)剛度和阻尼系數(shù)可以簡化為

(10)

2.3 補償環(huán)運動方程

在動環(huán)微擾運動的作用下，求解靜環(huán)的響應運動，就可以得到氣體端面密封的動態(tài)氣膜厚度分布規(guī)律，研究密封追隨性。由氣膜動態(tài)特性系數(shù)滿足的規(guī)律可知，靜環(huán)的軸向運動和角向運動是解耦的，其在動環(huán)微擾作用下沿3個方向的運動方程為

(11)

(12)

式中：z、α和β分別為靜環(huán)軸向響應以及擾x軸和y軸的角向響應；zr、αr和βr分別為動環(huán)軸向激勵以及擾x和y軸的角向激勵；m為靜環(huán)質量；Ix，Iy為靜環(huán)轉動慣量；ks為彈簧和O形圈綜合剛度;cs為O形圈阻尼；ksx和ksy為彈簧和O形圈的角向綜合剛度;csx和csy為O形圈的角向阻尼，計算公式為

(13)

為求解式(11)和式(12)，需要知道氣膜動態(tài)剛度和阻尼系數(shù)，氣膜剛度和阻尼系數(shù)又受到靜環(huán)微擾運動的頻率υ影響，所以需要聯(lián)立潤滑方程進行迭代求解。慶幸的是式(11)和式(12)屬于強迫運動的范疇，靜環(huán)微擾運動的頻率由動環(huán)激勵頻率決定，因此可以在給定動環(huán)激勵頻率的情況下，先求解潤滑方程得到氣膜動態(tài)剛度和阻尼系數(shù)，進而求解靜環(huán)微擾運動，研究追隨性。由于動環(huán)剛性地固連于轉軸，所以激勵頻率和轉軸角速度相等，故在此研究激勵頻率等于轉軸角速度ω的情形。

3 研究結果分析

在擾動頻率υ等于轉軸角速度ω的情況下，數(shù)值求解了三種密封的氣膜動態(tài)剛度和阻尼系數(shù)，如表2所示。

表2 密封剛度和阻尼系數(shù)

假定動環(huán)的激勵形式為簡諧振動zr=arzsin(ωt)，αr=arxcos(ωt)，βr=arysin(ωt)，其中arz、arx和ary分別為激勵的幅值，基于表1中的參數(shù)研究了彈簧和O形圈的剛度和阻尼對三種氣體端面密封追隨性的影響規(guī)律。忽略靜環(huán)響應運動與動環(huán)激勵運動的相位差，以靜環(huán)響應振動的幅值(分別為asz、asx和asy)大小來表征密封的追隨性，靜環(huán)響應振動的幅值越接近動環(huán)激勵的幅值，密封追隨性越好。

引入如下無量綱項

(14)

3.1 彈簧和O形圈綜合剛度對密封追隨性的影響

求解了彈簧和O形圈綜合剛度對靜環(huán)無量綱響應振動幅值的影響規(guī)律(由于密封端面的軸對稱性，靜環(huán)繞x軸和y軸的振動幅值相等)，如圖5所示。由圖可知當綜合剛度小于107N/m時，三種密封靜環(huán)的軸向和角向無量綱響應振動幅值基本保持不變；當綜合剛度超過107N/m時，三種密封靜環(huán)的各向無量綱響應振動幅值逐漸減小。這是由于當綜合剛度的數(shù)量級接近或者超過密封氣膜軸向剛度kzz的數(shù)量級時(動壓式密封和動靜壓混合式密封為108N/m，靜壓式密封為107N/m)，密封振動主要發(fā)生在密封靜環(huán)端面，而不是靜環(huán)背面與彈簧接觸處。因此，為了使密封具有良好的追隨性，綜合剛度不應該超過107N/m。通常情況下氣體端面密封用彈簧和O形圈的綜合剛度量級不會超過106N/m，因而密封不會因為剛度的選取不當而引起追隨性不好。

為了進一步研究綜合剛度對氣體端面密封追隨性的影響規(guī)律，分別求解了綜合剛度數(shù)量級低于和等于107N/m時靜環(huán)的響應運動，如圖6所示。圖6(a)～圖6(e)為綜合剛度等于5×106N/m時，靜環(huán)的軸向和角向響應運動，圖6(b)～圖6(f)為綜合剛度等于5×107N/m時，靜環(huán)的軸向和角向響應運動。由計算結果對比分析可知，當綜合剛度為5×106N/m時(數(shù)量級低于107N/m)，三種密封的靜環(huán)針對動環(huán)簡諧微擾在各個方向均可以提供良好的跟蹤響應，響應運動與微擾運動同步，沒有相位差，且簡諧運動的幅值也基本相等。

(a) 軸向振幅

(b) 角向振幅

Fig.5 Influence of stiffness on dimensionless amplitudes of stator

當綜合剛度為5×107N/m時(數(shù)量級等于107N/m)，三種密封的靜環(huán)針對動環(huán)簡諧微擾在各個方向均未能提供良好的跟蹤響應，響應運動與微擾運動出現(xiàn)了相位差，且靜環(huán)響應運動的幅值低于動環(huán)微擾運動的幅值。由計算結果亦可知，在同樣參數(shù)下三種密封追隨性的優(yōu)劣為：動壓式密封最好，動靜壓混合式密封次之，靜壓式密封最差。

3.2 密封圈阻尼對密封追隨性的影響

求解了密封圈阻尼對靜環(huán)無量綱響應振動幅值的影響規(guī)律，如圖7所示。由圖可知，三種密封靜環(huán)各向無量綱響應振動幅值的變化規(guī)律一樣，對于動壓式和動靜壓混合式密封而言，當密封圈阻尼小于5×104N·s·m-1時，靜環(huán)的各向無量綱響應振動幅值基本不變，隨著密封圈阻尼的繼續(xù)增大，無量綱響應振動幅值隨之減小，密封的追隨性變差。對靜壓式密封而言，當密封圈阻尼小于2×103N·s·m-1時，靜環(huán)的各向無量綱響應振動幅值基本不變，隨著密封圈阻尼的繼續(xù)增大，無量綱響應振動幅值隨之減小，密封的追隨性亦變差。因此，為了使氣體端面密封具有良好的追隨性，動壓式和動靜壓混合式密封補償環(huán)用輔助密封圈的阻尼不應超過5×104N·s·m-1，而靜壓式密封不應超過2×103N·s·m-1。通常情況下氣體端面密封用O形圈阻尼量級不會超過103N·s·m-1，因而密封不會因為O形圈阻尼選取不當而引起追隨性不好。

(a) ks=5×106 N/m

(b) ks=5×107 N/m

(c) ks=5×106 N/m

(d) ks=5×107 N/m

(e) ks=5×106 N/m

(f) ks=5×107 N/m

(a) 軸向振幅

(b) 角向振幅

Fig.7 Influence of secondary seal damping on dimensionless amplitudes of stator

同樣，為了進一步研究密封圈阻尼對氣體端面密封追隨性的影響規(guī)律，求解了不同數(shù)量級阻尼時靜環(huán)的響應運動，如圖8所示。圖8(a)～圖8(e)為密封圈阻尼等于2×103N·s·m-1時，靜環(huán)的軸向和角向響應運動，圖8(b)～圖8(f)為密封圈阻尼等于1×105N·s·m-1時，靜環(huán)的軸向和角向響應運動。由計算結果對比分析可知，當密封圈阻尼為2×103N·s·m-1時，三種密封的靜環(huán)針對動環(huán)簡諧微擾在各個方向均可以提供良好的跟蹤響應，響應運動與微擾運動同步，沒有相位差，且簡諧運動的幅值也基本相等。當密封圈阻尼等于1×105N·s·m-1時，三種密封的靜環(huán)針對動環(huán)簡諧微擾在各個方向均未能提供良好的跟蹤響應，響應運動與微擾運動出現(xiàn)了相位差，且靜環(huán)響應運動的幅值低于動環(huán)微擾運動的幅值。由計算結果亦可知，在同樣參數(shù)下三種密封追隨性的優(yōu)劣為：動壓式密封最好，動靜壓混合式密封次之，靜壓式密封最差。

4 結 論

(1) 當彈簧和O形圈綜合剛度低于某數(shù)量級時，三種密封在軸向和角向均能具有優(yōu)異的追隨性，隨著綜合剛度的增加，密封追隨性隨之變差。為了使氣體端面密封具有良好的追隨性，綜合剛度不應該超過107N/m。

(a) cs=2×103 N·s·m-1

(b) cs=1×105 N·s·m-1

(c) cs=2×103 N·s·m-1

(d) cs=1×105 N·s·m-1

(e) cs=2×103 N·s·m-1

(f) cs=1×105 N·s·m-1

(2) 當輔助密封圈阻尼低于某數(shù)量級時，三種密封在軸向和角向均能具有優(yōu)異的追隨性，隨著輔助密封圈阻尼的增加，密封追隨性隨之變差。為了使氣體端面密封具有良好的追隨性，動壓式和動靜壓混合式密封補償環(huán)用輔助密封圈的阻尼不應超過5×104N·s·m-1，而靜壓式密封不應超過2×103N·s·m-1。

(3) 在同樣的操作條件下，螺旋槽密封追隨性最好，動靜壓混合密封次之，靜壓密封最差。

[1] 張樹強, 李雙喜, 蔡紀寧, 等. 動靜壓混合式氣體密封追隨性及主動調控振動特性數(shù)值分析[J]. 航空學報, 2012, 33(7): 1336-1346.

ZHANG Shuqiang, LI Shuangxi, CAI Jining, et al. Numerical analysis for the tracking property and active regulation vibration characteristics of dynamic-hydrostatic hybrid gas seals[J]. Acta Aeronautica et Astronautica Sinica, 2012, 33(7): 1336-1346.

[2] 宋鵬云, 胡曉鵬, 許恒杰. 實際氣體對T槽干氣密封動態(tài)特性的影響[J]. 化工學報, 2014, 65(4): 1344-1352.

SONG Pengyun, HU Xiaopeng, XU Hengjie. Effect of real gas on dynamic performance of T-groove dry gas seal[J]. CIESC Journal, 2014, 65(4): 1344-1352.

[3] 劉向鋒, 徐辰, 黃偉峰. 基于半解析法的極端工況干氣密封動態(tài)特性研究與參數(shù)設計[J]. 清華大學學報(自然科學版), 2014, 54(2): 223-228.

LIU Xiangfeng, XU Chen, HUANG Weifeng. Analysis and parametric design of the dynamics of a dry gas seal for extreme operating conditions using a semi-analytical method[J]. Journal Tsinghua University (Natural Science), 2014, 54(2): 223-228.

[4] 丁雪興, 陸俊杰, 劉勇, 等. 熱耗散變形下干氣密封系統(tǒng)軸向振動穩(wěn)定性分析[J]. 振動工程學報, 2016, 29(1): 79-86.

DING Xuexing, LU Junjie, LIU Yong, et al. Stability analysis on axial vibration of dry gas seal system under the thermo-elastic deformation considering the thermal dissipation[J]. Journal of Vibration Engineering, 2016, 29(1): 79-86.

[5] GREEN I, BARNSBY R M. A simultaneous numerical solution for the lubrication and dynamic stability of noncontacting gas face seals[J]. Journal of Tribology, 2001, 123(2): 388-394.

[6] MILLER B A, GREEN I. Numerical formulation for the dynamic analysis of spiral-grooved gas face seals[J]. Journal of Tribology, 2001, 123(2): 395-403.

[7] MILLER B A, GREEN I. Semi-analytical dynamic analysis of spiral-grooved mechanical gas face seals[J]. Journal of Tribology, 2003, 125(2): 403-413.

[8] RUAN B. A semi-analytical solution to the dynamic tracking of non-contacting gas face seals[J]. Journal of Tribology, 2002, 124(1): 196-202.

[9] BLASIAK S, ZAHORULKO A V. A parametric and dynamic analysis of non-contacting gas face seals with modified surfaces[J]. Tribology International, 2016, 94: 126-137.

[10] RUAN B. Numerical modeling of dynamic sealing behaviors of spiral groove gas face seals[J]. Journal of Tribology, 2002, 124(1): 186-195.

[11] GREEN I, BARNSBY R M. A parametric analysis of the transient forced response of noncontacting gas coned face seals[J]. Journal of Tribology, 2002, 124(2): 151-157.

[12] ZIRKELBACK N, ANDRES L S. Effect of frequency excitation on force coefficients of spiral groove gas seals[J]. Journal of Tribology, 1999, 121(4): 853-863.

[13] ELROD H G, JR. MCCABE J T, CHU T Y. Determination of gas-bearing stability by response to a step-jump[J]. Journal of Lubrication Technology, 1967, 89(4): 493-498.

[14] MILLER B A, GREEN I. Numerical techniques for computing rotor-dynamic properties of mechanical gas face seals[J]. Journal of Tribology, 2002, 124(4): 755-761.

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

[16] 劉雨川, 徐萬孚, 陳國林, 等. 端面氣膜密封的高性能端面結構[J]. 航空學報, 2000, 21(2): 187-190.

LIU Yuchuan, XU Wanfu, CHEN Guolin, et al. Face configuration with high performance for gas film face seal[J]. Acta Aeronautica et Astronautica Sinica, 2000, 21(2): 187-190.