機車曲線通過的接觸分析

1、機車曲線通過的接觸分析吳昌華 孫麗萍 張 ?。ù筮B交通大學(xué)） 摘要 本文利用參變量變分原理和由其導(dǎo)出的有限元參數(shù)二次規(guī)劃法，結(jié)合多重多支子結(jié)構(gòu)技術(shù)，對機車通過曲線時的力學(xué)行為進行了研究。文章以東風(fēng)11型內(nèi)燃機車通過R1600的曲線為例，完全模擬輪軌實際接觸關(guān)系建立了計算模型，分析了輪軌各個組成部分的變形與應(yīng)力狀態(tài)以及他們隨著機車的運行速度和軸重等參數(shù)的變化而變化的規(guī)律。關(guān)鍵詞 機車曲線通過 輪軌接觸 有限元 參變量變分原理 二次規(guī)劃法 一 概述列車在線路上運行，或者通過直線，或者通過曲線，二者必居其一。因而列車曲線通過是列車運行的基本載荷工況，必須進行仔細的分析。列車通過曲線時主要的作用載荷是

2、其牽引力、離心力和曲線線路外軌超高引起的向心力。由于列車通過曲線的速度是變的，離心力和向心力，一般說，不會相等，加上牽引力的作用、鋼軌軌頭和車輪踏面的配合關(guān)系以及輪軌間摩擦力等的影響，于是就導(dǎo)致從機車到車輛各個轉(zhuǎn)向架在曲線線路上水平位置產(chǎn)生相對鋼軌不同角度的傾斜，從而也對各根車軸構(gòu)架力的大小產(chǎn)生影響。如果列車的運行速度、各個車輪的輪重和踏面形狀等與線路的曲率半徑、外軌超高以及鋼軌斷面等參數(shù)設(shè)計匹配得不好，列車通過曲線時會引起相當大的輪軌動力作用，降低線路和列車走行部分的強度和壽命，影響列車運行的安全，甚至造成脫軌。從1825年英國建成世界上第一條鐵路以來，鐵路工作者就一直沒有中止對列車曲線通過

3、的研究，但是由于問題太復(fù)雜，計算模型一直很粗糙。半個多世紀以前蘇聯(lián) 教授提出的轉(zhuǎn)心法是當時國際上機車曲線通過計算方法的代表，而且一直沿用至今。但該方法是以剛體力學(xué)作為基礎(chǔ)的，由于受到當時計算手段的歷史限制，不可能計入車輪和鋼軌的彈性變形，更不可能考慮輪軌接觸的非線性，至于扣件的作用則只能完全忽略了，顯然，這些都將對計算結(jié)果帶來比較大的誤差。最近幾年我國鐵路連續(xù)實行了五次大規(guī)模提速，很多干線客運列車的最高運行時速已經(jīng)從120km/h提高到了160km/h，有的試驗列車的時速甚至達到了200km/h以上。但由于多方面的原因提速后線路的變動不大，這就大大加劇了列車通過曲線時的輪軌動力作用。例如，在最

4、早實行提速的廣深線上，有的R600曲線鋼軌上已經(jīng)發(fā)現(xiàn)產(chǎn)生了多處裂紋，直接威脅到列車運行安全，亟需找出原因予以解決。這首先需要進行理論和數(shù)值分析。但這是輪軌多體彈塑性摩擦接觸問題，經(jīng)典的轉(zhuǎn)心法對其是無能為力的，號稱當前國際上在輪軌關(guān)系方面水平最高的Kalker的三維非赫茲彈性體滾動接觸理論，由于無法模擬輪軌的真實幾何形狀和邊界條件，對其也不能進行精細的分析。本文利用參變量變分原理1和由其導(dǎo)出的有限元參數(shù)二次規(guī)劃法，并結(jié)合多重多支子結(jié)構(gòu)技術(shù)，針對必然貼靠外軌的機車第一輪過曲線的情況建立了完全模擬實際的計算模型，對輪對、鋼軌、扣件、橡膠墊板、軌枕和道床等輪軌各個部件的變形和應(yīng)力狀態(tài)隨著機車的速度、軸

5、重等參數(shù)的變化而變化的規(guī)律進行了分析，對機車通過曲線時的輪軌關(guān)系的計算機模擬做了一個初步的嘗試。二 計算理論和算法本文采用中國科學(xué)院鐘萬勰院士提出的參變量變分原理及基于此原理的有限元參數(shù)二次規(guī)劃法來求解輪軌接觸問題1?？臻g彈性接觸力學(xué)邊界問題在區(qū)域（見圖1）上的基本方程為 平衡方程 (1) 應(yīng)變位移關(guān)系 (2) 邊界條件在上滿足在給定力的邊界上 (3) 圖1 兩物體接觸體系 在給定位移的邊界上 (4)在可能接觸的邊界上的邊界條件由接觸系統(tǒng)狀態(tài)方程確定，即 (5)本構(gòu)方程 (6)彈性接觸問題參變量最小勢能原理：在所有滿足幾何條件（2）、（4）的可能位移場中，真實解使總勢能泛函 (7)在接觸狀態(tài)方

6、程的控制下取總體最小值。其中，、分別為應(yīng)力、體積力、應(yīng)變、位移；、是微分算子矩陣；、為給定面力和給定位移；是原約束松弛變量；、分別為彈性矩陣和接觸面彈性矩陣；是常數(shù)矩陣，它表征單元發(fā)生滑動時的接觸彈性松弛力；而為不參加變分的參變量，其物理意義為滑動參數(shù)。將接觸系統(tǒng)離散再對總勢能泛函（7）求極值，可推導(dǎo)出彈性接觸問題有限元參數(shù)二次規(guī)劃方程。 （8）式中為總剛度矩陣；為接觸矩陣；為總約束矩陣；U為總強化矩陣；為總約束向量。標準二次規(guī)劃問題是一個凸規(guī)劃問題，其解法（如Lemke）是成熟的，這就從根本上保證了計算的高效率和高精度。三 計算模型為了真實地模擬輪軌接觸情況，按照機車車輪、60kg/m鋼軌、

7、扣件、橡膠墊板、軌枕和道床等實體三維尺寸建立有限元模型。計算以東風(fēng)11內(nèi)燃機車JM型磨耗型踏面的車輪為例。由于機車通過曲線時，第1輪對的外側(cè)車輪與鋼軌在輪緣處發(fā)生接觸，內(nèi)側(cè)車輪與鋼軌在踏面處發(fā)生接觸，這就造成整個輪軌系統(tǒng)相對軌道中心線不對稱，所以對輪對、軌道和道床建立整體計算模型。在模型中令鋼軌與車輪在軌枕正上方發(fā)生接觸，為考慮其它軌枕對輪軌接觸狀態(tài)的影響，軌道模型計入三排軌枕。軌道模型中，曲線半徑為1600m，外軌超高量為120mm。在輪軌系統(tǒng)的整體模型中，輪對和內(nèi)外側(cè)鋼軌的接觸按接觸狀態(tài)模擬，而鋼軌、扣件、橡膠墊板、軌枕和道床等的相互關(guān)系按連續(xù)體模擬，各個組成部分根據(jù)其材料性能賦予相應(yīng)的材

8、料特征參數(shù)。軌枕和道床離輪軌接觸區(qū)比較遠，它們對輪軌接觸狀態(tài)的影響較小，為了簡化計算采用普通的實體單元對其進行離散。整個模型采用八節(jié)點等參塊體元，共包括28個子結(jié)構(gòu)，44122個節(jié)點，27432個單元，126個接觸點對。本模型共有兩處接觸區(qū)域，即外側(cè)車輪輪緣與鋼軌接觸和內(nèi)側(cè)車輪踏面與鋼軌接觸，對這兩處接觸區(qū)域都進行精細的網(wǎng)格剖分，最小單元邊長為3mm。機車第1軸與外側(cè)鋼軌貼靠的輪軌系統(tǒng)整體有限元模型如圖2所示。機車通過曲線時，作用在輪對上的力主要是牽引電機產(chǎn)生的牽引力矩、由圓周運動產(chǎn)生的離心力和外軌超高產(chǎn)生的向心力以及輪軌之間的摩擦力。為了比較，計算時軸重分別取21t和23t，作用在兩個車輪的

9、輪心的節(jié)點上。牽引力矩由東風(fēng)11型內(nèi)燃機車的牽引力反推求出，并將其轉(zhuǎn)換成繞車軸表面的一圈切向力，均勻分布于牽引齒輪安裝部位的車軸表面各節(jié)點上2。牽引力根據(jù)東風(fēng)11型內(nèi)燃機車牽引曲線得出，共取120km/h和160km/h兩種速度對應(yīng)的牽引力。橫向力分別取120km/h和160km/h速度對應(yīng)的離心力，作用在車輪的輪心的節(jié)點上。輪軌間的摩擦系數(shù)取為0.45。圖2 輪軌系統(tǒng)整體有限元模型圖四 計算結(jié)果分析 4.1相同的軸重、不同曲線通過速度的影響圖3和圖4是23t軸重機車以120km/h通過曲線時內(nèi)外側(cè)車輪與鋼軌接觸應(yīng)力云圖，最大接觸應(yīng)力（Von mises 應(yīng)力）分別為749MPa和1137MP

10、a。內(nèi)外側(cè)車輪與鋼軌接觸的接觸斑上各接觸點對的最大法向接觸力分別為10909N和13287N，內(nèi)側(cè)車輪與鋼軌法向接觸力合力大約為112475N，接觸力方向與水平夾角為94.8o，外側(cè)車輪與鋼軌法向接觸力合力大約為104258N，接觸力方向與水平夾角為114.1o，內(nèi)外側(cè)車輪與鋼軌接觸斑的面積大約分別是162 mm2和135 mm2。 圖3 23t軸重機車以120km/h通過曲線 圖4 23t軸重機車以120km/h通過曲線時內(nèi)側(cè)車輪與鋼軌接觸應(yīng)力云圖 時外側(cè)車輪與鋼軌接觸應(yīng)力云圖圖5和圖6是23t軸重機車以160km/h通過曲線時內(nèi)外側(cè)車輪與鋼軌接觸應(yīng)力云圖，最大接觸應(yīng)力分別為723MPa和1

11、194MPa。內(nèi)外側(cè)車輪與鋼軌接觸的接觸斑上各接觸點對的最大法向接觸力分別為10704N和12550N，內(nèi)側(cè)車輪與鋼軌法向接觸力合力大約為109294N，接觸力方向與水平夾角為94.8o，外側(cè)車輪與鋼軌法向接觸力合力大約為108201N，接觸力方向與水平夾角為120.7o，內(nèi)外側(cè)車輪與鋼軌接觸斑的面積大約是162mm2和144mm2。 圖5 23t軸重機車以160km/h通過曲線 圖6 23t軸重機車以160km/h通過曲線時內(nèi)側(cè)車輪與鋼軌接觸應(yīng)力云圖 時外側(cè)車輪與鋼軌接觸應(yīng)力云圖分析以上數(shù)據(jù)可得知，機車以120km/h或 160km/h通過曲線時，內(nèi)外側(cè)輪軌接觸部位都發(fā)生嚴重的塑性變形，這會

12、造成內(nèi)外側(cè)輪軌接觸部位磨耗。以160km/h通過曲線時外側(cè)車輪輪緣接觸處最大接觸應(yīng)力比以120km/h通過曲線時大，這說明對于速度高的機車，其輪緣與鋼軌之間將產(chǎn)生更嚴重的塑性變形，這是鐵路提速之后輪緣與鋼軌之間磨耗加劇原因之一；另外，雖然160km/h通過曲線時外側(cè)車輪最大法向接觸力比120km/h時小，但輪軌之間接觸面積增加，所以160km/h通過曲線時外側(cè)車輪法向接觸力合力比120km/h時的合力大，對于速度高的機車，其輪緣與鋼軌之間將產(chǎn)生更多的摩擦功，這是鐵路提速之后輪緣與鋼軌之間磨耗加劇的另一個原因。4.2相同曲線通過速度、不同軸重的影響圖7和圖8是25t軸重機車以120km/h通過曲

13、線時內(nèi)外側(cè)車輪與鋼軌接觸應(yīng)力云圖，最大接觸應(yīng)力分別為775MPa和1168MPa。內(nèi)外側(cè)車輪與鋼軌接觸的接觸斑上各接觸點對的最大法向接觸力分別為11256N和13812N，內(nèi)側(cè)車輪與鋼軌法向接觸力合力大約為122607N，接觸力方向與水平夾角為94.8o，外側(cè)車輪與鋼軌法向接觸力合力大約為112917N，接觸力方向與水平夾角為114.1o，內(nèi)外側(cè)車輪與鋼軌接觸斑面積大約分別是171mm2和144mm2。分析以上數(shù)據(jù)可知，23t軸重和25t軸重機車以120km/h通過曲線時，內(nèi)外側(cè)輪軌接觸部位都發(fā)生嚴重的塑性變形，這將造成內(nèi)外側(cè)輪軌接觸部位磨耗。25t軸重機車的外側(cè)車輪輪緣接觸處最大接觸應(yīng)力比2

14、3t軸重機車更大，這說明當以同樣速度通過曲線時，對于軸重大的機車，其輪緣與鋼軌之間將產(chǎn)生更嚴重的塑性變形，這是鐵路開通重載牽引之后輪緣與鋼軌之間磨耗加劇的重要原因之一；不但如此，25t軸重機車的外側(cè)車輪最大法向接觸力也比23t軸重機車大，而且法向接觸力合力也比23t軸重機車大，所以對于軸重大的機車，其輪緣與鋼軌之間將產(chǎn)生更多的摩擦功，這也是鐵路開通重載列車引起輪緣與鋼軌之間磨耗加劇的原因之一。 圖7 25t軸重機車以120km/h通過曲線 圖8 25t軸重機車以120km/h通過曲線時內(nèi)側(cè)車輪與鋼軌接觸應(yīng)力云圖 時外側(cè)車輪與鋼軌接觸應(yīng)力云圖五 小結(jié)和今后的展望 在列車運行速度大幅度提高的今天，

15、提高機車曲線通過問題的力學(xué)分析精度，以確保列車運行的安全，已變得十分迫切。必須摒棄原來的基于剛體力學(xué)的分析方法，建立輪軌作為變形體的接觸模型進行分析。 用參變量變分原理和由其導(dǎo)出的有限元參數(shù)二次規(guī)劃法，結(jié)合多重多支子結(jié)構(gòu)技術(shù)，按彈性摩擦接觸模型求解機車曲線通過問題是行之有效的。 本文僅僅是用現(xiàn)代計算力學(xué)手段分析機車曲線通過問題的一個初步嘗試，還有許多問題在計算中沒有考慮，像扣件、軌枕、鋼軌等相互接觸關(guān)系對計算結(jié)果的影響，接觸斑局部產(chǎn)生塑性變形的問題，以及列車曲線通過時的動力響應(yīng)問題等，這些都有待于下一步研究解決。 任何理論分析和數(shù)值計算的基礎(chǔ)都是試驗。只有通過大量試驗，并與計算結(jié)果相對比，才能

16、真正建立起可靠的機車車輛曲線通過力學(xué)分析的理論方法和計算模型。因此，呼吁有關(guān)鐵路部門開展相應(yīng)的機車車輛曲線通過試驗研究，以滿足鐵路客運高速化和貨運重載化對設(shè)計的需要。主要參考文獻1鐘萬勰，張洪武，吳承偉. 參變量變分原理及其在工程中的應(yīng)用. 科學(xué)出版社，1997:251-2552張軍. 基于有限元法的輪軌蠕滑理論研究. 大連理工大學(xué)博士論文，2003:41-43吳昌華簡歷 吳昌華，男，1937年11月生于上海。大連交通大學(xué)教授，大連理工大學(xué)工程力學(xué)系客座教授。1961年畢業(yè)于前蘇聯(lián)哈爾科夫工學(xué)院動力機械系內(nèi)燃機車專業(yè)。19611986在大連機車車輛廠設(shè)計部門工作。1980年去聯(lián)邦德國進修。1987年調(diào)入大連鐵道學(xué)院。先后承擔了下列科研項目：鐵道部科研項目“機車和柴油機結(jié)構(gòu)的計算機輔助