摘 要：文章首先介紹了基于Hertz理論的輪軌滾動接觸理論，通過ABAQUS有限元軟件相結(jié)合從不同輪對橫移量對滾動接觸狀態(tài)下輪軌的滾動接觸特性進行分析。結(jié)果表明：輪對橫移量主要影響輪軌滾動接觸的法向特性；最大Mises應力以及最大法向接觸應力出現(xiàn)在無橫移量的工況下。

關(guān)鍵詞：輪軌滾動接觸特性；多體動力學仿真；有限元+

1 基于Hertz理論的輪軌滾動接觸理論

Hertz運用彈性力學理論對兩個彈性體接觸問題進行了研究，并提出了彈性接觸Hertz理論[1]。

Hertz還闡述了物體表面的法向位移必須滿足式（1）的條件：

（1）

式中，δ=δ1+δ2為兩個物體內(nèi)部距離較遠處點的移動位移；uz1，uz2為兩個物體的表面法向位移； Ax2+By2為接觸面初始的間隙。

輪軌接觸的曲面初始間隙的函數(shù)常數(shù)A和B可由下列方程確定：

（2）式與（3）式中，R1′-輪圓的半徑；R1″-踏面橫斷面的外形的半徑；R2′-鋼軌縱向曲率半徑值；R2″-鋼軌橫斷面外形的徑值。

2 基于ABAQUS有限元的輪軌滾動接觸分析

2.1 ABAQUS有限元模型

為了協(xié)調(diào)計算的效率和計算的精度，模型對鋼軌軌頭部分單元進行了細化。最后形成結(jié)點總數(shù)為37118、單元總數(shù)為24844的鋼軌模型和結(jié)點總數(shù)為8064、單元總數(shù)為5832的車輪模型。

采用庫倫摩擦模型定義切向關(guān)系，定義“硬”接觸作為法向關(guān)系；從面的選擇上，將輪圓面定義為主面，將軌頭頂面定義為從面；輪軌的接觸類型定義為面面接觸；將輪軌接觸對之間定義為“小滑移”且設置接觸滑動容差為0.006，以滿足輪軌滾動接觸而產(chǎn)生的彈性變形，會引起輪軌接觸對在接觸面上出現(xiàn)小滑動。

2.2 不同車輪橫移量的影響

橫移量是指輪對的垂向中心線偏離軌道中心線的距離。對輪對在30.920kN作用下、車輪角速度為129.2rad/s的輪軌接觸的特性隨輪對的橫移量變化的情況進行分析，此時摩擦系數(shù)取為0.2，輪對沖角取為0mard，橫移量取值分別為0、1、2、3、4、5mm。輪軌接觸斑內(nèi)單元節(jié)點上接觸斑形狀、接觸斑面積、最大應力和最大法向接觸應力如表2所示。

由表2我們可以得出，輪軌滾動接觸的特性與輪軌橫移量有較大的關(guān)系，隨著輪對橫移量的逐漸增大，接觸位置隨之遷移，接觸斑的面積逐漸增大，接觸區(qū)域上的最大Mises應力整體趨勢上隨著輪對橫移量的增加而減?。唤佑|區(qū)域上的最大法向接觸應力整體趨勢上隨著輪對橫移量的增加而減??；Mises應力、法向接觸應力、接觸斑內(nèi)橫向剪切應力、縱向剪切應力在無輪對橫移量時出現(xiàn)最大值。

從圖1、2和3中可以看出，最大Mises應力出現(xiàn)于鋼軌接觸面的表層，隨著輪軌橫移量的增大，最大Mises逐漸減小，Mises應力分布區(qū)域范圍也趨于減小。不難發(fā)現(xiàn)的是，輪軌發(fā)生橫移時，輪軌接觸斑的位置、尺寸以及輪軌滾動接觸的法向特性各應力、切向特性各應力的應力值都會發(fā)生變化。

3 結(jié)束語

文章利用ABAQUS有限元軟件，對輪軌滾動接觸中鋼軌表面接觸特性進行了分析，并實現(xiàn)了有關(guān)結(jié)果的三維立體顯示，最終主要獲得了以下結(jié)論：（1）最大Mises應力以及最大法向接觸應力出現(xiàn)在無橫移量下。（2）最大Mises應力出現(xiàn)于鋼軌接觸面的表層，且鋼軌接觸斑上的最大Mises應力隨著輪軌橫移量的增大逐漸減小，且應力分布區(qū)域范圍也趨于減小。

參考文獻

[1]Gladwell G M L.經(jīng)典彈性理論中的接觸問題（中譯本）[M].北京理工大學出版社，1991.

[2]郭俊.輪軌滾動接觸疲勞損傷機理研究[D].西南交通大學，2006.

[3]金學松，劉啟躍.輪軌摩擦學[M].中國鐵道出版社包鋼科技，2004.

[4]閻國臣，何慶復.車輪滾動接觸疲勞研究[J].鐵道機車車輛，2002，11（004）：18-19.

[5]金學松，沈志云.輪軌滾動接觸疲勞問題研究的最新進展[J].鐵道學報，2001，23（002）：92-100.