列車制動盤表面溫度分布演變的研究

農萬華，熊啟純，韓曉明，禤釗堅

列車制動盤表面溫度分布演變的研究

農萬華1，熊啟純1，韓曉明2，禤釗堅3

（1.柳州工學院 機械工程系，廣西 柳州 545616；2.大連交通大學 連續擠壓教育部工程研究中心，遼寧 大連 116028；3.東風柳州汽車有限公司，廣西 柳州 545616）

高速列車制動條件下制動盤的溫度不均勻分布是影響盤熱疲勞損傷的一個重要因素。本文針對銅基粉末冶金閘片與鍛鋼制動盤構成的摩擦副，利用1:1制動動力試驗臺，進行了最高速度達350km/h的制動試驗，結合紅外熱像與有限元分析技術，研究了制動速度和制動壓力對盤面溫度分布形態的影響。制動實驗結果顯示，在制動過程中盤面高溫區優先在摩擦中徑兩側形成，并隨制動時間的增加向摩擦中徑處偏移。有限元分析結果顯示，制動盤摩擦中徑兩側的接觸壓力明顯偏大。由此可推斷，閘片的燕尾結構致使摩擦中徑兩側壓力偏高，是盤面形成兩個高溫區的主要原因。摩擦速度隨摩擦半徑的增加而增加，導致摩擦半徑大于摩擦中徑的高溫區優先磨損，磨損的不均勻性改變了接觸壓力分布。摩擦速度和接觸壓力的綜合作用與高溫區演化規律密切相關。

制動盤；閘片結構；摩擦磨損；摩擦溫度

高速列車的基礎制動方式均為盤式摩擦制動，制動盤與閘片通過摩擦將列車的動能轉成熱能并加以吸收儲存，因此每次停車制動盤都要經歷一次熱負荷的作用。列車制動時，特別是高能制動，制動盤的溫度和溫度梯度均處于相當惡劣狀態[1]，容易致使制動盤熱疲勞損傷而失效。因此，熱疲勞損傷成為決定制動盤的壽命的主要因素之一。

為認識影響制動盤熱疲勞性能的基本規律，人們采用有限元法模擬制動過程，分析制動盤的溫度和熱應力分布規律，利用熱電偶[2]或紅外熱像儀[3-4]的測量制動盤的局部溫度變化或表面的溫度分布情況，研究了制動壓力[5-6]、制動盤材料[7-9]、制動模式[10-11]，閘片材料[12-13]等對制動盤的溫度影響，這些工作均為掌握盤表面的熱疲勞損傷機理奠定了良好的基礎。然而，列車制動盤溫度分布是由制動初速度、摩擦副接觸壓力、閘片結構和摩擦磨損等多因素綜合作用的結果，分析制動摩擦表面溫度分布與多種影響因素的綜合關系，對解決制動盤的熱疲勞損傷是有意義的。

本文利用高速列車制動盤與銅基粉末冶金閘片摩擦副，在1:1制動動力試驗臺上進行最高速度達350 km/h的制動試驗，通過紅外熱像技術檢測制動盤摩擦表面溫度，并從制動速度、制動壓力、閘片結構和摩擦磨損等方面，對影響制動盤表面溫度分布的因素進行分析。

1 材料與試驗方法

采用銅基粉末冶金閘片和鍛鋼制動盤組成摩擦副，在1:1制動動力制動試驗臺進行了120次制動測試。試驗參數如表1所示。每次制動前制動盤都進行充分冷卻，盤面溫度小于60℃后才開始。采用紅外熱像儀（flukeTi45）記錄制動盤面溫度，熱像儀的發射率取0.75[12]，熱像儀記錄視野窗口大小如圖1所示。

表1 試驗參數

圖1 熱像儀觀察位置示意圖

2 實驗結果及分析

2.1 速度、壓力對盤面溫度分布的影響

用熱像儀所記錄制動過程中不同時刻制動盤摩擦表面的紅外圖像，然后在制動盤相同周向位置，沿著徑向提取出制動盤表面的溫度值（圖2）。

圖2 盤面溫度采集位置示意圖

圖3是制動力為11 kN，制動初速度分別為120、200和350 km/h時，制動過程中盤摩擦表面溫度徑向分布情況。可見制動開始時，三種情況都分別在制動盤半徑為約300 mm和180 mm的兩個區域形成兩個高溫區，摩擦中徑處的溫度低于其兩側的溫度。在制動過程中，制動初速度為120 km/h（圖3a）時，高溫區的位置及寬度都沒有明顯變化，制動初速度升高（圖3b、c），高溫區的位置逐漸向摩擦中徑位置移動，寬度增大，但摩擦中徑處的溫度仍低于其兩側的溫度。這表明，制動初速度的增加對盤面高溫區形成初始位置影響不明顯，但對高溫區的分布演變有一定的影響。

圖3 不同速度下盤面溫度分布徑向演變情況（制動力為11 kN）

圖4為制動初速度250 km/h，制動力分別為5 kN、7.5 kN和11 kN時，制動過程中盤摩擦表面溫度徑向分布情況?？梢姡苿恿? kN時（圖4a），制動初始形成了一個位于盤半徑約為300 mm的高溫區，隨制動進行，溫度快速升高后逐漸降低，并向摩擦中徑方向有所偏移，制動時間超過60 s后，這個高溫區基本消失。而另一個處于盤半徑為170 mm的高溫區，隨制動時間的增加，向摩擦中徑方向偏移明顯，最后發展到盤半徑約為230 mm處。當制動壓力升高時（圖4b、c），盤面溫度的變化規律相似。這表明，制動壓力的增加對盤面高溫區形成初始位置及其分布演變的影響不明顯。

圖4 不同制動載荷下盤面溫度分布徑向演變情況（制動初速度為250 km/h）

此外，在制動過程中也會出現僅形成一個高溫區的情況，如圖5所示。一種是在盤半徑約為270 mm（圖5a）或230 mm（圖5b）處出現一個較明顯高溫區，在制動過程中，高溫區中心位置及其寬度未有明顯變化。另一種情況是在盤半徑約為250 mm處形成高溫區，高溫區隨制動的進行有逐漸發展為兩個的趨勢。

圖5 盤面溫度分布徑向演變

2.2 制動盤表面溫度分布及演變的影響因素分析

（1）高溫區初始形成位置影響因素分析

從實驗結果看，盤面高溫區的初始位置基本相同，制動力、制動初速度對高溫區初始形成位置的影響不明顯。圖7是試驗中采用的摩擦副及其連接件的軸向剖面圖，從摩擦副的結構看，閘片與閘片托采用燕尾方式連接，在燕尾槽處有安裝間隙。制動力通過閘片托從閘片燕尾的兩側作用到閘片背鋼板上，然后通過摩擦塊作用在摩擦副接觸面。

根據盤、閘片和閘片托三者的關系建立其3D模型（圖8），采用有限元方法計算制動盤表面與閘片的靜態接接觸力。由于閘片托厚度、剛度相對于閘片背鋼板大很多，在制動力的作用下發生變形可忽略不計，計算中假設閘片托為剛體，施加在閘片托的制動壓力均勻分布，制動盤與盤轂接觸一側的自由度個數為0。圖9是盤面接觸壓力分布分析結果，可見（基于閘片結構，沿徑向將制動盤接觸面分為三個部分：燕尾外側區Ⅰ，燕尾區Ⅱ和燕尾內側區Ⅲ），三個區域相比，燕尾區Ⅱ的接觸力分布處于最低值，其它兩區的接觸力相對較高，且分布位置與高溫區初始位置基本相同。由此可見，閘片結構形式決定了摩擦副接觸力分布，是影響制動過程中盤面高溫區優先形成位置的重要原因。

圖7 摩擦副軸向截面圖

圖8 有限元分析模型

圖9 盤面接觸壓力分布云圖

（2）制動盤的表面溫度分布演變影響因素分析

實驗結果表明，每次制動時隨著時間增加，盤面高溫區逐漸從初始位置向盤面中間位置移動，制動力、制動初速度增大，高溫區偏離初始位置越明顯，說明制動力、制動初速度對高溫區的移動方向無明顯影響，但對高溫區的移動快慢有明顯影響。盤面溫度是摩擦副相互作用的結果，與摩擦副相對運動、接觸力、摩擦副的接觸面磨損等有關。

根據Archard摩擦磨損模型[15,20]，磨損厚度與摩擦速度和接觸力的關系為：

由式（2）可知，表面磨損速度與接觸力和摩擦速度成正比，圖9表明，閘片燕尾兩側的接觸壓力高于燕尾區，顯然，從接觸力分布的角度，閘片兩側區域的磨損速度應該大于中部區域。另一方面，摩擦速度隨制動盤半徑增大而增大，因此閘片燕尾外側區的摩擦速度最快，磨損速度最大。隨著溫度的上升摩擦材料的硬度下降[21]，其磨損速度也會加快。綜合接觸力和速度的分布及溫度的變化情況，可知閘片燕尾外側區應該是磨損速度最高的區域。

閘片和閘片托的結構關系表明，若閘片內外兩側受力不平衡時，則閘片托具有繞連接銷孔中心轉動的趨勢。假設閘片出現轉動前后的幾何關系如圖10所示，軸表示制動盤的接觸面，制動力作用線與燕尾的中心線一致，與軸交與點，假設制動力方向始終與閘片鋼背平面法向保持一致。設制動盤摩擦中徑兩側接觸力中心分別為1、2，則根據力的平衡關系可知，在制動初始兩點間存在以下關系：

式中：1、2分別為制動盤摩擦中徑兩側接觸力中心的作用力，MPa；1、2分別為兩側接觸力中心到制動力作用線的距離，m。

隨著制動的進行，閘片燕尾外側磨損較快致使閘片發生傾斜，由于制動力方向始終與閘片鋼背平面法向保持一致，故制動力方向偏轉角。偏轉后制動力的作用線與軸交于點，如圖10（b）所示。根據力平衡原理可知，1、2兩點處的合力方向也偏轉角，則1、2點到制動力作用線的距離分別為：

故可得磨損后1、2兩點的力平衡關系為：

顯然，由式（4）、式（5）可知1＜1、2＜2，所以比較式（3）、式（6）可以得到1＞1、2＞2。這說明由于磨損導致燕尾外側區接觸力減小，摩擦功率密度降低，從而使外側高溫區在制動過程中優先消失。燕尾內側區盡管接觸力提高會增加磨損速度，但由于摩擦速度較外側低，磨損程度的增加不明顯，因此制動盤內側的高溫區具有較長的保持時間，如圖3、圖4所示。

圖10 摩擦副接觸壓力的演變示意圖

隨著這兩個區域磨損程度增加，造成閘片中部區域相對高度增加，接觸力增加，摩擦功率密度升高，從而使高溫區隨制動過程向制動盤摩擦中徑方向移動（圖3、圖4），因此摩擦副的磨損與高溫區移動密切相關[15-19]，是影響高溫區移動的主要因素。

3 結論

（1）在制動過程中，優先在制動盤摩擦中徑兩側出現兩個高溫區，而且隨著時間的增加高溫區逐漸向制動盤摩擦中徑移動。這說明制動初始，摩擦功率高密度區分布在制動盤摩擦中徑兩側區域，隨后摩擦功率高密度區向其它區域移動。

（2）有限元計算結果表明，由于閘片與閘片托采用的燕尾結構連接，盤面接觸力分布呈中間低兩側高的非均勻分布狀態，是造成高溫區優先出現在制動盤摩擦中徑兩側的主要原因。

（3）制動盤的外側處于高摩擦速度，內側則具有高的接觸力，使得各區域摩擦功率的不同，閘片接觸面各區域磨損快慢有差別，從而導致高溫區在制動過程中發生移動現象。

[1]Benseddiq N.， Seiderrnann J.. Optimization of Design of Railway Disc Brake Pads[J]. Proc Instn Mech Engrs，1996（210）：51-61.

[2]莊光山，王成國，王海慶，等. 盤形制動摩擦表面溫升研究[J].機械工程學報，2003，39（2）：150-154.

[3]Stéphane Panier，Philippe Dufrénoy，Pierre Brémond. Infrared characterization of thermal gradients on disc brakes[J]. Proceedings of SPIE，2003：295-302.

[4]P Litos，M Honner，V Lang，et al. A Measuring System for Experimental Research on the Thermomechanical Coupling of Disc Brakes[J]. Proc. IMechE Part D：J. Automobile Engineering，2008：1247-1257.

[5]趙文清. 高速列車“中華之星”制動盤溫度場及熱應力[J]. 兵工學報，2006，27（7）：132-136.

[6]Sung-Keun Cho，Jung-Hun Choi. Life Evaluation of a Disk Brake of Railway Vehicles considering Pressure Distributions at a Frictional Surface [J]. Key Engineering Materials，2007（353-358）：303-306.

[7]陳德玲，張建武，周平. 高速輪軌列車制動盤熱應力有限元研究[J]. 鐵道學報，2006，28（2）：39-43.

[8]劉金朝，卜華娜，劉敬輝，等. 整體制動盤熱應力有限元仿真分析[J]. 中國鐵道科學，2007，28（2）：80-83.

[9]M. Eltoukhy，S. Asfour et al. Thermoelastic Instability in Disk Brakes：Simulation of the Heat Generation Problem[C]. Boston：Excerpt from the Proceedings of the COMSOL Users Conference，2006.

[10]趙海燕，張海泉，湯曉華，等. 快速列車盤型制動熱過程有限元分析[J]. 清華大學學報，2005，45（5）：589-592.

[11]周凡華，吳光強，沈浩，等. 盤式制動器15次循環制動溫度計算[J]. 汽車工程，2001，23（6）：411-418.

[12]Gérard Degallaix，Philippe Dufrénoy，Jonathan Wong，et al. Failure Mechanisms of TGV Brake Disc [J]. Key Engineering Materials，2007（345）：697-700.

[13]朱愛強，劉佐民. 材料匹配性對盤式制動器摩擦溫度場的影響[J]. 潤滑與密封，2008，33（3）：63-65.

[14]Haruo SAKAMOTO，Kenji HIRAKWA. Fracture Analysis and Material Improvement of Brake Discs [J]. JSME International Journal，2005，48（4）：458-464.

[15]V. Hegadekatte，N. Huberb，O. Kraft, Modeling and Simulation of Wear in a Pin on Disc Tribometer[J]. Tribology Letters，2006，24（1）：51-60.

[16]K. Lee，J. R. Barber. An Experimental Investigation of Frictionally-Excited Instability in Automotive Disk Brakes under a Drag Brake Application[J]. ASME J. Tribol.，1994（116）：409-414.

[17]J.-J. Thomas，A. Bignonnet，T.M.-L. Nguyen-Tajan. Reliable Design of Automotive Brake Discs[C]. Lille：Proceedings of the European Conference on Braking JEF，2002：79-86.

[18]P. Dufrénoy，D. Weichert. Prediction of Railway Disc Brake Temperatures Taking the Bearing Surface Variations Into Account[J]. Proc. Instn Mech. Engrs, Part F，1995（209）：67-76.

[19]Dufrénoy，D. Weichert. A Thermomechanical Model for the Analysis of Disc Brake Fracture Mechanisms [J]. Journal of Thermal Stresses，2003（26）：815-828.

[20]S?ren Andersson，Wear Simulation. Advanced Knowledge Application in Practice. Igor Fuerstner (Ed.), ISBN: 978-953-307-141-1, InTech, Available from: http://www.intechopen.com/books/advanced-knowledge-application-in-practice/wear-simulation，2010：15-36.

[21]孫洪雨，馬元，陳輝，等. 高速列車制動材料高溫摩擦磨損行為研究[J]. 機械，2018，10（45）：5-10.

Research on the Evolution of Surface Temperature Distribution of Train Brake Disc

NONG Wanhua1，XIONG Qichun1，HAN Xiaoming2，XUAN Zhaojian3

(1.Department of Mechanical Engineering, Liuzhou Institute of Technology, Liuzhou 545616, China;2.Engineering Research Center of Continuous Extrusion, Ministry of Education, Dalian Jiaotong University, Dalian 116028, China; 3.Dongfeng Liuzhou Motor Co., Ltd., Liuzhou 545616, China)

The non-uniform distribution of the temperature of the brake disc during the braking of high-speed train is an important factor influencing the thermal fatigue of the brake disc. In this paper, a friction pair of copper-based powder metallurgy brake pad and steel disc is tested on a full-scale test bench at a maximum braking speed of 350 km/h.With the technology of infrared thermal imaging and finite element analysis, the influence of the braking speed and braking pressure on the evolution law of temperature distribution and influencing factors is investigated. The experimental results show that the high temperature zones first form at both sides of the median friction diameter of the disc surface, and shift to the median diameter with the increase of the braking time. The finite element analysis results show that the contact pressure on both sides of the median friction diameter of the brake disc is obviously at a higher level, which indicates that the high pressure on both sides of the median friction diameter caused by the dovetail structure of the brake pad leads to is the main reason for the formation of the two high temperate zones on the disk surface.. The friction velocity increases with the increase of the friction radius, which leads to the preferential wear in the high temperature zone where the friction radius is larger than the friction mid-diameter, and the non-uniform wear changes the contact pressure distribution. The combined effect of friction velocity and contact pressure is closely related to the evolution law of the high temperature zone.

brake disc；brake pad structure；frictional wear；friction temperature

U26O.351

A

10.3969/j.issn.1006-0316.2021.08.007

1006-0316 (2021) 08-0046-07

2020-12-18

廣西高校中青年教師基礎能力提升項目（2019KY1099）

農萬華（1980－），男，廣西貴港人，博士，講師，主要研究方向為軌道車輛閘片材料與結構、機械摩擦磨損，E-mail：328868608@qq.com。