陈尧兴, 李志刚, 晏鑫, 李军,2
(1.西安交通大学能源与动力工程学院, 710049, 西安; 2.先进航空发动机协同创新中心, 100191, 北京)
迷宫齿蘑菇形磨损时密封泄漏特性和转子动力特性系数研究
陈尧兴1, 李志刚1, 晏鑫1, 李军1,2
(1.西安交通大学能源与动力工程学院, 710049, 西安; 2.先进航空发动机协同创新中心, 100191, 北京)
采用基于转子多频椭圆涡动模型和动网格技术的Unsteady Reynolds-Averaged Navier-Stokes (URANS)方程求解方法,研究了迷宫齿蘑菇形磨损对密封泄漏特性和气流激振转子动力特性的影响,计算分析了未磨损以及磨损间隙分别为0.4、0.5、0.6 mm时的密封泄漏量,流场以及转子动力特性系数。结果表明:迷宫齿蘑菇形磨损使得密封间隙内缩流面积增加,导致迷宫密封泄漏量随磨损间隙增加而增加,且当磨损间隙大于0.4 mm时,泄漏量随磨损间隙增加而线性增加;迷宫齿蘑菇形磨损增加了密封的直接刚度,降低了密封交叉刚度的大小以及直接阻尼,但当磨损间隙超过0.5 mm时,直接阻尼不再发生改变;随着迷宫齿蘑菇形磨损间隙增加,迷宫密封有效阻尼降低,从而使得转子稳定性降低,但仍然处于稳定范围内。
迷宫密封;迷宫齿蘑菇形磨损;泄漏特性,转子动力特性;数值模拟
通常在透平机组稳定运行情况下,迷宫密封在其设计密封间隙下具有良好的封严性能,然而当透平机组运行在非稳定工况下,如开机、停机以及热重启阶段,迷宫密封静子件容易与转子件之间发生碰磨从而使得密封间隙增加,封严性能退化[1],同时引起密封气流激振转子动力特性发生改变。另一方面,迷宫齿的磨损形状主要与迷宫齿和其碰撞面的材料硬度有关。当迷宫齿材料硬度小于转子材料硬度时,迷宫齿在碰磨过程中则会形成蘑菇形结构。Ghasripoor等公开了实际运行过程中迷宫齿由于齿间磨损形成的蘑菇形迷宫齿结构[1],如图1所示。
图1 迷宫齿蘑菇形磨损的实际结构图[1]
目前,关于密封磨损的研究主要集中在泄漏和传热特性方面。Xu等研究了迷宫齿弯曲、蘑菇形磨损,以及蘑菇形磨损和弯曲同时存在的情况下密封的泄漏特性,指出迷宫齿磨损或者弯曲能显著增加密封泄漏量、降低透平效率[2]。Dogu等从几何结构(迷宫齿蘑菇形半径、密封间隙)和运行参数(压比、转速、齿数)方面对迷宫密封蘑菇形磨损后的泄漏特性进行了详细的研究[3]。Rhode等指出密封静子面磨损形成凹槽后密封泄漏量增加的原因是,迷宫齿间隙内有效间隙增加以及迷宫齿附近流场发生改变[4]。Yan等采用流固耦合传热的方法对迷宫齿蘑菇形磨损后密封传热特性进行了研究,指出静子面与转子面上的努塞尔数分布随迷宫齿蘑菇形半径增加而减小[5]。
然而,关于迷宫齿磨损后密封转子动力特性的研究较少。Childs等以压气机中典型层流缓冲油封为研究对象,开展了密封静子面磨损形成凹槽后的油封结构在不同偏心率下的转子动力特性研究,并指出当偏心率大于0.5时,静子面磨损后的油封结构稳定性低于光滑油封结构,但当偏心率低于0.3时,静子面磨损后的油封结构具有明显更低的交叉刚度[6]。Wang等采用单控制体以及有限扰动法对密封静子面凹槽磨损后的转子动力特性进行了研究,并分析了迷宫齿与凹槽中心的位置关系对转子动力特性的影响,指出当转子面上迷宫齿位于凹槽上游时,密封具有最大的直接阻尼[7]。
因此,本文以Ertas等的迷宫密封实验模型[8]为研究对象,开展迷宫齿蘑菇形磨损时密封泄漏特性和转子动力特性系数的研究,为迷宫密封磨损状态下的性能分析提供参考。
图2 迷宫齿未磨损时迷宫密封三维计算模型
图3 迷宫齿蘑菇形磨损前后密封几何尺寸
图2给出了迷宫齿未磨损时6齿直通式迷宫密封的三维计算模型,其中进出口部分通过加长处理以避免回流的产生。图3给出了迷宫齿间距、腔室高度以及磨损前后迷宫齿4种不同的齿顶尺寸,其中迷宫齿齿顶蘑菇形磨损采用等截面设计,即认为迷宫齿齿顶磨损前后密封质量或者密封横截面面积不变。图4给出了迷宫齿蘑菇形磨损时密封三维网格和齿顶局部网格图。表1给出了迷宫齿磨损前后部分几何参数和计算边界条件。
图4 迷宫齿蘑菇形磨损后密封三维网格示意图
参数数值蘑菇形半径Rm/mm0,0.15,0.19,0.23迷宫齿磨损间隙Cr/mm0.3,0.4,0.5,0.6迷宫齿数ni6转子直径D0/mm170进口总压Pt/kPa699.14进口总温Tt/℃14出口静压Ps/kPa388.07转速n/r·min-115000
采用基于动网格技术和多频椭圆涡动模型的非定常数值方法求解Unsteady Reynolds-Averaged Navier-Stokes (URANS)方程,计算分析了3种迷宫齿蘑菇形磨损结构下密封泄漏量、流场、转子动力特性系数、旋流强度、转子面气流激振力,并与未磨损状态的迷宫密封特性进行对比。表2给出了迷宫密封转子动力特性系数计算的数值方法和涡动模型参数。图5给出了转子椭圆涡动轨迹模型示意图。
表2 转子动力特性系数计算方法以及涡动模型参数
注:a和b分别为椭圆轨迹的长轴与短轴。
(a)x方向激励 (b)y方向激励图5 转子多频椭圆涡动轨迹模型示意图
根据小位移涡动理论可知,密封内力-位移方程满足如下关系
(1)
对x激励、y激励下力-位移方程进行快速傅里叶变换,并同时求解可得方程(1)的频率解
(2)
(3)
(4)
(5)
其中:Kij=Re(Hij);Cij=Im(Hij)/Ω,Ω为涡动频率。
图6给出了本文采用的多频椭圆涡动模型预测方法[9]下获得的迷宫密封转子动力特性系数,并与实验结果[8]进行比较。多频椭圆涡动模型预测的迷宫密封转子动力特性系数与实验数据整体吻合良好,其中刚度系数以及直接阻尼的预测结果较实验值偏低,当涡动频率小于120 Hz时,有效阻尼预测结果较实验值偏高,当涡动频率大于120 Hz时,有效阻尼的预测结果与实验值一致。
考虑到迷宫密封转子稳定性主要与切向力系数有关,因此选取直接阻抗Hxx的虚部ΩCxx以及交叉阻抗Hxy的实部Kxy进行网格无关性验证。图7给出了磨损间隙为0.6 mm时两种疏密网格下Im(Hxx)以及Re(Hxy)随涡动频率的变化曲线,其中疏网格和密网格的网格节点数分别是548.8万和823.6万。当网格节点数为548.8万时,密封切向力系数几乎不再变化;另外,疏网格与密网格计算得到的迷宫密封泄漏量分别为0.417 5和0.417 3 kg/s,表明当网格节点数为548.8万时,泄漏量基本不变。因此,综合考虑网格疏密程度对密封切向力系数、泄漏量的影响,本文迷宫密封的网格节点数选为548.8万。
(a)刚度系数
(b)阻尼系数图6 转子动力特性系数预测值与实验值比较
图7 不同网格数下迷宫密封直接阻抗虚部Im(Hxx)与交叉阻抗实部Re(Hxy)随涡动频率的变化曲线
图8给出了迷宫密封泄漏量随迷宫齿蘑菇形磨损间隙的变化曲线。可以发现:当磨损间隙大于0.4 mm时,密封泄漏量随磨损间隙增加而线性增加;另一方面,当迷宫齿磨损间隙分别从0.3 mm增加至0.4 mm、0.4 mm增加至0.5 mm、0.5 mm增加至0.6 mm时,泄漏量增幅分别为未磨损时的0.77倍、0.64倍、0.63倍,这种迷宫密封泄漏量的增幅差异主要与齿间隙内缩流面积和轴向速度变化有关。
图8 迷宫齿蘑菇形磨损间隙下密封泄漏量的变化曲线
图9给出了迷宫齿蘑菇形磨损前后第1个迷宫齿间隙内速度矢量分布。迷宫齿磨损前后间隙内缩流面积主要受到迷宫齿磨损间隙和流道压缩比(迷宫齿磨损间隙面积与缩流面积的比值)的影响。另一方面,气流通过磨损后的蘑菇形迷宫齿时接近齿顶表面流动,流道压缩比接近于1,从而使得磨损后密封的流道压缩比较未磨损时明显偏低。因此,当迷宫齿初始磨损时,即磨损间隙从0.3 mm增加至0.4 mm时,缩流面积受到磨损间隙增加和流道压缩比降低的共同影响;当磨损间隙从0.4 mm增加至0.5 mm或者从0.5 mm增加至0.6 mm时,流道压缩比不再存在明显差异,缩流面积主要受迷宫齿磨损间隙增加的影响。
图9 迷宫齿蘑菇形磨损前后第1个迷宫齿间隙内速度矢量分布
图10给出了不同迷宫齿蘑菇形磨损间隙下密封腔室内气流马赫数沿轴向分布。可以发现:马赫数沿气流流动方向逐渐增加;迷宫齿蘑菇形磨损增强了迷宫齿处节流加速过程,且随着磨损间隙增加,密封腔室内气流马赫数明显增加。值得注意的是,当磨损间隙为0.5 mm时,最后一个迷宫齿内气流最大马赫数为1.0,表明此时气流流动已经处于临界状态。
图10 迷宫齿蘑菇形磨损间隙下密封内马赫数沿轴向分布
图11给出了不同迷宫齿蘑菇形磨损间隙下密封腔室内气流静压沿轴向分布。静压沿气流流动方向呈梯形降低;随着迷宫齿蘑菇形磨损间隙增加,密封腔室内压力略微降低。需要指出的是,迷宫齿蘑菇形磨损改变了迷宫齿前后压差的大小,尤其是提高了第1个迷宫齿和最后1个迷宫齿前后压差,这是由于第1个迷宫齿和最后1个迷宫齿齿后更多的压力能转化为动能引起的。
图11 迷宫齿蘑菇形磨损间隙下密封内静压沿轴向分布
迷宫密封直接刚度主要与转子系统临界转速有关,正直接刚度会增加转子系统临界转速,负直接刚度则会降低转子系统临界转速。图12给出了不同迷宫齿蘑菇形磨损间隙下密封直接刚度随涡动频率的变化曲线。
图12 迷宫密封直接刚度随涡动频率的变化曲线
当迷宫齿蘑菇形磨损间隙从0.3 mm增加至0.4 mm时,密封直接刚度提高87.8%以上,而当磨损间隙从0.4 mm继续增加时,密封直接刚度变化幅度相对较小,其中,当磨损间隙从0.4 mm增加至0.5 mm时,密封直接刚度反而降低6.0%~14.9%,当磨损间隙从0.5 mm继续增加至0.6 mm时,密封直接刚度增加6.4%~18.3%。这种现象与文献[10]中密封间隙对直接刚度影响的变化趋势不同,主要是因为转子面上径向力受到迷宫齿蘑菇形磨损结构的影响。
迷宫密封交叉刚度与转子失稳的能力强弱有关,交叉刚度为正值表明易诱发转子失稳,交叉刚度为负值表明能够抑制转子失稳。图13给出了不同迷宫齿蘑菇形磨损间隙下密封交叉刚度随涡动频率的变化曲线。4种密封磨损结构下交叉刚度均呈现出明显的频率无关性;交叉刚度Kxy均为负值,且随着迷宫齿蘑菇形磨损间隙的增加,交叉刚度Kxy的大小逐渐减小,表明抑制转子失稳的能力同样随之减弱。迷宫密封交叉刚度随磨损间隙变化的趋势与Li等的研究结论[10]相符。
图13 迷宫密封交叉刚度随涡动频率的变化曲线
迷宫密封直接阻尼同样表征转子失稳的能力强弱,直接阻尼为正值表明能够抑制转子失稳,直接阻尼为负值则易诱发转子失稳。图14给出了不同迷宫齿蘑菇形磨损间隙下密封直接阻尼随涡动频率的变化曲线。4种磨损间隙下密封直接阻尼均随涡动频率增加而略微增加;直接阻尼均为正值,但随磨损间隙增加而略微减小,当磨损间隙从0.3 mm增加至0.5 mm时,密封直接阻尼仅降低11.1%~14.8%,而当磨损间隙继续增加至0.6 mm时,密封直接阻尼几乎不发生变化。Benckert等认为,迷宫密封直接阻尼与密封腔室内平均密度近似成线性相关[11]。因此,迷宫齿蘑菇形磨损使得腔室内压力略微降低,从而使得腔室内平均密度降低,进而引发直接阻尼随磨损间隙增加而略微降低的现象。
图14 迷宫密封直接阻尼随涡动频率的变化曲线
图15给出了不同迷宫齿蘑菇形磨损间隙下密封交叉阻尼随涡动频率的变化曲线。迷宫密封交叉阻尼存在明显的频率相关性;交叉阻尼值大小整体上随密封磨损间隙增加而逐渐减小。值得注意的是,当迷宫齿蘑菇形磨损间隙大于0.4 mm时,x方向交叉阻尼与y方向交叉阻尼的大小不再满足近似相等关系。
图15 迷宫密封交叉阻尼随涡动频率的变化曲线
迷宫密封的有效阻尼定义为
Ceff=Cxx-Kxy/Ω
(6)
综合考虑交叉刚度与直接阻尼对转子稳定性的影响,图16给出了不同迷宫齿蘑菇形磨损间隙下密封有效阻尼随涡动频率的变化曲线。4种迷宫齿蘑菇形磨损间隙下密封有效阻尼随转子涡动频率增加而逐渐减小,直至趋于稳定值;迷宫密封有效阻尼均为正值,但随着磨损间隙增加,不同涡动频率下迷宫密封有效阻尼均明显降低,表明迷宫密封转子稳定性随着迷宫齿蘑菇形磨损间隙增加而降低,但仍然处于稳定范围内。
图16 迷宫密封有效阻尼随涡动频率的变化曲线
图17给出了不同迷宫齿蘑菇形磨损间隙下密封内周向旋流强度沿轴向的分布。周向旋流强度定义为
β=60Ugas/(2πRn)
(7)
式中:Ugas是气流绝对周向速度。
由图17可知,由于转子面附近的周向黏性剪切力的作用,4种磨损间隙下迷宫密封周向旋流强度沿轴向逐渐增强;随着磨损间隙增加,迷宫密封周向旋流强度逐渐减弱,使得密封腔室内周向不平衡压力减弱,从而导致密封交叉刚度减小。
图17 迷宫齿蘑菇形磨损间隙下密封周向旋流强度沿轴向分布
表3给出了非定常计算T=0.1 s(此时x方向涡动位移最大、y方向涡动位移为0)时,不同迷宫齿蘑菇形磨损间隙下密封转子面上的径向力和切向力。可以发现,不同磨损间隙下转子径向力均为负值,即此时转子面径向力方向与转子半径增加方向相反,表明迷宫密封具有正的直接刚度(交叉阻尼过小,忽略交叉阻尼对径向力的影响);另一方面,除磨损间隙为0.4 mm以外,转子面径向力的大小整体随迷宫齿蘑菇形磨损间隙增加而逐渐增加,这一变化趋势与直接刚度随磨损间隙的变化趋势吻合一致。
表3 不同磨损间隙下转子面激振力径向分量、 切向分量
不同迷宫齿蘑菇形磨损间隙下转子切向力均为负值,即此时转子面切向力方向与转子涡动方向相反,表明迷宫密封有效阻尼为正值;转子切向力随磨损间隙增加而减小,表明迷宫密封有效阻尼随磨损间隙增加而减小,这与之前关于迷宫密封有效阻尼的分析一致。
本文数值研究了3种不同迷宫齿蘑菇形磨损间隙下密封泄漏和转子动力特性系数,并与迷宫齿未磨损时的特性进行了对比,得到如下结论。
(1)迷宫齿蘑菇形磨损使得密封间隙内缩流面积增加,且当磨损间隙大于0.4 mm时,泄漏量随迷宫齿蘑菇形磨损间隙增加而线性增加。
(2)迷宫齿蘑菇形磨损增加了密封的直接刚度,降低了密封交叉刚度以及直接阻尼,但当迷宫齿蘑菇形磨损间隙大于0.5 mm时,密封直接阻尼基本不再发生改变。此外,随着迷宫齿蘑菇形磨损间隙的增加,密封有效阻尼降低,从而转子稳定性降低,但仍然处于稳定范围内。
[1] GHASRIPOOR F, TURNQUIST N A. Wear prediction of strip seals through conductance [C]∥Proceedings of ASME Turbo Expo 2004. New York, USA: ASME, 2004: 331-337.
[2] XU J. Effects of operating damage of labyrinth seal on seal leakage and wheelspace hot gas ingress [D]. College Station, Texas, USA: Texas A&M University, 2006.
[3] DOGU Y, SERTÇAKAN M C, BAHAR A S, et al. Computational fluid dynamics investigation of labyrinth seal leakage performance depending on mushroom-shaped tooth wear [J]. ASME Journal of Engineering for Gas Turbines and Power, 2016, 138(3): 032503.
[4] RHODE D L, ADAMS R G. Computed effect of rub-groove size on stepped labyrinth seal performance [J]. Tribology Transactions, 2001, 44(4): 523-532.
[5] YAN X, LEI L, LI J, et al. Effect of bending and mushrooming damages on heat transfer characteristic in labyrinth seals [J]. ASME Journal of Engineering for Gas Turbines and Power, 2013, 136(4): 041901.
[6] CHILDS D W, RODRIGUEZ L E, CULLOTTA V, et al. Rotordynamic-coefficients and static(equilibrium loci and leakage) characteristics for short, la minar-flow annular seals [J]. ASME Journal of Tribology, 2006, 128(2): 378-387.
[7] WANG W Z, LIU Y Z, MENG G, et al. Influence of rub groove on rotordynamics associated with leakage air flow through a labyrinth seal [J]. Journal of Mechanical Science and Technology, 2010, 24(8): 1573-1581.
[8] ERTAS B H, DELGADO A, VANNINI G. Rotordynamic force coefficients for three types of annular gas seals with inlet preswirl and high differential pressure ratio [J]. ASME Journal of Engineering for Gas Turbines and Power, 2012, 134(4): 042503.
[9] LI Z, LI J, YAN X. Multiple frequencies elliptical whirling orbit model and transient RANS solution approach to rotordynamic coefficients of annual gas seals prediction [J]. Journal of Vibration & Acoustics, 2013, 135: 031005.
[10] LI Z, LI J, FENG Z. Labyrinth seal rotordynamic characteristics: part II Geometrical parameter effects [J]. Journal of Propulsion & Power, 2016, 32(5): 1-12.
[11] 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. College Station, Texas, USA: Texas A&M University, 1980: 189-212.
InvestigationontheLeakagePerformanceandRotordynamicCoefficientsofLabyrinthSealwithMushroom-ShapedToothWear
CHEN Yaoxing1, LI Zhigang1, YAN Xin1, LI Jun1,2
(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)
The influence of mushroom-shaped tooth wear on leakage performance and fluid excited rotordynamic characteristics of labyrinth seal was numerically investigated by means of the unsteady Reynolds-averaged Navier-Stokes (URANS) equation based on the multi-frequency elliptical orbit rotor whirling model and dynamic mesh technique. The leakage performance, flow fields and rotordynamic coefficients of labyrinth seal with unworn clearance of 0.3 mm and after-damage clearances of 0.4, 0.5 and 0.6 mm were calculated. The obtained results indicate that a mushroom-shaped tooth wear causes an increase of vena contracta area in sealing clearance, thus leakage flow rate is increased with the worn clearance and a linear increase is expected when after-damage clearance is over 0.4 mm. In addition, a mushroom-shaped tooth wear results in an obvious increase in direct stiffness and a drop of cross-coupling stiffness and direct damping, but direct damping will be no longer affected by worn clearance when the worn clearance is above 0.5 mm. Also, the stability is reduced as a result of a drop in effective damping with the increase of worn clearance, but is still within a stable range.
labyrinth seal; mushroom-shaped tooth wear; leakage performance; rotordynamic characteristics; numerical simulation
2017-06-07。 作者简介: 陈尧兴(1992—),男,博士生;李军(通信作者),男,教授,博士生导师。 基金项目: 国家自然科学基金资助项目(51376144)。
时间: 2017-10-18
网络出版地址: http:∥kns.cnki.net/kcms/detail/61.1069.T.20171018.1438.004.html
10.7652/xjtuxb201801007
V231.1
A
0253-987X(2018)01-0040-07
(编辑 荆树蓉 苗凌)