基于热传导角的机械密封环温度分布计算方法

本资料包含pdf文件1个，下载需要1积分

Calculation Method of Mechanical Seal Temperature DistributionBased on Heat Conduction AngleSong Pengyun Gao Jie Ma Fangbo Zhang Chun Mao Wenyuan(1.Faculty of Chemical Engineering，Kunming University of Science and Technology，Kunming Yunnan 650500，China；2.Kunming Tighscis Sealing Technology Co.，Ltd，Kunming Yunnan 650500，China)Abstract：The temperature distribution on the end faces of the mechanical seal is an important factor which afects thesealing performance.The approximate analytical calculation of the end face temperature can analyze the effect of tempera-ture on sealing perform ance easily.In the inner parts of the rotor and static rings，heat transfers in the form of heat conduc-tion，while in the contacting area of the tings and the medium，heat transfers in the form of convection.Considering theclearly defined heat conducting angle，and heat distribution in the rotor and static rings，a form ula calculating the innertemperature distribution of the rotor and static tings was derived based On the heat conducting theory.The calculating re-suhs were compared with the numerical simulation results.The resuhs show that as the end faces distance and the radiusincrease，the average errors of the temperature between the approximate analytical calculation and the numerical simulationresults will decrease。

Keywords：mechanical seal；face temperature；numerical calculation机械密封端面温度的确定与控制是机械密封设计和使用的-个关键技术。密封端面温度对机械密封性能有着明显的影响 ，端面温度分布的研究-直是机械密封研究领域-个重要内容。目前确定端面温度计算方法主要有以下几种 ：端面平均温度简算法，数值计算法和近似解析法等。其中解析法能从本质上揭示规律，其计算结果具有确定性和唯-性，因此确定机械密封端面温度的解析计算方法是机械密封研究-直致力追求的目标。1980年李克永等 推导了确定基金项目：昆明理工大学资助项目。

收稿 日期 ：2012-07-23作者简介：宋鹏云 (1963-)，男，博士，教授，主要研究方向为流体密封技术.E-mail：songpengyun###hotmail.com。

简单矩形断面密封环温度分布的解析计算方法。1986年李红等人 针对实际使用的密封环形状，推导了L型密封环的端面温度解析解。

1992年 Pascovici和 Etsion提出热传导角 的概念 ，1994年Etsion和Groper给出了热传导角 的明确定义 ，推导了全液体润滑状况下的端面温度分布计算公式，对比了解析计算与数值模拟计算结果。

1995年朱孝平等 ，2004年胡小云 推导了双端面机械密封混合摩擦状态下的端面温度分布计算公式，其中对于热传导角的计算，胡小云根据实际情况进行了计算。这几种解析法不足在于将静环表面 (包括与介质接触的表面和与空气接触的表面)作为绝热边界处理，忽略了热量在静环上的传递。2012年高2013年第 1期 宋鹏云等：基于热传导角的机械密封环温度分布计算方法 7杰等人 考虑热传导角、热量在动静环上的分配，推导了混合摩擦状况及全液体润滑状况的密封端面温度分布，对比分析了不同近似解析法，并分析了影响端面温度分布的因素及其规律。

但是前面几种方法仅仅只是分析了动静环接触端面的温度分布，对于动静环内部温度的分布 ，并没有给出计算表达式。因此，本文作者根据热传导理论，假定热流沿热传导角方向由密封间隙传至动静环的圆柱表面，并考虑热量在动静环上的分配，对密封环内部温度分布进行理论分析，并将理论分析结果与数值模拟结果进行对比。

1 理论计算1.1 计算模型对于动环或静环来说，假定热流是沿着固定角方向由密封间隙传至动环或静环圆柱表面，然后再与密封腔液体进行热量交换。因此热量在动静环内部的传递为热传导，在动静环表面与密封腔液体之间的热传递为对流换热。理论分析计算模型如图1所示，密封端面摩擦产生的热量沿固定的热传导角 传递到密封环的圆周表面，再与被密封液或冲洗液进行热量交换 (图中虚线为热流传递方向)。同时考虑到实际情况下，热量并不能完全由动环进行传递，静环也会传递-部分热量，即要考虑热量在动静环上传递时的分配 比图1 理论分析模型Fig 1 Theory analysis model图2示出了考虑热传导角的计算模型，其中 表示轴向方向，r为半径方向。表面1为动环或静环接触端面，表面2为动环或静环与介质接触表面，表面3为动环或静环内侧表面。

计算数据如下：外径 81.5 mm；内径 ri76mm；端面线速度 10 m/s；A179.6 W/(m·oC)；A 275.4 W/(m·oC)；k0.：599 W/(m·oC)； 20℃；弹簧比压为0.2 MPa；介质压力为1.5 MPa。

l I--J图2 考虑热传导角的计算模型Fig 2 The calculation model considered heat conduction angle1.2 传热效率分配比在实际模型中，不仅动环会传递热量，静环也会传递-部分热量，因此，在计算时要考虑热量在动静环上的分配。动静环热量分配系数的计算由文献[9]可得，动静环端面热量分配系数计算如下：tan 21- (1)r。HI，2kE ： (2) -A1cos 1A2cos 2 /E，：--坐 (3)- 、 J ，A1c0s 1A2cos 2 、式中：E 为动环的热量分配系数；E：为静环的热量分配系数。

从式 (1) ～ (3)能够看出，热量分配比与动静环材料的导热系数A、热传导角 动静环的对流换热系数日及几何模型等参数相关。由式 (1)可得动环 的热传导角为 54.06。，静环的热传导角为57.45。，同时考虑动静环材料的导热系数分别为79.6和75.4 W/(m·oC)，代人式 (2)， (3)中，计算可得热量分配比为：动：环E 0.535 2，静环E20.464 8。

1.3 考虑热传导角的端面温度解析法热流沿热传导角 从壁面1传至壁面2时，如图2中虚线所示，运用热传导理论中傅里叶定理，则不同径向位置r，和不同轴向位置 的密封环温度为： (f ( - )(4)8 润滑与密封 第 38卷壁面2与封液之间的热量交换为对流换热。则由传热学知识可得：～ A l：日( -rf) (5)dn l得 苦 (6)联立式 (4)，(6)可得) 吾 ( - simp) (7)其中，qfPgrto式中： 是径向位置 r，及轴向位置 的函数，对应不同的r，、 ，能够得到不同位置的温度分布。

由式 (7)可得不同轴向位置、不同径向位置的温度计算公式。对于动静环来说，密封环内部的温度分布可由分别由下两式确定。

) 十 ( --simp1)(8)) 噎 ( -志)(9)2 计算结果对比分析(10)按照上述方法，对动静环进行温度计算，并将计算结果与数值模拟结果进行对比。计算结果以动环为例进行分析。不同半径、不同轴向位置的动环内部温度分布结果如图3所示∩以看出，随着半径增大，tmm(a)距端面距离为0在-定的轴向位置时，温度是下降的，这是因为随着半径增大，对流换热效果逐渐增大，换热效果好，温度会逐渐降低。同时在径向位置固定，随着轴向位置的增加，温度也是增加的，这是因为随着轴向位置增大，远离密封端面的摩擦热源，故而温度逐渐降低。

908070呈605040300.Ol51图3 动环内部温度分布Fig 3 The temperature distribution of rotor ring2.1 动环径向温度计算结果对比图4示出了距端面距离分别为0、1.2、2 mm处的径向温度分布对比，结果表明径向温度分布都是随半径增加而逐渐减小，数值模拟结果表明温度分布减小的趋势近似为抛物线，而解析解减小的趋势为直线。3种端面距离处径向温度的解析解与数值模拟解平均误差为 19%、12.5%、10.3%。即离密封端面越远，径向温度的解析解与数值模拟解的平均误差越校fmm(b)距端面距离为1.2 mmrlmm(c)距端面距离为 2mm图4 动环径向温度解析解和模拟解对比Fig 4 Comparison of radial temperature of the rotor ring obtained by analytical calculation and simulationanalysis(a)0 from the end face；(b)1.2 mm from the end face；(C)2 mm from the end face2.2 动环轴向温度分布对比图5示出了不同半径处的轴向温度解析解与数值模拟结果的对比。不同半径处，随着轴向位置的增加，距离密封端面越远，端面温度逐渐降低，解析解为直线，降幅较小，数值模拟解近似为抛物线，降幅较大。图 (a) ～ (C)的解析解与数值模拟结果的2013年第 1期 宋鹏云等：基于热传导角的机械密封环温度分布计算方法 9对比表明，随着半径的增大，不同半径处的轴向温度平均误差是逐渐减小的。在图 (a) ～ (c)表示80706O50403020的不同半径处，轴 向温度的平均误差分别 为：20.5% 、13.5% 、8.2% 。

轴 向位置x/mm 轴向位置x/mm(a)r75.8 mm (b)r78.1 mm轴向位置x/mm(c)r80 mm图 5 动环轴向温度解析解和模拟解对比Fig 5 Comparison of axial temperature of the rotor ring obtained by analytical calculationand simulation analysis at r75.8 mm (a)；r78.1 mm (b)；r80 mm (c)3 结论(1)考虑热量在动静环上的分配、明确定义热传导角，推导了密封环内部温度分布计算公式，得到密封环内不同轴向位置及不同径向位置的温度分布，并将计算结果与数值模拟结果进行了对比分析。径向温度及轴向温度分布趋势如下：随着半径的增加，或随着距端面轴向距离的增加，温度均呈下降趋势↑似解析解表明，温度呈线性下降趋势，而数值模拟结果呈似抛物线下降趋势。不过，近似解析结果与数值模拟结果的平均误差均随端面距离的增加或随半径的增加呈下降趋势。

(2)近似解析计算对模型进行了部分简化，在密封环上的热量分配采用公式进行计算，但数值模拟计算时的模型与实际模型相同，采用整体建模方式，密封环上的热量分配比由系统自动分配。因此，两者计算的结果会出现-定的误差。有理由相信，数值模拟的结果更加接近实际情况，但解析计算方法更简洁和更实用。