摘 要:文章首先介紹了基于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.