Application of Good Lattice Point with Power Generator Method in Stochastic Dynamic Analysis of Vehicle-Bridge System
-
摘要:
针对难以精确选取具有代表性的超高维随机相位角问题,采用方幂好格子点法生成代表性轨道不平顺样本,并将该样本作用于车-桥系统,得到轮轨力的均值与标准差;然后,通过对比虚拟激励法、确定性时程法和蒙特卡洛法的计算结果来探究该方法的计算精度与计算效率;最后,采用线性与非线性轮轨接触关系研究考虑列车日运营量的脱轨系数阈值. 以和谐号通过桥梁为例,计算结果表明:与蒙特卡洛法相比,采用方幂好格子点法生成不同方向的轨道不平顺样本之间的均匀性较好;方幂好格子点法求得的随机动力响应的概率特征参数与其他方法相比,具有较高的计算精度,其计算效率较蒙特卡洛法提高了近5倍;分别采用线性与非线性轮轨接触关系时所得的脱轨系数阈值相差达4.68%,方幂好格子点法具有较广泛适用性.
Abstract:Since it is difficult to accurately select representative super-high-dimensional random phase angles, the good lattice point with power generator method (GLPPGM) was utilized to generate samples of representative track irregularities, which were applied to the vehicle-bridge system to obtain the mean and standard deviation of random dynamic responses. Then, the calculation accuracy and efficiency of this method were explored by comparing the results of the pseudo-excitation method, deterministic time history method, and Monte Carlo method (MCM). Finally, the threshold value of the derailment factor considering the daily operation volume of trains was studied by using linear and nonlinear wheel-rail contact relationships. The Harmony train passing over a bridge was studied, and the results show that compared with that by the MCM, the uniformity between the samples of track irregularities generated by the GLPPGM in different directions is better. The probability characteristic parameters of the random dynamic response obtained by the GLPPGM have higher calculation accuracy than different methods, and its calculation efficiency is nearly five times higher than that of the MCM. When linear and nonlinear wheel-rail contact relationships are considered, the threshold value of the derailment factor differs by 4.68%, and the GLPPGM has a wider applicability.
-
钢轨波磨是铁路行业尚未解决的难题,在地铁线路上尤为严重. 任何一条地铁线路上,只要曲线轨道半径低于350 m,则该轨道内轨上出现波磨的几率接近100%;然而在大半径曲线轨道(R≥650~800 m)和直线轨道上,很少出现钢轨波磨:这是地铁线路上钢轨波磨的普遍规律[1-2]. 当列车通过波磨轨道时,轮轨界面上的剧烈振动加剧了车辆和轨道部件的疲劳断裂,甚至危害行车安全[3-4].
目前钢轨波磨的形成机理可分为3类:1) 轮轨表面初始不平顺激励下,轮轨界面的瞬时动力作用引起摩擦功波动导致钢轨波磨[5-10];2) 轮轨系统的粘-滑振动引起了钢轨波磨[11-13];3) 饱和轮轨蠕滑力激发了轮轨摩擦耦合系统的自激振动导致钢轨波磨[14-19]. 第一类理论模型为它激振动模型,该理论可以较好地解释地铁线路上的部分波磨现象,但不足以解释地铁线路上波磨主要出现在小半径曲线内轨,在大半径曲线和直线轨道上波磨发生率较低这一普遍现象;第二类理论模型预测得出的波磨频率约为20~80 Hz[12-13],远小于地铁线路上绝大多数波磨的通过频率;第三类波磨理论由陈光雄教授提出,该理论认为当车辆通过小半径曲线轨道时,导向轮对上的轮轨蠕滑力饱和,引发了轮轨系统的摩擦自激振动,致使钢轨表面出现波磨,该理论可以合理地解释地铁线路上规律性的波磨现象,正逐渐被学者们认可[20-23].
梯形轨枕因其具有良好的减振降噪效果被铺设在减振要求高的地段[24]. 然而,在小半径曲线梯形轨枕轨道上出现了严重的钢轨波磨[25],目前对于梯形轨枕轨道波磨成因的研究还不充分. 本文基于轮轨系统摩擦自激振动导致钢轨波磨的观点对小半径曲线上的梯形轨枕轨道波磨展开研究,建立了导向轮对-梯形轨枕轨道系统有限元模型. 该模型采用实体单元对扣件系统进行建模,考虑了弹条的扣压作用以及轨距挡块的横向限位作用,建模方式更加符合实际工况. 应用复特征值分析和瞬时动态分析分别求解了该摩擦耦合系统的运动稳定性和时域动态响应,揭示了小半径曲线上梯形轨枕轨道波磨的成因.
1. 数值模型和理论方法
1.1 导向轮对-梯形轨枕轨道系统有限元模型
当车辆通过小半径曲线梯形轨枕轨道时,导向轮对与钢轨的接触位置如图1所示. 外轮的轮缘与高轨的轨距面接触,内轮的踏面与低轨的轨顶面接触. 图1中:δL、δR分别为外轮、内轮与钢轨间的接触角,FL和NL分别为外轮与高轨接触界面上的横向蠕滑力和法向接触力,FR和NR分别为内轮与低轨接触界面上的横向蠕滑力和法向接触力,FSVL和FSVR为作用在导向轮对左、右轴端的垂向悬挂力,FSLL和FSLR为作用在导向轮对左、右轴端的横向悬挂力,KC为弹条刚度,KDL和CDL分别为横向缓冲垫的刚度和阻尼,KDV和CDV分别为垂向减振垫的刚度和阻尼.
扣件系统由轨下垫板、弹条、轨距挡块、铁垫板等组成,轨距挡块约束着钢轨的横向运动,弹条产生的扣压力将钢轨紧紧地压在轨道上,钢轨的纵向运动通过钢轨与轨下垫板间的摩擦力来约束,螺旋道钉将铁垫板固定在梯形轨枕上. 为了简化模型,本文仅对扣件系统中的轨下垫板和铁垫板进行实体单元建模.
依托于图1所示的轮轨系统接触模型,使用ABAQUS有限元分析软件建立了导向轮对-梯形轨枕轨道系统有限元模型,如图2所示. 车轮廓形为LM磨耗型踏面,轮对滚动圆半径为0.42 m;钢轨廓形为60 kg/m标准钢轨截面,模型中单根钢轨的长度为36 m;轨下垫板离散地分布在钢轨底部,相邻轨下垫板/铁垫板间的距离为0.625 m;梯形轨枕由横向钢管和预应力混凝土纵梁组合而成,相邻横向钢管间的距离为2.5 m;导向轮对、钢轨、轨下垫板、铁垫板和梯形轨枕通过实体单元建模,整个有限元模型采用C3D8I实体单元划分,模型中共包含
422968 个单元,606942 个节点. 轮轨间的接触以及钢轨和轨下垫板间的接触采用面-面接触算法和有限滑移准则进行建模. 轮轨间的摩擦系数为0.40,钢轨与轨下垫板间的摩擦系数为0.75[26]. 通过垂向压缩弹簧来模拟弹条的扣压作用,通过“Equation”约束模拟轨距挡块的横向限位作用,通过“Tie”约束模拟轨下垫板和铁垫板以及铁垫板和梯形轨枕的连接关系. 模型中采用接地弹簧和阻尼单元模拟侧向缓冲垫和垂向减振垫对梯形轨枕的缓冲减振作用,相邻缓冲减振垫间的距离为1.25 m. 缓冲减振垫的刚度为2 × 107 N/m,阻尼为7 × 104 N·s/m[27]. 各部件的材料参数见表1.表 1 材料参数Table 1. Material parameters部件 密度/
(×103 kg·m−3)弹性模量/MPa 泊松比 钢轨/轮对 7.80 210000 0.30 轨下垫板 1.19 12 0.45 铁垫板 7.80 173000 0.30 横向钢管 7.80 210000 0.30 混凝土纵梁 2.50 36000 0.20 1.2 轮轨摩擦耦合系统自激振动的有限元分析方法
采用复特征值分析来求解导向轮对-梯形轨枕轨道系统的运动稳定性. 在小半径曲线轨道段,导向轮对上的轮轨蠕滑力达到饱和,约等于滑动摩擦力. 在此情形下,轮轨摩擦耦合系统的运动方程为
M¨u+C˙u+Ku=0, (1) 式中:M、C、K分别为系统的质量、阻尼和刚度矩阵,u为节点位移向量.
式(1)的通解为[28]
u(t)=z∑i=1φieλit=z∑i=1φie(αi+jωi)t, (2) 式中:z为拟提取的特征值阶数;λi=αi+jωi,为第i阶特征值;φi为特征向量;t为时间.
由式(2)可以看出,当特征值实部为正时,节点位移随时间指数式增长,这意味着轮轨系统发生了不稳定振动. 通常采用等效阻尼比ξ来判断系统发生自激振动的趋势,如式(3).
ξ=−2αi|ωi|. (3) 当ξ<0时,轮轨摩擦耦合系统发生了自激振动. 而且,ξ越小,轮轨系统发生自激振动的趋势越强,钢轨波磨更易发生.
复特征值分析是一种频域分析方法,而瞬时动态分析是一种时域非线性分析方法,可以获取摩擦自激振动发生时系统的动态响应. 本文使用ABAQUS/Standard隐式求解器计算了轮轨系统的动态响应. 该求解器采用Newmark方法计算每个增量步结束时刻系统的平衡方程,进而获取节点的位移和速度[29],如式(4)、(5).
u|t+Δt=u|t+Δt˙u|t+Δt2((12−β)¨u|t+β¨u|t+Δt), (4) ˙u|t+Δt=˙u|t+Δt((1−γ)¨u|t+γ¨u|t+Δt), (5) 式中:Δt为积分步长,β = (1–α2)/4,γ = 0.5–α,−0.5≤α≤0.
2. 计算结果分析
2.1 小半径曲线段导向轮对-梯形轨枕轨道系统摩擦自激振动分析
复特征值分析的步骤为:
步骤1 在导向轮对轴端施加垂向悬挂力和横向悬挂力;
步骤2 引入轮轨间的摩擦耦合作用;
步骤3 提取轮轨系统的固有频率;
步骤4 进行复特征值分析,提取轮轨系统的摩擦自激振动.
通过定义轮对的横向滑动速度引入导向轮对和钢轨间的摩擦耦合作用. 应用复特征值分析提取了50~
1200 Hz内导向轮对-梯形轨枕轨道系统的摩擦自激振动,如图3所示. 可以看到,在该频率范围内系统共存在7个自激振动,对应频率分别为150.35、163.07、422.16、447.81、471.05、493.97 Hz和519.94 Hz.图4为7个自激振动的振动模态,自激振动主要发生在低轨和内轮上,与钢轨波磨发生在小半径曲线梯形轨枕轨道内轨上相吻合. 根据等效阻尼比越小,自激振动发生趋势越强的原则,可知频率为150.35 Hz和493.97 Hz的自激振动更易发生.
图5为曲线半径为350 m的梯形轨枕轨道上的内轨波磨,其波长λ为60~80 mm. 列车通过该轨道的平均速度v为37 km/h,由λ=v/f可得波磨通过频率为130~180 Hz[27]. 综上可得:当车辆通过小半径曲线段梯形轨枕轨道时,发生在导向轮对内侧轮轨界面上频率为150.35 Hz的自激振动致使内轨上出现了波长为60~80 mm的波磨.
2.2 轮轨法向接触力变化
应用瞬时动态分析求解轮轨系统的动态响应分2个步骤完成:第一步,在导向轮对的轴端施加垂向悬挂力和横向悬挂力;第二步,设置轮对运行速度为40 km/h,运行时间为0.15 s,采样间隔为5 × 10–5 s. 由于未施加任何外部激励且钢轨表面无初始不平顺,因此,计算所得的轮轨系统的响应为自激振动响应. 图6所示为轮轨界面上的法向接触力的时间历程曲线. 从图中可以看出,内侧轮轨界面上的法向接触力波动幅值要大于外侧轮轨界面上的,其中外轮与高轨间法向接触力的波动幅值约为20 kN,而内轮与低轨间法向接触力的波动幅值约为30 kN.
图7为轮轨法向接触力的功率谱密度分析结果. 可以看到,内外侧轮轨界面上法向接触力的振动主频为160 Hz. 由λ=v/f可得波长约为69 mm,该预测结果与线路上梯形轨枕轨道波磨的波长(60~80 mm)接近.
2.3 侧向缓冲垫阻尼对轮轨系统摩擦自激振动的影响
不同侧向缓冲垫阻尼CDL,轮轨系统自激振动频率的分布如图8所示. 从图中可以看出:轮轨系统自激振动的数量不随侧向缓冲垫阻尼变化,轮轨系统产生了7个自激振动;侧向缓冲垫阻尼对400~600 Hz内的自激振动的影响很小,而对150~200 Hz内的2个自激振动有较大影响,2个自激振动的等效阻尼比随侧向缓冲垫阻尼值的增大而增大,这说明系统发生自激振动的趋势减弱. 因此,增大侧向缓冲垫阻尼可以在一定程度上抑制梯形轨枕轨道上的内轨波磨.
2.4 垂向减振垫刚度对轮轨系统摩擦自激振动的影响
图9所示为不同垂向减振垫刚度KDV,轮轨系统自激振动频率的分布. 从图9可以看出,当垂向KDV=10 kN/mm时,轮轨系统有6个自激振动,最小等效阻尼比为–
0.00789 ;当KDV=20 kN/mm时,轮轨系统有7个自激振动,最小等效阻尼比为−0.00506 ;当KDV=30 kN/mm时,轮轨系统有6个自激振动,最小等效阻尼比为−0.00679 ;当KDV=40 kN/mm时,轮轨系统有6个自激振动,最小等效阻尼比为−0.00798 . 由以上分析可知,当KDV=20 kN/mm时,最小等效阻尼比高于其余3种工况,这意味着系统发生自激振动的趋势减弱. 过小或过大的垂向减振垫刚度均不利于抑制梯形轨枕轨道上的内轨波磨.2.5 梯形轨枕结构对轮轨系统摩擦自激振动的影响
图10所示为4种不同结构的梯形轨枕. 每块轨枕上横向钢管的数量N=2根时,对应结构的梯形轨枕相邻横向钢管间的距离为5.000 m;N=3根时,对应结构的梯形轨枕相邻横向钢管间的距离为2.500 m;N=5根时,对应结构的梯形轨枕相邻横向钢管间的距离为1.250 m;N=9根时,对应结构的梯形轨枕相邻横向钢管间的距离为0.625 m. 需要说明的是,横向钢管间距的设计值为2.500 m.
图11所示为4个不同结构的导向轮对-梯形轨枕轨道系统自激振动频率的分布. 从图中可以看出:梯形轨枕结构对轮轨系统的摩擦自激振动有显著影响,并非每块梯形轨枕上的横向钢管数量越多,轮轨系统就越稳定. 如蓝色椭圆圈内所示,当每块梯形轨枕上安装9根横向钢管时,轮轨系统产生自激振动的等效阻尼比反而更小,轮轨系统更不稳定,钢轨波磨更易出现. 红色椭圆圈内所示为低频段的2个自激振动,当每块梯形轨枕上安装2根横向钢管时,等效阻尼比更低,轮轨系统更易趋于不稳定. 当每块梯形轨枕上安装5根横向钢管时,等效阻尼比最大,轮轨系统发生自激振动的趋势减弱. 因此,在地铁小半径曲线地段,铺设横向钢管间距为1.250 m的梯形轨枕可以一定程度抑制梯形轨枕轨道波磨.
3. 结 论
本文研究了地铁小半径曲线梯形轨枕轨道上的内轨波磨,建立了导向轮对-梯形轨枕轨道系统有限元模型,模型中采用实体单元对扣件系统进行建模,应用复特征值分析和瞬时动态分析分别求解了轮轨系统的运动稳定性和时域动态响应. 研究了缓冲减振垫参数及梯形轨枕结构对轮轨系统摩擦自激振动的影响,得到以下结论:
1) 饱和轮轨蠕滑力引起的频率为150 Hz的轮轨系统共振振动是地铁线路小半径曲线段梯形轨枕轨道上波长为60~80 mm内轨波磨的成因.
2) 增大侧向缓冲垫阻尼可以一定程度上抑制梯形轨枕轨道上的内轨波磨,过小或过大的垂向减振垫刚度均不利于抑制钢轨波磨.
3) 改变单块梯形轨枕上横向钢管的数量实际上是改变了导向轮对–梯形轨枕轨道系统的刚度矩阵,进而影响了轮轨系统发生自激振动的趋势以及钢轨波磨的发展趋势. 铺设横向钢管间距为1.250 m结构的梯形轨枕可以抑制梯形轨枕轨道上的钢轨波磨.
-
表 1 车辆参数
Table 1. Train parameters
名称 拖车 动车 名称 拖车 动车 转向架轴距/m 2.5 2.5 x 向转向架转动惯量/(t•m2) 3.2 1.6 车辆定距/m 18 18 y 向转向架转动惯量/(t•m2) 7.2 1.7 一系悬挂横向跨距/m 0.33 0.34 z 向转向架转动惯量/(t•m2) 6.8 1.7 转向架导一系悬挂/m 2.05 2.05 一系纵向阻尼(单侧)/(kN•s•m−1) 960 960 二系悬挂横向跨距/m 2.05 2.05 一系横向阻尼(单侧)/(kN•s•m−1) 960 960 车体中心到二系悬挂/m 0.36 0.83 一系竖向阻尼(单侧)/(kN•s•m−1) 1040 700 二系悬挂到转向架/m 0.24 0.15 二系纵向阻尼(单侧)/(kN•s•m−1) 240 210 轮对滚动园半径/m 0.4575 0.4575 二系横向阻尼(单侧)/(kN•s•m−1) 240 210 轮对质量/t 1.9 2.2 二系竖向阻尼(单侧)/(kN•s•m−1) 400 350 轮对转动惯量/(t•m2) 1.067 1.63 一系纵向阻尼(单侧)/(kN•s•m−1) 0 0 转向架质量/t 3.4 1.7 一系横向阻尼(单侧)/(kN•s•m−1) 0 0 x 向车体转动惯量/(t•m2) 101.5 74 一系竖向阻尼(单侧)/(kN•s•m−1) 30 38 y 向车体转动惯量/(t•m2) 1064.4 1370 二系纵向阻尼(单侧)/(kN•s•m−1) 120 150 z 向车体转动惯量/(t•m2) 867.2 1370 二系横向阻尼(单侧)/(kN•s•m−1) 30 15 车体质量/t 42.4 44 二系竖向阻尼(单侧)/(kN•s•m−1) 33 40 -
[1] 雷虎军,黄江泽. 考虑桩土相互作用的车-轨-桥系统地震响应分析[J]. 西南交通大学学报,2021,56(2): 229-237.LEI Hujun, HUANG Jiangze. Seismic responses analysis of train-track-bridge system considering pile-soil interaction[J]. Journal of Southwest Jiaotong University, 2021, 56(2): 229-237. [2] 何佳骏,向活跃,龙俊廷,等. 桥隧过渡段高速列车行车抗风安全分析[J]. 西南交通大学学报,2021,56(5): 1056-1064.HE Jiajun, XIANG Huoyue, LONG Junting, et al. Wind-resistant safety analysis of high-speed trains passing through bridge-tunnel transition[J]. Journal of Southwest Jiaotong University, 2021, 56(5): 1056-1064. [3] 王伟,张亚辉,欧阳华江. 考虑轨道不平顺随机性的车轨耦合系统动力可靠度分析[J]. 大连理工大学学报,2019,59(2): 111-117.WANG Wei, ZHANG Yahui, OUYANG Huajiang. Dynamic reliability analysis of vehicle-track coupled systems considering randomness of track irregularities[J]. Journal of Dalian University of Technology, 2019, 59(2): 111-117. [4] 林家浩,张亚辉,赵岩. 虚拟激励法在国内外工程界的应用回顾与展望[J]. 应用数学和力学,2017,38(1): 1-31. doi: 10.1007/s10483-016-2152-6LIN Jiahao, ZHANG Yahui, ZHAO Yan. The pseudo-excitation method and its industrial applications in China and abroad[J]. Applied Mathematics and Mechanics, 2017, 38(1): 1-31. doi: 10.1007/s10483-016-2152-6 [5] XU L, ZHAI W M. Probabilistic assessment of railway vehicle-curved track systems considering track random irregularities[J]. Vehicle System Dynamics, 2018, 56(10): 1552-1576. doi: 10.1080/00423114.2018.1424916 [6] XIA C Y, XIA H, DE ROECK G. Dynamic response of a train-bridge system under collision loads and running safety evaluation of high-speed trains[J]. Computers & Structures, 2014, 140: 23-38. [7] 朱艳,李小珍,强士中. 高速铁路简支梁桥车桥系统随机响应[J]. 西南交通大学学报,2011,46(4): 535-540. doi: 10.3969/j.issn.0258-2724.2011.04.002ZHU Yan, LI Xiaozhen, QIANG Shizhong. Stochastic responses of train-bridge coupling system for simply-supported high-speed railway bridge[J]. Journal of Southwest Jiaotong University, 2011, 46(4): 535-540. doi: 10.3969/j.issn.0258-2724.2011.04.002 [8] HE X H, SHI K, WU T. An efficient analysis framework for high-speed train-bridge coupled vibration under non-stationary winds[J]. Structure and Infrastructure Engineering, 2020, 16(9): 1326-1346. doi: 10.1080/15732479.2019.1704800 [9] 赵岩,张亚辉,林家浩. 列车-结构耦合作用的随机振动分析方法[M]. 北京:科学出版社,2018:156-190. [10] ZHANG N, ZHOU Z J, WU Z Z. Safety evaluation of a vehicle-bridge interaction system using the pseudo-excitation method[J]. Railway Engineering Science, 2022, 30(1): 41-56. doi: 10.1007/s40534-021-00259-6 [11] 李杰,陈建兵. 随机动力系统中的概率密度演化方程及其研究进展[J]. 力学进展,2010,40(2): 170-188. doi: 10.6052/1000-0992-2010-2-J2009-105LI Jie, CHEN Jianbing. Advances in the research on probability density evolution equations of stochastic dynamical systems[J]. Advances in Mechanics, 2010, 40(2): 170-188. doi: 10.6052/1000-0992-2010-2-J2009-105 [12] YU Z W, MAO J F. A stochastic dynamic model of train-track-bridge coupled system based on probability density evolution method[J]. Applied Mathematical Modelling, 2018, 59: 205-232. doi: 10.1016/j.apm.2018.01.038 [13] XU L, ZHAI W M, GAO J M. A probabilistic model for track random irregularities in vehicle/track coupled dynamics[J]. Applied Mathematical Modelling, 2017, 51: 145-158. doi: 10.1016/j.apm.2017.06.027 [14] 刘章军,张传勇,齐东春. 高速列车乘坐舒适度评价的概率密度演化方法[J]. 振动与冲击,2018,37(14): 149-155.LIU Zhangjun, ZHANG Chuanyong, QI Dongchun. Probability density evolution method for the ride comfort evaluation of high speed trains[J]. Journal of Vibration and Shock, 2018, 37(14): 149-155. [15] LIU Z J, LIU Z H, PENG Y B. Simulation of multivariate stationary stochastic processes using dimension-reduction representation methods[J]. Journal of Sound and Vibration, 2018, 418: 144-162. doi: 10.1016/j.jsv.2017.12.029 [16] ZHANG N, XIA H. Dynamic analysis of coupled vehicle–bridge system based on inter-system iteration method[J]. Computers & Structures, 2013, 114/115: 26-34. [17] 张楠,郭薇薇,夏禾. 高速铁路车桥耦合动力学[M]. 北京:北京交通大学出版社,2018:171-177. [18] 方开泰,刘民千,覃红,等. 均匀试验设计的理论和应用[M]. 北京:北京科学出版社,2018. [19] MARSAGLIA G. Ratios of normal variables[J]. Journal of Statistical Software, 2006, 16(4): 1-10. [20] ZHU Y, LI X Z, JIN Z B. Three-dimensional random vibrations of a high-speed-train–bridge time-varying system with track irregularities[J]. Proceedings of the Institution of Mechanical Engineers Part F Journal of Rail and Rapid Transit, 2015, 230(8): 1851-1876. [21] 翟婉明. 车辆-轨道耦合动力学(上册)[M]. 4版. 北京:科学出版社,2015:120-128. -