基于動柔度法的車?線?橋垂向耦合振動分析

石廣田，楊建近，楊新文，張小安

?

基于動柔度法的車?線?橋垂向耦合振動分析

石廣田1，楊建近1，楊新文2，張小安1

(1. 蘭州交通大學機電工程學院，甘肅蘭州，730070；2. 同濟大學交通運輸工程學院，上海，201804)

利用動柔度的思想建立高速鐵路車輛?軌道?橋梁垂向耦合動力學頻域模型。將車輛看作具有10自由度的多剛體系統，推導出車輛系統的動柔度；將鋼軌看作無限長Timoshenko梁、軌道板簡化為自由?自由的Euler梁、橋梁簡化為簡支Euler梁，扣件系統和CA砂漿層用線性彈性阻尼單元模擬，從而建立線?橋垂向動力相互作用系統，利用動柔度的思想求得線?橋垂向動力相互作用系統的動柔度；利用線性化Hertz接觸理論，將車輛系統與 線?橋垂向動力相互作用系統耦合成車?線?橋垂向耦合振動模型。分別求解單位簡諧激勵、輪軌幾何粗糙度激勵下的軌道系統和橋梁結構的動力學響應，最后結合虛擬激勵法求解軌道譜激勵下的軌道系統和橋梁結構的隨機振動響應。研究結果表明：利用動柔度思想建立的頻域車輛?軌道?橋梁垂向耦合動力學模型，具有邏輯清晰、求解快速的特點，能夠直接求解系統在頻域的動力學響應。

耦合振動；隨機振動；動柔度；車?線?橋耦合動力學；高速鐵路

由于列車運行速度高和無砟板式軌道的線路結構會使得輪軌相互作用急劇增加，導致高速鐵路車輛?軌道?橋梁耦合振動問題更突出。分析車輛?軌道?橋梁耦合系統的振動，對于確保高速鐵路運營安全、提高乘坐舒適性、降低高速鐵路運營引起的振動對沿線居民和結構建筑的影響有重要的理論價值和工程應用價值。松浦章夫[1]先后建立了半車模型和10自由度的整車模型分析車橋動力學問題。TANABE 等[2]開發了列車?軌道相互作用計算機程序，其中車輛模型具有31自由度并考慮了車輛懸掛系統的非線性，將軌道和橋梁作為一體用有限元法得出軌道?橋梁運動方程，用模態分解技術減小軌道橋梁結構自由度，以輪軌間不平衡力收斂為迭代終止條件，采用Newmark法求解列車、軌道與橋梁的動力響應結果。BHATTI等[3]和WANG等[4]在車橋空間動力分析中分別建立了21個自由度和23個自由度的車輛空間振動模型，并考慮了車輛懸掛的非線性。FRYBA[5]系統地研究了車輛參數、行車速度、軌道參數、橋梁參數等對橋梁振動的影響，并對橋梁結構的疲勞問題及隨機振動問題進行了較深入的研究。曾慶元等[6?7]利用勢能駐值原理和 “對號入座”法則，導出了車橋時變系統的動力學方程，其中機車車輛模型具有21自由度，橋梁利用有限元法進行離散。鄧子銘等[8]利用“對號入座”法則建立了列車與鋼桁梁橋耦合動力學模型，將軌道不平順作為系統的自激激勵源，地震作為外部激勵，分析地震對鋼桁梁橋車橋系統耦合振動的影響。夏禾等[9]將車輛視為多剛體系統，分析了車?橋?墩體系的耦合動力學問題，也研究了脈動風、地震荷載對列車過橋時的動力學響應和行車安全的影響。李小珍等[10?11]采用23個自由度的車輛模型，橋梁采用有限元模型，研究了車橋系統的空間動力響應。翟婉明等[12?13]建立了較為完備的車輛?軌道?橋梁動力相互作用理論研究，并開發了相應的計算機程序，為鐵路新橋設計、舊橋加固提供了一套理論方法和安全評估。以上建立的車?線?橋耦合動力學模型多為時域模型，一方面大多數學者考慮了系統的非線性因素，但同時會使得求解難度大大增加；另一方面，若要分析系統在頻域的動力學響應，還要進行時頻轉換分析，較為復雜。本文作者利用動柔度的思想建立高速鐵路車輛?軌道?橋梁垂向耦合動力學頻域模型[14]，能夠快速、準確地求解系統在頻域的振動響應。

1 車輛?軌道?橋梁耦合動力學模型

本文建立的車輛?軌道?橋梁耦合動力學模型如圖1所示。車輛為10自由度的多剛體系統，鋼軌、軌道板、橋梁分別用無限長Timoshenko梁、自由?自由Euler梁、簡支Euler梁模擬，扣件系統和CA砂漿層用線性彈性阻尼單元模擬，輪軌接觸關系采用線性化Hertz彈性接觸理論。

圖1 車輛?軌道?橋梁動力相互作用示意圖

單節高速列車車輛的振動幅值為

(2)

式中：[]=[][][]T為單節車輛在輪對處的柔度矩陣；為轉換矩陣；{()}為單節車輛受到的垂向輪軌作用力幅值。

鋼軌被視為無限長Timoshenko梁，其動柔度函數為[15]

β(12)為在位置2處施加單位力在1處引起的位移；1，2，1和2與振動波沿鋼軌的傳播有關。鋼軌的振動幅值為[15?16]

(4)

式中：x和x分別為第個激勵的坐標、第個扣件的坐標；Z為軌道板在坐標x處的垂向位移；K(Z?Z)為第個扣件的彈性恢復力。

多個軌道板的動柔度可表示為β(12)= diag[ [1]1，…，[β1]，…，[β1]]。其中，為軌道板數量，β1為一塊軌道板的兩端自由Euler梁模型的動柔度[16]，

W()為兩端自由的Euler梁的第階振型的陣型函數[17]。則整個軌道板的運動微分方程為

(6)

式中：N為軌道板下離散分布彈性支承的數量；x為軌道板下第個離散分布彈性支承點的坐標；K(Z?Z)為軌道板下第個離散支撐的彈性恢復力。

高架橋梁的運動微分方程為

式中：β(1，2)為高架橋梁的簡支Euler梁模型的動柔度，其表達式可通過將式(5)中的陣型函數換成橋梁簡支Euler梁模型的陣型函數得到；N為高架橋梁支撐的數量；x為高架橋梁第個支撐的坐標；Z為橋梁在x處的位移；KZ為高架橋梁第個支撐的彈性恢復力。

式(4)，(6)和(7)中均未明示阻尼，實際上梁的彎曲剛度以及彈簧的剛度都是包含損耗因子的復剛度。將式(4)，(6)和(7)合并，可得到矩陣形式的表達式：

[]主要由鋼軌、軌道板和橋梁結構的動柔度乘以剛度形成；{}由待求解的鋼軌、軌道板和橋梁結構的位移組成；{}為荷載矩陣。求解式(8)即可以得出激勵作用下軌道系統和橋梁結構的頻域位移響應。

車輛與軌道耦合關系通過輪軌相互作用力來實現。從下式可求得垂向輪軌力：

式中：[T]為軌道系統在輪軌接觸點處的原點垂向位移導納，由式(8)可求得；[]為單節車輛在輪對處的柔度矩陣，由式(2)可求得；[]為輪軌接觸彈簧的動柔度，β=1/k，k為線性化輪軌接觸剛度，以CRH2型列車為例，軸質量為14 t，LMA型磨耗踏面，滾動圓半徑為430 mm，根據Hertz彈性接觸理論，可求得k為1.1718×106kN/m；Δ()為輪軌不平順輸入。在得到動態輪軌力幅值后，按照式(1)和式(8)即可以計算車輛?軌道?橋梁耦合系統的動態響應。

2 系統的動力學響應

計算采用的車輛為CRH2型高速列車車輛。軌道結構為CRTS?I型板式無砟軌道，高架橋梁為我國高速鐵路單線32 m簡支箱梁橋，動力學參數如表1所示。

表1 模型的計算參數

2.1 單位簡諧荷載激勵

在鋼軌上施加單位簡諧荷載時，得到的軌道?橋梁相互作用系統的動力學響應就是軌道?橋梁相互作用系統的動柔度。輪軌力可以看作是一系列一定頻率的簡諧力的疊加，因此，分析軌道?橋梁相互作用系統在單位簡諧荷載作用下的響應，能夠在一定程度上反映軌道?橋梁相互作用系統的動力學特性。

在橋梁跨中截面的鋼軌上(扣件跨中)施加單位簡諧荷載，以1 Hz為步長計算系統的導納。由于橋梁跨中截面位于2個軌道板間的間隙中，故取靠近橋梁跨中截面扣件下的軌道板的垂向位移響應及橋梁跨中截面上鋼軌和橋梁的垂向位移響應為研究對象。軌道?橋梁相互作用系統的動柔度幅值和其相位分別如圖2和圖3所示。

1—鋼軌；2—軌道板；3—橋梁。

從圖2可知：在200 Hz以下頻率范圍內，軌道板的垂向位移響應和橋梁的變化趨勢幾乎是一致的，軌道板在此頻率范圍的位移響應幅值出現了3個峰值，也是伴隨著橋梁結構的位移響應幅值的峰值出現的，因此，可以推斷在垂向簡諧荷載作用下，200 Hz以下頻率范圍內軌道板是隨著橋梁在垂向產生振動的，在特別是在20 Hz以下，軌道板的位移幅值和橋梁的完全一致。在9 Hz左右鋼軌、軌道板和橋梁的垂向位移幅值同時達到最大值，而此頻率正是橋梁結構的一階豎彎頻率。

而頻率在25 Hz以上時，鋼軌的垂向振動位移幅值比軌道板的大；僅頻率在4~7 Hz范圍時，橋梁的垂向振動位移幅值稍比軌道板的大，而在其他頻率范圍時都比軌道板的小；隨著頻率的增加，鋼軌的垂向振動位移幅值大于軌道板的垂向振動位移幅值，軌道板的垂向振動位移幅值大于橋梁的垂向振動位移幅值，并且它們之間的差值越來越明顯。鋼軌、軌道板和橋梁都在軌道?橋梁相互作用系統的自振頻率處出現峰值，因此，可推斷這些系統自振頻率對鋼軌、軌道和橋梁振動都有較大的影響。

1—鋼軌；2—軌道板；3—橋梁。

從圖3可知：當頻率在6 Hz以下時，鋼軌的動柔度相位角比軌道板的大，而軌道板的動柔度相位角比橋梁的大，而在6 Hz左右時三者都急劇減小；鋼軌的動柔度相位在6~165 Hz范圍內幾乎保持不變，在一階pinned-pinned頻率(970 Hz左右)鋼軌的動柔度相位發生了突變；在300 Hz左右軌道板的動柔度相位幅值突變為最大；在100~450 Hz橋梁的動柔度相位出現4個大于零的峰值，并依此減小；整個頻率范圍內，鋼軌的動柔度相位角始終處于負值，而軌道板與橋梁的比較復雜，特別是在1 000 Hz以上，隨著頻率的增加，軌道板和橋梁的動柔度相位呈現較明顯的震蕩趨勢。

結構振動衰減率用來表示結構上隨著距離激勵點距離改變振動衰減的特性[16]，定義為

式中：0為激勵點處結構振動位移；z為距離激勵點處結構振動位移。

在橋梁跨中截面的鋼軌上(扣件跨中)施加單位簡諧荷載，分別在鋼軌、軌道板和橋梁上沿鋼軌方向取距離點=6.25 m處結構的振動，分析振動在高架軌道系統中的縱向衰減。圖4所示為振動在鋼軌、軌道板、橋梁沿縱向上衰減的情況。

從圖4可知：在77 Hz左右鋼軌的振動沿縱向的衰減達到最大值，而在此頻率軌道板和橋梁的振動沿縱向的衰減卻達到1個谷值。鋼軌的振動沿縱向的衰減在鋼軌?扣件系統的自振頻率(230 Hz左右)和第一個pinned-pinned頻率達到峰值，軌道板和橋梁的振動衰減率在其自身固有頻率處都會產生峰值或谷值。總體來看，鋼軌的振動沿縱向的衰減在1 000 Hz以下衰減較明顯，而在1 000 Hz以上高頻段鋼軌沿縱向的衰減較小。鋼軌振動沿縱向的衰減率變化較簡單，軌道板和橋梁的振動沿縱向的衰減率變化較復雜。

(a) 鋼軌；(b) 軌道板；(c) 橋梁

2.2 輪軌幾何粗糙度激勵

由于我國對高鐵的研究是近十幾年才開始的，尚沒有能反映我國高鐵輪軌組合粗糙度譜，計算中采用文獻[15]中給出的一組考慮輪軌接觸濾波效應的輪軌組合粗糙度譜，如圖5所示。利用建立的高速車輛?軌道?橋梁耦合動力學模型，求解車速為250 km/h時系統的垂向振動響應。圖6所示為車輪速度為250 km/h時輪軌組合粗糙度激勵下的4組垂向輪軌動態相互作用力幅值，圖6中的1，2，3和4分別對應圖1中垂向輪軌動態相互作用力的下標。取第三位輪對所在橫截面內的鋼軌、軌道板和橋梁為分析對象，車速為250 km/h時它們在輪軌組合粗糙度激勵下的垂向振動加速度幅值見圖7。

圖5 輪軌組合粗糙度譜(考慮接觸濾波的影響)[14]

圖6 輪軌組合粗糙度激勵下的垂向輪軌力幅值

(a) 鋼軌；(b) 軌道板；(c) 橋梁

由圖6可以看出：4個輪對上的垂向輪軌動態相互作用力幅值變化趨勢是一致的，總體上看，當頻率低于50 Hz時，隨著頻率的增加而增加；在50~210 Hz頻段內，隨著頻率的增加而減小；在210~970 Hz，隨著頻率的增加總體上增加；當頻率大于970 Hz時，隨著頻率的增加，4個輪對上的垂向輪軌動態相互作用力幅值總體上減小。4個輪對上的垂向輪軌動態相互作用力幅值的最大值為47 142 N。

由圖7可知：鋼軌垂向振動加速度幅值在低頻段很小，在500~3 500 Hz頻段特別是1 000 Hz左右，鋼軌垂向振動加速度幅值比較大，在一階pinned- pinned頻率(970 Hz左右)達到最大值963.6 m/s2。軌道板的垂向振動加速度幅值主要集中300~1 500 Hz的中高頻段，在一階pinned-pinned頻率附近達到最大值16 m/s2，主要是一階pinned-pinned頻率施加到軌道板上的扣件彈性力而導致。橋梁垂向振動加速度幅值主要集中在300 Hz以下的低頻段。

總體上看，低頻區域鋼軌的垂向振動加速度幅值較小，高頻區域鋼軌的垂向振動加速度幅值較大；橋梁的垂向振動加速度幅值主要集中在低頻區域，而在高頻區域較小。

2.3 軌道譜激勵

國內關于高速鐵路運用條件下的軌道譜尚沒有全國代表性的軌道譜分析式和標準，國外典型的高速鐵路軌道譜標準為德國高速軌道譜。德國高速軌道譜的波長范圍為1~100 m，當高速列車通過時激發的系統振動為低頻振動，而短波不平順功率譜能夠激發系統的高頻振動，常用的短波不平順功率譜為Sato譜。本節以德國低干擾譜高低不平順(波長范圍為1~100 m)和Sato譜(波長1 m以下)為激勵，利用建立的高速車輛?軌道?橋梁耦合動力學模型及本文作者在文獻[18]中采用的虛擬激勵法，求解車速為250 km/h時系統的隨機振動響應。Sato譜表達式為

式中：和是由軌道狀態決定的參數，本文選擇軌道狀態為好時的值，=0.065，=3.06。

圖8所示為車速250 km/h時垂向輪軌動態相互作用力幅值譜密度。取圖1所示車輛3號輪對所在平面內的鋼軌(輪軌接觸點處)、軌道板及橋梁的垂向隨機振動位移響應為研究對象，當車速為250 km/h時，德國低干擾譜激勵下車輛?軌道?橋梁相互作用系統的隨機振動垂向加速度功率譜見圖9。

從圖9可以看出：在車速為250 km/h時，波長為1~100 m的德國低干擾譜激發的車輛?軌道?橋梁相互作用系統的隨機振動為70 Hz以下的低頻振動，系統70 Hz以上的振動是由Sato譜激發的；隨著頻率增加，垂向輪軌力幅值譜在70 Hz左右和650 Hz左右有2個突出的峰值，因此，垂向輪軌力幅值譜整體上呈現先增加后減小再增加再減小的趨勢。鋼軌在650 Hz和970 Hz左右的垂向振動加速度功率譜達到最大值，后者是由于鋼軌振動的pinned-pinned頻率引起的。軌道板的垂向振動加速度功率譜在650 Hz左右達到最大值，在1 200 Hz以上軌道板的垂向振動加速度功率譜很小。橋梁的垂向振動加速度主要集中在300 Hz以下的低頻段，當頻率高于300 Hz時，橋梁的垂向振動加速度功率譜幾乎為零。

圖8 垂向輪軌力幅值譜

(a) 鋼軌；(b) 軌道板；(c) 橋梁

3 結論

1) 隨著頻率的增加，鋼軌的垂向動柔度幅值大于軌道板的動柔度幅值，軌道板的動柔度幅值大于橋梁的動柔度幅值。在整個頻率范圍內，鋼軌的動柔度相位角始終處于負值，而軌道板與橋梁的動柔度相位變化比較復雜，特別是在1 000 Hz以上隨著頻率的增加呈現較為明顯的震蕩趨勢。鋼軌振動沿縱向的衰減率變化較簡單，在1 000 Hz以下衰減較明顯而在1 000 Hz以上衰減較小，軌道板和橋梁的振動沿縱向的衰減率變化較復雜。

2) 在車速為250 km/h時，在輪軌幾何不平順激勵下或在德國低干擾譜及Sato譜激勵下，鋼軌垂向振動加速度主要集中在500~3 500 Hz頻段特別是1 000 Hz左右，鋼軌垂向振動加速度較大，而在低頻段很小；軌道板的垂向振動加速度主要集中300~1 500 Hz的中高頻段；橋梁垂向振動加速度主要集中在300 Hz以下的低頻段，在300 Hz以上幾乎為零。

3) 利用動柔度思想建立的頻域車輛?軌道?橋梁垂向耦合動力學模型，具有邏輯清晰、求解快速的特點，能夠直接求解系統在頻域的動力學響應。

[1] 松浦章夫. 高速鐵路における車輛と橋げたの動的舉動に關する研究[J]. 土木學會論文集(日), 1976, 258(12): 35?47. MATSUURA A. Study of dynamic behaviors of bridge girders for high-speed railway[J]. Journal of JSCE, 1976, 258(12): 35?47.

[2] TANABE M, WAKUI H, MATSUMOTO N, et al. Computational model of a Shinkansen train running on the railway structure and the industrial applications[J]. Journal of Materials Processing Technology, 2003, 140(1/2/3): 705?710.

[3] BHATTI M H, GARGV K, CHU K H. Dynamic interaction between freight train and steel bridge[J]. Journal of Dynamic Systems Measurement and Control, 1985, 107(1): 60?66.

[4] WANG T L, CHU K H. Railway bridge/vehicle interaction studies with new vehicle model[J]. Structural Engineering, 1991, 117(7): 2099?2116.

[5] FRYBA L. Vibration of solids and structures under moving loads[M]. Netherlands: Noordhoff International Publishing, 1972: 157?172.

[6] 曾慶元, 郭向榮. 列車橋梁時變系統振動分析理論及應用[M]. 北京: 中國鐵道出版社, 1999: 49?93.ZENG Qingyuan, GUO Xiangrong. Theory and application of train-bridge time variant system vibration analysis[M]. Beijing: China Railway Press, 1999: 49?93.

[7] 曾慶元, 周智輝, 赫丹, 等. 列車?軌道(橋梁)系統橫向振動穩定性分析[J]. 鐵道學報, 2012, 34(5): 86?90. ZENG Qingyuan, ZHOU Zhihui, HE Dan, et al. Study on stability of lateral vibration of train-track(bridge) system[J]. Journal of the China Railway Society, 2012, 34(5): 86?90.

[8] 鄧子銘, 郭向榮, 張志勇. 地震作用對鋼桁梁橋車橋系統耦合振動的影響分析[J]. 中南大學學報(自然科學版), 2011, 42(1): 184?191. DENG Ziming, GUO Xiangrong, ZHANG Zhiyong. Coupled vibration of train-bridge system of steel truss bridge with seismic effect[J]. Journal of Central South University (Science and Technology), 2011, 42(1): 184?191.

[9] 夏禾, 張楠. 車輛與結構動力相互作用[M]. 2版. 北京: 科學出版社, 2005: 37?198.XIA He, ZHANG Nan. Dynamic interaction of vehicle and structure[M]. 2nd ed. Beijing: Science Press, 2005: 37?198.

[10] 李小珍, 朱艷, 強士中. 高速列車作用下簡支梁車橋耦合振動隨機響應分析[J]. 振動與沖擊, 2012, 31(4): 168?172.LI Xiaozhen, ZHU Yan, QIANG Shizhong. Stochastic response analysis of train-bridge coupling system under high speed train loads excitation[J]. Journal of Vibration and Shock, 2012, 31(4): 168?172.

[11] 李小珍, 馬文彬, 強士中. 車橋系統耦合振動分析的數值解法[J]. 振動與沖擊, 2002, 21(3): 21?25.LI Xiaozhen, MA Wenbin, QIANG Shizhong. Coupling vibration analysis of vehicle-bridge system by iterative solution method[J]. Journal of Vibration and Shock, 2002, 21(3): 21?25.

[12] 翟婉明, 夏禾. 列車?軌道?橋梁動力相互作用理論與工程應用[M]. 北京: 科學出版社, 2011: 29?157.ZHAI Wanming, XIA He. Train?track?bridge dynamic interaction: thory and engineering application[M]. Beijing: Science Press, 2011: 29?157.

[13] 蔡成標. 高速鐵路列車?線路?橋梁耦合振動理論及應用研究[D]. 成都: 西南交通大學牽引動力國家重點實驗室, 2004: 21?122.CAI Chengbiao. Theory and application of train?track?bridge coupling vibration in high-speed railway[D]. Chengdu: Southwest Jiaotong University. State Key Laboratory of Traction Power, 2004: 21?122.

[14] 楊建近. 高速車輛-橋梁耦合振動與高架橋梁結構噪聲預測研究[D]. 蘭州: 蘭州交通大學機電工程學院, 2015: 8?23. YANG Jianjin. Study on prediction of high speed vehicle-bridge coupling vibration and viaduct bridge structure noise radiation[D]. Lanzhou: Lanzhou Jiaotong University. School of Mechatronic Engineering, 2015: 8?23.

[15] 劉海平. 高速鐵路輪軌滾動噪聲建模、預測與控制研究[D]. 上海: 上海交通大學機械系統與振動國家重點實驗室, 2011: 30?39.LIU Haiping. A study on modeling, prediction and its control of wheel/rail rolling noises in high speed railway[D]. Shanghai: Shanghai Jiaotong University. State Key Laboratory of Mechanical System and Vibration, 2011: 30?39.

[16] 楊建近, 楊新文. 鋼軌吸振器對高架結構垂向振動的影響[J]. 噪聲與振動控制, 2014, 34(4): 214?218.YANG Jianjin, YANG Xinwen. Effect of rail vibration absorbers on vertical vibration of elevated structures[J]. Noise and Vibration Control, 2014, 34(4): 214?218.

[17] 翟婉明. 車輛?軌道耦合動力學[M]. 4版. 北京: 科學出版社, 2015: 39?41.ZHAI Wanming. Vehicle?tack coupling dynamics[M]. 4th ed. Beijing: Science Press, 2015: 39?41.

[18] 楊新文. 高速鐵路輪軌噪聲理論計算與控制研究[D]. 成都: 西南交通大學牽引動力國家重點實驗室, 2006: 71?73.YANG Xinwen. Theorical analysis and control studies in wheel/rail noises of high speed railway[D]. Chengdu: Southwest Jiaotong University. State Key Laboratory of Traction Power, 2006: 71?73.

(編輯 楊幼平)

Vertical vehicle?track?bridge coupling vibration based on dynamic flexibility method

SHI Guangtian1, YANG Jianjin1, YANG Xinwen2, ZHANG Xiaoan1

(1. School of Mechat ronic Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China;2. School of Transportation Engineering, Tongji University, Shanghai 201804, China)

A frequency domain dynamic model of vehicle?track?bridge vertical coupling system of high-speed railway was established by using the dynamic flexibility method. The vehicle subsystem was treated as the multi-rigid-body system with 10 degrees of freedom, and its dynamic flexibility was deduced by using the dynamic flexibility method. The rail, slab and bridge were respectively considered as the infinite Timoshenko beam, the freedom-freedom Euler beam and the simply supported Euler beam, and the fastener system and CA mortar layer were modeled as linear elastic and damping elements. Therefore, a vertical track-bridge dynamic interaction subsystem was developed, and its dynamic flexibility was also deduced by using the dynamic flexibility method. Using linear Hertz contact theory, the vehicle system and the vertical track-bridge dynamic interaction subsystem were coupled into a vertical vehicle-track-bridge coupling dynamic system. The dynamic responses of the track and the bridge, respectively excited by the unit harmonic force and the wheel/rail geometric irregularity, were calculated. Finally, combined with the pseudo-excitation method, the random vibration responses of the track and the bridge excited by the track spectrum were also calculated. The results show that the vehicle-track-bridge coupling dynamic model proposed in this paper has the characteristics of clear logicality and fast computation, and it can directly solve the dynamic responses of system in the frequency domain.

coupling vibration; random vibration; dynamic flexibility; vehicle?track?bridge coupled dynamics; high speed railway

U211.5 U213

A

1672?7207(2017)04?1119?08

10.11817/j.issn.1672?7207.2017.04.036

2016?04?08；

2016?06?23

國家自然科學基金資助項目(51165017)；蘭州交通大學校青年基金資助項目(2015025)(Project(51165017) supported by the National Natural Science Foundation of China? Project (2015025) supported by the Youth Fund of Lanzhou Jiaotong University)

石廣田，博士，教授，從事軌道交通系統動力學研究；E-mail：shigt@mail.lzjtu.cn