三维非赫兹滚动接触理论在轮轨滚动接触中的应用

应用Ｋａｌｋｅｒ的三维弹性体非赫弱滚动接触理论，将运动于直线轨道上的单轮对蠕滑率／力关系，建立在ＴＰＬＲ（非赫兹轮轨蠕滑力数表）中，分析了不同正压力所对应的蠕骨力值以及用表中的数值修正的蠕滑力值的精度。简单介绍了该表的使用方法。     

轮对摇头运动对轮轨滚动接触蠕滑率／力的影响

用数值分析方法分析了单轮对的摇头运动对其左右轮轨滚动接触斑上蠕滑率／力的影响。在轮轨滚动接触蠕滑率／力关系分析方面，利用了Ｋａｌｋｅｒ的三维弹性体非赫兹滚动接触计算模型。通过分析计算可知，轮对摇头角运动参量是影响轮轮之间横向蠕滑力的主要因素。     

三维非赫兹滚动接触理论在 轮轨滚动接触中的应用* ——TPLR 的编制及应用

应用Kalker 的三维弹性体非赫弱滚动接触理论, 将运动于直线轨道上的单轮对蠕滑率 ?力关系, 建立在TPLR (非赫兹轮轨蠕滑力数表) 中。分析了不同正压力所对应的蠕 滑力值以及用表中的数值修正的蠕滑力值的精度。简单介绍了该表的使用方法。      

几种铁路轮轨蠕滑力计算方法的比较

在对各种轮轨蠕滑计算模型进行归类的基础上,详细地讨论了Ｋａｌｌｅｒ教授提出的ＵＳＥＴＡＢ插值表法的结构和应用的可能途径,并对几种典型的方法作了分析比较。结果表明不同模型之间的差异随不同的轮轨接触条件而变化。计算结果可以作为选择轨蠕滑模型的依据。     

轮轨接触的有限元研究

用有限元参数二次规划法来分析轮轨的接触问题，提出了研究轮轨关系问题的新方法；同时给出了用此方法计算轮轨接触区内力的结果，并指出了各接触点对的接触状态和受力情况．     

基于mixed Lagrangian/Eulerian方法的轮轨滚动接触特性分析

为了缩短有限元方法显式计算轮轨滚动接触的时间,采用mixed Lagrangian/Eulerian方法建立了轮轨滚动接触有限元模型.应用该模型对轮轨接触区的单元进行细化,计算了列车在启动、运行和制动工况下轮轨的接触特性.计算结果表明:不同工况下,轮轨滚动接触区最大Mises应力、最大接触应力和接触斑面积等法向特性变化幅度均在2%以内,但接触区纵向截面中Mises应力分布及纵向剪切应力分布有较大变化;启动和制动工况下,最大Mises应力和最大纵向剪切应力位置均比自由滚动时更接近于轮轨表面;不同工况下,摩擦力大小和方向发生变化,在列车牵引和制动工况中,摩擦力达到极限时轮轨间出现完全滑动,摩擦力方向与滑动方向相反,且不同速度等级下的纵向摩擦力变化幅度也在2%以内.      

短波波磨状态的轮轨纵向蠕滑力特性

为了对具有简谐波形的钢轨短波波磨进行分组与分析轮轨非稳态滚动接触的纵向蠕滑力特性, 引入了波磨深度指数与波长比, 采用Kalker三维滚动接触理论计算了车轮的纵向蠕滑力, 并与采用稳态滚动理论计算结果进行了对比, 使用频率响应的系统辨识法对纵向蠕滑力的波动分量进行了拟合, 在短波波磨等深度指数条件下, 用波长比的二阶传递函数描述了轮轨纵向蠕滑力的波动分量与稳态理论波动分量之间的关系, 使用传递函数, 由稳态纵向蠕滑力的波动分量计算了非稳态纵向蠕滑力的波动分量, 进而计算了非稳态的纵向蠕滑力。计算结果表明: 在小蠕滑条件下, 由Kalker三维滚动接触理论计算出的纵向蠕滑力的波动分量随着波长比的变化产生明显的幅值衰减和相位滞后, 波长比越大, 幅值衰减越大, 相位滞后越多, 而稳态滚动理论的计算结果与波长比无关。由传递函数和Kalker数值理论计算的纵向蠕滑力的时域波形、频域幅值谱和相位谱相同。     

一个非线性的非稳态蠕滑力计算模型

本文讨论了由轨道对车轮的简谐动态输入引起的线性化简谐蠕滑力的影响。用弹性理论并空间假设下的数值解法。同时考虑接触区内粘着区和滑动区的非稳态变化和切向面力的非稳态变化，反复迭代，计算了非稳态的蠕滑系数，得出线性的非稳态蠕滑率－力关系，并将其与沈志云－赫居里克－埃尔金斯理论结合，推广到了非线性的惰形。     

高速列车轮轨弹塑性滚动接触理论和数值方法的研究

课题组根据2004年度项目工作要点：在建立轮对和钢轨处于刚性约束条件下的运动方程和轮轨接触界面的滑动方程、建立二维弹塑性滚动接触理论模型、轮轨二维非稳态滚动接触条件下材料棘轮效应对轮轨滚动接触参与应力应变分析、二维滚动接触条件下滚滑温度分析计算、二维滚动接触条件下粗糙度和水膜对轮轨滚动接触影响分析等方面开展了研究。同时，     

基于轮轨接触的高速铁路钢轨磨耗量

为预测高速铁路钢轨的磨耗量,建立了轨道结构静力学有限元模型和动力分析模型,基于Archard磨耗理论从曲线半径、行车速度、轮轨横移量3个角度计算分析了钢轨磨耗量,利用垂直磨耗深度0.5mm的磨耗量为界反算出通过总质量.计算结果表明：曲线地段钢轨磨耗较为严重,垂直磨耗深度为0.5mm时,直线上通过的总质量为45.9～60.0 Mt,曲线上通过仅为22.9～29.9Mr;相同曲线半径条件下,单轮作用下的接触斑处钢轨磨耗量随着行车速度提高而增大;相同速度和曲线半径下,钢轨磨耗量随着轮轨横移量增大而增大.     

轮轨稳态滚动接触中的HellingerReissner变分原理

以ＨｅｌｌｉｎｇｅｒＲｅｉｓｓｎｅｒ 变分原理为基础,应用摄动拉格朗日技术,将轮轨视为有限弹性体,推导出轮轨稳态滚动接触的ＨｅｌｌｉｎｇｅｒＲｅｉｓｓｎｅｒ 变分原理。为进一步从数值方面精确计算轮轨中接触疲劳问题提供理论依据和计算基础     

轮轨接触应力与钢轨波磨分析

根据轮轨滚动接触中钢轨循环加载塑性累积和材料的Ratcheting效应，应用强化材料模型对钢轨内部的残余应力和累积变形进行了数值分析。分析结果表明钢轨材料的Ratcheting效应和轮轨接触应力的波动是钢轨表面剥离与压溃形波波磨产生的重要原因。     

货车采用不同轮径时轮轨接触应力分析

加大车辆轴重是提高铁路运输能力的重要方法，但这势必会引起轮轨接触应力的增加。本文从理论计算和实验分析两个方面，研究了加大车轮直径对改善轮轨接触应力的关系，计算及实验结果表明，将货车车轮直径由840mm加大至915mm是有利的。     

考虑表面微观粗糙度的轮轨接触弹塑性分析

采用接触单元方法,结合初应力法,采用测量获得的表面微观粗糙度,对轮轨弹塑性接触问题进行了研究,获得了轮轨的表面接触压力分布等结果。结果表明,轮轨表面微观粗糙度度使得接触区的峰值接触压力大大高于平整接触表面的接触压力,会造成轮轨表面出现塑性变形。对含表面粗糙度的轮轨接触问题,进行弹塑性分析是必要的。     

轮轨两点接触的简易计算方法

本文给出轮轨两点接触判断和计算的简化方法,编制了电子计算机程序。实例计算表明,该法简便直观、计算速度快,精度可满足机车车辆动力学研究的需要。     

一种轮轨接触几何算法

提出并实现一种轮轨接触几何算法，可以检测铁道车辆系统动力学仿真在线计算时轮对与钢轨的刚性单接触斑、多接触斑和跳离情形．所得刚性接触斑可以为Hertz理论提供刚性穿透量和曲率，为非Hertz理论提供接触区域的法向间隙，为动力学仿真提供接触中心位置和法向方向．轮对与钢轨的计算机三维图形显示表明该算法是有效的．该算法已用于开发空间耦合的铁道车辆系统动力学仿真软件．     

基于轮轨接触参数的固定辙叉设计方法

为了更好地设计固定辙叉结构,改进轮对与辙叉间的相互作用关系,提出了以接触参数为基础的设计方法,并对60 kg/m钢轨12号固定辙叉心轨及翼轨顶面的各部位尺寸进行优化.结果表明:将此固定辙叉心轨顶宽20 mm处降低值增大至4 mm,取消翼轨顶面1∶20横坡,并将翼轨顶面抬高范围改为从辙叉咽喉至心轨顶宽50 mm处,可较好的改善其轮轨关系,验证了本文设计方法的可行性和有效性.     

连续测量轮轨力的测力轮对技术

本文介绍了用标准的机车辆轮对作为力传感器来连续测量轮轨相互作用的测力轮对技术，并以作者研究的测力轮对为例，介绍了选择最佳贴片位置的方法、电桥设计的一般原则及试验结果处理等技术。