doi:10.16576/j.cnki.1007-4414.2015.05.007320t矿用车行走机构非线性有限元分析*王摇亭,马摇力,陈东发,朱祝英(武汉理工大学汽车工程学院,湖北武汉摇430070)摘摇要:对行走机构进行了整体的受力分析,利用HyperMash软件采用非线性有限元法对矿用车行走机构的强度和刚度进行分析校核,根据分析结果发现立柱部分位置应力偏大。在对强度和刚度不满足的部分进行改进设计后,对改进的结构进行校核,实践结果说明改进是有效的。关键词:矿用车;行走机构;非线性有限元中图分类号:U463.33摇摇摇摇摇摇文献标志码:A摇摇摇摇摇摇文章编号:1007-4414(2015)05-0018-004NonlinearFiniteElementAnalysisofRunningMechanismfor320TonMineTruckWANG摇Ting,MA摇Li,CHENDong-fa,ZHUZhu-ying(SchoolofAutomotiveEngineering,WuhanUniversityofTechnology,WuhanHubei摇430070,China)Abstract:Inthispaper,theoverallstressconditionofrunningmechanismwasanalyzed.Thenonlinearfiniteelementmethodwasusedforminecarrunninggear;thestiffnessandstrengthwereanalyzedbyusingtheHyperMashsoftware,andaccordingtotheanalysisresult,itwasfoundthatthestressofthepillarpositionistoolarge.Animproveddesignwasconductedtothestructurewhichdoesnotmeettherequirements,andthenthestiffnessandstrengthwerechecked,thepracticeresultsshowedthatimprovementswereeffective.Keywords:minetruck;runningmechanism;nonlinearfiniteelement0摇引摇言所研究的载重量320t矿用车为五轴线电动轮全轮驱动车辆,它的行走系统采用液压模块组合挂车行走机构的形式。车辆长约16m,宽约8m,载重量和横向尺寸均比较大,且每个车轮都有自己独立的驱动系统,因此行走机构的车轮采用宽距双胎结构。与传统车辆的行走机构相比,这种宽距双胎结构受力比较复杂,主体部件设计分析属于状态非线性问题。目前关于相关车辆的设计分析主要集中在液压模块组合挂车方面[1],涉及到车架轻量化和转向机构优化等问题[2-3],针对载重量200t以上的超重型矿用车有人进行了异形货厢的参数化和拓扑优化设计。但是针对多轴线电动轮全轮驱动矿用车的行走机构设计方面,目前国内还未见到相关报道,企业在产品设计方面遇到了困难,初始设计产品的行走机构遇到实际损坏。因此笔者研究了宽距双胎电动轮行走机构设计问题,对结构进行了受力分析,对典型结构采用状态非线性方法进行了整体仿真计算,找出结构损坏的原因,并进行了结构的改进设计。1摇行走机构的结构特点及受力分析1.1摇宽距双胎行走机构组成及特点五轴线电动轮全轮驱动矿用车共有10个车轮组,每侧各五个。图1为每个车轮组的行走机构的示意图,主要由回转盘、立柱、摆臂机构、驱动桥壳、轮胎和油气悬架和组成。回转盘通过轴承与车架回转连接,推动回转盘可实现车轮组的转向;立柱和回转盘焊接为一体,主要起支撑和传力的作用;摆臂左端通过销轴和立柱连接,通过绕销轴的转动来保证车轮的跳动和车身升降的运动关系[4-5];驱动桥壳内装有驱动电机、减速器以及制动器等部件,使得驱动桥轴向尺寸较大,两端轮胎相距较远,形成宽距双胎结构。驱动桥壳通过其纵向通孔与摆臂右侧的摆臂轴铰接在一起,使驱动桥壳和车轮可以绕摆臂轴横向偏转,保证宽距双胎结构中的两个轮胎都能有效接地,防止单胎过载爆胎;油气悬架上端与回转盘铰接,下端与摆臂中部铰接,实现减振和车身的升降功能。图1摇行走机构二维简图摇摇由此可见,行走机构主体部件均是接触连接关系,而且接触面比较大,呈现比较强的状态非线性特·81·研究与分析摇摇摇摇摇摇摇摇摇摇摇摇摇摇摇摇摇摇2015年第5期(第28卷,总第139期)·机械研究与应用·*收稿日期:2015-08-14作者简介:王摇亭(1991-),女,山东德州人,在读硕士,研究方向:汽车CAD/CAE。征。另一方面,由于宽距双胎结构受力比较复杂,因此针对单个部件的常规设计方法容易出现设计问题,需要面向整体结构利用非线性方法进行设计分析。1.2摇宽距双胎行走机构受力分析以车轮中心O为原点建立坐标系,车辆纵向为x方向,横向为y方向,垂直方向为z方向。在没有横向力作用的情况下,车轮组的两个轮胎受到地面的支撑反力与纵向力相等,其合力作用在轮胎与地面接触点P点,分别为Fz和Fx,如图2(a)所示,A、B、C、D、E点为各部件之间的连接点。制动情况下,制动力方向与图2(a)中Fx方向相同;驱动情况下,驱动力方向与图2(a)中Fx方向相反,这两种工况下受力情况与常规车辆的受力情况基本相同,车架与回转盘之间的相互作用力与这两个力互相平衡。显然有:Fx=Fz渍(1)式中:在制动情况下,渍为制动力系数,在驱动情况下,渍为附着系数。有侧滑力的情况下,宽距双胎结构与常规车辆的受力明显的不同。由于驱动桥壳3可以绕着摆臂4的纵向轴线转动,而且两轮胎的距离比较宽,约1.5m,因此两个轮胎上的垂直反力和横向力不相等,驱动或制动时的纵向力也不相同,如图2(b)所示,Fz1和Fz2分别为左右两轮胎上的支撑力,Fy1和Fy2左右两轮胎上的侧向力。设两轮胎间距离为L,P点和C点坐标分别为(xp,yp,zp)和(xc,yc,zc),则可得力和力矩的平衡方程:(Fy1+Fy2)(zc-zp)=(Fz2-Fz1L2)Fy1+Fy2=Fy;Fz1+Fz2=FzFy1=Fz1滋;Fy2=Fz2滋Fx1=Fz1渍;Fx2=Fz2渍式中:Fx1和Fx2分别为两轮胎上的纵向力;滋为侧向摩擦系数;Fy为侧滑力合力;fz为支撑力合力。Fy和Fz可由车辆的行驶状况和载荷状况求得。图2摇行走机构示意图1.油气悬架摇2.车轮摇3.驱动桥壳摇4.摆臂摇5.立柱摇摇转向时,回转盘上获得的转向力矩与地面车轮之间产生的转向阻力矩互相平衡。由此可见宽距双胎行走系统各部件受力情况要比常规车辆复杂一些。2摇有限元模型的建立建立的典型320t矿用车行走机构的几何模型如图3。由于油气悬架力学上是典型的二力弹性元件,且计算分析中不关注油缸和活塞的强度问题,只需考虑油气弹簧的刚度贡献,因此图3中没有必要给出油气悬架的三维模型,计算时采用弹性单元进行等效替换即可。这样可以减少计算规模,减少计算时间。主体结构采用三维实体单元划分网格。根据前面的结构特点分析可知,摆臂与立柱、摆臂与油气弹簧下端、油气悬架上端与回转盘的销轴连接以及摆臂与驱动桥壳通孔处的贯通铰接,理论上都属于有限元状态非线性问题。但是若全部连接关系都按非线性问题处理,接触面较多,计算规模大,分析时间较长,因此需要根据结构特点进行局部的线性化处理。摆臂与通孔接触长度大,受力复杂,初始设计时没有常规的设计可循,而且实际结构正是该部分出现设计不当早期失效的现象,因此该部分必须进行状态非线性问题处理。考虑到销轴连接的设计与计算都存在比较完善的常规计算方法,实际结构在这些部位没有存在设计不当出现损坏现象,故在满足静力学等效的前提下,将它们简化为线性问题,以减小计算规模和计算时间。根据圣维南原理,这样的简化不会影响稍远结构的应力水平和应力分布。定义摆臂右端的摆臂轴的外表面和驱动桥壳通孔的内表面为接触对,选择面接触单元。接触区域的网格大小取为10mm。大量的试算表明,由于通孔的接触长度长达0.5m,这样的网格密度能够保证接触区域具有足够的计算精度。利用HyperMesh里Spring单元模拟油气悬架。利用HyperMesh中的HyperBeam单元模拟销轴,油气弹簧与摆臂的销轴连接以及摆臂与立柱销轴连接的线性化处理示意图如图4所示。图3摇行走机构的几何模型摇图4摇销轴连接模拟1.摆臂摇2.立柱摇3.驱动桥克摇摇图4中销轴采用6个自由度的梁单元模拟,梁单元相关节点与摆臂和立柱内表面关联。这样既保证了各部件的连接关系和销轴的刚度贡献,又不会限制各部件绕销轴的转动趋势。只是注意要施加附加的·91··机械研究与应用·2015年第5期(第28卷,总第139期)摇摇摇摇摇摇摇摇摇摇摇摇摇摇摇摇摇摇研究与分析约束条件来限制梁单元即销轴绕自身转动的自由度。位移约束施加在回转盘与车架的连接处,地面载荷按各工况的实际值静力等效到车轮轮心位置,然后通过节点关联的方式等效传递到轮毂与车轮的安装连接端面。最后计算模型的规模为单元371850个,约30万个自由度。3摇有限元计算及结果分析结构材料为高强钢,其弹性模量取2.1伊1011N/m2,泊松比为0.3,材料的屈服强度极限约为785MPa,根据项目组的疲劳试验结果[6],材料的疲劳极限约为500MPa。计算时,悬架刚度在静平衡位置约为400伊103N/m2。3.1摇工况选择及载荷的计算根据320t矿用车行走机构的受力特点,行走机构在极限工况下更易造成部件的损坏,在这里选取三个典型工况对行走机构进行整体有限元分析。(1)满载静止或匀速工况摇模拟车辆静止或匀速行驶,考虑满载装载质量,行走机构之上的结构自重(行走机构自重不计)以及动载系数,计算得地面对车轮的载荷FZ大小为509600N。此时Fz1和Fz2均为FZ/2。(2)满载制动工况摇车轮的外径约为2.48m。取制动力系数渍为0.55,可得,行走机构在极限制动工况下,受到的力矩大小约为348248N·m。由于这种300t以上的重型车辆不允许在车速比较高的情况下出现自动抱死现象,因此按抱死计算制动力矩载荷偏大,但对设计而言,这样做是偏安全的。(3)满载转向工况摇满载转向工况时,转向阻力矩Mn按下式计算:Mn=渍3G3P式中:P为轮胎的压强,取0.7MPa,可求得Mn大小约为79714N·m。按偏安全设计,最大横向力取极限情况计算约为280280N。3.2摇计算结果分析在满载制动和满载转向工况下,整体结构应力水平较大,应力云图分别如图5和图6所示,最大应力分别为652MPa和639MPa,发生在图中1和2所指的立柱支撑板位置。在满载静止或匀速工况下整体应力水平不高,约为420MPa;在工况1和工况3下,摆臂的最大应力均在280MPa以下,工况2下的摆臂应力云图如图7所示,最大应力约为360MPa,发生在图中3所示位置,该位置出现了明显的应力集中;三种工况下的驱动桥壳的应力水平都在270MPa以下,有很大的轻量化空间,其设计另行描述。工况1下驱动桥壳的应力云图如图8所示,最大应力约为70MPa,发生在图中4所示位置。立柱、摆臂、驱动桥壳在三种工况下的最大应力的具体数据如表1所列。图5摇工况2下的整体应摇图6摇工况3下的整体应力云图力云图图7摇工况2下的摆臂应摇图8摇工况1下的驱动桥壳力云图应力云图表1摇三种工况下的最大应力/MPa工况立柱最大应力摆臂最大应力驱动桥壳最大应力满载静止或匀速42417170满载制动652364268满载转向639274208计算结果表明,摆臂和驱动桥壳的应力水平不高,符合要求。但是立柱局部位置应力水平达到了650MPa,应力水平比较高。高应力位置与早期损坏一致,需要改进设计。4摇改进设计4.1摇改进措施综上所述,摆臂和驱动桥壳应力水平比较低,暂不改进。这里的改进措施主要针对立柱结构,以解决存在的结构早期损坏的问题。(1)增加支撑板的厚度摇支撑板的厚度小,使得支撑板的整体应力水平比立柱的其它