连杆的接触有限元分析.docxVIP

连杆的接触有限元分析 接触有限元分析 可靠性是可靠性工程的一个重要方面。近年来, 利用有限元方法 (FEM) 对连杆的精细分析取得了一定的研究成果, 但由于连杆是一个较为复杂的装配体, 几何体复杂, 接触对繁多, 载荷也较为复杂, 这就决定了连杆接触有限元分析的复杂性。本文利用ABAQUS对连杆在额定工况下进行接触有限元分析。 一、 模型建立与网格划分 本文在UG中对连杆进行三维建模, 建模时, 删除了一些倒角和油孔, 由于本文接触分析的重点在连杆大头, 对活塞销将在ABAQUS采用分析刚体进行建模, 以简化计算。采用parasolid格式把建好的几何模型导入ABAQUS中, 网格处理采用混合网格划分, 定义材料属性等。建立的有限元模型如图1所示。 二、 定义接触对的类型 本文中的连杆是按实际连杆各个零部件进行建模装配的, 接触关系将从实际情况出发来定义。定义的接触对有:活塞销与连杆小头, 连杆身与连杆盖, 连杆身与上轴瓦, 连杆盖与下轴瓦, 上轴瓦与下轴瓦, 上轴瓦与曲柄销, 下轴瓦与轴柄销, 以及四个螺母与连杆盖, 共计11个接触对。对螺栓与螺纹孔将采用绑定约束进行简化处理。 三、 惯性力的特点 本文计算发动机在标定工况下的连杆载荷, 发动机参数见表1。 连杆在最大拉伸情况下的载荷包括连杆小头所受的活塞惯性力、气体压力、连杆惯性力载荷。连杆运动时的往复惯性力是一种分布的体积力, 由于连杆上各单元的位置不同, 各单元的往复惯性力也不同。在ABAQUS软件中, 可以通过Gravity来定义惯性载荷。惯性载荷不依赖于任何单元, 它通过指定的参考坐标系中加速度来定义。连杆小头载荷aP即活塞组件通过活塞销作用于连杆小头的力。 1. 气体效率 作用于活塞的气体作用力为: 2. 基于数据的4基于加速度的传动 在讨论活塞组件惯性力之前, 首先需计算出活塞的加速度j, 根据文献, 活塞加速度的计算公式如下: 活塞组件中的所有零件, 包括活塞、活塞环、活塞销和活塞销卡环, 均以活塞加速度j做变速往复直线运动。整个活塞组件的惯性力 3. 过连杆质心质量 (1) 因往复加速度而产生的惯性力 此力通过连杆质心C而平行于气缸中心线, 式中ml是连杆的质量。 (2) 因变速摆动中的向心加速度而产生的惯性力 式中ωl——连杆角速度。 4. 螺母提前紧张 参照装配经验, 螺栓预紧力的大小统一取15000N。 四、 单元应力及疲劳系数 本节研究连杆的各部件应力、接触面应力和连杆总体变形, 并对其疲劳强度进行了校核。 图2中, 可以清楚的看到连杆上应力的分布情况, 有三个地方应力比较突出, 第一个是连杆的连杆身部位, 第二个是连杆小头与活塞销接触的部位, 第三个是连杆大头与曲柄销接触的部位。 最大压缩工况下连杆各部件应力分布如图3~6所示。 连杆的材料为42Cr Mo, 其强度极限为1080Mpa;屈服极限为930Mpa。 从输出结果中, 可以找出连杆应力最大的单元为单元32577, 其应力情况如表2所示。 对单元32577来说,, 所以连杆完全能满足静强度设计的要求。 该单元在连杆受最大拉伸载荷时等效应力为86.9706MPa, 其应力情况如表3所示。 根据以下公式来计算其疲劳安全系数: 式中:σa——应力幅 (MPa) , σm——平均应力 (MPa) , σ-1z——对称拉压循环疲劳极限, 取σ-1z=0.41σb; kσ——应力集中系数, 取值为1; εσ——尺寸影响系数, 取值为1; β——工艺系数, 取值为0.6; 由文献, 对于材质均匀, 性能离散程度小, 而应力计算又比较准确的情况, 可取安全系数许可值[n]=1.3~1.5, 甚至更小些。计算得出, n=1.841.5。经比较知, 该点在额定工况下是安全的。 图7中, 可以清楚地看到连杆整体变形的分布情况, 连杆小头部分应变最大, 变形量为5.537E-4m, 沿连杆小头到大头应变逐渐减小, 这是符合应变理论的, 载荷较大的部位的变形往往较大。连杆盖上的变形大于与它接触的连杆身部位的变形。 五、 螺栓与连续剂misen应力平衡校核 最大拉伸工况时, 连杆身最大应力不超过150M Pa, 完全满足静强度的设计要求。与最大压缩工况不同, 最大拉伸工况时, 连杆最大应力不出现连杆身, 而是在连杆螺栓上以及螺栓与连杆接触绑定的部位, 其Mises应力值为284.876M Pa, 单元号为39。对单元号为39的单元进行疲劳强度校核。计算得:n==1.541.5。经比较知该点在额定工况下是安全的。 最大拉伸工况时, 连杆小头部分应变最大, 变形值为5.965E-4m, 沿连杆小头到大头应变逐渐减小。在螺栓连接连杆身与连杆盖的部位变形不一致, 原因就在于此力学模型的建立依据是:连杆是斜切口, 螺栓的变形较