熊军林,韩世昌
(1.昆明理工大学 机电工程学院, 昆明 650500;2.云南省先进装备智能控制及应用国际联合实验室, 昆明 650500)
列车车辆动力吸振器是利用质量块、弹簧以及阻尼器在列车主振系产生“反共振”现象来降低主振系振幅,并提升列车行驶安全性和舒适性的一种减振装置[1]。许多研究学者基于Voigt-DVA(voigt-dynamic vibration absorber)模型[2-3],从结构方式和实际应用中创新[4-8],并加强了DVA的减振性能,但是在实际能应用到传统DVA的工程中,DVA的过大振子重量会造成列车运输成本提高,且较大体积也会影响工程实现。
2002年,Smith提出了惯容器概念,一种能在运转过程中以较小的自身重量实现较大“虚质量”的物理元件,并设计出齿轮齿条式惯容器和滚珠丝杠式惯容器样机[9-11]。随后,Smith将惯容器与弹簧、阻尼结合应用于车辆的减振系统提升车辆的稳定性,验证了ISD(inerter spring damper)悬架的减振性能优于传统的SD(spring damper)悬架[12-13]。目前为止,由于结构的原因,滚珠丝杠式惯容器实现“虚质量”的能力大,且结构简单[14],易于运用在惯容特性需求大的工程系统中。随后,杨晓峰等[15]利用惯容器改进传统DVA的结构,解决了车辆悬架系统应用DVA的空间布置问题,实现从DVA结构转变为ISD悬架结构。葛正等[16]采用主动式惯容器对DVA进行了控制优化,提升了DVA的减振频带,提高了DVA的使用价值。综上,DVA能够实现对列车的减振作用,但DVA振子质量对实际工况会产生影响,并且惯容器与DVA的结合也具有一定减振效果。因此,从惯容器具有“虚质量”的特性入手,在性能相同的基础上通过惯容器介入列车DVA来替换振子实际质量的研究仍值得进一步探索。
本文中利用惯容器具有“虚质量”的物理特性,以DVA的减振性能为分析目标,建立了3种列车模型,并利用传统DVA参数经验公式以及智能优化算法得到DVA的结构参数,通过车体加速度均方根值对比分析了惯容器对DVA振子质量的影响。
ISD悬架在列车中的应用已经被广为证实,其中惯容器串联阻尼器后再与弹簧并联的结构较其他结构对列车的稳定性提升明显[17],本文中将在此基础上对含有ISD悬架列车进一步研究。图1为含ISD悬架的列车垂向振动模型,该模型将列车简化为基于轮对、转向架、车体的多刚体系统,并且包含了弹簧、阻尼器并联的一系隔振结构,ISD悬架的二系隔振结构,其中ISD悬架采用惯容器bs串联阻尼器cs后再与弹簧ks并联的结构。
图1 含ISD悬架的列车垂向振动模型
图1中,mt为转向架质量,mc为车体质量;kp、cp分别为一系垂向刚度和垂向阻尼系数;ks、cs、bs分别为二系垂向刚度、垂向阻尼系数和惯质系数,其中惯质系数为惯容器的惯性力与其两端点间相对加速度的比值;Zi为系统的垂向位移变量;Zwi为轨道不平度。
对简化的列车模型,取向下为正方向,根据牛顿第二运动定律推导出系统动力学方程。
前转向架沉浮运动方程:
(1)
级间垂向位移Z5沉浮运动方程:
(2)
后转向架沉浮运动方程:
(3)
级间垂向位移Z6沉浮运动方程:
(4)
车体沉浮运动方程:
(5)
在列车上安装DVA能够使得列车通过DVA自身振动给予的反方向力来实现减振效果,且可以定点安装。传统DVA选择Voigt-DVA模型[18],是由弹簧并联阻尼器后与DVA振子串联的结构,本文中选择安装于车体中间,如图2虚线框所示。
图2 传统DVA的ISD悬架列车模型
图2中,在含ISD悬架的列车垂向振动模型基础上,加入传统DVA,其中,k、c分别为传统DVA的弹簧刚度、阻尼系数;mV为DVA的振子质量。
安装了传统DVA后,列车系统的转向架沉浮运动方程和级间垂向位移运动方程与不安装DVA的模型方程一致,见式(1)—式(4)。根据牛顿第二运动定律推导出增加的系统动力学方程如下。
车体沉浮运动方程为
(6)
传统DVA振子沉浮运动方程为
(7)
传统DVA的最优参数采用定点理论设计,并且传统DVA的参数设计受到安装位置以及减振目标质量的影响,因此基于含ISD悬架的列车垂向振动模型,由文献[19]可得最优传统DVA最优参数公式为
(8)
(9)
(10)
式(8)—式(10)中:γo为DVA的频率比;μ为DVA振子质量与车体质量之比;ζo为DVA的阻尼比;mV为DVA振子的质量;n为正整数,对轨道列车可取值为3;v为列车行驶速度;lb为转向架定距之半;k、c分别为DVA的弹簧刚度和阻尼系数。
1.3.1列车模型
传统DVA振子质量一般为水泥和金属块,质量大、空间占用率多,制约着DVA的使用。因此本文中结合惯容器具有“虚质量”的物理特性,以及惯容器串联阻尼后再与弹簧并联的DVA结构,替换了传统DVA的减振结构,并提出了含惯容器式DVA的ISD悬架列车模型,如图3所示,虚线框内为含惯容器式DVA结构。
图3 含惯容器式DVA的ISD悬架列车模型
图3中,在含ISD悬架的列车垂向振动模型基础上,加入含惯容器式DVA,其中,k、c、b分别为含惯容器式DVA的弹簧刚度、阻尼系数和惯容器惯质系数。
安装了含惯容器式DVA后,列车系统的转向架沉浮运动方程和级间垂向位移运动方程与不安装DVA的模型方程一致,见式(1)—式(4)。根据牛顿第二运动定律推导出系统增加的动力学方程如下。
车体沉浮运动方程为
(11)
含惯容器式DVA振子沉浮运动方程为
(12)
级间垂向位移Z7沉浮运动方程为
(13)
1.3.2含惯容器式DVA中惯容器的设计
在含惯容器式DVA结构中,需要实现振子质量被惯容器的“虚质量”代替,且惯容器本身质量对DVA整体质量影响较小,则需要采用 “虚质量”与本身质量之比较大的惯容器,因此在本文中选择滚珠丝杠式惯容器。
滚珠丝杠式惯容器主要由惯容器两端固定的外壳(左)和外壳(右)、滚珠丝杠、螺母、滚动轴承以及飞轮组成,如图4所示。
图4 滚珠丝杆式惯容器
其动力学方程为
(14)
式(14)中:ωb为丝杠旋转角速度;vb为惯容器两端间的相对速度;P为滚珠丝杠副导程;T为丝杠驱动力矩;Jb为惯容器中飞轮的转动惯量;F为惯容器两端的受力;a为惯容器两端点间的相对加速度。
根据式(14)可得滚珠丝杠式惯容器的惯质系数b的计算式为
(15)
常见滚珠丝杠式惯容器中的飞轮为圆柱体,设飞轮质量为m,半径为r,由圆柱体转动惯量计算公式,则惯质系数还可以表示为
(16)
基于以上模型,以含ISD悬架的列车垂向振动模型为基础,分别对安装传统DVA和含惯容器DVA的模型进行定点理论计算和智能优化算法优化,最后对结果进行对比分析。
以含ISD悬架的列车垂向振动模型为对比分析基础,为了保证该模型的最优状态,采用智能优化算法实现参数优化,优化目标为一系kp、cp,二系ks、cs、bs。同样由于在含惯容器式DVA的ISD悬架列车模型中,不适用传统DVA设计方法,为实现优化设计,也采用智能优化算法,优化目标为DVA中的k、c、b。
粒子群优化算法(particle swarm optimization PSO)是一种模拟鸟群觅食行为的随机搜索算法,优化能力强,应用方便[20]。在粒子群算法中,粒子根据个体极值和群体极值更新自身速度和位置,在预设的迭代次数下,寻找出最优值,并且随着迭代次数的增加,适应度函数值趋于直线。流程如图5所示。
图5 粒子群优化算法流程框图
采用德国轨道谱高低不平度高干扰谱作为系统输入,速度为300 km/h。根据参考文献[17],得到时域计算公式为
(17)
利用Matlab/Simulink,考虑实际工况中前后轮对的间隔,得到德国谱位移时域信号。
基于列车的垂向振动模型,以乘坐的舒适性作为减振系统设计重点,在简化模型中,选择车体垂向加速度的均方根值作为性能评价指标。
车体加速度均方根值越小越能提升乘坐的舒适性。因此,优化目标为
(18)
2.3.1含ISD悬架的列车垂向振动模型参数优化
简化的列车垂向振动模型中,所需的基础参数如表1所示。
表1 列车-轨道垂向振动模型参数
利用粒子群优化算法,以含ISD悬架的列车垂向振动模型中车体加速度均方根值为目标进行参数优化,其中优化变量和算法约束条件如表2所示。
表2 悬架元件参数取值范围
在300 km/h的德国轨道谱激励路面下,通过粒子群优化算法得到优化结果如表3所示。
表3 含ISD悬架的列车垂向振动模型优化结果
2.3.2含惯容器式DVA的ISD悬架列车模型参数优化
由文献[21]可知,当传统DVA的振子与车体的μ为0.1时,系统能够同时满足最优同调和最优阻尼2个条件,使得系统在原来共振处的幅值大大减小,出现的2个新共振峰的峰值也很低,达到了定点理论参数优化的良好减振效果。当传统DVA的振子与车体的μ超过0.1时,会提高DVA的工程应用成本和降低DVA的减振效果。
因此在本文中,取传统DVA的振子与车体的μ等于0.1时为边界条件,在μ小于0.1时研究惯容器“虚质量”特性对DVA振子质量的替换范围,在μ等于0.1时研究含惯容器式DVA对传统DVA的减振优化。
为了在含ISD悬架的列车垂向振动模型最优参数基础上,实现惯容器“虚质量”替代DVA振子质量的目的,本文选择在含惯容器式DVA的振子质量mV从一定初始值递增到μ等于0.1的过程中,以及300 km/h的德国轨道谱激励路面下,通过粒子群优化算法得到每一次mV递增时含惯容器式DVA的元件最优参数,随后与传统DVA设计下的列车系统作对比分析。优化参数范围如表4所示。
表4 含惯容器式DVA的ISD悬架列车模型参数范围
在mV从1 000 kg开始步长为100 kg递增到2 600 kg(即μ达到0.1)的过程中,每一次mV递增时含惯容器DVA的元件最优参数如表5所示。
表5 含惯容器式DVA的ISD悬架列车模型优化结果
在经过粒子群算法参数优化后含ISD悬架的列车垂向振动模型基础上,利用上述优化数据以及传统DVA参数设计方法,对比分析含惯容器式DVA和传统DVA在含ISD悬架的列车垂向振动模型的性能影响,并借此讨论惯容器在DVA结构中对DVA振子的质量替代范围。
结合悬架的优化参数,根据式(8)—式(10),在μ为0.1时,可得传统DVA的ISD悬架列车模型的最优减振参数。如表6所示。
表6 传统DVA的ISD悬架列车模型参数取值
经过上述分析可以得到不同模型的车体加速度结果,汇总后如图6所示。其中,以含ISD悬架和未添加DVA的列车垂向振动模型的车体加速度均方根值作为对比基础,如图6中长实线所示,值为0.026 9 m/s2。
图6 惯容器对DVA的影响
在DVA振子质量从1 000 kg达到2 600 kg过程中,传统DVA的ISD悬架列车模型在传统DVA参数计算方法下得到传统DVA对含ISD悬架的列车垂向振动模型的减振效果,如图6中“传统DVA”所示。
在DVA振子质量从1 000 kg达到2 600 kg过程中,含惯容器式DVA的ISD悬架列车模型结合表5的优化结果,得到含惯容器式DVA对含ISD悬架的列车垂向振动模型的减振效果,如图6中“含惯容器式DVA”所示。
由图6可知,传统DVA和含惯容器式DVA对含ISD悬架的列车垂向振动模型都有性能提升。相较于没有DVA的模型,在所有振子质量取值情况下车体加速度均方根值都有所降低。同时,优化后的含惯容器式DVA较传统参数设计所得的传统DVA对车辆减振有更好的效果。
振子质量mV在2 100~2 600 kg内,含惯容器式DVA结构对传统DVA结构的性能提升如表7所示。
表7 含惯容器式DVA对传统DVA的优化
振子质量mV在2 100~2 200 kg内,含惯容器式DVA结构通过粒子群优化算法达到了传统DVA结构参数经验计算方法下mV为2 600 kg时给列车带来的减振效果,如图6中局部放大图所示,其中虚线处的值为0.023 98 m/s2,且mV在2 200 kg时含惯容器式DVA结构减振效果好于mV在2 600 kg时的传统DVA结构,综上表明了惯容器 “虚质量”的物理特性是等效了接近400 kg DVA振子的质量。并且由表7可知,振子质量mV在2 100~2 600 kg的变化过程中,含惯容器式DVA比传统DVA结构对含ISD悬架的列车垂向振动模型的减振性能提升要好。
为了进一步验证惯容器质量替换的可行性,从加速度均方根值接近的不同振子取值来分析时域和频域上的车体加速度效果。取传统DVA的ISD悬架列车模型在mV为2 600 kg时和含惯容器式DVA的ISD悬架列车模型在mV为2 200 kg时车体加速度时域和频域作对比,如图7、图8所示。
图7 车体加速度随时间变化曲线
图8 车体加速度功率谱变化曲线
图7为车体的加速度随时间变化的曲线,由图7可以看出,在仿真的10 s内,时域的变化基本一致;图8为车体加速度功率谱变化曲线,在仿真的70 Hz以内取0~15 Hz内对比,频域的变化也基本一致。说明了含惯容器式DVA的振子质量在2 200 kg时已经达到了传统DVA的振子质量在2 200 kg时对列车的减振效果,同时反映了惯容器能够替换DVA结构中振子的质量并且不损失结构的减振性能。
由表4可知,粒子群优化算法约束确定时,为实现更好的优化效果,b的优化约束范围设置较大,优化结果也表明了b的优化值需要采用 “虚质量”与本身质量之比较大的惯容器实现。
滚珠丝杠式惯容器在选择一定材料,调整r以及P的前提下,能够实现m到b较大的放大倍数。表5中能够用惯容器“虚质量”替换DVA振子质量的有效优化范围是mV从2 100~2 200 kg的中间值到2 600 kg,取该范围中惯质系数的最大值4 455 297.55kg,根据式(15)可得,滚珠丝杠式惯容器主要参数其中一种取值,如表8所示。表8的配置可以计算得到飞轮质量为91.84 kg,相较于4 455 297.55 kg的惯质系数,滚珠丝杠式惯容器实现了较小的飞轮质量得到较大 “虚质量”的物理特性。
表8 滚珠丝杠式惯容器主要参数取值
同时由表7可知,惯容器替换传统DVA振子的质量在400~500 kg,由表8可知,根据最大惯质系数实现出来的滚珠丝杠式惯容器在91.84 kg左右,粗略计算可得实际惯容器能替换传统DVA振子的质量在300~400 kg,且逼近400 kg。
本文中基于粒子群优化算法,对含ISD悬架的列车垂向振动模型及含惯容器式DVA结构进行了参数优化,并与传统DVA结构进行了对比分析,得到如下结论:
1) DVA的加入比没有DVA时候的车辆性能要有所提升,且含惯容器式DVA结构在振子质量相同的情况下比传统DVA结构减振效果好;
2) 在保证相同减振性能的条件下,含惯容器式DVA能够用惯容器的“虚质量”物理特性实现对DVA振子的质量替换,替换质量接近400 kg;
3) 滚珠丝杠式惯容器可根据具体的实物参数设计达到百万级的惯质系数,且自身体积质量比较小。