基于摩擦自激理论的单侧钢轨波磨机理分析

肖 宏 ，陈 鑫 ，赵 越

（北京交通大学土木建筑工程学院，北京 100044）

钢轨波浪形磨耗是指钢轨表面出现的周期性波浪形不平顺现象，又称作波磨[1].从已有文献来看，钢轨波磨在世界各国的铁路及城市轨道交通都普遍存在.中国在重载铁路上波磨尤其严重，据统计，约40%以上的重载铁路曲线地段均存在钢轨波磨[2]，特别是在小半径曲线（半径R≤ 800 m）地段，波磨现象出现早，发展快.由波磨引起的各种轨道结构病害也急剧增加，这些不仅增加了重载铁路养护维修工作量，也直接威胁运营安全.

钢轨波磨产生机理复杂，影响因素众多，国内外学者据此进行了大量的研究，至今未形成共识：英国学者Xie等[3]发现如果选择赫兹接触，各种磨损和法向力之间的相移近似相同，由此导致初始不平顺继续发展；国内Jin等[4-6]对初始不平顺导致钢轨波磨产生的机理进行了分析；张厚贵[7]针对北京地铁异常波磨现象进行了大量的研究，认为剪切型减振扣件导致的轨道结构阻尼过低，钢轨极易发生不稳定振动，从而导致钢轨波磨的形成；Chen等[8-9]提出了钢轨波磨的发生源于饱和蠕滑力诱发的轮轨系统摩擦自激振动；肖祥龙等[9-12]提到小半径曲线地段钢轨波磨首先出现在一侧，但是关于波磨首先出现在哪一侧尚不明确.

由上述国内外已有研究可以看出，轮轨间饱和蠕滑力引发的不稳定摩擦自激振动是导致轮轨间产生波磨的主要诱因之一.同时，钢轨波磨的不对称性普遍存在，但这种现象的产生原因以及演化机理尚未被探讨.

基于此，本文在现场钢轨波磨调研基础上，采用饱和蠕滑力诱发轮轨系统摩擦自激振动的理论，建立精细化的轮轨三维接触有限元模型；基于等效阻尼比进行复特征值计算，分析不同扣件刚度、摩擦系数和超高等因素对轮轨系统不稳定摩擦自激振动的影响，研究单侧钢轨波磨的产生原因；进一步结合瞬态动力学方法模拟钢轨表面波磨几何特征，分析单侧钢轨波磨的传递及演化过程，揭示不对称钢轨波磨产生的原因，为养护维修提供指导.

1 现场测试

1.1 钢轨波磨测试

神朔重载铁路是我国西煤东运大通道神黄铁路的一部分，主要开行C80、C70、C64货车，平均速度50 km/h，另外还开行一辆通勤客车，平均速度70 km/h.为提高旅客乘坐的舒适性，在曲线地段的超高均按过超高设置.采用钢轨波磨测量仪（CAT）在现场测试钢轨表面的粗糙度.

图1为用CAT钢轨波磨仪测得的钢轨表面垂向不平顺局部放大图.测试的波磨完整图形见补充材料图S1，对应的频谱图见补充材料图S2.从图1可以看出：内轨表面波长约为200 mm，波深约为0.5 mm，属于中波波磨不平顺，计算可知其对应的振动频率为69 Hz；外轨也检测出波长为200 mm的不平顺，但是波深较小.此外，通过现场调研及跟踪测试发现，神朔重载铁路小半径曲线地段钢轨波磨先在内侧钢轨出现，之后随着时间的推移，另一侧钢轨表面也会产生同波长的波磨不平顺，整个过程中内侧钢轨波磨始终比外轨大.

图1 钢轨波磨测量仪所测表面垂向不平顺Fig.1 Vertical surface irregularity measured by CAT

1.2 轮轨作用力测试

为验证后续模型，在现场对轮轨的垂向作用力进行了测试，测试结果如图2所示.

图2 轮轨垂向力时程曲线Fig.2 Time-history curves of vertical wheel-rail force

从图2可以看出：在峰值及非峰值处，外轨受力均略小于内轨，内、外轨轮轨垂向力分别为132.11 kN和126.72 kN，内轨受力较大，这主要是由于曲线地段过超高，以及内轨不平顺严重造成的.该试验结果将用于后续模型验证.

2 钢轨波磨仿真分析模型

2.1 轮轨摩擦自激振动及瞬态动力作用形成波磨机理

诸多学者认为，轮轨间摩擦功的变化会导致钢轨表面产生相应的磨耗与伤损[13-19]，根据Brockley等[20-23]提出的轨面磨损量的计算公式以及摩擦自激振动理论，在轮轨系统摩擦自激振动的影响下，法向接触力会产生相同频率的振动，进而导致摩擦功也产生同频率的周期性波动，这几种频率均来自整个轮轨系统，进而导致单位时间内钢轨表面磨损量亦发生相应的波动，由此产生钢轨波磨，详细内容见附加材料1.

2.2 复特征值分析法

复特征值分析法是一种能够判定系统产生摩擦自激不稳定振动频率及其对应复模态主导振型的线性方法.由于涉及弹性体之间的耦合作用，考虑摩擦的系统模型表现得非常复杂.一般认为，高频率的工况不适用无质量弹簧方法.这里采用无弹性接触面只约束接触节点的方法[8]，详细推导过程见附加材料2.系统对应第i阶次的等效阻尼比 ξi可以判断系统自激振动发生的趋势，如式（1）.

式中： αi为特征值实部，i为阶次；ωi为特征值虚部，即当 ξi为负值时，系统发生自激振动的可能性也会增大，ωi绝对值越大，系统越容易发生不稳定自激振动.

2.3 轮轨三维接触有限元模型

重载铁路小半径曲线段轮轨接触力学模型，通过在轮对两端施加垂、横向悬挂力以模拟悬挂系统对轮对的作用，如图3所示.图中：FVL、FVR分别为曲线外轨、内轨侧转向架作用在车轴上的垂向力；FLL、FLR分别为曲线外轨、内轨侧转向架作用在车轴上的横向力；NL、NR分别为曲线外轨、内轨侧轮轨法向力；FL、FR分别为曲线外轨、内轨侧轮轨横向力； δL、 δR分别为曲线外轨、内轨侧轮轨接触角；KK1、KK2分别为扣件系统横向、垂向刚度；CK1、CK2分别为扣件系统横向、垂向阻尼；KZ1、KZ2分别为枕下横向、垂向刚度；CZ1、CZ2分别为枕下横向、垂向阻尼，其中轮轨横向蠕滑力、法向力以及接触角由多体动力学软件Simpack计算得到.车辆选用神华号电力机车，计算工况地段线路形式为曲线，半径R=400 m，超高 105 mm，车速 50 km/h.其中，车辆参数见附加材料表S1.材料瑞利阻尼系数β取0.001 5，依此建立机车的轮轨三维接触有限元模型，其材料参数与详细模型参数设置见附加材料表S2、S3所示.

图3 小半径曲线段轮轨接触力学模型Fig.3 Contact model of wheelset-track system with small radius curve section

由于车辆悬挂装置和减振系统的作用，在复模态分析时忽略其簧上质量的影响，着重考虑轮轨接触系统处在平衡状态时轮轨产生相对纵向滑移发生的弹性自激不稳定振动情况，在复模态分析模型中仅考虑轮对、钢轨和轨枕等主要结构.车轮踏面采用JM3型，钢轨廓形为CHN75型.钢轨长度取36 m，由61根Ⅲ型轨枕和扣件系统分散支撑.利用Python脚本语言批量建立wire连接以模拟钢轨和轨枕之间的扣件系统.模型单元类型为C3D8I，单元总数为 234 552 个，节点总数为 312 609 个，模型效果示意图见附加材料图S3.

由于本模型是根据重载铁路小半径曲线段现场实际线路参数以及列车车型、速度等运营参数建立的，将模型仿真计算结果与现场实测结果进行对比，如表1所示，现场轮轨力实测结果与仿真结果较为接近，所建模型可靠性较好，可用于后续的研究分析.

表1 计算结果对比Tab.1 Comparison of calculation results kN

3 钢轨波磨产生原因分析

扣件刚度、轮轨摩擦系数、超高与速度等对钢轨波磨有直接影响.

3.1 扣件刚度对轮轨系统不稳定性的影响

钢轨扣件刚度是保证轨道结构安全、评价扣件性能最重要的参数之一.应用复特征值分析法，就扣件刚度对轮轨系统摩擦自激振动的影响进行分析，从而进一步了解扣件垂、横向刚度对钢轨波磨产生的影响.根据现场实际情况，设置如表2所示的组合工况.

表2 扣件刚度计算工况Tab.2 Working conditions of fastener stiffness calculation MN/m

由摩擦自激振动理论可知：负等效阻尼比的绝对值越大，系统越容易发生摩擦自激振动，导致钢轨波磨进一步发展；钢轨波磨的发生频率与该不稳定摩擦自激振动频率相对应；最小等效阻尼比出现的频率所对应的不稳定振型为轮轨系统的主导振型，对系统的影响最大，最容易发生不稳定振动.

由于重载铁路钢轨波磨频率范围一般在20～120 Hz[2,7,9]，所以本文在复模态分析时，为了节约计算资源，只提取20～120 Hz内的模态计算结果.9种工况的等效负阻尼比值计算结果如图4所示（红色、绿色投影部分分别表示不同扣件垂向、横向刚度下等效阻尼比的大小），由图4可知：整个轮轨系统的负等效阻尼比绝对值会随着钢轨扣件垂、横向刚度的增大而减小，说明一定程度地增大扣件垂向刚度和横向刚度能够显著地抑制自激振动，从而抑制钢轨波磨的发展.

图4 不同垂、横向扣件刚度组合计算等效阻尼比Fig.4 Calculation of equivalent damping ratio by different combination of vertical and transverse fastener stiffness

选取具有代表性的刚度最大、最小和一般工况下负等效阻尼比分布及主导自激振动振型图，见附加材料图 S4～S6，放大系数均为 10 000.3 种工况的最小等效阻尼比分别出现在 65、70、74 Hz处，各工况的主导不稳定自激振动振型均出现在65～75 Hz范围内，与现场情况较为吻合.

同时，3种工况的最小负等效阻尼比值分别为-0.021 89、-0.021 01 和-0.019 82，变化不大.在第一种工况下，当扣件垂向刚度和横向刚度都较小时，不稳定振动振型最为复杂.这也反映了增大扣件垂向刚度和横向刚度，能够显著地抑制自激振动，从而抑制钢轨波磨的发展.

另外，对比振型可以发现，各工况最小负等效阻尼比对应的不稳定振动均发生在轮轨系统的内轮-内轨上，与该重载铁路小半径曲线现场波磨产生情况较为符合.

为了便于观察，附加材料图S4～S6中列出了其横断面视角，图5为工况1的主导不稳定自激振动振型对应的三维视角图.9种工况的主导不稳定振动振型非常类似，且均发生在曲线的内轮-内轨处，故此处不再重复罗列其三维视角图.

图5 三维视角位移（工况 1）Fig.5 Three-dimensional displacement(working condition 1)

查阅文献可知，北京地铁波磨病害现场也有类似反映[7]，使用弹性扣件的线路比使用普通扣件的地段更容易产生波磨病害，这和本文的分析结果十分吻合，也从实际角度进一步证明了饱和蠕滑力会导致轮轨系统摩擦自激振动，导致波磨的发生和发展.

3.2 摩擦系数对轮轨系统不稳定性的影响

作为轮轨接触摩擦副中最为核心的参数，轮轨表面摩擦系数会随着养护作业状况、干湿情况等因素发生变化.一般而言，对钢轨表面做涂油处理后，轮轨表面摩擦系数会降至0.2，对钢轨做撒沙处理时，轮轨表面摩擦系数会增加至0.6[9].为了更加清晰地揭示摩擦系数对轮轨系统的影响，在本文研究中将摩擦系数扩大至0.8，故此处分析摩擦系数为0.2、0.4、0.6、0.8 的 4 种工况.进行系统复特征值分析，最小负等效阻尼比计算结果分别为-0.017 65、-0.021 01、-0.074 76、-0.192 18，对应的频率分别为68、70、74、77 Hz，如图6 所示.

图6中蓝色原点的大小与负等效阻尼比的绝对值正相关，即蓝色圆点越大，对应的轮轨系统更易产生不稳定振动，从而更易导致钢轨波磨发展，如图7所示.

图6 最小负等效阻尼比对应的频率分布Fig.6 Frequency distribution corresponding to the minimum negative equivalent damping ratio

由图7可知：摩擦系数对系统稳定性影响较大，且呈现抛物线式的增长，尤其是系统在轮轨间摩擦系数达到0.4之后，负阻尼比的绝对值急剧增长.在保证正常轮轨接触和牵引制动需求的前提下，控制轮轨表面摩擦系数在0.4以下，能够较好地避免钢轨波磨的发展.从神朔重载铁路现场情况来看，长大上坡地段钢轨波磨的产生和发展，与过度撒沙增大轮轨间的摩擦系数有直接关系，导致了系统产生摩擦不稳定自激振动，诱发钢轨波磨病害.

图7 摩擦系数对轮轨系统振动稳定性的影响Fig.7 Influence of friction coefficient on vibration stability of wheel-rail system

3.3 超高和速度的匹配对轮轨系统不稳定性的影响

根据超高计算方法，当列车运行速度v= 50 km/h，曲线半径R= 400 m 时，计算得到超高值为 75 mm.表3列出了5种工况下前3阶负等效阻尼比分布及最小对应的自激振动振型.附加材料图S7～S11展示了5种工况的详情.

表3 不同超高下系统最小负等效阻尼比分布及主导振动振型Tab.3 Distribution of minimum negative equivalent damping ratio and dominant vibration mode for different superelevation

由表3可知：无论是过超高、欠超高还是平衡超高，均出现了不稳定振型；当轨道为过超高状态时，最小负等效阻尼比对应的不稳定振动振型仍出现在内轮-内轨处，当轨道为欠超高状态时，最小负等效阻尼比对应的不稳定振动振型出现在了外轮-外轨处；系统最小负等效阻尼比的绝对值对超高状态更为敏感，超高状态的改变能够很明显地使得系统的稳定性发生变化.这也很好地解释了部分文献中提到在我国其他一些重载铁路小半径曲线地段发生的钢轨波磨首先出现在曲线外轨上的情况.上述分析表明，超高状态直接影响波磨的发展，也决定了波磨首先产生在曲线内轨或者外轨上，揭示了钢轨单侧先产生的原因.

综上可知，适当地提高扣件垂横向刚度，控制轮轨摩擦系数在0.4以下以及平衡轨道结构超高能够有效地降低系统发生不稳定振动的可能性，从而抑制波磨的发展.

4 单侧钢轨波磨的传递与演化分析

通过瞬态动力学数值分析方法模拟单侧钢轨波磨不平顺，研究一侧钢轨波磨产生后对另一侧钢轨的影响.通过对比分析揭示钢轨波磨的传递与演化规律.

4.1 单侧钢轨波磨不平顺的模拟

根据图1中现场实测钢轨波磨特征进行一般化和抽象化处理，在钢轨内轨表面设置波深为0.5 mm、波长为200 mm的余弦形波磨不平顺，外轨表面则设置为普通平顺状态.

基于重载铁路现场运营的神华号电力机车的实际参数建立瞬态动力学模型，相关参数见附加材料表S4.由于模型体量较大，车轮与钢轨接触区域网格划分较密，利用Python语言对钢轨表面节点位置进行批量调整，调整后的模型如附加材料图S12所示.

4.2 波磨激发振动的传递

根据波磨在轨道上施加的位置可知，车辆会在0.50～0.65 s经过该波磨段.图8 为车辆以 50 km/h的速度通过仅内轨表面有波磨不平顺区段时两侧轮轨垂、横向力瞬态计算结果.图9为有波磨一侧钢轨和无波磨一侧钢轨应力云图.由图8、9可以看出：两侧钢轨应力状态相差较大；内轨和外轨的垂、横向轮轨作用力均在波磨的波峰、波谷位置出现了幅值波动，但没有波磨一侧的外轨轮轨作用力波动幅值明显小于有波磨一侧内轨轮轨力，其波动幅值分别为36.2 kN和80.1 kN；轮轨横向作用力因其幅值较小，相应的波动幅值也表现地较小.

图8 轮轨力随时间的变化情况Fig.8 Variation of wheel force with time

图9 钢轨 Mises应力云图Fig.9 Mises stress nephogram of rails

观察0.50～0.65 s的波磨施加段还可以发现：在进入波磨地段后，有波磨侧和无波磨侧轮轨作用力的波动变化逐渐表现出了一定的相位差，外侧波动要略晚于内侧的波动，且随着在波磨地段的滚动距离增大相位差也会逐渐变大，因此在通过波磨地段时，无波磨一侧的垂向力整体波动出现相位滞后；同时，横向力波动相较于垂向力波动也发生了一个明显相位差；有波磨一侧的扰动会导致系统发生不稳定振动，从而进一步传递至无波磨一侧；模型中虽然只施加了垂向不平顺，但是横向力也发生了明显的扰动，说明垂向振动的不稳定，也会逐渐传递至横向振动.

为了更清楚地展示系统随波磨状态的振动变化，将时域坐标换算成空间坐标，如图10所示.由图10可以看出：有波磨一侧的轮轨振动变化和无波磨一侧有明显的相位差，垂向力的振动变化和横向力也有明显的相位差，但变化频率均与波磨几何变化频率相同；轮轨垂横向力的一个变化周期都对应着波磨几何的相应周期.由此，更好地验证了在小半径曲线地段，钢轨波磨由于不合理的超高设置导致蠕滑力饱和诱发轮轨系统摩擦自激振动，曲线一侧钢轨首先演化出波浪形磨耗，随着轮轨间的振动响应，再逐渐传递至另一侧钢轨.

图10 系统随波磨状态的振动变化Fig.10 Vibration of the system with rail corrugation

5 结 论

1）首次从钢轨波磨的传递角度揭示波磨演化的全过程，建立的有限元模型能够很好模拟小半径曲线地段轮轨三维接触状态，对轮轨系统进行复特征值计算，分析其不稳定摩擦自激振动状态，从而揭示钢轨波磨先于一侧产生的原因是超高和实际运行速度的不匹配，该结果与实际情况较为吻合.

2）过超高设置会导致系统在内轮-内轨处发生不稳定摩擦自激振动，欠超高设置则会导致系统在外轮-外轨处发生不稳定自激振动.该计算结果与我国很多铁路线路钢轨波磨首先产生在单侧钢轨的现象相吻合.因此，在重载铁路小半径曲线地段设置外轨超高时，应尽量根据该区段车辆实际行驶速度来确定，同时调整摩擦状态，控制摩擦系数在0.4以下.此外，适当提高扣件刚度，也会减缓钢轨波磨演化和发展.

3）波磨一侧的扰动会导致轮轨系统发生不稳定振动，并逐渐传递至无波磨一侧，同时垂向不稳定振动也会逐渐传递至横向振动，并进一步放大轮轨动态响应.该发现进一步揭示了钢轨波磨由一侧钢轨传递至另一侧钢轨的过程，说明不合理的超高设置会导致饱和蠕滑力诱发摩擦自激振动进而产生单侧钢轨波磨，之后向另一侧钢轨传递演化.

备注：附加材料在中国知网本文的详情页中获取.

致谢：神华科技创新项目（SHGF-16-50）.