40t軸重貨車制動熱負荷分析

陳 倩，李 芾，王軍平，張麗霞

（西南交通大學 機械工程學院，四川成都610031）

提高鐵路的運輸能力是改善鐵路工作運營效率，降低運輸成本的重要手段，其主要通過如下3種方法來實現：加大軸重、增加行車密度和提高行車速度。目前，我國鐵路主要干線上的行車密度已居于世界前列，而行車速度的提高又受到多重因素的限制，難度較大，潛力有限，相對而言，增加軸重對運能的提高意義更大，因此，發展重載運輸是提高運能的一條重要途徑［1］。

鐵路重載技術具有軸重大、載重大、自重輕、列車牽引噸位大的特點［2］。自20世紀60年代以來，發展大軸重運輸已成為重載運輸的一個主要發展方向。目前，較為成熟的國家有：美國（35.7～39t）、加拿大（35.7～39t）、澳大利亞（40t）、巴西（30t）、俄羅斯（27t）和南非（30t）等。我國鐵路貨車技術的自主發展已有50多年的歷史，自主設計的重載貨車已經出口到澳大利亞、巴西和委內瑞拉等國。與世界重載運輸大國相比，中國鐵路近年來軸重和載重雖然都有所提高，但仍存在一定差距。根據我國鐵路的中長期發展需求，在已有的大軸重貨車設計、制造的成功經驗基礎上，研制適應于中國鐵路的40t軸重貨車具有十分重要意義。

在鐵路運輸中，安全性能始終是最受關注的問題，在眾多影響列車運行安全的因素中，制動問題至關重要。為降低運營成本，貨車一般采用踏面制動，軸重的增加使得列車的動能增加，這就意味著制動時閘瓦產生的摩擦熱增加，車輪的溫度和熱應力也將隨之增加。一旦車輪的熱負荷惡化，就會形成熱裂紋和熱失效，踏面也會由于大量制動熱輸入而產生剝離，直接影響到車輪的疲勞壽命，對鐵路運輸的安全產生不利影響。

以我國重載專運線大秦鐵路為背景，分別對列車運行時的緊急制動和長大下坡道制動進行了模擬分析。其中，每一制動工況都包括加熱和冷卻兩個階段。根據常規試驗方法，緊急制動是在平坡道上進行。長大下坡道則選取大秦線重車方向為50km的長大下坡道（平均坡度為9.1‰）。根據實際運行經驗和相關規定，兩種制動的初速度（vo）選取為100，80km／h。

車輛在執行制動的過程中，熱源主要考慮來自閘瓦與踏面摩擦產生的熱量。由于車輪在制動時處于高速轉動狀態且屬于軸對稱結構，所以在計算踏面上熱量輸入和耗散時認為熱負荷是軸對稱且在周向上是沒有變化的。

1 模型建立

圖1 車輪有限元模型

以某40t軸重貨車為例，利用ANSYS有限元分析軟件建立了直徑為915mm的車輪1／2模型，如圖1所示。熱分析時，模型采用三維溫度單元solid70，結構分析時將其轉化為相應的三維結構單元solid45。輪軸采用過盈配合（過盈量為0.3mm），兩者接觸面通過TARGE170單元和CONTA174單元建立三維面—面接觸。閘瓦和車輪踏面接觸區域的實體單元承受熱量輸入，車輪外部（除輪轂輪孔和車輪截面外）建立表面效應單元SURF152，承受熱量輸出，同時考慮了熱輻射。

2 邊界條件確定

主要對緊急制動和長大下坡道制動兩種工況進行計算。文中緊急制動工況主要模擬列車在緊急狀況下采取一次制動，速度從最大衰減到零的情況。以我國目前的制動配置和2萬t重載列車的運用情況來看，其性能能夠滿足緊急制動1 400m停車的運用要求［2］。長大下坡道制動工況主要以我國大秦鐵路為背景，選取了其中平均坡度為9.1‰的線路區間（K275～K325），此線路長度為50km［3］。長大下坡道制動工況主要模擬列車恒速制動的情況，即控制列車坡道運行不增速，制動力和列車坡道下滑力達到平衡狀態。

制動過程中，閘瓦和踏面間的摩擦熱經踏面傳給車輪，車輪上熱量通過對流和輻射等方式耗散到周圍空氣中。本文所采用的工況參數如表1所示。

表1 制動參數

熱量輸入通過熱流密度表示，采用摩擦功率法［7］可以得到熱流密度。其中，閘瓦和踏面接觸區域的摩擦功率W如下：

式中η為輸入到車輪的熱流密度分配系數，取0.91［3］；F為滑動摩擦力，kN；N為閘瓦壓力，kN；μ為閘瓦和踏面間的摩擦系數；v（t）為車輪和閘瓦的相對滑動速度，km／h，t為制動時間。

假設車輪旋轉一周熱載荷均勻分布，則熱流密度q（t）如下：

式中Sf為閘瓦在踏面上掃過的面積，m2；D為車輪直徑，m；L為閘瓦寬度，m。

緊急制動速度變化急劇，對流換熱系數α與速度關系較大，沿用經驗公式［3-7］：

長大下坡道制動可看成圓柱擾流的對流換熱，則該對流換熱系數與車輪表面空氣的導熱系數、運動黏度、普朗特數、雷諾數、努謝爾特數相關，計算公式如下：

式中，λ為空氣導熱系數，W／（m·K）；l為固體壁面特征尺寸（本文取車輪圓周的一半），m；Nu為努謝爾特數。

式中，Pr為普朗特數，且Pr=0.687；Re為雷諾數。

雷諾數由下式決定：

其中，u∞為空氣流動速度，m／s；γ為空氣的運動黏度，m2／s。

在ANSYS有限元軟件中，輻射換熱主要通過輻射率來體現，本文輻射系數取ε=0.66［3］。

3 計算結果分析

3.1 溫度場

將邊界條件施加到有限元模型上，進行瞬態熱分析，得到40t軸重車輪制動時的最高溫度云圖如圖2～圖3。

圖2 緊急制動最高溫度

圖3 長大下坡道制動最高溫度

根據圖2、圖3可知，緊急制動和長大下坡道制動出現最大溫度時間分別在第36s和第2 250s，其最大溫度分別為233℃和231℃。前者出現最高溫度的地方為閘瓦與踏面接觸的中心，而后者則出現在偏離摩擦面中心線的位置。

緊急制動時踏面的最高溫度先急劇上升，然后緩慢下降，并在制動的第36s時達到最大。由于緊急制動開始的0～36s間，熱量在踏面上積聚，最高溫度曲線上升急劇，熱流密度隨著速度的降低而降低，36s后熱量輸入小于熱量輸出但并不為零，所以溫度梯度開始緩慢下降，且向輪輞內部轉移，在第82s時制動結束，此后的冷卻時間熱量輸入停止，熱量輸出通過自然對流的形式進行（緊急制動最高溫度時間歷程如圖4所示）。

圖4 緊急制動最高溫度時間歷程

長大下坡道制動踏面的最高溫度出現在制動的第2 250s，即制動結束時刻。從制動開始到結束，熱流密度和換熱系數都為恒定值，溫度一直處于上升的趨勢。0～750s期間，曲線上升比較急劇，之后曲線上升趨勢逐漸放緩，這是由于制動初期熱量迅速在踏面積聚，隨著制動時間的推移，熱量逐漸向輻板區域傳遞的緣故。制動結束后熱量停止輸入，熱量輸出以自然對流形式散熱，所以溫度迅速下降，表現為曲線中的尖角（長大下坡道制動最高溫度時間歷程如圖5所示）。

圖5 長大下坡道最高溫度時間歷程

3.2 應力場

進行應力場分析時，在軸端施加全位移約束，對稱面施加對稱約束，將溫度場作為載荷施加到各個節點上。分析得到的最高應力場云圖如圖6～圖7。

圖6 緊急制動最高熱應力

圖7 長大下坡道制動最高熱應力

根據圖6、圖7可知，緊急制動和長大下坡道制動最大熱應力分別出現在踏面和輻板外側，其值分別為348MPa和252MPa。

在緊急制動中，由于制動時間較短，產生的熱量在踏面區域積聚后無法及時向其他部位擴散，所以最大熱應力出現在踏面上，此后溫度慢慢向輪輞內部傳遞，溫度梯度減小，踏面熱應力逐漸減小，輻板熱應力增加。（緊急制動最高熱應力如圖6所示）。

在長大下坡道制動中，熱量輸入較大且持續時間較長，制動初期，踏面上的熱量積聚，溫度最大，熱應力也最大，雖然踏面一直有熱量輸入，但隨著時間的增加，熱量開始向輪輞內部轉移，最大熱應力由踏面逐漸轉移到了輻板外側［5-7］（長大下坡道制動最高熱應力如圖7所示）。

圖8、圖9分別為緊急制動和長大下坡道制動時最高熱應力時間歷程，將其與圖4、圖5進行比較可以發現，車輪最高熱應力的變化趨勢與相應的溫度變化趨勢基本一致，即溫度梯度越高，對應時刻的熱應力也就越大。同時，對不同溫度下的von Mises應力進行計算分析可得，最大應力值均未超出材料的屈服極限。圖9曲線上所出現的小波動，則是由最大熱應力位置的變化所致。

圖8 緊急制動最高熱應力時間歷程

圖9 長大下坡道制動最高熱應力時間歷程

本文考慮了溫度對材料力學性能的影響，特別是對膨脹系數和屈服應力的影響。制動過程中閘瓦與踏面接觸區域的材料因熱量輸入而迅速膨脹，由于受到周圍冷金屬基體的約束，踏面上會形成壓縮應力［6］。隨著向輪輞深度溫度梯度的變小，壓應力逐漸減小，甚至會變為拉應力，拉應力太大時，會使車輪發生塑性變形。故需要繼續對車輪的應力水平進行分析。

圖10為緊急制動時，踏面上最大熱應力部位的徑向、周向、軸向應力變化圖。根據圖10可知，踏面上熱應力的徑向分力為較小的拉應力且其值變化不大；軸向應力在制動結束前均為壓應力，其值先增大后減小，制動結束后，壓應力轉變為較小的拉應力；周向應力在制動期間和制動結束后均為壓應力。對比圖8和圖10發現，緊急制動時，最高熱應力為正值，周向應力為負值，但其最高熱應力跟周向應力變化規律一致。圖11為長大下坡道制動時，車輪輻板外側最大熱應力部位的三向應力變化圖。由圖可知，其徑向分力始終為較小的拉應力；在制動結束后，不再有熱量輸入，車輪受熱層開始冷卻，所以周向壓應力開始迅速減小，但并未變成拉應力；軸向壓應力轉變成了較小的拉應力。

圖10 緊急制動三向應力

圖11 長大下坡道制動三向應力

以上是對緊急制動工況和長大下坡道制動工況下車輪的各向應力水平的分析。從仿真結果來看，車輪的軸向壓應力在制動結束后雖都變成拉應力，但其應力水平較低，故該車輪發生塑性變形幾率較小。但是在實際運用過程中，列車制動的情況較為復雜，車輪不僅受到熱負荷的影響還受其他載荷的影響，如機械載荷。在熱—機械載荷耦合的情況下，車輪的應力變化情況將更加復雜，更精確的結果還有待進一步的深入研究。

4 結論

通過對軸重40t貨車車輪制動熱的分析，可以得出以下結論：

（1）兩種制動最高溫度均出現在踏面處，其幅值分別為233℃和231℃。緊急制動的最高溫度出現在第36s，此時熱量輸入和輸出達到平衡狀態。長大下坡道制動時列車速度恒定，持續時間長，輻板區域的溫度變化較明顯，最高溫度出現在制動結束時刻。

（2）兩種制動的最高熱應力變化趨勢與各自的最高溫度變化趨勢一致，出現最大熱應力的部位分別為踏面和輻板外側，其幅值分別為348MPa和252MPa。

（3）對不同溫度下的von Mises應力進行計算分析可得，最大應力值均未超出材料的屈服極限。（本文僅考慮熱應力）。

（4）考慮溫度對材料力學性能的影響，兩種制動的周向和軸向分力在制動結束前都為壓應力，制動結束后軸向分力轉變為拉應力，但其應力水平較低；徑向分力都為拉應力，但其應力水平也較低。所以，僅從熱負荷方面考慮，該車輪發生塑性變形的幾率較小。

本文僅從制動熱負荷方面探索了將貨車軸重增加到40t時的可行性。貨車軸重的增加勢必會使輪軌間作用力增加，制動距離增加，而且車輛的動力學性能也會隨著改變，同時還會對線路有很大的影響。所以還需對其動力學等方面進行深入研究。

［1］劉向東.淺析我國鐵路的重載運輸［J］.科技創新導報，2010，（18）：1.

［2］劉化宇，楊知猛，劉 華.30t軸重專用運煤貨車構想［A］.實踐 開拓 創新--2008年貨車技術發展學術研討會論文匯編［C］，2008：18.

［3］劉俊紅.重載貨車踏面制動熱響應分析研究［D］.成都：西南交通大學，2006，11：30-60.

［4］劉 云.提速貨車車輪溫度場及熱應力場的數值模擬［D］.北京：北京交通大學，2004.

［5］侯 耐.大軸重貨車車輪熱負荷下疲勞強度分析［J］.鐵道機車車輛，2011，31（01）：14-15.

［6］向 陽.30t軸重貨車車輪熱應力分析［J］.鐵道機車車輛，2009，30（06）：30-32.

［7］鄭紅霞，劉玉軍，張全忠，謝基龍.典型制動條件下貨車車輪溫度場和熱應力場的有限元仿真［J］.機械科學與技術，2011，28（04）：501-502.

［8］T Vernersson.Temperatures at railway tread braking Part 1：Modelling［J］.Proceedings of the Institution of Mechanical Engineers Part F：Journal of Rail and Rapid Transit，2007，221：170-175.