车辆可靠性

基于能量分布的轨道车辆动力学系统可靠性分析郑明亮(浙江理工大学杭州310018）摘要：为研究轨道车辆随机振动系统动力可靠性，根据分析动力学中的Hamilton正则方程，建立考虑轨道不平顺与结构自身参数随机因素共同作用的车轨悬挂2自由度复合随机振动模型。将系统的广义能量停留在安全域内的概率作为动力可靠性指标，运用拟不可积Hamilton系统随机平均法，求出系统动力响应Hamilton函数幅值的均值和方差。依据顺序统计理论和随机变量函数的矩法求出系统动力可靠度的均值和方差的计算公式。通过算例，分析轨道轨道不平顺和车辆物理参数的随机性对系统动力可靠度随机性的影响，为系统结构参数优化设计和最优控制策略等奠定基础，对保证高速列车更安全平稳运行提供参考。关键字：轨道车辆；Hamilton；均值和方差；动力可靠性ReliabilityanalysisofrailwayvehicledynamicssystembasedonenergydistributionZhengMingLiang(ZhejiangSci-TechUniversity,Hangzhou310018)Abstract:Inordertostudythedynamicreliabilityoftherandomvibrationsystemofthetrackvehicle,therandomvibrationmodelofthe2degreeoffreedomofthevehiclerailsuspensionisestablishedbyconsideringtherandomfactorsofthetrackirregularityandthestructuralparameters.Thegeneralizedenergyofthesystemisusedasadynamicreliabilityindex,andthemeanandvarianceoftheamplitudeofthesystemisobtainedbyusingthestochasticaveragingmethodoftheHamiltonsystem.Thecalculationformulaofthemeanandvarianceofthesystemdynamicreliabilitybasedontheorderstatisticstheoryandthemomentmethodoftherandomvariablefunction.Effectsofrandomthroughexamplesandanalysisoftrackirregularityandvehiclephysicalparametersofthepowersystemreliabilitydegreeofrandomness,asthesystemstructureparametersoptimizationdesignandoptimalcontrolstrategytolayafoundation,toensurethehighspeedtrainmoresafeandstableoperationofthereference.Keywords:Railvehicles;Hamilton;meanandvariance;dynamicreliability0引言：随着当前我国高速铁路和城市轨道交通的快速发展，车辆行驶安全可靠性是现代高速列车选择的基本指标之一，高速轨道车辆重载运行时，轮轨间作用力突增，车辆悬挂系统的垂向受迫振动加剧，当系统振动能量达到储存极限时，系统可靠性会被破坏，车辆发生首次穿越脱轨失效，严重会造成重大伤亡事故。因此，有必要对高速列车的可靠性展开深入研究，而车辆系统的参数分析对于改善其可靠性具有重要参考价值。目前轨道车辆系统的研究主要集中在列车-轨道各结构件的合理力学模型、系统耦合动力学方程的推导和高效数值计算方法以及应用虚拟样机的CAE技术建模来研究系统的性能动力学仿真计算，如ADAMs/Rail建模、SIMPACK建模、DADS、以及软件联合建模等。翟婉明[1]完成了《车辆—轨道耦合动力学》一书，该书从最初的垂向耦合动力学研究，到后来的横向耦合动力学研究，再到随机振动研究，逐渐形成了完整的理论体系。国内还有曾庆元[2]、陈果[3]、雷晓燕[4]等众多学者也从不同的角度研究了铁路车辆系统动力学问题，并取得了一系列丰富成果。在轨道车辆系统的可靠性分析中，刘伟渭等[5-6]主要研究了铁路车辆约束轮对的随机稳定和随机可靠性，研究内容忽略列车轨道耦合作用。王长一和张义民[7]进行了高速动车组转向架构架可靠性及可靠性灵敏度分析，该文用NASTRAN有限元和Monte-Carlo方法计算了转向架的可靠度。王洪礼等[8]研究了带有记忆合金弹簧的汽车半主动悬架系统在随机激励作用下，稳定性被打破后，系统发生的动力学行为以及相应需要采取的控制策略。含半主动悬架系统的车辆随机动力学和可靠性的研究还有文献[9-10]。目前同时考虑轨道车辆结构参数和轨道不平顺荷载随机性的系统动力可靠性分析研究成果较为鲜见。本文研究了轨道车辆系统在轨道不平顺多谐波和平稳高斯白噪声随机激励作用下，车轨悬架系统广义能量的随机响应过程统计特性和系统发生首次穿越失效的寿命分布问题，还考虑了车辆本身结构参数的随机变异性对系统的动力可靠性的影响，这样的分析更能反映系统的真实行为，同时也是满足理论发展的需要。1基于Hamilton原理的列车随机动力学模型研究表明，造成环境振动的主要原因是轨道车辆垂向荷载，而列车侧滚及横向振动荷载往往是忽略不计的。仅考虑列车系统的垂向动力响应，忽略碰撞及接触等非线性因素，并且假定钢轨是完全刚性的。通常车体在纵向和横向都是对称的，忽略轮轨之间的弹跳作用以及车体的摇摆和点头作用的影响[11]。建立轨道车辆系统随机几何模型见图1：图1车辆悬架系统垂向振动简化模型假定轮轨为绑定接触，建立系统Lagrange方程为：0))()(())(()()(0)()(12111221221112212222tWtryktrycyykyycymyykyycym（1）式中:)(tW为均值0，强度2D的高斯白噪声.)(tr为轨道不平顺多谐波激扰。式(1)的系统定义了一个具有弱阻尼随机激励的拟Hamilton系统。设系统广义坐标12211,yyqyq和广义动量12211,yypyp，式(1)整理变形如哈密顿正则方程为：))()()((2)()()(121211112111212212222222222121211112111212212111tWmktrmktrmcqmkpmcqmkpmcqmkpmcppqtWmktrmktrmcqmkpmcqmkpmcppq）（则系统的Hamilton函数(广义能量)可以表示为：22122221122221)(22)(2)(21)(21qmkmkqmkppH）（（3）2系统的随机动力响应统计特征分析根据拟不可积Hamilton系统的定义及性质，可知系统(2)依概率收敛到一维Itö扩散过程，以H表示这一极限扩散过程，则支配该随机过程的Itö随机微分方程为：)()()(tdBHHdH（4）其中，)(tB是标准Weiner过程，)(),(HH分别是随机过程的漂移系数与扩散系数。使用拟不可积Hamilton系统的随机平均法[12]，得到：212112122211212122222212211]/)2[()(1)(dpdpdqdqpHppHppHmpHmpHmHTHHpqqHppqqdpdpdqdqpHHTdpdpdqdqpHppHHTH),0,,(),,,()/1()(]/)[()(1)(221212121211212112122211(5)其中：1222221112111212212212121111,/,/)(,mcmcmpqmkpmcmpqmkpmcmmcm）（Dmk2121211，令,,11222212121mkwkmmmmw带入(5)计算可得漂移和扩散系数：0211212021021212121212)(]3432(8[2)(HmkwDwHHkcHmmcmmmkDwwHs）(6)其中),(min0cHHH为系统的初始能量值。3系统的动力可靠度分析车辆悬挂系统受随机参激不稳定和受轨道不平顺随机外激时，系统将在较大范围内作随机运动，系统广义能量状态停留在安全域内的概率就是可靠性。系统状态首次越出安全域就意味着损坏，系统状态首次穿越安全域边界的时间就是寿命。根据动态应力—强度干涉理论，动态可靠性可以定义为：CHHtHHtR),(Pr),(minXX（7）将系统广义能量随机过程),(tHX离散化为n个随机变量iH，其大小取为第i个时段的均值。这样就把设计期内系统广义能力状态响应函数离散为n个样本，随机能量载荷作用n次时系统的可靠度等价于n个载荷样本中最大载荷值所对应的可靠度。载荷样本中的最大载荷可定义为载荷作用n次时的等效载荷。由顺序统计量理论可知[13]，随机载荷作用n次时的最大能量载荷maxH实际上就是由载荷样本nHHH,...,,21所决定的最大顺序统计量)(nH。设单位时间内能量载荷出现一次，即nt，则),(),(),()()(),()]([)Pr()...Pr()Pr()(minmaxcmax21maxtRtRtHFHFtRzFzHzHzHzFHHtHnHXXXX（8）对于文中所考虑的悬挂系统，因其物理参数具有随机性，这将导致结构的动力可靠度亦具有随机性。利用求解随机变量函数数字特征的矩法[14]，可推得随机参数系统广义能量动力可靠度的均值和方差分别为：miiiHtHRHtHRxtRtR12),(),()]([]),([),(XXXXXXX（9）其中，m为结构物理随机变量的数目。全局灵敏度分析可以帮助了解影响系统功能各设计参数的相对重要程度，从而可以在设计和优化中优先或者重点考虑重要性程度高的设计参数或者忽略重要性程度低的设计参数，为系统工程设计和优化提供重要的指导。基于方差的重要性测度模型为：)()(RVarRVarSiiX(10)其中，Var
)(ixR为单独只在ix不确定作用下可靠度R的条件方差。4算例说明本文的轨道车辆系统算例参数参照下表1：表1：某轨道车辆相关参数名称单位均值方差车体质量m2t41.8/转向架质量m1t7.46/一系悬挂刚度k1MN/m1.050.05一系悬挂阻尼c1kN·s/m240.2二系悬挂刚度k2MN/m0.60.01二系悬挂阻尼c2kN·s/m100.1取地铁列车的速度v=60km/h，则多波简谐激扰位移输入模型为：ttttr44.209cos05.075.104cos15.047.10cos5.27.2)(计算时取安全域边界10,0,2.0cmin0HHH，得到不同白噪声强度下系统平均动力可靠性函数(图2)。00.20.40.60.811.21.41.61.8200.10.20.30.40.50.60.70.80.91D=0.2D=0.3D=0.4D=0.5D=0.2D=0.3D=0.4D=0.5H0=0.2tx102/无量纲可靠度函数R(t)图2不同噪声强度下可靠度函数由图中看出，随噪声强度系数D值的增加，可靠性函数)(tR减小越快；首次穿越时间概率密度峰值会变大，发生首次穿越的时间会提前;即系统可靠性能与轨道随机因素强弱有直接关系，轨道不平顺变化范围较大将使系统不易控制自身性能。本文四个物理参数2121,,,cckk基于方差的重要性测度为T]11.0,25.0,18.0,46