Simulation Calculation Method and Application for Relative Displacement of Heavy Hall Freight Suspension on Curved Track
-
摘要:
为准确求解曲线轨道上重载货车悬挂的相对位移,首先,建立曲线轨道数学模型,推导出曲线外轨超高、顺坡角、侧滚角和中心角随线路长度的变化公式,再根据车辆各刚体部件进出曲线的时间和所处曲线位置差异,编程计算悬挂点刚体间的超高及转角差;其次,以刚体质心为坐标原点建立本体坐标系,分别给出悬挂点在两刚体本体坐标系中的坐标表达式,通过坐标变换法将本体坐标转换到同一坐标系下,计算悬挂点瞬态相对位移;最后,结合车辆曲线动力学仿真程序计算,即可求出车辆曲线通过时各悬挂点的动态相对位移. 计算结果表明:车辆悬挂相对位移是车辆参数和曲线轨道参数综合作用的结果,当单独不计线路侧滚角差、顺坡角差、中心角差时,对应悬挂相对位移的最大偏差率可达42.85%、24.03%、71.42%;利用坐标变换结合动力学仿真计算的方法可全面考虑车辆和轨道参数,求解车辆悬挂相对位移更为准确.
Abstract:In order to accurately solve the relative displacements of heavy haul freight suspension on a curved track, firstly the mathematical model of curved track was established, and the formulas of outer rail superelevation, slope angle, roll angle and central angle with the length of curve were derived, then according to the time difference and the geometric position difference between the suspension components while the vehicle entering or leaving the curved track, programming to calculate the relative rail superelevation and angle difference of these rigid components. Secondly, the ontological coordinate system was established with the centroid of each rigid body as the origin, and the coordinate expressions of suspension points in the two ontology coordinate systems were given respectively, and converting these coordinates into a same coordinate system through coordinate transformation method, so the transient relative displacement of suspension points could be calculated. Finally, by combining the vehicle curving dynamics simulation, the dynamic relative displacements of suspension points could be calculated. The results indicate that the relative displacement of vehicle suspension is the result of the comprehensive action of vehicle parameters and curve track parameters, when the difference of roll angle difference, slope angle and central angle are excluded separately, the corresponding maximum deviation rate of relative suspension displacement can arrive at 42.85%, 24.03%, and 71.42%. The method of coordinate transformation combined with dynamic simulation can fully consider the vehicle and track parameters, and the relative displacement of vehicle suspension can be solved more accurately.
-
机车车辆的曲线通过问题是车辆动力学系统研究课题中的一个重要领域,深入探究和理解车辆曲线通过行为及其机理是合理设计车辆及其轨道参数,提高车辆曲线通过性能的理论基础. 对车辆曲线通过理论及计算模型的研究经历了早期的经典曲线通过理论、线性稳态曲线通过、非线性稳态曲线通过和现阶段的非线性动态曲线通过理论4个阶段[1-2],反映了对车辆曲线通过机理由定性到定量、由个体到整体、从局部到全面的系统认识过程,研究更趋实际和完善. 但动态曲线通过理论分析也更为复杂,不仅要考虑线路参数、轮轨接触关系变化,还要考虑车辆悬挂及其动态响应,即准确求解出车辆悬挂变化成为了车辆动态曲线通过计算的前提. 但查阅一般车辆悬挂计算的相关文献[1-14]发现,对车辆系统悬挂装置相对位移及悬挂力的计算并不多见,也不统一. 在常见车辆动态曲线通过计算中,大多仅列出各刚体的运动方程,而对车辆悬挂相对位移及悬挂力的具体计算却未给出详细的推导[1-3];翟婉明[4]在《车辆-轨道耦合动力学》中对我国铁路机、客、货等基本车型的悬挂力进行了推导,并给出了相应计算公式,但书中货车模型采用的转向架为传统的轴箱导框式结构,与现有重载货车一系悬挂结构存在较大差异,在悬挂力计算公式中,仅考虑了线路曲率和超高影响,而忽略了悬挂点刚体间的侧滚角差和顺坡角差;王开云等[5-6]通过坐标变换方法分析了车辆在曲线轨道上的悬挂相对位移及其对轮重特性的影响,并给出了车辆在不同曲线段(缓和曲线、圆曲线)上悬挂点相对平移和转角的变化表,但为简化计算,忽略了悬挂点不同刚体间的顺坡角差;多数使用的SIMPACK、UM等动力学仿真软件往往是直接利用软件自带的曲线模块进行参数化建模及计算,对于软件如何处理悬挂相对位移及如何计算悬挂力的内核却难以知晓[7-14].
在曲线轨道,因外轨超高及曲率变化,内外钢轨形成的轨道面实际为曲面,在该轨道面上运行的车辆将产生动态扭曲,使得车辆悬挂点各刚体间的参考坐标系不再平行,须将不同坐标系变换到同一坐标系下才可进行相应计算[5]. 而重载货车及其转向架在沿用传统的两级刚度弹簧 + 摩擦斜楔减振的基础上,又增加了一系轴箱橡胶垫和常接触弹性滚子旁承等弹性支撑悬挂,悬挂点更多,悬挂计算也更为复杂[15-18]. 但无论车辆系统各刚体如何变化,悬挂点不同刚体坐标系(即质心坐标系)间的关系都可认为是某一坐标系(OB-XBYBZB)是另一坐标系(OA-XAYAZA)经平动(a,b,c)和转动(α,θ,γ)后形成的新坐标系,其中:α为侧滚角;θ为顺坡角;γ为中心角. 则两坐标系的相互变化关系为
[XBYBZB]=[1γ−θ−γ1αθ−α1][XA−aYA−bZA−c]=ArA−Ar0, (1) [XAYAZA]=A−1[XBYBZB]+[abc]=A−1rB+r0, (2) 式中:
${\boldsymbol{A}}= \left[1γ−θ−γ1αθ−α1\right] ; {\boldsymbol{r}}_{\mathrm{A}}=\left[XAYAZA\right] $ ;rB=[XBYBZB];r0=[abc]. 所以,只要知道车辆悬挂点刚体间的相对平动(a,b,c)和转动(α,θ,γ)大小,通过式(1)或式(2)将不同刚体的坐标转换到同一参考坐标系下,即可进行车辆悬挂点相对位移及速度的计算.
1. 曲线轨道参数模型
铁路曲线一般包括水平曲线和竖曲线,本文只涉及水平曲线分析. 典型的曲线轨道如图1所示,由圆曲线连接前、后缓和曲线组成. 从直线经前缓和曲线到圆曲线时,线路曲率k和超高h由0分别增至1/R0(R0为圆曲线半径)和h0;圆曲线处,超高和曲率保持不变;从圆曲线经后缓和曲线到直线时,线路曲率k和超高h则分别由1/R0和h0减至0. 图1中坐标系定义如下(坐标轴符合右手法则):
1) O-XYZ绝对坐标系,固结于水平面轨道中心处. X轴沿行车方向水平向前,Y轴方向水平向内,Z轴方向垂直向下;
2) O1-X1Y1Z1质点运动坐标系,固结在轨道曲面上运动的刚体质心,相对绝对坐标系X、Y和Z轴的转角分别为(α,θ,γ).
为缓和因超高和曲率突变引起的冲击,在直线和圆曲线间需设置一条曲率和超高均连续变化的缓和曲线. 我国一般采用直线超高顺坡三次抛物线型缓和曲线[8],即外轨超高h和曲率k都随缓和曲线长度线性变化. 根据图1所示和三次抛物线型的特性,则某一质点沿曲线轨道中心线运行时的超高h和转角(α,θ,γ)与各曲线段走行距离l的关系为
h={l2lh1h0,前缓和曲线,l=0∼lh1,h02,圆曲线,l=0∼ly,lh2−l2lh2h0,后缓和曲线,l=0∼lh2. (3) α={l2lh1a0h0,前缓和曲线,l=0∼lh1,h02a0,圆曲线,l=0∼ly,lh2−l2lh2a0h0,后缓和曲线,l=0∼lh2, (4) 式中:
$ {a}_{0} $ 为内外钢轨中心间距的一半.θ={h02lh1,前缓和曲线,l=0∼lh1,0,圆曲线,l=0∼ly,−h02lh2,后缓和曲线,l=0∼lh2. (5) γ={l22R0lh1,前缓和曲线,l=0∼lh1,lh12R0+lR0,圆曲线,l=0∼ly,lh12R0+lyR0+l22R0lh2,后缓和曲线,l=0∼lh2. (6) 2. 悬挂点超高差和转角差的求解
重载货车基本采用一系和二系两系悬挂,主要悬挂点有轮对与侧架间的轴箱悬挂、侧架与摇枕间的中央悬挂以及车体与摇枕间的弹性旁承支撑悬挂及心盘支撑等. 图2是重载货车进入曲线轨道后,各刚体部件的位置关系示意,图中:lt为转向架固定轴距的一半;lc为车辆定距的一半;ds为二系悬挂横向跨距的一半;dw为一系悬挂横向跨距的一半;R为线路曲线半径;X、Y、Z分别为纵向、横向和垂向位移自由度,ψ为摇头自由度,下标w、B分别表示轮对、摇枕,下标tL、tR分别表示左、右侧转向架. 假定轮对、摇枕和车体中心均沿曲线轨道中心向前运动,并以车辆1位轮对作为起始参考点,根据式(3) ~ (6)和车辆各刚体前、后和左、右位置距离关系,可确定出车辆各悬挂点刚体部件对应的超高和转角.
2.1 轮 对
假定车辆沿曲线轨道纯滚动时,可认为轮对中心与轨道中心重合,则轮对的超高和相应转角可直接利用式(3) ~ (6)进行计算,但需将式中l换成lwi,lwi为1 ~ 4位轮对在不同曲线段的走行距离,其中i = 1, 2, 3, 4;若以1位轮对lw1作为起始点,则lw2 = lw1-2lt,lw3 = lw1-2lc,lw4 = lw1-2(lc + lt).
2.2 摇 枕
根据图2,可认为摇枕质心与曲线轨道中心线重合. 则相应的超高和转角也可直接按轨道曲线中心公式计算,只需将式中l换成lbj即可,lbj为前后摇枕在不同曲线段的走行距离,j = 1,2,分别表示前、后摇枕;若以1位轮对lw1作为起始点,则lb1 = lw1−lt,lb2 = lw1−lt−2lc.
2.3 侧 架
由于同一转向架的摇枕和左、右侧架的质心在横向共线,若忽略两者对应的曲率差异,则侧架的侧滚角和中心角与摇枕一致. 而左、右侧架与摇枕因横向位置差异,其超高和顺坡角与摇枕不同,但可根据几何位置关系推导,如式(7)、(8).
超高:
htLj=hbj+dsαbj. (7) 顺坡角:
θtLj={h02lh1+dsh02a0lh1,前缓和曲线,0,圆曲线,−(h02lh2+dsh02a0lh2),后缓和曲线, (8) 式中:
$ {h}_{\mathrm{b}j} $ 和$ {\alpha }_{\mathrm{b}j} $ (j = 1,2)分别为前、后摇枕的超高和中心角;下标L表示左侧架,当其改为R表示右侧架时,则式(7)、(8)右边的 + 号改为$ -\mathrm{号} $ .2.4 车 体
车体质心与曲线轨道中心线重合,则车体超高、曲率和各转角可直接利用式(3) ~ (6)计算,但需将式中l换成车体在各曲线段的走行距离lC,仍以1位轮对作为起始点,则lC = lw1−lt−lc.
计算出车辆各刚体的具体超高和转角后,即可计算出各悬挂点刚体间的相对超高差及转角差.
3. 悬挂点相对位移求解
计算车辆系统各悬挂点的相对位移必须要在同一参考坐标系下进行. 先需对车辆系统各刚体定义本体坐标系,如图1中质点坐标系O1-X1Y1Z1所示,坐标原点固结于各刚体质心;横向定义为刚体所在曲线位置的径向方向,即指向曲率中心;纵向为各刚体所在曲线位置的切向方向;而垂向定义为竖直向下方向.
求解悬挂点相对位移步骤如下:
步骤1 先建立各刚体的本体坐标系,并根据各悬挂点的几何位置和刚体运动自由度分析确定悬挂点在各刚体本体坐标系的坐标表达式;
步骤2 计算或直接给出悬挂点两刚体在曲线某位置的相对平动(a,b,c)和转动(α,θ,γ)值;
步骤3 按式(1)或式(2)进行坐标变换,将悬挂点的本体坐标转换到同一坐标系下,两者相减得出各悬挂点相对位移的表达式;
步骤4 结合车辆动力学计算求出的各刚体运动自由度值,代入相对位移表达式即可求出车辆各悬挂点的相对位移. 本文以求解车辆1位轮对和前左侧架悬挂点的相对位移作为算例来介绍整个过程,其他悬挂点的计算可参照执行.
令1位轮对的纵向、横向、垂向位移分别为Xw1、Yw1、Zw1,侧滚、点头、摇头分别为ϕw1、βw1、Ψw1;前左侧架的纵向、横向、垂向位移分别为XtL1、YtL1、ZtL1,侧滚、点头和摇头分别为ϕtL1、βtL1、ΨtL1;轮对质心与侧架质心垂向距离为htw. 则悬挂点在1位轮对和前左侧架本体坐标系的坐标分别如式(9)~式(11).
悬挂点在前左侧架本体坐标系的位置:
rtt=[XtYtZt]=[XtL1+htwβtL1+ltYtL1+ltψtL1ZtL1+htw−ltβtL1]. (9) 悬挂点在1位轮对本体坐标系的位置:
rww=[XwYwZw]=[Xw1+dwψw1Yw1−dwZw1−dwϕw1]. (10) 1位轮对质心相对前左侧架质心转动和平移为
{A=[1γ−θ−γ1αθ−α1],r0=[abc]≈[ltdw+ltγ2htw−Δhw1tL1], (11) 式中:
$\Delta {h}_{\mathrm{w}1\mathrm{t}\mathrm{L}1}$ 为1位轮对和前左侧架的超高差.将1位轮对悬挂点本体坐标向左侧架本体坐标系进行转换并简化(忽略小变量间的乘积项)得到
rtw=A−1rww+r0=[Xw1+dw(ψw1+γ)+ltYw1+ltγ2Zw1−dw(ϕw1+α)+htw−Δhw1tL1]. (12) 悬挂点相对位移为
Δrtw=rtt−rtw=[XtL1+htwβtL1+ltYtL1+ltψtL1ZtL1+htw−ltβtL1]−[Xw1+dw(ψw1+γ)+ltYw1+ltγ2Zw1−dw(ϕw1+α)+htw−Δhw1tL1]=[XtL1+htwβtL1−dw(ψw1+γ)−Xw1YtL1+ltψtl1−Yw1−ltγ2ZtL1−ltβtL1−Zw1+dw(ϕw1+α)+Δhw1tL1]. (13) 4. 仿真算例及应用
以我国27 t轴重C80E货车为例[17-18],基于交叉支撑转向架重载货车-轨道耦合模型[19],采用该方法对车辆通过曲线时各悬挂点的相对转角和相对位移进行动力学仿真. 车辆主要结构参数为:lc=4605 mm,lt=930 mm,dw和ds均为1003 mm,a0=746.5 mm. 曲线参数为:曲线半径800 m,最大外轨超高95 mm,前、后缓和曲线长度75 m,圆曲线长度50 m,前、后直线长度25 m. 车辆运行速度为80 km/h,即车辆以平衡最大超高通过曲线,且不考虑线路轨道不平顺谱的激扰.
图3为车辆通过该曲线时,车辆系统各悬挂点刚体间超高差的变化,由图3可见:车辆各刚体部件都在直线上时,超高差为0;在缓和曲线上时,因外轨超高随缓和曲线长度线性变化,而轮对质心与左、右侧架质心横向间距为dw,纵向间距为lt,故轮对与侧架以及侧架与摇枕悬挂点的超高差值也呈线性变化,但在前、后缓和曲线上的变化方向相反,且轮对与左、右侧架以及摇枕与左、右侧架超高差也刚好相反;而摇枕质心和车体质心的横向距离为0,纵向间距为lc,超高差由两者进出曲线前后不同步引起,所以其超高差值在缓和曲线上保持不变,其值也明显比轮对与侧架以及侧架与摇枕的超高差小;在圆曲线上,因外轨超高保持最大超高不变,悬挂点的超高差也相应不变;在直缓点和缓圆点附近,因悬挂点刚体质心同时处在不同的曲线段,超高差则会发生明显的转折.
悬挂点侧滚角差值的变化如图4所示,因同一转向架的左、右侧架质心和摇枕质心在径向共线,即它们在同一时刻的侧滚角相同,故车辆悬挂点的侧滚角差只需考虑各轮对与前、后侧架以及前、后摇枕与车体间的侧滚角差. 由图4可见:在直线和圆曲线段,悬挂点的侧滚角差值均为0;在缓和曲线段,各侧滚角差值最大并保持不变,其中前、后轮对与侧架以及前、后摇枕与车体间的侧滚角差方向相反,在前、后缓和曲线上的方向也刚好相反;在直缓点和缓圆点附近,悬挂点的侧滚角差会出现明显的突变,特别是轮对与侧架间的悬挂点,几乎呈直角变化,摇枕与车体间的侧滚角差变化也很快,且差值较大,最大为3.907 mrad.
图5为各悬挂点顺坡角差值的变化. 由图5可知:当车辆完全处于直线或圆曲线段时,各悬挂点顺坡角差值为0;在缓和曲线段时,轮对与侧架间、侧架与车体间的顺坡角差值最大并保持不变,且在前、后缓和曲线上的方向相反,而摇枕与车体因其质心在纵向共线,在缓和曲线段两者顺坡角相同,故顺坡角差值为0;在直缓点和缓圆点附近,因悬挂点两刚体质心处在不同的曲线段上,一个为0,一个为定值,所以会出现明显的跳跃现象,且左侧比右侧大(因左侧外轨超高大,故顺坡角更大).
悬挂点中心角差(或称摇头角差)变化如图6所示. 同理,因左、右侧架和摇枕质心在径向共线,忽略内、外侧轨道的曲率差异,同一转向架的侧架和摇枕的中心角可认为是同步变化,则只需考虑轮对与侧架以及摇枕与车体间的悬挂点中心角差. 由图6可见:在直线段,中心角差均为0;在前缓和曲线段,中心角差随轨道长度线性增加,后缓和曲线段,则线性减小,且前、后轮对与侧架和前、后摇枕与车体间的中心角差的正负值刚好相反,前为正,后为负;在圆曲线,各悬挂点的中心角差达到最大且保持不变.
求出各悬挂点刚体动态的相对平动和转动后,给定重载货车一二系悬挂刚度和阻尼参数,代入车辆系统运动微分方程,进行车辆曲线通过动力学仿真计算,即可求解出重载货车系统各悬挂点的相对位移. 图7是重载货车一系和二系悬挂点相对位移的计算结果(图中均是前项减后项). 由图7可见:因曲线轨道的曲率和超高变化,重载货车通过曲线时,其一系和二系悬挂相对位移均不断变化,车辆的前、后悬挂和左、右悬挂均存在不同程度差异,且与曲线线路的超高、侧滚角、顺坡角和中心角差值的单个变化特性也存在较大差异. 仿真发现:当单独不计线路侧滚角差、顺坡角差或中心角差时,与考虑全部线路偏差角时悬挂相对位移极值的最大偏差率分别可达42.85%、24.03%和71.42%. 这说明重载货车通过曲线时的悬挂相对位移是线路参数和车辆参数综合作用的结果,仿真计算中要准确确定车辆悬挂的相对位移,需要全面考虑车辆结构参数和轨道几何参数,忽略其中某些参数,将影响车辆悬挂相对位移及悬挂力,进而影响车辆系统动力学仿真计算结果的准确性.
5. 结 论
1) 悬挂点相对位移与线路侧滚角、顺坡角和摇头角差的变化特性均存在差异,说明车辆通过曲线时的悬挂相对位移是车辆参数或曲线轨道参数综合作用的结果. 但从各悬挂点相对转角和相对位移随曲线线路变化看,悬挂点相对位移在缓和曲线上的变化更为明显,且在前、后缓和曲线上的变化趋势相反,而在圆曲线上相对位移变化不大,甚至无变化.
2) 通过先确定车辆在曲线轨道上悬挂点各刚体间的平动和转角差,再通过坐标变换计算车辆悬挂点相对位移的方法,可以准确计算出车辆各悬挂点的相对位移. 利用该方法进行计算机编程,可推广至空间任意两刚体间相对位移及速度的仿真计算.
3) 该方法考虑了车辆和轨道结构相关参数,无需简化处理,使车辆在曲线轨道上悬挂力的计算更为准确,一定程度上丰富和完善了车辆曲线通过仿真计算方法. 基于此方法开展车辆及轨道相关参数的优化分析,可为铁路车辆及线路工程设计、维修等实际应用提供一定指导.
致谢:牵引动力国家重点实验室开放课题资助项目(TPL2209)、高校博士科研启动基金资助项目(MY2015B009).
-
-
[1] 王福天. 车辆系统动力学[M]. 北京: 中国铁道出版社, 1994: 144-149. [2] 姚建伟. 机车车辆动力学[M]. 北京: 科学出版社, 2014: 251-257. [3] 刘学毅, 王平. 车辆-轨道-路基系统动力学[M]. 成都: 西南交通大学出版社, 2010: 51-55. [4] 翟婉明. 车辆-轨道耦合动力学(上册) [M]. 4版. 北京: 科学出版社, 2015: 57-64. [5] 王开云,翟婉明. 直线与曲线轨道上车辆悬挂相对位移的计算[J]. 西南交通大学学报,2003,38(2): 122-126. doi: 10.3969/j.issn.0258-2724.2003.02.002WANG Kaiyun, ZHAI Wanming. Calculation of displacements of vehicle suspension on tangent and curved tracks[J]. Journal of Southwest Jiaotong University, 2003, 38(2): 122-126. doi: 10.3969/j.issn.0258-2724.2003.02.002 [6] 刘鹏飞,翟婉明,王开云,等. 机车车辆通过缓和曲线时悬挂系统及轮重的动态特性[J]. 中国铁道科学,2013,34(1): 67-74. doi: 10.3969/j.issn.1001-4632.2013.01.10LIU Pengfei, ZHAI Wanming, WANG Kaiyun, et al. The dynamic characteristics of suspension system and wheel load for rolling stock passing through transition curve[J]. China Railway Science, 2013, 34(1): 67-74. doi: 10.3969/j.issn.1001-4632.2013.01.10 [7] 李亨利,李芾,潘树平,等. 新型货车径向转向架低动力性能与连接可靠性[J]. 交通运输工程学报,2013,13(5): 39-46. doi: 10.3969/j.issn.1671-1637.2013.05.006LI Hengli, LI Fu, PAN Shuping, et al. Low track force dynamic performance and connection reliability of new type radial bogie for freight car[J]. Journal of Traffic and Transportation Engineering, 2013, 13(5): 39-46. doi: 10.3969/j.issn.1671-1637.2013.05.006 [8] 张建全,黄运华,李芾. 缓和曲线线型对铁道车辆动力学性能的影响[J]. 交通运输工程学报,2010,10(4): 39-44. doi: 10.3969/j.issn.1671-1637.2010.04.007ZHANG Jianquan, HUANG Yunhua, LI Fu. Influence of transition curves on dynamics performance of railway vehicle[J]. Journal of Traffic and Transportation Engineering, 2010, 10(4): 39-44. doi: 10.3969/j.issn.1671-1637.2010.04.007 [9] 高亮,王璞,蔡小培,等. 基于多车精细建模的曲线地段重载列车-轨道系统动力性能研究[J]. 振动与冲击,2014,33(22): 1-6,12.GAO Liang, WANG Pu, CAI Xiaopei, et al. Dynamic characteristics of train-track system in curved track sections based on elaborate multi-vehicle model[J]. Journal of Vibration and Shock, 2014, 33(22): 1-6,12. [10] TSUJIE M, MIURA M, CHEN H, et al. A study on the initiation of head check of the low rail using multibody dynamics[J]. Wear, 2019, 436/437: 202989.1-202989.10. doi: 10.1016/j.wear.2019.202989 [11] WU L, YAO X S, VANDERMAREL J, et al. Effects of curve radius and rail profile on energy saving in heavy haul achieved by application of top of rail friction modifier[J]. Wear, 2016, 366/367: 279-286. doi: 10.1016/j.wear.2016.07.003 [12] ZHAI W M, GAO J M, LIU P F, et al. Reducing rail side wear on heavy-haul railway curves based on wheel-rail dynamic interaction[J]. Vehicle System Dynamics, 2014, 52(S1): 440-454. [13] ZBOIŃSKI K, WOŹNICA P. Optimisation of the railway transition curves' shape with use of vehicle-track dynamical model[J]. Archives of Transport, 2010, 22(3): 387-406. [14] ZBOINSKI K, GOLOFIT-STAWINSKA M. Investigation into nonlinear phenomena for various railway vehicles in transition curves at velocities close to critical one[J]. Nonlinear Dynamics, 2019, 98(3): 1555-1601. doi: 10.1007/s11071-019-05041-2 [15] 邵文东,李立东,周国东,等. 我国铁路重载货车转向架技术及发展[J]. 铁道车辆,2012,50(12): 15-18. doi: 10.3969/j.issn.1002-7602.2012.12.005SHAO Wendong, LI Lidong, ZHOU Guodong, et al. Technologies and development of railway heavy haul freight car bogies in our country[J]. Rolling Stock, 2012, 50(12): 15-18. doi: 10.3969/j.issn.1002-7602.2012.12.005 [16] 刘宏友,邵文东,王宝磊,等. 货车转向架60年的回顾与展望[J]. 铁道车辆,2013,51(12): 42-48. doi: 10.3969/j.issn.1002-7602.2013.12.006LIU Hongyou, SHAO Wendong, WANG Baolei, et al. Review and prospects of the freight car bogies in the recent 60 years[J]. Rolling Stock, 2013, 51(12): 42-48. doi: 10.3969/j.issn.1002-7602.2013.12.006 [17] 段仕会,徐世锋,李立东,等. DZ1型转向架的研制[J]. 铁道车辆,2016,54(6): 31-34. doi: 10.3969/j.issn.1002-7602.2016.06.010DUAN Shihui, XU Shifeng, LI Lidong et al. Development of DZ1 bogies[J]. Rolling Stock, 2016, 54(6): 31-34. doi: 10.3969/j.issn.1002-7602.2016.06.010 [18] 张志彬,刘爽,李华,等. C80E(H、F)型通用敞车研制[J]. 铁道车辆,2018,56(6): 22-24. doi: 10.3969/j.issn.1002-7602.2018.06.009ZHANG Zhibin, LIU Shuang, LI Hua, et al. Development of C80E(H, F) general purpose gondola cars[J]. Rolling Stock, 2018, 56(6): 22-24. doi: 10.3969/j.issn.1002-7602.2018.06.009 [19] 杨春雷. 重载货车轴重与速度匹配关系研究[D]. 成都: 西南交通大学, 2013: 30-45. -