一种求解三维弹塑性轮轨法向接触问题的新方法  

A new approach for solving three-dimensional elasto-plastic wheel-rail normal contact problems

作　　者：赵吉中 伏培林 苗鸿臣 温泽峰 陈嵘 阚前华 Jizhong Zhao;Peilin Fu;Hongchen Miao;Zefeng Wen;Rong Chen;Qianhua Kan(Applied Mechanics and Structure Safety Key Laboratory of Sichuan Province,School of Mechanics and Aerospace Engineering,Southwest Jiaotong University,Chengdu 611756,China;State Key Laboratory of Traction Power,Southwest Jiaotong Universiy,Chengdu 610031,China;Key Laboratory of High-speed Railway Engineering,Ministry of Education,Southwest Jiaotong University Chengdu 610031,China)

机构地区：[1]Applied Mechanics and Structure Safety Key Laboratory of Sichuan Province,School of Mechanics and Aerospace Engineering,Southwest Jiaotong University,Chengdu 611756,China [2]State Key Laboratory of Traction Power,Southwest Jiaotong Universiy,Chengdu 610031,China [3]Key Laboratory of High-speed Railway Engineering,Ministry of Education,Southwest Jiaotong University Chengdu 610031,China

出　　处：《Acta Mechanica Sinica》2023年第9期126-140,共15页力学学报（英文版）

基　　金：supported by the National Natural Science Foundation of China(Grant Nos.U21A20167,12072295,12192214,12192210);Applied Basic Research Project of Sichuan Province(Grant No.2021YJ0050);Fundamental Research Funds for the Central Universities(Grant No.2682021ZTPY098).

摘　　要：提出了一种求解三维弹塑性轮轨法向接触问题的新方法,用于求解弹塑性轮轨接触下的接触斑和法向接触压力分布。首先,该方法通过矩形网格离散可能的接触区域;然后,采用共轭梯度法对每个矩形网格中的法向接触压力进行迭代求解,获取弹性情况下的接触斑和法向接触压力分布;进而,通过提出的弹塑性修正方法,对每个矩形网格的尺寸和法向接触压力进行修正,得到弹塑性情况下的轮轨接触斑和法向接触压力分布;最后,基于不同横移量和法向合力下的有限元计算结果,对该方法进行验证.结果表明,该方法具有较高的计算精度和计算效率,可为轮轨接触问题的研究提供参考.A new approach for solving three-dimensional elasto-plastic wheel-rail normal contact problems is presented,which is used to obtain the normal contact pressure distribution in the wheel-rail contact patch and its interior.First,the potential contact region dispersed by rectangular blocks is given.Then,the normal contact pressure in each block is iterated using the conjugate gradient algorithm,and the contact patch and normal contact pressure distribution under elastic conditions are obtained.After that,based on the proposed normal contact pressure correction equation and grid size correction equation,the size and normal contact pressure of each block are corrected,and so the normal contact pressure considering plastic deformation is obtained.Finally,the wheel-rail normal contact pressure under different lateral displacements and normal forces is calculated using the finite element method,and the accuracy of the proposed method is verified.Results show that the computational efficiency of the proposed approach is much higher than that of the conventional finite element method,which can provide new inspiration for railwaydynamics computation.

关 键 词：轮轨接触 矩形网格 接触斑 共轭梯度法 三维弹塑性 接触压力分布 横移量 有限元计算 

分 类 号：U211.5[交通运输工程—道路与铁道工程]