[摘? 要]：近年來(lái)，風(fēng)致行車(chē)事故較多。文章以橋上行駛車(chē)輛為研究對(duì)象，首先采用CFD數(shù)值模擬的方法計(jì)算車(chē)輛氣動(dòng)力，然后基于靜力簡(jiǎn)化的方法，計(jì)算得到車(chē)輛的行車(chē)臨界風(fēng)速，可以較快地對(duì)車(chē)輛的行駛安全進(jìn)行評(píng)價(jià)。研究結(jié)果表明：隨著風(fēng)偏角增大，車(chē)輛的側(cè)力系數(shù)、升力系數(shù)和側(cè)翻力矩系數(shù)呈現(xiàn)出先增大后減小的趨勢(shì)，側(cè)翻力矩系數(shù)相對(duì)其他兩項(xiàng)數(shù)值較小。隨著車(chē)速的增大，觀光巴士的側(cè)翻臨界風(fēng)速逐漸減小。

[關(guān)鍵詞]：斜風(fēng); 行車(chē)安全; 數(shù)值模擬; 靜力簡(jiǎn)化

U461.1A

風(fēng)致行車(chē)安全一直是國(guó)內(nèi)外的研究熱點(diǎn)，近年來(lái)，車(chē)輛在斜風(fēng)作用下發(fā)生的交通事故頻頻報(bào)道[1]。我國(guó)高速公路發(fā)展迅速，為跨越江河湖海及山區(qū)峽谷，修建了一大批大跨度橋梁。而大跨度橋梁的海拔通常較高，易受到大風(fēng)的侵襲，車(chē)輛在橋上行駛過(guò)程中容易發(fā)生危險(xiǎn)。因此，有必要對(duì)斜風(fēng)作用下橋上車(chē)輛的行駛安全進(jìn)行研究，以保證旅客的人身安全。韓萬(wàn)水等[2-3]建立了大跨度鋼桁梁風(fēng)-車(chē)-橋耦合系統(tǒng)。王少欽等[4]計(jì)算了鐵路懸索橋在列車(chē)及風(fēng)荷載作用下的振動(dòng)響應(yīng)，并分析了響應(yīng)極值的產(chǎn)生機(jī)理及變化規(guī)律。韓艷等[5]提出了一種可以考慮抖振力空間相關(guān)性的風(fēng)-車(chē)-橋耦合振動(dòng)分析方法，并將其制作成軟件程序。此外，還有很多學(xué)者都對(duì)風(fēng)-車(chē)-橋耦合振動(dòng)進(jìn)行了研究，并取得了一些研究成果[6-8]。但是這些研究通常要耗費(fèi)大量的計(jì)算資源，無(wú)法對(duì)行車(chē)臨界風(fēng)速進(jìn)行快速評(píng)價(jià)。本文以某橋上車(chē)輛為研究對(duì)象，采用CFD數(shù)值模擬計(jì)算車(chē)輛氣動(dòng)力，并基于靜力簡(jiǎn)化的方法，可以較快速地對(duì)車(chē)輛的行駛安全進(jìn)行研究。

1 工程背景

某橋梁斷面為流線型箱梁形式，如圖1所示，其寬32.0 m，高3.0 m。一輛觀光巴士在橋上行駛，車(chē)輛高度為3.3 m，寬度為2.5 m，長(zhǎng)度為10.4 m。為降低數(shù)值模型的建立難度，對(duì)觀光巴士的外形進(jìn)行了一定程度的簡(jiǎn)化，保留了車(chē)輛的主體外形，而忽略一些車(chē)輛的細(xì)部構(gòu)造，簡(jiǎn)化后的車(chē)輛幾何模型如圖2所示。

2 車(chē)輛氣動(dòng)力系數(shù)

車(chē)輛氣動(dòng)力系數(shù)是計(jì)算行車(chē)臨界風(fēng)速的基礎(chǔ)，可以通過(guò)CFD數(shù)值模擬計(jì)算出橋上車(chē)輛的氣動(dòng)力系數(shù)。在斜風(fēng)來(lái)流作用下，對(duì)車(chē)輛行駛安全影響最大的氣動(dòng)力主要有側(cè)向力FS、升力FL和側(cè)翻力矩MR，車(chē)輛受力示意如圖3所示。

車(chē)輛的氣動(dòng)力系數(shù)可定義為下式，其中CS、CL和CR分別代表車(chē)輛的側(cè)力系數(shù)、升力系數(shù)和側(cè)翻力矩系數(shù)，U0為來(lái)流風(fēng)速，Af為車(chē)輛正投影面積，L為車(chē)輛的長(zhǎng)度，ρ為空氣密度，本文中取1.225 kg/m3。

CS=FS0.5ρU20Af

CL=FL0.5ρU20Af

CR=MR0.5ρU20AfL

（1）

采用ICEM軟件進(jìn)行網(wǎng)格劃分，建立數(shù)值模型，整體計(jì)算區(qū)域如圖4所示，其中計(jì)算域尺寸為25B×15B，其中B為主梁寬度。主梁和車(chē)輛壁面設(shè)置多層貼壁層網(wǎng)格，首層網(wǎng)格高度為0.001 m，并對(duì)附近區(qū)域網(wǎng)格進(jìn)行加密?？拷诿娴木W(wǎng)格尺寸足夠小，可以保證了計(jì)算精度，遠(yuǎn)離壁面的網(wǎng)格尺寸逐漸增大，可以保證計(jì)算速度，網(wǎng)格總數(shù)控制在為280萬(wàn)左右，細(xì)部網(wǎng)格如圖5所示。迎風(fēng)側(cè)邊界設(shè)為Velocity-inlet;背風(fēng)側(cè)邊界設(shè)為Pressure-outlet;上下邊界設(shè)置為Symmetry;前后邊界視來(lái)流風(fēng)偏角而定;主梁表面和車(chē)輛表面設(shè)為Wall。選用SST k-ω湍流模型;用SIMPLEC算法解決動(dòng)量方程中速度分量和壓力的耦合問(wèn)題;動(dòng)量方程、湍動(dòng)能方程及湍流耗散率方程均采用二階離散格式。

鐵路與公路王妍： 基于靜力簡(jiǎn)化的斜風(fēng)作用下橋上車(chē)輛行駛安全研究

計(jì)算得到0～90°風(fēng)偏角來(lái)流下的車(chē)輛氣動(dòng)力系數(shù)如圖6所示。由圖6可知，車(chē)輛的氣動(dòng)力系數(shù)受風(fēng)偏角的影響明顯。隨著風(fēng)偏角增大，車(chē)輛的側(cè)力系數(shù)、升力系數(shù)和側(cè)翻力矩系數(shù)均呈現(xiàn)出先增大后減小的趨勢(shì)，側(cè)翻力矩系數(shù)相對(duì)其他兩項(xiàng)，數(shù)值較小。其中在0～45°范圍內(nèi)，側(cè)力系數(shù)的變化幅度很小。

3 行車(chē)臨界風(fēng)速

橋上車(chē)輛運(yùn)動(dòng)過(guò)程中同時(shí)受到風(fēng)荷載、重力和摩擦力的共同作用，車(chē)輛的受力示意如圖7所示。

風(fēng)致行車(chē)安全事故主要為側(cè)傾事故，主要與車(chē)輛受到的側(cè)力、升力、側(cè)翻力矩有關(guān)。值得注意的是，汽車(chē)在橋上行駛過(guò)程中，不僅受到自然風(fēng)引起的風(fēng)荷載，還有車(chē)輛運(yùn)動(dòng)引起的縱向風(fēng)荷載，二者的疊加才為車(chē)輛受到的實(shí)際風(fēng)荷載，記為Ure;相對(duì)風(fēng)向角記為θ，表達(dá)式如式（2）所示。車(chē)輛受到的合成風(fēng)荷載如式（3）所示。其中，U為自然風(fēng)的風(fēng)速，Uv為車(chē)輛行駛速度;CS（θ）、CL（θ）和CR（θ）分別為風(fēng)向角為θ時(shí)的車(chē)輛側(cè)力系數(shù)、升力系數(shù)和側(cè)翻力矩系數(shù)，可通過(guò)多項(xiàng)式擬合后進(jìn)行取值。

Ure=U2+U2v

θ=arctan（Uv/U）（2）

FS=12ρU2reAfCSθ

FL=12ρU2reAfCLθ

MR=12ρU2reAfLCRθ

（3）

根據(jù)D’Alembert’s原理引入慣性力將車(chē)輛運(yùn)動(dòng)學(xué)問(wèn)題簡(jiǎn)化為靜力學(xué)問(wèn)題處理，則車(chē)輛的橫向慣性力Fx與豎向慣性力Fy可表示為

Fx=-maxFy=-may（4）

式中：m代表車(chē)輛質(zhì)量，ax和ay分別代表車(chē)輛的橫向加速度和豎向加速度。

車(chē)輛受到的風(fēng)荷載的作用點(diǎn)與橋面有一定的距離，使得車(chē)輛頂部由迎風(fēng)側(cè)向背風(fēng)側(cè)傾斜，稱為總側(cè)翻力矩，它可分為氣動(dòng)側(cè)翻力矩、氣動(dòng)側(cè)向力和車(chē)輛橫向慣性力引起的力矩、氣動(dòng)升力與車(chē)輛豎向慣性力產(chǎn)生的力矩，如式（5）所示。而重力則產(chǎn)生一定的力矩與之平衡，使車(chē)輛保持安全行駛狀態(tài)，稱為抗側(cè)翻力矩，如式（6）所示。當(dāng)二者相等時(shí)，則處在車(chē)輛側(cè)翻的臨界狀態(tài)，此時(shí)自然風(fēng)風(fēng)速大小則為側(cè)翻臨界風(fēng)速Ur。

MR+（FS+max）hvcosφ+（FL+may）·0.5Bcosφ（5）

mg·（0.5Bcosφ+hvsinφ）

（6）

當(dāng)風(fēng)偏角為0°時(shí)，車(chē)輛向前方行駛，其相對(duì)風(fēng)向角為正值，對(duì)圖6中3種氣動(dòng)力系數(shù)進(jìn)行多項(xiàng)式擬合，不同車(chē)速與風(fēng)速下的氣動(dòng)力系數(shù)根據(jù)此公式進(jìn)行取值。

經(jīng)過(guò)計(jì)算，當(dāng)不同車(chē)輛以不同速度行駛時(shí)，行車(chē)臨界風(fēng)速如圖8所示。由圖可知，隨著車(chē)速的增大，觀光巴士的側(cè)翻臨界風(fēng)速逐漸減小。當(dāng)車(chē)速為100 km/h時(shí)，車(chē)輛的臨界風(fēng)速為29.6 m/s。在極端大風(fēng)天氣下，可能發(fā)生危險(xiǎn)。

4 結(jié)論

本文以斜風(fēng)來(lái)流下，在橋上行駛的觀光巴士為研究對(duì)象，研究了其行車(chē)安全性，得到結(jié)論：

（1）車(chē)輛的氣動(dòng)力系數(shù)受風(fēng)偏角的影響明顯。隨著風(fēng)偏角增大，車(chē)輛的側(cè)力系數(shù)、升力系數(shù)和側(cè)翻力矩系數(shù)均呈現(xiàn)出先增大后減小的趨勢(shì)，側(cè)翻力矩系數(shù)相對(duì)其他兩項(xiàng)，數(shù)值較小。

（2）隨著車(chē)速的增大，觀光巴士的側(cè)翻臨界風(fēng)速逐漸減小。當(dāng)車(chē)速為100 km/h時(shí)，車(chē)輛的臨界風(fēng)速為29.6 m/s。在極端大風(fēng)天氣下，可能發(fā)生危險(xiǎn)。

參考文獻(xiàn)

[1] Zhu L D， Li L， Xu Y L， et al. Wind tunnel investigations of aerodynamic coefficients of road vehicles on bridge deck[J]. Journal of Fluids and Structures. 2012， 30： 35-50.

[2] 韓萬(wàn)水，趙越，劉煥舉，等.風(fēng)-車(chē)-橋耦合振動(dòng)研究現(xiàn)狀及發(fā)展趨勢(shì)[J].中國(guó)公路學(xué)報(bào)，2018，31（7）：1-23.

[3] 韓萬(wàn)水，劉煥舉，包大海，等.大跨鋼桁梁懸索橋風(fēng)-車(chē)-橋分析系統(tǒng)建立與可視化實(shí)現(xiàn)[J].土木工程學(xué)報(bào)，2018，51（3）：99-108.

[4] 王少欽，馬骎，任艷榮，等.主跨1120 m鐵路懸索橋風(fēng)-車(chē)-橋耦合振動(dòng)響應(yīng)分析[J].鐵道科學(xué)與工程學(xué)報(bào)，2017，14（6）：1241-1248.

[5] 韓艷，陳浩，劉躍飛，等.橋梁抖振力空間相關(guān)性對(duì)風(fēng)-車(chē)-橋耦合動(dòng)力響應(yīng)的影響[J].湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2015，42（9）：82-88.

[6] 舒鵬.變高度雙層鋼桁結(jié)合梁風(fēng)-車(chē)-橋耦合振動(dòng)性能分析[J].四川建筑，2018，38（4）：134-136.

[7] 黃靜文. 基于風(fēng)-車(chē)-橋耦合振動(dòng)的大跨度懸索橋車(chē)輛行駛安全性和舒適性研究[D].長(zhǎng)沙：長(zhǎng)沙理工大學(xué)，2018.

[8] 蘇波，齊冠，唐詩(shī)思，等.基于CFD的風(fēng)-車(chē)-橋系統(tǒng)耦合氣動(dòng)力模擬研究[J].河南科學(xué)，2017，35（5）：778-786.