关 键 词：

起动 箱体 结构设计 计算 三维 ug cad 图纸 以及 说明书 仿单

资源描述：

内容简介：

编号： 毕业设计(论文)外文翻译（译文）院 （系）： 机电工程学院 专 业： 机械设计制造及其自动化 学生姓名： 伍荣展 学 号： 1000110131 指导教师单位： 桂林电子科技大学 姓 名： 高成 职 称： 助理研究员 2014 年 05 月 26 日桂林电子科技大学第 1 页 共 10 页后桥壳疲劳失效的有限元分析预测M.M. Topac, H. Gunal, N.S. Kuralay. Fatigue failure prediction of a rear axle housing prototype by using finite element analysisJ. Engineering Failure Analysis，摘 要对与在试验中，当施加循环垂直应力在后桥壳上，产生了过早的疲劳变形的现象 进行了研究。发现在这些试验中，裂缝主要出现在样品的同一区域。为了确定破坏的原因，建立了完整的后桥壳 CAD 模型。同时，桥壳的机械性能取决于其材料的拉伸性能。利用这些数据，运用有限元原理进行了应力及疲劳分析。确定了疲劳应变的发生位置以及不发生疲劳应变的最小循环垂直应力。将有限元分析的结果与实验的结果进行对比。设计提出了增强桥壳疲劳寿命的解决方案。关键字：后桥壳；应力集中；失效；有限元分析1 概述由于具有较高的承载能力，固体车桥通常用于重型商用车辆上1。固体车桥的结构见图 1。在车辆的使用中，车桥是主要承载部件，由路面不平产生的动态应力进而产生的动态压力导致了车桥产生疲劳破坏。因此，最重要的是进行桥壳抵抗疲劳破坏的寿命预测。在大规模生产前，有必要对桥壳模型在动态垂直应力作用下进行如图 2 所示的装载能力及疲劳寿命的有限元分析。在这些测试中，由液压机构提供的循环垂直载荷施加在样品上，直到样品出现疲劳破坏。根据承载标 准，桥壳必须能承载N=5X105 循环应力而不出现疲劳破坏。在对如图 3 所示不对称的桥壳模型进行垂直疲劳测试时，在应力达到极限前就有疲劳破坏出现在模型上。因此发现，不出现疲劳破坏的最小循环应力大约为 3.7X105。在这些测试中， 裂纹出现在班卓过渡区 E1 和E2。图 4 所示为一早期破坏的例子。为了找出早期失效的原因，运用 CATIA V5R15 商业软件建立了一个详细的桥壳三维模型。利用该模型，建立有限元模型。运用 ANSYS V11.0 商业有限元分析软件工作平台进行应力和疲劳分析。通过拉伸测试的有限元分析获得了桥壳的材料性能，运用RecurDyn 商业 CAE 软件进行车辆动力学模拟，获得了桥壳最大载荷。通过这些分析，找到应力集中部位。为了实现疲劳分析，引入疲劳强度修正系数建立了桥壳材料的估计 S-N 曲线。将分析获得的结果与垂直疲劳测试实验的结果进行比较。为了阻止早期破坏并获得增大的疲劳寿命，提出了一些解决方案。 桂林电子科技大学第 2 页 共 10 页图 1 商用车后桥壳总成图 2 桥壳模型乘直疲劳测试图 3 桥壳几何形状桂林电子科技大学第 3 页 共 10 页图 4 测试样品底部的疲劳开裂图 5 桥壳的完整 CAD 模型2 有限元模型2.1 CAD 和有限元模型分析用全尺寸车桥 CAD 模型如上图 5 所示。桥壳本质上由两个相同的薄壁壳组成，薄壁壳的厚度为 9.5mm 并沿着后桥壳的中性轴焊接。在前端面，一个用螺栓固定了差动齿轮装载器的曾环被焊接在桥壳上用来增强刚度。出于密封性的考虑，将一个圆盖焊接在后端面上。这里，元素 A 和 B 为下垂壁卡钳联接。支撑 C 和 D 代表轮与地面的接触。车桥支撑联接点之间的距离与后轴轮轨之间的距离相等。运用 CATIA V5R15建立桥壳三维模型。将桥壳的完整 CAD 模型导入 ANSYS V11.0 工作平台前置处理界面，建立分析所需的有限元模型。有限元模型用于图 6 所示的压力及疲劳分析。为了建立有限元模型，桥壳按照 SOLID187 进行网格划分。S0LID187 具有二次位移的三维实体单元并且适用于进行不规则网格划分。桥壳被定义为拥有 10 个节点，且每个节点拥有 3 个自由度。运用 CONTA174 和 TARGE170 元素建立桥壳各部件之间的联系。焊接表面的联接关系选择为完整的可靠联接。有限元模型由 779,305 个元素和 1,287,354个节点组成。 桂林电子科技大学第 4 页 共 10 页图 6 桥壳有限元模型表 1 S450N 的化学特性（ Wt%）表 2 抗拉测试结果2. 2 桥壳材料车桥壳是由 9.5mm 厚的微金属合金管壁经冲压焊接制成的，该管壁的材料为热成型标准钢铁 S460N （材料编号 1.8901，等同于 ISO 标准3中 E460） 。该材料的化学成分是从供应商获得的，具体见表 1 4。未加工的 S460N 的机械性能见参考文献5。然而，桥壳材料在制造过程中需经过若干道工序，包括退火至 800C 和 750C 热冲压。为了将工序对机械性能的影响引入有限元分析并确定加工后材料的精确机械性能，从后桥壳模型中抽取 5 个样本并进行拉伸试验。所有的试验均在室温下进行。从后桥壳模型中抽取的 5 个样本均在热影响区之外。表 2 给出的结果均为 5 个样本的最低值，并将这些结果引入有限元模型。将材料定义为显性各向同向性材料。2. 3 负荷条件有限元分析中的负荷条件是根据垂直疲劳测试中出现早期失效处的支撑区域确定的。测试是在如图 7 所示的可提供 80 吨载荷的装置上进行的。该装置是由两个具有承载单元的电动液压执行机构和伺服阀组成的，伺服阀安装在连接 A，B 的卡钳处。TS表示两个卡钳间的距离，T 因表示支撑 C，D 间的距离即真实后桥壳的轮距。车桥的模型是根据如图 8 所示的由两个空气弹簧支撑的真实桥壳设计的。因 为载荷施加在牵引臂的偏心轮上，所以弹簧的弹力也产生了弯曲应力，该应力在桥壳上产生了一个额外的弯曲 AM。测试样品中的额外弯曲影响由图 7 所示的液压驱动装置的偏距 c 提供。每个弹簧的最大设计载荷为 F = 2850kg。应力垂直的施加在弹簧底座的点。这导致了桂林电子科技大学第 5 页 共 10 页在卡钳 A, B 处产生了静态反应力 P=4550kg。因为路面不平使车身的集中质量产生的垂直加速度导致在每个卡钳处的最大动态载荷大约为 P 的两倍。由 ReoirDyn 商业 CAE软件进行的计算机路面模拟所得的载 荷变化范围为 182-9I00kg。垂直疲劳测试所得的载荷特性曲线如图 9 所示。有限元分析也考虑到了最大动载荷 9100kg 沿额外弯曲变形M 所产生的影响。如图 10 所示的车桥垂直应力模型是根据参考文献 6设计的。图 7 垂直疲劳测试原理图图 8 纵臂的偏心载荷桂林电子科技大学第 6 页 共 10 页3 有限元分析及结果有限元分析用于预测应力集中及疲劳寿命较低区域的准确位置。P 和M 施加在图 10 所示的卡钳连接处。运用装配 1.86GHz 因特尔至强四核处理器的 HPx 因 8400 工作站借助 ANSYS V11.0 工作平台进行压力分析。图 11 所示为有限元分析所得的等应力分布图。分析结果显示应力集中区域 F1、F 2 分布在桥壳承载区域底部的过渡区。从图 12 中可以看出疲劳失效区域与临界区域在同一位置。计算得出的最大分布应力为 max=388.7Mpa；是材料屈服应力点的 78. 1%。这说明桥壳在承受最大静载荷时符合安全条件。图 9 疲劳测试中的执行机构负荷特性曲线图 10 桥壳的外加负荷及弯矩图图 11 下壳体上的工作应力分布桂林电子科技大学第 7 页 共 10 页图 12 测试与分析结果比较4.疲劳寿命预测由于在使用中后桥壳承受动应力，也需要进行疲劳分析。压力寿命的疲劳极限估计值 为s（1）s=0.504钢材的强度极限小于 1400MPa7，8。这意味着疲劳强度的周期为 106 或更多。为了预测在 105 - 106 周期范围内的疲劳寿命，使用参考文献9中使用简单抗拉测试获得所需数据的方法作出桥壳材料的 S-N 曲线。代表理想实验样品的压力疲劳寿命。为了预测机械零件的真实疲劳强度 ， 需s s要乘上代表各种设计，制造和环境对疲劳强度影响的修正因子 10。S e 为 Se=kakbkckdke （2）s式中 ka 为根据下式得出的表面抛亮度得到的表面因数ka=a （3）由于桥壳表面的粗糙度与经过热冲压工艺的热轧钢板相似，所以推荐的标准为 a=57.7 和 b=-0.7187.经计算得出 ka=0. 564, Sut=629.9MPa。另外，喷丸工艺作为一种常见的爪于减少零件材料表面残留应力的方法，也用于增加热冲压后的桥壳表面的疲劳寿命。文献9 中给出这种方法可增加 70%的疲劳寿命。因此，在有限元分析中 ka 的取值为 0.959。因为桥壳为非圆形截面，根据横截面深度 h 远大于 50mm 假定尺寸因数kb 为 0.75。由于环境温度 T=0-250C，所以弯曲和环境因数 kd=1，进而确定负荷系数 kc=1。通过静态有限元分析，可得出应力集中区分布在班卓及横臂过渡区域。所以，除了上述修正因数外，疲劳强度修正因数 ke 必须引入分析，k e 可通过与应力集中系数 kf有关的应力集中系数 kt 得到。因此 ke 的计算式为ke=1/kf （4）桂林电子科技大学第 8 页 共 10 页出于安全考虑，k f 假设与 kt 相等 7。由于桥壳的大小及形状的复杂性，k t 无法从标准文献中查出。另一方面，k t 被定义为Kt= peak/ nominal （ 5 ）式中 peak 为凹口处得峰值应力， nominal 不出现应力集中时的常应力 p9, 12， peak 的使用数值可从 max=388. 7MPa 时的静有限元分析中得出。为了计算 nominal 将后桥壳简化为一简支梁，其沿纵轴 Y 的危险横截面 X1X1 都为矩形并适用于纯弯曲理论6。 nominal 按图 10 所给出的模型的计算公式为 nominal=M/Z （6）式中 M 为弯曲力矩，Z 为危险横截面的断面系数。M 的取值为 41.9xl06Nmm。断面系数 Z 取值为 127507mm3。因此计算得出 nominal 为 329MPa。发现kt kf=1.181，k e=0.846。运用 ANSYS V11.0 工作平台定义 S-N 曲线中标绘的修正因数。通过压力寿命决定桥壳材料的疲劳寿命。全部的疲劳分析都是以无限寿命进行的（N=10 6） 。用有限元分析得到的压力分布图进行疲劳寿命计算。由于载荷具有正弦波动特性（平均应力 m0） ，修正方法如下9（7）+=1式中，n 表示安全系数。振幅 为（8）=2（8）平均应力 am 可表示为（9）=+2式中，通过有限元分析得到 max 为最大值 9100kg， min 匹配的最小值为 182kg。壳体底部的分配系数 n 如图 13 所示。根据疲劳分析结果，估计在周期为 ca3.6X105 时，桥壳表面 F1 区域会发生裂纹开裂，该数值低于预测值为 5x105 周期的最小疲劳寿命。此处 n 的最小值为 0. 93。在桥壳的内表面，最大应力集中发生处 F2 区域的 n 值最小，计算结果为 0.767。这意味着，在垂直应力测试中区域 F1 和 F2 会在载荷周期 5X105 前发生疲劳幵裂。5.结构及讨论有限元分析显示在垂直疲劳测试中出现疲劳破坏的区域存在应力集中，该应力集中会导致在最小预测周期 5X105 前出现过早破坏。此结果与垂直疲劳试验中的结构相同。增大桥壳的疲劳寿命需减小应力集中。减小应力集中，增大疲劳寿命的最简单的方法是金属壁的厚度。然而，在 F1F2 区域外桥壳符合无限寿命周期条件。增加金属笔厚度导致了不必要的重量增加。例如，增加厚度 0.5mm，使得桥壳材料在临界区域的桂林电子科技大学第 9 页 共 10 页疲劳极限提高到了超过 5. 85X105 周期，此极限超过了设计的疲劳极限。另一方面，这也意味着提高了汽车非簧载质量 5%的重量。所以这并不是实用的解决方法。作为另一种解决方法，可从新设计过渡区域的几何形状。平整的过渡区几何形状可提高疲劳痔命而不增加重量。此外，加固环的形状也对应力集中产生影响。在所研究的该桥壳设计中，加固环的厚度为 20mm。为了预测加固环的影响，在没有加固环的情况下又进行了一次有限元分析。在临界区域 F2 处的最大分布应力为 428MPa。这意味着，实用加固环大约减少了 10%的应力集中。通过增加此部分的厚度，可能会增加硬度。在此设计中，由于动力系统外形的限制，增加的厚度为 5mm。根据此加固环的外形变化 进行静态疲劳分析。然而，分析显示疲劳强度的增加均为其自身的，因此桥壳的疲劳寿命不会增加到超过设计最小载荷周期 5X105 倍的程度。因此，増加加固环的厚度可与从新设计过渡区几何形状同时使用。图 13 下壳体安全系数分布6.总结运用有限元分析方法对卡车后桥壳模型的早期疲劳失效进行分析。在分析中，通过模拟垂直疲劳试验过程，预测应力集中区在班卓过渡区域。发生疲劳开裂的区域与分析所得结果相吻合。通过有限元分析可预测破坏发生的位置。通过稳态和循环张应力确定临界区域。裂缝导致破坏发生在桥壳的应力集中区域。尽管桥壳模型负荷最大垂直载荷静态忍耐条件，分析显示，如果为循环载荷，疲劳破坏可能在预测的最小周期 5X105 前发生。有限元分析同样可用于估计疲劳失效开始前的周期数。为了解决该问题，増加金属管壁的厚度因为会增加桥壳的重量，所以并不是实用的方法。重新设计班卓过渡区和增加加固环的厚度，这种符合最小设计准则的途径，也许是增强疲劳寿命的好方法。 桂林电子科技大学第 10 页 共 10 页感谢这篇论文在土耳其伊兹密尔市的 Ege Endustri ve Ticaret A.S.的帮助下完成。作者同时也对来自 Dokuz Eylul 大学的 E. Cmar Yeni 博士和 Pamukkale 大学的 Cemal Meran 博士的批评与建议表示感谢。译文原文出处：M.M. Topac, H. Gunal, N.S. Kuralay. Fatigue failure prediction of a rear axle housing prototype by using finite element analysisJ. Engineering Failure Analysis，(16)2009,:1474-1482.编号： 毕业设计(论文)外文翻译（译文）院 （系）： 机电工程学院 专 业： 机械设计制造及其自动化 学生姓名： 伍荣展 学 号： 1000110131 指导教师单位： 桂林电子科技大学 姓 名： 高成 职 称： 助理研究员 2014 年 05 月 26 日桂林电子科技大学第 1 页 共 9 页利用重型载货汽车的有限元应力分析的初步数据预测其疲劳寿命Roslan Abd Rahman, Mohd Nasir Tamin, Ojo Kurdi马来西亚工程大学机械工程系81310 UTM, Skudai,Johor Bahru摘 要本文对一重型货车底盘做了应力分析。应力分析能够确定零件的最大受力点，是分析零部件疲劳研究和寿命预测的重要手段。前人已有用商用有限元软件 ABAQUS 软件对底盘模型进行分析的。本次研究的底盘长 12.35 米，宽 2.45 米，材料是 ASTM 低合金钢 710（3 级），屈服极限 552MPa，抗拉强度 620MPa。分析结果显示，最大应力点出现在底盘与螺栓连接的空缺处，最大应力为 386.9MPa，底盘的疲劳破坏将会从最大应力点开始向车架各部位蔓延。关键字：应力分析；疲劳寿命预测；货车底盘1 简介在马来西亚，很多货车的车架寿命都有 20 多年，20 多年架就会有使用安全的问题。因此，为了确保底盘在工作期间的安全性能，就有必要对底盘作疲劳研究和寿命预测。利用有限元法作应力分析能够确定受最大应力的关键点，这个关键点是导致底盘疲劳损伤的因素之一。应力的大小能够预测底盘的寿命，所以可以根据应力分析的结果精确地预测底盘的寿命，应力分析越精确，底盘寿命预测的越合理。本文是用商用有限元软件 ABAQUS 软件完成底盘应力分析的。汽车工业（汽车总成及各部件）在马来西亚的工业中占据非常重要的地位。随着东盟自由贸易区的贸易自由化发展，当地的汽车制造商和供应商应该顺应汽车及其零部件的世界级标准要求，比如噪声和振动就有相应的标准。马来西亚的汽车工业主要是依赖于国外的技术，而底盘是实现汽车轻量化的关键结构，所以底盘大多从国外进口。为了改变这种趋势，有必要建立发展马来西亚自己的底盘设计产业，这是对底盘进行研究的目标。底盘车架是汽车的装配基体和承载基体，支承着汽车的各个总成及零部件，如车轴，悬架系统，传动系，驾驶室及拖挂部件等，并将它们整合成一部完整的汽车。货车的底盘经常受到静载荷，动载荷以及周期性载荷。静载荷主要是车厢质量、货物及桂林电子科技大学第 2 页 共 9 页乘客，底盘的动载荷是由于货车的运动产生的，而发动机的振动和路面的不平整将会产生周期载荷。现有的底盘设计通常是基于静载荷的分析，设计的重点是底盘的强度结构设计，以支承施加在底盘上的载荷。然而货车底盘的受力复杂，包括静载、动载和疲劳破坏方面。据估计，85%到 90%的货车底盘的结构破坏是由疲劳破坏引起的 1。因此，货车底盘的动态和疲劳分析是很重要的。为了获得底盘的动态和疲劳工况的情况，就要确定各个零部件，如发动机、悬架、变速器等的支承点，并对其优化。许多研究人员都曾研究过货车底盘。Karaoglu and Kuralay 曾用铆接的连接方式对底盘所有限元应力分析 2。研究数据表明，局部增大纵梁的厚度可以减小边梁的应力，如果不能增大变量的厚度，增加接触面的面积也可以减小应力。Fermer et al 用高级疲劳分析软件 MSC/Fatigue 软件对沃尔沃双燃料车 S10 做了疲劳寿命分析 3，Conle and Chu 对复杂的底盘结构的疲劳分析和局部的应力应变分布做了研究 4，Ferreira et al 研究了汽车零部件耐久性的结构优化问题 5，Fermr and Svensson 研究了工业上焊接的汽车结构的铁基寿命预测问题 6。Filho et.al.考虑到小规模生产的经济可行性，结合适当的动载荷和结构特性对一越野车底盘做了设计分析和优化设计 7。研究表明，增大底盘的抗扭刚度，维持车架重心位置不变可以用来优化越野车结构，这样，底盘车架结构的总质量得到优化，结构也跟简单，生产成本也少了。Cosme et al 利用计算机辅助设计和工程软件代码Pro/E,ADAMS and ANSYS 模拟了改变设计对货车车架的影响 8。Chiewanichakorn et al 用试验得到的有限元模型，将已破坏的混凝土桥面替换为FRP 钢板，分析了桁架桥 9。结果数据表明，修复过后，桥的疲劳寿命是修复前钢筋混凝土桥面的两倍，在货车交通研究数据的基础上，桥面载荷及 ERP 钢板系统的应力范围在无限疲劳寿命范围中，即在其使用期间不会有桁架和地板系统的疲劳破坏现象。Ye 和 Moan 已经用有限元分析法分析了铝制框加强筋的车架静态和疲劳特性 10，改变车架切割形状和相应的焊接过程，同时得到足够的疲劳强度，这样就能够减小制造成本，并且解决连接问题。利用铁的疲劳可以确定可能产生疲劳裂纹的关键点，并能预测门铰链系统的寿命 11。本次研究中，对重型载货汽车施加静载荷，对其做应力分析，确定产生疲劳裂纹的危险点位置，以此作为该车架的疲劳寿命预测的备用数据。2 货车车架的有限元分析2.1 有限元法基本概念有限元分析法是一种计算机辅助技术，用来获得工程中边值问题的近似解。简言之，边值问题是一个数学问题，其中一个或是多个应变量必须要满足一个自变量范围已知的微分方程，还要满足特定的边界条件 12。桂林电子科技大学第 3 页 共 9 页有限元法的通俗解释是将一个结构离散成无数个单元（结构碎片），用简单的方法描述每个单元，然后用节点加各个单元重新连接起来，就像这些点是针脚或者点滴粘贴在一起形成各个单元（如图 1 所示）。这样就会产生一系列的同步代数方程。在分析应力时，这些方程式是节点的平衡方程，这样就会有数百甚至数千个这种方程，那么电脑的硬件要求就较高 13。图 1 二维轮齿的网格，所有的节点和单元都在纸平面内2.2 有限元法一般步骤有限元法可以分析一些物理问题，包括结构分析、流体分析、热传递和其他问题，分析这些问题有些通用的步骤，这些步骤通常包括一些商用有限元分析软件。主要有三大步骤，即前处理模块、求解模块和后处理模块。前处理模块要建立模型，这是必要的，如果发生了错误，就不会有完美的计算机有限元求解结果。这一步骤包括：定义问题的几何域，所需的单元类型，单元的材料属性，单元的几何性质（长度、面积等等） ，单元的连通性（网格划分） ，物理约束（边界条件）和加载。接下来就是求解，在这一步骤中，以矩阵方式列出的控制代数方程和未知的主变量是合成的，用计算结果回带求得其他派生变量，如反应力，单元应力和热流量。这一步骤要进行矩阵计算，数值积分，方程求解，这些都是由软件自动解决的。最后是后处理模块，对结果进行分析和评估。在这一部中，可以完成的操作包括按单元应力的大小分类，检查平衡，计算安全因素，绘制结构的变形形状，以动画的形式显示模型，以不同的颜色显示温度的分布。大型软件都会有一个前处理模块和后处理模块来完成分析，这两个模块都可以和其他的软件相同。前期处理和后期处理根据不同的项目会有各自的程序。2.3 货车的定义和分类货车是一种重型机动车辆，是用来承载货物的。货车的另一种定义是用来牵引的桂林电子科技大学第 4 页 共 9 页激动车辆。对货车的其他定义将根据货车的类型变化，例如自动倾卸卡车的货物可以作清空处理，车前端的平台末端就可以有空气作用被升起，此时载荷通过重力施加。房车或拖车有两种分类，一种是根据重量分类，由美国政府定的从 1 级到 8 级，如表 1 和表 2 所示；第二种是更为广泛的分类：轻型载重汽车；中型载重汽车和重型载重汽车。表 1 货车分类及等级重量等级 最小总质量额定值/磅最大总质量额定值/磅VIUS 分类 一般分类1 级 6000 轻型 轻型2 级 6001 10000 轻型 轻型3 级 10001 14000 中型 轻型4 级 14001 16000 中型 中型5 级 16001 19500 中型 中型6 级 19501 26000 轻型-重型 中型7 级 26001 33000 重型-重型 重型8 级 33001 重型-重型 重型表 2 制造商的货车分类分类 等级 总质量额定值 参考车型1 0-27kN（0-6000 磅）2 27-45 kN（ 6001-10000 磅）轻型3 45-62 kN（ 10001-14000 磅）卡车，救护车，运钞车4 62-71 kN（ 14001-16000 磅）5 71-87 kN（ 16001-19500 磅）6 87-116 kN（ 19501-26000 磅）中型7 116-147 kN（ 26001-33000 磅）市运货车，饮料运货车，拖吊车，校车重型8 147 kN 及以上（33000 磅及以上）卡车拖拉机，水泥搅拌车，自动倾卸卡车，消防车，城桂林电子科技大学第 5 页 共 9 页市公交客车注：总质量额定值：制造商指定的质量作为一辆车的最大装载质量（货车加货物） 。2.4 货车车架模型该模型如图 2 所示。模型长 12.35m，宽 2,45m，材料为 ASTM 低合金钢 710（3级） ，屈服极限 552MPa，抗拉强度 620MPa。车架的其他属性见表 3。表 3 货车车架的材料属性弹性模量 E(Pa) 密度（kg/m 3）泊松比 屈服极限（MPa）抗拉强度（MPa）207109 7800 0.3 550 620图 2 货车底盘模型2.5 加载货车模型承受来自车身和货物的静载荷，该车的最大装载质量为 36000kg，假设由最大载质量求得一个总的压力，将这个压力平均的分配到货物和底盘上表面的接触面上，具体的加载如图 3 所示，底盘上表面的压强为 67564.6N/m2 。图 3 静载荷（压强为 67564.6N/m2）桂林电子科技大学第 6 页 共 9 页2.6 边界条件本模型有 3 个边界条件。第一个施加在底盘前端，第二和第三个边界条件在底盘的后端，如图 4 所示。第一个边界条件是固定的（约束有轴的平移自由度，释放所有轴的旋转自由度） ，底盘与驾驶室的接触条件如图 5（a）所示。车架与车轴间由弹簧连接，将货物和底盘的重量传递到车轴上，所以第二个边界条件施加在底盘与弹簧上端连接的地方。第二个边界条件如图 5 所示，平移自由度只约束在轴 2 上，所有轴的旋转自由度都释放。第三个边界条件施加在底盘孔的内表面和螺栓的外表面的接触面处，在ABAQUS 软件中，这种接触是相互作用的，本文中的相互作用是面与面之间的摩擦作用。此时，螺栓所在的轴的平移自由度和旋转自由度都为零，称为固定约束。假定螺栓都是刚性元件，故螺栓选用杨氏模量很高的材料。图 4 模型的约束桂林电子科技大学第 7 页 共 9 页图 5 实物的约束注：a第一个边界条件，b、c第二个边界条件，d第三个边界条件3 分析结果及讨论在等效应力云图中，最大应力点在底盘开孔的地方，即与螺栓接触的地方，如图6 所示，最大应力为 386.9MPa，最大应力点在 86104 单元和 16045 节点上。底盘开孔处的内表面与非常坚硬的螺栓接触。第三个边界条件也是固定约束，因此会产生一个很大的应力。基于静态安全系数理论，取安全系数为 1.43，由安全系数公式得：安全系数=极限应力 /许用应力 （1）图 6 等效应力云图及最大应力点Vidosic 建议根据结构的载荷和材料选取一些安全系数，对于一些常用的材料，当载荷很容易确定时，安全系数可以取 1.5 到 2。基于分析结果，为了得到底盘精确的安全系数值，有必要减小最大应力值，因此对底盘结构进行修改以提高安全系数，尤其是在临界点区。桂林电子科技大学第 8 页 共 9 页底盘的位移和最大位移点如图 7 所示，最大位移为 4.995mm，位于底盘中部，最大的偏转在第一个边界条件和第二个边界条件的中部。为了验证分析结果，最大应力发生在第一个边界条件和第二个边界条件之间，这一部分可以近似的简化为一维的简支梁结构，在其中点施加集中力载荷，用施加在中点的集中力代替均匀分布在梁上的压力，这一力的大小等于压强的大小乘以受到压力的所有面的面积，求得结果与分析结果近似。计算求得的结果表明，这个简支梁的应变点在梁的中部，大小为：(2)图 7 应变分布云图及最大应变点位置模拟结果的最大应变值为 4.99mm，比数值分析计算结果大 11.2%。4 结论从数值分析可以看出，应力关键点出现在与螺栓连接的底盘孔处，最大应力值是很重要的，因为安全系数低于推荐值。由于疲劳破坏是从最大应力点开始的，可以断定，这一关键点是一些破坏的起源。因此，要注意减少这一点上的应力值，这是很重要的。分析得到的最大挠度的位置与受均布载荷的简支梁的最大变形位置一致。桂林电子科技大学第 9 页 共 9 页参考文献1 MSC. Fatigue, 2003. Encyclopedia. Los Angeles (CA, USA): MacNeal,Schwendler Corporation.2 Karaoglu, C. and Kuralay, N.S., 2000. Stress Analysis of a Truck Chassis with Riveted Joints, Elsevier Science Publishers B.V. Amsterdam, the Netherlands,Vol. 38, 1115-1130.3 Fermer, M., McInally, G. and Sandin, G., 1999. Fatigue Life Analysis of Volvo S80 Bi-Fuel using MSC/Fatigue, Worldwide MSC Automotive Conference, Germany.4 Conle, F.A. and Chu, C.C., 1997. Fatigue Analysis and the Local Stress-strain Approach in Complex Vehicular Structures, International Journal of Fatigue.5 Ferreira, W.G., Martins, F., Kameoka, S., Salloum, A.S. and Kaeya, J.T.,2003. Structural Optimization of Automotive Components Applied to Durability Problems, SAE Technical Papers. 6 Fermr, M. and Svensson, H., 2001. Industrial Experiences of FE-based Fatigue Life Predictions of Welded Automotive Structures, Fatigue & Fracture of Engineering Materials and Structures 24 (7), 2001, 489-500.7 Filho, R.R.P., Rezende, J.C.C., Leal, M. de F., Borges, J.A.F., 2003. Automotive Frame Optimization, 12th International Mobility8 Cosme, C., Ghasemi, A. and Gandevia, J., 1999. Application of Computer Aided Engineering in the Design of HeavyDuty Truck Frames, International Truck & Bus Meeting & Exposition, Detroit, Michigan, November 15 17.9 Chiewanichakorn, M., Aref, A.J., Allampalli, S., 2007. Dynamic and Fatigue Response of a Truss Bridge with Fiber Reinforced Polymer Deck, International Journal of Fatigue, 29, 14751489. 10 Ye, N. and Moan, T., 2007. Static and Fatigue Analysis of Three Types of Aluminium Box-Stiffener/Web Frame Connections, International Journal of Fatigue, 29, 14261433. 11 Bekah, S., 2004. Fatigue Life Prediction in a Door Hinge System Under Uni- Axial and Multiaxial Loading Condition, Master Thesis, Ryerson University, Toronto, Ontario, Canada.12 Hutton, David, V., 2004. Fundamental of Finite Element Analysis, Mc Graw Hill, New York.13 Cook, Robert, D., 1995. Finite Element Modeling for Stress Analysis, John Willey & Sons, Inc, New York. 14 Juvinall, R.C. and Marshek, K.M., 2006. Fundamental Machine Component Design, John Willey & Son, Inc., USA.15 Vidosic, J.P., 1957. Machine Design Project, Ronald Press, New York. Jurnal Mekanikal December 2008, No. 26, 76 - 85 76 STRESS ANALYSIS OF HEAVY DUTY TRUCK CHASSIS AS A PRELIMINARY DATA FOR ITS FATIGUE LIFE PREDICTION USING FEM Roslan Abd Rahman, Mohd Nasir Tamin, Ojo Kurdi*Faculty of Mechanical Engineering, Universiti Teknologi Malaysia, 81310 UTM, Skudai, Johor Bahru ABSTRACT This paper presents the stress analysis of heavy duty truck chassis. The stress analysis is important in fatigue study and life prediction of components to determine the critical point which has the highest stress. The analysis was done for a truck model by utilizing a commercial finite element packaged ABAQUS. The model has a length of 12.35 m and width of 2.45 m. The material of chassis is ASTM Low Alloy Steel A 710 C (Class 3) with 552 MPa of yield strength and 620 MPa of tensile strength. The result shows that the critical point of stress occurred at the opening of chassis which is in contact with the bolt. The stress magnitude of critical point is 386.9 MPa. This critical point is an initial to probable failure since fatigue failure started from the highest stress point. Keyword: Stress analysis, fatigue life prediction, truck chassis 1.0 INTRODUCTION The age of many truck chassis in Malaysia are of more than 20 years and there is always a question arising whether the chassis is still safe to use. Thus, fatigue study and life prediction on the chassis is necessary in order to verify the safety of this chassis during its operation. Stress analysis using Finite Element Method (FEM) can be used to locate the critical point which has the highest stress. This critical point is one of the factors that may cause the fatigue failure. The magnitude of the stress can be used to predict the life span of the truck chassis. The accuracy of prediction life of truck chassis is depending on the result of its stress analysis. The more accurate result of stress analysis the more valid the predicted life of object. In this study, the stress analysis is accomplished by the commercial finite element packaged ABAQUS. The automotive industry (vehicles and components) represents a strategic and important business sector in Malaysia. With the eventual trade liberalization of ASEAN Free Trade Area (AFTA), local automotive manufacturers and vendors shall require cars and components of world class standard. Noise and vibration are *Corresponding author: E-mail: ojokurdi Jurnal Mekanikal, December 2008 77 key elements in such standard. The automotive industry in Malaysia is much relying on foreign technology. Truck chassis, which is important structure of lightweight commercial vehicle, is mostly designed and imported from foreign country. In order to change this trend, it is necessary to develop and built Malaysian own chassis design. Study and research on truck chassis is thus required to achieve this goal. The chassis of trucks is the backbone of vehicles and integrates the main truck component systems such as the axles, suspension, power train, cab and trailer. The truck chassis is usually loaded by static, dynamic and also cyclic loading. Static loading comes from the weight of cabin, its content and passengers. The movement of truck affects a dynamic loading to the chassis. The vibration of engines and the roughness of road give a cyclic loading. The existing truck chassis design is normally designed based on static analysis. The emphasis of design is on the strength of structure to support the loading placed upon it. However, the truck chassis has been loaded by complex type of loads, including static, dynamics and fatigue aspects. It is estimated that fatigue is responsible for 85% to 90% of all structural failures 1. The knowledge of dynamic and fatigue behavior of truck chassis in such environment is thus important so that the mounting point of the components like engine, suspension, transmission and more can be determined and optimized. Many researchers carried out study on truck chassis. Karaoglu and Kuralay 2 investigated stress analysis of a truck chassis with riveted joints using FEM. Numerical results showed that stresses on the side member can be reduced by increasing the side member thickness locally. If the thickness change is not possible, increasing the connection plate length may be a good alternative. Fermer et al 3 investigated the fatigue life of Volvo S80 Bi-Fuel using MSC/Fatigue. Conle and Chu 4 did research about fatigue analysis and the local stress-strain approach in complex vehicular structures. Structural optimization of automotive components applied to durability problems has been investigated by Ferreira et al 5. Fermr and Svensson 6 studied on industrial experiences of FE-based fatigue life predictions of welded automotive structures. Filho et. al. 7 have investigated and optimized a chassis design for an off road vehicle with the appropriate dynamic and structural behavior, taking into account the aspects relative to the economical viability of an initial small scale production. The design of an off-road vehicle chassis is optimized by increasing the torsional stiffness, maintenance of center of gravity, total weight of structure and simpler geometry for reduction of production cost. The integration of computer aided design and engineering software codes (Pro/Engineer, ADAMS, and ANSYS) to simulate the effect of design changes to the truck frame has been studied by Cosme et al 8. Chiewanichakorn et al 9 investigated the behavior of a truss bridge, where an FRP deck replaced an old deteriorated concrete deck, using experimentally validated finite element (FE) models. Numerical results show that the fatigue life of the bridge after rehabilitation would be doubled compared to pre-rehabilitated reinforced concrete deck system. Based on the estimated truck traffic that the bridge carries, stress ranges of the FRP deck system lie in an infinite fatigue life Jurnal Mekanikal, December 2008 78 regime, which implies that no fatigue failure of trusses and floor system would be expected anytime during its service life. Ye and Moan 10 have investigated the static and fatigue behavior of aluminium box-stiffener/web frame connections using Finite Element Analysis (FEA) to provide a connection solution that can reduce the fabrication costs by changing the cutting shapes on the web frame and correspondingly the weld process meanwhile sufficient fatigue strength can be achieved. FE based fatigue was used to locate the critical point of probable crack initiation and to predict the life in a door hinge system 11. In this study, stress analysis of heavy duty truck chassis loaded by static force will be investigated to determine the location of critical point of crack initiation as a preliminary data for fatigue life prediction of this truck chassis. 2.0 FINITE ELEMENT ANALYSIS OF TRUCK CHASSIS 2.1 Basic Concept of FEM The finite element method (FEM) is a computational technique used to obtain approximate solutions of boundary value problems in engineering. Simply stated, a boundary value problem is a mathematical problem in which one or more dependent variables must satisfy a differential equation everywhere within a known domain of independent variables and satisfy specific conditions on the boundary of the domain 12. An unsophisticated description of the FE method is that it involves cutting a structure into several elements (pieces of structure), describing the behavior of each element in a simple way, then reconnecting elements at nodes as if nodes were pins or drops of glue that hold elements together (Figure 1). This process results in a set of simultaneous algebraic equations. In stress analysis these equation are equilibrium equations of the nodes. There may be several hundred or several thousand such equations, which mean that computer implementation is mandatory 13. Figure 1: A coarse mesh, two-dimensional model of gear tooth. All nodes and elements lie in plane of the paper 13 Jurnal Mekanikal, December 2008 79 2.2 A General Procedure for FEA There are certain common steps in formulating a finite element analysis of a physical problem, whether structural, fluid flow, heat transfer and some others problem. These steps are usually embodied in commercial finite element software packages. There are three main steps, namely: preprocessing, solution and postprocessing. The preprocessing (model definition) step is critical. A perfectly computed finite element solution is of absolutely no value if it corresponds to the wrong problem. This step includes: define the geometric domain of the problem, the element type(s) to be used, the material properties of the elements, the geometric properties of the elements (length, area, and the like), the element connectivity (mesh the model), the physical constraints (boundary conditions) and the loadings 12. The next step is solution, in this step the governing algebraic equations in matrix form and computes the unknown values of the primary field variable(s) are assembled. The computed results are then used by back substitution to determine additional, derived variables, such as reaction forces, element stresses and heat flow. Actually the features in this step such as matrix manipulation, numerical integration and equation solving are carried out automatically by commercial software 13. The final step is postprocessing, the analysis and evaluation of the result is conducted in this step. Examples of operations that can be accomplished include sort element stresses in order of magnitude, check equilibrium, calculate factors of safety, plot deformed structural shape, animate dynamic model behavior and produce color-coded temperature plots. The large software has a preprocessor and postprocessor to accompany the analysis portion and the both processor can communicate with the other large programs. Specific procedures of pre and post are different dependent upon the program 12. 2.3 Truck definition and classification Generally, truck is any of various heavy motor vehicles designed for carrying or pulling loads. Other definition of the truck is an automotive vehicle suitable for hauling. Some other definition are varied depending on the type of truck, such as Dump Truck is a truck whose contents can be emptied without handling; the front end of the platform can be pneumatically raised so that the load is discharged by gravity. There are two classifications most applicable to Recreational Vehicle tow trucks. The first one is the weight classes, as defined by the US government, ranging from Class 1 to Class 8 as listed in Table 1 and Table 2. The second is classified into a broader category: Light Duty Truck Medium Duty Truck Heavy Duty Truck Jurnal Mekanikal, December 2008 80 Table 1: Classification and classes of truck Weight Class Minimum GVWR (lbs) Maximum GVWR (lbs) VIUS * Category Common Category Class 1 6,000 Light-duty Light Duty Class 2 6,001 10,000 Light-duty Light Duty Class 3 10,001 14,000 Medium-duty Light Duty Class 4 14,001 16,000 Medium-duty Medium Duty Class 5 16,001 19,500 Medium-duty Medium Duty Class 6 19,501 26,000 Light-heavy Medium Duty Class 7 26,001 33,000 Heavy-heavy Heavy Duty Class 8 33,001 Heavy-heavy Heavy Duty Table 2: Vehicle manufacturer truck classification Category Class GVWR1Representative Vehicles 1 0 - 27 kN (0 - 6,000 lbs.) 2 27 - 45 kN (6,001 - 10,000 lbs.) Light 3 45 - 62 kN (10,001 - 14,000 lbs.)pickup trucks, ambulances, parcel delivery 4 62 - 71 kN (14,001 - 16,000 lbs.)5 71 - 87 kN (16,001 - 19,500 lbs.)6 87 - 116 kN (19,501 - 26,000 lbs.)Medium 7 116 - 147 kN (26,001 to 33,000 lbs.)city cargo van, beverage delivery truck, wrecker, school bus Heavy 8 147 kN and over (33,000 lbs. and over)truck tractor, concrete mixer, dump truck, fire truck, city transit bus Notes: Gross Vehicle Weight Rating (GVWR): weight specifiedby manufacturer as the maximum loaded weight (truck pluscargo) of a single vehicle. 2.4 Model of Truck Chassis The model is depicted in Figure 2. The model has length of 12.35 m and width of 2.45 m. The material of chassis is ASTM Low Alloy Steel A 710 C (Class 3) with 552 MPa of yield strength and 620 MPa of tensile strength. The other properties of chassis material are tabulated in Table 3. Jurnal Mekanikal, December 2008 81 Table 3: Properties of truck chassis material 14 Modulus Elasticity E (Pa) Density (kg/m3)Poisson Ratio Yield Strength (MPa) Tensile Strength (MPa) 207 x 1097800 0.3 550 620 Figure 2: Model of truck chassis 2.5 Loading The truck chassis model is loaded by static forces from the truck body and cargo. For this model, the maximum loaded weight of truck plus cargo is 36.000 kg. The load is assumed as a uniform pressure obtained from the maximum loaded weight divided by the total contact area between cargo and upper surface of chassis. Detail loading of model is shown in Figure 3. The magnitude of pressure on the upper side of chassis is calculated as 67564.6 N/m2. Figure 3: Static load (pressure = 67564.5 N/m2) 2.6 Boundary Conditions There are 3 boundary conditions (BC) of model; the first BC is applied in front of the chassis, the second and the third BC are applied in rear of chassis. They are shown in Figure 4. The type of BC 1 is pinned (the displacement is not allowed in all axes and the rotation is allowed in all axes) that represent the contact condition Jurnal Mekanikal, December 2008 82 between chassis and cab of truck as shown in Figure 5(a). The BC 2 represents the contact between chassis and upper side of spring that transfer loaded weight of cargo and chassis to axle. The contact condition of BC 2 in the object is shown in Figure 5. In the BC 2, the displacement only occurred in axis 2 and the rotation respect to all axes is zero. In the position where the BC 3 applied, there is a contact between inside surface of opening chassis and outside surface of bolt. In ABAQUS, this contact is called interaction. In this case, the type of the interaction is frictionless surface to surface contact. In the BC3, the displacement and the rotation is zero in all axes on all of bolts body. This condition is called fixed constrain. The bolt in BC 3 was assumed perfectly rigid. This assumption was realized by choosing a very high Youngs Modulus value of the bolt properties. Figure 4: Boundary conditions representation in the model Figure 5: Boundary conditions representation in the object, 4(a) BC 1, 4(b) and 4(c) BC 2, 4(d) BC 3 Jurnal Mekanikal, December 2008 83 3.0 RESULTS AND DISCUSSION The location of maximum Von Misses stress is at opening of chassis which is contacted with bolt as shown in Figure 6. The stress magnitude of critical point is 386.9 MPa. This critical point is located at element 86104 and node 16045. The internal surface of opening of chassis was contacted with the very stiff bolt. The BC 3 is also a fixed constraint, thus it cause a high stress on it. Based on static safety factor theory, the magnitude of safety factor for this structure is 1.43. The formula of Safety Factor (SF) is defined by 11: SF = significant strength of material corresponding significant stress from normal load (1)Figure 6: Von Misses stress distribution and critical point location Vidosic 15 recommends some value of safety factor for various condition of loading and material of structures. He recommends the value of 1.5 to 2 for well known materials under reasonably environmental condition, subjected to loads and stresses that can be determined readily. Based on this result, it is necessary to reduce the stress magnitude of critical point in order to get the satisfy SF value of truck chassis. The truck chassis can be modified to increase the value of SF especially at critical point area. The displacement of chassis and location of maximum displacement is shown in Figure 7. The magnitude of maximum displacement is 4.995 mm and occurs at middle of chassis. Maximum deflection is occurred at the middle of BC 1 and BC 2. For validation purpose, the region between BC 1 and BC 2 of chassis where the highest stress occurred is approximated by one dimensional simple beam loaded by concentrated force at mid point. The uniform distributed pressure on this region is replaced by a single concentrated force at mid point. The magnitude of the single force is obtained by multiplying the magnitude of pressure with the total area where the pressure is applied. The result agrees well with this a