龚金科 卜艳平 钟超 胡辽平 陈少林
摘 要:基于耦合传热理论,建立了废气旁通阀式的涡轮箱模型,并将涡轮箱、隔热罩、轴承体作为一个装配体进行了传热分析,得到了轴承体的温度场和热应力分布.结果表明:温度从涡轮箱至轴承体逐渐降低,呈现明显的温度梯度,涡轮箱的最高温度比初始废气温度低71.5 ℃左右;由于隔热罩的隔热作用和冷却水、机油的双重冷却作用,轴承体的温度较低;由于轴承体内部和表面区域温差较大,其热应力较高.与实验对比,轴承体温度的仿真值与实验值误差最大为7.2%,说明该仿真方法具有较高的精度,能为增压器的设计和优化提供一定的理论依据.
关键词:涡轮箱;隔热罩;轴承体;耦合传热;热应力
中图分类号:TK411.8 文献标志码:A
Abstract:Based on coupled heat transfer theory, a model of turbine box with waste gate was built, and turbine box, heat shield and bearing body was tied as an assembly to be heat transfer analyzed. Temperature and thermal stress analysis of bearing body were then obtained. The results show that temperature is gradually reduced from turbine box to bearing body, which presents an obvious temperature gradient. And the temperature of turbine box is about 71.5℃ lower than original gas temperature. Because of good insulation of heat shield and dual cooling of cooled water and oil, the temperature of bearing body is low. However, the thermal stress of bearing body is high due to its temperature difference between inside and outside surface. Compared with the test results, the maximum deviation between the simulated and experimental value is 7.2%, which verifies that this simulation method has high precision and can provide a theoretical basis for the design and optimization of turbocharger.
Key words:turbine box; heat shield; bearing; coupled heat transfer;thermal stress
增压器在高温、高速的环境下工作,使得涡轮箱、涡轮叶片、轴承体、浮动轴承等关键零部件承受非常高的温度及热应力.日益严格的经济性、环保性要求使得增压器的强化程度不断提高,其热负荷、冷却散热问题严重影响了增压器结构强度和可靠性.且轴承和密封环受轴承体热负荷的影响,其使用寿命明显降低,成为增压器可靠性中的一个薄弱环节[1-3].
增压器冷却散热系统如果不能有效降低轴承体的温度,在长时间工作中就会使轴承体零部件产生热疲劳,同时还会进一步降低机油润滑效果和承载能力.长期以来对轴承体进行流固耦合传热分析时,采用的是轴承体涡轮端和压氣机端的温度边界条件[4-6].但在产品开发前期,以上边界条件未知,需采用涡轮箱的进气边界和轴承体压气机端的温度边界对涡轮箱、隔热罩、轴承体进行整体传热分析,以此分析轴承体的热负荷情况.
本文根据工程实际的需要,建立了废气旁通阀式的涡轮箱模型,考虑隔热罩隔热作用,分析涡轮箱、隔热罩和轴承体装配体的传热情况.并根据耦合传热理论,对高温废气和涡轮箱、油腔和轴承体、冷却水腔和轴承体分别进行耦合传热分析,使数值仿真更符合轴承体的实际工作状态,获得了精确的轴承体温度场及热应力分布情况,并进行了实验验证.
1 基本理论及计算流程
因为固体壁面和流体之间存在相互制约作用,两者的对流换热是非常复杂的过程,热边界条件不能预先规定[7].此时边界上的温度、换热系数不是已知条件,应看成计算结果的一部分[8].而耦合传热方法可以很好地解决这类由热量交换过程动态决定的热边界条件问题.
本文采用CFD软件计算得到各流体域的壁面温度T和换热系数h,并将其映射到流固交界面上,成为固体域的热交换边界条件.再将计算得到的固体域温度场映射到流体域,成为流体域计算的壁面温度条件,按此步骤重复迭代耦合.计算流程图如图1所示.
2 涡轮增压器仿真模型
2.1 涡轮增压器仿真网格模型
在数值仿真计算前,为避免计算出错、收敛更快,简化处理涡轮箱、轴承体外表面和细小结构,并加密涡轮箱舌口和轴承体水套等倒角处网格.通过几何处理后,对涡轮箱、隔热罩、轴承体进行体网格划分.各部件的网格信息如表1所示.
对废气旁通阀式的涡轮箱流体域划分体网格,首先将涡轮箱内腔壁面提取出来,并对进出口作封面处理,使得面网格转换成实体,再将涡轮旋转域分离出来.在CFD软件中建立交界面用于不同区域之间的能量及物质传递,并在涡轮旋转区域建立MRF模型以实现动网格的高速旋转,以更精确地模拟高温废气的流动.废气旁通阀式涡轮箱流体域网格模型如图2所示.
为保证流体计算精度,在轴承体中提取冷却水腔和机油腔的表面网格并导入相关软件中进行局部加密,生成流体3D网格.对边界层网格进行细化以精确地模拟壁面附近的流体流动情况.轴承体冷却水腔和机油腔的网格模型分别如图3(a)(b)所示.
2.2 物理模型材料
涡轮箱在高温环境下工作,采用抗热性好的片状石墨奥氏体铸铁XNi35Si5Cr2;隔热罩采用奥氏体型耐热钢06Cr19Ni9;轴承体则采用灰铸铁HT250.在仿真模拟中,对于这些材料通常需要定义其弹性模量、泊松比、热传导系数、热膨胀系数等,以获得准确的仿真结果.表2列出了各部件的材料属性.
对于高温废气和冷却水则需要定义其比热容、热导率、普朗特数等,且这些参数是温度的单值函数,根据传热学[9]设定.机油的牌号为5W40,其密度设为846.4 kg/m,比热容为2.265 J/(kg·K) , 热导率为0.141 6 W/(m·K),Pr为190,而粘度属性随温度变化很大,且影响机油的流动,机油粘度属性根据表3给定.
2.3 计算条件
在对轴承体进行仿真研究时,一般以测量得到的轴承体涡轮端和压气机端的温度边界作为计算边界条件.但在工程实际中,产品开发前期,上述边界条件往往未知,此时,应以废气作为热源,通过废气与固体壁面换热,并经过涡轮端及隔热罩的热传导,作为轴承体计算的涡轮端边界条件.热源改变时的轴承体计算模型如图4所示.
2.4 边界条件
在仿真过程中,认为流体为三维粘性不可压缩流体,采用kε湍流流动模型.计算涡轮箱流场时,涡前废气入口采用质量流量边界条件,入口质量为0.075 kg/s,温度为950 ℃,涡轮转子轴转速为250 000 r/min;出口则采用压力出口边界条件,出口压力为0.134 MPa,温度为25 ℃.
對于轴承体内的流场,冷却水和机油均采用压力出口边界条件,保持进出口压比为1.5.轴承体内流体的边界条件如表4所示.
3 计算结果分析
3.1 涡轮箱、隔热罩和轴承体装配体温度计算结果
在相关软件中预先将涡轮箱、隔热罩、轴承体进行装配,定义约束等条件,并模拟冷却水和机油在轴承体内的冷却和流动,最后将涡轮箱和轴承体内流场计算得到的换热系数和温度导入装配体中进行温度传热分析.图5和图6分别为装配体在机油腔剖面和水腔剖面的温度分布.
由图5和图6可知:温度从涡轮箱至轴承体逐渐降低,形成明显的温度梯度.涡流道、旁通阀流道以及出口端区域温度最高,这是由于涡轮箱的热量来源于发动机排气,在高温废气通过流道壁面时,产生耦合换热;但这些区域的温度比废气的初始温度低71.5 ℃左右,这是由于经过多次流固耦合换热计算,机油和冷却水降低了涡轮箱整体的温度.
为分析隔热罩的隔热作用,图7所示为隔热罩正反面的温度分布情况,隔热罩与涡轮箱接触的一侧温度明显高于与轴承体接触的一侧,可见在降低从涡轮箱传递至轴承体的热量上隔热罩起到了很重要的作用.
3.2 轴承体温度场分析
轴承体温度分布如图8所示.由图可知:从涡轮端至压气机端轴承体的温度逐渐降低,最高温度为349.58 ℃,主要分布在与隔热罩接触的部位.这是由于涡轮箱在传递热量的过程中,高温首先将热量传递至隔热罩,经过隔热罩的隔热作用,剩余的主要热量通过隔热罩的热传导传递至轴承体,故与隔热罩接触的部位轴承体温度最高.最低温度为135 ℃,主要分布在与压气机接触的边界上.
图9为轴承体平行于水腔的切片温度场,在轴承体冷却水腔的顶部,即与隔热罩连接处,其壁面最高温度达到了349.58 ℃.涡轮箱的热量沿着轴承座方向传递给轴承体,轴承体顶部热量未得到及时冷却,热负荷很高.但是经过冷却水的充分冷却,在靠近压气机一侧,温度已经下降至135 ℃.
图10为轴承体平行于机油腔的切片温度场,机油对轴承体不仅有润滑和清洗的作用,还能较好地对轴承体进行冷却.和图9类似,温度从机油腔顶部至底部依次降低,呈现明显的温度梯度.靠近压气机端一侧的轴承座附近,由于涡轮转子轴高速运转,机油摩擦生热导致部分区域温度偏高.
3.3 轴承体热应力分析
轴承体在工作过程中温度分布不均匀,且其两端被固定,必然会产生热应力,需对轴承体进行热应力分析.为了保证与温度场分析模型的一致,应将传热单元转化为结构分析单元[11].同时将流固耦合计算得到的轴承体温度场作为应力计算的预定义场,并限制轴承体两端的自由度.
图11和图12分别为轴承体平行于轴承座方向并分别通过冷却水腔和机油腔的剖面热应力图,从总体上看,轴承体油腔和水腔内表面的热应力较高,这主要是因为这些表面区域和流体进行了换热,而轴承体内部存在热阻,热量传递相对滞后,与表面形成了较大的温差,从而导致了较大的热应力.如果热应力过大,可以控制冷却液和轴承体冷却腔壁面的温差,在满足使用条件的前提下,适当提高进水温度和限制进水流量,来减少这部分热应力.
最大的热应力分布在涡轮端的轴承座内,这是因为根据轴承体的温度分布趋势,涡轮端的温度较高,且机油与高速旋转的转子轴摩擦也会产生热量,这些热量使得涡轮端的轴承座内热应力较大.对于该部位,可以适当减少轴承座内圈半径,以增加轴承座的壁厚.
4 实验验证
在对轴承体进行模拟计算时,传热边界条件的确定与仿真结果的精确性密切相关[12].为验证仿真结果是否准确,在轴承座上取9点测量其温度,测点位置标示在图13中.其中需要验证测点1、5模拟仿真的数据与实际边界数值的误差;同时,轴承和密封环的寿命直接受轴承体上涡轮端轴承座和密封档台部位的温度影响[13],因此有必要选取测点2和测点9对这两处的温度进行实际测量.
在增压器稳态运行工况下试验,使用燃烧室内加热后排出的气体模拟发动机废气驱动增压器,控制废气温度为950 ℃,流量为0.075 kg/s;冷却水温度为85 ℃,压力为0.25 MPa;机油温度为105 ℃,压力为0.45 MPa;涡轮转子轴转速为250 000 r/min.当增压器运行20 min并稳定后,测量所需要的温度值.
各测点温度仿真值与实验值的对比如表5所列,仿真计算结果与实测结果误差最大为7.2%.这是由于本文研究轴承体所需的边界条件是从涡轮箱的计算结果得到,累计的误差会比直接测得的轴承体边界条件大.作为涡轮增压器的前期设计,耦合传热方法可以减少产品开发的成本,缩短开发周期.
5 结 论
1)本文建立了废气旁通阀式的涡轮箱模型,并对涡轮箱、隔热罩、轴承体整体进行了耦合传热分析.结果表明,温度从涡轮箱到轴承体逐渐降低,最高温度分布在涡流道、旁通阀流道以及出口端区域;隔热罩与涡轮箱接触的一侧温度明显高于与轴承体接触的一侧,有效地减少了热量从涡轮箱传递至轴承体.
2)轴承体最高温度出现在与隔热罩接触的部位,经过冷却水和机油的双重冷却后,温度从349.58 ℃下降至135 ℃,充分体现了轴承体良好的冷却性能.但轴承体油腔和水腔内表面的热应力较高,可以控制冷却液和轴承体冷却腔壁面的温差,如提高进水温度和限制进水流量,来减少这部分热应力.对于热应力最高的轴承座,可增加其壁厚,以减少热应力.
参考文献
[1] LIN H C, CHANG Y T, LONG T G,et al. Oil coking prevention using electric water pump for turbocharge sparkignition engines[J]. Mathematical Problems in Engineering, 2014(6):1-8.
[2] 马国山,毕金光,王银凤,等.汽油机涡轮增压器的特点初谈[J].汽车实用技术,2016(8):85-87.
MA G S, BI J G, WANG Y F, et al. A preliminary discussion on the characteristics of turbocharger for gasoline engine[J]. Automobile Applied Technology, 2016(8):85-87.(In Chinese)
[3] BURKE R D, COPELAND C D, DUDA T,et al. Lumped capacitance and threedimensional computational fluid dynamics conjugate heat transfer modeling of an automotive turbocharger [J]. Journal of Engineering for Gas Turbines and Power, 2016, 138(9): 92602-92601.
[4] 龔金科,章滔,胡辽平,等.涡轮增压器水冷轴承体冷却性能仿真研究[J].汽车工程,2014,36(3):383-388.
GONG J K, ZHANG T, HU L P, et al. A research of the simulation on the cooling performance of watercooled bearing in turbocharger[J]. Automotive Engineering, 2014,36(3):383-388.(In Chinese)
[5] HEUER T, ENGELS B, WOLLSCHEID P. Thermomechanical analysis of a turbocharger based on conjugate heat transfer[C]//Proceedings of ASME Turbo Expo 2005. RenoTahoe, Nevada, USA, 2005:829-36.
[6] 龚金科,田应华,贾国海,等.涡轮增压器轴承体耦合传热的数值仿真[J].湖南大学学报(自然科学版),2015,42(4):7-12.
GONG J K, TIAN Y H, JIA G H, et al. Numerical simulation on coupled heat transfer of turbocharger bearing [J]. Journal of Hunan University(Natural Sciences),2015, 42(4):7-12. (In Chinese)
[7] 张艳岗,毛虎平,苏铁熊,等.基于整场离散的气缸盖热流固多物理场耦合状态数值模拟[J].北京理工大学学报,2016,36(6):563-568.
ZHANG Y G, MAO H P, SU T X. Numerical simulation based overall discrete method for thermofluidsolid coupling of cylinder head[J]. Transactions of Beijing Institute of Technology, 2016,36 (6):563-568. (In Chinese)
[8] 赵亚茹.发动机活塞组缸套冷却水耦合传热研究[D]. 太原:中北大学, 2016:43-45.
ZHAO Y R.The coupling heat transfer research on piston grouplinercoolant of engine[D]. Taiyuan:North University of China, 2016: 43-45.(In Chinese)
[9] 贾力,方肇洪,钱兴华.高等传热学[M]. 北京: 高等教育出版社, 2008:117-125.
JIA L, FANG Z H, QIAN X H. Higher heat transfer[M]. Beijing: Higher Education Press, 2008:117-125. (In Chinese)
[10]LAUNDER B E, SPALDING D B.Lectures in mathematical models of turbulence[M].London: Academic Press,1972.
[11]李发宗,童水光,袁铭鸿.热流固耦合下增压器涡壳多目标优化[J]. 中国机械工程,2016,27 (6) :754-766.
LI F Z, TONG S G, YUAN M H. Multiobjective optimization for turbocharger housing considering thermofluidsolid interaction[J]. China Mechanical Engineering, 2016,27 (6) :754-766. (In Chinese)
[12]吴新涛,何洪,徐思友,等.涡轮增压器轴承体传热计算方法研究[J].车用发动机,2010,186(1):11-14.
WU X T, HE H, XU S Y, et al. Heat transfer calculation method of turbocharger bearing[J].Vehicle Engine, 2010,186(1) : 11-14. (In Chinese)
[13]王本亮,刘石源,唐维新,等.汽油机涡轮增压器轴承体回热试验研究[J].车用发动机,2013, 209(6):40-43.
WANG B L, LIU S Y, TANG W X, et al. Heat soakback testing of turbocharger bearing housing for gasoline engine[J]. Vehicle Engine, 2013, 209(6): 40-43. (In Chinese)