馬浩原，張萬(wàn)福, 2

(1.上海理工大學(xué) 能源與動(dòng)力工程學(xué)院,上海 200093；2.上海理工大學(xué) 能源與動(dòng)力工程學(xué)院，上海市動(dòng)力工程多相流動(dòng)與傳熱重點(diǎn)實(shí)驗(yàn)室，上海 200093)

迷宮式汽封流體動(dòng)力特性研究

馬浩原1，張萬(wàn)福1, 2

(1.上海理工大學(xué) 能源與動(dòng)力工程學(xué)院,上海 200093；2.上海理工大學(xué) 能源與動(dòng)力工程學(xué)院，上海市動(dòng)力工程多相流動(dòng)與傳熱重點(diǎn)實(shí)驗(yàn)室，上海 200093)

建立了迷宮式汽封的全三維模型，采用CFD商業(yè)軟件Fluent對(duì)不同工況下汽封的流場(chǎng)進(jìn)行了計(jì)算分析，研究了轉(zhuǎn)子所受的氣流力與轉(zhuǎn)速、進(jìn)出口壓力比及偏心率的關(guān)系，得到了影響轉(zhuǎn)子穩(wěn)定性的密封動(dòng)力特性系數(shù)和渦動(dòng)系數(shù).結(jié)果表明：隨著轉(zhuǎn)速、進(jìn)出口壓力比和偏心率的增大，渦動(dòng)系數(shù)也相應(yīng)增大，轉(zhuǎn)子的穩(wěn)定性下降.

迷宮式汽封；氣流力；動(dòng)力特性；數(shù)值計(jì)算

汽封是汽輪機(jī)等透平機(jī)械中減少流體泄漏的關(guān)鍵部件.常見(jiàn)的汽封結(jié)構(gòu)有迷宮式[1]、刷式[2]、階梯式[3]、蜂窩式[4]、指尖式[5]、螺旋槽式[6]和混合式[7]等，而迷宮式汽封在透平機(jī)械中最為常見(jiàn)，已得到廣泛應(yīng)用[8].隨著機(jī)組向大容量和高參數(shù)方向發(fā)展，汽封氣流力對(duì)轉(zhuǎn)子穩(wěn)定性的影響更加顯著，汽封中流體周向壓力分布不均是引起氣流力的主要原因.由于存在偏心，轉(zhuǎn)子在自轉(zhuǎn)的同時(shí)還會(huì)發(fā)生渦動(dòng).在轉(zhuǎn)子發(fā)生渦動(dòng)時(shí)，轉(zhuǎn)子與汽封之間的間隙會(huì)發(fā)生周期性變化，進(jìn)而使其中的流體壓力產(chǎn)生周期性脈動(dòng)，而壓力脈動(dòng)滯后轉(zhuǎn)子位移一個(gè)相位角，在與轉(zhuǎn)子偏心方向垂直的方向產(chǎn)生一個(gè)加劇轉(zhuǎn)子渦動(dòng)的激振力，流體形成自激振動(dòng)[9].然而這種自激振動(dòng)無(wú)法通過(guò)動(dòng)平衡的方法消除[10]，輕則影響機(jī)組可靠性，重則使機(jī)組被迫停機(jī)，無(wú)法正常運(yùn)行.汽封流體動(dòng)力特性是設(shè)計(jì)時(shí)必須考慮的因素之一，以往的研究[11]主要集中在汽封靜特性方面.

筆者主要研究了汽封對(duì)轉(zhuǎn)子動(dòng)力特性的影響，采用計(jì)算流體力學(xué)(CFD)方法建立了迷宮式汽封的三維偏心渦動(dòng)數(shù)值分析模型，對(duì)汽封的偏心渦動(dòng)進(jìn)行數(shù)值模擬，研究了轉(zhuǎn)速、進(jìn)出口壓力比及偏心率對(duì)密封動(dòng)力特性系數(shù)和渦動(dòng)系數(shù)的影響，并分析了其對(duì)轉(zhuǎn)子動(dòng)態(tài)穩(wěn)定性的影響.

1 汽封流體動(dòng)力特性模型

1.1 汽封流體動(dòng)力特性模型

研究汽封中流體對(duì)轉(zhuǎn)子作用力的動(dòng)力學(xué)分析時(shí)一般采用線(xiàn)性化模型[12].如圖1所示，當(dāng)轉(zhuǎn)子在靜平衡位置受到位移或速度擾動(dòng)時(shí)，氣流作用在轉(zhuǎn)子上的氣流力會(huì)發(fā)生變化，當(dāng)擾動(dòng)量是微小量時(shí)，可將氣流力分解為水平和垂直2個(gè)分量，即Fx和Fy.

圖1 汽封流體動(dòng)力特性模型

定義剛度系數(shù)k為單位位移擾動(dòng)量引起的氣流力增量，阻尼系數(shù)c為單位速度擾動(dòng)量所引起的氣流力增量，剛度系數(shù)和阻尼系數(shù)統(tǒng)稱(chēng)為密封動(dòng)力特性系數(shù).則動(dòng)態(tài)氣流力的矩陣形式為

(1)

假設(shè)轉(zhuǎn)子為同心渦動(dòng)，即渦動(dòng)中心為靜子中心.在同心渦動(dòng)中，密封腔中靜態(tài)流動(dòng)參數(shù)沿周向分布均勻，密封動(dòng)力特性系數(shù)中直接剛度系數(shù)、阻尼系數(shù)項(xiàng)kii、cii(i=x,y)相等，交叉耦合剛度系數(shù)、阻尼系數(shù)項(xiàng)kij、cij(i=x,y;j=x,y)幅值相等，符號(hào)相反，即

則氣流力可表示為

(2)

式中：Fr為轉(zhuǎn)子所受徑向力，N；Ft為轉(zhuǎn)子所受切向力，N；δ為偏心率，mm；Ω為渦動(dòng)速度，r/min；K、k分別為直接剛度系數(shù)和交叉耦合剛度系數(shù)；C、c分別為直接阻尼系數(shù)和交叉耦合阻尼系數(shù)[13].

轉(zhuǎn)子的穩(wěn)定性可以用渦動(dòng)系數(shù)來(lái)判定，渦動(dòng)系數(shù)的絕對(duì)值越小，轉(zhuǎn)子的穩(wěn)定性越高[14].定義渦動(dòng)系數(shù)Ωf為

(3)

1.2 計(jì)算模型

采用前處理軟件Gambit來(lái)建立模型和劃分網(wǎng)格，采用CFD商業(yè)軟件Ansys Fluent16.0對(duì)汽封的流場(chǎng)和氣流力進(jìn)行模擬計(jì)算.目前，汽封流體動(dòng)力特性的CFD計(jì)算方法主要包括旋轉(zhuǎn)坐標(biāo)系法(MRF)和動(dòng)網(wǎng)格法(Dynamic Mesh).由于汽封間隙較小、流場(chǎng)復(fù)雜以及湍流度較大，故采用旋轉(zhuǎn)坐標(biāo)系法，將非穩(wěn)態(tài)問(wèn)題轉(zhuǎn)化成穩(wěn)態(tài)問(wèn)題進(jìn)行求解，可避免動(dòng)網(wǎng)格法產(chǎn)生網(wǎng)格變形、畸變的問(wèn)題，保證了計(jì)算結(jié)果的準(zhǔn)確性.以氣流激振實(shí)驗(yàn)臺(tái)中的迷宮式汽封為對(duì)象，建立其數(shù)值分析模型.圖2給出了汽封二維截面示意圖.表1給出了汽封各處的尺寸.

圖2 汽封二維截面示意圖

齒數(shù)n齒根間距a/mm齒根厚度b/mm汽封間隙e/mm齒尖寬度d/mm轉(zhuǎn)子半徑r/mm汽封半徑R/mm4520.50.29097.5

1.3 網(wǎng)格劃分及邊界條件設(shè)置

計(jì)算域?yàn)槠馔鞑糠旨扒昂髩毫M(jìn)出口，圖3為汽封的全三維模型計(jì)算域示意圖.圖4給出了汽封子午面的節(jié)點(diǎn)數(shù)和網(wǎng)格圖.表2給出了子午面及周向各段的網(wǎng)格節(jié)點(diǎn)數(shù).因齒頂與轉(zhuǎn)子間隙處的流動(dòng)變化劇烈，故在此處進(jìn)行網(wǎng)格的局部細(xì)化，經(jīng)驗(yàn)證網(wǎng)格的進(jìn)一步細(xì)化對(duì)計(jì)算結(jié)果影響不大，即通過(guò)了網(wǎng)格無(wú)關(guān)性檢驗(yàn).最終得到汽封全三維模型的網(wǎng)格節(jié)點(diǎn)數(shù)約為112萬(wàn).

n1n2n3n4周向4050151048

工質(zhì)為空氣，可視為理想氣體，采用k-ε湍流模型.進(jìn)出口邊界為壓力進(jìn)出口條件，固定出口壓力為1.0×105Pa，通過(guò)改變進(jìn)口壓力來(lái)模擬不同進(jìn)出口壓力比下的工況條件.壁面均設(shè)為絕熱邊界，氣體與壁面采用光滑無(wú)滑移邊界條件.轉(zhuǎn)子繞其軸心的轉(zhuǎn)速為N，渦動(dòng)速度為Ω.采用旋轉(zhuǎn)坐標(biāo)系法(見(jiàn)圖5)，設(shè)置靜子壁面的旋轉(zhuǎn)角速度為Ω，轉(zhuǎn)子旋轉(zhuǎn)壁面繞其軸心旋轉(zhuǎn)的相對(duì)角速度為N-Ω，汽封壁面繞轉(zhuǎn)子軸心的相對(duì)角速度為-Ω.

圖5 旋轉(zhuǎn)坐標(biāo)系法

通過(guò)控制變量，分別改變不同實(shí)驗(yàn)組別的偏心率、進(jìn)出口壓力比、轉(zhuǎn)速和渦動(dòng)速度來(lái)研究汽封流體動(dòng)力特性的影響因素，從而求出密封動(dòng)力特性系數(shù).不同實(shí)驗(yàn)組別的邊界條件設(shè)置見(jiàn)表3.

表3 邊界條件

2 計(jì)算結(jié)果與分析

2.1 轉(zhuǎn)速對(duì)氣流激振特性的影響

切向力是促使轉(zhuǎn)子渦動(dòng)的氣流力來(lái)源，也是導(dǎo)致轉(zhuǎn)子失穩(wěn)的主要因素.圖6給出了不同轉(zhuǎn)速下切向力與渦動(dòng)速度之間的關(guān)系.對(duì)圖6中的數(shù)據(jù)進(jìn)行一次線(xiàn)性擬合，將擬合出的系數(shù)代入式(2)，求出交叉耦合剛度系數(shù)k和直接阻尼系數(shù)C，再將其代入式(3)中求出渦動(dòng)系數(shù).

圖7給出了渦動(dòng)系數(shù)隨轉(zhuǎn)速的變化.從圖7可以看出，隨著轉(zhuǎn)速的增大，渦動(dòng)系數(shù)增大，因此轉(zhuǎn)速越高，轉(zhuǎn)子越容易失穩(wěn).

2.2 進(jìn)出口壓力比對(duì)氣流激振特性的影響

圖8給出了渦動(dòng)系數(shù)隨進(jìn)出口壓力比的變化.從圖8可以看出，當(dāng)轉(zhuǎn)速和偏心率一定時(shí)，隨著進(jìn)出口壓力比的增大，渦動(dòng)系數(shù)增大，即轉(zhuǎn)子的穩(wěn)定性隨著進(jìn)出口壓力比的增大而下降.

2.3 偏心率對(duì)氣流激振特性和泄漏量的影響

圖9給出了不同偏心率下的渦動(dòng)系數(shù).從圖9可以看出，隨著偏心率的增大，渦動(dòng)系數(shù)增大，轉(zhuǎn)子穩(wěn)定性下降，且與進(jìn)出口壓力比和轉(zhuǎn)速相比,偏心率對(duì)轉(zhuǎn)子穩(wěn)定性的影響更為顯著.

在對(duì)偏心模型進(jìn)行計(jì)算時(shí)，還可以得到不同偏心率下的迷宮式汽封泄漏量,如圖10所示.從圖10可以看出，當(dāng)其他條件相同時(shí)，泄漏量隨著偏心率的增大而增加，但總體增加量不大，在偏心率增大300%的情況下泄漏量只增加了2.1%.

3 結(jié) 論

(1) 隨著轉(zhuǎn)子轉(zhuǎn)速、進(jìn)出口壓力比和偏心率的增大，渦動(dòng)系數(shù)增大，轉(zhuǎn)子越容易失穩(wěn)，其中偏心率對(duì)轉(zhuǎn)子穩(wěn)定性的影響更為顯著.

(2) 偏心率增大會(huì)導(dǎo)致泄漏量增加，汽封的密封效果下降.

(a) N=3 000 r/min

(b) N=6 000 r/min

(c) N=9 000 r/min

(d) N=12 000 r/min

圖7 渦動(dòng)系數(shù)隨轉(zhuǎn)速的變化

圖8 渦動(dòng)系數(shù)隨進(jìn)出口壓力比的變化

圖9 渦動(dòng)系數(shù)隨偏心率的變化

圖10 泄漏量隨偏心率的變化

[1] 潘永密,鄭水英.迷宮密封氣流激振及其動(dòng)力系數(shù)的研究[J].振動(dòng)工程學(xué)報(bào),1990,3(2):48-58.

PAN Yongmi,ZHENG Shuiying.Study on the flow induced exciting forces and the dynamic characteristics of labyrinth seals[J].Journal of Vibration Engineering,1990,3(2):48-58.

[2] LELLI D,CHEW J W,COOPER P.Combined three-dimensional fluid dynamics and mechanical modeling of brush seals[J].Journal of Turbomachinery,2006,128(1):188-195.

[3] YüCEL U.Calculation of leakage and dynamic coefficients of stepped labyrinth gas seals[J].Applied Mathematics and Computation,2004,152(2):521-533.

[4] 張強(qiáng),何立東.蜂窩密封動(dòng)力特性系數(shù)的計(jì)算方法[J].中國(guó)電機(jī)工程學(xué)報(bào),2007,27(11):98-102.

ZHANG Qiang,HE Lidong.Study on calculation of the dynamic coefficients of honeycomb seals[J].Proceedings of the CSEE,2007,27(11):98-102.

[5] 蘇華.指尖密封結(jié)構(gòu)和性能的設(shè)計(jì)分析與試驗(yàn)研究[D].西安:西北工業(yè)大學(xué),2006.

[6] CHILDS D W,GANSLE A J.Experimental leakage and rotordynamic results for helically grooved annular gas seals[J].Journal of Engineering for Gas Turbines and Power,1996,118(2):389-393.

[7] LAOS H E,VANCE J M,BUCHANAN S E.Hybrid brush pocket damper seals for turbomachinery[J].Journal of Engineering for Gas Turbines and Power,2000,122(2):330-336.

[8] 何立東,袁新,尹新.刷式密封研究的進(jìn)展[J].中國(guó)電機(jī)工程學(xué)報(bào),2001,21(12):28-32,53.

HE Lidong,YUAN Xin,YIN Xin.The advance in the investegation of brush seal[J].Proceedings of the CSEE,2001,21(12):28-32,53.

[9] 徐悅,田愛(ài)梅.基于CFD的渦輪泵轉(zhuǎn)子密封流體激振研究進(jìn)展[J].火箭推進(jìn),2005,31(1):8-13.

XU Yue,TIAN Aimei.Progress of investigation on turbopump annular seal fluid-induced vibration based on CFD[J].Journal of Rocket Propulsion,2005,31(1):8-13.

[10] 沈耀陽(yáng).迷宮密封氣流激振的數(shù)值模擬與分析[D].保定:華北電力大學(xué),2008.

[11] 駱名文,丁學(xué)俊,祁小波,等.大型汽輪機(jī)汽流激振研究現(xiàn)狀及展望[J].熱能動(dòng)力工程,2006,21(6):551-555.

LUO Mingwen,DING Xuejun,QI Xiaobo,et al.The status quo and prospect concerning a study on steam flow excitations in large-sized steam turbines[J].Journal of Engineering for Thermal Energy and Power,2006,21(6):551-555.

[12] 張萬(wàn)福.汽輪機(jī)流體激振及抑制技術(shù)理論與試驗(yàn)研究[D].南京:東南大學(xué),2013:41-43.

[13] HIRANO T,GUO Zenglin,KIRK R G.Application of computational fluid dynamics analysis for rotating machinery-part II:labyrinth seal analysis[J].Journal of Engineering for Gas Turbines and Power,2005,127(4):820-826.

[14] CHILDS D.Turbomachinery rotordynamics:phenomena,modeling,and analysis[M].New York,USA:Wiley,1993.

Research on Flow Field Dynamic Characteristics in a Labyrinth Seal

MAHaoyuan1,ZHANGWanfu1,2

(1.School of Energy and Power Engineering,University of Shanghai for Science and Technology,Shanghai 200093,China; 2.Shanghai Key Laboratory of Multiphase Flow and Heat Transfer in Power Engineering,School of Energy and Power Engineering,University of Shanghai for Science and Technology,Shanghai 200093,China)

A full three-dimensional model was set up for a labyrinth seal,in which the flow field was calculated and analyzed under different working conditions using commercial CFD software Fluent,so as to study the relations of the fluid-induced force with the rotational speed,pressure ratio and the eccentricity,etc.,and to obtain the dynamic coefficient and whirl coefficient of the seal that may influence the stability of the rotor.Results show that with the rise of rotational speed,pressure ratio and eccentricity,the whirl coefficient increases accordingly,and the stability of the rotor is hence reduced.

labyrinth seal; fluid-induced force; dynamic characteristic; numerical simulation

2016-03-11

2016-04-18

國(guó)家自然科學(xué)基金資助項(xiàng)目(11402148)

馬浩原(1995-)，男，滿(mǎn)族，河北承德人，本科，研究方向?yàn)閯?dòng)力機(jī)械強(qiáng)度與振動(dòng).張萬(wàn)福(通信作者)，男，講師，博士，電話(huà)(Tel.)：18817882508；E-mail：wfzhang@usst.edu.cn．

1674-7607(2017)01-0021-05

TK263

A 學(xué)科分類(lèi)號(hào)：470.30