基于摩擦自激振動理論的地鐵線路梯形軌枕支承小半徑曲線輪軌系統(tǒng)鋼軌波磨研究

崔曉璐 何志強 程 梽 楊宗超 楊宏光

(1.重慶交通大學機電與車輛工程學院， 400074， 重慶; 2.重慶市軌道交通(集團)有限公司， 401120， 重慶;3.中鐵物總運維科技有限公司， 100036， 北京; 4.中車長春軌道客車股份有限公司， 130062， 長春∥第一作者， 副教授)

0 引言

鋼軌波磨作為常見的軌道損傷問題之一，是鋼軌頂面沿縱向出現(xiàn)的一種規(guī)律性、似波浪狀的不平順磨損現(xiàn)象[1]。其產(chǎn)生的噪聲不僅對乘客的乘車舒適性造成極大影響，而且還會使車輛和軌道結(jié)構(gòu)部件的使用壽命大為降低，極大增加了鐵路工務(wù)維修費用。目前，對于鐵路軌道上出現(xiàn)的鋼軌波磨問題尚未找到根本解決方法，但提出了一些相應(yīng)的控制措施[2]。

現(xiàn)階段關(guān)于鋼軌波磨的成因，大致可分成以下兩種觀點[3]：①鋼軌頂面的粗糙度導致輪軌間瞬時動態(tài)接觸特性波動，進而產(chǎn)生波磨病害; ②輪軌間黏滑振動導致輪軌摩擦自激產(chǎn)生，進而導致摩擦功波動和相應(yīng)的波磨病害。

軌道支承結(jié)構(gòu)常被作為研究鋼軌波磨病害的重要對象。基于上述兩種理論，研究者們以軌道支承結(jié)構(gòu)為著入點，開展了一系列的研究工作，并提出了相應(yīng)緩解波磨病害的措施。文獻[4]從頻域角度出發(fā)，調(diào)查了北京地鐵隧道、橋梁上減振器扣件軌道結(jié)構(gòu)上的鋼軌波磨問題，發(fā)現(xiàn)鋼軌波磨與扣件的剛度阻尼有關(guān)，通過優(yōu)化扣件參數(shù)可以減緩鋼軌波磨的發(fā)展。文獻[5]基于固定波長機理，建立了非線性鋼軌波磨預測模型，發(fā)現(xiàn)降低軌墊剛度可以降低磨耗速率，緩解鋼軌波磨發(fā)展。文獻[6]利用頻域波磨模型，對影響鋼軌波磨的主要軌道支承結(jié)構(gòu)參數(shù)進行了大量調(diào)查分析，發(fā)現(xiàn)軌枕質(zhì)量、軌枕間距、軌墊橫向剛度、軌墊垂向剛度、道床垂向剛度對鋼軌波磨的影響較大，通過優(yōu)化軌枕質(zhì)量、軌枕間距、軌墊垂向剛度參數(shù)可以減緩鋼軌波磨發(fā)生趨勢。文獻[7]基于輪軌磨損理論與瞬態(tài)動力學的觀點，利用線性鋼軌波磨的預測模型,發(fā)現(xiàn)輪軌頂面的粗糙度與鋼軌波磨病害的產(chǎn)生緊密相關(guān)，導致波磨產(chǎn)生的直接原因是車輛軌道結(jié)構(gòu)部分產(chǎn)生共振，通過采用連續(xù)的軌道支承結(jié)構(gòu)能有效避免波磨病害。

本文基于摩擦自激振動理論研究某地鐵線路梯形軌枕支承小半徑曲線軌道上出現(xiàn)的波磨病害，通過對梯形軌枕支承小半徑曲線輪軌系統(tǒng)有限元模型的復特征值分析，研究了軌道減振結(jié)構(gòu)中側(cè)向緩沖墊、底部橡膠墊板參數(shù)與扣件參數(shù)對輪軌摩擦自激振動發(fā)生趨勢的影響。同時，提出了相應(yīng)的緩解波磨病害的措施。

1 地鐵梯形軌枕支承小半徑曲線輪軌系統(tǒng)仿真模型及理論方法

1.1 梯形軌枕支承小半徑曲線輪軌系統(tǒng)接觸模型

梯形軌枕是地鐵線路上常用的一種軌道支承結(jié)構(gòu)，特別應(yīng)用于對減振要求較高的路段。梯形軌枕軌道主要由鋼軌、扣件系統(tǒng)、連接鋼管、道床、側(cè)向緩沖墊、底部橡膠墊板和支撐塊組成[8]。側(cè)向緩沖墊和底部橡膠墊板作為道床與支撐塊之間的減振緩沖結(jié)構(gòu)，又被統(tǒng)稱為軌道減振結(jié)構(gòu)。

根據(jù)梯形軌枕支承小半徑曲線輪軌系統(tǒng)的結(jié)構(gòu)形式，結(jié)合SIMPACK動力學仿真模型，建立梯形軌枕支承小半徑曲線輪軌系統(tǒng)的接觸模型，如圖1所示。通常而言，R(曲線半徑)≤350 m的曲線稱之為小半徑曲線[3]。當?shù)罔F車輛通過梯形軌枕支承小半徑曲線軌道時，輪對兩端外側(cè)受到的垂向懸掛力和橫向懸掛力分別為FSVL和FSLL。輪對兩端內(nèi)側(cè)受到的垂向懸掛力和橫向懸掛力分別為FSVR和FSLR。外輪法向力為NL，蠕滑力為FL；內(nèi)輪法向力為NR，蠕滑力為FR[3]?？奂到y(tǒng)、側(cè)向緩沖墊和底部橡膠墊板利用接地彈簧和阻尼單元模擬。外輪、內(nèi)輪與鋼軌間的接觸角分別為δL、δR?？奂拇瓜騽偠取⒋瓜蜃枘岱謩e為KRV和CRV。軌道減振結(jié)構(gòu)中側(cè)向緩沖墊的橫向剛度、橫向阻尼分別為KSL、CSL，底部橡膠墊板的垂向剛度、垂向阻尼分別為KSV、CSV。其中，KSL與KSV數(shù)值相同，CSL與CSV數(shù)值相同[8]。

圖1 梯形軌枕支承小半徑曲線輪軌系統(tǒng)接觸模型Fig.1 Contact model of the wheel-rail on a sharp curved track supported by trapezoid sleepers

1.2 梯形軌枕支承小半徑曲線輪軌系統(tǒng)有限元模型

根據(jù)標準型梯形軌枕軌道結(jié)構(gòu)以及輪軌接觸模型，利用ABAQUS有限元軟件建立梯形軌枕支承小半徑曲線輪軌系統(tǒng)的有限元模型，如圖2所示。

圖2 梯形軌枕支承小半徑曲線輪軌系統(tǒng)有限元模型Fig.2 Finite element model of the wheel-rail on a sharp curved track supported by trapezoid sleepers

該有限元模型中的各部件分別為輪對、鋼軌和軌道支承結(jié)構(gòu)，其材料參數(shù)取值如表1所示[3]。

表1 梯形軌枕支承小半徑曲線輪軌系統(tǒng)各部件參數(shù)

模型中各部件尺寸為實際尺寸，采用固體單元模擬。由于導向輪對對波磨病害的影響最大，因此模型中的輪對特指導向輪對。鋼軌為長30 m的標準鋼軌，其兩端邊界設(shè)置為6個自由度全約束。輪對與鋼軌的接觸細節(jié)如圖3所示，輪軌摩擦系數(shù)設(shè)置為0.45。為了使模型計算結(jié)果更好收斂，需要對輪對與鋼軌接觸區(qū)域的網(wǎng)格作細化處理。在軌道支承結(jié)構(gòu)中，扣件間距為625 mm，軌距為1 435 mm?？奂到y(tǒng)簡化為彈簧阻尼單元，由于單元數(shù)目眾多，因此采用Python腳本程序自動生成[3]。側(cè)向緩沖墊和底部橡膠墊板在軟件中利用點對面的形式直接設(shè)置生成。扣件系統(tǒng)與減振結(jié)構(gòu)參數(shù)如表2所示[8]。本文主要研究梯形軌枕支承小半徑曲線軌道的固有特性，有限元模型中不考慮支撐塊以簡化計算。

表2 梯形軌枕軌道支承結(jié)構(gòu)參數(shù)

1.3 復特征值法

通過對梯形軌枕支承小半徑曲線輪軌系統(tǒng)發(fā)生摩擦自激振動時的振動頻率和相應(yīng)的振動模態(tài)的提取，預測系統(tǒng)穩(wěn)定性的方法稱為復特征值法[3]。建立不考慮摩擦和外力作用影響的運動微分方程：

(1)

式中：

x——節(jié)點位移向量；

M、C、K——分別為質(zhì)量、阻尼、剛度矩陣,且均為對稱矩陣。

根據(jù)式(1)，可得到通解：

(2)

式中：

yk——第k階特征向量；

αk+jwk——第k階復特征值，其中，αk為特征值的實部，wk為特征值的虛部；

t——時間。

當αk>0時，輪軌系統(tǒng)可能會發(fā)生摩擦自激振動。

同時，通過引入評判摩擦自激振動發(fā)生趨勢的等效阻尼比ξk(ξk=-α/(π|wk|))亦可。當ξk為負值時，說明輪軌系統(tǒng)可能會發(fā)生摩擦自激振動，等效阻尼比越大，其發(fā)生的概率越小。

1.4 輪軌摩擦自激振動誘發(fā)波磨理論

通常鐵路領(lǐng)域認為輪對與鋼軌間的摩擦功波動是引發(fā)波磨病害的主要原因。文獻[9]提出的磨損公式如下：

w=K(H-C0)

(3)

式中：

w——單位時間內(nèi)鋼軌頂面的磨損量；

K——磨損系數(shù)；

H——摩擦功(H為飽和蠕滑力F與相對速度v0的乘積)；

C0——耐久性摩擦功。

在小半徑曲線軌道上F將達到飽和狀態(tài)時，F(xiàn)為：

F=μN

(4)

式中：

μ——輪軌間的動摩擦系數(shù)；

N——輪軌間的法向接觸力。

根據(jù)式(4)可知，當輪軌系統(tǒng)發(fā)生摩擦自激振動時，輪軌間的摩擦功也出現(xiàn)了同頻率的振動，進而導致波磨病害的產(chǎn)生[3]。

1.5 輪軌系統(tǒng)的復特征值分析

根據(jù)前期調(diào)查和現(xiàn)場測試，在曲線半徑為300 m的梯形軌枕支承小半徑曲線軌道上，鋼軌波磨的情況十分嚴重，其波長主要集中在60～80 mm[10]，且主要發(fā)生在低軌上，而通過此路段地鐵車輛的運營速度約為80 km/h。根據(jù)頻率f、速度v和波長λ的關(guān)系式f=v/λ，可以求得波磨的頻率為277.8～463.0 Hz。

本文針對梯形軌枕支承小半徑曲線輪軌系統(tǒng)有限元模型進行復特征值分析，提取相應(yīng)的摩擦自激振動模態(tài)和頻率，如圖4所示。

圖4 梯形軌枕支承小半徑曲線輪軌系統(tǒng)的等效阻尼比分布圖

由圖4可見，ξk分別為-0.013 430、-0.009 887、-0.023 590，其對應(yīng)的頻率分別為156 Hz、230 Hz、337 Hz，說明梯形軌枕支承小半徑曲線輪軌系統(tǒng)存在3個摩擦自激振動頻率。其中，ξk是用于評判輪軌摩擦自激振動的重要參數(shù)。當梯形軌枕支承小半徑曲線輪軌系統(tǒng)的等效阻尼比為負數(shù)時，系統(tǒng)才可能會發(fā)生摩擦自激振動，且數(shù)值越小，系統(tǒng)發(fā)生摩擦自激振動的趨勢越大。因此可認為當ξk=-0.023 590時，此時輪軌系統(tǒng)發(fā)生摩擦自激振動的情況最容易出現(xiàn)，其對應(yīng)的摩擦自激振動頻率為337 Hz。該頻率下對應(yīng)的振動模態(tài)如圖5所示。由圖5可知，輪軌系統(tǒng)摩擦自激振動發(fā)生在低軌上。綜上所述，當輪軌系統(tǒng)振動頻率為337 Hz時，其會產(chǎn)生摩擦自激振動且發(fā)生在低軌上，從而可能誘導波長為66 mm的波磨，這與實際測試結(jié)果相符，進一步證明了誘導梯形軌枕支承小半徑曲線軌道上波磨病害的成因可能是輪軌間的摩擦自激振動。

圖5 梯形軌枕支承小半徑曲線輪軌系統(tǒng)的摩擦自激振動模態(tài)

2 軌道減振結(jié)構(gòu)參數(shù)和總的扣件參數(shù)對輪軌摩擦自激振動的影響

2.1 軌道減振結(jié)構(gòu)參數(shù)對輪軌摩擦自激振動的影響

在梯形軌枕支承小半徑曲線軌道支承結(jié)構(gòu)中，梯形軌枕與支撐塊之間的底部橡膠墊板和側(cè)向緩沖墊是該軌道支承結(jié)構(gòu)中減振降噪的主要部件。因此，軌道減振結(jié)構(gòu)對于輪軌系統(tǒng)的振動具有重要影響[10]。采用參數(shù)化分析研究梯形軌枕支承小半徑曲線軌道減振結(jié)構(gòu)中側(cè)向緩沖墊與底部橡膠墊板的結(jié)構(gòu)參數(shù)呈倍數(shù)變化時對輪軌摩擦自激振動的影響。減振型軌道中減振結(jié)構(gòu)的靜態(tài)剛度一般在8×106～4×107N/m范圍內(nèi)，靜態(tài)阻尼一般在1.5×105Ns/m內(nèi)[4]。在此范圍內(nèi)，側(cè)向緩沖墊和底部橡膠墊板結(jié)構(gòu)參數(shù)的變化情況如表3所示[8]。文獻[8]發(fā)現(xiàn)側(cè)向緩沖墊中縱向、垂向參數(shù)以及底部橡膠墊板中縱向、橫向參數(shù)對波磨病害幾乎沒有影響。

表3 減振結(jié)構(gòu)參數(shù)隨其倍數(shù)變化值

利用復特征值法研究側(cè)向緩沖墊橫向剛度與阻尼、底部橡膠墊板垂向剛度與阻尼同時呈倍數(shù)變化時對輪軌摩擦自激振動頻率的影響，如圖6所示。由圖6可知，當側(cè)向緩沖墊和底部橡膠墊板的結(jié)構(gòu)參數(shù)同時呈倍數(shù)變化時，輪軌系統(tǒng)的摩擦自激振動頻率始終保持在337 Hz。綜上可知，側(cè)向緩沖墊橫向剛度與阻尼、底部橡膠墊板垂向剛度與阻尼同時呈倍數(shù)變化時，對輪軌系統(tǒng)的摩擦自激振動頻率基本沒有影響。

圖6 軌道減振結(jié)構(gòu)參數(shù)-輪軌摩擦自激振動頻率關(guān)系曲線

圖7為側(cè)向緩沖墊的橫向剛度、橫向阻尼呈倍數(shù)變化時對輪軌摩擦自激振動發(fā)生趨勢的影響。由圖7 a)可知，隨著側(cè)向緩沖墊橫向剛度的增加，等效阻尼比變化幅度不大，說明側(cè)向緩沖墊橫向剛度對輪軌摩擦自激振動發(fā)生趨勢影響很小。由圖7 b)可知，隨著側(cè)向緩沖墊橫向阻尼的增加，等效阻尼比基本不變，說明側(cè)向緩沖墊橫向阻尼的變化對輪軌摩擦自激振動發(fā)生趨勢的影響也很小。

a) 側(cè)向緩沖墊橫向剛度

b) 側(cè)向緩沖墊橫向阻尼圖7 側(cè)向緩沖墊參數(shù)-等效阻尼比關(guān)系曲線Fig.7 Relationship curve of parameters of the lateral cushion and equivalent damping ratio

圖8為底部橡膠墊板的垂向剛度、垂向阻尼呈倍數(shù)變化時對輪軌摩擦自激振動發(fā)生趨勢的影響。由圖8 a)可見，隨著底部橡膠墊板垂向剛度的增加，摩擦自激振動的等效阻尼比呈增大趨勢，但增幅極小。由圖8 b)可見，隨著底部橡膠墊板垂向阻尼的增大，等效阻尼比越大，即輪軌摩擦自激振動的發(fā)生趨勢就越小。綜上所述，軌道減振結(jié)構(gòu)中側(cè)向緩沖墊參數(shù)對輪軌摩擦自激振動的發(fā)生趨勢影響不大，所以通過增大其橫向剛度降低摩擦自激振動發(fā)生趨勢的意義不大。而底部橡膠墊板中垂向阻尼對摩擦自激振動發(fā)生趨勢有較大的影響，在一定范圍內(nèi)通過增加底部橡膠墊板的垂向阻尼，有助于降低輪軌摩擦自激振動的發(fā)生趨勢，當?shù)撞肯鹉z墊板中垂向阻尼接近于1.5×105Ns/m時，輪軌系統(tǒng)最不易發(fā)生摩擦自激振動。

a) 底部橡膠墊板垂向剛度

b) 底部橡膠墊板垂向阻尼圖8 底部橡膠墊板參數(shù)-等效阻尼比關(guān)系曲線Fig.8 Relationship curve of parameters of the bottom rubber pad and equivalent damping ratio

2.2 總的扣件參數(shù)對輪軌摩擦自激振動的影響

文獻[11-12]驗證了軌道支承結(jié)構(gòu)中扣件參數(shù)對波磨病害的影響較為明顯。本文采用參數(shù)化分析討論總的扣件參數(shù)呈倍數(shù)變化時對輪軌摩擦自激振動發(fā)生趨勢的影響。

現(xiàn)階段，國內(nèi)城市軌道結(jié)構(gòu)使用的扣件靜態(tài)剛度一般在5～60 MN/m，靜態(tài)阻尼一般小于1.5×105Ns /m[3]。在此范圍內(nèi)，扣件參數(shù)的變化如表4所示[11]。采用復特征值法研究不同扣件參數(shù)對梯形軌枕輪軌摩擦自激振動發(fā)生趨勢的影響，如圖9～10所示。

表4 扣件參數(shù)隨其倍數(shù)變化值

由圖9可見，隨著總的扣件剛度的增大，即扣件的垂向、橫向、縱向剛度均增大時，輪軌系統(tǒng)的負等效阻尼比增大，即發(fā)生摩擦自激振動的趨勢降低。因此，通過增加總的扣件剛度有助于降低輪軌系統(tǒng)摩擦自激振動的發(fā)生趨勢，即當總的扣件剛度接近60 MN/m時，輪軌系統(tǒng)最不易發(fā)生摩擦自激振動。

圖9 總的扣件剛度-等效阻尼比關(guān)系曲線Fig.9 Relationship curve of overall fastener stiffness and equivalent damping ratio

圖10 總的扣件阻尼-等效阻尼比關(guān)系曲線Fig.10 Relationship curve of overall fastener damping and equivalent damping ratio

由圖10可知，隨著總的扣件阻尼(即扣件垂向、橫向、縱向阻尼)的增大，輪軌系統(tǒng)負的等效阻尼比基本不變，即輪軌系統(tǒng)發(fā)生摩擦自激振動的趨勢基本不變。因此，通過調(diào)整總的扣件阻尼以降低摩擦自激振動的發(fā)生趨勢意義不大?；谀Σ磷约ふ駝拥挠^點，在一定范圍內(nèi)，通過改變總的扣件剛度對產(chǎn)生自激振動的等效阻尼比影響較大，而總的扣件阻尼變化對等效阻尼比的影響較小。由此可見，通過增加總的扣件剛度有助于緩解輪軌系統(tǒng)摩擦自激振動的發(fā)生趨勢。

3 結(jié)論

1) 在梯形軌枕支承小半徑曲線軌道上，使輪軌系統(tǒng)發(fā)生摩擦自激振動的頻率為337 Hz且主要發(fā)生在低軌上，在該頻率下可能會誘發(fā)波長為66 mm的波磨。

2) 在軌道減振結(jié)構(gòu)參數(shù)對輪軌摩擦自激振動的影響研究中，減振結(jié)構(gòu)參數(shù)的改變對該區(qū)段輪軌摩擦自激振動的誘導頻率的影響較小，而減振結(jié)構(gòu)參數(shù)對輪軌摩擦自激振動的發(fā)生趨勢的影響較為明顯。在一定范圍內(nèi)通過增加減振結(jié)構(gòu)中底部橡膠墊板的垂向阻尼，可以有效降低摩擦自激振動發(fā)生的趨勢，從而有助于減緩相應(yīng)區(qū)段的波磨病害。

3) 在一定范圍內(nèi)通過改變扣件參數(shù)，發(fā)現(xiàn)總的扣件剛度對輪軌系統(tǒng)摩擦自激振動的發(fā)生趨勢影響較大，而總的扣件阻尼對輪軌摩擦自激振動的發(fā)生趨勢影響較小。通過增加總的扣件剛度有助于減緩相應(yīng)區(qū)段內(nèi)的波磨病害。