基于离散元法的钢轨力学性能定性分析 Qualitative analysis of rail mechanical properties based on discrete element method 1 2 corresp first-author 陈宪麦，博士，副教授，从事轨道结构、轨道平顺性等研究；E-mail：chenxianmai@csu.edu.cn 陈宪麦，博士，副教授，从事轨道结构、轨道平顺性等研究；E-mail：chenxianmai@csu.edu.cn chenxianmai@csu.edu.cn 1 3 3 3 中南大学 土木工程学院，湖南 长沙，410075 中南大学 土木工程学院，湖南 长沙，410075 School of Civil Engineering, Central South University, Changsha 410075, China School of Civil Engineering, Central South University, Changsha 410075, China 重载铁路工程结构教育部重点实验室，湖南 长沙，410075 重载铁路工程结构教育部重点实验室，湖南 长沙，410075 Key Laboratory of Heavy Haul Railway Engineering Structure, Ministry of Education, Changsha 410075, China Key Laboratory of Heavy Haul Railway Engineering Structure, Ministry of Education, Changsha 410075, China 中国铁道科学研究院集团有限公司 基础设施检测研究所，北京，100081 中国铁道科学研究院集团有限公司 基础设施检测研究所，北京，100081 Infrastructure Inspection Research Institute, China Academy of Railway Sciences Group Co. Ltd., Beijing 100081, China Infrastructure Inspection Research Institute, China Academy of Railway Sciences Group Co. Ltd., Beijing 100081, China 陈灿华 为了更好地了解钢轨的本质特征，分析钢轨在列车荷载作用下的响应以及裂纹的萌生与变化规律，采用离散元法构建钢轨静力学分析模型并研究其静力学特性。通过抗拉强度试验仿真，验证离散元力学参数的可行性；在此基础上，对横纵断面钢轨的最不利状态进行静力加载分析，并进一步探究车轮垂向位移对钢轨裂纹演变规律的影响。研究结果表明：静力加载时，轮轨接触力主要由作用点附近的钢轨颗粒承担，距离作用点较远的钢轨内部颗粒承担的荷载较小；随车轮向下移动，初始阶段钢轨轨头出现小量倾斜，接触面增大而接触应力减小；随车轮位移增加，钢轨轨头和轨腰逐渐出现裂纹，并且裂纹不断扩展；在纵断面钢轨接头处，随车轮垂向位移增加，钢轨承担的荷载由轨头转向轨腰，轨底所承担的轮轨接触力较小；当车轮通过钢轨接头时，由于其对轨缝两侧钢轨分配的荷载不均匀，接头容易产生病害。 In order to better understand the essential characteristics of rails, analyze the response of rails under train load and the generation and change of cracks, the discrete element method was used to construct rail statics analysis model and study its statics characteristics. Through the tensile strength test simulation, the feasibility of the discrete element mechanical parameters was verified, the most unfavorable conditions of the transverse and longitudinal sections of the rail were analyzed by static loading, and the influence of the vertical displacement of the wheel on the crack evolution law of the rail was explored. The results show that the wheel-rail contact force during static loading is mainly borne by the rail particles near the point of application, and the internal particles of the rail farther from the point of application bear less load. As the wheel moves downwards in the initial stage, the rail head tiltes slightly to increase the contact surface and reduce the contact stress. With the increase of the wheel displacement, the rail head and rail waist gradually display cracks, and the cracks continue to expand. At longitudinal section rail joints, as the vertical displacement of the wheel increases, the load borne by the rail is shifted from the rail head to the rail waist, and the wheel-rail contact force borne by the bottom of the rail is small. When the wheel passes through the rail joint, the load is distributed to the rails on both sides of the rail seam unevenly, which makes the joints prone produce disease. 离散元 钢轨伤损 力学性能 裂纹扩展 discrete element rail damage mechanical properties crack growth 近年来，重载铁路、高速铁路迅速发展，年运量和列车运行频次也不断提高。钢轨作为轨道结构的重要组成部分，其功能主要是引导列车车轮前进，将承受的荷载传至轨枕，有时还兼有轨道电路的功能。随着列车荷载的往复作用，钢轨上逐渐出现裂纹、轨头核伤、轨头压溃、轨头剥离、折断等病害 [ 1 1 - 3 3 1-3 ] 。这些不同程度的钢轨伤损将会影响旅客的舒适性和列车运行的安全性，也会影响轮轨接触，产生更大的轮轨作用力，进一步增大钢轨裂纹、折断等伤损的可能性 [ 4 4 - 5 5 4-5 ] 。为了了解钢轨的性能，许多研究人员对钢轨进行了探究。马晓川等 [ 6 6 ] 运用近场动力学损伤理论研究了钢轨裂纹的萌生和发展规律，得出萌生于钢轨轨头的裂纹启裂于钢轨表面下约2 mm处；HU等 [ 7 7 ] 对不同硬度的钢轨磨耗情况进行了研究，分析了体积硬度、试验后硬度、硬化率等对钢轨磨损率的影响；胡文博等 [ 8 8 ] 采用深度学习方法对钢轨的伤损进行了探测，并对钢轨不同位置的伤损进行了分类；GUI等 [ 9 9 ] 研究了微观结构对贝氏体钢轨的滚动接触疲劳性能的影响；ZHAO等 [ 10 10 ] 研究了钢轨表面的球形凹痕对其微观结构和滚动接触的影响；PAN等 [ 11 11 ] 研究了钢材的微观结构对其服役性能的影响，指出钢材的使用寿命受微观结构的影响很大。 从宏观层面考虑，钢轨属于连续体，可以用有限元方法进行模拟 [ 12 12 - 15 15 12-15 ] ，这也是钢轨常用的数值模拟方式。有限元模型假设钢轨是一个各向同性的连续体，采用有限元法可以模拟钢轨在不同列车荷载作用下的受力情况，根据材料的属性，当应力达到钢轨屈服强度时，该处的应力不再增大，增加的荷载由其余未达到屈服强度的区域承担。而实际上，钢轨是由许多微细观颗粒组成的，当钢轨承受的荷载较大时，在其内部会产生细小裂纹，这些裂纹会随着外部荷载的增加而扩展，采用有限元模型无法模拟钢轨内部裂纹的扩展情况。为此，本文采用离散元法，将细小的颗粒通过平行黏结组成钢轨，并根据钢轨的力学性能对细小颗粒添加接触模型和相关的力学参数。当离散元模型钢轨承受的荷载超过屈服强度时，颗粒之间的接触力链会发生断裂，从而产生裂纹，这些裂纹会随着车轮垂向位移的增大而不断发展，直至钢轨失效。 对于钢轨微观性能的研究，已有学者运用元胞自动机的方法探究钢轨和车轮之间的磨损及其演变规律 [ 16 16 ] ，但是没有考虑钢轨内部的损伤机理。元胞自动机是一种时间和空间都离散的数值模拟方法，主要是以网格(包括三边、四边、六边网格)作为元胞单元，根据元胞的当前状态及其邻居状况确定下一时刻该元胞状态的动力学函数，下一步结果是以上一步结果为基础，可以预测复杂问题的发展趋势。而离散元法是一种在时间上连续、空间上离散的模拟方法，本文运用细小的球体颗粒作为模拟钢轨内部的小单元，通过对颗粒间施加平行黏结模型和黏结强度使其成为一个整体。本文拟构建基于离散元法的钢轨静力学分析模型，该模型由细小球体颗粒构成，并设置颗粒之间的接触模型和黏结强度，以保证模型与实际的钢轨性能接近，同时又可以从微观和细观的角度了解和分析钢轨内部的作用机理，对掌握和优化钢轨性能具有参考价值和指导意义。 1 钢轨离散元模型的建立 1 1.1 　物理力学参数的选定 2 钢材由金相晶体构成，其内部的晶体结构主要分为体心立方体、面心立方体和密排六方体3种 [ 17 17 - 18 18 17-18 ] 。为了近似分析钢轨内部晶体之间的荷载传递和相互作用规律，通过离散元法模拟钢轨，采用小粒径(如1.0 mm)球体构建钢轨，用小粒径球体颗粒近似模拟钢材中的晶体晶粒，并设置颗粒之间的接触模型为平行黏结模型，近似模拟晶体晶粒之间的相互作用关系，以便于分析颗粒间的荷载传递情况。 现代钢轨的生产流程主要为精炼、精轧、精整、质量自动检测、长尺化生产等。与传统的模铸加工普通孔型法轧制工艺相比，现代的生产工艺不仅可以保证钢轨的几何尺寸更加精确，而且可以使钢轨具有良好的内在质量和表面质量 [ 19 19 ] 。本文不考虑生产工艺导致的不利因素，以U71Mn钢轨为研究对象，参考既有的规范和文献[ 20 20 - 23 23 20-23 ]，选定的离散元模型细观力学参数如 表1 表1 所示。 参数 取值 钢轨的密度/(kg·m -3 ) 7 850 钢轨的抗拉强度/Pa 1.0×10 9 钢轨的屈服强度/Pa 8.0×10 8 钢轨的弹性模量/Pa 2.1×10 11 钢轨颗粒间的法向接触刚度/(N·m -1 ) 5.0×10 13 钢轨颗粒间的切向接触刚度/(N·m -1 ) 5.0×10 13 钢轨颗粒间的法向黏结刚度/(N·m -1 ) 5.0×10 13 钢轨颗粒间的切向黏结刚度/(N·m -1 ) 5.0×10 13 钢轨颗粒间的黏结强度/(N·m -2 ) 1.0×10 20 钢轨的摩擦因数 0.5 1.2 　模型参数的验证 2 为了保证离散元法建立的钢轨模型具有可靠性，对钢轨抗拉强度进行试验仿真。根据相关标准 [ 24 24 ] ，对宽度 d 为10 mm、长度为50 mm的钢材试样进行拉伸试验和离散元仿真，拉伸前后的钢材试样如 图1 图1 所示，经轴向拉伸后得到的钢材应力-应变曲线如 图2 图2 所示。 (a) 拉伸前的钢轨试样；(b) 拉伸后的钢轨试样； (c) 拉伸后的钢轨试样接触力分布 从 图1 图1 可以看出：拉伸前钢材内部接触力分布均匀，钢材各向同性；钢材拉伸后，在中间部位的受力最大，并且在该处出现了部分接触力链的断裂情况。 图1 图1 (b)中，蓝色球体颗粒表示构成钢轨的微观颗粒，绿色短线代表拉力，红色短线代表压力。从 图1 图1 (b)可以看出：在钢材中间一定区域内，钢材试样主要受到拉力的作用，而在其他区域受到拉力和压力的共同作用；随着轴向拉伸不断进行，试样的中间部位由于所受拉力过大而出现接触力链断裂的情况，此时，钢材试样达到屈服状态。 工程上规定 [ 25 25 ] ，将试样产生0.2%塑性变形时对应的应力作为该材料的屈服强度。由 图2 图2 可知：钢材试样应变为0.2%时的拉应力为810 MPa，与本文选定的屈服强度非常相近，并且根据文献[ 26 26 ]的应力-应变曲线可以验证本文选取的离散元力学参数是可信的，可以在此基础上继续进行模拟研究。 1.3 　钢轨的初始模型 2 为了探究钢轨在不同工况下的服役性能及其损伤特性，分别对钢轨横截面和钢轨接头纵断面进行2D离散元仿真。采用的主要建模方法是：将钢轨模拟为由ball组成的颗粒簇，并设定各个球体颗粒之间的接触属性；对于车轮的模拟方法，在分析静力作用时，为了便于加载，用clump单元模拟钢轨横截面上的车轮，用ball单元模拟钢轨纵断面上的车轮。在研究钢轨伤损时，为加快程序运行速率，用wall单元模拟车轮。本文以钢轨为主要研究对象，而不考虑车轮的变形，在钢轨横纵断面均将车轮视为刚体。 钢轨横截面与接头纵断面的仿真初始模型如 图3 图3 所示。其中，车轮选用LM磨耗型车轮踏面，据文献[ 27 27 ]，磨耗型车轮踏面对应的轮轨接触应力限值为1 600 MPa。 (a) 钢轨横截面；(b) 钢轨接头纵断面 通过对车轮施加荷载以实现对钢轨间接加载。由于施加的荷载是轮轨接触力，其中包含了车轮重力，所以，模拟时车轮的直径对加载力没有影响。为使图形显示更加直观， 图3 图3 (a)中的车轮仅展示部分区域， 图3 图3 (b)中的车轮直径缩小为原来的1/10。 2 钢轨横断面的受力特性 1 2.1 　静力分析 2 为了避免由于突然施加较大荷载而产生的冲击作用，对钢轨进行分级静力加载，分级加载的次数为100次，直至达到目标荷载为止。据文献[ 28 28 - 29 29 28-29 ]，并取轮轨最易受损的加载条件，最终选用的轮轨横向作用力为10 kN，垂向作用力为150 kN。钢轨横截面接触力分级静力加载曲线如 图4 图4 所示。 静力加载后钢轨横截面的接触力分布情况如 图5 图5 所示。 由 图5 图5 可知：在静力荷载作用下，轮轨接触力主要由轮轨作用点附近的颗粒承担，离接触点较远的钢轨内部颗粒受静载力作用较小。 2.2 　钢轨伤损 2 随着列车的往复运行，轮轨间的接触面由于受车轮踏面缺陷、钢轨表面不平顺、轮轨之间的冲击、道床沉降等不利因素的影响，钢轨顶面的平顺性逐渐下降。对车轮施加竖直向下的位移作用，直至垂向位移达到10 mm，以便于分析钢轨在车轮不同垂向位移下的伤损及裂纹分布，从而间接反映轨道不平顺对轮轨接触力的影响。横截面钢轨的伤损演变情况如 图6 图6 所示。建立轮轨离散元初始模型时，为了避免初始的轮轨接触力对后续车轮移动产生的轮轨接触力的影响，设置轮轨在初始阶段并没有绝对接触(即初始轮轨接触力为 0 kN)，即车轮与钢轨轨头之间有较小间隙；后续随车轮竖向移动，轮轨间才逐渐产生接触。 (a) 初始状态；(b) 车轮垂向移动2 mm； (c) 车轮垂向移动4 mm；(d) 车轮垂向移动6 mm； (e) 车轮垂向移动8 mm；(f) 车轮垂向移动10 mm 从 图6 图6 可以发现：随着车轮垂向位移增大，钢轨的伤损程度逐渐加大，在车轮垂向位移为2 mm时，主要在轮轨接触作用点附近产生较大的应力，几乎没有裂纹产生；当车轮垂向位移为4 mm时，钢轨轨头产生较小程度的倾斜，可增大轮轨接触面，减少接触应力；当垂向位移为6 mm时，轨头处产生明显的裂纹；当垂向位移为8 mm时，轨头上颚的裂纹即将贯穿至钢轨顶面，并且在轨腰位置出现部分伤损现象；当车轮垂向位移为10 mm时，轨头的裂纹数量没有明显增多，但裂纹宽度逐渐增大，并在轨头和轨腰出现掉块现象，钢轨基本丧失其承载能力。 轮轨之间接触力随车轮垂向位移的变化趋势如 图7 图7 所示，据 图7 图7 可对钢轨伤损的力学性能进行定量分析。 从 图7 图7 可见：由于在初始时轮轨之间没有接触，随着车轮的竖向移动逐渐产生接触力，所以，当车轮垂向位移小于1 mm时，轮轨间的接触力较小；当车轮垂向位移超过1 mm时，轮轨接触力随着车轮垂向位移增大而显著变大；当位移增大到 4 mm时，轮轨接触力达到最大，而后随着垂向位移增大，轮轨接触力逐渐减小；当位移达到8 mm时，轮轨接触力又出现了小幅度提高。这是因为当垂向位移达到4 mm后，钢轨轨头出现了细小裂纹，并且裂纹数量随着位移增大而增多，从而导致其承载力有所下降；当位移达到8 mm时，钢轨出现了较明显塑性变形，钢轨与车轮的接触面积增大，其后的轮轨接触力也有所增加。 3 钢轨接头纵断面的受力特性 1 3.1 　静力分析 2 钢轨接头是轨道结构的薄弱环节，故选取接头处一定长度范围内的钢轨进行研究。建立离散元钢轨纵向模型时，考虑到钢轨的轨头、轨腰、轨底3部分的横截面面积不相等，若直接采用钢轨的密度进行属性赋值则存在较大的偏差，因此，对钢轨密度赋值时，将轨头、轨腰、轨底3部分的颗粒密度按照相应部位的面积比例进行分配。钢轨纵向长度的选取以包含实际接头孔洞为原则，最终选取长度为600 mm的钢轨作为研究对象。为了用离散元法分析钢轨在压力作用下的力学性能，本文在纵向钢轨试样底部添加1个与其等长的wall单元，仅考虑钢轨在轮载作用下的压应力分布。 对于钢轨纵断面，分别取接头最左端、第1个螺栓孔、第2个螺栓孔、第3个螺栓孔、最右端5个位置添加车轮，以便于探究不同位置的车轮对钢轨受力性能的影响。车轮在纵断面钢轨的分布如 图8 图8 所示。 对钢轨纵断面进行静力加载时，采取与钢轨横截面静力加载类似的方式，均采用分级加载。静力加载后，车轮在不同纵向位置时的轮轨接触力见 表2 表2 。 车轮位置 轮轨接触力/kN Wheel_1 74.118 Wheel_2 34.799 Wheel_3 34.542 Wheel_4 34.477 Wheel_5 74.368 分析 表2 表2 可以发现：当车轮位于钢轨孔洞位置上方时，Wheel_2、Wheel_3、Wheel_4三者的轮轨接触力相差不大；而当车轮位于研究对象的两端时，轮轨接触力较大，约为孔洞位置轮轨作用力的2倍。出现该结果的原因是，位于钢轨试样两端的车轮由于其一边没有钢轨处于悬空状态，仅由车轮另一边钢轨承担传递的荷载，从而产生了较大的接触力。 由上述研究可知，接头处容易出现钢轨低头压溃等病害，其原因是：一方面，钢轨接头的螺栓孔减小了钢轨的截面，其承载力降低，同时，在孔洞处易产生应力集中现象；另一方面，钢轨接头存在轨缝，当列车运行到轨缝时，该处的钢轨由于部分悬空而将荷载不均匀地传递至轨缝两侧的钢轨，从而导致两侧钢轨由于受力不同而出现不同程度的沉降，使得轨道结构的平顺性及其服役性能进一步降低。 3.2 　钢轨接头伤损 2 当车轮的垂向位移达到10 mm时，不同工况下的钢轨内部颗粒接触力云图如 图9 图9 所示。 (a) Wheel_3位置的钢轨接触力云图；(b) 不同车轮位置的钢轨轨头接触力局部图 图9 图9 (a)所示为车轮位于Wheel_3时钢轨内部颗粒间的接触力分布云图，可见在轮轨作用点附近的接触力较大，距离作用点较远处的接触力均较小。当车轮作用在不同位置时，作用点附近的接触力分布情况见 图9 图9 (b)。从 图9 图9 (b)可见：当车轮位于3个螺栓孔上方时，Wheel_2、Wheel_3、Wheel_4三者的钢轨颗粒接触力分布情况类似，并且沿车轮中心轴对称分布；而位于接头两端的车轮由于车轮荷载仅由一边的钢轨承担，导致其实际承载力与中间钢轨承载力相比较大，并且钢轨两端均是一边处于自由状态，钢轨内部颗粒在较大荷载且不受约束的状态下容易出现掉粒、剥落现象。 在车轮垂向下移时，分别对不同车轮位置的钢轨应力进行监测，得到钢轨轨头、轨腰、轨底的应力随车轮垂向位移的变化情况。由于车轮在不同位置时的应力变化情况相似，因此，只列出车轮位于Wheel_3时，钢轨不同部位的应力随车轮垂向位移的变化，如 图10 图10 所示。 由 图10 图10 可知：随着车轮垂直向下移动，轨头由于与车轮直接作用而产生较大的应力并且远大于轨腰和轨底的应力，轨腰部位的应力次之，轨底处的内部颗粒应力最小。 当车轮位于不同位置时，对应的钢轨应力变化如 图11 图11 所示，据 图11 图11 可分析车轮位置对钢轨不同部位的影响。 (a) 轨头；(b) 轨腰；(c) 轨底 由 图11 图11 (a)可知：当车轮位于Wheel_2、Wheel_3、Wheel_4时，车轮位置对轨头应力随车轮垂向位移的变化几乎没有影响，3条曲线基本重合；当车轮位于Wheel_1、Wheel_5时的轨头应力变化曲线也几乎重合，并且均小于Wheel_2、Wheel_3、Wheel_4对应的轨头应力。轨腰应力曲线变化情况与轨头的应力曲线变化情况类似，但随着车轮垂向位移增大，轨头应力的增大速率逐渐减小，而轨腰应力的增大速率不断提高。这说明随着车轮向下移动，钢轨承受的荷载由以轨头承受为主转换为以轨腰承受荷载为主。 从 图11 图11 可见钢轨主要由轨头承受荷载，其次是轨腰，轨底对承载的贡献较小。虽然不同车轮位置的轨底应力不相同，但轨底应力均很小，并且相差不大。 4 结论 1 1) 轮轨接触力随车轮垂向位移增大，出现先增大后减小再增大的变化；接触力减小的原因是轨头产生细小裂纹，之后再增大的原因是钢轨发生了塑性变形。 2) 作用在钢轨轨缝的车轮由于其一边处于部分悬空状态，导致轨缝两侧的钢轨受力不均匀，容易产生钢轨低头、折角、台阶等现象。 3) 当车轮分别位于钢轨接头处的3个钢轨螺栓孔上方时，钢轨应力相差不大；当车轮位于钢轨端部时，应力曲线重合；不同车轮位置的钢轨应力均随着车轮垂向位移增大而增大。 4) 轮轨间接触力主要由钢轨轨头承担，轨腰分担的荷载次之，轨底分担的荷载最小。随着车轮向下移动，轨头承担的荷载逐渐向轨腰移动和传递。 参考文献 陈秀方 , 娄平 . 轨道工程 [M]. 2版 . 北京 : 中国建筑工业出版社 , 2017 : 11 - 36 . 轨道工程 CHEN Xiufang , LOU Ping . Track engineering [M]. 2nd ed . Beijing : China Architecture & Building Press , 2017 : 11 - 36 . Track engineering 姚锋涛 . 基于超声回波信号的钢轨伤损识别与分类研究 [D]. 西安 : 西安理工大学 , 2020 : 1 - 7 . 基于超声回波信号的钢轨伤损识别与分类研究 YAO Fengtao . Identification and classification of rail damage based on ultrasonic echo signals [D]. Xi'an : Xi'an University of Technology , 2020 : 1 - 7 . Identification and classification of rail damage based on ultrasonic echo signals 汤雪扬 , 蔡小培 , 彭华 , 等 . 基于轮轨系统耦合振动的地铁钢轨波磨研究 [J]. 中南大学学报(自然科学版) , 2022 , 53 ( 8 ): 3232 - 3244 . 基于轮轨系统耦合振动的地铁钢轨波磨研究 TANG Xueyang , CAI Xiaopei , PENG Hua , et al . Research on rail wave wear in subway based on coupled vibration of wheel rail system [J]. Journal of Central South University(Science and Technology) , 2022 , 53 ( 8 ): 3232 - 3244 . Research on rail wave wear in subway based on coupled vibration of wheel rail system 崔旭浩 , 马超智 , 金锋 . 重载铁路钢轨波磨对有砟道床动力特性的影响 [J]. 华中科技大学学报(自然科学版) , 2021 , 49 ( 1 ): 110 - 115 . 重载铁路钢轨波磨对有砟道床动力特性的影响 CUI Xuhao , MA Chaozhi , JIN Feng . The effect of rail corrugation on dynamic characteristics of heavy haul railway ballast bed [J]. Journal of Huazhong University(Science and Technology) , 2021 , 49 ( 1 ): 110 - 115 . The effect of rail corrugation on dynamic characteristics of heavy haul railway ballast bed 方修洋 , 张昊楠 , 马东文 . 轮轨型面匹配对钢轨亚表面抗疲劳裂纹扩展影响 [J]. 中南大学学报(自然科学版) , 2022 , 53 ( 5 ): 1824 - 1833 . 轮轨型面匹配对钢轨亚表面抗疲劳裂纹扩展影响 FANG Xiuyang , ZHANG Haonan , MA Dongwen . The effect of wheel rail profile matching on the fatigue crack propagation resistance of steel rail sub surface [J]. Journal of Central South University(Science and Technology) , 2022 , 53 ( 5 ): 1824 - 1833 . The effect of wheel rail profile matching on the fatigue crack propagation resistance of steel rail sub surface 马晓川 , 刘林芽 , 冯青松 , 等 . 铁路钢轨裂纹萌生的键型近场动力学预测模型 [J]. 交通运输工程学报 , 2021 , 21 ( 3 ): 228 - 237 . 铁路钢轨裂纹萌生的键型近场动力学预测模型 MA Xiaochuan , LIU Linya , FENG Qingsong , et al . Prediction model of rail crack initiation using bond-based peridynamics theory [J]. Journal of Traffic and Transportation Engineering , 2021 , 21 ( 3 ): 228 - 237 . Prediction model of rail crack initiation using bond-based peridynamics theory HU Y , WATSON M , MAIORINO M , et al . Experimental study on wear properties of wheel and rail materials with different hardness values [J]. Wear , 2021 , 477 ( SI ): 1 - 10 . Experimental study on wear properties of wheel and rail materials with different hardness values 胡文博 , 邱实 , 许馨月 , 等 . 基于深度学习的钢轨伤损超声检测与分类 [J]. 铁道学报 , 2021 , 43 ( 4 ): 108 - 116 . 基于深度学习的钢轨伤损超声检测与分类 HU Wenbo , QIU Shi , XU Xinyue , et al . Ultrasonic detection and classification for internal defect of rail based on deep learning [J]. Journal of the China Railway Society , 2021 , 43 ( 4 ): 108 - 116 . Ultrasonic detection and classification for internal defect of rail based on deep learning GUI Xiaolu , WANG Kaikai , GAO Guhui , et al . Rolling contact fatigue of bainitic rail steels: the significance of microstructure [J]. Materials Science and Engineering: A , 2016 , 657 : 82 - 85 . Rolling contact fatigue of bainitic rail steels: the significance of microstructure ZHAO X J , GUO J , LIU Q Y , et al . Effect of spherical dents on microstructure evolution and rolling contact fatigue of wheel/rail materials [J]. Tribology International , 2018 , 127 : 520 - 532 . Effect of spherical dents on microstructure evolution and rolling contact fatigue of wheel/rail materials PAN Jinzhi , CHEN Lin , LIU Chunpeng , et al . Relationship between the microstructural evolution and wear behavior of U 71 Mn rail steel [J]. Journal of Materials Engineering and Performance , 2021 , 30 ( 2 ): 1090 - 1098 . Relationship between the microstructural evolution and wear behavior of U71Mn rail steel 崔晓璐 , 包鹏羽 , 陈佳欣 , 等 . 高速铁路制动区间钢轨摩擦自激振动研究 [J]. 西南交通大学学报 , 2023 , 58 ( 1 ): 141 - 149 . 高速铁路制动区间钢轨摩擦自激振动研究 CUI Xiaolu , BAO Pengyu , CHEN Jiaxin , et al . Study on rail friction self-excited vibration in braking section of high-speed railway [J]. Journal of Southwest Jiaotong University , 2023 , 58 ( 1 ): 141 - 149 . Study on rail friction self-excited vibration in braking section of high-speed railway 刘佳璐 , 阚前华 , 汪渊 , 等 . 高速列车钢轨硌伤致损机理有限元分析 [J]. 四川轻化工大学学报(自然科学版) , 2021 , 34 ( 2 ): 60 - 69 . 高速列车钢轨硌伤致损机理有限元分析 LIU Jialu , KAN Qianhua , WANG Yuan , et al . Finite element analysis of damage mechanism of high-speed train rails [J]. Journal of Sichuan University of Light Chemical Technology(Science and Technology) , 2021 , 34 ( 2 ): 60 - 69 . Finite element analysis of damage mechanism of high-speed train rails 张友鹏 , 王东 , 赵斌 , 等 . 基于管状导体模型钢轨内阻抗计算 [J]. 中南大学学报(自然科学版) , 2020 , 51 ( 10 ): 2988 - 2997 . 基于管状导体模型钢轨内阻抗计算 ZHANG Youpeng , WANG Dong , ZHAO Bin , et al . Calculation of internal impedance of rail based on tubular conductor model [J]. Journal of Central South University(Science and Technology) , 2020 , 51 ( 10 ): 2988 - 2997 . Calculation of internal impedance of rail based on tubular conductor model 侯广友 . 基于压磨板的钢轨开式砂带打磨接触理论建模与仿真分析 [D]. 北京 : 北京交通大学 , 2020 : 1 - 40 . 基于压磨板的钢轨开式砂带打磨接触理论建模与仿真分析 HOU Guangyou . Contact theoretical modeling and simulation analysis of rail grinding using open belt based on grinding plate [D]. Beijing : Beijing Jiaotong University , 2020 : 1 - 40 . Contact theoretical modeling and simulation analysis of rail grinding using open belt based on grinding plate 尹波润 , 文永蓬 , 尚慧琳 . 基于元胞自动机方法的地铁车轮磨损动态建模与仿真 [J]. 机械工程学报 , 2019 , 55 ( 2 ): 135 - 146 . 基于元胞自动机方法的地铁车轮磨损动态建模与仿真 YIN Borun , WEN Yongpeng , SHANG Huilin . Dynamic modeling and simulation of metro wheel wear based on cellular automata method [J]. Journal of Mechanical Engineering , 2019 , 55 ( 2 ): 135 - 146 . Dynamic modeling and simulation of metro wheel wear based on cellular automata method XU Haiwe , YANG Wangyue , SUN Zuqing . Microstructure control of fine grained dual phase low carbon steel based on the dynamic transformation [J]. Acta Metallurgica Sinica , 2006 , 42 ( 10 ): 1101 - 1108 . Microstructure control of fine grained dual phase low carbon steel based on the dynamic transformation YEGANEH M , ESKANDARI M , ALAVI-ZAREE S R . A comparison between corrosion behaviors of fine-grained and coarse-grained structures of high-Mn steel in NaCl solution [J]. Journal of Materials Engineering and Performance , 2017 , 26 ( 6 ): 2484 - 2490 . A comparison between corrosion behaviors of fine-grained and coarse-grained structures of high-Mn steel in NaCl solution 周清跃 , 张银花 , 杨来顺 , 等 . 轨道的材质性能及相关工艺 [M]. 北京 : 中国铁道出版社 , 2005 : 9 - 29 . 轨道的材质性能及相关工艺 ZHOU Qingyue , ZHANG Yinhua , YANG Laishun , et al . Material properties and related technology of track [M]. Beijing : China Railway Publishing House , 2005 : 9 - 29 . Material properties and related technology of track 李闯 , 周韶博 , 王猛 , 等 . 中日韩高速铁路钢轨标准和维修规则对比分析 [J]. 铁道技术监督 , 2020 , 48 ( 12 ): 1 - 5 . 中日韩高速铁路钢轨标准和维修规则对比分析 LI Chuang , ZHOU Shaobo , WANG Meng , et al . Comparative analysis of rail standards and maintenance rules for high-speed railways in China, Japan and South Korea [J]. Railway Quality Control , 2020 , 48 ( 12 ): 1 - 5 . Comparative analysis of rail standards and maintenance rules for high-speed railways in China, Japan and South Korea ZHU Liqiang , DUAN XiangYu , YU ZuJun . On the identification of elastic moduli of in-service rail by ultrasonic guided waves [J]. Sensors , 2020 , 20 ( 6 ): 1 - 24 . On the identification of elastic moduli of in-service rail by ultrasonic guided waves MA Yuewei , MARKINE V L , MASHAL A A . et al . Effect of wheel-rail interface parameters on contact stability in explicit finite element analysis [J]. Proceedings of the Institution of Mechanical Part F-Journal of Rail and Rapid Transit , 2018 , 232 ( 6 ): 1879 - 1894 . Effect of wheel-rail interface parameters on contact stability in explicit finite element analysis SHARMA S , SANGAL S , MONDAL K . Wear behaviour of bainitic rail and wheel steels [J]. Materials Science and Technology , 2016 , 32 ( 4 ): 266 - 274 . Wear behaviour of bainitic rail and wheel steels GB/T 228.3—2019 . 金属材料-拉伸试验. 第3部分: 低温试验方法 [S]. GB/T 228.3—2019 . Metallic materials-tensile test. Part 3: low temperature test method [S]. 王磊 . 材料的力学性能 [M]. 3版 . 沈阳 : 东北大学出版社 , 2014 : 1 - 40 . 材料的力学性能 WANG Lei . Mechanical properties of materials [M]. 3rd ed . Shenyang : Northeast University Press , 2014 : 1 - 40 . Mechanical properties of materials CHOI W , SONG MinJi , LIM N H , et al . Microstructure and mechanical properties of continuous welded 50 N rail [J]. Korean Journal of Metals and Materials , 2019 , 57 ( 12 ): 755 - 763 . Microstructure and mechanical properties of continuous welded 50 N rail 翟婉明 . 车辆-轨道耦合动力学 [M]. 4版 . 北京 : 科学出版社 , 2015 : 210 - 254 . 车辆-轨道耦合动力学 ZHAI Wanming . Vehicle-track coupling dynamics [M]. 4th ed . Beijing : Science Press , 2015 : 210 - 254 . Vehicle-track coupling dynamics ZHAI Wanming , GAO Jianmin , LIU Pengfei , 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 . Reducing rail side wear on heavy-haul railway curves based on wheel-rail dynamic interaction 王少锋 . 重载线路钢轨裂纹研究 [M]. 成都 : 西南交通大学出版社 , 2018 : 87 - 127 . 重载线路钢轨裂纹研究 WANG Shaofeng . Study on rail crack of heavy-haul line [M]. Chengdu : Southwest Jiaotong University Press , 2018 : 87 - 127 . Study on rail crack of heavy-haul line 10.11817/j.issn.1672-7207.2023.05.034 U213.7 A 1672-7207（2023）05-2030-09 陈宪麦, 陈楠, 魏子龙, 等. 基于离散元法的钢轨力学性能定性分析[J]. 中南大学学报(自然科学版), 2023, 54(5): 2030-2038. 陈宪麦, 陈楠, 魏子龙, 等. 基于离散元法的钢轨力学性能定性分析[J]. 中南大学学报(自然科学版), 2023, 54(5): 2030-2038. CHEN Xianmai, CHEN Nan, WEI Zilong, et al. Qualitative analysis of rail mechanical properties based on discrete element method[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 2030-2038. CHEN Xianmai, CHEN Nan, WEI Zilong, et al. Qualitative analysis of rail mechanical properties based on discrete element method[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 2030-2038. 国家自然科学基金资助项目(51478482)；中国国家铁路集团有限公司科技研究开发计划重点课题(N2019G012)；中国国家铁路集团有限公司科技研究开发计划项目(K2019G043) Project(51478482) supported by the National Natural Science Foundation of China; Project(N2019G012) supported by the Key Program of Science and Technology Research and Development Plan of the China National Railway Group Co. Ltd.; Project(K2019G043) supported by the Science and Technology Research and Development Plan of the China National Railway Group Co. Ltd. 2022-07-10 2022-09-16 2023-05-26 2023-06-20 38909472 中南大学学报(自然科学版) 5 54