基于有限元分析方法的后桥壳疲劳失效的预测_外文文献doc

1、 毕业设计（论文）译文题目基于有限元分析方法的后桥壳疲劳失效的预测 学生姓名梅月媛 学号 2010800454 专业机械设计制造及自动化 班级20111053 指导教师 杜义贤 评阅教师 完成日期 2015 年1月3日基于有限元分析方法的后桥壳疲劳失效的预测M.M. Topac, H. Gunal, N.S. Kuralay关键字：后桥壳；应力集中；失效；有限元分析； 摘要当施加循环垂直应力在后桥壳上，会产生过早的疲劳失效的现象，之前就已经有了这方面的研究。在这些试验中，裂缝主要出现在样品的同一区域。为了找出这种失效的原因，我们建立了完整的后桥壳CAD模型，并通过材料的拉伸性能试验获得桥壳材料

2、的机械性能。利用这些数据，运用有限元原理进行了应力及疲劳分析。确定了疲劳失效初始发生的位置以及不发生疲劳应变的最小循环垂直应力。而且我们还将有限元分析的结果与实验的结果进行对比，提出了增强桥壳疲劳寿命的解决方案。1.概述由于具有较高的承载能力，固体车桥通常用于重型商用车辆上1。固体车桥的结构如图1所示.在车辆的使用过程中，车桥是主要承载部件，由路面不平产生的动态力进而产生的动态压力导致了车桥产生疲劳破坏。因此，对桥壳抵抗疲劳破坏的寿命进行预测是很重要的。在大规模生产前，需通过垂直力试验对桥壳模型在动态垂直应力作用下进行如图2所示的装载能力及疲劳寿命的分析。在这些测试中，由液压机构提供的循环垂直

3、载荷施加在样品上，直到样品开始出现疲劳破坏。根据承载标准，桥壳必须能承载N=5×循环力而不出现疲劳破坏。在对如图3所示不对称的桥壳模型进行垂直疲劳测试时，当力未达到极限前就有疲劳破坏出现在模型上。因此发现，不出现疲劳破坏的最小循环应力大约为3.7×。在这些测试中，裂纹出现在班卓过渡区E1和E2。疲劳失效的样品如图4所示。为了找出过早失效的原因，我们运用CATIA V5R15商业软件建立了一个详细的桥壳实体模型。并通过该模型，建立有限元模型。在ANSYS V11.0商业有限元分析软件工作平台进行应力和疲劳分析。通过拉伸测试的有限元分析获得了桥壳的材料性能，运用RecurDyn

4、 商业CAE软件进行车辆动态模拟，获得了桥壳最大动态载荷。通过以上分析，找到应力集中部位。为了进行疲劳分析，通过引入疲劳强度修正系数建立了桥壳材料的近似S-N曲线。将分析获得的结果与垂直疲劳测试实验的结果进行比较。为了阻止过早失效并获得增大的疲劳寿命的效果，我们提出了一些解决方案。 图1 商用车后桥壳总成 图2 桥壳模型垂直疲劳测试 图3 桥壳几何形状 图4 测试样品底部的疲劳开裂 图5 桥壳的完整CAD模型2.有限元模型2.1.CAD和有限元模型如图5所示，为用于分析的完整车桥CAD模型。桥壳本质上由两个相同的薄壁壳组成，薄壁壳的厚度为9.5mm并沿着后桥壳的中性轴焊接。在前端面，一个用螺栓

5、固定了差动齿轮装载器的装配环被焊接在桥壳上用来增强刚度。出于密封性的考虑，将一个圆盖焊接在后端面上。这里， A和B代表车辆纵臂纵向推力杆。支撑C和D代表轮与地面的接触。车桥支撑联接点之间的距离与后轴轮轨之间的距离相等。运用CATIA V5R15建立桥壳三维模型。将桥壳的CAD完整模型导入ANSYS Workbench V11.0工作平台前置处理界面，建立分析所需的有限元模型。有限元模型用于图6所示的压力及疲劳分析。为了建立有限元模型，桥壳按照SOLID187进行网格划分。具有二次位移的三维实体单元并且适用于进行不规则网格划分。桥壳被定义为拥有10个节点，且每个节点拥有3个自由度。运用CONTA

6、174和TARGE170元素建立桥壳各部件之间的联系。焊接表面的联接关系选择为完整的可靠联接。有限元模型由779,305个元素和1,287,354个节点组成。 图6 桥壳有限元模型表1S450N的化学特性（wt%）表2抗拉测试结果2.2.桥壳材料车桥壳是由厚的微金属合金管壁经冲压焊接制成的，该管壁的材料为热成型标准钢铁（材料编号1.8901，等同于ISO标准3中的E460）.该材料的化学成分是从供应商提供的的表中获得，具体见表14。未加工的S460N材料的机械性能见参考文献5。然而，在制造过程中，桥壳材料经过了若干道工序，包括退火至800和750热冲压。为了在有限元分析中将工序对机械性能的影响

7、考虑在内并确定处理后材料的精确机械性能，从后桥壳模型中抽取5个样本进行拉伸试验。所有的试验均在室温下进行。从后桥壳模型中抽取的5个样本均在热影响区之外。表2给出的结果为从5个样本中获得的最低值，将这些结果引入有限元模型。并将材料定义为显性各向同向性材料。2.3.负荷条件有限元分析中的负荷条件是根据垂直疲劳测试中出现早期失效处的支撑区域确定的。测试是在如图7所示的可提供80吨载荷的装置上进行的。该装置是由两个具有执行机构和伺服阀的电动液压组成的，伺服阀安装在连接A和B的卡钳处。Ts表示两个卡钳间的距离，Tw代表C与D间的距离，这个距离为真实后桥壳的轮距。车桥的模型是根据如图8所示的由两个空气弹簧

8、支撑的真实桥壳设计的。因为牵引臂的偏心力，所以弹簧的弹力也产生了弯曲应力，该应力在桥壳上产生了一个额外的弯曲变形M。测试样品中的额外弯曲影响由图7所示的液压驱动装置的偏距提供。每个弹簧的最大设计静载荷为F=2850kg。施加在弹簧底座上的力作用于点ZR和ZL。这样就在卡钳A和B处产生了P=4550kg的静态反力。因为路面不平使车身的集中质量产生的垂直加速度导致在每个卡钳处的最大动态载荷大约为P的两倍。由RecurDyn商业CAE软件进行的计算机路面模拟所得的载荷变化范围为182-9100kg。垂直疲劳测试所得的载荷特性曲线如图9所示。有限元分析也考虑到了最大动载荷9100kg沿额外弯曲变形M所

9、产生的影响。如图10所示的车桥垂直应力模型是根据参考文献6设计的。 图7 垂直疲劳测试原理图图8 纵臂的偏心载荷3.有限元分析及结果有限元分析用于预测应力集中及疲劳寿命较低区域的准确位置。P和M施加在图10所示的卡钳连接处。运用装配1.86GHz因特尔至强四核处理器的HPxw8400工作站借助ANSYS V11.0工作平台进行压力分析。图11所示为有限元分析所得的等应力分布图。分析结果显示应力集中区域分布在桥壳承载区域底部的过渡区。从图12中可以看出疲劳失效区域与临界区域在同一位置。计算得出的最大分布应力为max=388.7Mpa；是材料屈服应力点的78.1%。这说明桥壳在承受最大静载荷时符合

10、安全条件。 图9 疲劳测试中的执行机构负荷特性曲线图10 桥壳的外加负荷及弯矩图图11 下壳体上的工作应力分布图12 测试与分析结果比较4.疲劳寿命预测由于在使用中后桥壳承受动应力，也需要进行疲劳分析。压力寿命的疲劳极限估计值为 (1)钢材的强度极限Sut小于1400MPa7,8。这意味着疲劳强度的周期为或更多。为了预测在周期范围内的疲劳寿命，使用参考文献9中使用简单抗拉测试获得所需数据的方法作出桥壳材料的S-N曲线。Se代表理想实验样品的压力疲劳寿命。为了预测机械零件的真实疲劳强度，需要乘上代表各种设计，制造和环境对疲劳强度影响的修正因子10。为 (2)式中为根据下式得出的表面抛光度得到的表

11、面因数 (3)由于桥壳表面的粗糙度与经过热冲压工艺的热轧钢板相似，所以推荐的标准为a=57.7和b=-0.7187.经计算得出Ka=0.564，Sut=629.9MPa。另外，喷丸工艺作为一种常见的用于减少零件材料表面残留应力的方法，也用于增加热冲压后的桥壳表面的疲劳寿命。文献9中给出这种方法可增加70%的疲劳寿命。因此，在有限元分析中的取值为0.959.因为桥壳为非圆形截面，根据横截面深度h远大于50mm假定尺寸因数为0.75。由于环境温度T=0-250，所以弯曲和环境因数Kc=1,进而确定负荷系数Kd=1。通过静态有限元分析，可得出应力集中区分布在班卓及横臂过渡区域。所以，除了上述修正因数

12、外，疲劳强度修正因数必须引入分析，可通过与应力集中系数有关的应力集中系数得到。因此的计算式为 (4) 出于安全考虑，假设Kf与Kt相等7。由于桥壳的大小及形状的复杂性，Kt无法从标准文献中查出。另一方面，Kt被定义为 （5）式中为p凹口处得峰值应力，n为不出现应力集中时的常应力9,12，p的使用数值可从max=388.7MPa时的静有限元分析中得出。为了计算将后桥壳简化为一简支梁，其沿纵轴Y的危险横截面-都为矩形并适用于纯弯曲理论6。按图10所给出的模型的计算公式为 (6)式中M为弯曲力矩，Z为危险横截面的断面系数。M的取值为41.9× Nmm。断面系数Z取值为127507m。因此计

13、算得出为n=329MPa。发现KtKf=1.181，Ke=0.846。运用ANSYS V11.0工作平台定义S-N曲线中标绘的修正因数。通过压力寿命决定桥壳材料的疲劳寿命。全部的疲劳分析都是以无限寿命进行的（N=）。用有限元分析得到的压力分布图进行疲劳寿命计算。由于载荷具有正弦波动特性（平均应力m>0），修正方法如下9 (7)式中，n表示安全系数。振幅为amaxmin2 (8) （8）平均应力可表示为m(maxmin) 2 (9)式中，通过有限元分析得到max为最大值9100kg,匹配的min最小值为182kg。壳体底部的分配系数n如图13所示。根据疲劳分析结果，估计在周期为.3.6&#

14、215;时桥壳表面区域会发生裂纹开裂，该数值低于预测值为5×周期的最小疲劳寿命。此处n的最小值为0.93.在桥壳的内表面，最大应力集中发声处区域的n值最小计算结果为0.767。这意味着，在垂直应力测试中区域和会在载荷周期5×前发生疲劳开裂。5.结构及讨论有限元分析显示在垂直疲劳测试中出现疲劳破坏的区域存在应力集中，该应力集中会导致在最小预测周期5×前出现过早破坏。此结果与垂直疲劳试验中的结构相同。增大桥壳的疲劳寿命需减小应力集中。减小应力集中，增大疲劳寿命的最简单的方法是金属壁的厚度。然而，在区域外桥壳符合无限寿命周期条件。增加金属笔厚度导致了不必要的重量增加。例

15、如，增加厚度0.5mm，使得桥壳材料在临界区域的疲劳极限提高到了超过5.85×周期，此极限超过了设计的疲劳极限。另一方面，这也意味着提高了汽车非簧载质量5%的重量。所以这并不是实用的解决方法。作为另一种解决方法，可从新设计过渡区域的几何形状。平整的过渡区几何形状可提高疲劳寿命而不增加重量。此外，加固环的形状也对应力集中产生影响。在所研究的该桥壳设计中，加固环的厚度为20mm。为了预测加固环的影响，在没有加固环的情况下又进行了一次有限元分析。在临界区域处的最大分布应力为428MPa。这意味着，实用加固环大约减少了10%的应力集中。通过增加此部分的厚度，可能会增加硬度。在此设计中，由于动

16、力系统外形的限制，增加的厚度为5mm。根据此加固环的外形变化进行静态疲劳分析。然而，分析显示疲劳强度的增加均为其自身的，因此桥壳的疲劳寿命不会增加到超过设计最小载荷周期5×倍的程度。因此，增加加固环的厚度可与从新设计过渡区几何形状同时使用。 图13 下壳体安全系数分布6.总结运用有限元分析方法对卡车后桥壳模型的早期疲劳失效进行分析。在分析中，通过模拟垂直疲劳试验过程，预测应力集中区在班卓过渡区域。发生疲劳开裂的区域与分析所得结果相吻合。通过有限元分析可预测破坏发生的位置。通过稳态和循环张应力确定临界区域。裂缝导致破坏发生在桥壳的应力集中区域。尽管桥壳模型负荷最大垂直载荷静态忍耐条件，

17、分析显示，如果为循环载荷，疲劳破坏可能在预测的最小周期5×前发生。有限元分析同样可用于估计疲劳失效开始前的周期数。为了解决该问题，增加金属管壁的厚度因为会增加桥壳的重量，所以并不是实用的方法。重新设计班卓过渡区和增加加固环的厚度，这种符合最小设计准则的途径，也许是增强疲劳寿命的好方法。感谢这篇论文在土耳其伊兹密尔市的Ege Endustri ve Ticaret A.S.的帮助下完成。作者同时也对来自Dokuz Eylul大学的 E. Çnar Yeni博士和Pamukkale大学的Cemal Meran博士的批评与建议表示感谢。参考文献1 Reimpell J, Stol

18、l H, Betzler JW. The automotive chassis: engineering principles. Butterworth-Heinemann; 2002. p. 39.2 ANSYS Theory Reference. ANSYS Release 10.0. ANSYS, Inc.; 2005.3 Jiang Y, Hertel O, Vormwald M. An experimental evaluation of three critical plane multiaxial fatigue criteria. Int J Fatigue 2007;29:1490502.4 Erdemir Product Catalogue, Eregli Iron and Steel Works Co.; 2007. p.50 in Turkish.5 Hoffmeyer J, Döring R, Seeger T, Vormwald M. Deformation behaviour, short crack growth and fatigue lives under multiaxial nonproportional loading.Int J Fatigue