120吨重载冶金框架车悬挂推力系统研究.pdf
120 吨重载冶金框架车悬挂推力系统研究 ① 李 力1,2, 赵 楠1, 张卫东3 (1.中南大学 机电工程学院,湖南 长沙 410083; 2.深海矿产资源开发利用技术国家重点实验室,湖南 长沙 410012; 3.长沙凯瑞冶金机械有限公司,湖 南 长沙 410600) 摘 要 以 120 吨冶金重载框架车悬挂推力系统为研究对象,基于多刚体理论和有限元理论,建立了 120 吨重载框架车悬挂推力系 统的三维动力学模型、参数化模型及有限元模型,开展了悬挂推力系统在启动和制动工况下的动力学仿真、参数化优化及有限元分 析,获得了悬挂推力系统优化参数模型和优化前后悬挂推力系统应力云图。 仿真结果表明,优化后悬挂推力系统在路面起伏 50 mm 时其 V 型推力系统应力小于优化前,但二者均超过屈服极限,无法满足强度要求。 在此基础上,提出了悬挂推力系统梯形推力 机构的改进方案,并进行了有限元仿真与分析。 仿真结果表明,改进后悬挂推力系统应力大大降低,满足强度要求,为 120 吨重载 框架车悬挂推力系统研发及升级提供了理论基础。 关键词 冶金框架车; 悬挂系统; V 型推力机构; 参数化; 有限元; ADAMS 软件; ANSYS 软件 中图分类号 TF303;TP391.9文献标识码 Adoi10.3969/ j.issn.0253-6099.2014.06.034 文章编号 0253-6099(2014)06-0139-05 Suspension System with a Thrust Rod for 120t Heavy⁃Duty U⁃frame Truck Used in Metallurgy Industry LI Li1,2, ZHAO Nan1, ZHANG Wei⁃dong3 (1.School of Mechanical and Electrical Engineering, Central South University, Changsha 410083, Hunan, China; 2.The State Key Laboratory of Deep Sea Mineral Resources Development and Utilization Technology, Changsha 410012, Hunan, China; 3.Changsha Kairui Metallurgical Machinery Co Ltd, Changsha 410600, Hunan, China) Abstract Based on the multi⁃rigid⁃body theory and finite⁃element theory, three⁃dimensional kineticmodel, parameterization model and finite⁃element model were established for the suspension system with a suspension system with a thrust rod for 120 t heavy⁃duty u⁃frame vehicle used for transportation in metallurgy industry. After the dynamics simulation, parameter optimization and finite⁃element analysis for the suspension system with a thrust rod in the starting and braking conditions, the system′s optimized model, as well as stress nephograms before and after the optimization were obtained. The simulation shows that, while driving on the road with 50 mm hump height, the stress of the V⁃shaped thrust rod structure is lower than that before optimization. However, both of them cannot meet the strength requirement due to excessing the yield point. On this basis, an improved suspension system with a thrusting machine of trapezoidal shape was proposed, for which finite⁃element simulation and analysis were carried out correspondingly. Results show that the stress of the improved suspension system has decreased remarkably, meeting the requirement in strength. Thereby, this analysis can provide theoretical basis for the following R&D and upgrading of the suspension system. Key words u⁃frame vehicle used in metallurgy industry; suspension system; V⁃shaped thrust rod structure; parameterization; finite⁃element analysis; ADAMS; ANSYS 重载冶金框架车在冶金运输中的地位极为重要, 主要承担钢、铝、锌等冶金厂的金属液体包、锭、卷、板 等超重物的运输工作,其负重达上百吨。 目前,120 吨 重载冶金框架车是冶金运输的主力车型,它由前后驱 动桥、液压悬挂推力系统和调平系统、7 排实心橡胶轮 胎及车架等组成,如图 1 所示。 由于运输过程中的超 大负荷、不平地面及转弯等恶劣因素,120 吨重载冶金 框架车将承受多种极大动载荷的耦合作用,车体主要 零部件经常出现变形、断裂等破坏现象,直接威胁到车 体的正常运行,其中最容易出现损坏的关键部件之一 为悬挂推力系统,其经常出现连接杆断裂、铰接头磨 损、机座裂缝等现象,如图 2 所示。 ①收稿日期 2014-06-20 基金项目 国家重点基础研究发展计划(973 计划)资助项目(2014CB046305) 作者简介 李 力(1963-),女,湖南益阳人,教授,研究方向为机械电子工程、海洋采矿技术。 第 34 卷第 6 期 2014 年 12 月 矿矿 冶冶 工工 程程 MINING AND METALLURGICAL ENGINEERING Vol.34 №6 December 2014 图 1 120 吨重载冶金框架车 图 2 悬挂推力系统现场照片 悬挂推力系统在重载框架车中将前、后车桥分别 与两侧车架连接,主要用于传递与承受行驶过程中地 面的纵向力、横向力及力矩作用[1],以保持车身的稳 定性,缓冲和减轻由路面不平引起的强大冲击和振动, 减缓轮胎磨损。 基于悬挂推力系统在重载框架车中的 重要作用,有必要对其开展深入的研究,以提高悬挂推 力系统的可靠性与安全性,从而为重载框架车整车安 全性提供技术保证,对高质量和高水平的 120 吨重载 框架车产品研发与升级具有重要意义。 在悬挂推力系 统研究方面,已有诸多研究成果[2-7]。 本文以 120 吨 重载冶金框架车悬挂推力系统为研究对象,以 ADAMS 软件为平台,建立悬挂推力系统三维多刚体动力学模 型和参数化模型,开展悬挂推力系统仿真与优化,利用 ANSYS 软件,开展悬挂推力系统有限元仿真与分析, 提出悬挂推力系统改进方案,为 120 吨重载冶金框架 车悬挂推力系统设计提供理论支持。 1 悬挂推力系统参数化建模与优化 1.1 多刚体理论 ADAMS 是根据多体动力学理论开发的复杂机械 系统动力学仿真软件。 根据机械系统刚体质心坐标 系,运用 Lagrange(拉格朗日)算子,建立完整的系统动 力学方程[8-9] d dt ∂T ∂̇ q T - ∂T ∂q T + φqTρ + θ̇ q Tμ = Q 完整约束方程 φ(q,t) = 0 非完整约束方程 θ(q,̇ q,t) = 0 式中 T 为系统动能;q 为系统广义坐标矩阵;Q 为广义 力矩阵;ρ 为对应于完整约束的拉氏乘子矩阵;μ 为对 应于非完整约束的拉氏乘子矩阵。 1.2 悬挂推力系统动力学建模与仿真 1.2.1 悬挂推力系统动力学建模 120 吨重载冶金框 架车具有前、后悬挂推力系统两套。 将 120 吨重载冶 金框架车 Pro/ E 模型以 X_T 格式导入到 ADAMS 仿真 软件中,根据各部件间约束关系,合并零件,添加约束, 建立 120 吨重载框架车三维动力学模型,其中,悬挂推 力系统动力学模型如图 3 所示,V 型推力机构动力学 模型如图 4 所示,X 向为车体横向,Y 向为车体垂向,Z 向为车体纵向。 图 3 悬挂推力系统三维多刚体动力学模型 1 V 型推力机构;2 悬挂油缸;3 驱动桥;4 动臂; 5 车架 图 4 V 型推力机构三维多刚体动力学模型 由图 3 可知,悬挂油缸的一端以回转副铰接于车 架,另一端以回转副铰接于动臂,动臂的一端以回转副 铰接于车架,另一端以回转副铰接于驱动桥,驱动桥轴 041矿 冶 工 程第 34 卷 上安装车轮,V 型推力机构的双头端分别与车架以回 转副铰接,其单头与驱动桥以球副铰接。 由图 4 可知, 悬挂推力系统的关键部件 V 型推力机构由二连接杆 构成 V 字型,分别通过 3 个运动副与左、右车架和驱 动桥联接,其双头端 A 点、B 点为以车体横向 X 轴为 轴心的回转副,其单头端 C 点为球副。 悬挂推力系统工作原理为,驱动桥输送动力至两 端车轮而使车体运动,悬挂油缸通过动臂连接驱动桥 和车架,遇到不平路面时可上下浮动,缓和来自地面的 冲击力和车体振动,V 型推力机构通过连接左右车架 及驱动桥来承受车体前进和转向等所产生的力和 力矩。 1.2.2 悬挂推力系统动力学仿真 仿真工况车体 120 吨负载,直线轨迹,启动加速设置为在 0~8 s 内车速达 到 20 km/ h,在 8~22 s 内,前进车速均为 20 km/ h,制动 减速设置为从 22~30 s 内车速从 20 km/ h 至 0。 利用 ADAMS 软件平台,将仿真工况输入 120 吨 重载框架车三维动力学模型,完成其在启动与制动过 程中的动力学仿真,获得其动力学参数及动画等[10], 从而得到悬挂推力系统的动力学特征,其中,V 型推力 机构的双头 A 点和 B 点与左右车架的回转副处受力 曲线见图 5。 由图 5 可知,V 型推力机构 A 点和 B 点 处受力在启动加速至 8 s 时达到最大值 35 kN,制动减 速至 27.5 s 时达到 31.6 kN。 由此可知,120 吨重载框 架车在启动加速与制动减速时,其悬挂推力系统中 V 型推力机构受力状况非常恶劣。 图 5 V 型推力机构 A 点和 B 点处受力曲线 1.3 悬挂推力系统参数化建模与优化 将悬挂推力系统的 V 型推力机构 A、B 两点位置 作为设计变量,分别设 A 点和 B 点 X 向坐标为变量 DV_1 和-DV_1,A 点和 B 点 Y 向坐标为变量 DV_2,A 点和 B 点 Z 向坐标为变量 DV_3,如表 1 所示。 根据 V 型推力机构 A 点和 B 点的结构尺寸,设变化范围为 X 向( - 100 mm,0),Y 向( - 600 mm,100 mm),Z 向 (-100 mm,100 mm),A 点和 B 点原始坐标分别为 (367 mm,-803 mm,6 282 mm) 和(-367 mm,-803 mm,6 282 mm),则设计变量取值范围为 267≤DV_1≤ 367,-1 403≤DV_2≤-703,6 182≤DV_3≤6 382,同 时,设目标函数[11]为推力杆与固定臂铰接点 A 点和 B 点受力取最小。 表 1 V 型推力系统设计变量参数 铰接点Loc_XLoc_YLoc_Z A(.MODEL_PT1201.DV_1)(DV_2)(DV_3) B(SIGN(.MODEL_PT1201.DV_1,-1))(DV_2)(DV_3) 利用 ADAMS 软件平台,建立悬挂推力机构系统 参数化模型,其设计变量满足上述设置范围,创建 A 点和 B 点受力测量传感器,开展悬挂推力机构系统优 化仿真,获得 A 点和 B 点的优化参数结果,见图 6,A 点和 B 点受力变化曲线见图 7。 图 6 V 型推力机构优化结果 图 7 A 点和 B 点的受力变化曲线 分析图 6 可知,A 点坐标为(367 mm,-1 253 mm, 6 289 mm),B 点坐标为(-367 mm,-1 253 mm,6 289 mm),由此可知,A 点和 B 点 X 向与原来坐标基本一 致,Z 向坐标比原来坐标上移仅 7 mm,Y 向坐标比原 141第 6 期李 力等 120 吨重载冶金框架车悬挂推力系统研究 来坐标下移 450 mm。 分析图 7 可知,A 点和 B 点的受 力值为 23.39 kN,比优化前减少 25.7%,从而获得了悬 挂推力机构系统优化方案。 由此可知,根据 C 点坐标 (0,-1 253 mm,5 689 mm),V 型推力机构处于水平状 态时,其受力状况为最小。 2 悬挂推力系统有限元分析 2.1 有限元建模 根据优化前与优化后悬挂推力系统方案,建立二者 的三维实体模型,并导入 ANSYS 软件,经过修补与缝合 缝隙[12],采用 Solid186 和 Solid187 单元开展模型网格划 分,分别获得悬挂推力系统优化前和优化后有限元模 型,如图 8 所示,二者分别为网格节点数 200 460 和单 元数114272 及网格节点数181322 和单元数103068。 已知悬挂推力系统所用材料为低合金高强度结构钢 Q345B,弹性模量为 2.06105MPa,屈服极限为 345 MPa,泊松比为 0.3,密度为 7.8510 kg/ m3。 图 8 悬挂推力系统有限元模型 2.2 有限元仿真与分析 采用 120 吨重载框架车整车在启动与制动工况过 程中动力学仿真数据,将其加载到悬挂推力系统有限 元模型的驱动桥上,施加固定约束于车架,设置重力, 考虑由于路面不平引起悬挂油缸位移 50 mm,利用 ANSYS 软件平台,开展优化前和优化后悬挂推力系统 有限元仿真,获得了优化前后悬挂推力系统应力云图, 见图 9。 图 9 优化前后悬挂推力系统应力云图 分析图 9 可知,优化前悬挂推力系统 V 型推力机 构二连接杆双头及单头附近受力较大,此部分应力值 均大于其屈服极限,A 点和 B 点附近的连接杆处最大 应力值达到 2 179 MPa。 优化后 V 型推力机构二连接 杆最大应力值为 365.83 MPa,整体应力值大大下降, 但仍超过其屈服极限,不满足强度要求。 由此可知,此 方案的悬挂推力系统连接杆和铰接头处将出现断裂与 破坏情况,与其实际破坏情况相符,且有必要进一步开 展悬挂推力系统 V 型推力机构改进,以达到消除其断 裂与破坏状况的目的。 3 悬挂推力系统改进研究 根据悬挂推力系统结构空间[13-15]和其动力学仿 真所表明 V 型推力机构二连接杆受力随其夹角减小 而减小的特征,提出 V 型推力机构改进方案,见图 10。 V 型推力机构改进方案为将 C 点处增加杆 DE,减小 二连接杆间的夹角,V 型结构变为梯形,其各运动副的 型式不变。 从图 10 可知,以车体纵向为对称面, DE 在 X 向为 540 mm;A 点和 B 点 X 向距离为 694 mm,Y 向和 Z 向位置不变。 图 10 悬挂推力系统改进方案 将悬挂推力系统改进方案输入 ANSYS 软件,采用 上述同样的边界条件和设置,建立悬挂推力系统有限 元模型,如图 11 所示,获得了悬挂推力系统和梯型推 力机构有等效应力云图,分别如图 12~13 所示。 分析 图 12 和图 13 可知,改进后悬挂推力系统所受应力小 于屈服极限,最大应力位于驱动桥与动臂交界处,约为 180 MPa,比优化后悬挂推力系统最大应力减小了50.8%; 图 11 悬挂推力系统改进方案有限元模型 241矿 冶 工 程第 34 卷 图 12 悬挂推力系统改进方案应力云图 图 13 梯型推力机构应力云图 梯型推力机构最大应力位于二连接杆,为 63.5 MPa, 比优化前 V 型推力机构最大应力减小了82.6%。 由此 可知,经过改进后悬挂推力系统应力大大降低,梯型推 力机构受力状况得到有效改善,满足设计要求。 4 结 论 1) 基于多刚体理论,利用 ADAMS 平台,建立了 120 吨重载框架车三维动力学模型和虚拟样机,开展 了其在启动加速和制动减速工况下的动力学仿真研 究,获得了其悬挂推力系统动力学特性;建立了悬挂推 力系统参数化模型,对其开展了优化仿真,获得了悬挂 推力系统的优化模型。 2) 利用 ANSYS 软件,建立了优化前和优化后悬 挂推力系统有限元模型,以其在启动加速和制动减速 工况下的动力学特性、地面不平度 50 mm 为边界条 件,开展有限元分析,获得了优化前和优化后悬挂推力 系统的应力云图,分析可知,优化后悬挂推力系统应力 值小于优化前,但二者均不满足强度要求。 3) 提出了将悬挂推力系统 V 型推力机构改为梯 形推力机构的改进方案,开展了其有限元建模与分析, 获得了悬挂推力系统的应力云图,分析表明,改进后悬 挂推力系统最大应力约为 180 MPa,梯型推力机构二 杆应力均有效改善,满足强度要求。 本研究为 120 吨 重载框架车悬挂推力系统研发及升级提供了理论 基础。 参考文献 [1] 张长江,侯文英. 重卡悬架推力杆优化设计分析[J]. 商用汽车, 2013(6)80-82. [2] Leite VJS, Peres PLD. Pole location control design of an active suspension system with uncertain parameters [ J]. Vehicle System Dynamics, 2005,43(8)561-568. [3] Jun K J, Park T W, Lee S H, et al. Prediction of fatigue life and estimation of its reliability on the parts of an air suspension system [J]. International Journal of Automotive Technology, 2008,9(6) 741-747. [4] 张红红,陈昆山,张 建. 空气悬架 V 型推力杆静动态特性分析 [J]. 拖拉机与农用运输车,2009,36(3)40-42. [5] 谢万春,李英昊. 汽车空气悬架系统 V 型推力杆总成试验方法探 讨[J]. 汽车科技,2006(5)50-52. [6] 曾晶晶,卜继玲,刘 羽. 某型汽车推力杆结构与疲劳性能分析 [J]. 客车技术与研究,2011(4)14-16. [7] 江 洪,成中书,张 丛. 空气悬架导向机构纵向推力杆的优化设 计[J]. 机械设计与制造,2010(6)1-3. [8] 李增刚. ADAMS 入门详解与实例[M]. 北京国防工业出版社, 2006. [9] Stronge W J. Generalized Impulse and Momentum Applied to Multi⁃ body Impact with Friction[J]. Mechanics of Structures and Machines, 2001,29(2)239-260. [10] JIANG Libiao, WANG Dengfeng, XIE Dong, et al. Buildup and simulation of a truck model with rear air suspension using leaf spring as guiding rod[J]. Journal of Harbin Institute of Technology,2006, 13(6)654-657. [11] 陈立平. 机械系统动力学分析及 ADAMS 应用教程[M]. 北京 清华大学出版社, 2008. [12] 张银生. PT200 液压平板车结构有限元分析与试验测试研究 [D]. 长沙中南大学机电工程学院, 2008. [13] 候永坤. 悬架系统 V 型杆的结构分析[J]. 客车技术与研究, 2007,29(4)20-22. [14] 史文库,柯 俊,王 倩,等. V 型推力杆的极限载荷分析与球较 结构优化[J]. 西安交通大学学报,2013,47(10)132-136. [15] Antonelliv, Marissen R. FE⁃modelling of the rubber mould behavior during press forming of thermoplastic composites[J]. International Journal of Material Forming,2010,3(4)679-682. 341第 6 期李 力等 120 吨重载冶金框架车悬挂推力系统研究