基于瞬态动力学的航空活塞发动机连杆修理检测区域划分

第27卷第2期燃气涡轮试验与研究Vol 27. No. 22014年4月Gas Turbine Experiment and Researck基于瞬态动力学的航空活塞发动机连杆修理检测区域划分丁发军,魏武国.中国民用航空飞行学院飞机修理厂,四川广汉6183072.中国民用航空飞行学院航空工程学院,四川广汉618307)摘要:基于瞬态动力学基本理论,选取 Lycoming航空活塞发动机的普通连杆为分析对象,在通用有限元软件中建立其三维有限元模型。考虑在发动机额定转速工况下,连杆在压缩行程末期、膨胀行程初期气缸内气体压力作用下的动力学响应。分析发现:在所有载荷步下,最大等效应力和应变均出现在杆身和小头的过渡区域,且与缸内气体峰值压力同时出现;获得的最大等效应力随时间的变化曲线,可为预估连杆疲劳寿命、提高连杆疲劳可靠性提供数值依据;结合等效应力应变云图,划分出的连杆应力集中的危险区域,为修理前检测的重点区域。关键词:航空活塞发动机;连杄;瞬态动力学分析;有限元;等效应力应变;检测区域中图分类号:V234文献标识码:A文章编号:1672-2620(2014)02-0025-05Repair and Detection Area Dipartition of Aviation Piston EngineConnecting Rod Based on Transient Dynamic AnalysisDING Fa-jun, WEI Wu-guo(1. Aircraft Repair Overhaul Plant, CAFUC, Guanghan 618307, China;2. Aviation Engineering Institute, CAFUC, Guanghan 618307, China)Abstract: Based on the basic theory of transient dynamics, a general connecting rod from a Lycoming piston engine had been taken as analysis object, and a three-dimensional finite element model of the rod wasestablished in general finite element software according to initial design parameters. And then, the dynamicresponse of the rod to gas pressure in cylinder during final phase of compression stroke and initial phase ofexpansion stroke under engines rated speed conditions had been analyzed. Simulation results show that themaximum equivalent stress occurs at the transition zone between shaft and little head in all load steps, andwith the appearance of gas peak pressure in cylinder at the same time. And the maximum equivalent stressversus time curve was obtained, which provided numerical basis for estimating fatigue life and improving fatigue reliability of the rod Finally, hazardous areas of the rod had been identified according to equivalentstress and strain contours, and become the key areas of detection before repairKey words: aviation piston engine; connecting rod; transient dynamic analysis; finite elementequivalent stress and strain; detection area1引言易引起连杆疲劳破坏。文献[1]、[2]基于活塞发动机连杆是活塞发动机曲拐机构的重要连接件,工连杆故障统计数据,深入研究了连杆的失效规律和作中受力复杂,承受着急剧变化的动载荷括活塞机理,证明疲劳是连杆断裂的主要原因。传来的气体压力、往复运动质量惯性力及摆动运动目前,国内关于航空活塞发动机连杆疲劳失效质量惯性力)。这些载荷的大小和方向周期性变化,的文献较少,较多文献集中在车用活塞发动机连杆失效问题上。文献[8在分析连杆失效的内在机收稿日期:2013-12-04;修回日期:2014-04-10作者简介:丁发军(1973-),男,四川广汉人,高级工程师,硕士,主要从事航空器维修与适航、可靠性管理,航空活塞发动机及其部件维修、故障诊断与状态监控等方面的研究。燃气涡轮试验与研究第27卷理时,利用 Paris公式简单估算了失效连杆的疲劳寿式中:6a为 Newmark积分参数,△t为时间间隔。命;文献[9从表面工艺和力学角度讨论了连杆疲劳根据方程(1),tn+时刻的平衡方程为失效的影响因素,并基于S-N法、 Goodman平均应力Mln}+C1n+区K]1n,}={F"}(4)修正及 Miner线性累积损伤理论,估算了连杆在为求解un1,将方程(2)、(3代入方程(4中,得到近似正弦载荷下的疲劳寿命。以上分析基本都忽略了连杆实际受载情况,特别是活塞发动机每一工作(aol]+a, [C]+[K Dun+1=F.循环,气缸内燃气压力对连杆造成的冲击影响。在[MAo(un )+a2(i, )+a,(i,)连续进行进气、压缩、膨胀、排气四个行程时,气缸中C1a,un )+au)as(i)的燃气压力不断随曲轴转角变化,通过活塞作用在式中连杆上就形成了一个具有固定周期(固定转速)连续变化的载荷谱,对连杆的疲劳寿命起决定作用航空活塞发动机作为大部分通用飞机的动力裝利用方程(5)求出tn时刻的位移{un}后,就可置,其连杆结构更轻。为保证飞行安全,适航当局对连杆等受力关键件的可靠性有着严格的适航要求。利用方程2)、(3)求出tn+时刻的加速度{n,}和速度一方面,航空活塞发动机仍采用定时下发策略,部分{n}零件存在不到寿换件的情况。为节省维护费用,需根据 Zienkiewicz的理论,利用方程(2)和(3)得到过渡到视情下发策略上来,有必要对受力关键件如的 Newmark求解方法的无条件稳定必须满足连杆的疲劳寿命进行估算。另一方面,发动机大修α≥1(1+8),8≥1,1时,对连杆表面状况进行检查,对防止连杆疲劳断裂222++a×>0Newmark参数可根据下式输入也具有重要意义。但发动机厂家目前提供的资料中,没有连杆表面质量检查区域的规定。a=1(1+y)本文基于瞬态动力学理论,对航空活塞发动机式中:y为振幅衰减因子。由式(6和式(7)可知,只要连杆在燃气压力作用下的动力学响应进行分析,以y≥0,求解就是稳定的。期得到准确的应力时间历程,为准确预估连杆的疲2.2有限元分析模型的建立劳寿命提供数值依据;同时确定岀连杆受力的薄弱选取通航领域使用较多的某型 Lycoming水平环节,以便最终制定科学合理的检测修理策略。对置型发动机的普通连杆为分析对象,主要考虑压缩行程末期、膨胀行程初期,上死点前后气缸中气体2连杆瞬态动力学的有限元分析方法压力对连杆受载的影响,并在通用有限元软件AN2.1瞬态动力学分析基本理论SYS Workbench13.0中进行瞬态动力学分析。对于线性结构的瞬态动力学平衡方程(1),可使连杆有限元建模分两步进行:首先根据图1中用 Newmark积分法求解。的初始设计参数建立实体模型,在不影响计算精度}+Ca}+k]an}={F"}1)的情况下,忽略连杆小头和大头衬套,大头端孟,螺式中:[M]为系统质量矩阵,C]为系统阻尼矩阵,[K为系统刚度矩阵,{F}为外力函数向量,为节点位移向量,{i}为节点速度向量,{i}为节点加速度向SECTION A-A量171.5Newmark方法使用有限差分法,在一个时间间隔内有{n}={a}+[-6){an}+6{n1△n,=+1+[-a)1+a1图1连杆几何尺寸(单位:mm)Fig 1 Dimension of the rod(Unit: mm第2期丁发军等:基于瞬态动力学的航空活塞发动机连杆修理检测区域划分栓,简化分析模型;然后建立三维有限元模型。取连在上死点前20°时,缸内气体压力开始急剧攀升;上杆材料为40 CrMnMo,采用三维10节点高阶四面死点后曲轴转角15°位置为峰值压力点,缸内气体压体单元 Solid187和三维20节点高阶六面体单元Sol-力峰值达到3.923MPa;从峰值时刻往后,气体压力id186混合划分网格,单元尺寸取4mm;材料常数及又开始急剧下降有限元模型参数见表1。边界条件设置:连杆在工选取该型发动机连杆在额定转速2400r/min作中主要受压应力作用,在大头孔内表面施加固定时,压缩行程末期(上死点前20°)至膨胀行程初期(上约束,在小头圆心指向大头圆心方向上设定周期载死点后30)的曲轴转角范围,缸内气体压力通过活荷。连杆的有限元模型如图2所示。塞对连杆的动力学激励为求解对象。每隔5°曲轴转表1材料常数及有限元模型参数角取1个,共计11个载荷计算点,如表2所示。表Table 1 Material property and finite element model parameters中,上死点前曲轴转角为负,上死点后曲轴转角为正。由于该段曲轴转角范围内连杆受载的变化最剧材料常有限元模型参数烈,计算分析此段曲轴转角范围内连杆应力与变形弹性模盛泊松比密度/(km)抗拉强度屈服强单元数目节点数目的分布规律,能在一定程度上反映连杆在动载下的GPaMPH /MPa2100.37800>850>6502243639215工作状况表2载荷计算点Table 2 Computational load-10O上死点)p/MPa0.7350.83409321,0791.2261.7161015(峰值压力)a3.4323.923,20.22图2连杆的有限元模型312载荷步设置Fig 2 Finite element grid of the rod为对连杆进行瞬态动力学分析,需将表2中缸内气体压力随曲轴转角的变化转换为连杆受力随时3连杆的瞬态动力学分析间的变化,即连杆的载荷步与载荷值。取曲轴转速为额定转速2400r/min,设曲轴转3.1计算工况选择及载荷步设置角等于-20°时为初始时刻t=0,曲轴每转过5°需要时3.1.1计算工况选择间0.347ms(即△t=0.347ms),可依次算出11个载荷该型发动机为自然吸气式发动机,提前点火的计算点所对应的载荷步时刻,如表3所示。曲轴转角为20°。当在额定转速2400rmin时,缸内气体压力p随曲轴转角的变化曲线如图3所示。表3载荷步与载荷值Table 3 Load steps and load values03.432载荷P2N1133-12419-143632.942间/ms1.389(上死点)1.7362.0832.431峰值压力2.452载荷pN16315850-457535239寸间ms2.77833.472147I载荷p,N-4859143480-396790.980.490单个气缸中的曲拐机构组件,在上死点前后附上死点0.0近时的受力如图4所示2。根据 Lycoming发动机资料31,取曲臂半径R=55.5625mm,连杆小头圆心到图3缸内气体压力随曲轴转角的变化大头圆心距离L=171.5000mm,活塞顶面直径D=Fig 3 Gas pressure vs. crank shaft turn angle130.1750mm。在表2中各曲轴转角下,首先将已知燃气涡轮试验与研究第27卷的缸内气体压力乘以活塞顶面积,计算出作用在活N塞顶的气体作用力;然后进行如图4的受力分析,并利用三角函数公式,得到连杆在各曲轴转角下(对应各载荷步时刻)的载荷值,如表3所示,载荷为负表示连杆受压。上死上死点-4(a)等效应力云图下死点7图4上死点前后曲拐机构组件受力示意图Fig, 4 Force anad top dead center(b)等效应变云图3.2计算结果分析图6连杆在峰值压力点时的应力应变分布有限元计算分析发现,在所有载荷步下,最大等Fig6 Stress and strain distribution of the rod when(=2.431m效应力和应变均出现在杆身与小头的过渡区域,且死点和缸内气体出现峰值压力时的等效应力应变云应力集中比较严重。图5、图6分别示出了连杆在上图。可见,当曲拐机构运动到压缩行程末期上死点时,等效应力和应变均不是最大值;最大值出现在曲轴运动到上死点后峰值压力点时,此时的最大等效应力达32568MPa,但小于连杆材料的屈服强度,因此连杆不会发生静载塑性变形连杆在各载荷步下的最大等效应力随时间的变化曲线如图7所示。可见,其与气缸内气体压力随曲轴转角的变化曲线类似。在整个加载过程中,当气缸内气体压力增加到峰值压力时,杆身与小头的过渡区域同时出现最大等效应力。本曲线可为准确)等效应力云图预估连杆的疲劳寿命、分析连杆疲劳可靠性提供数值依据271.89218.10l64.3l10.5256.732(b)等效应变云图3.472t/ms图5连杆在上死点时的应力应变分布图7最大等效应力随时间的变化ig. 5 Stress and strain distribution of the rod when t=1. 389 msF第2期丁发军等:基于瞬态动力学的航空活塞发动机连杆修理检测区域划分3.3修理检测区域划分空活塞发动机的普通连杆为分析对象,在通用有限根据图5、图6的等效应力应变云图,结合连杆元软件中建立起连杆的三维有限元模型。在发动技术条件,可确定杆身和小头的过渡区域为危险区机额定转速工况下,考虑连杆在缸内气体压力作用域,如图8所示。图中红色区域为修理前检测的重下的动力学响应,经仿真分析得到以下结论:点区域,不允许有任何腐蚀及损伤坑点,否则连杆1)工作中,受缸内气体压力作用,最大等效必须报废。而杄身的大部分区域确定为可修理区应力和应变出现在杆身与小头的过渡区域,且最大域,允许一定程度的损伤。等效应力和缸内气体峰值压力同时出现。(2)计算得到的连杆最大等效应力随时间的变化曲线,为预估连杆疲劳寿命、提高疲劳可靠性提供了数值依据。(3)结合计算结果,划分出连杆的危险区域和可修理区域,为制定连杆修理标准打下了基础。参考文献图8连杆应力集中的危险区域(红色区域)王银燕,张鹏奇,王善.柴油机连杆杆身疲劳强度可Fig 8 The stress concentration area of connecting rod (red靠性分析[哈尔滨工程大学学报,2001,221):67-2]古莹奎,唐淑云,周志博.车用发动机连杆的失效分析另外,连杆在发动机工作中除了承受高频交变机械科学与技术,2011,302):256-25应力外,还要承受冲击载荷、热循环载荷作用,极易阝3]陆有,吕晓霞,安健,等.汽车连杆失效分析汽导致连杆弯曲变形和扭曲变形。如果其变形超过车工艺与材料,2003,3:31-32定极限,则连杆只能报废,不能采取校正的办法4叶明.柴油机连杆断裂失效分析热处理,2002,17修复。在发动机翻修时,般采用图9所示的平行度和垂直度检查方法,来检测连杆的弯曲变形和扭5高学敏,于志伟,许晓磊42CMoA钢连杆断裂失效分曲变形。析J金属热处理,2008,335):100-1016朱红娟.16V280ZJ型柴油机连杆断裂失效分析J内燃机车,2008,7:34-367]王荣滨.连杆失效分析与影响调质钢力学性能的因素热加工工艺,2006,35(12):73-758]张进,赵晓博,陶世杰,等.发动机连杆失效分析铸造设备与工艺,2010,4:394219]吴波,谭建松,胡定云,等.柴油机连杆疲劳失效的影响因素分析[机械制造,2010,48(9):60-62.I温秉权,王宾,路学成,等.金属材料手册M!北京电子工业出版社,2009图9连杆的平行度和垂直度检查!]徐燕茹.基于有限元法的柴油机连杆动力学响应分析与研究[D太原:中北大学,2012Fig9 Checking the parallelism and perpendicularity2唐庆如.活塞发动机[M北京:兵器工业出版社,2007[13 Textron Lycoming Engine Overhaul Manual 60294-7MI4结论2013.基于瞬态动力学基本理论,以某型 Lycoming航1411cextron Lycoming Engine Operators Manual 60297-12[M.2013.