液体火箭发动机涡轮泵机械密封磨损机理研究

赵伟刚，张鹏鹏，任姗姗，董光能

（1.西北工业大学航天学院，陕西西安710072；2.西安交通大学现代设计及转子轴承系统教育部重点实验室，陕西西安710049）

液体火箭发动机涡轮泵机械密封磨损机理研究

赵伟刚1，张鹏鹏2，任姗姗2，董光能2

（1.西北工业大学航天学院，陕西西安710072；2.西安交通大学现代设计及转子轴承系统教育部重点实验室，陕西西安710049）

对比分析了机械密封静环端面原始表面和磨损表面形貌，得到了磨损形式和内外径磨损差异；仿真分析了密封端面接触应力及温度场，分析了转速、介质压力等条件对端面接触应力和温度场的影响；探明了端面接触应力和温度变化对磨损性能的影响，阐述了机械密封的磨损机理。

表面形貌；仿真分析；机械密封；磨损机理

0 引言

当前，接触型机械密封是液体火箭发动机涡轮泵中轴端密封的常见形式[1－2]，它在我国的航天工业中已经得到了成功的应用，具有高振动、高压力、高转速及密封介质特殊的特点[3－4]。密封环配对材料多为软-硬组合，且在特殊介质环境和工况条件下，密封端面易出现磨损[5]。表面形貌、密封环变形及温度场是影响机械密封环磨损的重要因素[6-8]。Shimomura和Hirabayashi通过实验探讨了密封端面的摩擦磨损状况随表面形貌参数变化的规律[9]。Mayer按照摩擦系数对摩擦类型进行分类，包括边界摩擦、混合摩擦及流体摩擦等[10]。顾永泉认为机械密封的其他参数通过临界工况参数来确定，有利于实现机械密封的低摩擦磨损运行[11]。国际上对密封环端面变形的研究较早，1999年Baheti运用有限元分析了机械密封的变形，认为刚度和介质压力影响密封环的变形[12]。Lebeck分析了各种载荷对端面变形的影响，比较了几种计算端面变形的方法；国内许多学者对机械密封温度场及热变形开展了大量的研究工作[13-15]。

目前针对密封磨损的研究较多，但对特殊介质下的机械密封磨损机理的研究较少。随着我国对长寿命、高可靠性及可重复使用的先进推进系统研究目标的提出，接触式机械密封要适应高承载能力、高工作转速及高工作寿命的要求。亟需研究特殊介质条件下的磨损机理，将为高可靠性、高稳定性的泵用机械动密封的设计提供基础，对于提高航天安全可靠性也具有重要意义。

1 有限元模型的建立与实验过程

1.1 有限元基本模型的建立

本文以低温推进剂所用接触型机械密封动静环为研究对象。利用ANSYS建立轴对称二维实体模型，并在静环和动环接触区域 (S1S2)建立接触模型，如图1所示。

1.2 接触应力分析边界条件的确定

1.2.1 力学分析

作用在密封环上的力可以分为2种：密封面闭合力和密封面开启力。受力分析如图2所示。

图2中：p0为介质压力；Fsp为弹簧载荷；Ff为辅助密封圈与静环之间的摩擦力；F0为轴承对动环的支撑力；pc为密封端面接触比压；pe为密封介质压力；pm为平均液膜压力。

1.2.2 参数计算及边界条件的确定

在启动和停车过程中，随着泵转速增大，介质压力随之增大，可近似线性变化。转速达到10 000 r/min左右时，动静环分离；在启动至分离的过程中介质压力的变化范围为0.3～2.5 MPa。密封环材料参数如表1所示，转速和介质压力对应数值如表2所示。

根据受力分析，确定位移约束及载荷等边界条件。对应图1中模型，在S8S9施加x方向零位移边界条件，在S7S8添加弹簧支撑力，在S6S7，S4S5，S3S4，S2S3，S2R2及R2R4施加介质压力，在R2R4施加y方向零位移边界条件。

表1 密封副材料的基本参数Tab.1 Basic parameters of seal pair materials

表2 启动过程中的转速及介质压力数值Tab.2 Values of rotating speed and medium pressure in start process

1.3 温度场有限元分析边界条件的确定

1.3.1 基本假设及微分方程

机械密封工况复杂，模拟过程中难以考虑到所有因素，需做出基本假设：

1）密封副是轴对称结构，密封副材料均匀、各向同性，物理性能不受力和温度的影响；

2） 接触副摩擦产生的热量全部在接触副传递，热流密度的分布均匀；搅拌产生的热量较小，忽略；

3） 运行时，摩擦系数不变，材料磨损的影响不计，不计少量被泄漏带走的摩擦热；

4） 不考虑材料的特性随温度变化的情况，材料的密度、导热系数、热膨胀系数及弹性模量等不变。

机械密封环的热量传递简化为二维轴对称稳态热传导，热传导微分方程为：

式中：T为物体的温度；x为径向坐标；y为轴向坐标。

1.3.2 计算参数及边界条件的确定

1.3.2.1 计算参数

密封介质的物理性能参数如表3所示。

表3 介质的基本参数Tab.3 Basic paameters of medium

1.3.2.2 摩擦热的计算与分配

对于混合摩擦机械密封，摩擦热的计算公式如下：

式中：Q为摩擦产生的总热量；f为摩擦系数；A为端面接触面积；p为端面压力；v为平均线速度；t为摩擦接触时间。

则热流密度的计算公式为：

式中：Dm为密封环的平均内径；n为转速。

本文按照Golubiev推导的计算公式进行计算[16]。对密封端面有：

式中：η为摩擦热分配系数；qs为静环端面摩擦热；qw为动环端面摩擦热；q为端面产生的总摩擦热，λ为导热系数；h为轴向厚度。

1.3.2.3 对流换热系数的确定

密封环外边界与密封介质接触，密封环内边界与大气接触。针对不同的边界有不同的半经验计算公式。对于动环外边界的对流换热系数α计算公式如下[17]：

其中

式中：Nu为努赛尔数；λ为流体的导热系数；Dr为动环外径；Rec，Rea分别为反映介质搅拌影响和横向绕流影响的雷诺数；Pr为普朗特数；μ为介质的动力粘度；cp为介质的比热容；ω为主轴角速度；U为动环周围介质的轴向流速；v为介质的运动粘度；ρ为介质密度。

对于动环内边界和静环外边界，对流换热系数的计算公式如下：

式中：Re=2Uδ/v；δ为动环与轴 (或轴套)之间或静环与密封腔内壁之间的间隙；ε为修正系数，一般取1.2～2.0；U为静环或动环周围介质的轴向流速。

对于静环内边界，对流换热系数：

1.3.2.4 边界条件的确定

根据计算的参数，确定的温度场有限元分析边界条件。对应图1中模型，在S1S2添加接触区域热流密度 (Heat Flux)，在S2S4S6添加静环外边界热对流 (Convection)，在S10S11添加静环内边界热对流 (Convection)，在S2R2R4动环外边界热对流，在S1R1动环内边界热对流 (Convection)。

1.4 静环磨损形貌分析实验

静环工作面材料为铜基石墨材料。取原始静环工件与静环磨损工件，利用MM-2金相显微镜，按照从密封环外径到内径的顺序分别对两静环工作面进行表面形貌分析，分析对比表面形貌差异，进而验证仿真分析结果。

2 结果与分析

2.1 摩擦副接触应力仿真结果分析

对不同转速和介质压力下的模型进行接触应力分析，选择接触区域，提取接触区域线段上各介质压力条件下的接触区域应力数值，并建立如图3（a）所示的坐标系，x坐标代表距离S，y坐标代表接触应力σ。输出的应力值如图3（b）所示。

图3（b）中，接触应力由内径向外径不断增大，应力集中在外径区域；转速越高，介质压力越大，密封端面外径区域的接触应力越大。应力突变值发生在距离原点4.0 mm附近处，静环在运行时存在偏载。

在10 000 r/min条件下，接触区域外径处应力最大，在此条件密封环发生变形，将其变形情况放大30倍后，其变形如图4所示。

密封环发生微小变形，端面间形成楔形空间。由于密封环的变形，密封端面外径区域接触紧密，外径表面承受的载荷较大。在高速转动条件下，外径区域易产生较严重的磨损。

2.2 温度场仿真结果分析

选取启动过程中的4个速度值，分析该条件下的温度场，如图5所示。

由温度场分析结果可知，在整个运行过程中，接触区域出现高温，温度最高情况出现在瞬间启动时，接触区域最高温升在极短时间内约达到180℃，较高温升会加剧密封环的磨损[15]。

2.3 静环磨损形貌结果分析

取原始静环工件，按照从密封环外径到内径的顺序分析原始表面进行表面形貌。静环原始表面平整，材料均匀，各处表面形貌没有明显差异，有大量的石墨填充物，同时也存在轻微的加工划痕。其中一处放大300倍后的表面形貌图如图6(a)所示。

磨损后的静环，按照从密封环外径到内径的顺序进行表面形貌分析，其放大300倍后的表面形貌图如图6(b)，图6(c)及图6(d)所示。

与图6(a)相比，图6(b)中的静环外径区域表面存在材料的粘着、剥落及犁沟，平整端面被破坏，磨损严重；图6(c)中，磨损表面有明显的分界，靠近外径区域磨损严重，靠近内径区域磨损轻微，表面较为完整；环内径区域表面完整，磨损轻微，如图6（d）。接触副摩擦生热，铜石墨与密封介质发生化学反应，铜发生氧化，使磨损件表面存在一层氧化物。

3 磨损机理分析

3.1 静环材料转移机理分析

所研究的密封形式为接触型密封，具有高速、高载荷的特点，且静环所用的铜基石墨材料强度较低。在高PV值的作用下，在滑动接触区域表面温度急剧升高，影响材料的摩擦磨损性能，摩擦系数增大，因铜基石墨材料较软，表面材料被撕裂，发生材料的转移，部分铜基石墨材料附着在动环上。发生材料转移时，配副间出现磨屑，使动环出现轻微磨损。

最终，在动环和静环上都有明显的材料转移，使表面凸起，在静环上有不规则微坑和犁沟产生。

3.2 内外径磨损情况差异的机理分析

由接触应力和温度场分析可知：在介质压力、弹簧支撑力等因素的共同作用下，引起材料的轻微变形，形成如图4所示的楔形，外径区域的接触应力较大，外侧易产生磨损。在高速转动过程中，密封接触面摩擦产生大量热，影响静环的材料性能，加剧材料转移的发生。表面的磨损，必然引起密封环系统的振动，将加剧磨损情况。

综合静环磨损形貌分析、接触应力及温度场有限元分析，摩擦过程中，主要发生静环的材料转移。接近外径区域接触应力大，磨损严重，导致密封的泄漏；接近内径区域接触应力较小，磨损轻微；密封环变形和摩擦产生的大量热也会对密封端面的磨损性能产生影响。

4 结论

1）该机械密封磨损的主要形式是粘着磨损，静环端面外径区域磨损严重，内经区域磨损轻微。

2)密封环接触区域外径区域接触应力大，内径处接触应力小；密封环发生微变形，接触端面形成楔形，外径区域较内径区域接触紧密。综合导致外径区域磨损情况较内径区域严重。

3)启动过程中，高速摩擦产生较大温升，摩擦产生的大量热影响静环材料性能，加剧密封环端面的磨损。

4)综合表面形貌分析和ANSYS有限元分析及磨损机理分析，为所研究的液体火箭发动机涡轮泵机械封环的性能优化提供了思路。

[1]张国渊,袁小阳,赵伟刚,等.螺旋槽端面密封脱开转速的理论及试验[J].机械工程学报，2008,44(8)：55-60.

[2]白东安,段增斌,张翠儒.涡轮泵端面密封性能与漏气量影响研究[J].火箭推进,2010,36(1)：38-42. BAIDong'an,DUAN Zenbin,ZHANG Cuiru.Effects of leaking rate on turbopump end-face sealing performance [J].Journalof RocketPropulsion,2010,36(1)：38-42.

[3]张树强,王良,赵伟刚.液体火箭发动机涡轮泵用机械密封温度场及热载变形研究[J].火箭推进,2014,40(5)：92-98. ZHANG Shuqiang,WANG Liang,ZHAO Weigang. Research on temperature field and heat deformation of mechanical seal in liquid rocket engine turbopump[J]. Journalof rocketpropulsion,2014,40(5)：92-98.

[4]黄智勇,胡钟兵,李惠敏.大功率、高转速、高扬程涡轮泵振动分析与减振研究[J].火箭推进,2005,31(6)：1-6. HUANG Zhiyong,HU Zhongbing,LIHuim in.Analysisof vibration and vibration reduction for turbopumps w ith high-power,high-rotation speed and high-delivery-head [J].Journalof rocketpropulsion,2005,31(6)：1-6.

[5]WANG J,Jia Q,YUAN X,et al.Experimental study on friction and wearbehaviourofamorphouscarbon coatings for mechanical seals in cryogenic environment[J]. Applied surface science,2012,258(24)：9531-9535.

[6]WEIL,ChANGX Z,ZhANGPG.Fractalcharacterization of changing of surface topography of the end face of soft ring in running-in procedure formechanical seal[J].China surfaceengineering,2011,24(5)：78-82.

[7]张书贵,顾永泉.机械密封变形的研究[J].中国石油大学学报(自然科学版),1992(2)：48-53.

[8]刘伟,彭旭东,白少先,等.流体静压型机械密封的三维传热数学模型及端面温度分析 [J].摩擦学学报,2010, 30(1)：57-63.

[9]TAKAO S,HIROSHI H,TOSHIKATSU N.A study of the relationship between sealing performance and surface condition of mechanical face seals[J].Tribology transactions,1992,35(4)：659-666.

[10]E.Mayer.机械密封[M].化学工业出版社,1981.

[11]顾永泉.机械密封实用技术 [M].北京：机械工业出版社,2001．

[12]BATHEISK,Kirk RG.Analysis of high pressure liquid seal ring distortion and stability using finite element methods[J].Journalof tribology,1999,121(4)：921-926.

[13]LEBECK A O.A mixed friction hydrostaticmechanical face seal model with thermal rotation and wear[J]. Tribology transactions,1980,23(4)：375-387.

[14]顾伯勤,周剑锋,陈晔等.机械密封端面间液膜摩擦热的传热规律[J].中国科学(E辑：技术科学),2008,(1)：137-147.

[15]刘伟,彭旭东,白少先等.流体静压型机械密封的三维传热数学模型及端面温度分析[J].摩擦学学报,2010, (1)：57-63.

[16]GOLUBIEV A I.Studies on seals for rotating shafts of high-pressure pumps[J].Wear,1965,8(4)：270-288.

[17]YAN G P,Liu Z L,Zhu X M,etal.Numericalanalysisof the thermal-field of ship stern-shaft mechanical sealed faces under the variationalworking conditions[J].Journal of shipmechanics,2008,12(3)：483-489.

（编辑：马 杰）

Research on wear mechanism of mechanical seal for turbopump in liquid rocket engine

ZHAOWeigang1,ZHANG Pengpeng2,REN Shanshan2,DONGGuangneng2
（1.Schoolof Astronautics,Northwestern PolytechnicalUniversity,Xi’an 710072,China; 2.MOEKey Laboratory forModern Design and Rotor-bearing System, Xi’an Jiaotong University,Xi’an 710049,China）

The original surface morphology and the worn surface morphology of static ring end-face of mechanical seal was analyzed by metallographic m icroscope.The wear form and the differences between inner diameter and outer diameter were obtained.The contact stress and temperature field distribution on the sealend-faceareanalyzed bymeansofsimulation.The influences of rotating speed,medium pressure and other conditions on end-face contact stress and temperature field arealso analyzed bymeansof ANSYSsoftware platform.Theeffectofend-face contactstressand temperature variation on the propertiesofwearwas verified.Thewearmechanism ofmechanical seal wasexpounded.

surface topography;simulation analysis;mechanicalseal;wearmechanism

V433-34

A

1672-9374(2017)03-0010-05

2016-10-25；

2017-01-06

民用航天产品高可靠长寿命专项（〔2014〕618号）

赵伟刚（1977—），男，高级工程师，研究领域为机械密封