高速列車—無砟軌道—橋梁耦合系統(tǒng)振動特性分析

利遠(yuǎn)翔,雷曉燕,張 斌

(華東交通大學(xué)鐵路環(huán)境振動與噪聲教育部工程研究中心,江西南昌330013)

國內(nèi)外對于車輛—軌道系統(tǒng)動力學(xué)及車輛—橋梁系統(tǒng)動力學(xué)的研究已經(jīng)相當(dāng)深入,但是由于各自不同的目的,將車輛、軌道和橋梁作為一個整體系統(tǒng)加以研究的卻不多見。從上世紀(jì)80年代起,國內(nèi)外學(xué)者對快速列車和重載列車引起的地面振動及對周圍環(huán)境的影響進(jìn)行了較多研究,提出了一些計算方法和計算模型。如在1994年,日本的Fujikake[1]和松浦章夫[2]分別就交通車輛引起的結(jié)構(gòu)振動發(fā)生機(jī)理、振動波在地下和地面的傳播規(guī)律及其對周圍居民的影響進(jìn)行了研究,提出了環(huán)境振動水平的預(yù)測方法。日本鐵道技術(shù)研究所還就新干線高速列車對環(huán)境振動的影響進(jìn)行了現(xiàn)場測試,分析了車輛、軌道、橋梁等不同部位的振動特點(diǎn),并研究了不同類型屏障的隔振效果。而對于高架軌道系統(tǒng)誘發(fā)環(huán)境振動特性的研究,源自于1999年開始動工修建的臺灣高速鐵路,經(jīng)過臺南高科技園區(qū)的一段高架軌道振動超標(biāo),從而引起了國內(nèi)外學(xué)者對高架軌道誘發(fā)環(huán)境振動的關(guān)注并開始研究。臺灣成功大學(xué)教授朱圣浩等[3-4]采用有限元的方法,在建立車輛—橋梁—大地耦合振動的三維有限元模型時,將列車以動輪單元模擬,大地采用三維塊體單元模擬,分析了高架軌道引起的周邊建筑物的振動。本文建立了車輛單元模型、軌道—橋梁單元模型,推導(dǎo)了車輛—軌道—橋梁耦合系統(tǒng)的振動方程,并在Matlab軟件中編制了程序,最后通過一個實(shí)例分析了高速列車—無砟軌道—橋梁耦合系統(tǒng)垂向振動特性。

1 車輛單元模型

本文的模型是在雷曉燕教授的車輛—軌道—地基耦合系統(tǒng)模型[5-7]的基礎(chǔ)上建立的。圖1為車輛單元模型,在整車模型中,車體和轉(zhuǎn)向架考慮沉浮振動和點(diǎn)頭振動,車輪考慮沉浮振動,每節(jié)鋼軌上兩個節(jié)點(diǎn)考慮豎向位移和轉(zhuǎn)角,整車模型有26個自由度。

圖1中,Mc,Jc為車體的質(zhì)量與轉(zhuǎn)動慣量;Mt,Jt為轉(zhuǎn)向架的質(zhì)量與轉(zhuǎn)動慣量;ks1,ks2為車輛一、二系懸掛剛度;cs1,cs2為車輛一、二系懸掛阻尼;Mwi(i=1,2,3,4)表示第 i個車輪的質(zhì)量;kc為輪軌間接觸剛度;vc,θc為車體沉浮振動的豎向位移、車體點(diǎn)頭振動的角位移;vi,θi(i=t1,t2)表示前、后轉(zhuǎn)向架沉浮振動的豎向位移、點(diǎn)頭振動的角位移;vi(i=w1,w2,w3,w4)為第i個車輪的豎向位移;vci(i=1,2,3,4)為第i個輪軌接觸處鋼軌的豎向位移。

圖1 車輛單元模型

從圖1中我們可以定義車輛單元的節(jié)點(diǎn)位移向量為

為了建立車輛單元的有限元方程,可以運(yùn)用Hamilton原理,即

其中:L為Lagrange函數(shù),L=T-Π;T為動能;Π為勢能;R為耗散能。

通過列出車輛單元的彈性勢能、動能和耗散能可推導(dǎo)出車輛單元的剛度矩陣、質(zhì)量矩陣和阻尼矩陣

2 軌道—橋梁單元模型

四層軌道—橋梁單元模型如圖2所示,將軌道結(jié)構(gòu)和橋梁簡化為四層梁模型,從上到下分別為鋼軌、軌道板、混凝土支承層和橋梁。模型僅考慮沉浮和點(diǎn)頭振動,共有16個自由度。圖2中v1,v5表示鋼軌的豎向位移;θ1,θ5表示鋼軌的轉(zhuǎn)角;軌下墊層的支承彈性系數(shù)和阻尼系數(shù)分別用ksy1和csy1表示;v2,v6表示軌道板的豎向位移;θ2,θ6表示預(yù)制軌道板的轉(zhuǎn)角;軌道板下瀝青水泥砂漿層的支承彈性系數(shù)和阻尼系數(shù)分別用ksy2和csy2表示;v3,v7表示混凝土板支承層的豎向位移;θ3,θ7表示混凝土板支承層的轉(zhuǎn)角;底座混凝土板支承層下橋梁的支承彈性系數(shù)和阻尼系數(shù)分別用ksy3和csy3表示;ν4,ν8表示橋梁單元的豎向位移;θ4,θ8表示橋梁單元的轉(zhuǎn)角。

從圖2中我們可以定義軌道—橋梁單元的節(jié)點(diǎn)位移向量為

同樣為了建立軌道—橋梁單元的有限元方程,利用Hamilton原理,通過列出軌道—橋梁單元的彈性勢能、動能和耗散能可推導(dǎo)軌道—橋梁單元的剛度矩陣e、質(zhì)量矩陣e和阻尼矩陣e。

圖2 板式軌道—橋梁單元模型

2.1 軌道—橋梁單元的彈性勢能

軌道單元的勢能包括鋼軌的彎曲勢能,預(yù)制軌道板的彎曲勢能,底座混凝土板支承層的彎曲勢能,橋梁的彎曲勢能,離散支承彈簧的勢能以及連續(xù)支承彈簧的勢能。

(1)鋼軌的彎曲勢能

鋼軌上任意一點(diǎn)的位移可表示為

其中位移插值函數(shù)為

其中

2.2 軌道—橋梁單元的動能

軌道單元的動能包括鋼軌的彎曲動能,預(yù)制軌道板的彎曲動能,底座混凝土板支承層的彎曲動能以及橋梁單元的彎曲動能。

(1)鋼軌的彎曲動能

其中

同理可以推導(dǎo)出預(yù)制軌道板單元的質(zhì)量矩陣、底座混凝土板支承層單元的質(zhì)量矩陣和橋梁單元的質(zhì)量矩陣

2.3 軌道—橋梁單元的耗散能

(1)比例阻尼

(2)離散支承彈簧的耗散能

2.4 軌道—橋梁單元的運(yùn)動方程

得板式軌道—橋梁單元四層梁的有限元方程如下

其中

3 車輛—軌道—橋梁耦合系統(tǒng)有限元方程

考慮車輛—無砟軌道—橋梁耦合系統(tǒng)的有限元方程包含兩種單元,即無砟軌道—橋梁單元和車輛單元,相應(yīng)的單元剛度矩陣、質(zhì)量矩陣和阻尼矩陣分別為Ku,Mu,Cu見(3)和KT,MT,CT見(19)。利用標(biāo)準(zhǔn)有限元的集成方法可形成車輛—無砟軌道—橋梁耦合系統(tǒng)的總剛度矩陣、總質(zhì)量矩陣和總阻尼矩陣。最后得到車輛—無砟軌道—橋梁耦合系統(tǒng)的動力有限元方程如下:

其中

考慮橋梁支承處的邊界條件,即橋墩沿豎直方向的位移為零,車輛—無砟軌道—橋梁耦合系統(tǒng)動力有限元方程數(shù)值解的實(shí)現(xiàn),可通過直接積分法,如Newmark數(shù)值積分法。

4 車輛—軌道—橋梁耦合系統(tǒng)算例以及分析

對上述車輛—無砟軌道—橋梁耦合系統(tǒng)有限元模型運(yùn)用Matlab編制了計算程序,并進(jìn)行了實(shí)例計算?，F(xiàn)利用編制的程序分析高速列車—無砟軌道—橋梁耦合系統(tǒng)垂向振動特性。考慮8跨連續(xù)梁的情況,橋梁跨度32 m,整個橋長256 m,列車時速300 km?h-1,混凝土支承塊間距0.65 m,時間步長0.001 s。車輛、軌道和橋梁結(jié)構(gòu)的參數(shù)見表1和表2。為簡化計算,不計軌道隨機(jī)不平順。

通過仿真分析可以得出輪軌接觸力、鋼軌豎向位移、橋梁豎向位移等時程曲線圖,分別見圖3～圖7。

表1 和諧號高速動車CRH3車輛結(jié)構(gòu)參數(shù)

表2 CRTSⅡ型板式無砟軌道和橋梁結(jié)構(gòu)參數(shù)

從時程曲線圖可以看出:輪軌接觸力保持在70 kN左右,從鋼軌豎向位移時程曲線圖中可以明顯的看出四個輪對與鋼軌接觸時的四個峰值,而且在鋼軌豎向加速度時程曲線圖中也可以看出四個峰值。還可看出鋼軌位移時程曲線與橋梁時程曲線圖相比較,前者的振動位移比后者的大,在軌道—橋梁單元模型中的四層梁中,從上到下位移呈遞減的趨勢。鋼軌和橋梁的振動加速度也呈現(xiàn)相同的規(guī)律。

圖3 輪軌接觸力時程曲線

圖4 鋼軌豎向位移時程曲線

圖5 橋梁豎向位移時程曲線

圖6 鋼軌豎向加速度時程曲線圖

5 結(jié)論

通過建立車輛單元模型和軌道—橋梁四層梁單元模型,利用拉格朗日方程推導(dǎo)出整個耦合系統(tǒng)的剛度、質(zhì)量和阻尼矩陣,最后得出耦合系統(tǒng)的振動方程。相對于其它復(fù)雜的車輛—軌道—橋梁耦合動力學(xué)模型,該方法思路清晰,借助于Matlab軟件使得程序編制容易實(shí)現(xiàn),計算時效快且能反映軌道—橋梁振動的基本規(guī)律和特點(diǎn)。

圖7 橋梁豎向加速度時程曲線圖

[1]FUJIKAKE T.A predictionmethod for the propagation of groundvibration from railway trains[J].Journalof Sound and Vibration,1986,111(2):357-360.

[2]松浦章夫.高速鐵路車輛與橋梁相互作用[J].鐵道技術(shù)研究資料,1974,31(5):14-17.

[3]JU SHENHAW,LIN HUNGTA.Analysis of train-induced vibrations and vibration reduction schemes above and below critical Rayleigh speeds by finite element method[J].Soil Dynamics and Earthquake Engineering,2004,24(4):993-1 002.

[4]JU SHENHAW,LIN HUNGTA.A simple finite element model for vibration analyses induced by moving vehicles[J].Soil Dynamics and Earthquake Engineering,2006,68(12):1 232-1 256.

[5]雷曉燕現(xiàn)代鐵路理論研究[M].北京:中國鐵道出版社,2006.

[6]雷曉燕,張斌,劉慶杰.列車—軌道系統(tǒng)豎向動力分析的車輛軌道單元模型[J].振動與沖擊,2009,29(3):168-173.

[7]雷曉燕,張斌,劉慶杰.軌道過渡段動力特性的有限元分析[J].中國鐵道科學(xué),2009,30(5):15-21.