DOI： 10.11817/j.issn.1672-7207.2020.06.026

地铁车轮凹陷磨耗对踏面接触应力的影响

周志军,李伟,陈斌,温泽峰

（西南交通大学 牵引动力国家重点实验室，四川 成都，610031）

摘要:针对地铁列车运营后出现的车轮踏面滚动接触疲劳现象，利用有限元软件ABAQUS建立考虑轮轨真实几何关系的三维弹塑性轮轨接触数值仿真模型，结合现场调查车轮磨耗结果和轮轨接触几何关系分析，计算分析车轮不同凹陷磨耗状态对踏面材料应力应变状态的影响。研究结果表明：车轮踏面主要凹陷磨耗区域为车轮踏面横向位置-30~50 mm处，轮轨接触几何关系呈强非线性特性，其轮轨接触点位置集中在车轮踏面横向位置20～32 mm或-32 ~ -20 mm。车轮不同凹陷磨耗状态下的轮轨接触状态差异显著，在磨耗突变区(-30 ~ -20 mm)轮轨接触斑呈狭长椭圆分布，导致相同载荷下轮轨接触应力较大。当轮对向外轨横移时，车轮凹陷磨耗接触区域材料易进入屈服状态，此时车轮踏面沿接触斑深度方向0~3.6 mm范围内Von-Mises等效应力最大，踏面表层材料等效塑性应变最大。车轮踏面出现凹陷磨耗后，车辆频繁地通过小半径曲线时易在磨耗突变区造成较高的等效应力和等效塑性应变，从而导致轮缘根部表面材料产生剥离损伤。

关键词:地铁；车轮凹陷磨耗；剥离损伤；弹塑性；Von-Mises等效应力；塑性应变

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

文章编号:1672-7207（2020）06-1712-09

Influence of metro wheel hollow-worn on contact stress of wheel tread

ZHOU　Zhijun, LI　Wei, CHEN　Bin, WEN　Zefeng

(State Key Laboratory of Traction Power, Southwest Jiaotong University, Chengdu 610031, China)

Abstract: In view of the phenomenon that wheel tread fatigue damage occurs after metro train operation, a three-dimensional simulation model of wheel-rail elastic-plastic contact considering the true geometry of wheel-rail was established by using the finite element software ABAQUS. Combined with the test results of wheel wear and the analysis of wheel-rail contact geometry, the effects of wheel hollow-worn on the stress and strain states of the tread material were investigated. The results show that the main wear area of hollow-worn wheel tread is from -30 to 50 mm in the lateral position of the wheel thread. The position of wheel-rail contact point is concentrated in from 20 to 32 mm or from -32 to -20 mm in the lateral position of the wheel thread, and the contact geometry is strongly non-linear. The shape of wheel-rail contact area varies obviously under different wear conditions. In abrupt wear region(from -32 to -20 mm), the shape of wheel-rail contact patch shows a narrow ellipse, which leads to greater wheel-rail contact stress under the same load. When the wheelset moves to the outer rail, the material in the contact area of wheel hollow-worn is easy to enter the yield state. At this time, the Von-Mises equivalent stress of wheel reaches the maximum in the range of 0-3.6 mm below the contact patch, and the equivalent plastic strain of tread surface material reaches the maximum. When the vehicle with hollow-worn wheels passes through the small radius curve frequently, high equivalent stress and equivalent plastic strain in the abrupt wear region occur, resulting in shelling damage of the wheel flange root.

Key words: metro; wheel hollow-worn; shelling damage; elastic-plastic; Von-Mises equivalent stress; plastic strain

地铁车辆在运行过程中受线路条件、轮轨参数匹配、运行速度和频繁牵引制动等因素的影响，车轮踏面不可避免地发生磨耗和疲劳损伤等现象。车轮凹陷磨耗是典型的踏面磨耗现象，表现为轮缘磨耗轻微，名义滚动圆附近磨耗严重，两侧磨耗量逐次递减。车轮出现凹陷磨耗后，轮轨接触区域产生强非线性接触特性，轮轨相互作用关系恶化。车轮踏面在长期反复滚动接触载荷作用下，会产生局部永久性累积损伤，导致接触表面产生剥离等疲劳损伤现象，影响列车安全运行。国内外学者针对车轮凹陷磨耗和踏面剥离损伤现象，结合现场试验和数值仿真法展开了研究。SAWLEY等^[1]调查发现北美货车6 700个车轮踏面廓形中21%的车轮凹陷磨耗量大于2 mm，模拟结果表明，当凹陷磨耗量大于2 mm时，轮轨滚动阻力增加，轮轨接触应力明显增加。黄照伟等^[2]的研究表明当高速列车凹陷磨耗车轮与标准钢轨匹配时，在多个平衡位置间跳跃造成“轮缘到假轮缘”的冲击振动，影响车辆运行性能。孙宇等^[3]基于传统的Kik-Piotrowski方法，考虑轮对摇头和多点接触问题，仿真分析了具有磨耗车轮的轮轨动力相互作用关系，发现车轮凹陷磨耗会增加轮轨高频垂向振动，使轮轨横向力的振动幅值显著增加，影响行车平稳性和安全性。陶功权等^[4-5]基于车辆系统动力学仿真模型，建立了考虑安定图和损伤函数的车轮滚动接触疲劳分析模型，研究了车辆通过3种不同半径曲线时轮轨磨耗状态对车轮滚动接触疲劳损伤的影响。陆文教^[6]研究了轮轨磨耗对车辆动力学性能及轮轨接触损伤特性的影响，发现车轮磨耗后(踏面磨耗均匀，轮缘磨耗明显)，表面疲劳指数比标准轮轨时的大，出现滚动接触疲劳的可能性增大。陶功权等^[7]对比分析了不同磨耗车轮与标准CN60钢轨匹配时，利用三维弹性体非Hertz 滚动接触理论及其数值程序CONTACT以及三维轮轨接触有限元模型计算轮轨接触应力的算法，发现CONTACT程序计算得到的轮轨接触斑面积小于有限元计算结果，但CONTACT计算得到的轮轨间最大接触压力和最大等效应力大于有限元结果。针对车轮磨耗后容易出现的多点接触问题，2种算法误差较大，采用有限元计算结果更可靠。由于有限元方法计算轮轨接触应力时不受弹性无限半空间和材料弹性假设限制，赵腾等^[8]建立了不同轮缘磨耗程度的轮轨弹塑性接触模型，发现车轮轮缘厚度由32 mm磨耗至27 mm过程中，与磨耗钢轨匹配关系较好，接触斑面积较大，轮缘部分的接触等效应力较小。结合现场试验和有限元方法，FROHLING等^[9]研究了车轮踏面凹陷磨耗导致的滚动接触疲劳损伤危害问题，讨论了车辆曲线通过性能和凹陷磨耗限值。LIU等^[10]基于轮对纵向振动研究车轮踏面剥落机理，发现轮对纵向振动改变了轮轨接触特性、力学特性及接触面上主应力方向，导致轮轨接触面疲劳，最终导致车轮踏面剥落。结合全尺寸制动装置试验和基于有限元法，ESMAEILI等^[11]考虑车轮踏面滚动接触疲劳和反复制动同时受到的机械载荷和热载荷问题，对踏面制动进行数值模拟，研究了制动温度和滚动接触疲劳对踏面裂纹的影响。我国某城市地铁线路车辆运营后，车轮普遍出现了凹陷磨耗和踏面接触疲劳剥离现象。经现场调查发现，车轮踏面剥离位置位于轮缘根部附近，呈不规则的网状或鱼鳞状裂纹，沿裂纹处有层状剥离掉块特征。在车辆实际运营过程中，由于车轮凹陷磨耗导致其型面发生显著变化，恶化了轮轨接触相互作用关系。因此，本文作者开展车轮磨耗状态下对轮轨接触状态影响的研究，基于地铁列车运营后出现的车轮踏面疲劳损伤现象，利用有限元软件ABAQUS建立考虑轮轨真实几何关系的三维轮轨弹塑性接触数值仿真模型，结合轮轨接触几何关系分析，计算车轮不同凹陷磨耗状态下踏面材料应力、应变状态及其对踏面剥离损伤的影响。

1 轮轨接触几何关系分析

轮轨接触几何是研究轮轨关系的基本参数，轮轨型面的匹配程度直接影响列车运行的稳定性、平稳性、曲线通过性能及轮轨接触力和磨耗等。

1.1　凹陷磨耗车轮廓形及磨耗分布

现场试验调查发现全列车轮中80%以上的车轮踏面磨耗量大于2 mm，50%以上的车轮踏面凹陷量(车轮踏面范围内，凸起与凹陷部分的差值)大于0.4 mm，但钢轨磨耗比较轻微。定义T1，T2和T3为3种不同踏面磨耗量轮对，其中T1为标准LM型车轮踏面。

图1所示为3种不同磨耗状态的车轮廓形和磨耗分布情况，车轮磨耗量统计见表1。由图1和表1可知：车轮轮缘磨耗轻微，踏面磨耗较严重。由图1可见：车轮踏面主要凹陷磨耗区域为车轮踏面横向位置-30 ~ 50 mm，且在-30 ~ -20 mm处存在较明显的磨耗突变。

图1　车轮廓形及踏面磨耗分布

Fig. 1　Wheel profile and tread wear distribution

表1　车轮磨耗量统计

Table 1　Wheel wear statistics

1.2　轮轨接触点位置

车轮出现凹陷磨耗后会严重影响轮轨接触点位置，选取不同踏面磨耗状态(3种工况)的车轮与60 kg·m^-1标准钢轨进行轮轨接触几何关系的数值计算(不考虑摇头角)。其中，车轮名义滚动圆半径为420 mm，轨距为1 435 mm，轨底坡为1/40，轮对横移量为-9~9 mm，每隔0.5 mm计算1个接触点。

图2所示为不同工况下轮轨接触点对分布图。图2中，正值表示轮对向左侧横移，负值表示轮对向右侧移动，0表示轮轨对中接触。图2(a)中，T1轮对发生横移时轮对接触点位置随横移量基本均匀分布，接触状态良好。图2(b)和(c)中，车轮凹陷磨耗后轮轨接触几何呈强非线性特性。T2左侧车轮横移-8~0 mm和右侧车轮横移0.5~8.0 mm时，轮对接触点位置分别集中在车轮踏面横向位置-30 ~ -20 mm或20~30 mm。T3左侧车轮横移-8.0 ~ -0.5 mm和右侧车轮横移0~8.0 mm时，轮对接触点位置位于车轮踏面横向位置-32 ~ -20 mm或20~32 mm。另外，轮轨对中接触时，在多个平衡位置处轮轨接触点出现了特别明显的跳跃，导致在跳跃处易产生两点接触，轮对运行会伴随有一定的冲击振动。

若考虑线路不平顺时轮轨接触点位置将会在不考虑不平顺时的轮轨接触点位置附近波动，因此凹陷磨耗车轮的轮轨接触点分布位置都较接近。横移量为0~8.0 mm或-8.0 ~ -0.5 mm时，主要接触区域位于车轮车轮踏面横向位置-32 ~ -20 mm或20~32 mm。当车辆频繁通过小半径曲线时，易造成该区域轮轨接触频次较多，由于踏面磨耗突变导致的轮轨接触应力较大，踏面材料易发生疲劳剥离损伤。这与现场调查发现的车轮踏面剥离位置基本吻合。

图2　不同工况下的轮轨接触点对图

Fig. 2　Wheel-rail contact point diagrams under different conditions

1.3　滚动圆半径差

列车在运行过程中受外界激扰和轮轨系统本身的影响容易产生蛇行运动，伴随出现轮对不同横移现象^[12]。带有凹陷磨耗车轮的车辆在实际线路运行时，滚动圆半径差在不同横移量时会发生明显的变化。图3所示为T1，T2和T3轮对滚动圆半径差与轮对横移量之间的关系。由图3可知：T1在横移量为-9~9 mm范围内滚动圆半径差变化较平缓；T2和T3在横移量为0~2 mm和-2~0 mm范围内滚动圆半径差增大，在横移量为2~8 mm和-8 ~ -2 mm的滚动圆半径差值变化相对缓慢。当凹陷磨耗车轮通过小半径曲线时，轮对横移量增大，车轮滚动圆半径差增大，轮轨纵向蠕滑力随之升高，车轮踏面磨耗和疲劳损伤增大^[12]。

图3　滚动圆半径差与横移量之间的关系

Fig. 3　Relationship between radius difference of rolling circle and transverse displacement

以上轮轨接触几何关系表明车辆通过小半径曲线时，轮轨在磨耗突变区容易集中接触，而现场观测到的车轮滚动接触疲劳裂纹就发生在此区域，在车轮磨耗突变区域，轮轨接触半径较小，轮轨接触应力较高，材料容易发生塑性变形，很难用解析法对其分析^[13]。虽然KULKARNI等^[14]建立了弹塑性滚动接触问题的分析方法，但方法过于简化，且其前提假设是满足赫兹接触条件，无法反映真实轮轨关系，需考虑非赫兹接触。下面基于轮轨接触几何分析结果，通过建立考虑轮轨真实几何和材料弹塑性本构关系的三维轮轨接触数值模型，进一步探明车轮凹陷磨耗对轮轨接触区域的踏面接触应力及踏面剥离的影响。

2 三维弹塑性轮轨接触有限元模型

图4　三维弹塑性轮轨接触有限元模型

Fig. 4　Finite element model of three-dimensional elastic-plastic wheel-rail contact

利用有限元软件ABAQUS建立三维弹塑性轮轨接触有限元模型，如图4所示。模型考虑了不同凹陷磨耗状态车轮、60 kg·m^-1标准钢轨和轨道板实体模型，具体参数见表2。

表2　三维弹塑性轮轨接触有限元模型结构参数

Table 2　Structure parameters for finite element model of three-dimensional elastic-plastic wheel-rail contact

张宝安等^[15]的研究表明当轮对横移量超过5 mm时，轮轨的结构柔性会导致轮轨接触参数发生明显变化，因此在该模型中将轮对考虑为柔性是必要的。模型中未将轨枕和扣件实体考虑在内，两侧离散支撑的扣件系统由间距为50 mm的20组弹簧阻尼单元模拟，且考虑了垂横纵方向刚度阻尼特性。将车辆簧上质量M等效为一质量点，并通过一系悬挂(弹簧阻尼单元)与车轴相连。轮轨接触类型采用“面－面”接触，法向接触关系定义为“硬”接触，切向接触关系采用罚函数的摩擦公式定义，轮轨摩擦因数为0.5。为了分辨接触斑内接触压力的分布状态，将轮轨接触区域网格细化，其他区域进行过渡稀疏网格划分，最小单元尺寸为1 mm。车轴对称面横向(x)对称约束，端面横向位移约束；钢轨端面纵向(y)对称约束；轨道板对称面横向对称约束，端面纵向对称约束，侧面横向位移约束。实体模型采用C3D8R八节点六面体单元离散，共673 076个节点和585 688个单元。因考虑到轨道、车轮和轨道板的对称性，本模型只取车轮、钢轨和轨道板的一半建模，以提高求解效率^[16]。模型中车轮沿横向平移施加横移量，沿垂向(z)调整接触主从表面的初始位置，且在轮轨接触时设置误差限度使得车轮与钢轨法向接触，确保轮轨接触求解完全收敛。

考虑车轮磨耗突变接触区域的局部塑性变形，轮轨材料采用双线性随动强化弹塑性材料模型，材料本构关系如图5所示。钢轨材料为U75V，车轮材料为CL60，屈服条件为Von-Mises屈服准则。线性强化塑性材料模型的应力-应变方程^[17]为

(1)

式中：σ为应力；ε为应变；σ_s为屈服极限；E_e为弹性模量；ε_s为屈服点处总应变；E_p为应变强化模量；钢轨和车轮E_p为17.254 GPa。

图5　轮轨材料本构关系

Fig. 5　Constitutive relationship of wheel-rail material

3 结果分析

利用建立的轮轨接触有限元模型，计算分析轮对不同横移量下(-9~9 mm)轮轨接触参数(包括轮轨接触斑面积、最大接触压力、Von-Mises等效应力、等效塑性应变)及其对应的变化趋势。

车轮发生凹陷磨耗后，轮轨主要接触点位于车轮踏面横向位置20~32 mm或-32 ~ -20 mm。图6所示为典型的轮对横移量为-6 mm时轮轨接触斑情况。图中，a和b分别为接触斑长轴和短轴，A为接触斑面积，x为接触斑中心与钢轨中心线的相对距离。由图6可知：3种工况轮轨接触斑形状有显著变化(轮轨接触斑面积分别为107.3，57.3和57.4 mm²)，T1轮轨接触斑近似圆形，T2和T3轮轨接触斑为沿纵向狭长椭圆形。值得注意的是，凹陷磨耗车轮相对标准车轮的轮轨接触斑中心位置向钢轨内侧偏移约6.7 mm，接触位置在钢轨轨距角或轮缘根部附近。凹陷磨耗车轮狭长的轮轨接触斑面积使得相同载荷下轮轨接触应力较高，导致在该接触位置中心处容易出现滚动接触疲劳。对于T1车轮(标准车轮踏面)，轮轨最大接触应力为960 MPa；对于凹陷磨耗T2和T3车轮，轮轨最大接触应力分别为1 465和1 490 MPa，与T1车轮相比，增加了约53%。

图6　轮轨接触斑面积

Fig. 6　Area of wheel-rail contact patches

图7所示为轮对不同横移量时轮轨接触Von-Mises等效应力变化趋势。由图7可知：当轮对横移量为正时，T2和T3的Von-Mises等效应力变化趋势基本相同；当轮对对中接触时，T1的Von-Mises等效应力最大(605 MPa)，T3的最小(460 MPa)，这是因为凹陷磨耗车轮在横移量为0 mm时，轮轨呈曲面接触，凹陷磨耗量较大的车轮接触斑面积较大；当轮对横移量为负时，3种工况下轮轨接触Von-Mises等效应力差异明显，在凹陷磨耗车轮的轮轨频繁接触区域内(轮对横移量为-0.5~-8.0 mm)，T2和T3轮轨Von-Mises等效应力剧增。

图7　Von-Mises等效应力随轮对横移量变化趋势

Fig. 7　Trend of maximum Von-Mises equivalent stress of wheel-rail with wheelset transverse shift

图8所示为车轮横移量为-6 mm时轮轨接触区域Von-Mises等效应力云图。由图8可知：3种工况下轮轨接触区域均出现了较明显的等效应力集中区域，凹陷磨耗车轮相比标准车轮的等效应力集中点向钢轨内侧偏移倾向，其位置靠近轮缘和钢轨轨距角。为了表征Von-Mises等效应力影响范围，分析Von-Mises等效应力沿车轮踏面接触斑深度方向的变化趋势，如图9所示。由图9可知：凹陷磨耗车轮在接触斑表面以下0~3.6 mm深度范围内Von-Mises等效应力最大，而标准车轮Von-Mises等效应力较小。另外，与T1轮对相比，T2和T3车轮对Von-Mises应力出现的位置更接近车轮表面；随着深度的增大，3种工况下Von-Mises等效应力逐渐减小。

图8　轮轨Von-Mises等效应力云图

Fig. 8　Von-Mises equivalent stress distribution of wheel-rail

图9　Von-Mises等效应力沿接触斑深度方向的变化

Fig. 9　Variation of Von-Mises equivalent stress along depth of contact patch

图10所示为轮对不同横移量时轮轨等效塑性应变变化情况。由图10可知：当轮对横移量为负时，凹陷磨耗车轮相比标准车轮的轮轨等效塑性应变急剧增大，且在横移量为-9~-6 mm时其值达到最大。这主要是因为横移量为负时轮轨接触点位置发生较大偏移，轮轨接触斑面积变小和应力集中导致的。图11所示为等效塑性应变最大值处的轮轨等效塑性应变云图。由图11可知：T1轮轨等效塑性应变较小，塑性变形区主要集中在车轮次表层。T2和T3轮轨等效塑性应变较大，T3等效塑性应变约为T2的4.35倍，为T1的151倍，变形区集中在车轮表面且沿车轮轮缘圆周狭长分布。凹陷磨耗车轮较大的塑性应变容易引起接触表面材料的塑性流动和残余应力的产生。车轮踏面出现凹陷磨耗后，车辆频繁地通过小半径曲线，在长期反复滚动接触载荷作用下，易在磨耗突变区造成较高的等效应力和等效塑性应变，产生局部永久性累积损伤，导致轮缘根部表面材料剥离损伤，这与现场观察到的踏面剥离基本吻合。

图10　轮轨等效塑性应变随轮对横移量的变化

Fig. 10　Change of equivalent plastic strain of wheel-rail with wheelset transverse shift

图11　轮轨等效塑性应变云图

Fig. 11　Equivalent plastic strain distribution of wheel-rail

4 结论

1) 地铁列车车轮踏面主要凹陷磨耗区域为车轮踏面横向位置-30~50 mm，全列车轮中80%以上的车轮踏面磨耗量大于2 mm，钢轨磨耗轻微。轮轨接触几何关系呈强非线性特性，其接触点位置主要集中在车轮踏面横向位置20~32 mm或-32 ~ -20 mm，与现场调查发现的车轮踏面剥离位置吻合。且当轮轨对中接触时，在多个平衡位置处轮轨接触点位置出现了明显的跳跃。

2) 车轮不同凹陷磨耗状态下的轮轨接触状态差异显著。尤其在磨耗突变区(-30 ~ -20 mm)轮轨集中接触，接触斑面积呈狭长椭圆分布，导致相同载荷下轮轨接触应力较大。

3) 当轮对向外轨横移时，车轮凹陷磨耗接触区域材料易进入屈服状态，此时车轮踏面沿接触斑深度方向0~3.6 mm范围内Von-Mises等效应力最大，踏面表层材料等效塑性应变最大。当车轮踏面出现凹陷磨耗后，车辆频繁地通过小半径曲线时易在磨耗突变区造成较高的等效应力和等效塑性应变，从而导致轮缘根部表面材料产生剥离损伤。

参考文献：

[1] SAWLEY K, URBAN C, WALKER R. The effect of hollow-worn wheels on vehicle stability in straight track[J]. Wear, 2005, 258(7/8): 1100-1108.

[2] 黄照伟, 崔大宾, 杜星, 等. 车轮偏磨对高速列车直线运行性能的影响[J]. 铁道学报, 2013, 35(2): 14-20.

HUANG Zhaowei, CUI Dabin, DU Xing, et al. Influence of deviated wear of wheel on performance of high-speed train running on straight tracks[J]. Journal of the China Railway Society, 2013, 35(2): 14-20.

[3] 孙宇, 朱胜阳, 翟婉明. 车轮踏面凹形磨耗对轮轨相互作用的影响研究[J]. 机械工程学报,2018, 54(4): 109-116.

SUN Yu, ZHU Shengyang, ZHAI Wanming. Influence of tread hollow-worn wheel on wheel/rail interaction[J]. Journal of Mechanical Engineering, 2018, 54(4): 109-116.

[4] 陶功权, 王衡禹, 赵鑫, 等. 基于轮轨关系的车轮踏面损伤机理研究[J]. 机械工程学报, 2013, 49(18): 23-29.

TAO Gongquan, WANG Hengyu, ZHAO Xin, et al. Research on wheel tread damage mechanism based on interaction of wheel and rail[J]. Journal of Mechanical Engineering, 2013, 49(18): 23-29.

[5] 陶功权. KKD客车车轮踏面剥离成因初探[D]. 成都: 西南交通大学牵引力国家重点实验室, 2013: 9-10.

TAO Gongquan. Preliminary study on the cause of wheel tread shelling of KKD passenger car[D]. Chengdu: Southwest Jiaotong University. State Key Laboratory of Traction Power. 2013: 9-10.

[6] 陆文教, 陶功权, 王鹏, 等. 地铁车轮磨耗对轮轨接触特性及动力学性能的影响[J]. 工程力学, 2017, 34(8): 222-231.

LU Wenjiao, TAO Gongquan, WANG Peng, et al. Influence of wheel wear on wheel-rail contact behavior and dynamic performance of metro vehicle[J]. Engineering Mechanics, 2017, 34(8): 222-231.

[7] 陶功权, 李霞, 温泽峰, 等. 两种轮轨接触应力算法对比分析[J]. 工程力学, 2013, 30(8): 229-235.

TAO Gongquan, LI Xia, WEN Zefeng, et al. Comparative analysis of two algorithms for wheel-rail contact stress[J]. Engineering Mechanics, 2013, 30(8): 229-235.

[8] 赵腾, 张军, 孙传喜. 机车轮缘磨耗对轮轨接触状况的影响[J]. 工程力学, 2012, 29(10): 308-312.

ZHAO Teng, ZHANG Jun, SUN Chuanxi. Effects of flange wear on the wheel/rail contact status of locomotive[J]. Engineering Mechanics, 2012, 29(10): 308-312.

[9] FROHLING R, EKBERG A, KABO E. The detrimental effects of hollow wear-field experiences and numerical simulations[J]. Wear, 2008, 265(9/10): 1283-1291.

[10] LIU Wei, MA Weihua, LUO Shihui, et al. Research into the problem of wheel tread spalling caused by wheelset longitudinal vibration[J]. Vehicle System Dynamics, 2015, 53(4): 546-567.

[11] ESMAEILI A, WALIA M S, HANDA K, et al. A methodology to predict thermomechanical cracking of railway wheel treads: From experiments to numerical predictions[J]. International Journal of Fatigue, 2017, 105: 71-85.

[12] 王晨, 马卫华, 罗世辉, 等. 轮径差对地铁车辆轮对磨耗的研究[J]. 铁道机车车辆, 2013, 33(S1): 87-91.

WANG Chen, MA Weihua, LUO Shihui, et al. Study on wheel wear of the subway vehicle based on radius difference[J]. Railway Locomotive & Car, 2013, 33(S1): 87-91.

[13] 张立民. 轮轨接触应力与钢轨波磨分析[J]. 西南交通大学学报, 2003, 38(1): 34-37.

ZHANG Limin. Analysis of wheel/rail contact stresses and rail corrugation[J]. Journal of Southwest Jiaotong University, 2003, 38(1): 34-37.

[14] KULKARNI S M, HAHN G T, RUBIN C A, et al. Elasto-plastic finite element analysis of repeated three-dimensional, elliptical rolling contact with rail wheel properties[J]. Journal of Tribology, 1991, 113(3): 434-441.

[15] 张宝安, 陆正刚, 王恒亮. 柔性轮对的轮轨接触参数分析[J]. 同济大学学报(自然科学版), 2013, 41(4): 577-582.

ZHANG Bao'an, LU Zhenggang, WANG Hengliang. Analysis on wheel-rail contact parameters of elastic wheelset[J]. Journal of Tongji University(Natural Science), 2013, 41(4): 577-582.

[16] 朱强强. 地铁先锋扣件轨道钢轨波磨形成机理初探[D]. 成都: 西南交通大学牵引力国家重点实验室, 2018: 38-39.

ZHU Qiangqiang. A preliminary investigation into the cause of rail corrugation on subway track with vanguard fasteners[D]. Chengdu: Southwest Jiaotong University. State Key Laboratory of Traction Power. 2018: 38-39.

[17] 侯传伦, 翟婉明, 邓锐. 曲线磨耗状态下轮轨弹塑性接触有限元分析[J]. 中国铁道科学, 2009, 30(5): 28-33.

HOU Chuanlun, ZHAI Wanming, DENG Rui. Finite element analysis of the elastic-plastic contact of the worn wheels and rails on curve[J]. China Railway Science, 2009, 30(5): 28-33.

（编辑 伍锦花）

收稿日期： 2019 -10 -29; 修回日期： 2019 -12 -29

基金项目(Foundation item)：国家重点研发计划(2016YFB1200501-005)；四川省科技计划项目(2019YFH0053)；中央高校基本科研业务费专项资金资助项目(A0920502051904-43) (Project(2016YFB1200501-005) supported by the National Key Research and Development Program of China; Project(2019YFH0053) supported by the Science and Technology Plan of Sichuan Province; Project(A0920502051904-43) supported by the Fundamental Research Funds for the Central Universities)

通信作者：温泽峰，博士，研究员，从事轮轨关系及减振降噪研究；E-mail: zefengwen@126.com