代　迪，劉美紅

(昆明理工大學(xué) 機(jī)電工程學(xué)院，云南 昆明 650500)

在工業(yè)燃?xì)廨啓C(jī)和航空發(fā)動(dòng)機(jī)中，使用密封技術(shù)能明顯降低油耗、提高相關(guān)構(gòu)件的工作壽命和耐久性，并能提升整機(jī)的工作效率[1]。然而，航空發(fā)動(dòng)機(jī)的工作環(huán)境具有高環(huán)境溫度、高界面滑速和高邊界溫差的特點(diǎn)。并由于熱力變形及劇烈振動(dòng)，發(fā)動(dòng)機(jī)轉(zhuǎn)子存在較大的動(dòng)態(tài)位移[2-5]。傳統(tǒng)的迷宮密封技術(shù)、刷式密封和端面密封技術(shù)，難以適應(yīng)上述條件[6]。因此，研究可以實(shí)現(xiàn)較大柔性的氣膜密封技術(shù)變得越發(fā)重要。

目前，柔性氣膜密封技術(shù)主要包括端面氣膜和柱面氣膜[7-8]。端面氣膜由靜環(huán)和動(dòng)環(huán)兩個(gè)部分組成，常用的端面氣膜密封結(jié)構(gòu)包括“吸氣式”靜壓端面氣膜密封技術(shù)、反轉(zhuǎn)軸間端面氣膜密封技術(shù)和瑞利階梯型端面氣膜密封技術(shù)[9]。然而，端面氣膜密封技術(shù)存在較大的局限性，難以適應(yīng)航空發(fā)動(dòng)機(jī)的工作環(huán)境[10-11]。柔性支承柱面氣膜采用密封浮動(dòng)構(gòu)件和大柔性支承結(jié)構(gòu)制造[12]，在確保最大氣膜剛度和最小氣膜厚度的條件下，能抵消熱力變形和轉(zhuǎn)子的劇烈振動(dòng)。其相比于傳統(tǒng)的端面氣膜較易實(shí)現(xiàn)，且具有較大的徑向柔性[13]。典型的柱面氣膜密封結(jié)構(gòu)包括柔順箔密封、薄片板密封和直筒式柱面氣模密封[14]。

對(duì)柔性支承柱面氣膜的定性和定量分析，是設(shè)計(jì)更加合理的密封系統(tǒng)的基礎(chǔ)。然而建立合理的柱面氣膜密封系統(tǒng)的理論分析模型仍較為困難，現(xiàn)有的研究手段主要包括數(shù)值模擬和實(shí)驗(yàn)研究。本文首先在Proe[15]中建立柱面氣膜的密封參數(shù)化模型，然后基于有限元分析的Linear Soid劃分方法，使用ANSA網(wǎng)格劃分軟件[16]來(lái)劃分網(wǎng)格，并分析柔性支承柱面的靜壓、速度和剪應(yīng)力。同時(shí)，研究柱面結(jié)構(gòu)參數(shù)與操作參數(shù)對(duì)柔性支承柱面氣膜密封性能的影響。

1　柔性支承柱面氣膜密封原理與建模

柔性支承柱面氣膜密封系統(tǒng)通常由軸承座、波箔片、平箔和軸頸組成。柱面氣膜轉(zhuǎn)子的高速旋轉(zhuǎn)使得浮環(huán)和軸頸間產(chǎn)生動(dòng)壓密封氣膜，而靜子與浮環(huán)間具有較大的柔性支撐可抵消變形和振動(dòng)位移對(duì)密封系統(tǒng)的不良影響。同時(shí)，支承與氣膜間的共同作用可以避免密封表面間的磨損及接觸摩擦，并能實(shí)現(xiàn)更優(yōu)的密封效果。

為了分析柱面氣膜的各種性能，本文首先用Proe建立了柱面氣膜模型，如圖1所示。使用ANSA中的CFD版塊對(duì)所建立的模型進(jìn)行網(wǎng)絡(luò)劃分，并使用FLUENT的有限元計(jì)算程序設(shè)置邊界條件與求解器。最后，使用FLUENT的后處理軟件繪制模型的流體圖、壓力云圖和速度矢量場(chǎng)等。

圖1　柔性支承柱面氣膜模型

圖2　整體與局部網(wǎng)格劃分結(jié)果

流體網(wǎng)格劃分的質(zhì)量將直接影響CFD數(shù)值計(jì)算的速度和精度。本文使用Linear Solid劃分方法，首先通過(guò)控制外圓面的網(wǎng)格數(shù)劃分外圓面；然后將寬度方向設(shè)置為5層，并劃分為92 484個(gè)網(wǎng)格，以既能模擬出模型的實(shí)際流動(dòng)，又能節(jié)省運(yùn)算成本和時(shí)間。圖2所示為本文劃分結(jié)果的整體和局部展示。

本文為了便于計(jì)算與分析，忽略了氣膜體積力與慣性力，并假設(shè)動(dòng)環(huán)和靜環(huán)與氣膜間無(wú)相對(duì)滑動(dòng)，氣膜運(yùn)行穩(wěn)定、無(wú)振動(dòng)和擾動(dòng)。設(shè)置邊界條件為：進(jìn)口壓力pr=0.113 325 MPa，出口壓力ps=0.101 325 MPa；動(dòng)環(huán)壁面為旋轉(zhuǎn)Wall Motion，靜環(huán)壁面為Stationary Wall，動(dòng)環(huán)轉(zhuǎn)速為8 000 r/min；軟件求解器設(shè)置為SIMPLE壓力修正法。下面分別仿真分析柔性支承柱面的靜壓、速度和剪應(yīng)力分布。

(1)靜壓分析。如圖3所示為FLUENT軟件得到的柔性支承柱面氣膜靜壓分布圖。從圖中可以看出，氣膜最大壓力為0.116 MPa，且壓力分布具有流體憋油升壓效應(yīng)。其在氣膜從較大膜厚處流向最小膜厚處時(shí)，壓力快速增大，并形成一個(gè)負(fù)壓區(qū)。

圖3　柔性支承柱面氣膜靜壓分布圖

(2)速度分析。如圖4所示為FLUENT軟件得到的柔性支承柱面氣膜速度矢量圖。圖中兩種顏色分別為動(dòng)環(huán)內(nèi)壁與靜環(huán)內(nèi)壁處速度矢量?？梢钥闯?，動(dòng)環(huán)外壁處速度明顯大于靜環(huán)內(nèi)壁處速度。

圖4　柔性支承柱面氣膜速度矢量圖

(3)剪應(yīng)力分析。如圖5所示為FLUENT軟件得到的柔性支承柱面氣膜剪應(yīng)力分布圖。從圖中可以看出，在最小膜厚處剪應(yīng)力最大為1.54×10-4MPa，且能產(chǎn)生更強(qiáng)的剪切效應(yīng)。

圖5　柔性支承柱面氣膜剪應(yīng)力分布圖

2　密封性能分析

2.1　不同結(jié)構(gòu)參數(shù)下柱面氣膜密封性能分析

本文比較了在不同偏心率和間隙下柔性支承柱面氣膜的軸向靜壓、泄漏量、浮升力、摩擦轉(zhuǎn)矩和氣膜剛度等密封性能。

由流體動(dòng)壓原理可知，當(dāng)壓差和轉(zhuǎn)速等邊界條件恒定時(shí)，偏心率的大小將嚴(yán)重影響流體動(dòng)壓效應(yīng)。如圖6所示為偏心率 分別取0.2、0.4、0.6和0.8時(shí)，柱面氣膜的靜壓分布圖。從圖中可以看出，ε為0.2時(shí)，柱面氣膜下端面壓力出口處壓力明顯小于上端面壓力入口處壓力，形成軸向泊肅葉壓差流；ε為0.4、0.6和0.8時(shí)，柱面氣膜在較薄處壓力值較大，依此為0.113 MPa、0.116 MPa和0.128 MPa，并出現(xiàn)流體動(dòng)壓效應(yīng)，且 越大流體動(dòng)壓效應(yīng)越明顯。

圖6　不同偏心率時(shí)柱面氣膜的靜壓分布圖

當(dāng)ε取0.6時(shí)，密封間隙從0.01 mm變化到0.03 mm時(shí)，柱面氣膜的泄漏量、浮升力、摩擦轉(zhuǎn)矩和最大壓強(qiáng)等密封性能變化如圖7所示。從圖中可以看出，隨著密封間隙的增加，柔性支承柱面氣膜的泄漏量隨之增加，浮生力、摩擦轉(zhuǎn)矩和最大壓強(qiáng)隨之減小，柱面氣膜的密封性能逐漸降低。

圖7　不同密封間隙下柱面氣膜的泄漏量、浮升力、摩擦轉(zhuǎn)矩和最大壓強(qiáng)等密封性能變化

2.2　不同操作參數(shù)下柱面氣膜密封性能分析

本文使用控制變量法，分析不同操作參數(shù)下柱面氣膜的密封性能。如圖8所示為轉(zhuǎn)速分別為2 000 rpm、4 000 rpm、6 000 rpm和8 000 rpm時(shí)，柱面氣膜的靜壓分布圖。從圖中可知，隨著轉(zhuǎn)速的增加，流體動(dòng)壓效應(yīng)越明顯。

圖8　不同轉(zhuǎn)速下柱面氣膜靜壓分布圖

如圖9所示為氣體粘度從0.5×10-5Pa·s增加至3.5×10-5Pa·s時(shí)，柱面氣膜的泄漏量、浮升力、摩擦轉(zhuǎn)矩和最大壓強(qiáng)等密封性能變化情況。從圖中可以看出，隨著氣體粘度的增加，泄漏量和摩擦轉(zhuǎn)矩在逐漸減小。而浮升力先增加后不變，最大壓強(qiáng)在逐漸增加。由此表明，氣體粘度越大，柱面氣膜密封性能越好。

如圖10所示為邊界壓差從0.01 MPa增加到0.1 MPa時(shí)，柱面氣膜的泄漏量、浮升力、摩擦轉(zhuǎn)矩和最大壓強(qiáng)等密封性能變化情況。從圖中可以看出，隨著壓差增加，各項(xiàng)指標(biāo)均在增長(zhǎng)。表明壓差越大，柱面氣膜密封性能越好。

圖9　不同氣體粘度時(shí)柱面氣膜的泄漏量、浮升力、摩擦轉(zhuǎn)矩和最大壓強(qiáng)等密封性能變化

圖10　不同壓差時(shí)柱面氣膜的泄漏量、浮升力、摩擦轉(zhuǎn)矩和最大壓強(qiáng)等密封性能變化

3　結(jié)束語(yǔ)

本文使用CFD數(shù)值分析技術(shù)分析柔性支承柱面氣膜的密封性能。首先建立了柔性支承柱面氣膜的參數(shù)化模型，并仿真分析氣膜的靜壓、速度和剪應(yīng)力分布；然后，比較結(jié)構(gòu)參數(shù)與操作參數(shù)下軸向靜壓、泄漏量、浮升力、摩擦轉(zhuǎn)矩和氣膜剛度等密封性能的變化情況。仿真與測(cè)試結(jié)果表明，隨著轉(zhuǎn)速的增加，柱面氣膜表面流體動(dòng)壓效應(yīng)越明顯；隨著壓差和氣體粘度的增加，浮升力及泄漏量呈線性增加，柱面氣膜密封性能越好。

[1]Lattime S B,Braun M J,Choy F K.Design considerations towards the construction of Hybrid Floating Brush Seal (HFBS)[J].Tribology International, 2004,37(2):159-167.

[2]王虹,洪杰,馬艷紅.用于轉(zhuǎn)子系統(tǒng)的氣膜密封阻尼結(jié)構(gòu)穩(wěn)態(tài)特性[J].航空動(dòng)力學(xué)報(bào),2012,27(1):190-197.

[3]Pugachev A O,Deckner M.CFD prediction and test results of stiffness and damping coefficients for Brush-Labyrinth gas seals[C].Sydney:ASME Turbo Expo 2010: Power for Land, Sea and Air,2010.

[4]Wang H,Hong J,Ma Y H.Steady-state properties of gas film seal damper used in a rotor system[J]. Journal of Aerospace Power,2012,27(1):190-197.

[5]劉曉玉,王之櫟,丁蕾,等.反轉(zhuǎn)軸間雙端面氣膜密封性能分析[J].北京航空航天大學(xué)學(xué)報(bào),2017,43(3):608-614.

[6]孫軍鋒,劉美紅,李遇賢,等.基于CFD柱面氣膜密封靜態(tài)力學(xué)性能分析[J].排灌機(jī)械工程學(xué)報(bào),2017, 35(11):968-974.

[7]謝靜,白少先.高速氣流阻塞效應(yīng)對(duì)傾斜微孔端面密封動(dòng)壓特性的影響[J].摩擦學(xué)學(xué)報(bào),2017,37(6):806-813.

[8]陳志,高昀皞,趙鵬,等.高壓干氣密封流場(chǎng)數(shù)值模擬[J].四川大學(xué)學(xué)報(bào):工程科學(xué)版,2017,49(2):254-261.

[9]曹飛梅,姚平喜.不同結(jié)構(gòu)閥芯的滑閥流場(chǎng)CFD分析[J].液壓與氣動(dòng),2017(8):40-43.

[10] 張勇,谷正氣,劉水長(zhǎng),等.基于風(fēng)洞試驗(yàn)與CFD分析的汽車前后軸氣動(dòng)升力計(jì)算的研究[J].汽車工程,2017,39(3):275-280.

[11] 于騰祥.離心式空氣壓縮機(jī)級(jí)內(nèi)流動(dòng)特性及葉輪參數(shù)的數(shù)值分析[D].上海:華東理工大學(xué),2014.

[12] 張偉建.基于CFD的屋頂全開(kāi)型溫室自然通風(fēng)流場(chǎng)分析和降溫調(diào)控[D].鎮(zhèn)江:江蘇大學(xué),2016.

[13] 黃苗苗,吳乘勝,卜淑霞.重力和自由面影響下的氣泡尾流擴(kuò)散特性CFD分析[C].長(zhǎng)春:全國(guó)水動(dòng)力學(xué)研討會(huì),2017.

[14] 高宏閣.單缸風(fēng)冷汽油機(jī)缸頭熱負(fù)荷CFD分析[J].內(nèi)燃機(jī)與配件,2017(2):35-36.

[15] 張羽玲,陳允春,邵明月,等.基于MRF模型的冷卻輥冷卻性能的CFD分析[J].包裝工程,2017(7):11-14.

[16] 馬綱,孫曉軍,羅先海,等.端柱面組合密封氣膜穩(wěn)態(tài)特性數(shù)值仿真分析[J].北京航空航天大學(xué)學(xué)報(bào),2014,40(4):439-443.