基于傳遞矩陣法的車輛-軌道垂向耦合系統(tǒng)動力分析方法

朱志輝， 王盈瑩， 龔 威， 馮乾朔

(1.中南大學(xué) 土木工程學(xué)院，長沙 410075； 2.高速鐵路建造技術(shù)國家工程實驗室，長沙 410075)

鐵路車輛-軌道耦合系統(tǒng)由移動車輛和連續(xù)軌道結(jié)構(gòu)組成。列車沿軌道高速運行將引起軌道結(jié)構(gòu)振動，強(qiáng)烈的振動不僅會直接影響軌道的工作狀態(tài)和使用性能，還會降低列車的行車平穩(wěn)性和安全性，影響乘客乘坐舒適度[1]。因此，準(zhǔn)確和有效地評估車輛和軌道間的相互作用是至關(guān)重要的。

建立準(zhǔn)確的軌道結(jié)構(gòu)動力學(xué)模型是車輛-軌道耦合動力學(xué)理論研究的重要部分。有限元直接剛度法(direct stiffness method, DSM)可直接基于軌道結(jié)構(gòu)有限元模型建立其動力方程[2-4]。這種方法可準(zhǔn)確計算軌道結(jié)構(gòu)的局部與高頻振動[5]，但軌道模型的自由度數(shù)目會隨著軌道長度的增長而增加，當(dāng)涉及較長軌道結(jié)構(gòu)動力分析時，不可避免會產(chǎn)生大量結(jié)構(gòu)自由度，使得動力方程系數(shù)矩陣規(guī)模龐大，導(dǎo)致建模和計算的困難[6]。為了克服這個問題，翟婉明[7]假設(shè)車輛相對于軌道結(jié)構(gòu)只在一個固定點振動，提出通過加載車輪上的軌道隨機(jī)不平順性得到車輛響應(yīng)的近似辦法，具有較高的計算效率，但忽略了軌道結(jié)構(gòu)變化的影響。Guo等[8-10]采用模態(tài)疊加法建立軌道結(jié)構(gòu)動力方程，縮減了軌道自由度數(shù)目，但由于軌道結(jié)構(gòu)的局部高頻振動特性，導(dǎo)致其主導(dǎo)振型難以選取[11]。移動軌道技術(shù)[12-13]通過刪除非荷載作用處軌道結(jié)構(gòu)有效降低了自由度數(shù)目，但新增軌道連接處存在位移不協(xié)調(diào)問題，新增單元處會出現(xiàn)動力響應(yīng)突變[14]。

傳遞矩陣法(transfer matrix method，TMM)的基本思想是將周期性結(jié)構(gòu)分割為若干元胞結(jié)構(gòu)，將元胞的力學(xué)特性用傳遞矩陣表示，建立元胞結(jié)構(gòu)狀態(tài)矢量傳遞模型，實現(xiàn)整個周期性結(jié)構(gòu)的求解。傳遞矩陣法基于元胞結(jié)構(gòu)進(jìn)行分析，無需建立結(jié)構(gòu)整體動力方程，避免了大型線性方程組的求解[15-16]，降低了計算過程中的內(nèi)存占用。Mead等[17-21]采用傳遞矩陣法分別針對周期性的彈性支承梁、框架結(jié)構(gòu)、薄壁構(gòu)件、拱橋和鋼橋進(jìn)行了研究，取得了良好的效果。對于鏈?zhǔn)较到y(tǒng)，狀態(tài)矢量的傳遞僅相當(dāng)于元胞單元傳遞矩陣相乘，傳遞矩陣連續(xù)相乘即可得到系統(tǒng)總傳遞矩陣。然而，由于單元傳遞矩陣連續(xù)相乘，累積誤差被加劇，可能產(chǎn)生數(shù)值不穩(wěn)定現(xiàn)象[22]。針對該問題，Horner等[23]結(jié)合Riccati變換與傳遞矩陣法，提出了一種Riccati傳遞矩陣法，通過傳遞矩陣的遞推建立系統(tǒng)傳遞關(guān)系，有效提高了傳遞矩陣法的數(shù)值穩(wěn)定性。Xue[24]在元胞結(jié)合面處引入剛度方程假定，將傳遞矩陣連續(xù)相乘演變?yōu)樵Y(jié)合面剛度方程的遞推傳遞，有效減小了因矩陣連續(xù)相乘產(chǎn)生的累積誤差。

軌道結(jié)構(gòu)具有明顯的周期特性，Ma等[25-26]利用這一特性，有效地解決了軌道結(jié)構(gòu)的頻域響應(yīng)分析問題，但該類方法無法考慮軌道的非線性特征和軌道構(gòu)成差異性。因此，本文利用傳遞矩陣方法建立軌道結(jié)構(gòu)傳遞模型并求解其動力響應(yīng)，從而實現(xiàn)不增加軌道系統(tǒng)計算自由度數(shù)目而可以增長軌道結(jié)構(gòu)計算長度的目的。同時，當(dāng)軌道結(jié)構(gòu)縱向構(gòu)成有差異，如出現(xiàn)軌枕空吊、扣件失效等情況時，只需將該處元胞更換為包含差異構(gòu)成的元胞，不影響整體傳遞關(guān)系的建立。

本文將車輛-軌道耦合系統(tǒng)分解為車輛子系統(tǒng)和軌道子系統(tǒng)，兩子系統(tǒng)間通過輪軌相互作用力實現(xiàn)耦合。其中，車輛子系統(tǒng)采用多剛體動力學(xué)方法建立動力方程；對于具有顯著周期性特征的軌道子系統(tǒng)，以有砟軌道和CRTSⅡ型板式無砟軌道為例，基于TMM將軌道結(jié)構(gòu)的周期性重復(fù)部分劃分為若干類元胞結(jié)構(gòu)，建立元胞結(jié)構(gòu)的動力方程，并基于引入剛度方程假定的Riccati傳遞矩陣法，確定元胞內(nèi)部及相鄰元胞間的傳遞關(guān)系，從而建立周期性軌道結(jié)構(gòu)傳遞矩陣模型。通過對比TMM與DSM計算的CRH2型動車通過無砟軌道結(jié)構(gòu)引起的車輛和軌道結(jié)構(gòu)動力響應(yīng)結(jié)果，驗證了本文方法的可行性和高效性。

1 軌道結(jié)構(gòu)傳遞矩陣模型

根據(jù)鋼軌下部有無道砟支撐層，現(xiàn)有軌道結(jié)構(gòu)可分為有砟和無砟軌道兩類，其結(jié)構(gòu)形式和力學(xué)特征均有所不同，因此采用傳遞矩陣法對其進(jìn)行動力分析時，應(yīng)針對不同類型的軌道結(jié)構(gòu)，提出相應(yīng)的元胞劃分方案。CRTSⅡ型板式無砟軌道是我國廣泛應(yīng)用的無砟軌道結(jié)構(gòu)形式之一，京津城際、京滬高鐵等均采用了此種軌道結(jié)構(gòu)[27]。故本文以有砟軌道和CRTSⅡ型縱連板式無砟軌道為例，提出元胞劃分方案，推導(dǎo)元胞間的動力傳遞關(guān)系，建立基于傳遞矩陣法的軌道結(jié)構(gòu)元胞模型。軌道系統(tǒng)左右對稱，為簡便起見，取半結(jié)構(gòu)進(jìn)行介紹，且只考慮在豎向平面內(nèi)的動力響應(yīng)。

1.1 元胞結(jié)構(gòu)的劃分

軌道結(jié)構(gòu)由重復(fù)的單跨軌道構(gòu)成，因此應(yīng)用傳遞矩陣法時，僅需對單跨軌道進(jìn)行劃分。為便于建立元胞內(nèi)部傳遞關(guān)系，元胞結(jié)構(gòu)的劃分應(yīng)遵循僅元胞端部存在自由度，且相鄰元胞交界面自由度數(shù)目相同的原則。

有砟軌道垂向整體模型如圖1(a)所示，其中鋼軌視作連續(xù)彈性離散點支撐的垂向Timoshenko梁模型，考慮豎向位移Zr及轉(zhuǎn)動角θr，長度為L；軌枕(及道砟質(zhì)量)等效為集中質(zhì)量ms，考慮豎向位移Zs，相鄰兩軌枕間距為l；扣件和道砟等效為彈簧-阻尼器，考慮豎向剛度系數(shù)kf，kb及阻尼系數(shù)cf，cb。為對鋼軌進(jìn)行較為精細(xì)的劃分，如圖1(b)所示，將有砟軌道單跨軌道劃分為3類元胞結(jié)構(gòu)：元胞A為左側(cè)支撐結(jié)構(gòu)，由鋼軌、軌枕、扣件和道砟組成；元胞B為中間連接結(jié)構(gòu)，由鋼軌組成；元胞C為右側(cè)支撐結(jié)構(gòu)，構(gòu)成與元胞A相同，由于相鄰兩元胞包含同一支承結(jié)構(gòu)，故二者支承結(jié)構(gòu)參數(shù)均為整體的1/2(端部元胞除外)。

CRTSⅡ型板式無砟軌道如圖2(a)所示，鋼軌、軌道板和底座板被視為連續(xù)彈性離散點支撐的Timoshenko梁模型，長度為L，考慮豎向位移Zr、Zt、Zc及轉(zhuǎn)動角θr、θt、θc,下標(biāo)r、t、c分別表示鋼軌、軌道板、底座板；扣件、CA砂漿和基礎(chǔ)等效為周期性離散分布的黏滯阻尼和線性彈簧，間距為l，考慮豎向剛度系數(shù)kf、kc、ks和阻尼系數(shù)cf、cc、cs，下標(biāo)f、c、s分別表示扣件、CA砂漿、基礎(chǔ)。類似于有砟軌道，將CRTSⅡ型板式無砟軌道單跨軌道結(jié)構(gòu)劃分為3類元胞結(jié)構(gòu)，如圖2(b)所示。其中，元胞A和C均由鋼軌、軌道板、底座板、扣件、CA砂漿和地基組成，由于相鄰兩元胞包含同一支承結(jié)構(gòu)，故二者支承結(jié)構(gòu)參數(shù)均為整體模型的1/2(端部元胞除外)；元胞B由鋼軌、軌道板和底座板組成。

(a) 整體垂向模型

(b) 單跨軌道劃分示意圖圖1 有砟軌道垂向模型劃分示意圖Fig.1 Structure division of ballast track

(a) 整體垂向模型

(b) 單跨軌道劃分示意圖圖2 CRTSⅡ型板式無砟軌道垂向模型劃分示意圖Fig.2 Structure division of CRTSⅡ slab ballastless track

1.2 元胞結(jié)構(gòu)動力方程

任意取整體結(jié)構(gòu)中第i個軌道元胞，其運動方程可表示為

(1)

采用Newmark-β直接積分法[29]對元胞動力方程進(jìn)行求解，式(1)可寫為

(2)

圖3 元胞i輸入端到輸出端傳遞關(guān)系示意圖Fig.3 The transmission of vectors from the input to the output of cell i

對于元胞i

(3)

如圖2所示，元胞i輸入和輸出端自由度對應(yīng)的位移和內(nèi)力可表示為

(4)

式中:Z、θ、Q和M分別表示位移、轉(zhuǎn)角、剪力和彎矩;下標(biāo)r、t和c分別表示鋼軌、軌道板和底座板。

(5)

1.3 傳遞模型的建立

由位移和內(nèi)力的連續(xù)條件[31]，元胞結(jié)合面處的狀態(tài)矢量傳遞關(guān)系可表示為

Ui-1,O=Ui,I,Ni-1,O=-Ni,I

(6)

引入類似于狀態(tài)矢量的廣義Riccita變換，假設(shè)元胞i左側(cè)結(jié)合面處內(nèi)力向量與位移向量的廣義剛度方程為

Ni,I=SiUi,I+Ei(i≥2)

(7)

式中:Si是結(jié)合面處剛度方程的系數(shù)矩陣;Ei是結(jié)合面處的等效外力向量。

展開式(5)，可得：

(8)

(9)

將式(6)和式(7)代入式(8)，可得由第i個元胞輸出端位移響應(yīng)求解第i-1個元胞輸出端響應(yīng)的傳遞關(guān)系式

(10)

將式(10)代入式(9)，可得元胞i+1輸入端內(nèi)力與位移向量的廣義剛度方程為

Ni+1,I=Si+1Ui+1,I+Ei+1

(11)

其中：

(12)

(13)

式(12)和式(13)即為相鄰剖面間，剛度方程的系數(shù)矩陣和等效外荷載向量的傳遞關(guān)系。

式(11)中，令i=2，有：

N2,I=S2U2,I+E2

(14)

將式(6)代入式(14)，則上式可表示為

N1,O=-S2U1,O-E2

(15)

對于元胞1，擴(kuò)展式(5)可得：

(16)

(17)

軌道結(jié)構(gòu)兩端無內(nèi)力作用，N1,I是已知的，因此由式(16)可得元胞1的輸入端位移為

(18)

將式(18)代入式(17)可得元胞1的輸出端內(nèi)力為

(19)

對照式(15)可知

(20)

(21)

將S2和E2代入式(12)和式(13)，可完成Si和Ei在整個結(jié)構(gòu)中由輸入端至輸出端作剛度方程的連續(xù)傳遞，進(jìn)而建立整體結(jié)構(gòu)的狀態(tài)矢量傳遞模型?？紤]結(jié)構(gòu)輸出端邊界條件，通過式(7)可得結(jié)構(gòu)輸出端的位移和內(nèi)力響應(yīng)，將其代入輸出端響應(yīng)傳遞關(guān)系式(10)，由輸出端至輸入端依次求解，可得各元胞的位移響應(yīng)，并通過Newmark-β方法解得其速度和加速度響應(yīng)。

由上述推導(dǎo)過程可知，采用傳遞矩陣法對軌道子系統(tǒng)進(jìn)行動力分析時，無需建立軌道結(jié)構(gòu)整體模型，僅需拆分其單跨軌道為若干元胞結(jié)構(gòu)，并依據(jù)元胞結(jié)構(gòu)間的傳遞關(guān)系，建立軌道結(jié)構(gòu)的整體傳遞模型即可。對于垂向有砟/無砟軌道模型，無論軌道整體長度如何變化，計算過程中所涉及的最大矩陣階數(shù)始終為3/6。

2 車輛-軌道耦合系統(tǒng)動力分析模型

本文在車輛-軌道動力分析模型中，將整個耦合系統(tǒng)劃分為車輛子系統(tǒng)和軌道子系統(tǒng)，通過輪軌間的相互作用力來建立兩個子系統(tǒng)的聯(lián)系。其中，軌道子系統(tǒng)模型通過軌道結(jié)構(gòu)傳遞矩陣模型建立。

2.1 車輛子系統(tǒng)

車輛視作多剛體模型，單節(jié)車輛由1個車體、2個轉(zhuǎn)向架、4個輪對、4個一系懸掛系統(tǒng)和2個二系懸掛系統(tǒng)組成，如圖4所示。車體、轉(zhuǎn)向架和輪對采用剛體模擬，懸掛系統(tǒng)通過線性彈簧-阻尼器模擬。本文僅關(guān)注車輛在豎向平面內(nèi)的運動，考慮車體的沉浮zc和點頭βc、轉(zhuǎn)向架的沉浮zt1、zt2和點頭βt1、βt2，以及4個輪對的沉浮zwj1～zwj4，共10個自由度。圖中：mc、mt和mw別為車體、轉(zhuǎn)向架和輪對的質(zhì)量；kpz、cpz和ksz、csz分別為車輛一系和二系懸掛系統(tǒng)的剛度和阻尼系數(shù)；Icy和Ity分別為車體和轉(zhuǎn)向架的點頭轉(zhuǎn)動慣量；lc為車輛定距之半；lt為車輛軸距之半；v為車輛運行速度。車輛主要參數(shù)參見文獻(xiàn)[32]。

車輛子系統(tǒng)的動力方程可表示為

(22)

圖4 車輛模型Fig.4 Vehicle model

2.2 軌道子系統(tǒng)

采用TMM建立軌道子系統(tǒng)模型，無需建立軌道結(jié)構(gòu)整體模型，僅需將周期性軌道結(jié)構(gòu)劃分為若干元胞結(jié)構(gòu)，建立元胞結(jié)構(gòu)的動力方程，并依據(jù)元胞結(jié)構(gòu)間的傳遞關(guān)系，建立軌道結(jié)構(gòu)傳遞模型進(jìn)行求解。對于元胞m，其動力方程可表示為

(23)

(24)

式中:Brvm(1×(Nm,INm,O))為布爾矩陣，該矩陣中與元胞結(jié)構(gòu)鋼軌節(jié)點對應(yīng)自由度的元素為1，其余元素為0;Nm,I和Nm,O分別為元胞m輸入端和輸出端的自由度數(shù)目;Fr為軌道子系統(tǒng)的外荷載向量。

2.3 輪軌間相互作用力

車輛子系統(tǒng)和軌道子系統(tǒng)之間的動力相互作用通過輪軌接觸實現(xiàn)，本文采用基于切線斜率法的線性Hertz接觸模型模擬輪軌接觸關(guān)系。切線斜率法是指過非線性赫茲接觸曲線中靜態(tài)輪軌力P0對應(yīng)的點作切線，切線的斜率即為輪軌接觸彈簧剛度kh，可根據(jù)下式求解

(25)

式中：G為輪軌接觸常數(shù),m/N2/3；P0為靜輪重,N。

在式(22)和(24)中，車輛子系統(tǒng)和軌道子系統(tǒng)通過軌道不平順激勵和輪軌間壓縮引起的輪軌力實現(xiàn)耦合。需要注意的是，假設(shè)車輛在平面內(nèi)運動，軌道為三維空間模型。因此，每個輪對由兩個相同的赫茲彈簧連接在兩條鋼軌上。

式(22)中的車輛外荷載向量Fv可寫為

Fv=Fvg+Fvi+Fvr

(26)

其中，

Fvgj=[mcg0mtg0mtg0mwgmwgmwgmwg]T

(27)

Fvij=[0 0 0 0 0 0fvij1fvij2fvij3fvij4]T

(28)

fviji=2khr(xji)

(29)

式中:Nv為車輛子系統(tǒng)的列車節(jié)數(shù)；Fvgj(10×1)為第j節(jié)列車所受的重力荷載，由各剛體所受重力組成；Fvij(10×1)為車輛所受的軌道不平順激勵；fviji為第j節(jié)車輛第i個輪對受到的軌道不平順激勵；r(x)表示軌道不平順取值，由輪對在軌道上的位置決定。

Fvr(Nv×1)為車輛所受的由輪軌間壓縮引起的輪軌力，需結(jié)合輪下鋼軌進(jìn)行求解，如圖5所示。

圖5 車輛和軌道子系統(tǒng)耦合方式Fig.5 Vehicle and track subsystem coupled mode

對于第j節(jié)車的第i個輪對，fvrji可表示為

fvrji=kh(zrwji-zwji)

(30)

式中:(zrwji-zwji)表示第j節(jié)車的第i個輪對所在位置的輪軌壓縮量;zwji為輪對豎向位移;zrwji為輪下鋼軌的豎向位移。可由輪下元胞中鋼軌單元節(jié)點的位移插值求得

zrwji=XmNT(x)|x=xji

(31)

其中，

Xm=Brvm[Um-1,OUm,O]T

(32)

式中:NT(x)為鋼軌單元的形函數(shù)向量轉(zhuǎn)置；Xm為輪下鋼軌單元的位移；Brvm同上；Um-1,O和Um,O分別為位移矩陣中與軌下元胞的輸入端和輸出端對應(yīng)的元素。

式(24)中的軌道外荷載向量Fr為

Fr=Fri+Frv

(33)

其中，

(34)

(35)

式中:Fri為軌道所受的不平順激勵；Frv為軌道所受的由輪軌間壓縮引起的輪軌力；Njis(1×Nr)為第j節(jié)車第i個輪對所在的第s個鋼軌單元(1個輪對與2個鋼軌單元存在接觸)的形函數(shù)。

2.4 計算流程

圖6給出了基于本文方法的車輛-軌道耦合系統(tǒng)動力響應(yīng)計算分析流程圖，具體計算步驟為：①采用多剛體動力學(xué)方法建立車輛動力方程；②將軌道結(jié)構(gòu)劃分為元胞結(jié)構(gòu)，建立元胞有限元模型，并采用TMM建立軌道結(jié)構(gòu)傳遞矩陣模型；③令系統(tǒng)初始響應(yīng)為零；④根據(jù)列車位置確定軌下元胞，計算輪對和鋼軌位移差，并通過輪軌接觸關(guān)系和軌道不平順求解式(22)和(24)中輪軌間相互作用力；⑤施加輪軌力至對應(yīng)的輪對和鋼軌自由度，通過數(shù)值積分法和傳遞矩陣法求解當(dāng)前時間步車輛及軌道的動力響應(yīng)；⑥判斷收斂性；⑦存儲響應(yīng)，并將其作為下一時間步的初始狀態(tài)，重復(fù)步驟④和⑤直至計算結(jié)束。

圖6 基于TMM的車輛-軌道耦合系統(tǒng)動力響應(yīng)計算流程圖Fig.6 The calculation flow of vehicle-track coupled system dynamic response based on TMM

3 方法驗證

為驗證上述關(guān)于車輛-軌道耦合系統(tǒng)傳遞矩陣算法的正確性，以高速列車作用下的無砟軌道結(jié)構(gòu)為例，對比了TMM和DSM計算所得的輪軌力、車輛位移和軌道動力響應(yīng)結(jié)果。并通過兩種方法計算耗時的對比，進(jìn)一步驗證本文方法的高效性。

CRTSⅡ型板式無砟軌道結(jié)構(gòu)如圖7所示，軌道長度為60 m，扣件間距為0.6 m；鋼軌通過扣件直接與軌道板(內(nèi)置軌枕的預(yù)制板)相連，對稱分布于距軌道板中心0.745 m處；軌道板長6 m，寬2.55 m，高0.2 m，其下鋪設(shè)有CA砂漿層；底座板長6 m，寬3.25 m，高0.3 m，通過板下的滑動層與下部結(jié)構(gòu)相連。假定單節(jié)CRH2動車以300 km/h的速度沿鋼軌從左至右移動，計算時考慮我國新規(guī)范高速鐵路無砟軌道不平順譜樣本，結(jié)構(gòu)阻尼比取為0.02，數(shù)值積分步長取1/10 000 s。

圖7 CRTSⅡ型板式無砟軌道結(jié)構(gòu)示意圖Fig.7 CRTSⅡslab ballastless track

對于空間的軌道結(jié)構(gòu)，其元胞劃分方案與平面軌道結(jié)構(gòu)類似，其中元胞A、B和C長度均為0.1 m。采用有限元法建立軌道結(jié)構(gòu)的元胞模型，并導(dǎo)出其質(zhì)量、剛度和阻尼矩陣進(jìn)行元胞動力方程求解。元胞有限元模型中，鋼軌采用空間梁單元Beam188模擬，軌道板和底座板采用空間板單元Shell181模擬，扣件、砂漿層和滑動層均通過彈簧-阻尼器單元模擬，具體建模參數(shù)參考文獻(xiàn)[34]。

圖8、圖9和圖10給出了TMM和DSM計算所得的第一輪對輪軌力、第一節(jié)車體響應(yīng)及軌道跨中截面處鋼軌響應(yīng)時程曲線。其中，虛線表示采用TMM的計算結(jié)果，實線表示DSM的計算結(jié)果，從圖中可以看出，兩種方法的計算結(jié)果有很好的一致性。以DSM計算結(jié)果為標(biāo)準(zhǔn)，分析兩種方法計算所得各項響應(yīng)的最大值及相對誤差，如表1所示。其中最大值誤差均小于1%，表明基于本文方法建立的軌道結(jié)構(gòu)傳遞矩陣模型，應(yīng)用于車-軌耦合系統(tǒng)動力響應(yīng)分析時具有較高的計算精度。

圖8 第一輪對輪軌力時程曲線Fig.8 Time history curve of vertical wheel-rail contact force of the first wheel set

圖9 第1節(jié)車體豎向加速度時程曲線Fig.9 The dynamic response time history curve of the vehicle

(a) 豎向位移時程

(b) 豎向加速度時程圖10 跨中截面處鋼軌動力響應(yīng)時程曲線Fig.10 The dynamic response time history curve of midpoint of rail

表1 兩種方法計算精度對比Tab.1 Comparison of calculation accuracy

為進(jìn)一步說明本文方法具有更高的計算效率，對比了TMM和DSM兩種方法的計算耗時，如表2所示。其中，采用DSM求解軌道結(jié)構(gòu)動力響應(yīng)時，需要建立無砟軌道結(jié)構(gòu)的整體模型，動力方程最大矩陣階數(shù)為50 484；采用TMM時，僅需建立元胞結(jié)構(gòu)模型，極大地減少了軌道自由度數(shù)目，動力方程最大矩陣階數(shù)僅為84，為整體模型的1/601。對比TMM和DSM的計算耗時，當(dāng)無砟軌道長度為60 m時，TMM計算耗時為DSM的17.1%。表明基于本文方法求解車輛-軌道耦合系統(tǒng)模型動力響應(yīng)時，僅需對各個軌道元胞結(jié)構(gòu)的動力方程進(jìn)行求解，無需建立和求解整體軌道結(jié)構(gòu)的動力方程，有效地降低了計算過程中動力方程系數(shù)矩陣的規(guī)模，提升了計算效率。

表2 兩種方法計算耗時對比

4 結(jié) 論

本文針對車輛-軌道耦合系統(tǒng)動力分析問題，利用軌道結(jié)構(gòu)的周期性特征，基于傳遞矩陣法提出了一種便捷的建模及求解方案。針對有砟軌道和CRTSⅡ型無砟軌道系統(tǒng)結(jié)構(gòu)，分別提出了對應(yīng)的元胞劃分方案和傳遞矩陣模型，并應(yīng)用于車輛-軌道耦合系統(tǒng)動力響應(yīng)分析。以高速動車通過無砟軌道結(jié)構(gòu)為例，對其正確性進(jìn)行了驗證。得到了以下結(jié)論：

(1) 本文方法計算結(jié)果與直接剛度法計算結(jié)果吻合度極高，表明基于本文方法應(yīng)用于車-軌耦合系統(tǒng)動力響應(yīng)分析時具有較高的計算精度。

(2) 本文方法在建立車輛-軌道耦合系統(tǒng)模型時應(yīng)用了軌道結(jié)構(gòu)的周期特性，只需建立其單跨軌道劃分而成的元胞模型，無需建立其整體模型，使得建模更為便捷。

(3) 基于傳遞矩陣法求解車輛-軌道耦合系統(tǒng)模型動力響應(yīng)時，僅需對各個軌道元胞結(jié)構(gòu)的動力方程進(jìn)行求解，無需組裝整體軌道結(jié)構(gòu)的動力方程，有效地降低了計算過程中涉及矩陣的規(guī)模，提升了計算效率。