車載作用下公路橋梁耦合振動精細化建模及驗證分析

殷新鋒 晏萬里 任厚乾 劉揚

摘? ?要：現(xiàn)有車-橋耦合振動分析中車輛模型不能精確考慮車輛動力特性和柔性輪胎對車橋耦合振動響應的影響.為了進一步研究充氣輪胎胎壓對車-橋耦合振動的影響，基于LS-DYNA程序，采用線彈性橡膠材料模擬輪胎并定義輪胎內(nèi)氣壓，結(jié)合常用重載三軸汽車的結(jié)構(gòu)參數(shù)，運用彈簧阻尼單元及梁、殼單元模擬車輛懸架系統(tǒng)的動力特性，建立可分析車輪氣壓的三維車輛模型;并基于實橋試驗結(jié)果及響應面法得到高精度有限元橋梁模型;通過顯式求解程序LS-DYNA內(nèi)置的接觸算法，將車輛子系統(tǒng)和橋梁子系統(tǒng)聯(lián)立耦合起來，形成顯式的車-橋耦合振動分析模型.計算結(jié)果與實測結(jié)果對比分析驗證了該方法的正確性，并分析了輪胎胎壓對橋梁振動的影響.

關(guān)鍵詞：響應面法;有限元模型修正;連續(xù)剛構(gòu)橋;LS-DYNA;動力響應

中圖分類號：U441.2? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼：A

Fine Modeling of Coupled Vibration of Highway Bridge under

Vehicle Loading and Verification Analysis

YIN Xinfeng1，YAN Wanli1，REN Houqian1，LIU Yang1，2

（1. College of Civil Engineering，Changsha University of Science & Technology，Changsha 410114，China;

2. College of Civil Engineering，Hunan University of Technology，Zhuzhou 412007，China）

Abstract：In the existing coupling vibration analysis of vehicle bridge，the influence of vehicle dynamic characteristics and flexible tire on the vehicle bridge coupling vibration response can not be accurately considered. In order to further study the influence of pneumatic tire pressure on vehicle bridge coupling vibration，the linear elastic rubber material is used to simulate the tire and define the air pressure in the tire based on LS-DYNA program. Combined with the structural parameters of the commonly used heavy-duty three-axle vehicle，the dynamic characteristics of the vehicle suspension system are simulated by using the spring damping element，beam and shell element，and a three-dimensional vehicle model is established. Based on the test results of the real bridge and the response surface method，a high-precision finite element bridge model is obtained. The vehicle subsystem and the bridge subsystem are coupled together to form the explicit vehicle bridge coupling vibration analysis model through the contact algorithm built in the explicit solution program LS-DYNA. The correctness of the method is verified by comparing the calculation results with the measured results，and the influence of tire pressure on bridge vibration is analyzed.

Key words：response surface method;finite element model updating;continuous rigid-frame bridge;LS-DYNA;dynamic response

車輛在橋梁上行駛時產(chǎn)生的激勵導致車輛和橋梁產(chǎn)生相互振動，該振動稱為車橋耦合振動[1]. 國內(nèi)外學者對車橋耦合振動問題進行了大量研究，并取得非凡成果.夏禾[2]將車輛簡化為懸掛振動模型，研究了車-橋-墩相互作用的動力響應. Cai等[3-4]采用兩軸車輛模型，基于功率譜密度函數(shù)生成隨機橋面不平整度，分析了車輛荷載作用下橋梁振動響應. Huang等[5]采用三軸車輛有限元模型研究了簡支梁的車橋隨機振動. 韓萬水等[6]結(jié)合實測數(shù)據(jù)對橋梁模型進行修正以得到橋梁基準模型，并采用梁格法對車橋振動進行分析. 李奇等[7]考慮車體柔性的影響，分析了高速列車和簡支梁橋相互作用的影響. 鄧露等[8]運用動態(tài)稱重技術(shù)（BWIM）實時監(jiān)測車輛總重和軸重，進行中小型跨徑混凝土梁橋的振動研究.趙越等[9]基于等參映射與改進折半法在傳統(tǒng)車橋耦合分析的基礎(chǔ)上進一步提升其分析精度及計算效率，從而進行公路車橋耦合分析.但關(guān)于充氣輪胎胎壓對車-橋耦合振動影響的研究較少.主要原因為現(xiàn)有車-橋耦合振動分析中車輛模型多為簡化的質(zhì)量-彈簧-阻尼多自由度振動模型[10-13]，該模型常將車輪簡化為點或者等效線面接觸，故不能精確考慮車輛動力特性和柔性輪胎對車橋耦合振動響應的影響.因此，建立精確車輛模型和橋梁模型是至關(guān)重要的，這樣才能提高數(shù)值模擬精度，為橋梁結(jié)構(gòu)的健康運營提供有利建議.

本文以一座連續(xù)剛構(gòu)箱梁橋為工程背景.首先，基于LS-DYNA程序以車輛實際構(gòu)造及動力特性為基準，建立車輛精細三維有限元模型;然后，結(jié)合響應面法并依據(jù)實橋試驗結(jié)果對橋梁模型進行修正，以得高精度的橋梁有限元模型;最后，通過LS-DYNA程序?qū)蛄耗Ｐ秃蛙囕v模型進行耦合，求解車橋的振動響應，并通過改變車輪氣壓，分析該參數(shù)對橋梁振動響應的影響.

1? ?車輛模型的建立及校驗

1.1? ?車輛模型的建立

在車橋耦合振動模型中，影響橋梁振動響應的主要因素為車輛模型的動力特性和荷載分布. 車輛模型需著重模擬懸架系統(tǒng)、車輪和軸重. 本文參照東風牌三軸載重自卸貨車，建立車輛有限元模型，其相關(guān)參數(shù)如表1所示.前中軸距為3.5 m，中后軸距為1.4 m，后輪距為1.8 m. 基于CAD軟件SolidWorks分別建立各部位的幾何模型并進行網(wǎng)格劃分，最后使用梁、殼、實體單元及離散單元等賦予不同部位網(wǎng)格不同的屬性，從而構(gòu)建車輛模型，如圖1所示.

車輪由輪胎、輪盤和輪轂組成.輪胎采用線彈性橡膠材料，輪盤和輪轂采用線彈性鋼材材料，車輪模型中各部位連接均為剛性連接，邊界為剛性固態(tài)約束，有限元模型如圖2所示. 根據(jù)輪胎內(nèi)實際氣壓，使用LS-DYNA程序關(guān)鍵字*AIRBAG_SIMPLE_AIRBAG_MODEL定義由輪胎、輪轂組成的封閉體內(nèi)的氣壓，其氣體壓力值為0.6 MPa.

由于車輛懸架系統(tǒng)結(jié)構(gòu)復雜，本文采用殼單元、剛體、彈簧阻尼單元和多點約束來模擬懸架系統(tǒng).前后懸架的彈簧和阻尼器的參數(shù)參考文獻[14]，前懸架的彈簧剛度取800 N/mm，阻尼系數(shù)取20 Ns/mm;后懸架的彈簧剛度取1 200 N/mm，阻尼系數(shù)取25 Ns/mm. 前后懸架實體構(gòu)造及有限元模型如圖3所示.

為了使車輪正常轉(zhuǎn)動，在輪盤與車軸之間設置旋轉(zhuǎn)鉸，再通過定義*INITIAL_VELOCTIY_GENERATION關(guān)鍵字設置車輪的轉(zhuǎn)動和平動速度，實現(xiàn)車輪滾動向前的狀態(tài)，如圖4所示.

1.2? ?車輛模型的校驗

在車輛模型的3個軸上選擇6個節(jié)點，約束豎向位移. 然后對車輛施加重力荷載使車輛產(chǎn)生瞬時振動，再通過定義關(guān)鍵字*DAMPING_GLOBAL對車輛模型施加全局阻尼，使其快速達到穩(wěn)定狀態(tài)，計算出約束反力，并與實測車輛軸重進行比較，最終得到車輛模型的軸載曲線如圖5所示.

由表2知，車體總重偏差為0.87%，說明實測車輛軸重分布特性能體現(xiàn)于建立的車輛模型.

通過對比車輛模型和實測車輛的自振頻率，驗證車輛模型動力特性的有效性. 如圖6所示，該車實測自振頻率為1.635 Hz，車輛有限元模型的自振頻率為1.647 Hz，顯然，兩者基頻相差很小.

2? ?橋梁模型

2.1? ?橋梁概況

該橋主橋為三跨變截面預應力混凝土連續(xù)剛構(gòu)結(jié)構(gòu)，跨徑布置為（65+120+65） m，其立面及斷面示意圖分別如圖7和圖8所示.

2.2? ?實橋試驗

利用有限元軟件ANSYS建立初始有限元模型如圖9所示. 主梁和墩身都用實體單元模擬，墩底采用固結(jié)連接. 混凝土密度取為2 500 kg/m3，主梁和墩的混凝土彈性模量分別為34.5 GPa和32.5 GPa.

現(xiàn)場對橋梁進行靜力試驗，測試中加載車輛總數(shù)為6輛，每輛總重為35 t，車輛前軸重7 t，中軸和后軸各重14 t，車輛照片如圖10所示.

在正式試驗前先進行預加載，以消除非彈性變形，確保試驗及設備處于良好工作狀態(tài). 正式試驗中將加載車分3級加載，每一級加載持荷2 ～ 3 min，待實測應變及撓度數(shù)據(jù)穩(wěn)定后進行數(shù)據(jù)的采集工作，再進行下一級的加載工作. 限于篇幅，僅就其中一個工況做詳要概述. 為了使加載截面承受最大正彎矩，在中跨跨中截面中心加載，相應的荷載效率系數(shù)為0.967.撓度測點布置如圖11所示.

實橋試驗的撓度采用水準儀進行測量，現(xiàn)場測試照片如圖12所示.脈動試驗采用多通道數(shù)據(jù)采集分析系統(tǒng)NI公司的PXI系統(tǒng)進行試驗，采用8330B3型超低頻加速度傳感器進行數(shù)據(jù)采集，最低采樣頻率從0 Hz開始，采集主梁橫向和豎向振動數(shù)據(jù)，再經(jīng)信號分析得到全橋的各階固有振動特性，現(xiàn)場采集照片如圖13所示.

2.3? ?基于響應面法的有限元模型更新

基于響應面法的有限元模型修正是用響應面函數(shù)來模擬實際結(jié)構(gòu)的響應函數(shù)，將試驗設計與數(shù)理統(tǒng)計相結(jié)合，通過樣本選取、方差分析參數(shù)選取、響應面的擬合及采用優(yōu)化算法尋求響應面模型中的最優(yōu)解來進行有限元模型修正.

據(jù)文獻[15-16]知，影響有限元模型與實測橋梁結(jié)構(gòu)差別的主要因素為結(jié)構(gòu)混凝土密度、主梁和墩混凝土彈性模量，因此取這三個因素作為變量. 為獲得響應與所選定的三參數(shù)之間的聯(lián)系，首先需要參數(shù)設計. 根據(jù)參數(shù)取值的變化規(guī)律，假設三參數(shù)的單位長度值為10%，則可得該三參數(shù)值的變化范圍見表3.

以橋梁模態(tài)和靜力變形為目標來更新橋梁模型，選取橋梁結(jié)構(gòu)的第一階自振頻率（R1）、測點5撓度值（R2）和測點6撓度值（R3）作為目標函數(shù).這3個目標函數(shù)充分利用了現(xiàn)場的實橋試驗條件，且包含了橋梁靜、動力性能指標，可較為全面、準確地反映橋梁的力學性能.

根據(jù)上述試驗設計，應用回歸分析技術(shù)對樣本數(shù)據(jù)進行響應面擬合，并用R2準則和R2adj準則進行響應面擬合精度的檢驗，其表達式見方程（1）（2）.

式中：SST = SSE + SSR表示模型的總方差;dT表示模型的總自由度.如果R2和R2adj的值都接近1且兩者差值很小，則表示響應面方程擬合得很好.

采用二次多項式對上述樣本結(jié)果進行響應面擬合，并用R2和R2adj準則對響應面的精度進行檢驗，檢驗結(jié)果如表4所示.

從表4知，R2和R2adj值均接近1，且兩者差值很小，說明模型擬合精度較高，能很好地表示各參數(shù)與目標函數(shù)的關(guān)系.

響應面函數(shù)采取省略交差項影響的二次多項式，則響應面擬合方程如方程（3）～（5）所示：

R1 = 0.821 9 - 0.000 3X1 + 0.012 7X2 +

0.012 0X3 - 0.000 1X22 - 0.000 1X23? ? ? ?（3）

R2 = 9.764 8 + 0.005 0X1 - 0.632 7X2 +

0.067 9X3 + 0.008 7X22 - 0.001 5X23? ? ? ?（4）

R3 = 16.557 2 + 0.008 6X1 - 1.050 3X2 +

0.091 6X3 + 0.014 4X22 - 0.002 3X23? ? ? ?（5）

該橋的第一階實測頻率為0.867 Hz;實橋跨中試驗加載的情形下，測點5、6的撓度值分別為5.687 mm和9.679 mm. 根據(jù)實測值和數(shù)值模擬值可以建立一個目標函數(shù)：

式中：R（i）為根據(jù)響應面法計算的響應;M（i）為實測值;coef（i）為權(quán)重系數(shù)，此處取為[1.5 1 1].目標函數(shù)可以利用Matlab中的遺傳算法（Genetic Algorithm）工具箱優(yōu)化.更新后的上述三參數(shù)分別如表5所示.

從表5可得，主梁和墩的混凝土彈性模量和密度都有所增加，說明依托工程的施工質(zhì)量較為可靠，同時參數(shù)的變動范圍仍在工程實踐的常見范圍之內(nèi)，說明了修正結(jié)構(gòu)的可靠性.利用所更新的參數(shù)計算得到大橋的第一階振動頻率及上述工況的試驗加載測點撓度如表6所示.

3? ?振動響應對比分析

3.1? ?現(xiàn)場測點布置及荷載工況

該連續(xù)剛構(gòu)橋主要對主箱梁在如圖14所示截面1、2處（即邊跨跨中、中跨跨中）的動位移與應變變化情況進行了分析，全橋共計布設動位移測點5個，應變測點18個.

現(xiàn)場布置了3種荷載工況：工況1，一輛35 t重車置于橋縱向中軸線上;工況2，兩輛35 t重車并排對稱置于橋縱向中軸線兩側(cè)，兩車橫向間距為1.3 m;工況3，四輛35 t重車以橋縱向中軸線對稱放置，車輛橫向間距為1.3 m，縱向間距為4 m. 所有工況，車輛均以40 km/h的速度從橋外駛?cè)氩⑼ㄟ^橋梁.限于篇幅，僅列出工況3的車輛布置圖，如圖15所示.

3.2? ?結(jié)果對比分析

3.2.1? ?動位移對比

將上述各個工況進行對比分析，邊/中跨跨中的位移實測值與計算值對比如圖16和圖17所示.

由圖16和圖17可知，計算值與實測值隨時間變化的波動規(guī)律一致，但存在同一時刻計算值與實測值吻合不是很好的情況，且計算值的曲線波峰稍大于實測值，兩者靜態(tài)豎向位移差值在1.71%～6.55%之間，振動位移（消除靜態(tài)位移值）幅值差值在3.71%～6.35%之間.其誤差主要原因可能是：1）橋梁模型雖經(jīng)過合理的參數(shù)更新，但與實測橋梁結(jié)構(gòu)特性還有一定差別;2）實橋試驗受車輛速度控制、行駛軌跡、儀器誤差等因素的影響，測量數(shù)據(jù)精度有較大的誤差.

3.2.2? ?動應變對比

將上述各個工況進行對比分析，邊、中跨跨中的動應變實測值與計算值對比如圖18和圖19所示.

由圖18和圖19可知，計算值與實測值兩者的變化規(guī)律一致，且數(shù)值模擬結(jié)果與實測數(shù)據(jù)吻合得較好.但計算值的曲線波峰稍大于實測值，兩者動應變的差值在3.35%～6.94%之間，振動應變（消除靜態(tài)應變值）幅值差值在2.83%～5.96%.

從上述對比情況可知，豎向振動位移與動應變的計算值與實測值隨時間的變化趨勢一致，且整體吻合較好，峰值接近，說明該法建立的橋梁模型和車輛模型是有效的，能應用于車橋耦合計算中.

4? ?輪胎氣壓對橋梁振動的影響

輪胎氣壓的大小與輪胎變形成幾何關(guān)系，因此，輪胎氣壓對跨中振動響應影響反映了輪胎變形對梁跨中振動的影響.以上述工況2為例，取3種分別為P1=0.4 MPa、P2=0.6 MPa、P3=0.8 MPa的不同氣壓，來分析車輪氣壓變化對中跨跨中動位移和動應變的影響.

從圖20和圖21可得，當車輪氣壓由0.6 MPa減少到0.4 MPa時，車輪柔性變大，與橋面的接觸面增大，橋梁中跨跨中動位移峰值減小3.6%，動應變減小2.7%;當車輪氣壓由0.6 MPa增加到0.8 MPa時，車輪剛度變大，與橋面的接觸面減少，橋梁中跨跨中動位移峰值增大4.4%，動應變增大3.85%.

5? ?結(jié)? ?論

本文基于LS-DYNA程序建立了車-橋相互作用模型.設立多種荷載工況，通過計算值與實測值的對比，驗證了車輛模型的有效性及該法建立車橋耦合模型的可行性，并進一步分析了車輪氣壓變化對車橋振動的影響.結(jié)果表明：1）本文方法可較真實再現(xiàn)車-橋耦合振動情形，無需自編復雜的車橋耦合分析程序，只需定義LS-DYNA軟件自帶的面面接觸關(guān)鍵字，就可以實現(xiàn)車橋間的耦合作用;2）可方便獲得橋梁任意部位的應力、應變及其他響應數(shù)據(jù);3）可高精度分析車輪氣壓變化及車輛其他參數(shù)對耦合系統(tǒng)振動影響.

參考文獻

[1]? ? 殷新鋒. 汽車荷載作用下梁式橋與斜拉橋的動態(tài)響應分析[D].長沙：湖南大學，2010：6—20.

YIN X F. Dynamic response analysis of girder bridge and cable-stayed bridge under vehicle load [D]. Changsha：Hunan University，2010：6—20. （In Chinese）

[2]? ? 夏禾. 車輛與結(jié)構(gòu)動力相互作用[M]. 北京：科學出版社，2005：154—164.

XIA H. Vehicle-structure interaction [M]. Beijing：Science Press，2005：154—164. （In Chinese）

[3]? ? CHEN S R，CAI C S. Accident assessment of vehicles on long-span bridges inwindy environments[J]. Journal of Wind Engineering and Industrial Aerodynamics，2004，92（12）：991—1024.

[4]? ? CAI C S，SHI X M，ARAUJO M. Effect of approach span condition on vehicle-induced dynamic response [J]. ASCE Journal of Bride Engineering，2007，14（4）：440—449.

[5]? ? HUANG D Z，WANG T L. Impact analysis of cable-stayed bridges[J]. Computers and Structures，1992，43（5）：897—908.

[6]? ? 韓萬水，王濤，李永慶，等. 基于模型修正梁格法的車橋耦合振動分析系統(tǒng)[J]. 中國公路學報，2011，24（5）：47—55.

HAN W S，WANG T，LI Y Q，et al. Analysis system of vehicle-bridge coupling vibration with grillage method based on model updating[J]. China Journal of Highway and Transport，2011，24（5）：47—55. （In Chinese）

[7]? ? 李奇，吳定俊，邵長宇. 考慮車體柔性的車橋耦合系統(tǒng)建模與分析方法[J]. 振動工程學報，2011，24（1）：41—47.

LI Q，WU D J，SHAO C Y. Modeling and dynamic analysis method of vehicle-bridge coupling system considering car-body flexibility[J]. Journal of Vibration Engineering，2011，24（1）：41—47. （In Chinese）

[8]? ? 鄧露，李樹征，淡丹輝，等. 橋梁動態(tài)稱重技術(shù)在中小跨徑混凝土梁橋上的適用性研究[J]. 湖南大學學報（自然科學版），2020，47（3）：89—96.

DENG L，LI S Z，DAN D H，et al. Study on applicability of bridge weigh-in-motion technology in short-to medium-span concrete girder bridges[J].Journal of Hunan University （Natural Sciences），2020，47（3）：89—96. （In Chinese）

[9]? ? 趙越，黃平明，劉修平，等. 基于等參映射與改進折半法的公路車橋耦合分析系統(tǒng)[J]. 湖南大學學報（自然科學版），2020，47（11）：29—37.

ZHAO Y，HUANG P M，LIU X P，et al. Highway vehicle-bridge interaction analysis system based on isoparametric mapping and improved binary search [J]. Journal of Hunan University （Natural Sciences），2020，47 （11）：29—37. （In Chinese）

[10]? 李小珍，劉德軍，晉智斌. 大跨度鐵路懸索橋車-線-橋耦合振動分析[J]. 鋼結(jié)構(gòu)，2010，25（12）：6—12.

LI X Z，LIU D J，JIN Z B. Analysis of train-track-bridge coupled vibration of a railway long-span suspension bridge[J]. Steel Construction，2010，25（12）：6—12. （In Chinese）

[11]? 朱志輝，王力東，龔威，等. 多種垂向輪軌關(guān)系的對比及改進的車-線-橋系統(tǒng)迭代模型的建立[J].中南大學學報（自然科學版），2017，48（6）：1585—1593.

ZHU Z H，WANG L D，GONG W，et al. Comparative analysis of several types of vertical wheel/rail relationship and construction of an improved iteration model for train-track-bridge system[J]. Journal of Central South University （Science and Technology），2017，48（6）：1585—1593.（In Chinese）

[12]? 譚國金，劉寒冰，程永春，等. 基于車-橋耦合振動的簡支梁橋沖擊效應[J]. 吉林大學學報（工學版），2011，41（1）：62—67.

TAN G J，LIU H B，CHENG Y C，et al. Analysis of impact of vehicle to simply supported beam bridge based on vehicle-bridge coupled vibration[J]. Journal of Jilin University （Engineering and Technology Edition），2011，41（1）：62—67. （In Chinese）

[13]? 李巖，吳志文，蔡明，等. 一種隨機車流與橋梁耦合振動的分析方法[J]. 哈爾濱工業(yè)大學學報，2018，50（3）：46—52.

LI Y，WU Z W，CAI M，et al. An analysis method for coupled vibration random traffic flow and bridge[J]. Journal of Harbin Institute of Technology，2018，50（3）：46—52. （In Chinese）

[14]? 陳水生，趙輝，桂水榮，等. 基于LS -DYNA公路橋車橋耦合的車輛模型研究[J]. 計算力學學報，2019，36（6）：747—756.

CHEN S S，ZHAO H，GUI S R，et al. Vehicle model research on vehicle bridge coupling of highway bridge based on LS -DYNA[J]. Chinese Journal of Computational Mechanics，2019，36（6）：747—756. （In Chinese）

[15]? DENG L. System identification of bridge and vehicle based on their coupled vibration[D]. Baton Rouge：Louisiana State University，2009：7—23.

[16]? YIN X F，F(xiàn)ANG Z，CAI C S，et al. Non-stationary random vibration of bridges under vehicles with variable speed[J]. Engineering Structures，2010，32（8）：2166—2174.