30t轴重重载铁路轮轨滑动接触引起的钢轨热相变分析

30t轴重重载铁路轮轨滑动接触引起的钢轨热相变分析

杨新文;顾少杰;周顺华;练松良

【摘要】轮轨接触热损伤是重载铁路轮轨结构的主要失效形式之一.当重载列车机车加速牵引时,轮轨接触界面间发生相对滑动,轮轨相对滑动会导致轮轨接触区出现较高的接触温升,这种接触温升积累到车轮和钢轨材料的相变时会发生相变,从而导致轮轨接触表面的擦伤、裂纹和剥离等.本文将钢轨看作是弹性半空间体.首先利用轮对切片投影法计算轮轨接触几何关系,然后利用接触弹性方法计算了轮轨接触斑和法向接触应力,最后将得到的轮轨法向接触应力施加到钢轨弹性半空间模型中,利用Matlab软件编制了轮轨接触热效应分析程序,计算分析了轮轨接触的热效应,并研究了钢轨材料马氏体层产生的厚度、位置及其分布.针对不同的轮对横移量、车轮移动速度、摩擦因数和相对滑动速度进行了钢轨的温度场研究.结果表明:轮轨滑动接触热效应剧烈,温升可以达到奥氏体化临界点,形成脆硬的马氏体.计算结果也表明,利用弹性半空间假设研究轮轨接触热效应问题是可行的.

【期刊名称】《铁道学报》

【年(卷),期】2016(038)007

【总页数】7页(P84-90)

【关键词】重载铁路;轮轨接触;热效应;温升;相变

【作者】杨新文;顾少杰;周顺华;练松良

【作者单位】同济大学道路与交通工程教育部重点实验室,上海201804;同济大学道路与交通工程教育部重点实验室,上海201804;同济大学道路与交通工程教育部

重点实验室,上海201804;同济大学道路与交通工程教育部重点实验室,上海201804

【正文语种】中文

【中图分类】U255

轮轨接触热损伤是重载铁路轮轨结构的主要失效形式之一。重载列车加速牵引时，车轮和钢轨接触面间容易发生滑动。轮轨间瞬时接触滑动会急剧增加接触区的温升，甚至可突破材料的临界相变点，从而产生局部奥氏体化，滑动摩擦过后，轮轨接触区的温度会快速下降淬火而形成马氏体组织。在摩擦热应力和荷载应力的共同作用下，脆硬的马氏体层容易产生裂纹，将成为引起钢轨裂纹与疲劳伤损的主要根源，见图1。

国内外学者对轮轨摩擦温升引起的热应力等问题做了大量的研究。在滑动体热传导的基本问题方面，Blok[2]和Jaeger[3]做出了首创性的工作。早期的学者主要利用拉普拉斯变换法以及格林函数法研究轮轨接触温升[4-6]，利用热源移动来模拟轮

轨相对滑动，他们的热源模型为瞬态静热源。由于大部分情况下轮轨接触不是椭圆接触，所以热源为非静热源。近些年借助计算机的发展，很多学者利用有限元法求解轮轨摩擦温升问题[7-13]，增加工程的可行度，但是有限元法同样存在着一些难以解决的问题，比如在模拟车辆运行速度、轮轨空间位置关系、轮轨蠕滑状态等方面还不能很好的实现。

铁路车辆行驶的过程中车轮是边滚边滑的状态，因此本文研究了不同车速、不同相对滑动速度下的钢轨热效应。钢轨表层的温升分布在接触表面下大约1.6 mm的

范围内[11,14]，热扩散影响的范围有限且轮轨接触是小变形问题，接触斑的大小

相对于轮轨曲率半径小很多，故本文将钢轨型面看作是弹性半空间体。首先利用轮对切片投影法计算轮轨接触几何关系[14]，然后利用虚拟渗透法计算轮轨接触斑，计算了非椭圆斑的接触应力，计算结果能更好的反应轮轨真实的接触状态。最后将真实的轮轨接触应力加在钢轨弹性半空间模型中，计算轮轨接触的热效应并分析了钢轨马氏体层产生的厚度、位置及其分布。

1 轮轨法向接触应力的计算

1.1 轮对切片投影法

钢轨是一个无限长柱状体，所有钢轨横截面所得的轮廓线都是一样的。因此钢轨顶面的形状与竖坐标z和横坐标y有关，而与纵坐标x无关[15]。而车轮是一个回转体，其踏面是一个回转面，它的形状与x、y、z的坐标位置都有关[16]。在轮对坐标系Gwx′y′z′中，作平面x′=0切割车轮，该平面与车轮周边的交线即为主轮廓线。再用平面分别切割车轮(见图2)，便可以得到n个轮廓线[15]。由于车轮是回转体，在轮廓线上的任意点的竖坐标可以根据横坐标相同的主轮廓线上的对应点的竖坐标求得

( 1 )

逐点求出轮廓线上各点的坐标，从而构成轮廓线。故由主轮廓线及所有x′轮廓线

构成的车轮踏面的轮廓面。再将轮廓面上的各点的坐标转换到轨道坐标系Oxyz上，再结合钢轨型面离散点进行扫描，得到轮轨法向间隙曲面[15]。根据轮轨的法向间隙找到极值点，比较轮轨法向间隙的最小极值点的大小，若相等，则为接触点的位置，否则旋转侧滚角，经过迭代直到极值点相等，这样就计算得到轮轨接触点[17]。

1.2 法向接触应力的计算

由上得到的轮轨接触点的坐标，引入轮轨接触实际嵌入量h0，J. Piotrowski[18]

经过计算研究得出h0=0.55h，h为轮轨法向弹性压缩量的最大值，在这里，利用这一嵌入量可得到实际的接触斑。对于出现多点接触的情况，其它接触点同时也分担一部分轮轨力，嵌入量h0就不够准确，此时我们需要对其进行修正。

在根据法向间隙获得接触斑大致范围，如果接触斑是关于车轴与主轮廓线形成的面对称的(忽略摇头角对其接触斑形状的影响)，则接触斑几个半轴的长度分别为

OA=b1，OB=b2，OC=OD=a，见图3，面ABCD为接触斑是关于AB轴对称，故将接触斑看作左右两个半椭圆的连接[15]。

根据Hertz轮轨接触理论，接触应力呈椭圆分布[15]，假设最大压缩变形处的接触应力为p0。将Hertz理论推广到非椭圆斑上，得到接触斑上各点的接触应力为

p(x,y)=

( 2 )

若轮轨法向力N已知，则接触斑上的最大法向接触应力p0为

p0=

( 3 )

把p0代入式( 2 )中就可以得到任意点的接触应力。

2 钢轨接触热效应分析

2.1 热流密度

轮轨间相对滑动摩擦热流密度确定如下

q(x,y)=μp(x,y)v

( 4 )

式中：q(x,y)为热流密度,J/(m2·s)；μ为摩擦系数；p(x,y)为接触区法向压力分布；