2. 东北大学 机械工程与自动化学院, 辽宁 沈阳 110819;
3. 奇瑞捷豹路虎汽车有限公司, 上海 210103;
4. 沈阳工业大学 后勤集团, 辽宁 沈阳 110870
2. School of Mechanical Engineering & Automation, Northeastern University, Shenyang 110819, China;
3. Chery Jaguar Land Rover Automotive Co., Ltd., Shanghai 210103, China;
4. Logistics Group, Shenyang University of Technology, Shenyang 110870, China
车辆悬架系统连接车身与车轮, 对缓和路面冲击、提升车辆行驶稳定性具有重要意义[1-2].对车辆悬架系统减振性能的改善, 始终是车辆工程领域研究的热点问题[3-5].传统的车辆悬架结构由减振元件弹簧和阻尼器并联构成, 具有结构简单、性能可靠、造价低等优点, 但是, 传统车辆悬架因其独特的结构限制了减振性能的提高.添加控制系统与能量输入的主动、半主动悬架, 可以有效提高悬架系统的减振性能, 同时也带来了能耗大、控制系统时滞等问题.惯容器的发明为车辆悬架系统减振性能的提升提供了新的可能.
Smith等[6]将减振元件惯容器应用于车辆悬架系统, 研究了应用惯容器的多种悬架拓扑结构的性能改善.Zhang等[7]对应用惯容器的车辆悬架结构进行分析, 研究了不同控制策略对悬架减振性能的影响.Wang等[8]探究了摩擦与弹性效应等非线性因素对含惯容器的车辆悬架系统减振性能的影响, 研究结果表明,受惯容器非线性因素影响的车辆悬架系统的减振性能仍然优于传统车辆悬架.Kuznetsov等[9]为改善乘坐舒适性, 建立了含惯容器的车辆悬架系统模型, 以加权加速度为评价指标对参数进行优化.李小彭等[10]提出了一种每一级悬架中均含有惯容器的两级汽车悬架, 进一步强化了惯容器在车辆悬架系统中起到的作用.Hu等[11]探究了惯容器、弹簧、阻尼器不同组合形式对减振性能改善的影响.
多数研究中, 通过调整减振元件惯容器、弹簧、阻尼器之间的组合关系来改善悬架的减振性能, 对综合考虑车辆多种性能评价指标的优化设计还有进一步探究的空间.本文基于含惯容器的车辆悬架系统进行优化设计, 通过对减振元件的参数进行优化及对衬套进行优化选型, 进一步改善车辆悬架系统的减振性能.研究结果将为含减振元件惯容器的车辆悬架系统的实际应用提供理论参考.
1 理论基础 1.1 惯容器减振原理减振元件惯容器具有质量属性, 同时惯容器突破了质量块单端子接地的限制, 是一种双端子减振元件[12].滚珠丝杠式惯容器具有结构简单、成本低、便于安装等优点, 是一款被广泛应用的惯容器.
滚珠丝杠式惯容器的理想动力学模型可以表示为
(1) |
式中:F为惯容器两端所受到的力; m为飞轮的质量; r为飞轮的半径; p为滚珠丝杠副的导程; v为惯容器两端的相对运动速度; b为惯容器的惯质系数.
1.2 振动功率流理论功率可表示为
(2) |
式中:F(t)为力的瞬时值; v(t)为速度的瞬时值.
(3) |
(4) |
稳态功率流的表达式为
(5) |
式中:ω为振动圆频率; F为力的幅值; V为速度的幅值; *表示共轭.
根据机械阻抗理论, 机械阻抗Z可表示为
(6) |
由式(5)与式(6)可以看出, 机械阻抗可决定流经该点的功率, 系统本身的动态特性直接影响到振动能量的传递.
2 车辆悬架系统建模为了分析含减振元件惯容器的车辆悬架系统的动力学响应, 建立如图 1所示的车辆悬架系统模型.
图 1所示的车辆悬架系统模型的动力学方程可表示为
(7) |
式中:ms为簧载质量; mu为非簧载质量; F为惯容器两端受到的力; ks为主弹簧刚度; cs为主阻尼系数; kt为轮胎刚度; b为惯质系数; ki为副弹簧刚度; ci为副阻尼系数; xs为车身垂向位移; xb为惯容器垂向位移; xu为轮胎垂向位移; xt为路面垂向位移.
采用滤波白噪声法生成的路面垂向位移输入为
(8) |
式中:G0为路面不平度系数; v为车辆行驶速度; w(t)为均值为零的高斯白噪声信号; f0为下截止频率.
将式(7)转换为式(9)所示的状态空间方程的一般表示形式:
(9) |
式中:A,B,C为系数矩阵; U为系统输入变量; X为系统状态变量; Y为系统输出变量.
选取系统的状态变量:
(10) |
其中
系统的状态方程组:
(11) |
系数矩阵:
(12) |
式中,
(13) |
(14) |
式中,
系统输入变量:
(15) |
系统输出变量:
(16) |
车辆悬架系统的模型参数如表 1所示.
为进一步提高车辆悬架的综合减振性能,采用遗传算法对惯质系数b、副弹簧刚度ki、副阻尼系数ci进行参数优化.
综合考虑车身加速度、悬架动行程、轮胎动载荷指标的适应度函数为
(17) |
式中:J为适应度函数值; A(x)为车身加速度的均方根值; L(x)为悬架动行程的均方根值; T(x)为轮胎动载荷的均方根值.
设定的约束条件为
(18) |
式中:A(x)past,L(x)past,T(x)past分别为车身加速度、悬架动行程、轮胎动载荷在优化前参数下的均方根值.
优化后的参数如下:
含衬套的车辆悬架系统机械网格如图 2所示.
从图 2可以得到左侧的等效阻抗为
(19) |
右侧的等效阻抗为
(20) |
则连接点A处的力为
(21) |
式中:Fc为恒力源,取Fc为简谐单位力.
连接点A处的速度为
(22) |
流经A点的功率流为
(23) |
为了验证所建立的车辆悬架系统模型的正确性, 在仿真软件ADAMS/car中建立车辆悬架系统的多体动力学模型, 并选取原始设计参数kb=450 N/mm, cb=0.50 N·s·mm-1.通过理论计算与仿真数据得到的功率流对比如图 3所示.
从图 3可以看出, 随着频率的增大, 通过理论计算和仿真数据得到的功率流谱的变化趋势基本一致, 证明了模型的正确性.
4 数值仿真 4.1 悬架系统的动力学响应为表征采用遗传算法优化后得到的悬架参数对车辆悬架系统综合减振性能的影响, 以优化前设计参数下的车辆悬架系统动力学响应为比较对象, 将优化前后车辆悬架系统的动力学响应进行对比分析.
参数优化前后车辆悬架系统的动力学响应如图 4所示.
从图 4可以看出, 与优化前相比, 优化后的车身加速度、悬架动行程、轮胎动载荷指标均有较大程度的改善, 优化后的车辆悬架系统的综合减振性能得到了显著提升.
4.2 衬套参数对功率传递的影响由于橡胶件的特殊性, 衬套的动态特性需要通过测试来获得, 必须模拟其实际工况, 然后采用专用设备进行测试.因此一味地按照优化方法对衬套的刚度、阻尼进行优化设计, 优化得到的参数可能不适合生产或生产周期过长, 优化设计缺乏实际意义.因此, 在实际应用过程中针对现有的衬套参数进行优化选型, 从而完成悬架衬套参数的优化匹配.
市场上主流的3种车辆悬架的衬套参数如表 2所示.
依据表 2中不同类型衬套的刚度、阻尼系数, 基于振动功率流理论探究衬套参数对系统振动功率的影响.不同刚度下的功率流谱如图 5所示:在0~1 Hz频段内, 衬套的刚度变化对传入车身的振动功率流的影响较小; 在0~50 Hz的仿真频段内, kb=300 N/mm时的功率流谱始终低于其他两种类型的功率流谱.
不同阻尼系数下的功率流谱如图 6所示:随着衬套阻尼系数的增大, 车辆悬架系统功率流谱的振动幅值逐渐减小.在0~50 Hz的仿真频段内, 衬套的阻尼系数取cb=2.00 N·s·mm-1时的功率流谱始终低于其他两种类型的功率流谱.
综上所述, 相比于衬套1与衬套2, 衬套3更加适合减小经车辆悬架系统传递至车身的能量, 使车辆悬架系统达到更好的减振效果.
5 结论1) 优化后得到的参数明显优于优化前的设计参数, 优化后车辆悬架系统的综合减振性能得到了显著提升.
2) 通过理论计算与仿真数据两种方式得到的车辆悬架系统的功率流谱的变化趋势基本一致, 验证了所建模型的正确性.
3) 衬套3的减振性能优于衬套1与衬套2, 更加适合应用于车辆悬架减振系统中.
[1] |
Mahmoodabadi M J, Farhadi F, Sampour S. Firefly algorithm based optimum design of vehicle suspension systems[J]. International Journal of Dynamics and Control, 2019, 7(1): 134-146. DOI:10.1007/s40435-018-0453-8 |
[2] |
李武杰, 陈从根, 郭立新. 基于微分几何法的主动悬架鲁棒H∞控制[J]. 东北大学学报(自然科学版), 2019, 40(5): 716-721. (Li Wu-jie, Chen Cong-gen, Guo Li-xin. Robust H∞ control of active suspension based on differential geometry[J]. Journal of Northeastern University (Natural Science), 2019, 40(5): 716-721.) |
[3] |
Sun S S, Tang X, Yang J, et al. A new generation of magnetorheological vehicle suspension system with tunable stiffness and damping characteristics[J]. IEEE Transactions on Industrial Informatics, 2019, 15(8): 4696-4708. DOI:10.1109/TII.2018.2890290 |
[4] |
Zheng X Y, Zhang H, Yan H C, et al. Active full-vehicle suspension control via cloud-aided adaptive backstepping approach[J]. IEEE Transactions on Cybernetics, 2019, 49(Early Access): 1-12. |
[5] |
段亮, 杨树凯, 宋传学, 等. 平衡悬架钢板弹簧动态特性的研究[J]. 机械工程学报, 2016, 52(6): 153-158. (Duan Liang, Yang Shu-kai, Song Chuan-xue, et al. Research on dynamic characteristics of the tandem suspension leaf spring[J]. Journal of Mechanical Engineering, 2016, 52(6): 153-158.) |
[6] |
Smith M C, Wang F C. Performance benefits in passive vehicle suspensions employing inerters[J]. Vehicle System Dynamics, 2004, 42(4): 235-257. DOI:10.1080/00423110412331289871 |
[7] |
Zhang X J, Ahmadian M, Guo K H. On the benefits of semi-active suspensions with inerters[J]. Shock and Vibration, 2012, 19(3): 257-272. DOI:10.1155/2012/640275 |
[8] |
Wang F C, Su W J. Impact of inerter nonlinearities on vehicle suspension control[J]. Vehicle System Dynamics, 2008, 46(7): 575-595. DOI:10.1080/00423110701519031 |
[9] |
Kuznetsov A, Mammadov M, Sultan I, et al. Optimization of improved suspension system with inerter device of the quarter-car model in vibration analysis[J]. Archive of Applied Mechanics, 2011, 81(10): 1427-1437. DOI:10.1007/s00419-010-0492-x |
[10] |
李小彭, 李凡杰, 曹洲, 等. 含惯容器的两级汽车悬架振动性能[J]. 东北大学学报(自然科学版), 2019, 40(10): 1448-1453. (Li Xiao-peng, Li Fan-jie, Cao Zhou, et al. Vibration performance of two-stage vehicle suspension with inerter[J]. Journal of Northeastern University (Natural Science), 2019, 40(10): 1448-1453. DOI:10.12068/j.issn.1005-3026.2019.10.015) |
[11] |
Hu Y, Chen M Z Q. Performance evaluation for inerter-based dynamic vibration absorbers[J]. International Journal of Mechanical Sciences, 2015, 99: 297-307. DOI:10.1016/j.ijmecsci.2015.06.003 |
[12] |
Smith M C. Synthesis of mechanical networks:the inerter[J]. IEEE Transactions on Automatic Control, 2002, 47(10): 1648-1662. DOI:10.1109/TAC.2002.803532 |