李志剛,李軍,2,豐鎮(zhèn)平

(1.西安交通大學能源與動力工程學院,710049,西安;2.先進航空發(fā)動機協(xié)同創(chuàng)新中心,100191,北京)

?

進口防旋板對孔型密封非定常氣流激振特性的影響

李志剛1,李軍1,2,豐鎮(zhèn)平1

(1.西安交通大學能源與動力工程學院,710049,西安;2.先進航空發(fā)動機協(xié)同創(chuàng)新中心,100191,北京)

為了展示透平機械密封系統(tǒng)設計中進口防旋板在改善密封轉(zhuǎn)子動力特性和增強轉(zhuǎn)子系統(tǒng)穩(wěn)定性的止旋方面的性能,采用基于轉(zhuǎn)子多頻橢圓軌跡渦動模型和動網(wǎng)格技術的非定常數(shù)值方法,研究了進口防旋板對孔型密封非定常氣流激振力和轉(zhuǎn)子動力特性系數(shù)的影響,計算分析了進口預旋比分別為0和0.6時2種進口防旋板結(jié)構(gòu)的孔型密封轉(zhuǎn)子動力特性系數(shù)、非定常氣流激振力和流場特性,并與實驗結(jié)果進行了對比。研究表明:非定常數(shù)值方法能夠準確預測與頻率相關的孔型密封轉(zhuǎn)子動力特性系數(shù);進口預旋可顯著減小孔型密封的有效阻尼、增大有效阻尼項的穿越頻率,易誘發(fā)轉(zhuǎn)子渦動失穩(wěn);進口防旋板能夠有效減小孔型密封的進口旋流速度、孔型密封的交叉剛度和有效阻尼穿越頻率,增大孔型密封的有效剛度和有效阻尼,進而提高轉(zhuǎn)子系統(tǒng)的穩(wěn)定性。

孔型密封;進口預旋;防旋板;氣流激振;轉(zhuǎn)子動力特性

透平機械旋轉(zhuǎn)動密封在控制泄漏流動的同時,還會產(chǎn)生流體激振力,影響轉(zhuǎn)子的穩(wěn)定性[1]。目前透平機械主要通過采用先進阻尼密封結(jié)構(gòu)(蜂窩密封、孔型密封和袋型密封)、調(diào)整動靜間隙、安裝止旋裝置(分流注射器和防旋板)等,來抑制和解決密封氣流激振力引起的轉(zhuǎn)子失穩(wěn)問題[2-3]。安裝在密封進口的一系列與密封面成一定角度的導流板,能夠阻滯和引導密封進口的周向旋流。

迷宮密封進口周向旋流速度對密封交叉剛度具有決定性的作用,而氣流激振力是誘發(fā)轉(zhuǎn)子失穩(wěn)的主要因素。Benckert等通過實驗證明了在迷宮密封進口安裝防旋板結(jié)構(gòu)能夠減小密封交叉剛度[4-5],而Childs等研究表明防旋板對直接剛度和直接阻尼影響很小且實驗測量了迷宮密封在無防旋板、直角防旋板和反旋流角防旋板3種情況下的轉(zhuǎn)子動力特性系數(shù)[5-6]。Nielson等發(fā)現(xiàn)反旋流角防旋板結(jié)構(gòu)并沒有比直角結(jié)構(gòu)表現(xiàn)出更優(yōu)的性能[7]。圖1給出了迷宮密封進口直角防旋板和反旋流角防旋板結(jié)構(gòu)[7]。

工程應用中單獨采用防旋板結(jié)構(gòu)改善來流條件或采用阻尼密封代替迷宮密封的方法往往達不到完全消除轉(zhuǎn)子失穩(wěn)的目的。Scharrer報道了采用蜂窩密封代替迷宮密封無法完全消除航天飛機高壓氧渦輪泵轉(zhuǎn)子的亞同步渦動,但通過在蜂窩密封上游安裝防旋板達到了此目的[8],Zeidan等也是通過采用帶進口防旋板的蜂窩密封使多級離心壓氣機的轉(zhuǎn)子非同步渦動降低到了安全水平[9]。Childs等研究表明,減小密封進口正預旋或增大負預旋能夠顯著提高孔型密封的有效阻尼,增強轉(zhuǎn)子系統(tǒng)的穩(wěn)定性[10-11]。

(a)直角防旋板 (b)反旋流角防旋板圖1 迷宮密封進口防旋板結(jié)構(gòu)

針對旋轉(zhuǎn)阻尼密封轉(zhuǎn)子動力特性計算,晏鑫等提出了雙控制體Bulk Flow模型的方法和基于單頻渦動模型的非定常計算流體動力學(CFD)方法[12-13],李志剛等提出了基于多頻單向和橢圓渦動模型的非定常CFD方法[14-15]。

本文采用基于多頻橢圓軌跡渦動模型和動網(wǎng)格技術的非定常CFD方法,計算了2種進口預旋和防旋板結(jié)構(gòu)下的孔型密封轉(zhuǎn)子動力特性系數(shù)、非定常氣流激振力和流場,分析了進口預旋和防旋板對孔型密封非定常氣流激振特性的影響,以期為阻尼密封防旋板的設計提供理論基礎和技術數(shù)據(jù)。

1 計算模型和數(shù)值方法

1.1 計算模型

圖2給出了孔型密封的幾何結(jié)構(gòu)[10]。表1給出了孔型密封的幾何參數(shù)和實驗工況[10]。本文選取了2種進口預旋比u0(u0=0,0.6)并定義了密封進口的流動方向。u0為密封進口氣體周向旋流速度與轉(zhuǎn)子面線速度的比值,即

(1)

式中:Uθ為密封進口周向旋流速度;n為轉(zhuǎn)子轉(zhuǎn)速。式(1)可表征密封進口旋流速度的大小。根據(jù)圖1所示的2種迷宮密封進口防旋板結(jié)構(gòu),本文設計了2種孔型密封進口防旋板結(jié)構(gòu),如圖3所示,其由周向等弧度排列的36個導流片組成,且具有與密封相同的徑向間隙。

圖2 孔型密封實驗件幾何結(jié)構(gòu)

圖4給出了孔型密封和防旋板的三維計算模型和網(wǎng)格?？仔兔芊?、進口防旋板的計算模型和網(wǎng)格均為360°整周全三維結(jié)構(gòu)。采用ANSYS ICEM CFD軟件生成多塊結(jié)構(gòu)化網(wǎng)格,每個密封孔和防旋板導流片為O型網(wǎng)格,壁面區(qū)域采用加密網(wǎng)格,在密封間隙處沿徑向布置了20個節(jié)點。根據(jù)文獻[15]網(wǎng)格無關性分析,當網(wǎng)格節(jié)點數(shù)達到519萬時,本文計算可獲得網(wǎng)格無關解。因此,計算孔型密封的網(wǎng)格節(jié)點數(shù)為592萬,直角防旋板和反旋流角防旋板結(jié)構(gòu)的網(wǎng)格節(jié)點數(shù)分別為248萬和297萬。

表1 孔型密封實驗件幾何參數(shù)和實驗工況[10]

(a)直角防旋板

(b)反旋流角防旋板圖3 2種進口防旋板幾何結(jié)構(gòu)

圖4 孔型密封、進口防旋板計算模型和網(wǎng)格

1.2 數(shù)值方法

本文采用基于多頻橢圓軌跡渦動模型和動網(wǎng)格技術的非定常數(shù)值方法、采用ANSYS-CFX數(shù)值求解的Unsteady Reynolds-Averaged Navier-Stokes(URANS)方程,計算了孔型密封的非定常流場和流體激振力。圖5給出了轉(zhuǎn)子橢圓軌跡渦動模型。表2給出了具體的數(shù)值方法和渦動模型參數(shù)。多頻橢圓軌跡渦動模型中的轉(zhuǎn)子渦動位移的數(shù)學表達式如下。

x方向激勵

(2)

y方向激勵

(3)

式中:a和b分別為橢圓軌跡的長軸和短軸;N為渦動頻率數(shù);Ω0=2πf0為渦動基準頻率,f0為渦動頻率。

(a)x方向激勵 (b)y方向激勵圖5 轉(zhuǎn)子橢圓軌跡渦動模型示意圖

表2 數(shù)值方法和渦動模型參數(shù)

(4)

對于x方向激勵和y方向激勵,在對2種轉(zhuǎn)子渦動位移下的密封流場進行獨立的非定常求解后,通過快速傅里葉變換(FFT)可獲得式(4)在頻率內(nèi)的表達式,即

(5)

式中:Hij=Kij+iωCij為密封阻抗系數(shù)。通過非定常CFD數(shù)值求解和FFT變換可獲得頻率內(nèi)的密封流體激振力Fij和轉(zhuǎn)子渦動位移Dij,進而通過式(5)獲得密封轉(zhuǎn)子動力特性系數(shù)。本文非定常數(shù)值方法、多頻轉(zhuǎn)子橢圓軌跡渦動模型和轉(zhuǎn)子動力特性系數(shù)提取方法的詳細介紹和驗證可參考文獻[14-15]。

2 結(jié)果分析與討論

2.1 轉(zhuǎn)子動力特性系數(shù)

圖6給出了在2種進口預旋和防旋板結(jié)構(gòu)下,孔型密封轉(zhuǎn)子動力特性系數(shù)隨渦動頻率的變化。本文引入有效剛度系數(shù)Keff和有效阻尼系數(shù)Ceff來評價孔型密封的轉(zhuǎn)子動力特性,即

Keff=Kxx+CxyΩ

(6)

Ceff=Cxx+Kxy/Ω

(7)

為驗證數(shù)值方法的準確性,對進口無預旋(u0=0)、無進口防旋板時孔型密封的轉(zhuǎn)子動力特性系數(shù)進行了計算,并與實驗結(jié)果進行了比較。如圖6所示,除直接阻尼和有效阻尼在80～280Hz范圍內(nèi)略微偏小外(偏差小于20%),其余數(shù)值結(jié)果與實驗值吻合良好,表明本文的非定常數(shù)值方法能夠準確分析孔型密封的轉(zhuǎn)子動力特性系數(shù)。

(a)直接剛度

(b)直接阻尼

(c)交叉剛度

(d)交叉阻尼

(e)有效剛度

(f)有效阻尼圖6 孔型密封轉(zhuǎn)子動力特性系數(shù)隨渦動頻率的變化

如圖6所示:進口預旋會使孔型密封的直接剛度和交叉剛度增大,但使直接阻尼和交叉阻尼減小;相比于直接動力系數(shù),進口預旋會顯著影響交叉動力系數(shù);在80～280 Hz范圍內(nèi),進口預旋會使有效剛度減小、有效阻尼顯著減小、穿越頻率顯著增大。可見,進口預旋不但會減小轉(zhuǎn)子系統(tǒng)穩(wěn)定運行的頻率范圍(密封有效阻尼為正的頻率范圍)、引起轉(zhuǎn)子失穩(wěn),還會降低轉(zhuǎn)子系統(tǒng)的一階臨界轉(zhuǎn)速。因此,工程中應采取措施來減小密封進口預旋。

圖6還顯示:進口防旋板對孔型密封的直接剛度和直接阻尼影響較小,但使交叉剛度顯著減小、交叉阻尼顯著增大;進口防旋板使孔型密封的有效剛度和有效阻尼均顯著增大,有效阻尼項的穿越頻率顯著減小;2種進口防旋板結(jié)構(gòu)下的孔型密封具有幾乎相同的轉(zhuǎn)子動力特性系數(shù)。因此,采用進口防旋板能夠有效抑制轉(zhuǎn)子系統(tǒng)的振動,消除密封進口預旋引起的轉(zhuǎn)子失穩(wěn)問題,其中直角防旋板和反旋流角防旋板對孔型密封轉(zhuǎn)子動力特性具有同等的改善效果。

2.2 非定常氣流激振力和流場

圖7給出了進口預旋比u0=0.6時不同進口防旋板下孔型密封作用于轉(zhuǎn)子面上的非定常氣流激振力的變化。如圖7所示,進口防旋板能有效減小孔型密封的非定常氣流激振力。

(a)x方向氣流激振力

(b)y方向氣流激振力圖7 孔型密封內(nèi)非定常氣流激振力(u0=0.6,x方向激勵)

進口防旋板之所以能夠改變密封的非定常氣流激振特性,是因為其改變了密封進口流體的流動參數(shù),特別是進口旋流速度。表3給出了不同進口防旋板下密封處的流動參數(shù)。進口防旋板使密封處的旋流速度顯著減小,使泄漏量略微增大(小于4.5%),對靜壓和軸向速度影響很小。

表3 密封處的流動參數(shù)(u0=0.6)

圖8給出了密封通道內(nèi)平均旋流速度沿軸向的分布。圖9給出了3種防旋板下防旋板通道內(nèi)的三維流場。相比于無進口防旋板情況,2種防旋板結(jié)構(gòu)均能有效阻斷流體的周向流動,減小流體的周向速度。受轉(zhuǎn)子面的黏性剪切作用和孔型靜子面的阻滯耗散作用,防旋板的止旋效果僅在密封上游區(qū)域有效。2種防旋板結(jié)構(gòu)止旋效果相近,反旋流角防旋板結(jié)構(gòu)并未達到有效控制流體的流動方向、使流動反向旋轉(zhuǎn)的目的。這主要是相鄰導流片間的弧度間隔太大、流體未發(fā)生有效偏轉(zhuǎn)的緣故。因此,工程中需要對導流板的結(jié)構(gòu)參數(shù)(導流片數(shù)目、安裝角度、軸向長度等)進行優(yōu)化,使其能夠有效控制流體的流動方向。

圖8 密封通道內(nèi)平均旋流速度沿軸向的分布(u0=0.6)

(a)無防旋板 (b)直角防旋板

(c)反旋流角防旋板圖9 防旋板通道內(nèi)的三維流場(u0=0.6)

3 結(jié) 論

本文采用基于多頻橢圓軌跡渦動模型的非定常CFD方法,研究了進口預旋、進口導流板對孔型密封非定常氣流激振特性的影響。

進口預旋會使孔型密封交叉動力系數(shù)的絕對值顯著增大、密封有效剛度和有效阻尼顯著減小、有效阻尼項的穿越頻率明顯增大;進口預旋不但會減小轉(zhuǎn)子系統(tǒng)穩(wěn)定運行的頻率范圍、引起轉(zhuǎn)子失穩(wěn),同時還會降低轉(zhuǎn)子系統(tǒng)的一階臨界轉(zhuǎn)速。

進口防旋板能夠有效阻斷密封進口流體的周向流動、減小密封進口旋流速度,進而減小密封交叉剛度、增大有效剛度和有效阻尼、減小有效阻尼項的穿越頻率。直角防旋板和反旋流角防旋板對孔型密封轉(zhuǎn)子動力特性具有同等的改善效果,因此在密封進口處安裝防旋板是一種提高轉(zhuǎn)子系統(tǒng)阻尼、消除轉(zhuǎn)子渦動失穩(wěn)的有效措施。

[1] VANCE J M. Machinery vibration and rotordynamics [M]. New York, USA: John Wiley & Sons, 2010: 271-272.

[2] CHILDS D W, VANCE J M. Annular seals as tools to control rotor dynamic response of future gas turbine engines, AIAA 94-2804 [R]. Reston, VA, USA: AIAA, 1994.

[3] LI J, CHOUDHURY P D, KUSHNER F. Evaluation of centrifugal compressor stability margin and investigation of antiswirl mechanism [C]∥Proceeding of the 32nd Turbomachinery Symposium. Texas, USA: Texas A&M University, 2003: 49-57.

[4] BENCKERT H, WACHTER J. Flow induced spring coefficients of labyrinth seals for application in rotordynamics [C]∥Proceedings of a Workshop on Rotordynamic Instability Problems in High-Performance Turbomachinery. Texas, USA: Texas A&M University, 1980: 189-212.

[5] CHILDS D W, RAMSEY C. Seal rotordynamic coefficient test results for a model SSME ATD-HPFTP turbine interstage seal with and without a swirl brake [J]. ASME Journal of Tribology, 1991, 113(1): 198-203.

[6] CHILDS D W, MCLEAN J E, ZHANG M, et al. Rotordynamic performance of a negative-swirl brake for a tooth-on-stator labyrinth seal, GT2014-25577 [R]. New York, USA: ASME, 2014.

[7] NIELSON K K, CHILDS D W, MYLLERUP C M. Experimental and theoretical comparison of two swirl brake designs [J]. ASME Journal of Turbomachinery, 2001, 123(2): 353-358.

[8] SCHARRER J. Test results for leakage and rotordynamic coefficients;comparisons to labyrinth and smooth configurations [J]. ASME Journal of Triboloby, 1989, 111(2): 300-301.

[9] ZEIDAN F Y, PEREZ R X, STEPHENSON E M. The use of honeycomb seals in stabilizing two centrifugal compressors [C]∥Proceedings of 22nd Turbomachinery Symposium, Turbomachinery Laboratory. Texas, USA: Texas A&M University, 1993: 3-15.

[10]CHILDS D W, WADE J. Rotordynamic-coefficient and leakage characteristics for hole-pattern-stator annular gas seals-measurements versus predictions [J]. ASME Journal of Tribology, 2004, 126(2): 326-333.

[11]BROWN P, CHILDS D. Measurement versus predictions of rotordynamic coefficients of a hole-pattern gas seal with negative preswirl [J]. ASME Journal of Engineering for Gas Turbines and Power, 2012, 134(12): 122503.

[12]晏鑫, 李軍, 豐鎮(zhèn)平. BULK FLOW方法分析孔型密封轉(zhuǎn)子動力特性的有效性 [J]. 西安交通大學學報, 2009, 43(1): 24-28. YAN Xin, LI Jun, FENG Zhenping. Validation of two-control-volume BULK FLOW method for rotordynamic characteristics of hole-pattern seals [J]. Journal of Xi’an Jiaotong University, 2009, 43(1): 24-28.

[13]YAN X, LI J, FENG Z. Investigations on the rotordynamic characteristics of a hole-pattern seal using transient CFD and periodic circular orbit model [J]. ASME Journal of Vibration and Acoustics, 2011, 133(4): 041007.

[14]李志剛, 李軍, 豐鎮(zhèn)平. 袋型阻尼密封轉(zhuǎn)子動力特性的多頻單向渦動預測模型 [J]. 西安交通大學學報, 2012, 46(5): 13-18. LI Zhigang, LI Jun, FENG Zhenping. Multiple frequencies one-dimensional prediction model of rotordynamic characteristics for pocket damper seals [J]. Journal of Xi’an Jiaotong University, 2012, 46(5): 13-18.

[15]LI Zhigang, LI Jun, YAN Xin. Multiple frequencies elliptical whirling orbit model and transient RANS solution approach to rotordynamic coefficients of annual gas seals prediction [J]. ASME Journal of Vibration and Acoustics, 2013, 135(3): 031005.

(編輯 苗凌)

Effect of Swirl Brakes on Unsteady Flow Excitation Characteristics of Hole-Pattern Seal

LI Zhigang1,LI Jun1,2,FENG Zhenping1

(1. School of Energy & Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China; 2. Collaborative Innovation Center of Advanced Aero-Engine, Beijing 100191, China)

To reveal the anti-swirl performance of the swirl brake on improving the seal rotordynamic performance and enhancing rotor stability in turbomachinery, the effects of swirl brakes on the transient flow excitation characteristics and rotordynamic coefficients of a hole-pattern seal are numerically investigated. With computational fluid dynamics (CFD) method based on the multi-frequency elliptical orbit rotor whirling model and mesh deformation technique, the rotordynamic coefficients, transient fluid response forces and flow pattern of the hole-pattern seal are solved and analyzed for two types of swirl brakes at inlet preswirl ratios of 0 and 0.6. The predicted rotordynamic coefficients are compared with the experiment data. The numerical results show that the present transient numerical method enables to accurately predict the frequency-dependent rotordynamic coefficients of hole-pattern seal. Inlet preswirl significantly decreases effective damping in the hole-pattern seal and increases in the crossover frequency of the effective damping term, which usually induces rotor instability vibration. The swirl brake is very useful to eliminate rotor instability vibration because it can obviously reduce the cross-coupling stiffness and the crossover frequency, and increase the hole-pattern seal effective stiffness and effective damping by reducing the inlet preswirl velocity.

hole-pattern seal; inlet preswirl; swirl brake; flow excitation; rotor dynamic

2015-03-30。 作者簡介:李志剛(1986—),男,講師。 基金項目:國家自然科學基金資助項目(51406144);博士后科學基金資助項目(2014M550490)。

時間:2015-10-23

網(wǎng)絡出版地址:http:∥www.cnki.net/kcms/detail/61.1069.T.20151023.1106.012.html

10.7652/xjtuxb201601003

TK474

A

0253-987X(2016)01-0016-06