基于mixed Lagrangian/Eulerianし椒ǖ穆止旃龆接触特性分析

肖乾　徐红霞　王成国　张海

在使用有限元方法分析轮轨滚动接触动态问题时，可采用中心差分法显式计算[1]，将车轮行进过程中轮轨接触部分全部进行单元细化，且钢轨模型较长，这不仅使得有限元模型单元规模庞大，也给计算带来困难，还将带来不可估量的显式计算时间问题[23].因此，为了降低轮轨滚动接触有限元模型的计算难度，国内外学者采取子结构法[4]、移动分布载荷等方法[56]分析轮轨接触问题，但这又使计算精度受到影响.本文采用mixed Lagrangian/Eulerian方法[7]建立轮轨滚动接触三维有限元模型，进而完成其有限元分析.这种方法分析轮轨滚动接触动态问题可真实反应轮轨接触状态，克服了显式计算的弊端，只需要对轮轨接触区进行网格细化，采用隐式迭代方法求解，能高效准确的完成轮轨不同工况下的滚动接触特性分析.