盧 勇，王成國(guó)，李 球

（1 廣州地鐵運(yùn)營(yíng)事業(yè)總部 車輛中心，廣東廣州510380；2 中國(guó)鐵道科學(xué)研究院 鐵道科學(xué)技術(shù)研究發(fā)展中心，北京100081）

隨著計(jì)算機(jī)技術(shù)發(fā)展，有限元技術(shù)廣泛應(yīng)用于輪軌接觸力學(xué)研究領(lǐng)域。為了了解廣州地鐵A1型車輪軌接觸特性，本文建立三維輪軌滾動(dòng)接觸有限元模型，計(jì)算分析牽引力和橫移量對(duì)輪軌接觸特性的影響。

1 接觸本構(gòu)模型

接觸本構(gòu)模型包括接觸法向應(yīng)力的定義和接觸切向應(yīng)力的定義。在某一點(diǎn)兩個(gè)接觸面間的接觸壓力p可以定義成相互滲透量h的函數(shù)：

式中K為接觸剛度。

接觸表面切向面力的模型稱之為摩擦模型。目前應(yīng)用最廣泛的Coulomb摩擦模型，該模型源于剛體的摩擦模型，當(dāng)Coulomb摩擦模型用于連續(xù)體時(shí)，它們應(yīng)用在接觸界面的每一點(diǎn)，如果兩接觸面在x處接觸，則

式中tT為接觸表面的切向面力；兩個(gè)物體在某一點(diǎn)處接觸時(shí)，則法向力p≤0，因此－μFp總是正值。式（2）條件（a）表示：當(dāng)接觸點(diǎn)處的切向面力小于臨界值，接觸面之間無(wú)相對(duì)切向運(yùn)動(dòng)，即兩個(gè)物體在該點(diǎn)處是黏著狀態(tài)。式（2）條件（b）表示當(dāng)接觸點(diǎn)處切向力等于臨界值時(shí)，接觸面間將發(fā)生切向相對(duì)滑動(dòng)。切向摩擦的方向和相對(duì)切向速度的方向相反，摩擦力起防止相對(duì)滑動(dòng)的作用。

拉格朗日乘子法或罰函數(shù)法均可將法向和切向接觸條件引入由虛位移原理建立的泛函方程中，并對(duì)其進(jìn)行求解。拉格朗日乘子法可使約束條件得到精確的滿足，但該方法增加了求解方程的自由度，并使求解方程的系數(shù)矩陣含有零對(duì)角元素，給求解帶來(lái)不便。罰函數(shù)法在求解靜力接觸問(wèn)題時(shí)可避免由于系數(shù)矩陣非正定導(dǎo)致的迭代過(guò)程不收斂問(wèn)題。因此，本文有限元模型采用罰函數(shù)法進(jìn)行求解。

2 三維輪軌滾動(dòng)接觸有限元模型

采用輪軌實(shí)際型面建立三維有限元模型，車輪直徑840mm，磨耗型踏面，踏面型面如圖1所示，鋼軌采用國(guó)內(nèi)60kg／m2標(biāo)準(zhǔn)軌（CHN60）。有限元模型共有135 064個(gè)C3D8R六面體單元，在接觸區(qū)網(wǎng)格細(xì)化為1mm，為了減小計(jì)算量，接觸區(qū)采用彈塑性材料，遠(yuǎn)離接觸區(qū)的材料設(shè)置為線彈性材料。在低速、干燥條件下的輪軌摩擦系數(shù)為0.3～0.5，本文摩擦系數(shù)取0.4，泊松比為0.3，輪對(duì)內(nèi)側(cè)距為1 353mm，軌距為1 435mm，軌底坡1／40。鋼軌底部和兩端面全約束。軸重Q為16t，根據(jù)UIC505－5規(guī)定的計(jì)算載荷，車輪受到的垂向作用力為Fz：

式中Q為單個(gè)車輪承受的重力；g為重力加速度，取值為9.81m／s2。

圖1 輪軌接觸有限元模型

起動(dòng)力矩M0是由列車平均起動(dòng)牽引力F和車輪半徑反推得到，M0＝6.23×106N·mm。在車軸上施加垂向載荷FZ和扭矩M0計(jì)算坐標(biāo)如圖1所示，X為橫向，Y為縱向，牽引方向與Y的正方向相同，Z為垂向。

3 計(jì)算結(jié)果分析

3.1 不同牽引力下的輪軌滾動(dòng)接觸特性分析

為了分析牽引力對(duì)輪軌接觸特性的影響，本文分別取牽引扭矩為0、0.25M0、0.5M0、1M0進(jìn)行計(jì)算，計(jì)算結(jié)果數(shù)據(jù)見(jiàn)表1。數(shù)據(jù)顯示，最大縱向摩擦力隨著牽引力增大而顯著增大，不同牽引力作用下接觸斑縱向摩擦力分布見(jiàn)圖2，在沒(méi)有牽引力作用時(shí)，接觸斑內(nèi)有兩個(gè)主要的方向相反的縱向摩擦力中心區(qū)，隨著牽引力的增大，縱向摩擦力都變?yōu)檎^對(duì)值也逐漸增大。圖3為輪軌接觸Mises應(yīng)力分布縱向截面圖，從圖中可看出，在沒(méi)有牽引力時(shí)，輪軌接觸Mises應(yīng)力左右對(duì)稱分布，隨著牽引力的施加，車輪的Mises應(yīng)力區(qū)向前移動(dòng)，鋼軌的Mises應(yīng)力向后偏移，在牽引力為0和M0時(shí)輪軌接觸最大Mises應(yīng)力均出現(xiàn)在接觸面以下約3mm的軌頭內(nèi)。

表1 不同牽引力時(shí)輪軌接觸計(jì)算數(shù)據(jù)

圖2 縱向摩擦力分布

圖3 輪軌接觸Mises應(yīng)力分布縱向截面圖

圖4 黏滑區(qū)分布

根據(jù)公式（2），輪軌接觸斑可分為黏著區(qū)和蠕滑區(qū)。不同牽引力作用下接觸斑黏滑區(qū)分布見(jiàn)圖4，輪軌接觸斑形狀為橢圓形，在沒(méi)有牽引力時(shí)，接觸斑為全黏著區(qū)，隨著牽引力的增大，在接觸斑后部出現(xiàn)月牙形蠕滑區(qū)，接觸斑中前部為黏著區(qū)，牽引力越大，黏著區(qū)越小，蠕滑區(qū)越大。

3.2 不同橫移量下的輪軌接觸特性分析

為了進(jìn)一步了解輪軌匹配情況，本文分別取橫移為－7.5，－5，－2.5，0，2.5，5，7.5mm 進(jìn)行輪軌接觸計(jì)算分析，其中輪緣靠近軌頭為正橫移，遠(yuǎn)離軌頭為負(fù)橫移，輪對(duì)中心線與軌道中心線重合時(shí)輪對(duì)橫移量為0，計(jì)算結(jié)果見(jiàn)表2。數(shù)據(jù)顯示，橫移量在－5～5mm范圍內(nèi)最大接觸壓力、接觸面積和最大Mises應(yīng)力都變化不大。當(dāng)橫移為7.5mm時(shí)，最大接觸壓力顯著減小，接觸面積顯著增大，最大Mises應(yīng)力顯著減小，從圖5（a）可以看到此時(shí)已經(jīng)出現(xiàn)了兩點(diǎn)接觸。除了橫移為7.5 mm時(shí)的接觸斑外，其他工況接觸斑都落在傾斜角為1／46的斜面上，并非圓弧區(qū)與鋼軌接觸，從這個(gè)角度看，該踏面與國(guó)內(nèi)CHN60鋼軌匹配并沒(méi)有完全發(fā)揮磨耗型踏面增大接觸面積，減小接觸應(yīng)力的作用。

當(dāng)兩個(gè)滾筒接觸面存在一個(gè)傾斜角的時(shí)候?qū)l(fā)生自旋蠕滑，本文輪軌模型在橫移量為0時(shí)正是這種接觸狀態(tài)，圖6（a）是Kalker依據(jù)其簡(jiǎn)化理論求得自旋蠕滑情況下接觸斑變形圖，圖6（b）是橫移量為0時(shí)有限元法計(jì)算得到的接觸斑變形圖，二者十分吻合。

表2 不同橫移量時(shí)輪軌接觸計(jì)算數(shù)據(jù)

4 結(jié)論

（1）牽引力改變Mises應(yīng)力在縱向的分布。隨著牽引力增大，車輪Mises應(yīng)力向前偏移，鋼軌Mises應(yīng)力向后偏移。

（2）隨著牽引力的增大接觸斑縱向摩擦力分量明顯增大，接觸斑后部最先出現(xiàn)蠕滑區(qū)，牽引力越大接觸斑黏著區(qū)越小，蠕滑區(qū)越大。

（3）在橫移為－5～5mm范圍內(nèi)，接觸斑分布在車輪踏面傾斜角為1／46的斜面上，并非圓弧區(qū)與鋼軌接觸。從這個(gè)角度看，該踏面與國(guó)內(nèi)CHN60鋼軌匹配并沒(méi)有完全發(fā)揮磨耗型踏面增大接觸面積，減小接觸應(yīng)力的作用。

圖5 車輪踏面接觸位置

圖6 接觸斑變形圖

［1］陳澤深，王成國(guó).鐵道車輛動(dòng)力學(xué)與控制［M］.北京：中國(guó)鐵道科學(xué)出版社，2004.

［2］常崇義.有限元輪軌滾動(dòng)接觸理論及其應(yīng)用研究［D］.北京：中國(guó)鐵道科學(xué)研究院，2010.

［3］肖 乾，車宇翔，周新建，等.CRH2010A綜合檢測(cè)車測(cè)力車輪與鋼軌滾動(dòng)接觸彈塑性分析.［J］.中國(guó)鐵道科學(xué).2012，33（3）：56－60.

［4］XiaoQian Wang Chenguo GuoGe.The Research of Parallel Computing for Large－scale Finite Element Model OF Wheel／Rail Rolling Contact［J］.Proceedings of 2010 3rd IEEE International conference on Computer Science and Information Technology，Chengdu，pp：254－257 ，July.2010.

［5］T Telliskivi，U Olofsson.Contact mechanics analysis of measured wheel－rail profiles using the Finite element method［J］.Proc Instn Mech Engrs，2000，215：65－72.

［6］Xin Zhao，Zili Li and Rolf Dollevoet.Solution of the wheelrail rolling contact in elasticity and elasto－plasticity using a transient finite element model［J］.8th International Conference on Contact Mechanics and Wear of Rail／Wheel Systems（CM2009）.Italy：September 15－18，2009，1－9.

［7］石亦平，周玉蓉.ABAQUS有限元分析實(shí)例詳解［M］.北京：機(jī)械工業(yè)出版社，2008.

［8］國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局.GB 2585－2007.鐵路用熱軋鋼軌［S］.北京：中國(guó)標(biāo)準(zhǔn)出版社，2008.

［9］大野熏.輪軌間滾動(dòng)摩擦及其控制［J］.國(guó)外鐵道車輛，2001，38（1）：36－40.

［10］International Union of Railway.UIC 505－5Techniacal Approval of Solid Wheels［S］.Paris：International Union of Railway，2003.

［11］Kalker J.J..Simplified Theory of Rolling Contact［R］.Delft Progress Report 1，1－10.