, 王爽1, 胡广阳2, 张立静2 1. 东北大学 机械工程与自动化学院, 辽宁 沈阳 110819;
2. 中国航空发动机集团 航空发动机动力传输重点实验室, 辽宁 沈阳 110015
收稿日期：2018-10-22
基金项目：国家自然科学基金资助项目(51875095)；航空科学基金资助项目(20140450001)。
作者简介：闫玉涛(1970-), 男, 吉林磐石人, 东北大学副教授。
摘要：根据偏摆角圆周密封分析, 建立受力模型和有限元分析模型.通过结构分析与热-结构耦合分析, 得到了不同偏摆角下主、辅助密封面应力和变形分布规律及泄漏间隙和泄漏量变化规律.主、辅助密封面应力和变形分布均匀, 最大应力和最大变形位于凸舌处, 主密封面最大变形大于辅助密封面, 最大应力和最大变形随偏摆角增加而增大.热-结构耦合分析应力和变形与结构分析分布规律基本一致, 但均大于结构分析.主密封面间隙从轴承腔向气腔增大, 辅助密封面间隙从中心向两端增大, 接头处间隙最大.泄漏量随偏摆角增加而增大.
关键词：偏摆角圆周密封结构分析热-结构耦合分析泄漏
Effect of Deflection Angles on Structural Characteristics of Circumferential Sealing and Leakage
YAN Yu-tao1
, WANG Shuang1, HU Guang-yang2, ZHANG Li-jing2 1. School of Mechanical Engineering & Automation, Northeastern University, Shenyang 110819, China;
2. Key Laboratory of Power Transmission Technology on Aero-Engine, Aero-Engine Corporation of China, Shenyang 110015, China.
Corresponding author: YAN Yu-tao. E-mail:ytyan@mail.neu.edu.cn
Abstract: Based on the circumferential sealing analysis of deflection angles, the force model and the finite element analysis model were established. Through structural analysis and thermo-structural coupling analysis, the distribution of stress and deformation of the main and auxiliary sealing surfaces were obtained under different deflection angles, as well as the variation law of leakage gap and leakage amount. The stress and deformation of the main and auxiliary sealing surfaces are evenly distributed. The maximum stress and maximum deformation are located at the tongue. The maximum deformation of the main sealing surface is larger than that of the auxiliary sealing surface. The maximum stress and maximum deformation increase with the increasing of deflection angles. The thermal-structural coupling analysis stress and deformation are basically consistent with the structural analysis distribution, but they are larger than that in the structural analysis. The gap of the main sealing surface increases from the bearing chamber to the air chamber, and the gap of the auxiliary sealing surface increases from the center to both ends, and the gap at the joint is the largest. The amount of leakage increases with the increasing of deflection angles.
Key words: deflection anglecircumference sealingstructural analysisthermal-structural coupling analysisleakage
航空发动机密封是发动机动力传输系统的重要组成部分, 密封性能的退化、结构故障及其导致的发动机结构故障是更换发动机主要原因之一.先进密封技术能给发动机带来意想不到的收益, 且在满足发动机性能, 如燃油消耗率、推重比、寿命等方面具有十分重要的作用[1-2].航空发动机圆周密封属于接触式动密封, 具有较低的泄漏率, 在燃气涡轮发动机中取代了迷宫式密封, 能够有效限制泄漏, 降低润滑油的氧化以获得最大效率[3].文献[4-5]以低泄漏、高耐久性密封条件为目标, 对涡轮泵中的分段圆周密封环进行了实验研究, 给出了密封装置的磨损特性和考虑气体压强时的泄漏率变化规律.文献[6]通过对主轴承腔密封研究, 提出了针对提高圆周密封工作能力和耐久性的改进研究方向.圆周密封环处于高压差、高滑动速度及高温工作环境, 齿轮驱动涡扇发动机中, 由于齿轮传动系统受到惯性和陀螺力影响, 密封环会产生极高偏摆角[7].文献[7-8]针对四段石墨密封环组成的圆周密封进行了角度和径向失调试验, 发现密封环在达到径向失调极限为1.016 mm和偏摆角极限为0.5°时, 工作100 h后超过了密封环磨损极限, 并发生了不同程度失效.文献[9]研究了高压气体密封环热行为, 证明在高速旋转情况下应考虑密封环热变形.文献[10]基于ANSYS, 对机械密封稳态和瞬态热-结构耦合进行了研究, 获得了温度场和变形分布规律.文献[11-12]针对指尖密封泄漏和磨损, 通过碳/碳复合材料的改善和建立等效动力学模型对其进行了分析.
圆周密封要求具有大幅度角向偏摆动适应性能, 本文采用ANSYS WORKBENCH, 对某型在役圆周密封环工况下考虑密封环偏摆角开展结构分析、热-结构耦合分析及泄漏分析, 为圆周密封结构优化, 提升其密封性能和使用寿命提供基础.
1 考虑偏摆角圆周密封力学模型圆周密封装置主要包括密封座、密封跑道、密封环、周簧和波簧等, 密封环由沿周向等分5段结构相同石墨环组成.圆周密封为接触式动密封, 工作时密封跑道高速转动为转子, 密封环嵌在密封座中由周簧和波簧定位.由于磨损及安装偏差等, 工作中存在一定偏摆角, 主、辅助密封面与无偏摆角时受力和密封性能存在差异, 导致密封环局部受力不均, 造成更严重磨损和泄漏, 影响其密封性能和使用寿命, 考虑密封环偏摆角圆周密封受力如图 1所示.其中, P1为气腔压力, MPa; P2为轴承腔压力, MPa; R1为密封跑道外半径, mm; R2为密封座内半径, mm; R3为密封环端面密封凸缘外半径, mm; t1为气腔温度, ℃; t2为轴承腔温度, ℃; l1为主密封缝隙流动长度, mm; l2为辅助密封缝隙流动长度, mm; β为偏摆角; l′为接头处密封缝隙长度, mm; b为接头处凸缘宽度, mm; c为密封环接头间隙, mm.
图 1(Fig. 1)
 | 图 1 考虑偏摆角圆周密封环受力分析Fig.1 Force analysis of circumferential sealing ring considering deflection angles (a)—圆周密封结构；(b)—密封环接头. |
1.1 气体载荷1) 密封缝隙中气体压力形成的单位长度推开力为
 | (1) |
2) 径向气体单位长度不平衡载荷为
 | (2) |
3) 轴向气体单位长度不平衡载荷为
 | (3) |
4) 密封环接头间隙形成径向气体不平衡载荷为
 | (4) |
1.2 总接触载荷1) 密封环与密封跑道表面之间径向单位长度总接触载荷为
 | (5) |
2) 密封环与密封座端面之间轴向单位长度总接触载荷为
 | (6) |
2 考虑偏摆角圆周密封结构特性分析2.1 有限元分析模型由圆周密封结构知, 每段密封环受力情况相同, 同时考虑到圆周密封结构对称性, 取单段密封环和四分之一密封跑道与密封座建立有限元分析模型, 通过ANSYS WORKBENCH进行分析.根据工程应用和性能提升需要, 分析中偏摆角β分别取0°, 0.1°, 0.2°, 0.3°, 0.4°, 0.5°, 0.7°, 0.9°和1.0°共9种情况, 考虑偏摆角圆周密封结构如图 2所示.密封环为石墨材料, 密封跑道和密封座为金属材料, 材料属性见文献[13], 结构参数见表 1.
图 2(Fig. 2)
 | 图 2 圆周密封有限元模型Fig.2 Finite element model of circumferential sealing |
表 1(Table 1)
 表 1 结构参数Table 1 Structure parameters
| 表 1 结构参数 Table 1 Structure parameters | ||||||||||||||||||||||||||||||
压差≤0.5 MPa, 封严气体温度为320~400 ℃, 线速度为120~140 m/s, 偏摆角为0~1.0°, 周簧力为2.9~3.5 N, 波簧力为24~30 N, 摩擦因数为0.15, 满足密封环强度要求.
2.2 考虑偏摆角结构分析2.2.1 网格划分与接触建立网格划分选用尺寸控制, 设定单元尺寸为0.6 mm, 密封环采用10节点四面体SOLID187单元, 密封座与密封跑道采用20节点六面体SOLID186单元, 如图 3a所示, 该模型共有544 821个节点, 195 307个单元.
图 3(Fig. 3)
 | 图 3 网格划分与密封环接触面Fig.3 Meshing and contact surface of sealing ring (a)—网格划分；(b)—接触面. |
密封环与密封跑道和密封座均采用面-面接触方式, 选择密封环主、辅助密封面为接触面, 采用8节点高阶四边形CONTA174单元, 密封跑道外表面与密封座凸缘端面为目标面, 使用TARGE170单元, 如图 3b所示.接触类型设置为标准接触类型, 允许接触分离, 分离后法向接触力为零.
2.2.2 加载和约束根据实际工况, 在柱坐标系下对密封环施加约束和载荷, X为径向, Y为周向, Z为轴向.密封座保持固定, 密封跑道在周向无位移约束.Y方向对密封环施加位移约束, 保证Y方向密封环固定.X方向对密封环周向弹簧槽面施加固定约束, 固定密封环.密封环在X方向受周簧产生径向载荷、波簧产生径向载荷、径向气体不平衡载荷和密封环接头间隙径向气体不平衡载荷作用, 故将接头间隙处产生的径向气体不平衡载荷施加到主密封面两端, 其余施加于主密封面.密封环在Z方向受周簧产生轴向载荷、波簧产生轴向载荷、轴向气体不平衡载荷和密封环接头间隙产生的轴向气体不平衡载荷.将接头间隙产生的轴向气体不平衡载荷施加于辅助密封面两端, 其余施加到辅助密封面.
2.2.3 结果分析图 4为偏摆角β=0.1°时，密封环主、辅助密封面应力分布云图.主密封面应力基本呈均匀分布, 密封环凸舌处应力最大为3.955 1 MPa(图 4a).辅助密封面应力大部分呈均匀分布, 最大应力出现在密封环凸舌处(图 4b).分析认为:最大应力出现在主、辅助密封面交界，是由于密封环接头处的结构, 导致接头处有较大气体不平衡载荷作用, 密封环存在偏摆角后, 导致主、辅助密封面交界与密封跑道存在一定程度线面接触, 加之气体不平衡载荷作用, 致使最大应力出现在主、辅助密封面交界处.
图 4(Fig. 4)
 | 图 4 密封面应力云图Fig.4 Stress nephogram of sealing surface (a)—主密封面；(b)—辅助密封面. |
图 5为偏摆角β=0.1°时，密封环主、辅助密封面变形云图.主密封面变形基本呈均匀分布, 密封环接头处变形最大为5.742×10-3mm(图 5a).辅助密封面变形分布均匀, 密封环接头凸舌处变形最大, 其值为5.174×10-3mm(图 5b).各密封面最大变形均出现在靠近密封环接头处, 原因是由于密封环接头凸舌处较大的气体不平衡载荷作用所致.主密封面最大变形量大于辅助密封面是由于径向气体不平衡载荷大于轴向气体不平衡载荷所致.
图 5(Fig. 5)
 | 图 5 密封面变形云图Fig.5 Deformation nephogram of sealing surface (a)—主密封面；(b)—辅助密封面. |
图 6为密封环主、辅助密封面应力随偏摆角变化规律.主、辅助密封面最大应力值随偏摆角增加而增大, β=0.1°, 主密封面最大应力增大了2.1 MPa, 比β=0°时增加了1倍.辅助密封面最大应力增大了3.08 MPa, 比β=0°时增加了近3倍.但随偏摆角增大, 主、辅助密封面最大应力增大梯度逐渐变小, 由式(5)和式(6)可得偏摆角的存在导致径向总接触载荷增大, 接触载荷的增大将加剧密封环磨损和摩擦热量增加, 降低密封环使用寿命.
图 6(Fig. 6)
 | 图 6 应力随偏摆角β的变化Fig.6 Variation of stress with deflection angle β |
密封环主、辅助密封面变形随偏摆角变化规律如图 7所示.主、辅助密封面最大变形量随偏摆角增加而增大.偏摆角从0增至0.1°, 主、辅助密封面最大变形值变化显著, 分别增大了3.83 μm和4.42 μm, 偏摆角超过0.3°后, 主、辅助密封面最大变形量增大梯度减小, 辅助密封面变形量小于主密封面变形量.分析认为, 由于径向总接触载荷比无偏摆角大, 随偏摆角增大而增大; 轴向总接触载荷比无偏摆角时小, 且随偏摆角增大而减小, 所以, 存在偏摆角时, 最大变形量均大于无偏摆角情况.径向总接触载荷大于轴向总接触载荷, 导致主密封面最大变形量大于辅助密封面最大变形量.
图 7(Fig. 7)
 | 图 7 变形量随偏摆角β变化Fig.7 Variation of deformation with deflection angle β |
2.3 考虑偏摆角热-结构耦合分析温度对石墨材料特性影响显著, 在高温环境下, 考虑圆周密封结构和工作特性, 圆周密封主、辅助密封面均存在摩擦生热, 必须考虑热效应对圆周密封的密封性能及使用寿命的影响, 因此, 热-结构耦合分析非常关键[14-16].
2.3.1 稳态热分析密封环采用10节点三维实体SOLID87单元, 密封座与密封跑道选用20节点三维实体SOLID90单元.SOLID87与SOLID187都是三维10节点实体单元, SOLID90与SOLID186都是三维20节点实体单元, 所以结构和热分析单元有相同节点号, 网格划分后, 分析模型与结构分析相同, 设定材料导热系数为定值.主、辅助密封面摩擦生热量求解与分配方式见文献[13, 17].以热流密度形式分别加载到主、辅助密封面, 非接触面温度设定为338 ℃.
图 8为偏摆角β=0.1°时，圆周密封环温度场分布云图.密封环温度场分布比较均匀, 密封环接头非凸舌处温度最高, 处于主密封面, 最高温度为392.75 ℃.原因是该处接触面积相对较大, 为高速滑动摩擦，生热量大, 散热条件较差, 所以温升最高.
图 8(Fig. 8)
 | 图 8 密封环温度云图Fig.8 Temperature nephogram of sealing ring |
2.3.2 热-结构耦合分析将热分析温度场以载荷施加于热-结构耦合分析模型, SOLID87单元与SOLID187单元等效, SOLID90单元与SOLID186单元等效, 网格划分、加载和约束与结构分析相同.
图 9为偏摆角β=0.1°时，主、辅助密封面应力云图.可以看出, 主密封面应力分布比较均匀, 最大应力出现在凸舌处, 最大值为4.063 4 MPa(图 9a).辅助密封面应力分布也较均匀, 最大应力也位于凸舌处, 与主密封面最大应力值相同(图 9b).与结构分析对比发现, 考虑热应力作用, 应力场分布规律基本一致, 但是应力值增大了3 %, 这是因为在热效应作用下, 密封跑道、密封座和密封环热膨胀使接触压力增大, 导致应力值增大, 应力值的增加会加剧密封环摩擦磨损, 影响密封性能和使用寿命, 因此, 热效应对密封性能的影响不可忽视.
图 9(Fig. 9)
 | 图 9 密封面应力云图Fig.9 Stress nephogram of sealing surface (a)—主密封面；(b)—辅助密封面. |
偏摆角β=0.1°时，主、辅助密封面变形云图如图 10所示.主密封面变形分布总体比较均匀, 最大变形位于凸舌区域, 与高应力区对应, 最大值为6.382 μm(图 10a).辅助密封面变形分布较均匀, 最大变形位于凸舌区域, 变形沿径向外法线方向递减, 最大变形量为5.763 μm, 为主密封面最大变形量的90 % (图 10b).对比结构分析, 热-结构耦合分析主密封面变形量增大了11 %, 辅助密封面最大变形量为结构分析的1.11倍, 热效应对密封面变形量也存在较大影响.热效应作用下圆周密封组件发生热膨胀, 热膨胀变形和接触应力作用下变形出现耦合作用, 热效应出现导致密封环接触压力增大, 热膨胀可使得密封性能得到更好的保证, 但将增大密封环更大摩擦生热和磨损, 降低使用寿命.
图 10(Fig. 10)
 | 图 10 密封面变形云图Fig.10 Deformation nephogram of sealing surface (a)—主密封面；(b)—辅助密封面. |
3 考虑偏摆角圆周密封泄漏特性分析3.1 密封面间隙图 11为偏摆角β=0.1°时，热-结构耦合分析主、辅助密封面间隙云图.可以看出, 主密封面间隙沿周向基本封严, 但在密封面两端及接头处存在间隙, 最大间隙值为0.017 62 mm(图 11a).辅助密封面间隙在单环周向表现为从中心向两端逐渐增大变化规律, 密封接头处的间隙最大为0.026 09 mm(图 11b).
图 11(Fig. 11)
 | 图 11 密封面间隙云图Fig.11 Gap nephogram of sealing surface (a)—主密封面；(b)—辅助密封面. |
密封面间隙随偏摆角变化规律如图 12所示.主密封面与密封跑道间的间隙表现为突然增大后趋于平缓, 不考虑偏摆角的间隙为0.002 5 mm, 当偏摆角增大到0.2°时, 间隙达到0.240 8 mm, 是0.1°偏摆角的14倍, 偏摆角超过0.2°, 主密封面间隙变化趋于平缓.辅助密封面间隙随偏摆角表现为突然增大, 趋于平稳而后又近似呈线性增大, 偏摆角为0°时, 间隙为2.996×10-5 mm, 当偏摆角达到0.2°时, 间隙达到了0.240 7 mm, 是0.1°偏摆角下的9倍.圆周密封最大偏摆角要求不超过0.5°, 可以看出, 在0.2°到0.5°之间, 辅助密封面间隙变化放缓, 偏摆角超过0.2°后, 偏摆角对辅助密封面间隙影响非常显著, 而对主密封面间隙影响较小.主、辅助密封面间隙变化规律与文献[8]实验研究规律基本一致, 说明本文分析结果是正确有效的, 为圆周密封设计和性能提升提供了基础.
图 12(Fig. 12)
 | 图 12 密封面间隙随偏摆角β变化Fig.12 Variation of gap with deflection angle β |
3.2 圆周密封泄漏量圆周密封总泄漏量由两部分泄漏量组成[18], 主密封面泄漏量G1可表示为
 | (7) |
 | (8) |
 | (9) |
 | (10) |
总泄漏量G为
 | (11) |
图 13(Fig. 13)
 | 图 13 泄漏量随偏摆角β的变化Fig.13 Variation of leakage with deflection angle β |
4 结论1) 通过圆周密封分析, 建立了考虑密封环偏摆角圆周密封力分析模型和计算模型.
2) 通过考虑偏摆角圆周密封结构分析, 得到了主、辅助密封面应力和变形分布及最大应力和最大变形随偏摆角变化规律.主、辅助密封面应力和变形分布较均匀, 最大应力和最大变形出现在凸舌处, 主密封面最大变形大于辅助密封面变形.最大应力和最大变形随偏摆角增加而增大.
3) 由考虑偏摆角圆周密封热-结构耦合分析得到了主、辅助密封面应力与变形分布, 应力和变形分布比较均匀, 最大应力和最大变形出现在凸舌处, 与结构分析相比, 最大应力和最大变形均显著增大, 热效应对圆周密封性能有显著影响.
4) 通过考虑偏摆角圆周密封泄漏特性, 得到了主、辅助密封面间隙分布规律, 主密封面间隙随偏摆角增加先增大后平稳变化, 辅助密封面间隙先增大, 再趋于平缓, 而后又呈线性增加的变化规律, 泄漏量随偏摆角度的增加而增大.
参考文献
| [1] | Steinetz B M, Hendricks R C.Advanced seal technology role in meeting next generation turbine engine goals[R].Cleveland: National Aeronautics and Space Administration Lewis Research Center, 1998. |
| [2] | Steinetz B M, Hendricks R C. Engine seal technology requirements to meet NASA′s advanced subsonic technology program goals[J].Journal of Propulsion & Power, 2015, 12(4): 786–793.DOI:10.2514/3.24102 |
| [3] | Swikert M, Johnson R. Wear of carbon-type seal materials with varied graphite content[J].ASLE Transactions, 1958, 1(1): 115–120. |
| [4] | Oike M, Nosaka M, Watanabe Y, et al. Experimental study on high-pressure gas seals for a liquid oxygen turbopump[J].Tribology Transactions, 1988, 31(1): 91–97.DOI:10.1080/10402008808981803 |
| [5] | Kikuchi M, Oike M, Nosaka M, et al.Durability of a carbon segmented circumferential seal for a liquid oxygen turbopump[R].Tokyo: Technical Report of National Aerospace Laboratory, 1992: 1-15. |
| [6] | 林基恕, 张振波. 21世纪航空发动机动力传输系统的展望[J].航空动力学报, 2001, 16(2): 108–114. ( Lin Ji-shu, Zhang Zhen-bo. Prospects of aeroengine power transmission system in the 21st century[J].Journal of Aerospace Power, 2001, 16(2): 108–114.DOI:10.3969/j.issn.1000-8055.2001.02.003) |
| [7] | Shaughnessy D, Dobek L.High misalignment carbon seals for the fan drive gear system technologies[R].Hanover: NASA Center for Aerospace Information, 2006. |
| [8] | Shaughnessy D, Dobek L.High misalignment carbon seals for the fan drive gear system technologies[R].Hanover: NASA Center for Aerospace Information, 2005. |
| [9] | Suganami T, Masuda T, Oishi N, et al. A study on thermal behavior of large seal-ring[J].Journal of Tribology, 1982, 104(4): 449–454.DOI:10.1115/1.3253248 |
| [10] | Liu L, Liu Z L, Wang J, et al. Coupled thermal-structural analysis of mechanical seal[J].Advanced Materials Research, 2012, 479: 1110–1114. |
| [11] | 路菲, 陈国定, 苏华, 等. 2.5D碳/碳复合材料指尖密封动态性能分析[J].航空学报, 2013, 34(11): 2616–2625. ( Lu Fei, Chen Guo-ding, Su Hua, et al. Analysis of dynamic performance for 2.5D carbon-carbon composite finger seals[J].Acta Aeronautica et Astronautica Sinica, 2013, 34(11): 2616–2625.) |
| [12] | 马晓林, 陈国定. 基于等效模型的流体动压指尖密封动力学及泄漏量分析[J].航空学报, 2008, 29(5): 1356–1363. ( Ma Xiao-lin, Chen Guo-ding. Dynamics and leakage analysis of padded finger seal based on equivalent model[J].Acta Aeronautica et Astronautica Sinica, 2008, 29(5): 1356–1363.DOI:10.3321/j.issn:1000-6893.2008.05.042) |
| [13] | 闫玉涛, 张博, 胡广阳, 等. 石墨圆周密封热-结构耦合分析[J].航空动力学报, 2018, 33(2): 273–281. ( Yan Yu-tao, Zhang Bo, Hu Guang-yang, et al. Analysis of thermal-structure coupling for graphite circumferential seal[J].Journal of Aerospace Power, 2018, 33(2): 273–281.) |
| [14] | 闫玉涛, 李雪娟, 胡广阳, 等. 石墨密封材料高温摩擦磨损行为及预测[J].航空动力学报, 2014, 29(2): 314–320. ( Yan Yu-tao, Li Xue-juan, Hu Guang-yang, et al. Friction/wear behaviors and predication of graphite seal material under high temperature[J].Journal of Aerospace Power, 2014, 29(2): 314–320.) |
| [15] | Soltani R, Heydarzadeh S M, Ansari M, et al. Effect of APS process parameters on high-temperature wear behavior of nickel-graphite abradable seal coatings[J].Surface & Coatings Technology, 2017, 321: 403–408. |
| [16] | Rajaram G, Kumaran S, Rao T S, et al. Studies on high temperature wear and its mechanism of Al-Si/graphite composite under dry sliding conditions[J].Tribology International, 2010, 43(11): 2152–2158.DOI:10.1016/j.triboint.2010.06.004 |
| [17] | B布尚.摩擦学导论[M].葛世荣, 翻译译.北京: 机械工业出版社, 2007. ( Bhushan Bharat.Introduction to tribology[M].Translated by Ge Shin-rong.Beijing: China Machine Press, 2007.) |
| [18] | 王伟, 赵宗坚, 张振生. 单环圆周密封装置设计和应用研究[J].航空发动机, 2009, 35(4): 6–11. ( Wang Wei, Zhao Zong-jian, Zhang Zhen-sheng. Design and application of single circumferential seal[J].Aeroengine, 2009, 35(4): 6–11.) |
