轮轨关系
轨道车辆和线路的作用问题是铁路轮轨接触式运输的基本问题。发展重载运输必须解决好轮轨之间的动力作用，努力减轻重载列车与线路的动态作用。但由于轮轨关系自身的复杂性，目前的研究理论和模型仍然基于一些假设[1]： 1）法向接触满足Hertz 接触条件； 2）轮轨接触副视为弹性半空间；
3） 接触表面是理想光滑连续的,而接触表面之间的 “第三介质 ,’ 如水、油和其它污染物的影响被忽略；
4) 轮轨接触斑以外边界支撑和约束条件对轮轨接触行为的影响被忽略； 5) 高速轮轨滚动接触时的惯性力被忽略； 6) 不考虑温度的影响。
上述几点假设是不符合实际但是理论的前提。
轮轨关系的主要研究内容为轮轨接触几何的确定和轮轨滚动接触理论的应用。实际接触参数计算和列解微分方程的过程可简述如下：
在某一瞬时位置确定轮轨接触点是关键，之后就可以在确定了接触点的基础上利用几何推导出各个重要的接触几何参数，如左右轮/轨在接触点的接触角?L、?R，左右轨在接触点处的钢轨顶面曲率半径?RR、?RL，左右轮在接触点处的踏面曲率半径?WR、?WL，左右轮实际滚动半径rR、rL，轮轨接触时的侧滚角?k，轮对中心的上下位移zk，其中变量为轮对相对轨道的横移量和摇头角yw、?w。利用已求得的接触参数和Hertz接触理论公式计算出接触椭圆的长短半轴，从而确定轮轨接触斑。然后利用接触椭圆的长短半轴长和查表得到的kalker系数及材料常数计算得到蠕滑系数，之后再通过实际速度和纯滚动速度计算出蠕滑率，将二者带入蠕滑力公式求得蠕滑力。最后就可以列出含有蠕滑力，悬挂力，惯性力的运动微分方程，利用计算机求解得到位移、速度、加速度和相关模态值。 最初进行轮轨接触几何关系研究并确定接触参数的实用方法有两种：一种是圆弧形截面模型，一种是任意截面模型。前者可直观的用数学解析的方法确定其几何关系，后者是数值方法，需编程实现。前者在综述中提到；现重点论述后者，它是一种通用性很强的求解轮轨接触几何的数值方法。
任意截面轮廓形状轮轨接触几何求法：（以二维模型为例）
假设及准则：1.轮轨为刚性，互不嵌入；2两侧轮轨同时接触；3.轮轨上的接触点具有相同的空间位置，且接触点处轮轨具有公切面。
任何轮轨外形都可以用一系列离散坐标来表示轮轨外形，然后用一条通过各坐标点的拟合曲线来表示其外形的描述函数。如果要求精度高，则离散点可以密集些。下面具体介绍求解任意截面轮廓形状轮轨接触几何的方法和步骤： I.轮轨外形拟合曲线
如图所示，过轮对中心线做一个垂向平面切割轮对，车轮外部轮廓与切割面得交线称为主轮廓线，即车轮踏面外形线。主轮廓线可以通过实测或计算轮廓线上各点坐标用三次样条函数拟合，得到描述外形线的样条函数WR(y)和WL(y)。同理，找到切割面与轨顶交线上的离散点可拟合成RR(y)和RL(y)样条函数。这样就得到了左右轮、轨的四个主轮廓线的样条函数WR(y)、WL(y)、RR(y)、RL(y)。 注：三次样条函数曲线S(?)的性质：
i.在每个子区间[?k?1,?k] (k=1,2….n) , S(?)不超过三次； ii. S(?)?yk (k=0,1....,n)
kiii. S(?)在[a,b]区间上两次连续可微。
图 确定轮轨接触点原理
II.确定轮轨接触点的原理和接触参数表达式：（无摇头的二维情况）
已知轮轨轮廓曲线的四条样条函数，就确定了轮轨轮廓上的任一点坐标。在某横
移值yi下，可依据样条函数求得轮轨轮廓的垂直距离zi；变化yi的值,增量取决于切割面（图示垂直于纸面）的间距，得到不同的zi（i=1,2,3….m）。在zi（i=1,2,3….m）中，找到其中最小的zmin。比较左右轮轨的zminL和zminR，若二者相等，车轮垂直平行下落时，两最小距离点zminL和zminR即为轮轨接触点。若不等，则要调整车轮侧滚角?k一个微小的角度?i重复以上各步，直到zminL=zminR确定接触点。最终可以通过左轮左轨，右轮右轨是否分别具有公切面验证接触点的真假。
接触点确定，接触点的横坐标yLmin和yRmin就确定了。之后就可以确定出各个接触参数如下：
1.侧滚角?k，直接由迭代求得?k??k0???i；
2.实际车轮滚动圆半径rR和rL，可由yLmin和yRmin确定； 3.轮对中心垂向位移量zk，可由rR、rL和?k确定； 4.钢轨接触角（接触点切面与水平面夹角）?L, ?R:
ddyddy左轨：?L?arctg[RL(yLmin)] 右轨：?R?arctg[RR(yRmin)]
5.车轮接触角（接触点切面与轮对中心线间夹角）?L，?R
ddyddy左轮：?L?arctg[WL(yLmin)]??k 右轮：?R?arctg[WR(yRmin)]??k
6.车轮接触点处踏面曲率半径?W，?W：
LR{1?（ddyd22[WL(yLmin)]）}[WL(yLmin)]23/2{1?（ddyd22[WR(yRmin)]）}[WR(yRmin)]23/2左轮：?W?L 右轮：?W?R
dydy7.钢轨在接触点处轨面曲率半径?R,?R
LR{1?（ddyd22[RL(yLmin)]）}[RL(yLmin)]23/2{1?（ddyd22[RR(yRmin)]）}[RR(yRmin)]23/2左轨：?R?L 右轨：?R?R
dydy轮轨三维几何接触关系参考文献[2]，当轮对有摇头时，轮对上的接触点已不在主轮廓线所在的垂平面内，因此要在车轮的整个踏面的轮廓面上来寻找轮轨接触点。
III .轮轨滚动接触及蠕滑力计算（轮轨动力作用部分）：
为了求得运动微分方程中蠕滑力项，首先必须获知轮轨接触几何关系，在给定的轮轨截面参数下，由轮对横移yw及摇头角?w马可确定接触点位置及接触点处的
?w，曲率半径等（上步已经求得）。第二步根据左右两个接触点的位移yw及速度y车辆实际速度和纯滚动速度等求出蠕滑率，根据Hertz接触理论求出接触斑尺寸，并结合材料属性利用Carter滚动接触理论或Kalker滚动接触理论计算蠕滑系数。最后代入蠕滑力公式求出蠕滑力项。下面具体介绍第二步内容： 车轮在牵引力的作用下要向前运动，如果轮轨间粘着力不足，则车轮在钢轨上滑行，如果粘着力充分，则车轮在钢轨上滚动。滚动的车轮在钢轨上走行的距离要小于没有弹性变形的纯理论滚动所走的距离，也就是车轮的实际表面速度低于理论表面速度，这种现象称为轮轨蠕滑。轮轨间出现蠕滑的三个条件：轮轨为弹性体，车轮和钢轨之间有一定的正压力和切向力使车轮沿钢轨滚动。缺少三者中的任何一个，均不会出现蠕滑。车轮实际走行速度与理论纯滚动速度之差称为蠕滑速度，蠕滑速度与纯滚动速度之比称为蠕滑率；使轮轨产生蠕滑的切向力称为蠕滑力。蠕滑力由蠕滑系数和蠕滑率求得。具体的：
图 蠕滑力F和蠕滑率?之间关系
先要依据Hertz接触理论[3]确定轮轨接触斑的尺寸，即接触椭圆的长短半轴a和b的大小。
图 椭圆接触斑示意图
Hertz推导出的接触椭圆长短半轴公式为：
3a?m3?N(k1?k2)4k33 b?n3?N(k1?k2)4k3 参数k1、k2、k3、k4，?定义为：
1??w2k1??Ew121rw , k2?1??R2?ER1 , k3?11111[???], 2rw?wrR?R11rR1k4?[(?1?w)?(21rR??R)?2(21rw??w)(??R)cos2?]
?=arccos（k4/k3）;
其中：N为轮轨接触斑处法向力，rw车轮实际滚动圆半径，?w车轮踏面轮廓在接触点处的曲率半径，rR钢轨在接触点处的滚动滚动圆半径，?R钢轨轮廓面在接触点处曲率半径，?w、?R车轮和钢轨材料的泊松比，Ew、ER车轮和钢轨材料的杨氏弹性模量；m、n是Hertz理论给出的系数，可由?查表知，?为rw、rR法平面夹角即轮对摇头角（实际中cos2?约等于1）。