北京理工大学 机械与车辆学院, 北京 100081
收稿日期: 2015-06-10
基金项目: 国家部委基金资助项目(1030020220707).
作者简介: 王振峰(1987-), 男, 河南新乡人, 北京理工大学博士研究生;
顾亮(1957-)男, 山东淄博人, 北京理工大学教授, 博士生导师。
摘要: 为了研究扭杆弹簧扭转精度对其强扭过程的影响，利用弹塑性理论对扭杆弹簧强扭过程进行了理论探讨，获得扭杆弹簧在屈服与塑性强扭工况下的断面应力、扭矩及残余应力等变化的数学模型.应用ABAQUS有限元软件对不同型号、强扭工况的扭杆弹簧进行数值模拟，分析了扭转精度对扭杆弹簧强扭效果的影响.结果表明，在强扭工艺允许范围内，扭转误差对扭杆弹簧强扭效果影响不超过5%.试验结果不仅验证了有限元模型的正确性，也得出扭杆弹簧强扭误差在不超过±1°工况下的强扭效果可近似等同的结论，此研究为扭杆弹簧强扭工艺的完善提供理论依据.
关键词：扭杆弹簧弹塑性强扭工艺数值模拟强扭试验
Effect of Torsion Precision on Torsion Bar Spring Pre-setting Process
WANG Zhen-feng, DONG Ming-ming, ZHAO Kai, GU Liang
School of Mechanical and Vehicle Engineering, Beijing Institute of Technology, Beijing 100081, China
Corresponding author: GU Liang, E-mail: guliangbit@gmail.com
Abstract: In order to study the effect of torsion precision on pre-setting process, the elastic-plastic theory was used to discuss the pre-setting process of torsion bar spring. Then, a mathematical model changing with stress, torque and residual stress was obtained under yield pre-setting and plasticity pre-setting situation. The ABAQUS software was used to simulate torsion bar spring with different models and pre-setting situations, and the effect of pre-setting precision on torsion bar spring was analyzed. The results showed that the effect of torsion bar spring pre-setting process in the range of torsion error on pre-setting is less than 5%. The test data not only validates correctness of the finite element model, but shows that torsion bar spring with no more than±1° pre-setting error can be seen approximately equivalent. The results provide a theoretical reference for improving the process of torsion bar spring pre-setting.
Key Words: torsion bar springelastic-plasticpre-setting processnumerical simulationpre-setting test
扭杆悬挂作为履带车辆的缓冲和储能元件，具有结构可靠、质量轻、空间小、装配维修方便等特点而得到了普遍关注.同时，扭杆弹簧作为扭杆悬挂系统的关键部件之一，具有支撑车体，传递车体与地面作用力，缓和车轮冲击力及保证车辆行驶平稳性和乘员舒适性等作用.因此扭杆弹簧工作时，所承受的许用剪应力越大越易设计出负重轮与车体间较少刚性碰撞的悬挂装置[1-3].
由于扭杆弹簧强扭工艺的特殊性，强扭过大或者不足都会降低扭杆弹簧的力学性能，而在具体试验前强扭量往往不可预知，需要进行大量的试验验证，造成强扭试验周期过长及成本较高的问题[2].
针对这一问题，本文采用弹塑性理论分析了扭杆弹簧强扭过程中材料非线性的应力应变情况，结合有限元分析方法数值模拟了不同型号扭杆弹簧在不同强扭工艺要求下的残余变形情况；通过与试验结果对比，得出强扭精度对扭杆弹簧的影响及扭杆弹簧强扭工艺的合理性，为扭杆弹簧强扭工艺的发展提供了参考.
1 扭杆弹簧强扭理论1.1 扭杆弹簧强扭原理扭杆弹簧强扭工艺工作简图如图 1所示.其中，扭杆弹簧的直径为d，长度为L；且一端固定一端施加转矩M，其产生相对转动的转角为φ；如果在整个断面上没有产生塑性变形，则沿半径的应力遵循虎克定律[3].
图 1(Fig. 1)
 | 图 1 扭杆弹簧及其工作示意图Fig.1 Schematic of torsion bar spring and force condition |
1)当加载转矩M未超过材料的屈服极限转矩Ms时，有
 | (1) |
2)当加载转矩M超过Ms时，扭杆弹簧某一半径以外的轴表层，材料剪切应变超过屈服极限应变，但此时剪切应变随扭杆弹簧截面半径的线性变化关系仍可假定不变[3-6].如图 2所示，扭杆未超过屈服极限卸载工况时，应力沿SO直线变化；在剪切应变γs≤γ≤γd区域，剪应力并不沿着原来的ASO曲线下降，而是沿着与SO平行的直线族AB，A′B′…下降.
图 2(Fig. 2)
 | 图 2 扭杆弹簧强扭示意图Fig.2 Schematic of torsion bar spring pre-setting |
3)完全卸载时，扭杆弹簧芯部保持了原加载方向相同的残余应力，其数值可根据完全卸载时扭杆弹簧内外力矩平衡条件求得.
1.2 扭杆弹簧屈服强扭断面应力与扭矩、残余应力分析由于扭杆弹簧并非理想弹塑性材料，故在屈服阶段变形过程中表面会出现硬化现象[3].文献[3-6]提出扭杆弹簧屈服与塑性硬化过程中，断面半径方向的切应力与屈服极限τs存在近似线性关系；同时，硬化变形区的应力变化为
 | (2) |
由此可知，屈服强扭的扭矩是断面上三个部分[4](弹性变形区、塑性流动区和硬化变形区)应力作用的扭矩总和，以M′表示：
 | (3) |
 | (4) |
 | (5) |
1)弹性变形区(0≤r≤ρ0)的残余应力为
 | (6) |
 | (7) |
 | (8) |
 | (9) |
因强扭卸荷后扭杆弹簧处于平衡状态，此时扭杆弹簧内力矩也应处于平衡状态，即将式(6)~式(8)代入式(9)中可得Mr=0.
同时，扭杆弹簧在卸载后，其横断面的残余应力在靠近表层为压应力，而靠近中心为正应力，其断面必有一层残余应力为零，从而使扭杆弹簧横断面上的切应力变成如图 3所示的a-b-e-o分布，这样充分地发挥了断面内部材料的性能，提高了扭杆弹簧的承载能力[3, 7-8].
图 3(Fig. 3)
 | 图 3 扭杆弹簧不同阶段断面的应力分布Fig.3 Stress distribution at different stages of torsion bar spring section (a)-强扭卸载后；(b)-加载应力; (c)-加载后的合成应力. |
2 扭杆弹簧有限元分析2.1 扭杆弹簧参数材料选用45CrNiMoVA，且通过查阅机械设计手册可以得到低温回火处理后扭杆弹簧材料45CrNiMoVA的属性如表 1所示[3, 6].
表 1(Table 1)
 表 1 45CrNiMoVA材料参数Table 1 Material parameters of 45CrNiMoVA
| 表 1 45CrNiMoVA材料参数 Table 1 Material parameters of 45CrNiMoVA |
根据表 1参数可以计算出45CrNiMoVA弹性系数E=19.8 GPa.
本文应用不同型号扭杆弹簧进行强扭工艺加载，具体参数如表 2所示.
表 2(Table 2)
表 2 不同型号扭杆弹簧Table 2 Different models of torsion bar springs
| 表 2 不同型号扭杆弹簧 Table 2 Different models of torsion bar springs |
2.2 扭杆弹簧几何模型扭杆弹簧在履带车辆上装配的三维模型如图 4所示.由于仿真的要求，需要对模型进行简化处理和几何清理以便于网格划分等[8-10]，故可得扭杆弹簧的几何模型如图 5、图 6所示.
图 4(Fig. 4)
 | 图 4 履带车辆三维模型示意图Fig.4 Schematic of tracked vehicles three-dimensional model |
图 5(Fig. 5)
 | 图 5 扭杆弹簧实体三维几何模型示意图Fig.5 Schematic of torsion bar spring three-dimensional solid geometry |
图 6(Fig. 6)
 | 图 6 扭杆弹簧模型加载与约束情况Fig.6 Model of loaded and restrained torsion bar spring |
2.3 扭杆弹簧模型网格划分本文应用ABAQUS分析软件中的八节点六面体单元C3D8I，因为该单元在分析各向同性率相关和速率无关的大塑性变形问题中，可以较好地避免剪切自锁的现象.
考虑实际情况，此次采用8 mm单元尺寸共划分单元网格总数为11 924个，节点总数为12 980个，划分网格后的扭杆弹簧有限元模型如图 7所示.
图 7(Fig. 7)
 | 图 7 扭杆弹簧二维有限元网格划分示意图Fig.7 Schematic of two dimensions finite element meshing for torsion bar spring |
3 扭杆弹簧强扭数值模拟分析3.1 扭杆弹簧强扭原理在ABAQUS分析程序中提供了多种塑性材料选项，其中在分析材料非线性时，双线性随动强化模型最为常用.所以，在扭杆弹簧强化处理的数值模拟中，采用上述双线性模型较宜[3, 6].
3.2 数值模拟结果在ABAQUS分析程序中设置边界条件，依据表 2扭杆弹簧工艺参数进行数值模拟；此处以表 2的序号1为例进行说明.
当加载到扭转角度大于103°时，扭杆弹簧的外表层与内层进入屈服状态；通过数值模拟可得扭杆弹簧即其断面剪切应力分布如图 8所示.
图 8(Fig. 8)
 | 图 8 扭杆弹簧横截面处的应力分布图Fig.8 Stress distribution of torsion bar spring cross section |
从图 8仿真结果得出，扭杆弹簧工作部分在最大扭转角处其表面及一定内部区域应力值为1.369 GPa，大于屈服应力1.350 GPa而发生屈服，但因未大于其强度极限1.470 GPa而不发生断裂；通过以上的仿真可以得出扭杆弹簧在此工况下所达到的屈服深度.扭杆弹簧强扭后的残余应力分布如图 9所示.
图 9(Fig. 9)
 | 图 9 扭杆弹簧强扭卸载后的应力状态Fig.9 Unload stress distribution of torsion bar spring pre-setting |
当经过6次反复强扭完成后及卸载完全时，扭杆弹簧整体塑性应变状态，以及其表面和一定深度处产生塑性变形如图 10所示.
图 10(Fig. 10)
 | 图 10 扭杆弹簧强扭6次结束后横截面Fig.10 Plastic deformation distribution of cross section after six times torsion bar spring pre-setting |
根据以上ABAQUS仿真分析的强扭后塑性变形情况，可知型号1样件在6次强扭结束工况下，沿半径方向塑性变形(PEEQ)值为正且具有逐渐增大的趋势，说明表面发生了塑性变形.
利用以上方法，当扭转误差在±1°变化时，其具体的数值分析情况如表 3所示.
表 3(Table 3)
表 3 扭杆弹簧强扭仿真数据Table 3 Simulation data of torsion bar spring pre-setting
| 表 3 扭杆弹簧强扭仿真数据 Table 3 Simulation data of torsion bar spring pre-setting |
利用相同方法，可以得到其他型号扭杆弹簧的模拟数据.
4 扭杆弹簧强扭试验4.1 试验仪器扭杆试验主要应用扭杆弹簧预扭机进行，且其主要参数如下：最大扭矩45 kN·m，电机功率11 kW，电机转速1 460 r/min，最大扭转速度1.21 r/min，最大转角±180°.
4.2 试验数据此强扭试验主要针对45CrNiMoVA钢低温回火扭杆弹簧，利用此扭杆弹簧预设参数进行强扭试验.从试验结果与仿真结果数据可知，不同型号扭杆弹簧的强扭工艺参数与有限元仿真强扭工艺数据最终存在一定的误差，但总体有较好的一致性，且最大相对误差不超过15%.同时，随着扭转角度的变化，其仿真分析数值之间变化不超过5%，这表明所建模型的可靠性.
5 结论本文研究了扭杆弹簧强扭工况下的受力与变形情况，利用弹塑性理论与非线性有限元分析方法，数值模拟了扭转精度对扭杆弹簧强扭效果的影响，得出在强扭工艺允许范围内的扭转误差对扭杆弹簧强扭效果的影响不超过5%.试验研究证明：扭杆弹簧具有强扭误差在不超过±10工况下的强扭效果可近似等同的特点.
参考文献
| [1] | 丁法乾, 杨景义. 扭力轴强扭处理[J].兵工学报, 1987, 10(1) : 1–10. ( Ding Fa-qian, Yang Jing-yi. The process of torsion bar spring pre-setting[J].[J] Acta Armamentarii, 1987, 10(1) : 1–10.) |
| [2] | 周流芳, 钱士文, 曹玉成. 扭杆弹簧硬化强扭理论的研究[J].锻压技术, 2009, 5(34) : 25–131. ( Zhou Liu-fang, Qian Shi-wen, Cao Yu-cheng. Study on hardening pre-setting theory of torsion bar spring[J].Forging & Stamping Technology, 2009, 5(34) : 25–131.) |
| [3] | Paul A.Models for the dynamic simulation of tank track components[D].Bedfordshire:University of Cranfield, 2006.https://dspace.lib.cranfield.ac.uk/bitstream/1826/1106/1/Paul Allen thesis.pdf |
| [4] | Nguyen T V T.Automated torsion bar design for vehicle suspension system[D].Louisville:University of Louisville, 1999. |
| [5] | Dhir A, Sankar S. Ride dynamics of high-speed tracked vehicles[J].Vehicle System Dynamics, 1994, 23(6) : 379–409. |
| [6] | Liu C J, Zhu B S. Development of a measuring system for stiffness evaluation of micro torsion bar[J].Applied Mechanics and Materials, 2010, 3(3) : 190–194. |
| [7] | Juan M M G, Daniel F C, Víctor R C, et al. Eneralized spiral torsion spring model[J].Mechanism and Machine Theory, 2012, 1(51) : 110–130. |
| [8] | Lee H C, Choi J H. Spatial dynamics of multi-body tracked vehicles Part Ⅱ:contact forces and simulation results[J].Vehicle System Dynamics, 1998, 29(1) : 113–137. |
| [9] | Rubinstein D, Hitmn R. A detailed multi-body model for dynamic simulation of off-road tracked vehicles[J].Journal of Terramechanics, 2004, 41(5) : 163–173. |
| [10] | Han S R, Kun S H, Dae S B, et al. Development of multibody dynamics simulation tool for tracked vehicles[J].JSME International Journal, Series C, 2003, 46(2) : 540–549.DOI:10.1299/jsmec.46.540 |
