DOI： 10.11817/j.issn.1672-7207.2021.02.030

铁道车辆车轮跳轨全过程计算方法

向俊¹，陈林¹，苏玮¹，杨海明¹，龚凯²，彭子祥¹

(1. 中南大学 土木工程学院，湖南 长沙，410075；

2. 华东交通大学 土木建筑学院，江西 南昌，330013)

摘要：为研究提速及高速列车跳轨脱轨机理，提出跳轨车轮位移及轮轨力的基本特征，强调考虑轮轨冲击作用的必要性。在具体计算时，应包含3种完全不同的动力学控制方程，分别是轮轨密贴时的方程、轮轨分离时的方程和轮轨冲击时的方程。研究结果表明：本文方法从本质上揭示了列车车轮跳轨全过程，能够描述车轮跳轨的主要信息；速度越大，车轮跳轨高度越高，跳轨持续时间越长；是否考虑轮轨冲击作用，对于车体加速度、钢轨速度、钢轨加速度、轨枕速度及轨枕加速度等车-轨系统振动响应影响较大。

关键词：铁道工程；行车安全；跳轨全过程；轮轨冲击；跳轨高度；跳轨时间

中图分类号：U211.5            文献标志码：A              开放科学(资源服务)标识码(OSID)

文章编号：1672-7207（2021）02-0614-11

A method of calculating the whole process of wheel-jumping of railway vehicle

XIANG　Jun¹, CHEN　Lin¹, SU　Wei¹, YANG　Haiming¹, GONG　Kai², PENG　Zixiang¹

(1. School of Civil Engineering, Central South University, Changsha 410075, China;

2. School of Civil Engineering and Architecture, East China Jiaotong University, Nanchang 330013, China)

Abstract: To study the mechanism of train derailment on speed-up and high-speed railway, the basic characteristics of wheel displacement and wheel-rail interaction force during the wheel-jumping process were presented, and the necessity of considering wheel-rail impact was emphasized.In actural calculation, three completely different dynamic governing equations namely the equations of contact, separation and wheel-rail impact stage, should be included. The results show that this developed method essentially reveals the entire process of wheel-jumping of vehicle, and can obtain valuable information of wheel-jumping. The higher the vehicle speed is, the larger the wheel-jumping height, and whether to consider the effect of wheel-rail impact has a great influence on the vibration response of the vehicle-track system such as the acceleration of car body, velocity and acceleration of rail and sleeper.

Key words: railway engineering; running safety; whole process of wheel-jumping; wheel-rail impact; wheel-jumping height; wheel-jumping duration

列车脱轨问题一直是铁路工作者的研究重点与难点，此问题很复杂，导致列车脱轨问题一直未能得到有效解决。曾庆元等^[1]突破传统的研究方法，提出了列车脱轨能量随机分析方法，并在最近10几年来，于高速与重载铁路上得到应用。龚凯等^[2-4]应用此方法，研究了诸多因素对货物列车运行安全性的影响规律，并提出了一些良好的工程措施。YU等^[5-6]应用此方法，研究了无砟轨道诸多病害对高速列车运行安全性的影响规律，并提出了相应措施及无砟轨道维修标准建议方案。尽管如此，此方法也只是针对爬轨脱轨进行了研究，并未对跳轨脱轨进行研究。试验结果表明，随着列车速度不断提高，列车车轮更易发生跳轨(即车轮与钢轨分离)现象，严重时将导致列车脱轨事故发生，这对列车运行安全性直接构成威胁，因此，需进一步研究考虑跳轨脱轨的列车脱轨能量随机分析方法。FR?BA等^[7]研究了车轮跳轨现象及二轴车辆跳离钢轨的车桥振动问题，指出了考虑跳轨现象的必要性。LEE^[8-9]采用车轮与钢轨刚性接触模型，研究了单个车轮跳离钢轨的车桥振动问题。LIU等^[10-11]研究了车辆在桥梁上的跳轨问题。ST?NCIOIU等^[12-13]研究了单轴及两轴车辆模型在简支梁上的跳轨问题。CHENG等^[14]采用车轮与钢轨刚性接触模型，研究了单个轮对跳离钢轨的车桥振动问题。BAEZA等^[15]研究了单轴车辆在简支桁梁上的跳轨问题。翟婉明^[16]建立了单轮对-轨道相互作用模型，对单轮对爬轨和跳轨两种工况下的脱轨行为进行了计算分析。娄平^[17]基于轮轨刚性接触模型，研究了二轴车辆跳离桥梁结构的振动问题，模拟了车轮与钢轨接触、跳离、再次接触的动态过程。曾京等^[18-19]分析了轮对横向冲击对车轮跳轨的影响规律。肖新标等^[20-21]研究了复杂环境状态下的轮轨分离现象，建立了判断轮轨是否分离的指标。孙丽霞等^[22]研究了车辆蛇形运动对于跳轨脱轨的影响。XU等^[23]在空间车轨模型中建立了时变耦合矩阵，简化了轮轨分离时的振动方程的形式。ZHANG等^[24]通过考虑分离和路面不平顺的车桥相互作用的非线性多弹簧模型发现，无论是否从桥上分离，其耦合的控制方程式都不会改变。JU^[25]建立了一种非线性轮轨接触单元模型，可以模拟轮轨密贴、滑动和分离。由此可见，国内外对铁道车辆车轮跳轨全过程计算方面的研究太少，并未取得令人满意的结果。文献[7-12, 24]没有考虑轨道振动及轮轨冲击作用。文献[13-15]虽然考虑了轮轨冲击作用，但是没有考虑轨道振动。文献[16, 17, 23]虽然考虑了轨道振动，但在轮轨再次密贴时，没有考虑轮轨冲击作用。文献[18-22]主要考虑了轮轨横向振动与冲击对跳轨脱轨的影响。

本文提出一种铁道车辆车轮跳轨全过程的计算方法，同时考虑轨道振动与轮轨冲击作用，可以反映车轮跳轨全过程的完整信息，以便为进一步研究考虑跳轨脱轨的列车脱轨能量随机分析方法打下良好基础。

1  车轮跳轨的基本特征

跳轨车轮的位移及轮轨力可以直接表征车轮跳轨的基本特征，如图1所示。

图1　时程曲线示意图

Fig. 1　Schematic diagram of time histories

跳轨车轮的位移基本特征表现为车轮经历了从轮轨密贴—车轮上升、车轮下降—轮轨密贴的全过程。图1(a)所示为跳轨车轮的跳轨高度时程曲线。从图1(a)可见：在列车正常行驶过程中，所有车轮与钢轨之间保持密贴关系；当列车中的某车轮在跳轨诱因作用下发生轮轨分离时，车轮就会从开始起跳而上升，直到上升到最高点，车轮跳轨高度达到最大值，之后开始下降，直到落在钢轨上，又维持了轮轨密贴状态(这里假定车轮没有发生跳轨脱轨)。

跳轨车轮的轮轨力基本特征表现为4个阶段，即轮轨密贴、轮轨分离、轮轨冲击及轮轨再次密贴。图1(b)所示为跳轨车轮的轮轨力时程曲线。从图1(b)可见：在列车正常行驶过程中，车轮与钢轨间的相互作用力处于轮轨密贴状态下的幅值水平；当列车中的某车轮在跳轨诱因作用下发生轮轨分离时，从至，轮轨力保持为0；在时，轮轨本应保持紧贴状态，但由于跳轨车轮与钢轨之间存在相对速度，从而产生轮轨冲击力，并在时达到最大值(其中，为小于1的正数，并在计算时根据精度要求加以确定)。在经历了冲击时间后，又恢复到轮轨密贴状态下的幅值水平。

由图1(a)可见：从至计算终点时刻，轮轨保持密贴状态。由图1(b)可见：从至，虽然轮轨处于密贴状态，但轮轨之间产生了冲击力，之后才恢复到正常密贴状态下的幅值水平。根据跳轨车轮的轮轨力基本特征，可将车轮跳轨全过程划分为轮轨密贴(第Ⅰ阶段)、轮轨分离(第Ⅱ阶段)、轮轨冲击(第Ⅲ阶段)及轮轨再次密贴(第Ⅳ阶段)共4个阶段。

由此可见，除了第Ⅰ与Ⅳ阶段具有共同动力学特征外，其他各阶段具有完全不同的动力学特征，尤其体现在它们的动力学控制方程之中。这样，在计算列车车轮跳轨全过程时，必须分别建立3种不同的动力学控制方程，分别为轮轨密贴时方程、轮轨分离时方程和轮轨冲击时方程。

2  计算方法

2.1　车辆振动分析模型

采用二轴车辆竖向振动分析模型，如图2所示。将车辆视为由1个车体、2个轮对以及车体和轮对之间的悬挂组成的多刚体系统，其中，车体视为质量为和绕其质心转动惯量为的刚体，其自由度为质心处的竖向位移和转角；车体与轮对之间的悬挂弹簧刚度与阻尼系数分别为和；每个轮对质量为，后、前轮对的竖向位移分别为和；为车辆定距之一半；v为车速。假定车辆向下的位移和顺时针转角为正向，其位移从各自静平衡位置处开始测量。

图2　二轴车辆及轨道振动分析模型

Fig. 2　Two-axle vehicle and track vibration analysis model

2.2　轨道振动分析模型

考虑两层支承的轨道竖向振动分析模型如图2所示，其中，钢轨视为弹性点支承的Bernoulli-Euler梁，总长为，并将相邻两支承点之间划分为1个梁单元，单元长为；轨枕视为质量为的单自由度刚体，其竖向位移为；钢轨与轨枕之间的线性弹簧刚度和黏滞阻尼系数分别为和；道床的线性弹簧刚度和黏滞阻尼系数分别为和。从轨道左端到轮对的水平距离为(=1，2)，钢轨表面的竖向不平顺为。假定轨道向下的位移为正，转角顺时针方向为正向。

轨道节点与轨枕位移从左端依次排序，钢轨单元的形函数采用三次Hermite插值函数，如果以表示单元形函数矩阵，表示单元节点位移矢量，表示轮对与所在钢轨单元左节点之间的距离，那么，单元内任意一点在时刻的竖向位移可表示为

(1)

式中：；；；；。

2.3　振动方程的建立及其求解

运用弹性系统动力学总势能不变值原理及形成矩阵的“对号入座”法则^[1]，建立车辆-轨道系统振动方程。采用Wilson-θ逐步积分法，基于Fortran编制程序求解系统振动方程。

2.3.1　轮轨密贴阶段的控制方程　

轮轨密贴时，轮对的位移不独立，后、前轮对的竖向位移和可表示为

；(2)

轮轨密贴时的振动方程可表示为

(3)

式中，各分矩阵表达如下。

1) 车辆的位移列阵、钢轨的位移列阵及轨枕的位移列阵分别表示为：

(4)

(5)

(6)

其中：中奇数项为节点位移，偶数项为节点转角；为钢轨总自由度数；为轨枕总自由度数。

车辆、钢轨及轨枕的速度列阵分别为，和，车辆、钢轨及轨枕的加速度列阵分别为，和。

2) 车辆质量矩阵可表示为

(7)

钢轨质量矩阵可表示为

(8)

其中：；

；。

钢轨本身质量矩阵由钢轨单元质量矩阵(阶数为4×4)组集而成，为钢轨单位长度质量；受轮对影响引发的钢轨质量矩阵中的表示轮对所在钢轨单元的形函数矩阵；为轮对与所在钢轨单元左节点之间的距离；下标，分别对应后轮对和前轮对；除轮对所在钢轨单元对应的4个元素外，其余元素为0。

轨枕质量矩阵可表示为

(9)

3) 车辆刚度矩阵可表示为

(10)

钢轨刚度矩阵可表示为

(11)

其中：；；

钢轨本身刚度矩阵由钢轨单元刚度矩阵(阶数为4×4)组集而成，和分别为钢轨弹性模量和钢轨截面对水平轴的惯性矩，和分别为对局部坐标的一次导数和二次导数；和()分别为车辆和扣件刚度引起的钢轨刚度矩阵，后者对角线上除第1个和最后1个奇数外，每个奇数位置对应的元素为。

车辆与钢轨相互作用刚度矩阵和的表达式为：

(12)

(13)

其中：

；

；

;

;

和分别为后轮对与钢轨相互作用刚度矩阵；和分别为前轮对与钢轨相互作用刚度矩阵。

轨枕的刚度矩阵可表示为

  (14)

钢轨与轨枕相互作用刚度矩阵为与。，刚度矩阵中除第行、第列外，位置对应的元素为，其余均为0。

4) 车辆阻尼矩阵可用代替中的获得。

钢轨阻尼矩阵可表示为

(15)

其中：；和分别为车辆和扣件阻尼引起的钢轨刚度矩阵，后者对角线上除第1个和最后1个奇数外，每个奇数位置对应的元素为。

车辆与钢轨相互作用阻尼矩阵和的表达式为：

(16)

(17)

其中：；

;

和分别为后轮对与钢轨相互作用阻尼矩阵；和分别为前轮对与钢轨相互作用阻尼矩阵。

轨枕的阻尼矩阵可表示为

(18)

钢轨与轨枕相互作用阻尼矩阵为和。，用代替中的可得。

5) 车辆的荷载列阵可表示为

(19)

钢轨的荷载列阵可表示为

(20)

其中：

；

；轨枕的荷载列阵元素全为0。

2.3.2　轮轨分离阶段的控制方程　

对于车轮跳轨阶段，车轮与钢轨已经分离，共有以下3种情况：前轮对跳轨，后轮对密贴；后轮对跳轨，前轮对密贴；前、后轮对同时跳轨。以下介绍不同情况下的轮轨分离阶段的车辆-轨道系统振动控制方程。

1) 对于一个轮对跳轨、另一个轮对密贴的情况，其振动方程相似，这里仅以前轮对单独跳离钢轨为例介绍。当前轮对跳离钢轨且后轮对密贴时，前轮对的位移独立，后轮对的竖向位移不独立，系统振动方程在密贴的基础上发生变化，可表示为

(21)

式中，车辆位移列阵、质量矩阵和刚度矩阵可表示为：

(22)

(23)

(24)

注意，将中的替换为，即得车辆振动阻尼矩阵。

钢轨质量矩阵、刚度矩阵和阻尼矩阵可表示为

(25)

其中：；

(26)

;

(27)

。

车辆与钢轨相互作用刚度矩阵和及阻尼矩阵和可表示为：

(28)

(29)

(30)

车辆荷载列阵和钢轨荷载列阵可表示为

(31)

(32)

2) 当前轮对与后轮对同时跳离钢轨时，前、后轮对的竖向位移和均独立，车辆系统振动方程可表示为

(33)

式中：车辆位移列阵、质量矩阵、刚度矩阵及车辆荷载列阵可表示为：

(34)

(35)

(36)

(37)

同时，车辆振动阻尼矩阵可通过将中的替换为获得。

轨道振动方程可表示为

(38)

钢轨质量矩阵、刚度矩阵和阻尼矩阵可表示为：；；。

2.3.3　轮轨冲击阶段的控制方程　

由图1(a)可见：在时，轮轨从分离阶段进入密贴阶段，此时，车轮位移等于钢轨位移与轨道不平顺之和，轮轨相对位移量值为零。然而，此时，并不能保证轮轨相对速度也为零，这意味着轮轨之间会出现速度差值，必然引起轮轨冲击，如图1(b)所示。这里仅以前轮对冲击、后轮对密贴为例来介绍轮轨冲击阶段的计算方法，此时，振动方程可表示为

(39)

车辆荷载列阵和钢轨荷载列阵的表达式如下：

(40)

(41)

其中：为冲击力。

鉴于轮轨冲击力的计算比较复杂，本文把它视为优化问题并采用迭代算法完成计算，具体计算步骤如下。

第一步：假定冲击力函数如图1(b)所示，具体表达式为

(42)

式中：和为假定的作用力函数，可以是线性的，也可以是非线性的；为最大冲击力幅值；为冲击终端时的幅值，其特点是应与后续轮轨正常密贴时的轮轨力幅值相当；为跳轨起点；为冲击起点；为冲击时间；为冲击终点；为小于1的正数。

第二步：假定控制变量为，和。

第三步：建立目标函数，并以时刻的轮轨相对速度表示。建立目标函数如下：

(43)

式中：为车轮速度；为钢轨随时间变化的速度；为钢轨随距离变化的速度；为时刻轮轨接触点距离轨道计算起点的距离；为轨道不平顺随距离变化的速度。

第四步：分别建立冲击终点时刻的轮轨相对位移、轮轨相对速度及轮轨力的约束条件^[14]，即轮轨相对位移为零或为很小的数值、轮轨相对速度为零或为很小的数值、轮轨力与轮轨正常密贴时的轮轨力相当。

轮轨相对位移约束条件如下：

(44)

当，，时，

(45)

当，，时，

(46)

当，，时，

(47)

关于轮轨相对速度约束条件如下：

(48)

(49)

另外，轮轨力约束条件如下：

(50)

式(44)~(50)中的为较小的正数，用以控制计算精度；为前轮对实际轮轨力，可表达为

 (51)

第五步：满足所有的约束条件且使目标函数达到最小值时的3个控制变量，和即为所求，从而最终可以确定轮轨冲击力。

3  计算实例及其分析

二轴车辆模型参数如表1所示，轨道模型参数如表2所示。选取波长为0.8 m、波幅a为2 mm的轨面单一谐波型不平顺作为跳轨诱因，并且设置于轨道中间部位，具体表达式如下：

(52)

轮轨冲击力函数假定如下：

(53)

式中：取0.5。同时，用以控制计算精度的参数取值为，。

取轨道计算长度=72 m，分别计算车辆以速度80，90和100 km/h通过轨道时的车轮跳轨全过程。为了符合习惯，图3(a)和图4中符号取向上为正，其他图中的符号与前述模型图中的坐标方向一致。

表1　二轴车辆模型参数

Table 1　Two-axle vehicle model parameters

表2　轨道模型参数

Table 2　Track model parameters

图3(a)所示为当车辆以速度90 km/h通过该轨道不平顺时，前轮对车轮跳轨高度时程曲线。由图3(a)可知：当车辆运行至=0.403 s时，轮轨开始分离，车轮开始起跳，并且一直上升；当=0.408 s时，车轮跳轨高度达到最大值=5.52 mm，之后，车轮开始下降；当=0.412 s时，车轮回落到钢轨之上。由此可见，车轮跳轨高度时程曲线，反映了跳轨车轮位移的基本特征，即车轮开始起跳、车轮上升、车轮下降、轮轨再次密贴。

图3(b)所示为当车辆以速度90 km/h通过式(52)表示的轨道不平顺时，前轮对轮轨力时程曲线。由图3(b)可知：当车辆运行至=0.403 s时，轮轨力为零，表示轮轨开始分离；当=0.412 s时，车轮经过跳轨阶段后重新回落到钢轨之上，但是，由于此时车轮竖向速度=0.411 m/s，钢轨竖向速度=0.004 m/s，二者之间出现了相对速度，必然引起轮轨冲击作用，冲击时间=9.0 ms，冲击过程中的最大值=386.2 kN；当=0.421 s时，轮轨冲击结束，重新开始正常状态下的轮轨密贴。由此可见，轮轨力时程曲线反映了跳轨车轮轮轨力的基本特征，即轮轨密贴、轮轨分离、轮轨冲击及轮轨再次密贴。特别地，从至，虽然轮轨处于密贴阶段，但轮轨之间存在冲击作用，在列车车轮跳轨全过程计算中，必须考虑到这一因素，否则，对时刻之后的车-轨系统振动响应将会失真。

图3　车速为90 km/h时的时程曲线示意图

Fig. 3　Schematic diagram of time histories when vehicle velocity is 90 km/h

由此可见，本文提出的计算方法从本质上揭示了列车车轮跳轨全过程，能够描述车轮跳轨的关键信息，如跳轨高度、跳轨时间、跳轨姿态、轮轨冲击力幅值、冲击时间等。

图4所示为车辆分别以速度80，90，100 km/h通过轨道不平顺(见式(52))时，前轮对跳轨高度时程曲线，其中，与分别表示车速80 km/h时跳轨时间起点与终点，与分别表示90 km/h时的跳轨时间起点与终点，与分别表示车速100 km/h时的跳轨时间起点与终点，，与分别表示车速80，90和100 km/h时的跳轨高度最大值；同一车速下的跳轨终点时刻与起点时刻的差值表示车轮跳轨持续时间，其中，=0.363 s，=0.374 s，=11 ms；=0.855 mm，=0.403 s，=0.412 s，=9 ms；=0.552 m，=0.454 s，=0.459 s，=5 ms，=0.075 mm。>>，即速度越高，车轮跳轨持续时间越长；>>，即速度越高，车轮跳轨高度越高。

图4　不同速度条件下的车轮跳轨高度时程曲线

Fig. 4　Time history of wheel-jumping height under different speeds

为了探究是否考虑轮轨冲击作用对于车-轨系统振动响应的影响程度，以车速90 km/h为例，当车辆通过式(52)表示的轨道不平顺时，分别考虑轮轨冲击和不考虑轮轨冲击，计算车体加速度、钢轨速度、钢轨加速度、轨枕速度和轨枕加速度等振动响应时程曲线，如图5~7所示，其中，红色实线表示考虑轮轨冲击作用时的结果，蓝色虚线表示不考虑冲击时的结果。

图5所示为车体加速度时程曲线局部图。由图5可见：在=0.412 s至=0.421 s期间，考虑冲击时的车体加速度最大幅值为4.25 ，不考虑冲击时的车体加速度最大幅值为0.71 ，并且考虑冲击时的车体加速度曲线幅值普遍远大于不考虑冲击时的车体加速度曲线幅值，同时，考虑冲击时的车体加速度振动频率也远高于不考虑冲击时的车体加速度振动频率。可见，是否考虑轮轨冲击对于车体振动加速度影响很大。

图5　=90 km/h、考虑与不考虑冲击时的车体加速度时程曲线

Fig. 5　Time histories of car body acceleration considering wheel-rail impact or not when =90 km/h

图6(a)所示为前轮对作用点下的钢轨速度时程曲线局部图。由图6(a)可见：从=0.412 s开始，是否考虑轮轨冲击作用导致钢轨速度时程曲线发生了很大差异，主要表现在振动幅值和相位，其中，考虑冲击时的钢轨最大速度为0.410 m/s，不考虑冲击时的钢轨速度为0.225 m/s，前者是后者的约1.82倍；另外，二者曲线的相位几乎为反相位状态。

图6(b)所示为前轮对作用点下的钢轨加速度时程曲线局部图。由图6(b)可见：从=0.412 s开始，考虑冲击时钢轨最大加速度为-362.188 ，不考虑冲击时的钢轨最大加速度为78.205 ，前者是后者的约4.63倍，同时，在=0.421 s之后的二者曲线相位几乎呈反相位状态。

由此可见，图6充分表明是否考虑轮轨冲击，对于钢轨速度及加速度影响甚大。

图6　=90 km/h、考虑冲击与不考虑冲击前轮对作用点下钢轨的时程曲线

Fig. 6　Time history of rail under action point of front wheelset considering impact or not when =90 km/h

图7(a)所示为计算长度中点处的轨枕速度时程曲线局部图。由图7(a)可见：从=0.412 s开始，考虑冲击时轨枕最大速度为0.113 m/s，不考虑冲击时的轨枕最大速度为0.043 m/s，前者是后者的约2.63倍；同时，在=0.421 s之后的二者曲线相位几乎为反相位状态。

图7(b)所示为计算长度中点处的轨枕加速度时程曲线局部图。，由图7(b)可见：从=0.412 s开始，考虑冲击时钢轨最大加速度为-97.490 ，不考虑冲击时的钢轨最大加速度为18.045 ，前者是后者的约5.40倍；同时，在=0.421 s之后的二者曲线相位也有了很大差异。

由此可见，图7也充分表明是否考虑轮轨冲击对于轨枕速度及加速度影响很大。

图7　=90 km/h、考虑冲击与不考虑冲击计算长度中点处轨枕的时程曲线

Fig. 7　Time histories of sleeper at midpoint of calculated length considering impact or not when =90 km/h

4  结论

1) 在提出了跳轨车轮位移及轮轨力基本特征的基础上，建立了列车车轮跳轨全过程计算方法。此方法能够描述车轮跳轨的主要信息。

2) 随着车速增大，车轮跳轨高度越高，跳轨持续时间越长。

3) 是否考虑轮轨冲击作用对于车体加速度、钢轨速度、钢轨加速度、轨枕速度、轨枕加速度及轮轨力等车-轨系统振动响应影响较大，因此，计算列车车轮跳轨时，必须考虑轮轨冲击作用。

参考文献：

[1] 曾庆元, 向俊. 列车脱轨分析理论与应用[M]. 长沙: 中南大学出版社, 2006: 125-143.

ZENG Qingyuan, XIANG Jun. The theory and application of train derailment analysis[M]. Changsha: Central South University Press, 2006: 125-143.

[2] 龚凯, 向俊, 毛建红, 等. 洪涝灾害引起的货物列车脱轨全过程分析[J]. 中南大学学报(自然科学版), 2015, 46(10): 3954-3960.

GONG Kai, XIANG Jun, MAO Jianhong, et al. Analysis of derailment course of freight train in floods[J]. Journal of Central South University(Science and Technology), 2015, 46(10): 3954-3960.

[3] 龚凯, 向俊, 毛建红, 等. 有砟及无砟轨道结构对货物列车运行安全性的影响[J]. 中南大学学报(自然科学版), 2017, 48(8): 2152-2161.

GONG Kai, XIANG Jun, MAO Jianhong, et al. Influence on freight train operation safety for ballasted and ballastless track[J]. Journal of Central South University(Science and Technology), 2017, 48(8): 2152-2161.

[4] 龚凯, 向俊, 毛建红, 等. 重载铁路桥上货物列车脱轨预防措施[J]. 中南大学学报(自然科学版), 2017, 48(12): 3406-3414.

GONG Kai, XIANG Jun, MAO Jianhong, et al. Derailment precautions of freight train on bridge in heavy haul railway[J]Journal of Central South University(Science and Technology), 2017, 48(12): 3406-3414.

[5] YU Cuiying, XIANG Jun, MAO Jianhong, et al. Influence of slab arch imperfection of double-block ballastless track system on vibration response of high-speed train[J]. Journal of the Brazilian Society of Mechanical Sciences and Engineering, 2018, 40(2): 109-123.

[6] 余翠英, 向俊, 林士财, 等. 双块式无砟轨道离缝对高速行车安全性影响及维修标准研究[J]. 铁道科学与工程学报, 2019, 16(8): 1865-1874.

YU Cuiying, XIANG Jun, LIN Shicai, et al. Influence of double-block ballastless track seam on high-speed train operation safety and maintenance standard research[J]. Journal of Railway Science and Engineering, 2019, 16(8): 1865-1874.

[7] FR?BA L. Impacts of two-axle system traversing a beam[J]. International Journal of Solids and Structures, 1968, 4: 1107-1123.

[8] LEE U. Revisiting the moving mass problem: onset of separation between the mass and beam[J]. Journal of Vibration and Acoustics, 1996, 118(3): 516-521.

[9] LEE U. Separation between the flexible structure and the moving mass sliding on it[J]. Journal of Sound and Vibration,1998, 209(5): 867-877.

[10] LIU Xuewen, XIE Jun, WU Chuan, et al. Semi-analytical solution of vehicle-bridge interaction on transient jump of wheel[J]. Engineering Structures, 2008, 30(9): 2401-2412.

[11] ZHU D Y, ZHANG Y H, OUYANG H. A linear complementarity method for dynamic analysis of bridges under moving vehicles considering separation and surface roughness[J]. Computers & Structures, 2015, 154: 135-144.

[12] ST?NCIOIU D, OUYANG Huajiang, MOTTERSHEAD J E. Dynamics of a beam and a moving two-axle system with separation[J]. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, 2008, 222(10): 1947-1956.

[13] STANCIOIU D, OUYANG Huajiang, MOTTERSHEAD J E. Vibration of a beam excited by a moving oscillator considering separation and reattachment[J]. Journal of Sound and Vibration, 2008, 310(4/5): 1128-1140.

[14] CHENG Y S, AU F T K, CHEUNG Y K, et al. On the separation between moving vehicles and bridge[J]. Journal of Sound and Vibration, 1999, 222(5): 781-801.

[15] BAEZA L, OUYANG Huajiang. Dynamics of a truss structure and its moving-oscillator exciter with separation and impact-reattachment[J].Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2008, 464(2098): 2517-2533.

[16] 翟婉明. 车辆-轨道耦合动力学[M]. 北京: 科学出版社, 2015: 434-442.

ZHAI Wanming. Vehicle-track coupled dynamics[M]. Beijing: Science Press, 2015: 434-442.

[17] 娄平. 列车-轨道(桥梁)系统竖向振动分析[D]. 长沙: 中南大学土木工程学院, 2007: 27-50.

LOU Ping. Vertical vibration analysis of train-track/bridge system[D]. Changsha: Central South University. School of Civil Engineering and Architecture, 2007: 27-50.

[18] 曾京, 胡松. 轮轨摩擦碰撞及脱轨的研究[J]. 振动工程学报, 2001, 14(1): 1-6.

ZENG Jing, HU Song. Study on frictional impact and derailment for wheel and rail[J]. Journal of Vibration Shock, 2001, 14(1): 1-6.

[19] 关庆华, 曾京. 轮轨横向碰撞引起的脱轨研究[J]. 振动与冲击, 2009, 28(12): 38-42.

GUAN Qinghua, ZENG Jing. Study on derailment induced by lateral impact between wheel and rail[J]. Journal of Vibration and Shock, 2009, 28(12): 38-42.

[20] 肖新标. 复杂环境状态下高速列车脱轨机理研究[D]. 成都: 西南交通大学牵引动力国家工程实验室, 2013: 21-88.

XIAO Xinbiao. Study on high-speed train derailment mechanism in severe environment[D]. Chengdu: Southwest Jiaotong University. State Key Laboratory of Traction Power, 2013: 21-88.

[21] 金学松, 郭俊, 肖新标, 等. 高速列车安全运行研究的关键科学问题[J]. 工程力学, 2009, 26(A2): 8-22, 105.

JIN Xuesong, GUO Jun, XIAO Xinbiao, et al. Key scientific problems in the study on running safety of high speed trains[J]. Engineering Mechanics, 2009, 26(A2): 8-22, 105.

[22] 孙丽霞, 姚建伟. 高速铁道车辆蛇形脱轨安全性评判方法研究[J]. 中国铁道科学, 2013, 34(5): 82-92.

SUN Lixia, YAO Jianwei. Hunting derailment safety evaluation method of high speed railway vehicle[J]. China Railway Science, 2013, 34(5): 82-92.

[23] XU Lei, ZHAI Wanming. A three-dimensional dynamic model for train-track interactions[J]. Applied Mathematical Modelling, 2019, 76: 443-465.

[24] ZHANG Yao, ZHAO Haisheng, LIE S T. A nonlinear multi-spring tire model for dynamic analysis of vehicle-bridge interaction system considering separation and road roughness[J]. Journal of Sound and Vibration, 2018, 436: 112-137.

[25] JU S H. A frictional contact finite element for wheel/rail dynamic simulations[J]. Nonlinear Dynamics, 2016, 85(1): 365-374.

(编辑  陈灿华)

收稿日期： 2020 -06 -12; 修回日期： 2020 -08 -22

基金项目(Foundation item)：国家自然科学基金资助项目(52068028)；国家自然科学基金委员会与神华集团有限公司联合资助项目(U1261113)；江西省教育厅科学技术项目(GJJ170392)；华东交通大学科研启动基金资助项目(10012003416034) (Project(52068028) supported by the National Natural Science Foundation of China; Project(U1261113) jointly supported by the National Natural Science Foundation Committee of China and Shenhua Group Corporation Limited; Project(GJJ170392) supported by Jiangxi Provincial Department of Education; Project(10012003416034) supported by Research Startup Fund of East China Jiaotong University)

通信作者：向俊，博士，教授，从事列车脱轨控制、列车-轨道(桥梁)系统空间振动及铁路轨道结构等研究；E-mail：jxiang@mail.csu.edu.cn

引用格式： 向俊, 陈林, 苏玮, 等. 铁道车辆车轮跳轨全过程计算方法[J]. 中南大学学报(自然科学版), 2021, 52(2): 614-624.

Citation: XIANG Jun, CHEN Lin, SU Wei, et al. A method of calculating the whole process of wheel-jumping of railway vehicle[J]. Journal of Central South University(Science and Technology), 2021, 52(2): 614-624.