Influence of Track Parameters on Wheel/Rail Contact Behavior of High-Speed Turnout
-
摘要:
为研究60N钢轨350 km/h 18号高速道岔合理的轨距和轨底坡,利用60N钢轨高速道岔关键断面和实测LMA磨耗车轮,基于迹线法原理和Kalker三维非赫兹滚动接触理论,分析不同轨距和轨底坡参数下的轮轨接触几何和力学特性,并与CHN60钢轨高速道岔计算结果进行对比. 结果表明:在保证安全的前提下适当将轨距加宽可改善轮轨匹配关系,提升列车过岔平稳性,减小轮对横移量大于8 mm时的轮轨接触应力和表面滚动接触疲劳因子,延长尖轨使用寿命;轨底坡为1/30、1/40和1/50时,轮轨接触参数相差较小,匹配性能较优;轨底坡为1/10和1/20时,横向不平顺和轮轨滚动接触疲劳因子普遍较大,且1/10轨底坡对车轮磨耗的适应性较差;与CHN60钢轨高速道岔相比,60N钢轨高速道岔的等效锥度普遍更小,列车过岔平稳性更优;车轮磨耗易导致车轮在轮轨过渡区段空转,引起尖轨伤损.
Abstract:In order to investigate reasonable rail gauge and rail cant of No.18 turnout with 60N rail under the straight passing speed 350 km/h, based on the track line method and Kalker’s 3D non-Herz rolling contact theory, key sections of 60N rail turnout and measured LMA worn wheels were used to analyzed the wheel-rail contact geometric and mechanical characteristics under different rail gauges and different rail cant parameters. The calculation results were compared with those of CHN60 rail turnout. The results show that under the premise of ensuring safety, properly widening the rail gauge can improve the wheel-rail matching relationship, improve the stability of the train went through the turnout, reduce the wheel-rail contact stress and surface rolling contact fatigue factor when the wheelset-lateral displacement is greater than 8 mm, and extend the service life of the rail. When the rail cant is 1/30, 1/40 and 1/50, the wheel-rail contact parameters have little difference and the matching performance is better. When the rail cant is 1/10 and 1/20, the lateral irregularity and wheel-rail rolling contact fatigue factors are generally large, and the adaptability of 1/10 rail cant to wheel wear is poor. Compared with CHN60 rail turnout, the equivalent conicity of 60N rail turnout is generally smaller, and the stability of the train went through the turnout is better. Wheel wear is easy to lead to wheelspin in the wheel-rail transition section, resulting in switch rail damage.
-
Key words:
- rail gauge /
- rail cant /
- 60N rail /
- worn wheel /
- high-speed turnout /
- wheel-rail contact
-
道岔是机车车辆从一股轨道转入或越过另一股轨道的线路设备,是铁路轨道的重要组成部分,也是线路上的薄弱环节[1]. 轨距和轨底坡是高速道岔重要的设计技术条件,其参数的设置直接影响行车安全和行车品质. 60N钢轨在区间线路上的使用取得了成功,为提升列车过岔的平稳性和轮轨接触力学特性,在高速道岔区同样需要应用60N钢轨. 列车运行过程中,车轮始终处于动态磨损状态,对轮轨接触行为影响较大.
致谢:京沪高速铁路股份有限公司科技研究项目(京沪科研-2020-11).
针对线路上的轨距和轨底坡参数,国内外学者做了大量的研究. 杜星等[2]建立了LMD车轮和CHN60钢轨匹配的动力学模型,分析同一轨道在不同轨底坡条件下的动力学行为发现,轨底坡变化对列车直线运行时的平稳性、舒适性影响很大. 钱瑶等[3]对比分析了不同轨底坡下60N钢轨和高速车轮LMA、XP55、S1002G匹配时的轮轨接触行为,匹配较优的轨底坡是1/20和1/30. 陈嵘等[4-5]研究了我国地铁线路常用的LM型面与CHN60钢轨在不同轨距和非对称轨底坡下的轮轨接触特性,发现增大轨距和改变轨底坡可改善轮轨匹配关系. Cui等[6]提出一种优化轮轨廓形的“正向求解法”,车辆动力学行为结果表明,轨底坡1/30、轨距1435 mm是优化后廓形的最优轨道参数. Sánchez等[7]开发了一个严谨的测量轨距和轨底坡的程序确保测量过程的准确性. 李超等[8]分析了转辙器区采用动态轨距加宽技术的轮轨接触作用,该技术可减轻尖轨磨耗和滚动接触疲劳. Ye等[9]通过建立多体动力学模型,分析轨道参数对车轮磨耗的影响,发现轨道参数对车轮磨耗和脱轨安全性有较大影响. 闫正等[10]分析了高速动车组车轮踏面和高速60N钢轨道岔断面的静态接触特性,发现适当地增大轨距和轨底坡有利于改善轮轨接触状态. 上述研究表明,改变轨距和轨底坡参数对轮轨接触行为影响较大,且现行的1/40轨底坡和1435 mm轨距往往不是最优的轨道参数,而国内的研究主要集中在区间线路或是针对CHN60钢轨道岔,因此,有必要对新型350 km/h 60N钢轨18号高速道岔的合理轨距和轨底坡展开研究.
本文根据实测LMA磨耗车轮型面和60N钢轨高速道岔关键断面,基于迹线法原理和三维非赫兹滚动接触理论,建立道岔区轮轨滚动接触模型,计算不同轨距和轨底坡下的接触几何参数和静力学指标,并与CHN60钢轨高速道岔进行对比.
1. 道岔关键断面和磨耗车轮
1.1 60N钢轨关键断面
利用CAD导出各关键断面离散数据,将尖轨和基本轨分开,并线性插值,得到如图1所示60N钢轨18号高速道岔转辙器区钢轨模型,X为沿钢轨的纵向里程坐标,Y为钢轨横坐标,Z为钢轨竖向坐标,尖轨顶宽35 mm钢轨断面位于轮载过渡段,其接触行为较为复杂,轮轨相互作用剧烈,对研究道岔区的轮轨接触行为具有代表性,因此,本文选取尖轨顶宽35 mm关键断面进行计算. 60N与CHN60钢轨高速道岔在尖轨顶宽35 mm处钢轨廓形见图2.
1.2 LMA磨耗车轮
跟踪记录某线路上运行的CRH2型动车组上LMA车轮型面演变[11],从标准LMA磨耗车轮到运营里程15万km,再到25万km,磨耗量增大较为明显,本文取这三种磨耗车轮进行分析. 车轮凹型磨耗对轮对接触具有较大影响[12],随着列车运营里程的增加,凹型磨耗车轮的磨耗量大致呈线性增大,如图3所示. 对磨耗车轮廓形通过三次样条函数进行插值、拟合.
2. 计算原理
2.1 轮轨接触几何计算原理
利用二维迹线法原理,不考虑车轮摇头角,由最小距离法搜索轮轨接触点. 接触点位置的求解有两个等价几何条件:1) 轮轨接触点处轮轨垂直距离为0,非接触点轮轨垂直距离大于0;2) 轮轨接触点处轮轨的轮廓线具有相同的斜率. 本文采用条件1求解轮轨接触点,再用条件2对结果加以验证.
2.2 UIC 519标准计算等效锥度
EN 15302标准[13]中,等效锥度的计算采用轮对周期运动的假设,这个方法称之为UIC 519标准[14]等效锥度. 自由轮对在轨道上的运动用微分方程表示为
¨y+v2er0Δr=0, (1) 式中:
y 为轮对在钢轨上的横向位移;v 为车辆前进速度;e 为接触点跨距;r0 为名义滚动圆半径;Δr 为右侧和左侧车轮滚动半径差.最后应用Klingel公式计算等效锥度
γe :γe=2er0(π λ)2 (2) 式中:
λ 为轮对在钢轨上做正弦运动的波长,其求解见文献[14].2.3 轮轨滚动接触力学原理
本文在接触力学部分采用Kalker的三维非赫兹滚动接触理论及其数值程序CONTACT,这是目前为止最为完善的滚动接触理论[15]. 该理论将轮轨接触问题转化为数学规划问题,利用Bossinesq-Cerruti公式可得轮轨滚动接触离散模型为:
{minCpJj=12pIiAIiJjpJj+[(g0J−q)pJz+(WJτ−uJτ)pJτ]A0,pJz⩾0,|pJτ|⩽bJ, ∀x∈Ac, (3) 式中:
C 为接触问题的余能;I、J 为接触斑上离散单元的编号;i、j、z 为X、Y、Z 坐标轴的方向;τ 为切向,对应XOY 平面;pIi 为单元I 沿方向i 的接触应力;pJj、pJz、pJτ 分别为单元J上沿Y轴、Z轴及方向τ 的接触应力;AIiJj 为影响函数,表示单元J方向j上单位力引起的单元I方向i上位移;g0J 为单元J中心处轮轨变形前接触面间的法向间隙;q为渗透量;WJτ 为单元J中心处从前时刻t′ 到当前时刻t总的刚性滑动量;uJτ 为轮轨滚动接触过程中时刻t时单元J中心处的弹性位移差;A0 为单元J的面积;bJ 为是单元J中心处Coulomb极限摩擦力;x为轮轨在X轴上坐标;Ac 为接触斑内坐标集合;车轮和尖轨的接触为异型接触,容易产生接触疲劳现象. 本文采用基于安定图模型的表面滚动接触疲劳因子预测不同工况下的轮轨表面滚动接触疲劳伤损. 轮轨接触斑中任一点的轮轨表面滚动接触疲劳系数
fI(•) 定义为{fI(x,y)=ft−kpz(x,y),ft=√(px(x,y))2+(py(x,y))2pz(x,y), (4) 式中:
ft 为牵引系数;y、z 分别为接触斑在Y、Z 轴上坐标值;px(•) 、py(•) 、pz(•) 分别为轮轨纵向、横向和法向接触应力;k 为材料自身的剪切屈服强度,取350 MPa.式(4)是根据赫兹接触理论得到,对于非赫兹接触问题,依据条带法,认为条带中间位置符合赫兹假设条件,从而将式(4)的应用扩展至非赫兹接触的范围,求得接触斑内任一单元的轮轨表面滚动接触疲劳系数,轮轨表面滚动接触疲劳因子定义为接触斑内滚动接触疲劳系数的最大值.
3. 计算结果与分析
3.1 轮轨接触点
轮轨几何接触是解释轮轨接触关系的基础[4]. 利用迹线法原理计算轮轨接触点,计算参数有:轮背距1353 mm,名义滚动圆半径460 mm,轮背到名义滚动圆处水平距离为70 mm,以Y轴正向为正,轮对横移量取−12 ~ 12 mm,横移量步长取0.5 mm. 以往研究表明,加宽轨距和改变轨底坡往往能够改善轮轨匹配关系[3-5],因此,对轨距1433、1435、1437、1439 mm,轨底坡1/10、1/20、1/30、1/40、1/50进行计算分析.
不同轨距和轨底坡参数下的接触点分布如图4和图5所示. 由图4可见:随着轨距的变化,轮轨接触点分布存在明显的差异;随着轨距的增加,轮轨接触点更多地分布在尖轨顶部,不易发生轮缘接触;随着轨距的增大,轮轨接触点由基本轨转移到尖轨所需要的横移量增大,也即轮载过渡延后,有利于减小尖轨受力,但会增大轮载过渡时轮轨接触点的跳跃,横向不平顺增大;在相同轨距下,随着车轮的磨耗,轮轨接触点更多地分布在基本轨,轮载过渡位置延后,横向不平顺增大.
由图5可见:不同轨底坡条件下,轮轨接触点的分布存在较大差异. 当轨底坡为1/10、1/20时,发生轮缘接触需要的轮对横移量最大,但和其余轨底坡相差较小;当轨底坡为1/10、1/20时,轮载过渡延后,横向不平顺增大;不同轨底坡下,发生轮载过渡时所需要的横移量最大为1/10,其次为1/20和1/30,最后为1/40和1/50.
3.2 等效锥度
等效锥度作为轮轨接触线性化指标,被广泛用于表征轮轨接触几何特征. 计算等效锥度的方法有简化法、谐波法和UIC 519法,UIC 519法采用轮对周期运动的假设,计算更准确[16],本文通过UIC 519标准[14]计算等效锥度.
在不同轨距下,磨耗车轮和60N关键断面匹配时的等效锥度见图6. 增大轨距有利于减小车轮踏面的等效锥度,从而提升列车过岔的平稳性;当轨距为1439 mm时,等效锥度基本上小于0.05,横移量相同条件下,其等效锥度是轨距1435 mm时的20%左右,极大减小了车轮等效锥度,提升了列车过岔平稳性;轨距变化对运营里程为25万km车轮的等效锥度影响较小.
图7为不同轨底坡和不同磨耗车轮条件下的等效锥度. 对于标准车轮和运营里程25万km车轮,1/30、1/40、1/50轨底坡条件下的等效锥度相差较小;当轨底坡为1/10、1/20,在横移量小于6 mm时,和轨底坡为1/30、1/40、1/50相比,等效锥度普遍较大,车辆过岔平稳性较差,横移量大于8 mm时,结果相反;在车轮运营里程为15万km时,1/30轨底坡条件下的等效锥度较小,列车过岔平稳性较好.
3.3 接触斑面积
相同法向轮轨力作用下,轮轨接触斑面积越大,其接触应力越小. 轮轨接触应力是影响轮轨磨耗和接触疲劳的重要因素. 利用非赫兹接触理论,计算轮轨接触斑面积和滚动接触疲劳因子,分析轨道参数取值对轮轨静力学接触行为的影响. 不同轨道参数工况下的参数取值:轴重14 t,单侧车轮轮心施加一半轴重,剪切模量82 GPa,泊松比0.28,摩擦系数0.3,划分网格单元0.2 mm × 0.2 mm,轮对横移量取0 ~ 12 mm. 选取右轮轨作为分析对象.
不同轨距和不同磨耗车轮条件下的接触斑面积如图8所示. 当横移量小于9 mm时,轨距和接触斑面积参数呈负相关,而当轮对横移量大于9 mm时,轨距和接触斑面积参数大致呈正相关;轮对横移量小于8 mm时,不同轨距条件下的轮轨接触斑面积相差较小,但当轮对横移量大于8 mm时,轨距越大,轮轨接触斑面积普遍越大;在轮对横移量大于8 mm时,轨距加宽有利于增大轮轨接触斑面积,减小轮轨接触应力.
不同轨底坡和不同磨耗车轮条件下的接触斑面积如图9所示. 1/10轨底坡下的轮轨接触斑面积普遍较小,且随车轮磨耗量的增大,其接触斑面积远小于其余轨底坡条件下的,说明1/10轨底坡下的轮轨接触力学性能较差,且会随着车轮磨耗变得更差. 由标准轮轨条件下,轨底坡和接触斑面积大致呈负相关. 随车轮磨耗量增大,1/30 ~ 1/50轨底坡的接触斑面积相差较小.
3.4 滚动接触疲劳因子
利用式(4)计算轮轨滚动接触疲劳因子,分析轮轨在不同轨道参数下的接触疲劳现象. 不同轨道参数下的表面滚动接触疲劳因子分布如图10、11所示. 由图10可知:在轮对横移量小于7 mm时,不同轨距下的表面滚动接触疲劳因子相差较小,在横移量大于 7 mm时,增大轨距可延缓轮轨表面进入滚动接触疲劳区;车轮运营里程达到15万km,在横移量为3 ~ 7 mm时,易发生轮轨间的两点接触,导致轮轨表面材料易进入疲劳区,从而产生疲劳破坏;增大轨距有利于减少轮轨材料出现接触疲劳现象,延长轮轨服役寿命.
由图11可知:标准轮轨条件下,轨底坡和滚动接触疲劳因子大致呈正相关;车轮磨耗导致轮轨滚动接触疲劳因子减小,原因是轮轨过渡位置延后,有利于减少尖轨磨耗,在横移量大于8 mm时,1/10和1/20轨底坡能延缓轮轨材料进入滚动接触疲劳区,但在横移量较小时,其滚动接触疲劳因子明显较大;1/10和1/20轨底坡下的轮轨滚动接触疲劳因子普遍较大,易引起轮轨材料进入滚动接触疲劳区,降低轮轨材料的使用寿命,轨底坡为1/30、1/40、1/50时,滚动接触疲劳因子相差较小.
随着车轮磨耗加深,道岔区的轮轨接触行为变得更为复杂. 运营里程15万km车轮和60N钢轨高速道岔在轨距1439 mm条件下发生轮轨两点接触时的接触斑分布如图12所示,接触斑中箭头指向表示切向应力合力的方向,箭头长短代表合力大小,滑动区轮轨表面间发生了相对滑移,黏着区轮轨间存在滑动趋势,但没有相对滑移. 由图12(a)可见,轨距1439 mm,轮对横移量5 ~ 7 mm时,发生了轮载过渡行为,随着横移量的增大,轮载更多由尖轨承载. 由图12(b)可见:在轮对横移量5 ~ 6 mm时,尖轨上存在较大的滑动区,在轮对横移量7 mm时,基本轨上存在较大的滑动区,两点接触导致了较大滑动区的存在,易导致车轮在钢轨上空转,引起道岔钢轨的磨损,从而缩短道岔区钢轨服役寿命. 需要指出的是,法向接触应力较小时同样可能存在较大滑动区,如图12中轮对横移量为5、7 mm时所示.
3.5 60N和CHN60钢轨高速道岔轮轨接触几何对比
为对比60N和CHN60钢轨高速道岔的轮轨接触几何行为,在轨距1439 mm和轨底坡1/30、1/40、1/50条件下,计算60N、CHN60钢轨高速道岔35 mm顶宽关键断面与三种磨耗车轮的等效锥度,如图13.
由图13可知:相同轨底坡和轮对横移量条件下,相比于CHN60钢轨,60N钢轨高速道岔区的等效锥度更小,列车过岔平稳性更优;车轮运营里程达到25万km时,当车轮横移量小于5 mm时,60N钢轨高速道岔区等效锥度普遍更小,具有较好的过岔平稳性,车轮横移量大于5 mm时,60N钢轨高速道岔区等效锥度普遍更大,轮对的对中性能更好;随车轮磨耗量增加,等效锥度整体呈增加趋势.
4. 结 论
1) 轨距对轮轨接触行为影响较大. 轨距加宽有利于减少轮缘接触,较大程度减小等效锥度,提升列车过岔的平稳性;轨距加宽可减小轮对横移量大于8 mm时的轮轨接触应力和表面滚动接触疲劳因子,减少轮轨材料发生接触疲劳,延长尖轨使用寿命. 轮对横移量小于8 mm时,轨距加宽对运营里程25万km车轮的接触性能影响较小.
2) 轨底坡对轮轨接触行为影响较大. 标准轮轨条件下,轨底坡和接触斑面积大致呈负相关,与滚动接触疲劳因子大致呈正相关. 轨底坡为1/10和1/20时,轮载过渡位置延后,横向不平顺增大,车轮横移量小于6 mm时,等效锥度普遍较大;1/10和1/20轨底坡下的接触斑面积普遍较小,轮轨滚动接触疲劳因子普遍较大,较易引起轮轨材料进入滚动接触疲劳区,降低轮轨材料的使用寿命,且1/10轨底坡对车轮磨耗的适应性较差. 轨底坡为1/30、1/40、1/50时,轮轨接触参数相差较小,匹配性能较优.
3) 和CHN60钢轨高速道岔相比,60N钢轨的等效锥度普遍更小,列车过岔平稳性更优;车轮运营里程为25万km时,当轮对横移量小于5 mm时,60N钢轨的等效锥度普遍更小,当轮对横移量大于5 mm时,结果相反.
4) 车轮磨耗易引起道岔区轮轨间的两点接触,在较小轮轨法向接触应力下,接触斑上也易出现较大滑动区,导致车轮空转,引起钢轨伤损. 随车轮磨耗量增加,轮轨间等效锥度整体呈增加趋势.
致谢:京沪高速铁路股份有限公司科技研究项目(京沪科研-2020-11)的支持.
-
[1] 王平,陈嵘,徐井芒,等. 高速铁路道岔系统理论与工程实践研究综述[J]. 西南交通大学学报,2016,51(2): 357-372. doi: 10.3969/j.issn.0258-2724.2016.02.015WANG Ping, CHEN Rong, XU Jingmang, et al. Theories and engineering practices of high-speed railway turnout system: survey and review[J]. Journal of Southwest Jiaotong University, 2016, 51(2): 357-372. doi: 10.3969/j.issn.0258-2724.2016.02.015 [2] 杜星,陶功权,杨城,等. 轨底坡变化对高速车辆运行行为的影响[J]. 西南交通大学学报,2022,57(2): 286-294.DU Xing, TAO Gongquan, YANG Cheng, et al. Influence of different rail cants change on dynamical characteristics of high-speed railway vehicles[J]. Journal of Southwest Jiaotong University, 2022, 57(2): 286-294. [3] 钱瑶,王平,苏谦,等. 轨底坡对高速铁路轮轨接触行为影响分析[J]. 铁道工程学报,2018,35(3): 18-25. doi: 10.3969/j.issn.1006-2106.2018.03.004QIAN Yao, WANG Ping, SU Qian, et al. Effect analysis of rail cant on the wheel-rail contact behavior of high-speed railway[J]. Journal of Railway Engineering Society, 2018, 35(3): 18-25. doi: 10.3969/j.issn.1006-2106.2018.03.004 [4] 陈嵘,温静,于浩,等. 地铁线路轨距对轮轨接触行为的影响[J]. 中南大学学报(自然科学版),2020,51(3): 824-831. doi: 10.11817/j.issn.1672-7207.2020.03.028CHEN Rong, WEN Jing, YU Hao, et al. Influence of rail gauge on wheel-rail contact behavior of metro line[J]. Journal of Central South University (Science and Technology), 2020, 51(3): 824-831. doi: 10.11817/j.issn.1672-7207.2020.03.028 [5] 陈嵘,温静,李博,等. 考虑非对称轨底坡的轮轨滚动接触应力分析[J]. 铁道工程学报,2019,36(5): 13-19,70. doi: 10.3969/j.issn.1006-2106.2019.05.003CHEN Rong, WEN Jing, LI Bo, et al. Analysis of wheel/rail rolling contact stress considering asymmetric rail cant[J]. Journal of Railway Engineering Society, 2019, 36(5): 13-19,70. doi: 10.3969/j.issn.1006-2106.2019.05.003 [6] CUI D B, LI L, JIN X S, et al. Wheel-rail profiles matching design considering railway track parameters[J]. Chinese Journal of Mechanical Engineering, 2010, 23(4): 410-417. doi: 10.3901/CJME.2010.04.410 [7] SÁNCHEZ R A, SANJUÁN E L, BRAVO J L. Experimental validation of track inspection trolley using a rigorous self-checking procedure[J]. Journal of Surveying Engineering, 2020, 146(3): 0000315.1-0000315.8. [8] 李超,张军,李霞,等. 动态轨距优化技术在重载道岔上的应用[J]. 大连交通大学学报,2015,36(S1): 55-61.LI Chao, ZHANG Jun, LI Xia, et al. Research on dynamic gauge optimization of heavy haul switch rail[J]. Journal of Dalian Jiaotong University, 2015, 36(S1): 55-61. [9] YE Y G, SUN Y. Reducing wheel wear from the perspective of rail track layout optimization[J]. Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multi-Body Dynamics, 2021, 235(2): 217-234. doi: 10.1177/1464419320956831 [10] 闫正,陈嘉胤,徐井芒,等. 不同车轮踏面与高速60N钢轨道岔静态接触特性研究[J]. 中南大学学报(自然科学版),2021,52(4): 1358-1370. doi: 10.11817/j.issn.1672-7207.2021.04.032YAN Zheng, CHEN Jiayin, XU Jingmang, et al. Study of static contact properties of diverse wheel treads and high-speed 60N rail turnout[J]. Journal of Central South University (Science and Technology), 2021, 52(4): 1358-1370. doi: 10.11817/j.issn.1672-7207.2021.04.032 [11] 陈嵘,方嘉晟,汪鑫,等. 车轮型面演变对高速道岔区轮轨接触行为影响分析[J]. 铁道学报,2019,41(5): 101-108. doi: 10.3969/j.issn.1001-8360.2019.05.012CHEN Rong, FANG Jiasheng, WANG Xin, et al. Influence of wheel profile evolution on wheel-rail contact behavior in high-speed turnout area[J]. Journal of the China Railway Society, 2019, 41(5): 101-108. doi: 10.3969/j.issn.1001-8360.2019.05.012 [12] 王晨,罗世辉,邬平波,等. 动车组踏面凹型磨耗对车辆稳定性的影响[J]. 西南交通大学学报,2021,56(1): 84-91.WANG Chen, LUO Shihui, WU Pingbo, et al. Effect of hollow worn tread of electric multiple units on vehicle stability[J]. Journal of Southwest Jiaotong University, 2021, 56(1): 84-91. [13] British Standards Institution. BS EN 15302: 2008 Railway applications—method for determining the equivalent conicity[S]. London: British Standards Institution, 2008. [14] International Union of Railways. UIC code 519 method for determining the equivalent conicity[S]. Paris: International Union of Railways, 2004. [15] 马晓川,王平,徐井芒,等. 铁路道岔轮轨非赫兹滚动接触算法对比与分析[J]. 机械工程学报,2019,55(18): 95-103. doi: 10.3901/JME.2019.18.095MA Xiaochuan, WANG Ping, XU Jingmang, et al. Analysis and comparison of different wheel-rail non-hertzian rolling contact approaches in railway turnout[J]. Journal of Mechanical Engineering, 2019, 55(18): 95-103. doi: 10.3901/JME.2019.18.095 [16] 干锋,戴焕云,高浩,等. 铁道车辆不同踏面等效锥度和轮轨接触关系计算[J]. 铁道学报,2013,35(9): 19-24. doi: 10.3969/j.issn.1001-8360.2013.09.004GAN Feng, DAI Huanyun, GAO Hao, et al. Calculation of equivalent conicity and wheel-rail contact relationship of different railway vehicle treads[J]. Journal of the China Railway Society, 2013, 35(9): 19-24. doi: 10.3969/j.issn.1001-8360.2013.09.004 -