隨機車流作用下橋梁沖擊系數分析

殷新鋒，鄧 露

（1．長沙理工大學 土木與建筑學院，湖南 長沙 410004； 2．湖南大學 土木工程學院，湖南 長沙 410082）

單一移動車輛對橋梁結構產生的振動響應問題研究已經較為成熟．以往的車－橋耦合振動的分析中，大多數學者選擇單個的三維車輛模型［1－4］，或者僅僅考慮一個確定的車列荷載［3－5］．近年來，雖少數研究者開始將研究重點轉移到隨機車流作用下橋梁振動響應分析，如韓萬水，Chen S．R．等，但還未形成統一結論［5－9］，且針對沖擊系數的研究仍然相對較少．沖擊系數作為一個重要參數被廣泛運用在橋梁設計中，因此，隨機車流作用下橋梁的沖擊系數分析顯得尤為重要．

由于隨機車流本身存在較大的不確定性，使得沖擊系數的研究較為困難，主要表現為：首先，由于多車道、多向行駛以及各個車輛的行駛速度不同導致車輛間的相對位置不斷變換，使得最大靜撓度的確定十分困難；其次，隨機車流中車輛數較多且類型復雜，如車輛均采用三維模型模擬，計算最大動撓度時，則計算量將較大．

本文提出了簡便且實用的隨機車流下橋梁結構沖擊系數計算方法：首先根據響應面分析方法，擬合影響面的函數表達式，將確定最大靜位移時隨機車流布載位置的問題轉變為求函數的最大值問題，從而獲得隨機車流作用下橋梁結構最大靜撓度；其次在計算最大動撓度時，根據已編制的車－橋耦合振動分析程序和隨機車流中各類車型的三維模型，獲得各類車型單獨作用下車－橋耦合接觸力并等效代替各類車輛三維模型作用在橋梁結構上，這樣便可將三維車輛模型作用等效轉化為隨機接觸力作用在橋梁結構，從而獲得隨機車流作用下橋梁結構的動力效應，基于統計分析，求得隨機車流作用下橋梁結構沖擊系數．

1 隨機車流模擬

1．1 公路橋梁隨機車流分類

按照全國高速公路管理部門制定的收費標準，根據車輛裝載貨物或乘客的能力大小將所有車輛劃分為12小類，具體分類見文獻［7］．

1．2 隨機車流數據分析

公路車流擁有較強的隨機性，通過對實測數據進行分析整理，下面僅給出隨機車流主要參數按時間段進行分類的統計分析結果．限于文章篇幅，僅給出車型、車道及車重分析．

1．2．1 車型分析

本次調查歷時24h，共采集到13 874輛的車流數據，將實測車流數據按本文的車型分類并進行分類統計，具體數據見表1．

表1 車型占有率Tab．1 Occupancy of vehicles

1．2．2 車道分析

車輛作用在橋梁橫向上的不同位置會得到不同的動力荷載效應，所以很有必要考慮車流在車道上的分布．根據實際測量的數據，選取某一時刻下，C1車型到C12 車型在超車道、第一行車道、第二行車道上的分布比例進行分類統計，如圖1所示．

圖1 24hC1～C12車型車道比例Fig．1 Ratio of from C1to C12vehicle models in 24hours

1．2．3 車重分析

限于文章篇幅，僅對C1車型，9：00～10：00在不同車道下的車重數據進行分析．通過橋梁動態測試儀，利用車橋耦合程序對車輛的實際車重進行識別，結合分類后的車型進行數據的篩選．通過K－S檢驗法，將統計所得結果，對常見的幾類分布函數進行判別，結果如圖2～圖3所示．從上述數據可知C1車型的車重在各車道下服從極值I型分布．

1．3 隨機車流模擬

通過Matlab軟件編寫出隨機車流程序分析并進行抽樣，即可獲得如圖4所示隨機車流模擬圖．圖4中不同標志代表不同車型，x坐標表示縱向位置，y坐標表示車道，z坐標表示車重，每點上的數值表示車速．

圖2 C1車型在超車道時車重分析Fig．2 Weight analysis of C1vehicle in overtaking lane

圖3 C1車型在行車道一時車重分析Fig．3 Weight analysis of C1vehicle in lane one

圖4 隨機車流模擬Fig．4 Simulation of random traffic flow

2 隨機車流下最大靜撓度分析

求解隨機車流下最大靜撓度的問題關鍵在于怎樣獲得最大靜撓度下對應車輛的布載情況．主要思路為：根據響應面分析方法，擬合影響面的函數表達式，將確定隨機車流布載位置的問題轉化為求對應函數表達式的最大值問題．下面以斜拉橋為例，詳細說明過程．

2．1 有限元軟件建模及影響面求解

本文選取某斜拉橋為研究對象，該橋為斜拉橋，全長840m（210m＋420m＋210m）．

采用機動法進行影響面的求解，選擇跨中截面點作為研究對象，將包含位置節點在內的單元作為強迫應變盒［10］，對單元上的4個節點各施加豎直向上的單位強制位移，為了突出豎向變形形狀，將顯示比例調整到25倍，形成了如圖6所示曲面圖形．由圖6可知，豎向最大位移為1．025m，出現在強制位移處，由于受到橋梁支座約束的影響，最小位移為－0．177m，出現在邊跨跨中位置．將相應強迫應變產生的節點反力反向施加在對應的節點上，形成了如圖7所示變形曲面．對比兩圖的形狀以及數據可知，變形形狀十分接近，最大值相差0．001，僅為變形量的0．975%；最小值相差則完全吻合．

圖6 強制位移Fig．6 Forced displacement

圖7 等效荷載位移Fig．7 Equivalent load displacement

2．2 影響面的擬合

本文采用多項式擬合法進行擬合，在擬合中，其正規方程組往往是病態的，而且正規方程組系數矩陣的階數越高，病態越嚴重；擬合節點分布的區間［x0，xm］偏離原點越遠，病態越嚴重；xi（i＝0，1，…，m）的數量級相差越大，病態越嚴重．針對上述問題，對模型數據做出以下調整：將縱、橫橋向坐標進行平移處理；根據影響面的形狀特點，將橋梁結構劃分成（50m＋23m＋24 m＋23m＋50m）共5個區段；選擇x，y均為2階函數進行擬合．利用Matlab軟件分段擬合，如圖8～圖13所示．擬合公式如下：

圖8 第一部分擬合曲面Fig．8 The first part of the fitting surface

圖9 第二部分擬合曲面Fig．9 The second part of the fitting surface

圖10 第三部分擬合曲面Fig．10 The third part of the fitting surface

圖11 第四部分擬合曲面Fig．11 The fourth part of the fitting surface

圖12 第五部分擬合曲面Fig．12 The fifth part of the fitting surface

圖13 分段函數圖形Fig．13 Piecewise graph

由圖13可知，橋面節點與擬合后曲面吻合較好，由確定系數可知，所有的擬合曲面都大于0．96，部分曲面甚至達到0．99．所以認為擬合效果能夠滿足計算精度的要求．

2．3 最大靜撓度求解

通過上述擬合處理獲得影響面表達式，布載位置的確定就轉變為函數求解最大值的問題，從而使得復雜的問題簡單明了．具體表達式的推導如下：

1）t時刻，第i輛車的第j輪對應的f（x，y）值表達如下：

式中：Li為第i＋1輛車與第i輛車之間的間距．

2）t時刻，第i輛車對影響面點的位移效應值如下式：

式中：k為i車的車軸數；Fj為第i輛車對應的輪重．Fj通過隨機生成的車重數據，按照車重對車軸的分配比率分配到各個軸，然后將軸重平均值賦予對應位置的輪胎重力Fj．

式中：Ll為車輪間距；Lz為車軸間距；F1＝F2，F3＝F4．t時刻，橋梁某車道上車輛對影響面點的效應之和為：

考慮到實際橋梁結構上存在車輛多向行駛的問題，因此需要對式（5）進行調整，當某一車道上車輛往順橋向行駛時；當車輛為逆橋向行駛時．這樣就可解決橋梁結構上雙向行駛的問題，而且不受車道數的影響，可以逐個累加．其中：S為橋梁跨徑，vi為第i輛車的速度，t時刻，橋梁上車輛對影響面點的效應之和為：

通過Matlab對P值進行循環計算，計算步長Δt＝0．1s，然后對所求的P值進行篩選，獲得Pmax．

3 隨機車流下的車 －橋耦合振動分析程序實現

本文考慮隨機車流中車輛數較多，如果所有車輛均模擬成三維車輛模型，程序的計算過程及計算難度會大大增加．為簡化，計算過程分兩步：首先得到車流樣本中各類型車輛各自單獨作用下與橋梁耦合振動的接觸力；然后在隨機車流橋梁振動分析時，將各自振動接觸力代替各類型車輛作用在橋上．鑒于篇幅原因，下面僅以重為20t的雙軸貨車（第三類型）為例，介紹得到單一車輛作用下與橋梁耦合振動的接觸力時程的計算過程．

3．1 單一車輛模型及其運動方程

采用7自由度兩軸汽車的空間整車模型模擬20 t的大型雙軸貨車（第三類型），自由度包括車體的豎向位移、點頭及翻滾轉動，前后車軸的豎向位移及翻滾轉動．該車輛模型雖然為最簡單的整車模型，卻能全面地分析車輛與橋梁的耦合振動特性，能比較真實地模擬實際車輛的動力特性．圖14所示為整車模型的正視圖和側視圖．

圖14 整車模型Fig．14 Vehicle model

根據文獻［1］可知，該車輛模型的運動微分方程矩陣形式為：

式中：Mv，Cv和Kv分別為車輛的質量矩陣、阻尼矩陣和剛度矩陣；Yv為車輛的位移向量；FG為車輛的重力向量；Fv－b為車輪與路面接觸面作用于車輪的反力．

3．2 橋梁運動方程

根據文獻［1］可知，橋梁的運動方程可以表示為：

式中：Mb，Cb和Kb分別表示橋梁的質量矩陣、阻尼矩陣和剛度矩陣；Yb表示橋梁全部自由度的位移向量；˙Yb和¨Yb分別表示Yb關于時間的一階和二階導數；Fv－b是作用于橋梁的所有外力向量．

一般將路面的不平度作為空間頻率函數，視為零均值的平穩隨機過程［8－9］．通過傅立葉變換，由此譜密度函數應用三角級數模擬得到路面不平度函數［1］：

式中：Gq（n）＝Gq（n0）n20／n2，n為空間頻率，n0＝0．1（1／m）是標準空間頻率，Gq（n0）為路面不平度系數；θi為0～2π之間均勻分布的隨機數．

3．3 單一車輛作用下橋梁振動響應分析

以第三類Ⅵ型為例，總重為20t的雙軸貨車距離橋頭15m 以速度10m／s行駛通過該橋．通過車橋耦合程序，得到車輛通過橋梁時橋梁的振動響應．

由圖15可知，車輛上橋后行駛在邊跨位置時，跨中位移僅出現較小的波動．15s之后逐漸增加．在22．88s時刻，即車輛行駛在跨中位置時位移達到最大值－4．09cm，隨著車輛離開跨中位置，豎向位移逐漸減小，44s位置車輛離開橋梁，跨中位移仍然在小范圍內發生波動．

圖15 單車作用下車－橋耦合系統振動響應Fig．15 Vibration response of vehicle－bridge coupled system

跨中截面豎向加速度則是隨著車輛上橋逐漸增大．行駛到22．92s時刻，豎向加速度達到最大值0．131 7m／s2．隨著車輛離開跨中位置，加速度逐漸減?。嬎愕煤筝S重為12．1t，在考慮車橋耦合作用下，后軸的接觸力如圖15（c）所示，最大值出現在30．02s位置，達到12．83t．

3．4 隨機車流作用下橋梁振動響應分析

如圖15（c）所示，基于車橋耦合程序可得考慮單車－橋耦合振動效應后的單車作用在橋梁時車軸接觸力時程數據值．在隨機車流橋梁振動分析時，為簡化計算過程，可將各類車單獨作用在橋上時所得時程接觸力數據值代替三維車輛模型直接作用在橋梁結構上［9］．為與單車作用下橋梁振動對比研究，車流總車重約為20t，則車流樣本作用下的振動響應如圖16所示．

圖16 車流樣本作用下橋梁振動響應Fig．16 Vibration response of the bridge under the traffic loads

跨中豎向位移在1～10s期間有較明顯的增加趨勢，由于受到的是車流影響，豎向位移與單一車輛相比規律性不明顯．在考慮60s內的豎向位移時，最大位移為3．647cm，發生在22．24s位置；對比單車作用下最大位移值可知，車流作用下最大位移值要小，這可能是由于車流中部分車輛作用在跨中截面的負彎矩區域，使跨中位置產生了反向位移，從而減小了車流作用下跨中最大位移．

4 隨機車流作用下橋梁的沖擊系數統計分析

本文通過大量的計算獲得一系列分散數據，再對數據進行統計分析從而獲得對應的分布函數和統計參數．借助Matlab工作平臺，采用K－S檢驗法對數據進行常規的5種分布函數檢驗．若分布函數與經驗函數差值的絕對值的最大值小于臨界值，則表示假設成立．由表2可知，經過K－S檢驗得出其服從極值Ⅰ型分布和正態分布，但是比較觀測值和臨界值可以發現，數據與極值Ⅰ型分布的吻合程度要高于正態分布（如圖17所示），所以認為沖擊系數服從極值Ⅰ型分布．

本文根據獲得的分布函數求解0．9 分位數［10］的沖擊系數值為0．161 5．與按照04規范求解得到的沖擊系數值0．149相比，隨機車流作用下所得沖擊系數要大8．38%，雖然04 規范與本文計算值有差別，但若從工程設計角度出發，04規范還是可運用于斜拉橋的設計中．

表2 沖擊系數檢驗Tab．2 Test of impact factor

圖17 沖擊系數擬合圖Fig．17 Fitted curve of impact factor

5 結 論

現有沖擊系數的研究較少考慮到交通車流的隨機特征，這與實橋上作用的交通荷載不符．因此，本文提出了簡便且實用的隨機車流作用下橋梁結構沖擊系數計算方法，結果表明：

1）基于響應面分析方法，提出了隨機車流作用下橋梁結構最大靜撓度的計算方法，算例表明該方法具有足夠的計算精度，滿足計算精度的要求．

2）根據隨機車流中各類車型的三維模型和已編制的車－橋耦合振動分析程序，提出了隨機車流作用下橋梁結構沖擊系數的計算方法，算例表明所提出的方法能有效地計算隨機車流下橋梁的沖擊系數．

3）沖擊系數統計分析表明，隨機車流作用下沖擊系數服從極值Ⅰ型分布，統計值比04規范求解所得的0．149大8．38%．

［1］殷新鋒．汽車荷載作用下梁式橋與斜拉橋的動態響應分析［D］．長沙：湖南大學土木工程學院，2010：20－30．

YIN Xin－feng．Dynamic response analysis of beam bridge and cable－stayed bridge under vehicle loads［D］．Changsha：College of Civil Engineering，Hunan University，2010：20－30．（In Chinese）

［2］夏禾．車輛與結構動力相互作用［M］．北京：科學出版社，2005：154－164．

XIA He．Vehicle－structure interaction［M］．Beijing：Science Press，2005：154－164．（In Chinese）

［3］DENG L，CAI C S．Development of dynamic impact factor for performance evaluation of existing multi－girder concrete bridges［J］．Engineering Structures，2010，32（1）：21－31．

［4］DENG L，CAI C S，BARBATO M．Reliability－based dynamic load allowance for capacity rating of pre－stressed concrete girder bridge［J］．Journal of Bridge Engineering，2011，16（6）：872－880．

［5］姚成釗．考慮橋面不平順影響的公路車流－橋梁耦合振動分析研究［D］．長沙：中南大學土木工程學院，2008：40－50．

YAO Cheng－zhao．Consider deck irregularity affecting road traffic－bridge coupled vibration analysis［D］．Changsha：Col－lege of Civil Engineering，Central South University，2008：40－50．（In Chinese）

［6］韓萬水，陳艾榮．隨機車流下的風－汽車－橋梁系統空間耦合振動研究［J］．土木工程學報，2008，41（9）：97－102．

HAN Wan－shui，CHEN Ai－rong．Three dimensional coupling vibration of wind－vehicle－bridge systems under random traffic flow［J］．China Civil Engineering Journal，2008，41（9）：97－102．（In Chinese）

［7］彭獻，殷新鋒，方志．變速車輛與橋梁的耦合振動及其TMD控制［J］．湖南大學學報：自然科學版，2006，33（5）：61－66．

PENG Xian，YIN Xin－feng，FANG Zhi．Vibration and TMD control of coupled system of girder bridge and vehicle with variable speeds［J］．Journal of Hunan University ：Natural Sciences，2006，33（5）：61－66．（In Chinese）

［8］周勇軍，趙煜，賀全海，等．剛構－連續組合橋梁沖擊系數多因素靈敏度分析［J］．振動與沖擊，2012，31（3）：97－101．

ZHOU Yong－jun，ZHAO Yu，HE Quan－hai，etal．Muti－parameters sensitivity analysis of impact factors for rigid－continuous combined bridge［J］．Journal of Vibration and Shock，2012，31（3）：97－101．（In Chinese）

［9］CHEN Su－ren．Dynamic performance of bridges and vehicles under strong wind［D］．Baton Rouge：Louisiana State University，2004：60－74．

［10］肖汝誠．空間結構內力影響面的有限元計算方法［J］．計算結構力學及其應用，1990（4）：56－63．

XIAO Ru－cheng．Spatial structure of internal force influence surface of the finite element method［J］．Computational Structural Mechanics and Applications，1990（4）：56－63．（In Chinese）