轮轨接触应力的有限元计算.pdfVIP

　试验研究 铁道车辆 　第 45 卷第 6 期 2007 年 6 月 　 文章编号 (2007) 轮轨接触应力的有限元计算 1 2 张 澎 湃 , 井 秀 海 ( 1. 中国北车集团四方车辆研究所 研究试验部 ,山东 青岛 26603 1 ; 2 . 北京市地铁运营有限公司车辆厂 筹备组 ,北京 100077) 摘 　要 :为了得到精确的数值解 ,应用弹塑性理论及大型有限元软件 AN SYS 对轮轨接触应力进行了有限元计算 , ( ) 同时对建模 、单元选择 、网格划分 、参数确定等计算过程进行了说明 ,得到了不同轴重 1 t ～30 t 时轮轨塑性区尺寸 、接 触斑形状 、Mi ses 当量应力最大值及轮轨接触面压应力随轴重的变化规律 。 关键词 :轮轨关系 ;接触应力 ;有限元法 ;弹塑性变形 中图分类号 :U2 11. 5 　　　文献标识码 :A 　　对轮轨间的相互作用 , 国内外学者进行过大量的 元算法的目的是克服一阶单元所特有的剪切 自锁效 研究 ,并取得了不少研究成果 ,但其中涉及的计算多数 应 ,提高计算精度 。 ( ) [ 14 ] 车轮踏面和钢轨的相互作用可以看作柔体柔体 基于 非 赫兹接触理论 ,计算结果为弹性解析解 。 实际上 ,轮轨接触时轨头和踏面已经发生塑性变形 ,因 的接触问题 。AN SYS 面面接触算法支持有大滑动和 ( ) 摩擦的柔体柔体大变形接触问题的计算 ,能够为工程 此基于 非 赫兹接触理论的计算结果不全面 。也有根 据弹塑性理论和有限元法进行计算的实例[ 58 ] ,不过模 计算提供很好的计算结果 ,通过设置接触对建立不同 型均进行过简化 。 物体之 间的接触关系 。本文选用接触单元 CON 为了真实地反映轮轨接触状态 ,加深对踏面剥离 、 TA 174 和 目标单元 TA R GE170 组成接触对 。接触算 轮轨接触疲劳等问题的研究 ,本文应用弹塑性理论及 法选用扩展的拉格朗日算法 ,接触对间的摩擦和运动 大型通用有限元软件 AN SYS 对轮轨弹塑性接触应力 状态根据库仑摩擦模型确定 ,轮轨接触表面作用模式 进行了有限元计算 ,并对计算结果进行了分析 。 为法向单边接触 。 1 . 2 　网格划分和模型参数确定 1 　轮轨接触应力的有限元计算 1 . 2 . 1 　网格划分 1 . 1 　单元及接触算法选取 由于轮轨接触面积远小于接触表面的曲率半径 , 轮轨接触时轨头和踏面位置发生了塑性变形 , 由 可知接触区应力远大于非接触区 , 即接触区内存在明 于金属塑性变形具有不可压缩性 , 因此模拟金属塑性 显的应力集中。在这种情况下 ,计算结果的精度高度 将增加对单元的运动学约束 ,并要求在单元积分点处 依赖于网格划分情况 , 即接触区网格要细划到一定程 的体积要保持常数 。