緊急制動(dòng)下重載貨車車輪溫度場(chǎng)和應(yīng)力場(chǎng)分析

袁征　趙越超　張?chǎng)螘?/p>

摘 要：本文研究了25噸軸重的重載貨車在緊急制動(dòng)條件下車輪的瞬態(tài)溫度和應(yīng)力變化情況，建立了HESA型貨車車輪的三維有限元模型，分析了貨車車輪的溫度場(chǎng)和應(yīng)力場(chǎng)。通過(guò)踏面部分加載熱流密度，整體車輪加載對(duì)流換熱、熱輻射，得到了緊急制動(dòng)過(guò)程中踏面的溫度和應(yīng)力的變化趨勢(shì)，同時(shí)得到了車輪整體最高溫度、最高應(yīng)力及其各自出現(xiàn)的時(shí)間，為車輪的使用壽命和疲勞分析提供了理論基礎(chǔ)。

關(guān)鍵詞：鐵路貨車；緊急制動(dòng)；車輪；溫度場(chǎng)；應(yīng)力場(chǎng)

DOI：10.16640/j.cnki.37-1222/t.2016.09.002

車輪是鐵路貨車走行部的重要組成部件之一，它不但承受著貨車的自重與載重，而且還與鋼軌、閘瓦及環(huán)境介質(zhì)有著非常復(fù)雜的作用，承受著極其復(fù)雜的應(yīng)力和踏面制動(dòng)所產(chǎn)生的熱負(fù)荷[1-3]。隨著貨物列車的提速和軸重的增加，特別是25t軸重的貨車車輪的服役運(yùn)行和貨車運(yùn)行速度達(dá)到120km/h，列車的動(dòng)能也會(huì)急劇上升。由此，依靠閘瓦和車輪踏面的機(jī)械摩擦來(lái)實(shí)現(xiàn)制動(dòng)，必然在車輪踏面部位產(chǎn)生大量的熱量，這種熱載荷和機(jī)械載荷的聯(lián)合作用，將使得貨車車輪踏面受到較為嚴(yán)重的破壞[4]。因此，對(duì)于貨車車輪在踏面制動(dòng)情況下的溫度場(chǎng)和應(yīng)力場(chǎng)進(jìn)行深入研究具有重要的經(jīng)濟(jì)價(jià)值和理論意義。

有學(xué)者曾經(jīng)對(duì)貨車車輪熱損傷問(wèn)題進(jìn)行了研究，并對(duì)21t和23t軸重車輪的溫度場(chǎng)和熱應(yīng)力場(chǎng)進(jìn)行了模擬計(jì)算。王京波曾使用有限元分析軟件MARC建立快速貨車車輪二維有限元模型，并分別用數(shù)值方法和實(shí)驗(yàn)方法研究合成閘瓦對(duì)車輪的熱影響[5]。劉云曾對(duì)提速貨車車輪的溫度場(chǎng)及熱應(yīng)力場(chǎng)進(jìn)行了數(shù)值模擬，但其模型為二維模型，且只考慮了熱傳導(dǎo)[6]。侯耐在2011年曾對(duì)重載貨車踏面制動(dòng)熱負(fù)荷進(jìn)行了分析研究，對(duì)緊急制動(dòng)工況和長(zhǎng)大下坡道制動(dòng)工況下的溫度場(chǎng)以及熱應(yīng)力場(chǎng)進(jìn)行了模擬[7]。等等。

本文從我國(guó)鐵道車輛的實(shí)際情況出發(fā)， 選取了轉(zhuǎn)K6轉(zhuǎn)向架所使用的HESA型輾鋼車輪進(jìn)行建模，全面考慮了現(xiàn)今鐵路貨車不斷向“高速重載”發(fā)展的趨勢(shì)。通過(guò)實(shí)際緊急制動(dòng)過(guò)程的模擬分析，得到了車輪所承受的機(jī)械載荷和熱載荷，為車輪的疲勞分析提供了理論依據(jù)，有利于預(yù)防車輪損害，增長(zhǎng)車輪使用壽命，減少經(jīng)濟(jì)損失。

1 車輪邊界條件的確定

本文對(duì)軸重為25t，制動(dòng)初速度為120km/h，緊急制動(dòng)距離為1400m工況下重載貨車踏面制動(dòng)時(shí)車輪溫度場(chǎng)與應(yīng)力場(chǎng)進(jìn)行分析。設(shè)定車輪的初始整體溫度為20℃。

1.1 熱流密度計(jì)算

在本文中，采用能量轉(zhuǎn)換法。假設(shè)列車制動(dòng)過(guò)程中重載貨車的動(dòng)能全部轉(zhuǎn)化為熱能，忽略輪軌摩擦熱量輸入，該熱量全部被摩擦的閘瓦和踏面吸收。傳入到踏面的熱量作為熱流密度來(lái)處理，即在與閘瓦摩擦的踏面上形成一個(gè)周向移動(dòng)的面熱源。則每輛車在制動(dòng)過(guò)程中閘瓦和車輪踏面間摩擦產(chǎn)生的熱量Q（ t）為：

2 車輪模型的建立

以直徑840mm的 HESA形整體輾鋼車輪為研究對(duì)象，假設(shè)熱量在整個(gè)與閘瓦摩擦的車輪踏面周向上均勻輸入，故車輪受到的載荷是軸對(duì)稱的，可建立 1/2 車輪模型，進(jìn)行簡(jiǎn)化計(jì)算。之后采用 8 節(jié)點(diǎn)六面體單元（ Solid70） 對(duì)其進(jìn)行網(wǎng)格劃分，最終得到的有限元模型單元數(shù)量為122430個(gè)，節(jié)點(diǎn)數(shù)量為138080個(gè)。

3 車輪溫度場(chǎng)仿真計(jì)算

仿真的總時(shí)間為265s，其中緊急制動(dòng)時(shí)間為65s，制動(dòng)后冷卻散熱時(shí)間為200s。 從曲線可以看出，制動(dòng)初始時(shí)，與閘瓦接觸的踏面溫度急劇上升，并達(dá)到整個(gè)仿真過(guò)程的最高溫度，其他位置，包括踏面兩側(cè)及踏面內(nèi)部的溫度在摩擦面之后達(dá)到峰值；之后隨著制動(dòng)的結(jié)束，各部分溫度開(kāi)始逐漸下降，由溫度云圖可以看出，最高溫度逐漸向輪輞內(nèi)部偏移。

分析溫度隨時(shí)間變化曲線，可見(jiàn)圖（a）中▲線與圖（b）中◆線溫度相同，此處為踏面表面與閘瓦接觸部分的正中位置，距離踏面外側(cè)57mm。一旦制動(dòng)開(kāi)始，則此處溫度急劇上升，在30秒左右達(dá)到最高溫度160.528℃，對(duì)應(yīng)溫度云圖中踏面紅色的區(qū)域。而（a）圖中的●線和╳線對(duì)應(yīng)溫度云圖中踏面橙色的區(qū)域，由于他們位于（a）中▲線的兩側(cè)，溫度的峰值比最高溫度略低；因?yàn)槎吲c中心線距離相同，所以溫度曲線大致重合。（a）中◆線和■線分布在踏面兩側(cè)，故溫度峰值更低，由溫度云圖可見(jiàn)，制動(dòng)結(jié)束后最高溫度逐漸向輪輞內(nèi)部偏移，故距踏面外側(cè)近的◆線較晚達(dá)到溫度峰值。分析（b）圖中的各線，因最高溫度逐漸向車輪內(nèi)部偏移，所以隨著深度的增加，各位置的溫度峰值出現(xiàn)的時(shí)間逐漸延后。

在40～60秒，也就是制動(dòng)過(guò)程趨近結(jié)束時(shí)，車輪各部分大致達(dá)到溫度峰值。隨著制動(dòng)的結(jié)束，不同的位置峰值過(guò)后，溫度大多有一個(gè)急劇下降的過(guò)程，最終在接近200秒時(shí)整體車輪溫度大致相同，為50℃左右，但摩擦面仍保持整體最高溫度。

4 車輪應(yīng)力場(chǎng)仿真計(jì)算

通過(guò)間接耦合法， 將上述溫度載荷施加到重載貨車車輪的模型上，并在對(duì)稱面上施加對(duì)稱約束，與軸配合處施加固定位移約束，建立相應(yīng)的載荷步。從曲線可以看出，在制動(dòng)過(guò)程中車輪踏面的熱應(yīng)力最大，與閘瓦直接接觸的摩擦面出現(xiàn)最大熱應(yīng)力，踏面兩側(cè)以及隨著輪輞深度增加各處的熱應(yīng)力在摩擦面之后達(dá)到峰值；隨著制動(dòng)過(guò)程的結(jié)束，車輪整體熱應(yīng)力迅速減小，由溫度云圖可看出，冷卻過(guò)程中最大熱應(yīng)力逐漸出現(xiàn)在s型輻板的兩個(gè)凹槽處。

分析應(yīng)力隨時(shí)間變化曲線，可見(jiàn)圖（a）中■線與圖（b）中■線應(yīng)力相同，此處為踏面表面與閘瓦接觸部分的正中位置。一旦制動(dòng)開(kāi)始，則此處溫度急劇上升，導(dǎo)致該處的熱應(yīng)力也急劇上升，迅速達(dá)到最高應(yīng)力，時(shí)間是30秒，為193.529Mpa，對(duì)應(yīng)應(yīng)力云圖中踏面紅色的區(qū)域。（a）圖中的▲線和◆線由于高溫的影響也能達(dá)到較高的熱應(yīng)力，但由于其分布在摩擦面兩側(cè)，所以距最高熱應(yīng)力值仍有差距。而 （a）中╳線和●線離摩擦面更遠(yuǎn)，故應(yīng)力峰值更低，距踏面內(nèi)側(cè)較近的●線甚至差不多無(wú)峰值，而且整條線有抖動(dòng)的趨勢(shì)，說(shuō)明應(yīng)力的變化并不平緩。對(duì)于（b）圖中的各線，由于最高溫度逐漸向車輪內(nèi)部偏移，所以熱應(yīng)力峰值的分布也逐漸向輪輞內(nèi)移動(dòng)，隨著深度的增加，各位置的熱應(yīng)力峰值出現(xiàn)的時(shí)間延后，但對(duì)于●線，由于距踏面最深，應(yīng)力變化并不平緩且有凹槽，應(yīng)力的峰值較早較小。

最終在接近200秒時(shí)整體車輪應(yīng)力大致相同，為40MPa左右，但整體最高熱應(yīng)力卻由摩擦面轉(zhuǎn)移到了s型輻板的凹槽處，可見(jiàn)此處為強(qiáng)度薄弱處。

5 結(jié)論

（1）車輪的最高溫度值出現(xiàn)在30s，為160.528℃，位于閘瓦與踏面的摩擦面處。

（2）車輪的最大熱應(yīng)力同樣出現(xiàn)在30秒，位于閘瓦與踏面的摩擦面處，為193.529Mpa，小于該溫度下材料的屈服強(qiáng)度。

（3）制動(dòng)開(kāi)始后，踏面溫度和熱應(yīng)力急劇上升，最高溫度和最大熱應(yīng)力同時(shí)出現(xiàn)在30秒左右的閘瓦與踏面的摩擦面處，可見(jiàn)應(yīng)重視踏面的疲勞損傷。之后隨著制動(dòng)結(jié)束，車輪整體溫度和熱應(yīng)力下降，最高溫度向踏面內(nèi)側(cè)偏移，最高熱應(yīng)力出現(xiàn)在輻板上，但其值都較小。

參考文獻(xiàn)：

[1]李芾，安琪，劉俊紅等.重載貨車車輪溫度場(chǎng)與應(yīng)力場(chǎng)研究[J]. 系統(tǒng)仿真學(xué)報(bào)，2010，22（02）：344-347，369.

[2]宋志坤，謝基龍，張勵(lì)忠.鑄鋼輪材料在700℃～20℃熱循環(huán)下斷裂機(jī)制分析[J].鐵道學(xué)報(bào)，2007，29（02）：141-144.

[3]趙陽(yáng)春.貨車運(yùn)用常見(jiàn)故障分析與處理[M].北京：中國(guó)鐵道出版社，2000.

[4]藍(lán)春紅，吳萌嶺，王勇等.重載貨車踏面制動(dòng)時(shí)車輪溫度場(chǎng)與應(yīng)力場(chǎng)研究[J]. 鐵道車輛，2011，49（07）：1-5.

[5]王京波.合成閘瓦對(duì)車輪熱影響的研究[J].鐵道機(jī)車車輛，2003，23（02）：74-79.

[6]劉云.提速貨車車輪溫度場(chǎng)及熱應(yīng)力場(chǎng)的數(shù)值模擬[D].北京：北京交通大學(xué)，2004.

[7]侯耐.重載貨車車輪踏面制動(dòng)熱負(fù)荷研究[D].成都：西南交通大學(xué)，2011.

[8][美]R.西格爾，J.R.豪厄爾.熱輻射傳熱[M].北京：科學(xué)出版社，1990.

項(xiàng)目：本文由國(guó)家大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計(jì)劃項(xiàng)目資助。