活塞热疲劳分析

柴油机设计与制造Design Manufacture of Diesel Engine2013年第1期第19卷(总第142期)do:10.3969/jssn.1671-06142013.01003活塞热疲劳分析谢琰',席明智2,刘晓丽3(1.长安汽车动力研究院,重庆4001;2内蒙古工业大学能源与动力工程学院,呼和浩特014010;3.渤海船舶职业技术学院,葫芦岛市1250000摘要用Pro/E建立活塞几何模型,在 ANSYS单元库里选取热结构耦合单元,对模型网格进行优化,并对活塞温度场进行标定,然后进行热机耦合分析计算,得到活塞温度场、热应力场和变形。计算结果表明,在低频热疲劳下,活塞循环次数最少约是112000次大,这为活塞的结构改进和优化提供了依据。关键词:活塞热应力热变形热疲劳analysis of Piston Thermal FatigueXie Yan, Xi Mingzhi?, Liu Xiaoli(1.CHANGAN Automotive Motor Engine R&D Institute, Chongqing 400021, China2. School of Energy and Power Engineering, Inner Mongolia University of Technology, Hohhot 010051,China:3. Bohai Shipbuilding Vocational College. Huludao 125000, China)Abstract: A 3-D model of piston was built with Pro/E and thermal structure couple unit was obtainedKey words: piston, thermal stress, thermal deformation, thermal fatigue 65Bom ANSYS unit base. Then optimization of model mess and calibration of piston thermal field were madeThermal couple analysis was carried out and piston thermal field, thermal stress and thermal deformationwere obtained. The analysis result shows that the number of piston life circulation is 1 12000 at least underlow frequent fatigue, which provides a base for the improvement and optimization ofn structure1引言负荷状况,首先严格按照图纸建立活塞几何模型,某些情况下,在发动机的受热零部件中,活塞对所建活塞模型温度场与实测特征点进行了温度值的热应力可能比燃气爆发压力造成的机械应力要大的对标,最终通过 ANSYS软件计算得到活塞的热好几倍。随着活塞温度的提高,对于铝合金活塞,应力和热变形,分析低频热疲劳下活塞的寿命。当温度达到300℃时;材料的抗拉强度会下降2活塞材料参数和模型建立22%,当超过350℃时,则抗拉强度将下降到原来ZHl105W型柴油机缩口四角ω燃烧室活塞采的一半。如果活塞的热应力和热变形过大,就会造用硅铝合金材料ZL09G,其常温下的弹性模量成发动机不能正常运转,所以活塞的热负荷仍然是E=7100MPa,泊松比p=031,导热系数A=124W一个不容忽视的问题叮。由于热应力是活塞总应力(m2K),比热c=902J/(kgK,密度p=2700的主要来源,热膨胀变形在活塞总变形中占绝对主kgm3,20-300℃时的材料线形膨胀系数B=20.96导地位,而机械负荷的作用仅使活塞边缘向内弯10°/℃,材料的抗拉强度an=2682MP,抗压强度曲、抵消边缘向外的热膨胀变形,贡献很小四。因a=260.7MP。此,针对活塞进行热疲劳分析。活塞模型的建立,首先采用Pro/E软件建立活为了取得精确的缩口四角ω燃烧室活塞的热塞的三维几中国煤仁活塞模型通过来稿日期:2012-09-26CNMHG作者简介:谢琰(1979-),男,硕士,主要研究方向为动力机械结构CAE分析及设计方甚。PoE和 ANSYS软件的接口导入到 ANSYS中,从稳态计算。计算表明,这样的约束不使活塞产生刚而得到活塞在 ANSYS中的几何模型,然后进行活体位移,也没有引人附加载荷,是合理的塞模型的后续处理工作,图1是活塞几何模型。41活塞热应力与热变形的计算结果分析在温度场的基础上进行的有限元热应力分析结果如图3和图4所示。八N热应力/P0(AVG)e0aNE=。684E+08300图1活塞三维实体几何模型84E+08图3标定工况下活塞 von mises热应力3活塞温度场由于活塞在标定工况下热负荷最为严重,因D-D AN此,本研究选择在标定工况下对活塞热应力和热变热应力/PaN=2477形进行分析。活塞的热应力与热变形的有限元分析6648+07实质上是活塞温度场和活塞结构的一种热结构耦合266E+0分析,它是在活塞温度场分析的基础上进行的。因此活塞温度场对低频热疲劳有着决定性的影响。该活塞的温度场分析,最高温度为311℃,分布在活塞燃烧室喉口,最低温度为120℃,分布在活塞裙图4热负荷下活塞截面 von mises热应力部下端。温度从上到下呈下降趋势,活塞温度分布图3和图4分别是活塞在标定工况下 von mises热趋势合理,如图2所示。应力整体图和截面图,从图中可以看出:(1)标定工况下活塞总体热应力不高,最高热品度/℃应力为596MPa,出现在排气一侧的回油孔顶部:1:52销座外侧销孔正上方第3环岸处热应力也较大,达到42.1MPa。主要原因是,该处有明显的尖角和棱183.19角,使得热流传递过程中热阻增大,出现热应力集中。0.672(2)活塞内腔顶部热应力较高,计算结果显示图2活塞温度场活塞内腔顶部最大热应力为439MPa。因此造成内4活塞应力场腔顶部中心温度高、温差大、热应力集中。本研究采用有限元间接法进行热应力分析,因(3)燃烧室进气侧旁部分底圈出现热应力集此热应力分析所用的活塞模型必须是温度场分析用中,达到462MPa,这是因为低温进气与高温燃气的模型。在进行热应力计算时,先进行温度场计算,交接碰撞而产生;燃烧室周面排气口侧出现应力集然后转换温度单元soid87到结构单元soid187;活中,达到386MPa;活塞其他部位,热应力不高,塞载荷直接从温度场以体载荷的方式读入,并对活基本都在30MPa以下;活塞销座和裙部的应力较塞进行约束,然后开始求解计算。小,基本都在18MPa以下。活塞的约束:将活塞一边销座中心上方内侧点图5是热负荷下的第一主应力场。由图可以看的y,z两个方向约束,将另一边销座同一位置点出,在热负荷下,活塞最大拉应力为689MPa,出的y,z两个方向约束,将活塞内腔上面中心点的现在刮油槽的V凵中国煤化工布在活塞头x,z两个方向约束。其中,x轴与销座孔轴线平部边缘和环槽CNMHG上缘也比较行,y轴是活塞中心轴线。该约束对标定工况进行大,其余地方都不大。最大压主应力出现在油环槽的回油孔处,达到10.3MPa。其中燃烧室底圈和活产生局部的残余变形,反复循环产生的热疲劳最终塞销座内侧上缘最大压主应力为10MPa,其余地将导致材料的破坏A。柴油机在起动-运行-停方应力不高。车的过程中造成的损伤最为严重。强度分析可以归结为预测热疲劳寿命68。第一主应力/Pa发动机的热负荷基本分为稳定热负荷、低频热89E+0日负荷和高频热负荷三类。稳定热负荷对应于发动机14:1:0在稳定工况运行时各受热件除受热表层以外绝大部337+08分结构所处的温度状况,即本文前面分析过的稳态601E+08温度场属于稳定热负荷的范畴;低频热负荷对应于669E+D8发动机在反复变换工况运转的过程中各受热件内部图5热负荷下活塞整体的第一主应力场温度的反复变化,以及由于各受热件内部的温度变4.2活塞热变形分析化滞后所造成的短时间改变的温度分布状况;高频图6是活塞在温度载荷下放大50倍的热变形热负荷对应于发动机运转过程中因缸内燃气温度周图。从图中可看出:活塞最大热变形量出现在活塞期变化所造成的受热件表层温度循环波动。从内燃顶面边缘排气口侧,达到0328mm。主要原因是机的可靠性和耐久性出发,稳定热负荷是设计者主排气温度比较高,热辐射能力强,气流速度较髙,要考虑的问题;高频热负荷可能引起附加的高频循对流换热加剧,致使这部分温度很高,变形量大。环热应力,这种热应力只在活塞表层,一般情况其活塞顶面边缘以及整个活塞头部的变形量都比较数值较小,可以不作专门考虑;低频热负荷加大了大,燃烧室底圈和凸台的变形量不大。活塞第1环活塞稳定热负荷造成的高温和温差,因而引起了附槽最大变形量达到025mm,主要是沿着活塞半径加增大的热变形和热应力,而且活塞热应力的反复方向膨胀;活塞第1环槽轴向变形量最大为0.07变化在持续一段比较长的时间后会导致活塞材料的mm,没有超出活塞环设计侧隙。活塞中间裙部和热疲劳损坏。柴油机在起动-运行-停车的过程销座变形量较小,在0.18mm以下;内腔顶部中心中负荷通常并不高,但经常作变工况运行。因此,变形量最小,在0.087mm以下。活塞裙部下端变设计人员在考虑稳定热负荷的同时,还应关注低频形量较大,达到0.23mm,是由于活塞结构和膨胀热负荷是十分必要的4s,国外也将低频热疲劳作为双重影响的结果。活塞整个变形呈两头大中间小的内燃机性能指标之-。趋势。研究热疲劳强度时要考虑一个非常重要的关系即一个循环的塑性应变与达到破坏重复次数的关热变形/mUBUM(AV系。目前针对低频热疲劳的寿命估算还没有一个统x=.328E-0的方法,但在工程和学术界中一般采用曼森和科62E-04芬提出的方法:即高温疲劳和蠕变交互作用的过程33E=03是消耗材料塑性的过程,当材料的塑性耗竭时就发241E-0270E=03生破坏叫。28E图6热负荷下放大50倍的活塞热变形C5活塞的低频热疲劳Er IC(1)热疲劳是由高温燃气周期性变化的温度作用下产生的。热疲劳源于材料内部为抵消物体热膨胀和式中收缩之差而产生的循环热应变,而且材料的延性与C—标志材料塑性大小的量;热应力强度密切相关。由于材料的延性存在,当热Δsp—循环塑性应变范围(全振幅);应力超过了材料本身的屈服点,即使尖峰应力值超N—断裂循环数过屈服点好几倍,在局部区域产生的塑性变形也不中国煤化工会立刻破坏材料,而在周边环境的影响下仍能压回CNMHG或拉回到原状,但当热应力超过屈服点太大时就要q—静拉伸断裂颈缩率。对于活塞低频热疲劳的分析,本文取活塞材料回油孔顶部;最大热变形为0328mm,出现在活ZL0G的材料颈缩率φ=50%,代入公式(1),得塞顶面边缘排气口侧。活塞的回油孔、活塞内腔顶到标志材料塑性大小的量为C=0.347。则计算低频部中心、销座外侧销孔正上方和燃烧室进气侧旁底热疲劳公式简化为圈部分地方出现不同程度的热应力集中现象,这些地方在设计活塞时也要重点考虑。0.347(2)计算结果还表明,该柴油机活塞进行低频热疲劳分析,得出柴油机经历起动-运行-停车图6是活塞仅在温度载荷下产生的热变形。由循环的最危险点的低频热疲劳寿命是112×10次图可知,活塞热变形最大出现在活塞顶面排气口这对柴油机来说寿命足够。侧,达到0.328mm;活塞温度主要集中在头部位置。现选取几何变形比较大的区域点进行计算,结参考文献果如表1所示1肖永宁等.内燃机热负荷和热强度[M].北京:机从表1可以看出在柴油机经历起动-运行-停械工业出版社,1988车的循环次数最少大约是1.12x10次,这对2冯立岩,高希彦,夏惠民等.8E160柴油机活塞组1105W型柴油机的低频热疲劳寿命是非常足够的。热负荷及机械负荷耦合分析.内燃机学报,表1活塞关键点的热疲劳寿命2002,20(5):441-4463谢琰,席明智,刘晓丽.基于 ANSYS的活塞温度关键点位置应变/mm循环次数N场数值模拟研究叮柴油机设计与制造,2009(4)活塞顶面排气侧边缘l.12x10°4雷基林.增压柴油机活塞三维有限元分析及温度0.32796场试验研究D].昆明:昆明理工大学,2005:1-8.活塞顶面进气侧边缘0.305921.29×105 Rodriguez M P, Shammas NY A Finite element燃烧室喉口0.25566184×10°simulation of thermal fatigue in multiplayer and燃烧室凸台0.16766structures:thermal and mechanical approach Ul4.28×10°Microeletronics Reliability. 2001, 41(4): 517-523第1环槽(排气侧)0.25042192×1056张文孝,郭成壁.船舶柴油机活塞的热疲劳强度第1环岸(排气侧)0235212.18×10°分析内燃机学报,2001,18(3):258-2627张卫正,魏春源等.内燃机铝合金活塞疲劳寿命第2环槽(排气侧)0.216722.56×10预测研究.中国机械工程,2003,14(10)第2环岸(排气侧)|0206112.83×10°865-867.活塞裙部(中间)0.172184.06×10°8张文孝.应用当量应变法预测柴油机活塞的多维疲劳寿命、内燃机工程,2002,(23)4活塞裙部下缘0.229612.28×1059万德玉编.柴油机试验测试与分析实用手册[M]活塞内腔顶部中央0.0672226.65×10°北京:学苑出版社,2000活塞销座内侧上缘009762212.63×1010《中国内燃机工业年鉴》编委会.中国2001年内燃机工业年鉴M].上海:上海交通大学出版6结论社,2001.(1)计算结果表明,在标定工况下,活塞最大11[日平修二.热应力与热疲劳[M],郭延玮,李安on mises热应力为684MPa,出现在排气一侧的定译.北京:国防工业出版社,1984中国煤化工CNMHG14-(14)