(完整版)总结轮轨关系

轮轨关系轨道车辆和线路的作用问题是铁路轮轨接触式运输的基本问题。发展重载运输必须解决好轮轨之间的动力作用，努力减轻重载列车与线路的动态作用。但由于轮轨关系自身的复杂性，目前的研究理论和模型仍然基于一些假设［1］：1)法向接触满足Hertz接触条件；轮轨接触副视为弹性半空间；接触表面是理想光滑连续的,而接触表面之间的“第三介质,'如水、油和其它污染物的影响被忽略；轮轨接触斑以外边界支撑和约束条件对轮轨接触行为的影响被忽略；高速轮轨滚动接触时的惯性力被忽略；不考虑温度的影响。上述几点假设是不符合实际但是理论的前提。轮轨关系的主要研究内容为轮轨接触几何的确定和轮轨滚动接触理论的应用。实际接触参数计算和列解微分方程的过程可简述如下：在某一瞬时位置确定轮轨接触点是关键，之后就可以在确定了接触点的基础上利用几何推导出各个重要的接触几何参数，如左右轮/轨在接触点的接触角5.5,LR左右轨在接触点处的钢轨顶面曲率半径p、P，左右轮在接触点处的踏面曲RRRL率半径P、P，左右轮实际滚动半径r、r，轮轨接触时的侧滚角9，轮对WRWLRLk中心的上下位移z，其中变量为轮对相对轨道的横移量和摇头角y、屮。利用kww已求得的接触参数和Hertz接触理论公式计算出接触椭圆的长短半轴，从而确定轮轨接触斑。然后利用接触椭圆的长短半轴长和查表得到的kalker系数及材料常数计算得到蠕滑系数，之后再通过实际速度和纯滚动速度计算出蠕滑率，将二者带入蠕滑力公式求得蠕滑力。最后就可以列出含有蠕滑力，悬挂力，惯性力的运动微分方程，利用计算机求解得到位移、速度、加速度和相关模态值。最初进行轮轨接触几何关系研究并确定接触参数的实用方法有两种：一种是圆弧形截面模型，一种是任意截面模型。前者可直观的用数学解析的方法确定其几何关系，后者是数值方法，需编程实现。前者在综述中提到；现重点论述后者，它是一种通用性很强的求解轮轨接触几何的数值方法。任意截面轮廓形状轮轨接触几何求法：(以二维模型为例)假设及准则：1•轮轨为刚性，互不嵌入；2两侧轮轨同时接触；3•轮轨上的接触点具有相同的空间位置，且接触点处轮轨具有公切面。任何轮轨外形都可以用一系列离散坐标来表示轮轨外形，然后用一条通过各坐标点的拟合曲线来表示其外形的描述函数。如果要求精度高，则离散点可以密集些。下面具体介绍求解任意截面轮廓形状轮轨接触几何的方法和步骤：I•轮轨外形拟合曲线如图所示，过轮对中心线做一个垂向平面切割轮对，车轮外部轮廓与切割面得交线称为主轮廓线，即车轮踏面外形线。主轮廓线可以通过实测或计算轮廓线上各点坐标用三次样条函数拟合，得到描述外形线的样条函数W和W。同理，TOC\o"1-5"\h\zR(y)L(y)找到切割面与轨顶交线上的离散点可拟合成R和R样条函数。这样就得到R(y)L(y)了左右轮、轨的四个主轮廓线的样条函数W、W、R、R。R(y)L(y)R(y)L(y)注：三次样条函数曲线S的性质：(1)在每个子区间［1,1］(k=1,2・・・.n),S不超过三次；k-1k(1)S=y(k=0,1....,n)(1kS在［a,b］区间上两次连续可微。(i)图确定轮轨接触点原理II.确定轮轨接触点的原理和接触参数表达式：(无摇头的二维情况)已知轮轨轮廓曲线的四条样条函数，就确定了轮轨轮廓上的任一点坐标。在某横

移值y下，可依据样条函数求得轮轨轮廓的垂直距离z；变化y的值，增量取决ii于切割面（图示垂直于纸面）的间距，得到不同的z（i=l,2,3・・・.m）。在zii（i=l,2,3・・・.m）中，找到其中最小的z。比较左右轮轨的z和z，若二minminLminR者相等，车轮垂直平行下落时，两最小距离点z和z即为轮轨接触点。若minLminR不等，则要调整车轮侧滚角e—个微小的角度丫重复以上各步，直到z=zkiminLminR确定接触点。最终可以通过左轮左轨，右轮右轨是否分别具有公切面验证接触点的真假。接触点确定，接触点的横坐标y和y就确定了。之后就可以确定出各个接LminRmin触参数如下：1•侧滚角e，直接由迭代求得e+sY；kkk0i实际车轮滚动圆半径r和r，可由y和y确定；RLLminRmin3•轮对中心垂向位移量z，可由r、r和e确定；kRLk4•钢轨接触角（接触点切面与水平面夹角）A,A:LR左轨：A左轨：A=arctg—[R]LdyL(yLmin)右轨：A=arctg—[R]RdyR(yRmin)5•车轮接触角（接触点切面与轮对中心线间夹角）5，5LR左轮：5左轮：5=arctg—[W]一0LdyL(yLmin)k右轮：5=arctg—[W]+eRdyR(yRmin)k6•车轮接触点处踏面曲率半径p，p:WLWR左轮：{1+Q—[W])}左轮：{1+Q—[W])}3/2dyL(yLmin)d2dy2[W]L(yLmin)右轮：{1+Cd[W])2}3/2dyR(yRmin)d2dy2[WR(yRmin)]7•钢轨在接触点处轨面曲率半径p,pRRr左轨:{1+C—[R])2左轨:{1+C—[R])2}3/2dL(y£mm)二[R]dy2L(yLmJ右轨:Rr{1+C—[R)}3/2dyR(yRmJd2dy2[RR(yRmin)轮轨三维几何接触关系参考文献[2],当轮对有摇头时，轮对上的接触点已不在主轮廓线所在的垂平面内，因此要在车轮的整个踏面的轮廓面上来寻找轮轨接触点。III•轮轨滚动接触及蠕滑力计算（轮轨动力作用部分）：为了求得运动微分方程中蠕滑力项，首先必须获知轮轨接触几何关系，在给定的轮轨截面参数下，由轮对横移y及摇头角屮马可确定接触点位置及接触点处的ww曲率半径等（上步已经求得）。第二步根据左右两个接触点的位移y及速度y，ww车辆实际速度和纯滚动速度等求岀蠕滑率，根据Hertz接触理论求岀接触斑尺寸，并结合材料属性利用Carter滚动接触理论或Kalker滚动接触理论计算蠕滑系数。最后代入蠕滑力公式求岀蠕滑力项。下面具体介绍第二步内容：车轮在牵引力的作用下要向前运动，如果轮轨间粘着力不足，则车轮在钢轨上滑行，如果粘着力充分，则车轮在钢轨上滚动。滚动的车轮在钢轨上走行的距离要小于没有弹性变形的纯理论滚动所走的距离，也就是车轮的实际表面速度低于理论表面速度，这种现象称为轮轨蠕滑。轮轨间岀现蠕滑的三个条件：轮轨为弹性体，车轮和钢轨之间有一定的正压力和切向力使车轮沿钢轨滚动。缺少三者中的任何一个，均不会岀现蠕滑。车轮实际走行速度与理论纯滚动速度之差称为蠕滑速度，蠕滑速度与纯滚动速度之比称为蠕滑率；使轮轨产生蠕滑的切向力称为蠕滑力。蠕滑力由蠕滑系数和蠕滑率求得。具体的：

图蠕滑力F和蠕滑率£之间关系先要依据Hert接触理论［3］确定轮轨接触斑的尺寸，即接触椭圆的长短半轴a和b的大小。33n33nN(k+k)12—4k3图椭圆接触斑示意图Hertz推导出的接触椭圆长短半轴公式为:4k参数k、k、k、k，1234P定义为:11111111.=t[+++],2rprpwwRR71—G2k=w-i兀Ew

71—G2k=r-2兀ER111111111k=—[(+)2+(+)2+2(—)(—)cos刘]42rprprprpwwRRwwRRP=arccos(k/k);43其中：N其中：N为轮轨接触斑处法向力,r车轮实际滚动圆半径，p车轮踏面轮廓在wwp钢轨轮廓面在R接触点处的曲率半径，rp钢轨轮廓面在RE、E车轮和钢轨材wRE、E车轮和钢轨材wRwR料的杨氏弹性模量；m、n是Hertz理论给出的系数，可由P查表知，屮为r、rwR法平面夹角即轮对摇头角（实际中cos2屮约等于1）。

以上求得了接触椭圆长短半轴a,b。为以后求蠕滑系数奠定了基础。在最后的蠕滑力计算中将详列蠕滑系数的表达式。求解蠕滑力的另外一个参数项蠕滑率可用下图表达，坐标C-123是以接触中心C在接触斑上建立的坐标系，如图所示。车轮上接触斑沿C1,C2轴的平动速度为v、v，绕C3轴的转动速度为0；钢轨上接触斑沿C1,C2轴的平动速度为v、w1w2w3R1TOC\o"1-5"\h\zv，绕C3轴的转动速度为0。定义蠕滑率：车轮实际前进速度与车轮纯滚动R2R3前进速度差与车轮实际前进速度（有的书用纯滚动速度［5］）的比值。具体的：沿C1轴纵向蠕滑率2二2（舄1_寫）v+vR1w1沿C2轴横向蠕滑率2二2（vR2_vw2）v+vR1w1沿C3轴自旋蠕滑率2二2（0w3_0R3）v+vR1w1其中，^1+vw1为车轮沿钢轨运行的平均速度，实用时用运行速度代替纯滚动速2度。图蠕滑率计算示意图有Kalker线性蠕滑理论知,蠕滑力（力矩）与蠕滑率成正比且反向:图蠕滑率计算示意图有Kalker线性蠕滑理论知,蠕滑力（力矩）与蠕滑率成正比且反向:纵向蠕滑力F=-f21111横向蠕滑力F=-f-fg222233自旋蠕滑力矩M二fg-fg322333了=EabCTOC\o"1-5"\h\z1111f=EabC公式中的蠕滑系数f的确定：ij<2222，其中a公式中的蠕滑系数f的确定：ijf=f=E(ab)3/2C233223f=E(ab)2C3333长短半轴长，E杨氏模量，C是查表得到的Kalker系数，与材料泊松比和a/bij的值有关。大量实验证明，轮轨蠕滑力和蠕滑率之间只有在蠕滑率不大情况下二者才成线性关系。确定轮轨非线性蠕滑力和蠕滑率的关系最简便实用的方法为经验修正法。它是利用线性Kalker理论计算出F和F求二者的合力F=JF2+F2，修正F'得1