铁道车辆平稳性分析

1、铁道车辆平稳性分析1.车辆平稳性评价指标1.1 sperling平稳性指标欧洲铁路联盟以及前社会主义国家铁路合作组织均采用平稳性指数来评定车辆的运行品质。等人在大量单一频率振动的实验基础上提出影响车辆平稳性的两个重要因素。其中一个重要因素是位移对时间的三次导数,亦即z=a（加速度变化率）。若上式两边均乘以车体质量mc,并将之积改写为f，则mcz=f。由此可见,z在一定意义上代表力f的变化率的增减变化引起冲动的感觉。如果车体的简谐振动为z=z0sint,则z=-z03sint，其幅值为：zmax=z0(2f)3（1）影响平稳性指数的另一个因素是振动时的动能大小,车体振动时的最大动能为：12mcz

2、2=12mcz02=12mcz02f2=ed（2）所以：z02f2=2edmc（3）sperling在确定平稳性指数时,把反映冲动的z0(2f)3和反映振动动能z02f2的乘积25z03f5作为衡量标准来评定车辆运行平稳性。车辆运行平稳性指数的经验公式为：w=2.710z03f5ff=0.89610a3fff（4）式中 z0振幅（cm）；f振动频率（hz）；a加速度cm/s2,其值为：a=z02f2；ff与振动频率有关的加权系数。ff对于垂向振动和横向振动是不同的，具体情况见表1。表1 振动频率与加权系数关系对于垂向振动的加权系数对于横向振动的加权系f的取值范围（hz）ff公式f的取值范围（h

3、z）ff公式0.55.9ff=0.325f20.55.5ff=0.8f25.920ff=400f25.42.6ff=650f2大于201大于261以上的平稳性指数只适用一种频率一个振幅的单一振动,但实际上车辆在线路上运行时的振动是随机的,即振动频率和振幅都是随时间变化的。因此在整理车辆平稳性指数时,通常把实测的车辆振动加速度按频率分解,进行频谱分析,求出每段频率范围的振幅值,然后对每一频段计算各自的平稳性指数,然后再求出全部频率段总的平稳性指数：w=w110+w210+wn100.1（5）sperling平稳性指标等级一般分为5级，sperling乘坐舒适度指标一般分为4级。但在两级之间可按要

4、求进一步细化。根据w值来评定平稳性等级表见表2表2 车辆运行平稳性及舒适度指标与等级w值运行品质w值乘坐舒适度（对振动的感觉）1很好1刚能感觉2好2明显感觉3满意2.5更明显但无不快4可以运行3强烈，不正常，但还能忍受3.25很不正常4.5运行不合格3.5极不正常，可厌，烦恼，不能长时忍受5危险4极可厌，长时忍受有害我国也主要用平稳性指标来评定车辆运行性能，但对等级做了简化，见表3。表3 车辆运行平稳性指标与等级平稳性等级评定平稳性指标客车机车货车1优2.52.753.52良好2.52.752.753.103.54.03合格2.753.03.103.454.04.25对sperling评价方法

5、的分析：1.该评价方法仅按照某一个方向的平稳性指标等级来判断车辆的性能是不全面的,需要同时考虑垂向与横向振动对人体的生理及心理的相互影响,因为有时根据垂向振动确定的平稳性指标等级与根据横向振动确定的平稳性指标等级存在较大的差异。2.该评价方法不够灵敏。由于人体对不同振动频率的反应不同,当对应某一频率范围的平稳性指标值很大值大于,在该窄带中的振动已超出了人体能够承受的限度,但在其它频带中值都很小,由于该方向总的平稳性指标是不同振动频率的平稳性指标求和,因而可能该方向总的砰值并不大,从而认为该车辆的平稳性能符合要求是不正确的。1.2 iso2631标准1.2.1 iso2631标准概述iso263

6、1是有关人体承受振动评价的国际标准,它是由iso/tc108，即国际标准化组织机械振动与冲击标准化技术委员会的sc4一一人体承受的机械振动与冲击技术委员会指定的权威性标准,得到世界的公认。首次颁布工标准,该标准的目的是量化人体受到从固体表面传到人体过程中主要频率范围在一振动暴露极限值。应用于预测在特定频率范围内随机或非周期振动信号的频谱,自从此标准颁布以来经历了几次的修改。这些极限标准的制定是根据三条普遍公认的认知准则而来的保持舒适性、工作效率和安全或健康,三个标准分别依据三条准则定义为“减少舒适性界限”“疲劳降低工效界限”和“暴露极限”。其具体是在180hz频率范围内定义了三条区域界限（is

7、o 2631/1）.1.疲劳降低效率界限:这个界限确定了人体暴露于振动的时间极限,如果超过该极限,人们的工作就视为进行一项危险的损害工作效率的工作,特别是那些受时间影响较大的工作,如车辆驾驶员等。2.暴露时间极限:这个极限值与人体的健康和安全保护有关,在没有特别理由和事先警告,我们一般是不建议在暴露极限范围外进行工作,甚至没有任何工作任务允许在暴露极限范围外完成。3.减少舒适界限:此界限涉及到人体的舒适性保护,它是有关人在乘坐交通运输工具时,人们进行诸如吃饭、阅读和写作行为的难易程度问题。1.2.2 iso2631的几种评价方法1.2.2.1 iso2631的总的加权值评价法总的加权值评价法是

8、在某一方向上所有加速度均方根值分量的方值和根值作为评价指标。pw=i=120pwi2（6）但是,这种评价方法是建立在把人体作为一个整体接受带宽随机振动的基础上的,这样就会导致在某窄带中加速度均方根值远远超过了允许值,但在其他频带中加速度均方根值较小,由于补偿作用,使总的加权值不大。并且没有考虑不同振动方向对人体的共同影响以及可能出现某些车辆在不同车速段及不同的运行线路人体所承受的振动时间有较大差异时,导致对车辆平稳性能产生误判。1.2.2.2 采用三分之一倍频带法进行评价该方法将人体受振敏感频率080hz用三分之一倍频程法分为20个频段,较倍频带、窄带分析能更准确诊断信号。三分之一倍频带法认为

9、许多三分之一倍频带中对人体产生影响最大的,主要是由人体感觉的振动强度最大的(折算到人体敏感频带范围以后)那一个三分之一倍频带所造成。将算得的值与标准曲线对照从而得到各参数的评价值。方法简洁,便于操纵。按照这种评价方法,人能够承受的时间均为4h,由此而认为这两种车辆的平稳性能相同,显然是不合理的。因此,三分之一倍频程评价法的缺陷在于没有考虑不同频率加速度均方根值对人体的总体主观感觉的影响及不同方向振动的影响。 1.2.2.3加权加速度单值评价法 is02631标准指出振动频谱包含多个振动分量或是一个宽频带的振动时，使用加权加速度有效值方法更合适。据此，结合铁道车辆的实际振动情况，从该标准推荐的几

10、种数据处理法中选取频率分析或加权滤波网络的单值评价法。加权加速度有效值a定义为：a=i=1nkiai2（7）式中：ai一个频率组的振动加速度（m/s2）;ki频率加权函数；见表4。n频率分组数，与频谱带宽相关。表4 频率加权函数频率hz（1/3倍频程中心频率）加权系数垂向振动纵向振动1.00.51.01.250.561.01.60.631.02.00.711.02.50.80.83.150.90.634.01.00.55.01.00.46.31.00.3158.01.00.2510.00.80.212.50.830.1816.00.500.125当采用加权滤波网络时，aw等效地定义为：a=1t

11、a2tdt（8）式中：a(t)经加权滤波后的加速度时间历程。1.3 uic513r（欧洲铁路联盟标准）欧洲规范eurocode对客车车体垂向振动加速度的评定标准见表5，车体横向振动加速度没有考虑。我国均采用最大振动加速度amax和司机室振动加速度有效值aw来评定，标准如表6 表5 “eurcode”关于车体垂向振动加速度评定标准评价等级优秀良好合格车体垂向加速度cms2100130200表6 我国机车振动加速度平稳性评定等级评价等级amaxmm2awmm2垂向横向垂向横向优2.451.470.3930.273良好2.951.960.5860.407合格3.632.450.8400.5801.4

12、 gb5595-85标准与tb/t-2360-93标准分析我国制定的gb5595-85铁道车辆动力性能评定和试验鉴定规范标准基本上与平稳性指标评价法相同,因此存在的问题也类同这里的tb/t-2360-93标准主要从机车的振动加速度方面对平稳性进行侧面的反映。因为当振动加速度增大,列车的平稳性就会降低。但是,这种评价方法比较死板,不能直接反映出机车的平稳性,而且各加速度级间差距比较大,这使得测试的准确性降低。1.5 本文评价指标的选择车辆的运行平稳性是评价车辆系统动力学性能的重要指标，本文采用 sperling 运行平稳性指标，该指标基于大量试验而制定，用于评价车辆本身的运行品质和乘客乘坐舒适度

13、。其指标的大小与车辆的振动加速度和振动频率有关，横向与垂向的计算方法不同其计算方法及评定标准见第一节。2. 车辆模型建立2.1 车辆详细参数如图1所示的车辆系统动力学模型中，轴箱簧上质量被分成车体质量和构架质量。该模型得到的结果更接近于车辆的实际振动特性。需要说明的是，模型中第一悬挂刚度为车辆各轴箱弹簧刚度之和，第二系悬挂刚度为车体与构架各弹簧刚度之和。mb为两转向架构架质量之和，mc为车体质量。当然该模型也可以理解为半车模型，即一个转向架与半个车体之间的垂向振动关系，此时第一系悬挂刚度为单个转向架轴箱弹簧刚度之和，第二系悬挂刚度为车体与构架间两个弹簧刚度之和，mc为车体质量一半。铁道客车刚柔

14、模型参数含义及原始数值见表7。表7 火车转向架crh2参数含义及原始数值参数单位整备状态数值含义mct31.6车体质量lckgm21548400车体点头转动惯量mbkg3200构架质量lbkgm217652构架点头转动惯量mwkg2000轮对质量kskn/m114.65二系垂向刚度（每转向架）cskns/m120二系垂向阻尼系数（每转向架）kpkn/m1200一系垂向刚度（每轴箱）cpkns/m25一系垂向阻尼系数（每轴箱）lbm9转向架定距之半lwm1.25轴距之长lm24.5车体总长w3rad/s8.52车体垂向第一阶弯曲频率3%1.5第一阶车体弹性振型阻尼比w4rad/s172车体垂向第

15、二阶弯曲频率4%1.5第二阶车体弹性振型阻尼比dmm790磨耗型车轮直径2.2 车辆系统垂向动力学模型图1 铁道车辆垂向动力学模型对于图1所示的车辆系统模型，当车辆处于平衡状态时（此时重力与弹簧力平衡）且自由振动时，其运动微分方程为：（9）式中：zc车体垂向位移；zc构架垂向位移；zc车轮垂向位移，即路面输入位移。将上式化成矩阵形式，上式可改写为：mz+cz+kz=0（10）式中：x为状态向量：z=（zc zb zw）m为质量矩阵：m=mc000mb000mwc为阻尼矩阵：c=cs-cs0-cscs+cpcp0-cp-cpk为刚度矩阵：k=ks-ks0-ks-ks+kp-kp0-kpkp在实际

16、车辆运行过程中，客车会受到路面对其的激励，其轮轨界面存在外加激扰，将会对系统产生影响因此上式的方程右端将不完全为0，通常我们将其表示为：mz+cz+kz=f（11）在不考虑钢轨振动，可一定程度上将轨道垂向不平顺视为钢轨位移，轮轨垂向作用力可由赫兹非线性弹性接触理论确定：fwt=1gxw(t)32（12）式中：fwt轮轨垂向作用力；xw(t)路面不平顺时域激励信号。g轮轨接触常数，g=3.68r-0.11510-8mn23根据参考文献23轮轨之间的接触线性化刚度为：kh=321gp013p0为近似轴重：p0=mc+2mb+4mwg/8故激励矩阵为：f=kh00xwt2.3车辆路面激励模型对于上式

17、，为系统在零路面激励下的车辆模型方程，为了求出随机输入下的振动响应谱，首先建立时域上的路面激励模型，在此我们根据根据我国p50钢轨无缝轨道不平顺的实测数据,推荐供设计计算时用的轨道不平顺功率谱密度函数设计路面激励时域模型。路面垂向不平顺功率谱密度函数如下式：ssvfs=2.75510-9fs2+8.87810-1fs4+2.52410-2fs2+9.6110-7（13）式中：fs空间频率，cycle/m;ssvfs空间域内路面功率谱函数；将上述空间域内功率谱密度函数，转换为时间域内功率谱密度函数，再利用白噪声通过滤波器的方法实现对对路面模型的数数学模型的建立。根据参考文献，利用遗传算法得到列车

18、运行速度为200 km/h时对应于轨道垂向不平顺的优化的滤波器方程为：gvs=1.135710-46434.68s2+478.45s+168.20s3+60.44s2+121.54s（14）将白噪声函数通过上述滤波器可得到路面功率谱的时域模型的近似表示。利用matlab建立模型如下：图2 轨道谱时域模型图3 模拟轨道的时间序列3.在轨道谱激励下的车辆稳定性分析3.1车辆模型时域分析在建立基于轨道功率谱的路面激励时域模型，结合第二节中车辆在路面模型激励下的运动微分方程我们可以利用simulink建立车辆仿真模型如下图：图4 车辆仿真模型在前两节给定参数下得出车辆车体及构架振动加速度仿真曲线：图5

19、 车体加速度图6 构架加速度3.2车身振动加速度频域分析下面我们求车辆的平稳性指标，由于sperling指标需要车辆车体垂向加速度的幅频特性，故在上面模型的基础上我们将车体加速度数据利用simulink中的to workplace模块导入matlab工作空间，在将数据通过matlab编程处理进行傅立叶变换后我们得到了车体加速度的幅频特性图如下：图7 车体加速度幅频特性根据上图我们可以得到车辆加速度及其频率的数据，并计算其平稳性评价指数w见下表：表8 垂向平稳性指数w计算参数及结果加速度幅值cm/s2加速度频率hz平稳性指数2.96111.10902.63151.25731.80271.0857

20、1.30790.91442.693111.06951.860190.81234.166231.00472.496310.83631.813400.64911.314570.7407根据各频率段w值，利用下式可以计算出车辆全频率段的平稳性指数w=w110+w210+wn100.1=1.3404由表2及表3可以看出该车的运行品质为好，平稳性等级为1级。参考文献1 杨岳,张晓峰,张兆丰,徐勇. 面向运行平稳性的铁道车辆悬挂参数灵敏度分析j. 铁道科学与工程学报,2014,(02):116-120.2 宫岛,周劲松,孙文静,沈钢. 基于格林函数法的铁道车辆弹性车体垂向振动分析j. 机械工程学报,201

21、3,(12):116-122.3 王辉,钟晓波,沈钢. 基于simulink的柔性车体垂向振动特性分析j. 噪声与振动控制,2012,(05):79-82+140.4 孙文静,周劲松,宫岛. 弹性车体垂向运行平稳性一系最优控制研究j. 振动与冲击,2012,(12):150-154+164.5 张学铭,谢晓波,杨少彬. 铁道车辆平稳性指标对比分析j. 佳木斯大学学报(自然科学版),2011,(05):661-666.6 贾璐,曾京,池茂儒. 车辆系统横向运动稳定性评判的数值仿真研究j. 铁道车辆,2011,(09):1-7.7 廖英英,刘金喜,刘永强,杨绍普. 半主动控制与时滞对高速铁道车辆平

22、稳性、稳定性及安全性的影响j. 振动与冲击,2011,(06):53-578 宫岛,周劲松,孙文静,谢维达. 铁道车辆弹性车体垂向运行平稳性最优控制j. 同济大学学报(自然科学版),2011,(03):416-420.9 周劲松,宫岛,任利惠. 铁道车辆弹性车体被动减振仿真分析j. 同济大学学报(自然科学版),2009,(08):1085-1089.10 周劲松,宫岛,孙文静,任利惠. 铁道客车车体垂向弹性对运行平稳性的影响j. 铁道学报,2009,(02):32-37.11 周劲松,李大光,张祥韦,沈钢. 平稳性快速算法及其在高速铁道车辆动力学分析中的运用j. 铁道学报,2008,(06):

23、36-39.12 池茂儒,张卫华,曾京,金学松,朱旻昊. 蛇行运动对铁道车辆平稳性的影响j. 振动工程学报,2008,(06):639-643.13 周劲松,张洪,任利惠. 模态参数在铁道车辆运行平稳性研究中的运用j. 同济大学学报(自然科学版),2008,(03):383-387.14 刘李莹,周文祥,徐娜. 铁路车辆运行平稳性指标的测试精度分析j. 电力机车与城轨车辆,2008,(01):38-40+56.15 张洪,周劲松,任利惠,杨国桢,沈钢. 基于运行模态参数辨识的客车运行平稳性研究16 刘转华. 铁道车辆运行平稳性评价方法研究d.西南交通大学,2007.17 万鹏. 考虑轮对弹性时车辆系统动力学建模与仿真分析d.西南交通大学,2008.18 段