磨耗鋼軌經(jīng)濟(jì)性打磨型面研究

林鳳濤，胡偉豪, 2

磨耗鋼軌經(jīng)濟(jì)性打磨型面研究

林鳳濤1，胡偉豪1, 2

(1. 華東交通大學(xué) 機(jī)電與車輛工程學(xué)院，江西 南昌 330013；2. 中鐵物總運(yùn)維科技有限公司，北京 100000)

磨耗鋼軌的恢復(fù)原型打磨方式極易造成鋼軌材料的過(guò)度去除，合理的磨耗鋼軌打磨型面對(duì)于改善列車的輪軌匹配性能、提高列車運(yùn)行安全及延長(zhǎng)鋼軌壽命具有重要的意義。基于NURBS曲線理論建立鋼軌型面曲線的重構(gòu)方法，結(jié)合鋼軌型面特性，選擇標(biāo)準(zhǔn)60 kg/m鋼軌(CN60)型面曲線上的離散化13個(gè)型值點(diǎn)作為設(shè)計(jì)變量，以減少打磨去除材料和減小脫軌系數(shù)為目標(biāo)，以鋼軌型面幾何特性和降低鋼軌磨耗為約束條件，基于Archard磨耗計(jì)算模型，提出貨運(yùn)鐵路鋼軌經(jīng)濟(jì)性打磨型面的設(shè)計(jì)方法。采用優(yōu)化算法進(jìn)行求解，建立輪軌接觸模型和車輛-軌道耦合動(dòng)力學(xué)模型進(jìn)行計(jì)算。研究結(jié)果表明：設(shè)計(jì)的gri60打磨型面具有良好的輪軌接觸特性，與LMa車輪匹配中，軌距角與車輪輪緣部分接觸性能變好；與標(biāo)準(zhǔn)軌面打磨方法相比，材料的去除量減少了12.35%；脫軌系數(shù)的最大值降低了5.63%。同時(shí)，與CN60/LMa輪軌匹配相比，gri60/LMa匹配的左右鋼軌磨耗分別減少28.19%和33.48%；導(dǎo)向輪對(duì)的左右車輪磨耗分別降低11.38%和33.88%。與恢復(fù)標(biāo)準(zhǔn)打磨型面方法相比，本文設(shè)計(jì)的gri60打磨型面具有更好的動(dòng)態(tài)特性和緩磨特性。

鋼軌磨耗；鋼軌打磨；型面優(yōu)化；NURBS

鐵路運(yùn)營(yíng)中出現(xiàn)的輪軌磨耗問題越來(lái)越突出，成為影響列車安全性、運(yùn)行穩(wěn)定性和鋼軌維護(hù)成本的主要因素。設(shè)計(jì)合理的鋼軌打磨型面是保障良好輪軌匹配特性的一種有效手段，能夠改善輪軌接觸關(guān)系，減少輪軌間接觸傷損，進(jìn)而合理控制車輪璇修周期、軌道打磨周期，進(jìn)而控制運(yùn)營(yíng)成本降低運(yùn)維人員工作壓力[1]。對(duì)于鋼軌打磨型面的設(shè)計(jì)，國(guó)內(nèi)外學(xué)者進(jìn)行了相關(guān)研究。Sato[2]設(shè)計(jì)的型面廣泛應(yīng)用在新干線，降低了磨耗，取得了良好的經(jīng)濟(jì)效益。Choi等[3]采用遺傳算法對(duì)地鐵曲線段鋼軌磨耗量進(jìn)行優(yōu)化，使曲線段偏磨現(xiàn)象有所改善。王軍平等[4]提出了個(gè)性化鋼軌型面打磨和設(shè)計(jì)方法，提高打磨效率。王文健等[5]根據(jù)現(xiàn)場(chǎng)調(diào)研情況設(shè)計(jì)了一種新型鋼軌打磨型面，解決了廣深線鋼軌的斜裂紋等傷損問題。周駿等[6]提出了基于軌道質(zhì)量指數(shù)的打磨方法，對(duì)后期鋼軌打磨優(yōu)化具有指導(dǎo)意義。鋼軌型面的描述方法和求解方法國(guó)內(nèi)外專家也進(jìn)行了相關(guān)研究。丁軍君等[7]利用Savitzky-Glay法進(jìn)行平滑處理，并且推薦使用算術(shù)平均法來(lái)擬合鋼軌廓形。王璞等[8]建立輪軌磨耗演化的數(shù)值仿真模型對(duì)75 kg/m鋼軌型面進(jìn)行優(yōu)化設(shè)計(jì)。Gerlici等[9]提出了一種通過(guò)迭代方法，改變剖面圓弧半徑的鋼軌型面設(shè)計(jì)方法。Ingemar等[10]基于輪軌動(dòng)力學(xué)理論和數(shù)值模擬理論，提出用遺傳算法對(duì)曲線鋼軌進(jìn)行型面優(yōu)化。Choi等[11]通過(guò)遺傳算法優(yōu)化非對(duì)稱軌頭型面，降低了小曲線半徑鋼軌磨耗。LIU等[12]提出了基于NURBS曲線的型面優(yōu)化設(shè)計(jì)方法，針對(duì)確定的滾動(dòng)半徑差得到優(yōu)化的鋼軌型面。以上方法雖然對(duì)鋼軌型面進(jìn)行了優(yōu)化，但是大多型面設(shè)計(jì)方法中存在擬合階次過(guò)高、型面曲線局部特性不明顯、平滑性描述困難、多段圓弧擬合誤差大等缺點(diǎn)，同時(shí)未考慮到材料去除量等問題。本文基于NURBS擬合算法，優(yōu)化鋼軌廓形打磨區(qū)域曲線，實(shí)現(xiàn)少量型值點(diǎn)描述鋼軌型面曲線的方法，并驗(yàn)證該方法的有效性。同時(shí)以磨耗后鋼軌型面為研究對(duì)象，以鋼軌型面幾何特性和降低鋼軌磨耗為約束條件，建立鋼軌打磨型面曲線的多目標(biāo)優(yōu)化模型，同時(shí)建立相關(guān)的約束函數(shù)，采用優(yōu)化算法對(duì)模型進(jìn)行求解計(jì)算，以獲得良好的打磨型面。并通過(guò)Archard磨耗模型計(jì)算磨耗，在減少鋼軌型面材料去除量的同時(shí)，獲得良好的輪軌匹配、減緩鋼軌的磨損磨耗。

1 3次NURBS曲線模型

鋼軌型面曲線是一種特殊形式的曲線，除了有平滑連續(xù)性要求外，其磨耗形式的局部變化特性也極其重要。NURBS擬合算法能夠更好地控制物體表面的曲線度，同時(shí)具有局部修改性、強(qiáng)凸包性、曲線段連接點(diǎn)處具有C1連續(xù)等特性。

1.1 NURBS曲線的定義

由+1個(gè)多邊形控制頂點(diǎn)定義的一條次NURBS曲線，可以表示為一條分段有理多項(xiàng)式矢函數(shù)，形式如下：

按照上述定義，=3就是一條3次NURBS曲線的數(shù)學(xué)表達(dá)式。

1.2 3次NURBS曲線控制點(diǎn)反求算法

采用鋼軌型面上的有限離散點(diǎn)(型值點(diǎn))建立鋼軌型面曲線的NURBS描述方法，即已知個(gè)型值點(diǎn)p及其相應(yīng)的2個(gè)權(quán)因子(,…,1)，需進(jìn)行控制點(diǎn)的反算設(shè)計(jì)。設(shè)控制點(diǎn)權(quán)因子w與型值點(diǎn)權(quán)因子h有下述關(guān)系：

根據(jù)上式可以列出個(gè)方程，本文針對(duì)鋼軌型面的曲線特性，補(bǔ)充2切矢邊界條件：

這樣，+2個(gè)控制頂點(diǎn)即可確定。

1.3 鋼軌型面3次NURBS描述方法驗(yàn)證

根據(jù)鋼軌型面曲線特點(diǎn)，采取與曲率正比關(guān)聯(lián)的方式選取型值點(diǎn)，集成反求樣條曲線的控制定點(diǎn)坐標(biāo)及權(quán)因子計(jì)算，建立鋼軌型面曲線3次NURBS描述的通用方法。

1.3.1 確定型值點(diǎn)

描述曲線的精確程度直接受型值點(diǎn)個(gè)數(shù)及型值點(diǎn)的確定方法影響。對(duì)于鋼軌型面曲線的規(guī)律探究，發(fā)現(xiàn)鋼軌型面曲線的軌距角部分斜率變化較大，軌頂部分斜率變化較小，所以在軌距角部分選取密集的型值點(diǎn)有利于鋼軌型面曲線整體精確的表達(dá)?？紤]到鋼軌廓形軌距角附近區(qū)域斜率變化比較大，軌頂部分曲率變化比較小，從控制計(jì)算量的角度出發(fā)密集軌距角附近區(qū)域的型值點(diǎn)。

對(duì)鋼軌橫坐標(biāo)=[?36,36] mm段的曲線采用離散化取點(diǎn)的方法進(jìn)行取點(diǎn)。取的型值點(diǎn)與鋼軌型面曲線的曲率進(jìn)行正比關(guān)聯(lián)，相對(duì)增加軌距角曲線部分的型值點(diǎn)的采樣數(shù)量和采樣頻率，同時(shí)取點(diǎn)過(guò)還需要考慮NURBS有限型值點(diǎn)描述復(fù)雜曲線的特性，對(duì)鋼軌型面曲線分別進(jìn)行9，13和17個(gè)型值點(diǎn)3種情況進(jìn)行對(duì)比分析，從而確定型值點(diǎn)適宜個(gè)數(shù)。

1.3.2 型值點(diǎn)權(quán)因子

型值點(diǎn)在曲線中的重要程度通過(guò)型值點(diǎn)的權(quán)因子反映出來(lái)，某個(gè)性質(zhì)點(diǎn)的權(quán)因子相對(duì)來(lái)說(shuō)比較大，則說(shuō)明此型值點(diǎn)對(duì)附近點(diǎn)和曲線的影響程度比較大，反之亦然。由于鋼軌型面軌距角部分對(duì)于輪軌匹配有著重要意義，軌距角部分=[?34,?23] mm和=[23,34] mm內(nèi)型值點(diǎn)權(quán)因子設(shè)定范圍為[0.9,1.0]，軌頂部分=[?23,23] mm內(nèi)型值點(diǎn)權(quán)因子值設(shè)定范圍為[0.8,0.9]，其余部分型值點(diǎn)權(quán)因子設(shè)定范圍為[0.5,0.7]。

對(duì)CN60鋼軌型面進(jìn)行=9，13，17型值點(diǎn)的3種情況進(jìn)行NURBS擬合設(shè)計(jì)，在檢驗(yàn)該方法可行性的同時(shí)對(duì)標(biāo)準(zhǔn)鋼軌型面CN60曲線與3次NURBS擬合設(shè)計(jì)曲線f60進(jìn)行相關(guān)性分析，采用相關(guān)系數(shù)參數(shù)rail評(píng)價(jià)該擬合方法的有效性，結(jié)果如表1所示。

表1 3次NURBS設(shè)計(jì)前后對(duì)比結(jié)果

從結(jié)果分析可知，型值點(diǎn)個(gè)數(shù)越多，所得差值擬合曲線與原曲線越接近，但是其擬合曲線所需時(shí)間和計(jì)算成本越高。綜合考慮曲線的精確性及型值點(diǎn)個(gè)數(shù)對(duì)計(jì)算的效率影響，選擇=13個(gè)型值點(diǎn)進(jìn)行差值擬合。CN60鋼軌型面與重構(gòu)型面對(duì)比如圖1所示?？梢钥闯鰯M合型面f60與CN60磨耗型面基本吻合，可采取NURBS曲線設(shè)計(jì)方法，實(shí)現(xiàn)鋼軌磨耗型面描述良好的描述效果。

式中，CN60鋼軌型面曲線擬合主要參數(shù)型值點(diǎn)P，權(quán)因子，及反算求得的控制點(diǎn)d值如表2所示。

圖1 Rf60與原CN60磨耗型面對(duì)比

表2 CN60鋼軌型面3次NURBS參數(shù)

2 經(jīng)濟(jì)性磨耗鋼軌型面設(shè)計(jì)模型

2.1 設(shè)計(jì)變量

2.2 目標(biāo)函數(shù)

2.2.1 磨耗鋼軌打磨去除材料最少目標(biāo)函數(shù)

2.2.2 脫軌系數(shù)目標(biāo)函數(shù)

根據(jù)公式，建立脫軌系數(shù)目標(biāo)函數(shù)為：

式中：和分別輪軌橫向力和輪軌垂向力；1為輪緣角，對(duì)脫軌系數(shù)取均方根值。

2.3 約束函數(shù)

在描述過(guò)程中，為了使描述曲線光滑，沒有奇異點(diǎn)和波浪點(diǎn)，選擇磨耗鋼軌型面的縱坐標(biāo)統(tǒng)計(jì)量、曲線單調(diào)性、凹凸特性作為鋼軌型面曲線的幾何約束條件。同時(shí)選取降低軌道磨耗作為約束條件。

2.3.1 型值點(diǎn)的縱坐標(biāo)范圍約束條件

選取某運(yùn)營(yíng)線路磨耗鋼軌統(tǒng)計(jì)型面和標(biāo)準(zhǔn)CN60鋼軌型面作為設(shè)計(jì)變量的上下邊界條件:

2.3.2 鋼軌型面曲線的單調(diào)性約束條件

2.3.3 鋼軌型面曲線凹凸性約束條件

基于鋼軌型面凹凸特性分析，設(shè)定型面曲線整體為凸的形式，約束條件為：

2.3.4 降低軌道磨耗約束條件

基于列車輪軌磨耗特性，考察左右鋼軌累計(jì)磨耗量，約束條件為：

式中：W()，W()分別為左右鋼軌的磨耗功；整個(gè)仿真的運(yùn)營(yíng)線路長(zhǎng)度。

3 Archard磨耗模型

鋼軌和車輪磨耗問題，大多采用摩擦磨損領(lǐng)域廣泛使用的滑動(dòng)磨損模型Archard磨耗模型計(jì)算，這是一種與材料密切相關(guān)的磨耗模型，其主要目的就是解決滾動(dòng)和滑動(dòng)接觸的問題。模型中材料磨耗的體積與法向力、滑動(dòng)距離成正比，與材料的強(qiáng)度成反比，如式(16)所示。

式中：wear為材料的磨損體積；為無(wú)量綱的磨耗系數(shù)，其值由單元上的滑動(dòng)量和正壓力決定；為材料的硬度；F為輪軌法向力；為滑動(dòng)距離，其值由式(17)確定。

對(duì)長(zhǎng)度為Δ和Δ的接觸斑內(nèi)的任一單元，其磨耗深度Δ為：

式中：Δ為單位時(shí)間增量Δ內(nèi)的滑動(dòng)距離，p為單元上的正壓力密度。

由上式進(jìn)一步得到磨耗深度Δ：

V表示某一點(diǎn)通過(guò)接觸斑的速度。忽略因彈性變形產(chǎn)生的速度分量，式(22)可寫為：

式中：V為車輪滾動(dòng)速度。

在磨耗表達(dá)式中引入磨耗系數(shù),該系數(shù)通過(guò)實(shí)驗(yàn)獲得，3個(gè)區(qū)域的數(shù)值如圖2所示。

圖2 Archard模型磨耗區(qū)的磨耗系數(shù)

4 打磨型面計(jì)算與分析

4.1 動(dòng)力學(xué)模型建立及優(yōu)化流程

建立車輛?軌道耦合動(dòng)力學(xué)模型，模型考慮輪軌接觸幾何、輪軌蠕滑特性、車輛懸掛等非線性環(huán)節(jié)。其動(dòng)力學(xué)模型拓?fù)鋱D如圖3所示，車輛模型主要由1個(gè)車體、2個(gè)搖枕、4個(gè)側(cè)架和4個(gè)輪對(duì)以及一系和二系懸掛組成。經(jīng)過(guò)驗(yàn)證，確認(rèn)動(dòng)力學(xué)模型各部分模擬無(wú)誤[13]。

根據(jù)上文所述的鋼軌型面描述方法，形成以鋼軌型值點(diǎn)縱坐標(biāo)為設(shè)計(jì)變量，以鋼軌幾何特性以及磨耗特性為約束條件，以減小鋼軌打磨去除量和脫軌系數(shù)為目標(biāo)函數(shù)的鋼軌型面設(shè)計(jì)方法。通過(guò)優(yōu)化算法對(duì)目標(biāo)函數(shù)求解，然后通過(guò)動(dòng)力學(xué)模型以及輪軌接觸特性對(duì)求解結(jié)果進(jìn)行驗(yàn)證，最終得到符合條件的鋼軌型面[14]。

圖3 車輛動(dòng)力學(xué)拓?fù)鋱D

4.2 計(jì)算結(jié)果

最終求得的gri60曲線與原CN60標(biāo)準(zhǔn)型面對(duì)比如圖4所示，可見在鋼軌側(cè)面∈[?20,20]內(nèi)，gri60與CN60曲線形式一致。在橫坐標(biāo)位置在[?35,?30] mm和[30,35] mm內(nèi)的軌距角范圍內(nèi)，gri60型面的縱坐標(biāo)較小，曲線斜率變化稍快，二者縱坐標(biāo)最大差值為1.2 mm。在橫坐標(biāo)在[?15.2,15.2] mm的軌頂范圍內(nèi)，兩曲線變化都平緩，gri60曲線縱坐標(biāo)值稍小，斜率基本相同。

圖4 鋼軌打磨型面Rgri60曲線

4.3 材料去除量

對(duì)gri60和CN60與磨耗鋼軌型面圍成的截面面積即打磨去除材料面積進(jìn)行分析，如圖5所示，式中因受到輪軌力作用使得鋼軌worn60的軌距角部分發(fā)生磨耗變形，若按照CN60型面進(jìn)行打磨，為了還原軌距角部分，需要打磨掉軌頂豎直高度H2為2.25 mm，打磨型面gri60的打磨豎直高度H1為1.64 mm，減小了0.61 mm。采用gri60型面打磨后，材料去除量由原來(lái)145.6 mm2降低到111.2 mm2，截面的材料去除量減少了23.64%?？傮w來(lái)說(shuō)，gri60打磨型面有利于保留和利用原磨耗鋼軌剩余型面部分。

圖5 打磨量對(duì)比分析

4.4 輪軌接觸幾何

CN60型面和gri60型面輪軌接觸分析，設(shè)置軌底坡為1/40，輪對(duì)直徑860 mm，軌距為1 435 mm，如圖6所示。gri60/LMa匹配接觸點(diǎn)在分布的寬度和范圍上小于CN60/LMa匹配，在橫移量為[?10,10] mm的范圍內(nèi)，與輪軌接觸良好，且接觸點(diǎn)均勻的分布在軌頂部分，無(wú)明顯的跳躍點(diǎn)。而CN60/LMa匹配在橫移量為[8,10] mm范圍內(nèi)發(fā)生輪緣接觸，且存在輪軌跳躍點(diǎn)。整體上，gri60/LMA匹配性能優(yōu)于CN60/LMa匹配。

(a) Rgri60/LMa匹配；(b) CN60/LMA匹配

等效錐度參數(shù)對(duì)于車輛的動(dòng)力學(xué)性能有重要的影響，gri60/LMa匹配和CN60/LMa匹配的等效錐度曲線如圖7所示。對(duì)比發(fā)現(xiàn)，在橫移量[2,9] mm的范圍內(nèi)，gri60/LMa匹配比CN60/LMa匹配的等效錐度值略大，橫移量為6 mm時(shí)，二者最大差值為0.01，總體上，等效錐度曲線變化趨勢(shì)相同，都在橫移量為[0,7] mm的范圍內(nèi)，等效錐度都小于0.07，且呈現(xiàn)“緩慢”上升，這一趨勢(shì)有利于列車高速運(yùn)行。在橫移量為[7,10] mm內(nèi)，等效錐度呈現(xiàn)快速上升趨勢(shì)，有利于輪對(duì)自動(dòng)對(duì)中，提高曲線通過(guò)性。

4.5 動(dòng)力學(xué)分析

貨運(yùn)線路多數(shù)為小半徑曲線，設(shè)定=800 m的右曲線工況，計(jì)算gri60/LMa匹配和CN60/LMa匹配下的貨車橫向平穩(wěn)性指標(biāo)、一位輪對(duì)輪軌橫向力、脫軌系數(shù)及輪重減載率，計(jì)算結(jié)果如圖8所示。根據(jù)GB5599—85動(dòng)力學(xué)標(biāo)準(zhǔn)對(duì)計(jì)算得到的動(dòng)力學(xué)指標(biāo)進(jìn)行評(píng)價(jià)。計(jì)算結(jié)果表明，對(duì)于CN60/LMa來(lái)說(shuō)，gri60/LMa的脫軌系數(shù)最大值由0.176減小到0.166，減少了5.68%；輪重減載率最大值由0.249減小到0.218，減少了12.45%；車體橫向振動(dòng)加速度的最大值由1.37減小到1.16，減少了15.33%；一位輪對(duì)輪軌橫向力的最大值由20.6 kN減小到17.1 kN，減少了17%；整體來(lái)說(shuō)，gri60型面動(dòng)力學(xué)性能相較于CN60型面有一定提升，可以保證列車運(yùn)行的安全性和穩(wěn)定性。

圖7 等效錐度對(duì)比

4.6 磨耗分析

采用上述Archard磨耗模型分別計(jì)算gri60/ LMa和CN60/LMa下，通過(guò)次數(shù)為20萬(wàn)次和60萬(wàn)次時(shí)，左右鋼軌的磨耗量對(duì)比情況。通過(guò)計(jì)算曲線與坐標(biāo)軸所圍面積表示磨耗面積；曲線大于零的位置即為磨耗位置；曲線最大值為磨耗最深位置。根據(jù)圖9可以看出左軌通過(guò)20萬(wàn)次時(shí)，gri60/LMA和CN60/LMA的磨耗面積分別為5.81 mm2和8.08 mm2，面積減少了28.19%；磨耗寬度由[?9,32] mm減小到[?8,29] mm，且磨耗深度最深的位置由27 mm移動(dòng)到14 mm處，有效減小軌距角磨耗和側(cè)面磨耗，有利于延長(zhǎng)軌道的使用壽命。在通過(guò)次數(shù)為60萬(wàn)次時(shí)，gri60/LMA和CN60/LMA下磨耗面積分別為12.11 mm2和15.46 mm2，面積減少了21.68%；磨耗寬度由原來(lái)的[?11,35] mm減小到[?10,33] mm。右軌在通過(guò)量次數(shù)為20萬(wàn)次時(shí)gri60/LMA和CN60/LMA下磨耗面積分別為2.73 mm2和4.11 mm2，面積減少了33.48%；在通過(guò)次數(shù)為60萬(wàn)次時(shí)磨耗面積分別為5.49 mm2和8.04 mm2，面積減少了31.7%。累計(jì)磨耗量和磨耗寬度都有明顯的優(yōu)化。

(a) 脫軌系數(shù)對(duì)比；(b) 輪重減載率對(duì)比；(c) 車體橫向振動(dòng)加速度對(duì)比；(d) 輪軌橫向力對(duì)比

(a) 左側(cè)鋼軌磨耗；(b) 右側(cè)鋼軌磨耗

在對(duì)比分析鋼軌磨耗的同時(shí)，計(jì)算2萬(wàn)和5萬(wàn)km運(yùn)營(yíng)里程下，gri60/LMA匹配和CN60/LMA匹配的車輪磨耗情況，對(duì)于導(dǎo)向輪左輪在2萬(wàn)km下累計(jì)磨耗面積分別為5.475 mm2和6.157mm2，面積減少11.38%；5萬(wàn)km的累計(jì)磨耗面積分別為10.834 mm2和13.162 mm2面積減少17.69%；導(dǎo)向輪右輪在2萬(wàn)km下累計(jì)磨耗面積分別為2.525 mm2和3.819 mm2，面積減少33.88%；5萬(wàn)km的累計(jì)磨耗面積分別為5.965 mm2和8.993 mm2面積減少33.68%；gri60型面相較于CN60標(biāo)準(zhǔn)型面有較大改善。

(a) 一位輪對(duì)左側(cè)車輪磨耗圖；(b) 一位輪對(duì)右側(cè)車輪磨耗圖

5 結(jié)論

1) 對(duì)CN60鋼軌型面型值點(diǎn)選取采用與曲率相關(guān)聯(lián)的方法，分成10，14和18段弧長(zhǎng)，取得9，13和17個(gè)型值點(diǎn)，采用非均勻有理B樣條(NURBS)曲線理論重構(gòu)CN60型面曲線，與原曲線相關(guān)系數(shù)為0.89，0.94和0.98，對(duì)比分析后選取14段分割法即取=13個(gè)型值點(diǎn)，能夠滿足鋼軌型面良好 描述。

2) 建立以鋼軌型面上13個(gè)型值點(diǎn)縱坐標(biāo)為設(shè)計(jì)變量，以磨耗鋼軌打磨去除材料最少和減小脫軌系數(shù)為目標(biāo)函數(shù)，以軌道曲線幾何和鋼軌累計(jì)磨耗為約束條件的鋼軌型面多目標(biāo)設(shè)計(jì)函數(shù)。

3) 采用優(yōu)化算法求解，建立輪軌接觸模型和車輛?軌道耦合動(dòng)力學(xué)模型計(jì)算，結(jié)果表明：設(shè)計(jì)打磨型面gri60具有良好的輪軌接觸特性，采用gri60型面后，脫軌系數(shù)的最大值降低了5.63%；鋼軌線路磨耗比標(biāo)準(zhǔn)型面降低了28.19%；與標(biāo)準(zhǔn)軌面打磨方法相比，gri60打磨型面的材料去除量減少了12.35%。

4) 打磨型面gri60與標(biāo)準(zhǔn)的CN60型面的對(duì)比，gri60鋼軌型面在通過(guò)小曲率半徑情況下鋼軌磨耗和車輪磨耗較大幅度的降低；磨耗位置有明顯的改善；同時(shí)能夠減小列車車輪和鋼軌的損傷，延長(zhǎng)鋼軌使用壽命。

[1] 黃育斌, 王文健, 郭俊, 等. 鋼軌非對(duì)稱打磨對(duì)車輛運(yùn)行性能影響[J]. 鐵道科學(xué)與工程學(xué)報(bào), 2015, 12(4): 264?268. HUANG Yubin, WANG Wenjian, GUO Jun, et al. Effect of asymmetrical rail grinding on running performance of vehicles[J]. Journal of Railway Science and Engineering, 2015,12 (4): 264?268.

[2] Sato Yoshihiko. Historial study on designing Japanese rail profile[J]. Wear, 2005, 258(7?8): 1064?1070.

[3] Choi Ha Young, LEE Dong Hyong, Song Chang Yong, et al. Optimization of rail profile to reduce wear on curved track[J]. International Journal of Precision Engineering and Manufacturing, 2013, 14(4): 619?625.

[4] 王軍平, 單連琨, 丁軍君, 等. 個(gè)性化鋼軌廓形打磨方法分析[J]. 鐵道建筑, 2015, 55(11): 131?133.WANG Junping, SHAN Liankun, DING Junjun, et al. Analysis of personalized rail profile grinding methods[J]. Railway Engineering, 2015, 55(11): 131?133.

[5] 王文健, 陳明韜, 郭俊, 等. 高速鐵路鋼軌打磨技術(shù)及其應(yīng)用[J]. 西南交通大學(xué)學(xué)報(bào), 2007, 42(5): 574?577. WANG Wenjian, CHEN Mingtao, GUO Jun, et al. Rail grinding technique and its application in high-speed railway[J]. Journal of Southwest Jiaotong University, 2007, 42(5): 574?577.

[6] 周駿, 劉林芽, 毛順茂. 基于軌道質(zhì)量指數(shù)的鋼軌打磨型面研究[J]. 鐵道建筑, 2014(8): 86?89. ZHOU Jun, LIU Linya, MAO Shunmao. Study on rail profile grinding based on track quality index[J]. Railway Engineering, 2014(8): 86?89.

[7] 丁軍君, 吳朋朋, 王軍平, 等. 基于輪軌關(guān)系的鋼軌打磨代表廓形計(jì)算方法研究[J]. 鐵道學(xué)報(bào), 2019, 41(7): 135?140. DING Junjun, WU Pengpeng, WANG Junping, et al. Study on algorithm of representative profile for rail grinding based on wheel-rail relationship[J]. Journal of the China Railway Society, 2019, 41(7): 135?140.

[8] 王璞, 高亮, 蔡小培. 重載鐵路鋼軌磨耗演變過(guò)程的數(shù)值模擬[J]. 鐵道學(xué)報(bào), 2014, 36(10): 70?75. WANG Pu, GAO Liang, CAI Xiaopei. Numerical simulation of rail wear evolution of heavy-hual railways[J]. Journal of the China Railway Society, 2014, 36(10): 70?75.

[9] Gerlici J, Lack T. Railway wheel and rail head profiles development based on the geometric characteristics shapes[J]. Wear, 2011, 271(1): 246?258.

[10] Ingemar Persson, Rickard Nilsson, Ulf Bik, et al. Use of a genetic algorithm to improve the rail profile on Stockholm underground[J]. Vehicle System Dynamics, 2010, 48(Supp1): 89?104.

[11] Choi H Y, Lee Dong Hyong. Optimization of rail profile to reduce wear on curved track[J]. International Journal of Precision Engineering & Manufacturing, 2013, 14(4): 619?625.

[12] LIU Y, LIU X. Railway wheel profile optimization design based on nurbs curve[C]// International Conference on Computer Modeling and Simulation. IEEE, 2010: 331?335.

[13] LIN Fengtao, ZHOU Shuang, DONG Xiaoqing, et al. Design method of LM thin flange wheel profile based on NURBS[J]. Vehicle System Dynamics, 2019:1?16.

[14] 林鳳濤. 高速列車車輪磨耗及型面優(yōu)化研究[D]. 北京:中國(guó)鐵道科學(xué)研究院, 2014. LIN Fengtao. Research on wheel wear and wheel profile optimization of high speed train[D]. Beijing: China Academy of Railway Sciences, 2014.

Study on the economical grinding surface of wear rail

LIN Fengtao1, HU Weihao1, 2

(1. School of Mechatronics and Vehicle Engineering, East China Jiaotong University, Nanchang 330013, China;2. China Railway Mateirals Operations and Maintenance Technology Co., Ltd., Beijing 100000, China)

The way to restore the prototype of worn rail will lead to excessive removal of rail material. Therefore, reasonable grinding profile of worn wear is of great significance to improve wheel-rail matching, improve the safety of train operation and extend rail life. Based on the theory of non-uniform rational B spline (NURBS) curve, the reconstruction method of rail profile curve was established. In combination with the characteristics of rail profile, 13 data points on the rail profile curve of CN60 were selected as design variables. Aiming at reducing grinding and removing materials and reducing derailment coefficient, taking geometric characteristics of rail profile and reducing rail wear as constraint conditions, the design method of economical grinding rail profile of heavy load railway was established by Archard wear model. The wheel-rail contact model and vehicle-track coupling dynamic model were established. The results show that the designedgri60 grinded profile has good wheel-rail contact characteristics and when matching with LMa wheels, the contact performance between track gauge angle and wheel rim is improved. Compared with the standard rail surface grinding method, the material removal is decreased by 32.35% less. The maximum derailment coefficient is decreased by 5.63%. Besides that, when comparing with the CN60/LMa rail-wheel match, the left and right rail wear ofgri60/LMa match is decreased by 28.19% and 33.48% respectively and the left and right guide wheel wear ofgri60/LMa match is decreased by 11.38% and 33.88% respectively. Compared with the method of recovering standard grinding profile, the designed Rgri60 grinding profile has better dynamic characteristics and slow wear characteristics.

rail wear; profile optimization; particle swarm optimization; NURBS

U271.91

A

1672 ? 7029(2020)10 ? 2493 ? 10

10.19713/j.cnki.43?1423/u.T20191152

2019?12?19

國(guó)家自然科學(xué)基金資助項(xiàng)目(51665014)

林鳳濤(1977?)，男，內(nèi)蒙古赤峰人，副教授，博士，從事輪軌關(guān)系研究；E?mail：46473697@qq.com

(編輯 涂鵬)