肖乾 徐紅霞 王成國 張海
在使用有限元方法分析輪軌滾動接觸動態問題時,可采用中心差分法顯式計算[1],將車輪行進過程中輪軌接觸部分全部進行單元細化,且鋼軌模型較長,這不僅使得有限元模型單元規模龐大,也給計算帶來困難,還將帶來不可估量的顯式計算時間問題[23].因此,為了降低輪軌滾動接觸有限元模型的計算難度,國內外學者采取子結構法[4]、移動分布載荷等方法[56]分析輪軌接觸問題,但這又使計算精度受到影響.本文采用mixed Lagrangian/Eulerian方法[7]建立輪軌滾動接觸三維有限元模型,進而完成其有限元分析.這種方法分析輪軌滾動接觸動態問題可真實反應輪軌接觸狀態,克服了顯式計算的弊端,只需要對輪軌接觸區進行網格細化,采用隱式迭代方法求解,能高效準確的完成輪軌不同工況下的滾動接觸特性分析.