大攻角下典型主梁斷面顫振臨界風(fēng)速數(shù)值模擬

劉志文 陳岳飛 陳政清

摘? ?要：針對(duì)大攻角下主梁斷面顫振穩(wěn)定性問題，基于計(jì)算流體動(dòng)力學(xué)軟件ANSYS FLUENT用戶自定義函數(shù)UDFs和動(dòng)網(wǎng)格技術(shù)，結(jié)合Newmark - β法建立了橋梁主梁斷面二維流固耦合分析方法，采用該方法在不同攻角下（0°、±3°、±5°、±8°）對(duì)薄平板和流線型箱梁斷面顫振穩(wěn)定性進(jìn)行數(shù)值模擬研究，并將數(shù)值模擬結(jié)果與風(fēng)洞試驗(yàn)研究結(jié)果進(jìn)行了比較. 結(jié)果表明：薄平板和流線型箱梁斷面的顫振臨界風(fēng)速數(shù)值模擬結(jié)果與風(fēng)洞試驗(yàn)結(jié)果吻合較好，驗(yàn)證了本文所建立的橋梁主梁斷面二維流固耦合分析方法的精度. 薄平板斷面的顫振臨界風(fēng)速隨攻角的增大顯著降低;流線型箱梁斷面在正攻角范圍內(nèi)顫振臨界風(fēng)速隨著攻角的增大而降低，在負(fù)攻角范圍內(nèi)顫振臨界風(fēng)速隨攻角絕對(duì)值的增加先增大后降低;當(dāng)攻角較大時(shí)，薄平板斷面和流線型主梁斷面均表現(xiàn)出“鈍體”特征，氣流繞過斷面前緣時(shí)發(fā)生分離，沿?cái)嗝嫔舷戮壆a(chǎn)生較大的渦，從而導(dǎo)致主梁斷面顫振臨界風(fēng)速降低.

關(guān)鍵詞：典型主梁斷面;大攻角;顫振臨界風(fēng)速;流固耦合;數(shù)值模擬;風(fēng)洞試驗(yàn)

中圖分類號(hào)：U441.3? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1674—2974（2021）01—0082—10

Abstract：Aiming at the flutter stability of the main girder section at high attack angles， this paper studies the flutter stability of a thin plate and a streamlined box girder section at different attack angles（0°，±3°，±5°，±8°）by numerical simulation method. Two-dimensional fluid-structure interaction（FSI） method combined with Newmark-β method is developed based on ANSYS FLUENT user-defined function（UDFs） and dynamic mesh technology. Meanwhile， the numerical simulation results are compared with the wind tunnel test data. The results show that the critical flutter wind speeds of the thin plate section and the streamlined box girder section obtained from numerical simulation method are in good agreement with those from wind tunnel tests， which verifies the accuracy of the two-dimensional FSI analysis method for bridge deck section. The flutter critical wind speed of the thin plate section decreases significantly with the increase of attack angle. The flutter critical wind speed of the streamlined box girder section decreases with the increase of attack angle in the range of positive attack angles， and increases firstly and then decreases with the increasing of absolute values of attack angles in the range of negative attack angles. When the attack angle is large enough， the thin plate section and the streamlined deck section all characterize bluff body， and therefore the incoming flow passes around the leading edge of the section separates and generates large vortices along the upper and lower edges of the section， which leads to a decrease of the critical flutter wind speeds of the bridge girder section.

Key words：typical main deck sections;high attack angle;critical flutter wind speed;fluid-structure interaction;numerical simulation;wind tunnel tests

大攻角對(duì)顫振穩(wěn)定性的影響最早在機(jī)翼斷面顫振穩(wěn)定性研究中受到關(guān)注. 當(dāng)攻角增大到一定程度時(shí)，機(jī)翼斷面的氣動(dòng)力與飛行狀態(tài)參數(shù)之間呈現(xiàn)明顯的非線性關(guān)系，線性理論模型已經(jīng)不再適用. 國(guó)內(nèi)外學(xué)者[1-2]對(duì)此進(jìn)行了深入研究，目前大攻角對(duì)機(jī)翼顫振性能影響的研究已經(jīng)較為成熟. 近年來，隨著橋梁跨度的不斷增大，橋梁主梁結(jié)構(gòu)在風(fēng)作用下的靜力變形會(huì)進(jìn)一步增大，導(dǎo)致主梁在高風(fēng)速條件下附加風(fēng)攻角增大;另外復(fù)雜山區(qū)地形及峽谷地形橋位風(fēng)場(chǎng)由于受局部地形影響，主梁高度處存在較大的來流風(fēng)攻角. 因此，大攻角下主梁斷面的顫振穩(wěn)定性問題逐漸受到關(guān)注.

關(guān)于大攻角下橋梁主梁斷面抗風(fēng)性能的研究主要集中在大攻角下主梁斷面靜氣動(dòng)力系數(shù)、顫振性能和渦振性能等方面的研究. 大攻角下主梁斷面靜氣動(dòng)力系數(shù)研究方面，Reinhold等[3-4]分別采用風(fēng)洞試驗(yàn)和數(shù)值模擬方法研究了薄平板和流線型箱梁斷面三分力系數(shù)隨攻角的變化，結(jié)果表明攻角對(duì)氣動(dòng)升力系數(shù)的影響最為顯著，升力系數(shù)在一定的攻角范圍內(nèi)會(huì)出現(xiàn)先增加后減小的趨勢(shì).

隨著橋梁跨度的增加，高風(fēng)速條件下主梁斷面附加攻角逐漸增大，部分學(xué)者開展了附加風(fēng)攻角對(duì)主梁斷面顫振性能的影響研究. 許福友等[5]在進(jìn)行蘇通大橋全橋氣彈模型的三維顫振分析過程中，發(fā)現(xiàn)接近顫振失穩(wěn)臨界狀態(tài)時(shí)主梁斷面的附加攻角達(dá)到了1.5°，若不考慮附加攻角的影響會(huì)高估該橋的顫振臨界風(fēng)速. 朱青等[6]通過彈性懸掛風(fēng)洞試驗(yàn)研究了附加攻角對(duì)象山港大橋鋼箱梁節(jié)段模型顫振臨界風(fēng)速的影響，從氣動(dòng)導(dǎo)數(shù)的角度分析了顫振臨界風(fēng)速降低的原因. 朱樂東等[7]研究了附加風(fēng)攻角對(duì)扁平箱梁顫振的影響，近似10%的攻角增量也會(huì)引起顫振風(fēng)速的顯著變化. 張宏杰等[8]采用風(fēng)洞試驗(yàn)的方法研究了附加攻角對(duì)鋼箱梁斜拉橋顫振穩(wěn)定性影響，結(jié)果表明，考慮附加攻角影響的顫振分析得到的結(jié)果更為合理. 歐陽克儉等[9]采用強(qiáng)迫振動(dòng)試驗(yàn)裝置研究了附加攻角對(duì)洞庭二橋桁架加勁主梁斷面初步設(shè)計(jì)方案顫振臨界風(fēng)速的影響，結(jié)果表明，附加攻角效應(yīng)會(huì)降低桁架加勁梁的顫振臨界風(fēng)速. 劉志文等[10]采用自由振動(dòng)風(fēng)洞試驗(yàn)結(jié)合數(shù)值模擬的方法研究了伶仃洋大橋箱梁斷面的不同方案在0°、±3°和±5°下的顫振臨界風(fēng)速，結(jié)果表明各方案下的顫振臨界風(fēng)速在正攻角范圍內(nèi)均隨攻角的增大而降低，而在負(fù)攻角范圍內(nèi)隨攻角的增大而增大. 同時(shí)，風(fēng)攻角的增加還會(huì)直接改變橋梁結(jié)構(gòu)的顫振性能. 當(dāng)風(fēng)攻角增加時(shí)，主梁斷面形狀相對(duì)來流風(fēng)變鈍，顫振性能會(huì)由小攻角下的非線性較弱的發(fā)散型振動(dòng)（硬顫振）轉(zhuǎn)變?yōu)榇蠊ソ窍碌姆蔷€性自限幅極限環(huán)振動(dòng)（LCO），即軟顫振現(xiàn)象[11-15].

顫振導(dǎo)數(shù)是表征橋梁顫振性能的重要參數(shù)，通過強(qiáng)迫振動(dòng)風(fēng)洞試驗(yàn)，研究大攻角下不同斷面的顫振導(dǎo)數(shù)，對(duì)顫振穩(wěn)定性分析具有重要意義. 劉慕廣等[16]采用強(qiáng)迫振動(dòng)裝置分別對(duì)矩形斷面和H型斷面氣動(dòng)導(dǎo)數(shù)的識(shí)別方法進(jìn)行了試驗(yàn)研究，結(jié)果表明在大攻角下有必要采取與結(jié)構(gòu)實(shí)際振動(dòng)一致的模態(tài)來識(shí)別氣動(dòng)導(dǎo)數(shù). Tang等[17-18]采用數(shù)值模擬方法分別研究了中央開槽平板斷面和桁架加勁梁斷面在不同攻角下的顫振性能，結(jié)果表明，中央開槽平板斷面的顫振導(dǎo)數(shù)A*2在6°攻角下（桁架加勁主梁斷面在8°攻角下）開始出現(xiàn)由負(fù)變正的趨勢(shì)，顫振形式由彎扭耦合顫振變?yōu)榕まD(zhuǎn)顫振;隨著攻角的增加，斷面的氣動(dòng)特性接近鈍體，顫振臨界風(fēng)速顯著降低. 王騎等[19]采用強(qiáng)迫振動(dòng)和自由振動(dòng)風(fēng)洞試驗(yàn)的方法，對(duì)薄平板在不同攻角下的顫振特性進(jìn)行了研究，結(jié)果表明在7°攻角下，當(dāng)折算風(fēng)速Vred > 15時(shí)，薄平板的顫振導(dǎo)數(shù)A*2變?yōu)檎担冶∑桨逶?°攻角內(nèi)采用Scanlan線性顫振導(dǎo)數(shù)計(jì)算的顫振臨界風(fēng)速仍是可靠的.

大攻角對(duì)橋梁渦振性能的影響也受到關(guān)注. 楊陽等[20]采用風(fēng)洞試驗(yàn)對(duì)大攻角下的渦激振動(dòng)性能進(jìn)行了研究，結(jié)果表明大攻角下扁平箱梁的渦振振幅和范圍明顯增大，對(duì)于橋址位于山區(qū)等容易發(fā)生大攻角地區(qū)的橋梁應(yīng)進(jìn)行大攻角試驗(yàn). 劉小兵等[21]研究了單箱梁在-5° ～ 5°攻角范圍內(nèi)渦激振動(dòng)性能的變化情況，研究表明隨著風(fēng)攻角的逐漸減小，單箱梁渦激振動(dòng)的風(fēng)速鎖定區(qū)間逐漸變短，最大振幅值逐漸變小，渦激振動(dòng)性能逐漸變好.

綜上可知，目前關(guān)于大攻角下橋梁主梁斷面的顫振性能研究主要采用風(fēng)洞試驗(yàn)方法進(jìn)行研究，且研究的攻角范圍大多數(shù)為-5° ～ 5°. 本文在已有研究文獻(xiàn)的基礎(chǔ)上采用數(shù)值方法進(jìn)行典型主梁斷面在大攻角條件下的顫振穩(wěn)定性研究，并對(duì)大攻角對(duì)主梁斷面顫振性能影響機(jī)理進(jìn)行研究.

1? ?典型斷面簡(jiǎn)介

分別針對(duì)寬高比B/D=22.5的薄平板（B為薄平板斷面寬度，D為薄平板斷面高度）和幾何縮尺比為1/50的大帶東橋（Great Belt East Bridge，GBEB）主梁斷面進(jìn)行大攻角下顫振性能研究. 主梁斷面具體幾何尺寸分別見圖1和圖2，對(duì)應(yīng)的結(jié)構(gòu)動(dòng)力特性參數(shù)見表1.

2? ?數(shù)值模擬方法介紹

2.1? ?流體控制方程

不可壓黏性流體的動(dòng)量方程（即N-S方程）表示為：

式中：t為時(shí)間;ρ為空氣密度，ρ = 1.225 kg/m3;ui、uj分別為速度分量;p為壓力;υ為運(yùn)動(dòng)黏性系數(shù).

2.2? ?湍流模型

大渦模擬（Large Eddy Simulation，LES）的基本假設(shè)為：1）動(dòng)量、能量、質(zhì)量及其他標(biāo)量主要由大渦輸送;2）流動(dòng)的幾何和邊界條件決定了大渦的特性，且流動(dòng)特性主要由大渦體現(xiàn); 3）小尺度旋渦受幾何和邊界條件影響較小，并且各向同性. 大渦模擬過程中，直接求解大渦，而利用亞網(wǎng)格尺度模型模擬小尺度紊流運(yùn)動(dòng)對(duì)大尺度紊流運(yùn)動(dòng)的影響.

大渦模擬的控制方程通過對(duì)N-S方程進(jìn)行空間濾波得到，即

另需要指出的是，由于旋渦具有空間三維特性，LES模擬需要進(jìn)行三維分析來反應(yīng)實(shí)際流動(dòng)，但陳艾榮等[22]和詹昊等[23]采用二維LES模擬分別對(duì)方柱和大帶東橋進(jìn)行研究后發(fā)現(xiàn)，二維和三維大渦模擬數(shù)值結(jié)果均與試驗(yàn)結(jié)果吻合較好，三維大渦計(jì)算結(jié)果整體略好于二維大渦模擬結(jié)果，但二維大渦模擬計(jì)算效率遠(yuǎn)高于三維大渦模擬. 綜上，本文采用二維LES模擬對(duì)薄平板和典型流線型箱梁斷面的顫振穩(wěn)定性進(jìn)行直接計(jì)算分析.

2.3? ?流固耦合分析方法

采用LES模型求解不可壓縮黏性流體N-S方程，結(jié)合動(dòng)網(wǎng)格技術(shù)將Newmark-β算法通過UDF嵌入Fluent中，建立主梁斷面二維流固耦合分析方法（Flutter solid interaction，F(xiàn)SI）[24].

以結(jié)構(gòu)豎向振動(dòng)為例來介紹具體的流固耦合分析方法，其中結(jié)構(gòu)振動(dòng)方程為：

式中：m為結(jié)構(gòu)單位長(zhǎng)度質(zhì)量（kg/m）;c為阻尼系數(shù)（N·s/m）;k為剛度（N/m）;F（t）為結(jié)構(gòu)單位長(zhǎng)度受到的豎向氣動(dòng)力，通過求解流體動(dòng)力學(xué)方程獲得.

采用Newmark-β法求解結(jié)構(gòu)振動(dòng)方程：

2.4? ?計(jì)算域及網(wǎng)格

取薄平板計(jì)算域?yàn)?6B × 10B，示意圖如圖3所示. 對(duì)應(yīng)的邊界條件為：計(jì)算域左側(cè)為速度入口邊界（Velocity inlet），計(jì)算域右側(cè)為速度出口邊界（Pressure outlet），計(jì)算域上、下側(cè)為對(duì)稱邊界（Symmetry），斷面邊界為固定壁面邊界條件（Wall）. 薄平板斷面計(jì)算風(fēng)攻角范圍（0°～8°）對(duì)應(yīng)的阻塞率范圍為0.44%～1.65%. 流線型箱梁斷面計(jì)算域?yàn)?6B × 10B，邊界條件設(shè)置與薄平板斷面計(jì)算邊界條件一致，在計(jì)算風(fēng)攻角范圍內(nèi)（-8°～+8°）對(duì)應(yīng)的阻塞率范圍為1.40%～2.31%，滿足阻塞率小于5%的要求. 為計(jì)算方便，通過改變速度入口邊界參數(shù)，模擬來流的不同攻角.

計(jì)算域網(wǎng)格劃分采用“剛體運(yùn)動(dòng)網(wǎng)格+動(dòng)網(wǎng)格+靜態(tài)網(wǎng)格”方案，其中剛體運(yùn)動(dòng)網(wǎng)格區(qū)域和靜態(tài)網(wǎng)格區(qū)域采用結(jié)構(gòu)化網(wǎng)格，動(dòng)網(wǎng)格區(qū)域采用非結(jié)構(gòu)網(wǎng)格. 薄平板與流線型箱梁斷面的局部網(wǎng)格劃分分別見圖4和圖5. 對(duì)應(yīng)的首層網(wǎng)格高度均為0.000 1 m，邊界層內(nèi)網(wǎng)格增長(zhǎng)率分別為1.2和1.018，對(duì)應(yīng)的無量綱高度Y+均小于1. 相應(yīng)斷面網(wǎng)格數(shù)量分別約為50 380和99 470.

2.5? ?無關(guān)性驗(yàn)證

1）網(wǎng)格無關(guān)性驗(yàn)證. 主梁斷面三分力系數(shù)定義如下：

式中：CH、CV和CM分別為體軸坐標(biāo)系下主梁斷面水平力系數(shù)、豎向力系數(shù)和力矩系數(shù);FH、FV和MT分別為主梁斷面水平力、豎向力和扭矩;U為來流風(fēng)速（m/s）;長(zhǎng)度L取1 m. 綜合考慮網(wǎng)格精度與計(jì)算時(shí)間，首先進(jìn)行了薄平板斷面和流線型箱梁斷面網(wǎng)格無關(guān)性驗(yàn)證.

圖6（a）給出了薄平板斷面不同網(wǎng)格方案，即TP-Mesh-1（網(wǎng)格數(shù)為34 650）、TP-Mesh-2（網(wǎng)格數(shù)為50 380）和TP-Mesh-3（網(wǎng)格數(shù)為92 680）的網(wǎng)格無關(guān)性驗(yàn)證結(jié)果;圖6（b）為流線型箱梁斷面不同網(wǎng)格方案，即GBEB-Mesh-1（網(wǎng)格數(shù)為68 420）、GBEB-Mesh-2（網(wǎng)格數(shù)為99 470）和GBEB-Mesh-3（網(wǎng)格數(shù)為148 520）的網(wǎng)格無關(guān)性驗(yàn)證結(jié)果. 由圖6可知，兩種斷面采用不同網(wǎng)格方案得到的體軸坐標(biāo)系下的三分力系數(shù)十分接近，為了兼顧求解效率和計(jì)算精度，分別采用TP-Mesh-2和GBEB-Mesh-2的網(wǎng)格方案進(jìn)行后續(xù)計(jì)算.

2）時(shí)間步無關(guān)性檢驗(yàn). 薄平板和流線型箱梁斷面分別采用TP-Mesh-2和GBEB-Mesh-2的網(wǎng)格方案，在0°攻角下進(jìn)行斷面時(shí)間步無關(guān)性驗(yàn)證. 圖7分別給出了兩種斷面在不同時(shí)間步長(zhǎng)下流固耦合計(jì)算的扭轉(zhuǎn)位移（α/（°））時(shí)程曲線和無量綱豎向位移（Y/D）時(shí)程曲線. 薄平板斷面（0.000 2 s、0.000 5 s和0.001 s）和流線型箱梁斷面（0.000 1 s、0.000 2 s、0.000 5 s和0.001 s）在不同時(shí)間步長(zhǎng)下計(jì)算得到的位移時(shí)程曲線基本一致. 薄平板在不同時(shí)間步長(zhǎng)下計(jì)算的扭轉(zhuǎn)頻率分別為2.48 Hz、2.59 Hz和2.56 Hz，時(shí)間步為0.000 5 s和0.001 s的結(jié)果更為接近，因此取時(shí)間步長(zhǎng)為0.000 5 s進(jìn)行計(jì)算. 流線型箱梁在不同時(shí)間步長(zhǎng)下計(jì)算的扭轉(zhuǎn)頻率分別為4.09 Hz、4.08 Hz、4.03 Hz和4.09 Hz，因此取時(shí)間步長(zhǎng)為0.000 2 s進(jìn)行計(jì)算.

2.6? ?參數(shù)設(shè)置

采用LES中的Smagorinsky-Lilly湍流模型，Smagorinsky常數(shù)Cs = 0.10;速度-壓力耦合求解采用SIMPLEC格式;空間離散梯度、壓力和動(dòng)量分別采用Green-Gause Node Based、Second Order和Bounded Central Differencing;時(shí)間離散采用Second Order Implicit格式.

3? ?節(jié)段模型顫振穩(wěn)定性試驗(yàn)

為驗(yàn)證數(shù)值模擬的精度，需進(jìn)行薄平板斷面和流線型箱梁斷面顫振臨界風(fēng)速的風(fēng)洞試驗(yàn)測(cè)試. 薄平板與流線型箱梁斷面節(jié)段模型斷面幾何尺寸分別見圖1和圖2，模型的長(zhǎng)度均為L(zhǎng) = 1.54 m. 模型骨架采用鋁合金制作而成，外衣分別采用優(yōu)質(zhì)木板和ABS板制作，以保證幾何外形和剛度的設(shè)計(jì)要求. 為了確保節(jié)段模型周邊流動(dòng)的二元性，在薄平板節(jié)段模型和流線型箱梁節(jié)段模型的端部設(shè)置了尺寸分別為780 mm×160 mm×2 mm和270 mm×103.5 mm×2 mm的端板（端部設(shè)圓弧倒角）.

采用彈性懸掛主梁節(jié)段模型進(jìn)行不同攻角條件下薄平板斷面和顫振臨界風(fēng)速試驗(yàn)研究. 根據(jù)表1所示的薄平板斷面和流線型斷面節(jié)段模型動(dòng)力特性參數(shù)，進(jìn)行節(jié)段模型風(fēng)洞試驗(yàn)參數(shù)設(shè)計(jì). 模型振動(dòng)位移響應(yīng)采用激光位移計(jì)測(cè)試，采用東華DH5920多通道信號(hào)采集系統(tǒng)進(jìn)行數(shù)據(jù)采集，采樣頻率為200 Hz. 分別進(jìn)行了不同攻角下的薄平板斷面（α = 0°，+3°，+5°，+8°）和流線型箱梁斷面（α=-8°，-5°，-3°，0°，+3°，+5°，+8°）的顫振穩(wěn)定性試驗(yàn). 圖8為薄平板斷面和流線型箱梁斷面節(jié)段模型風(fēng)洞試驗(yàn)照片.

4? ?結(jié)果分析

4.1? ?0°攻角下的顫振臨界風(fēng)速結(jié)果分析

首先將0°攻角下薄平板顫振臨界風(fēng)速Vcr和顫振臨界狀態(tài)扭轉(zhuǎn)頻率fα的流固耦合數(shù)值模擬結(jié)果（Num.）與顫振臨界風(fēng)速理論解及試驗(yàn)結(jié)果（Exp.）進(jìn)行對(duì)比，以驗(yàn)證數(shù)值模擬方法的精度，對(duì)比結(jié)果如表2所示. 0°攻角下薄平板斷面顫振臨界風(fēng)速和顫振發(fā)生扭轉(zhuǎn)頻率的數(shù)值模擬結(jié)果與理想平板顫振臨界風(fēng)速和顫振扭轉(zhuǎn)頻率的理論計(jì)算結(jié)果較為接近，相對(duì)誤差分別為-0.14%和-0.39%;且數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果也較為吻合，驗(yàn)證了本文所建立的主梁流固耦合數(shù)值模擬方法的計(jì)算精度. 圖9為薄平板斷面0°攻角、不同風(fēng)速下結(jié)構(gòu)豎向和扭轉(zhuǎn)振動(dòng)位移響應(yīng)時(shí)程曲線.

采用相同的方法得到0°攻角下GBEB橋流線型主梁斷面顫振臨界風(fēng)速結(jié)果匯總，具體如表3所示. 0°攻角下GBEB橋流線型主梁斷面數(shù)值模擬結(jié)果和風(fēng)洞試驗(yàn)結(jié)果分別為18.0 m/s和20.7 m/s，對(duì)應(yīng)的扭轉(zhuǎn)頻率分別為4.00 Hz和3.90 Hz，換算到實(shí)橋的風(fēng)速結(jié)果分別為36.0 m/s和41.4 m/s（風(fēng)速比為1 ∶ 2），對(duì)應(yīng)的扭轉(zhuǎn)頻率分別為0.160 Hz和0.156 Hz（頻率比為1 ∶ 25）. 由表3可知，GBEB橋流線型主梁斷面數(shù)值模擬結(jié)果與風(fēng)洞試驗(yàn)結(jié)果、有限體積法計(jì)算結(jié)果[25]、Walther離散渦計(jì)算值[26]和Poulsen風(fēng)洞試驗(yàn)結(jié)果[27]之間的相對(duì)誤差分別為-13.04%、-11.11%、-4.26%、0.0%.

4.2? ?不同攻角下顫振臨界風(fēng)速結(jié)果分析

分別采用流固耦合方法和風(fēng)洞試驗(yàn)方法，對(duì)兩種典型斷面的彈性懸掛節(jié)段模型進(jìn)行了不同攻角、不同風(fēng)速下風(fēng)致振動(dòng)響應(yīng)計(jì)算和試驗(yàn)研究. 表4給出了兩種典型斷面的數(shù)值模擬與風(fēng)洞試驗(yàn)之間的對(duì)比結(jié)果. 由表4可知，薄平板斷面顫振臨界風(fēng)速均隨攻角的增加而逐漸降低，不同攻角下薄平板斷面顫振臨界風(fēng)速的數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果總體較為接近，最大相對(duì)誤差為-7.8%;流線型主梁斷面在正攻角范圍內(nèi)的顫振臨界風(fēng)速隨攻角的增加而逐漸減小，負(fù)攻角范圍內(nèi)在-3°攻角時(shí)最大，-5°攻角時(shí)次之，-8°攻角時(shí)最低，接近0°攻角時(shí)的顫振臨界風(fēng)速，不同攻角下流線型箱梁斷面顫振臨界風(fēng)速的數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果最大相對(duì)誤差為-18.6%，且數(shù)值模擬計(jì)算結(jié)果普遍低于風(fēng)洞試驗(yàn)結(jié)果，這一現(xiàn)象及誤差范圍與蘇通大橋的研究結(jié)果相近[15].

圖10和圖11分別給出了兩種斷面在不同攻角下扭轉(zhuǎn)位移幅值根方差和豎向位移幅值根方差隨風(fēng)速變化的數(shù)值模擬結(jié)果，圖中箭頭代表該風(fēng)速下的位移響應(yīng)已經(jīng)發(fā)散，位移會(huì)繼續(xù)發(fā)展. 由圖10和圖11可以看出，薄平板在+5°攻角范圍內(nèi)的顫振特性均表現(xiàn)為“硬顫振”的特點(diǎn)，而在+8°攻角時(shí)卻表現(xiàn)為“軟顫振”的特點(diǎn). 流線型箱梁斷面在不同攻角下的顫振特性都表現(xiàn)出“軟顫振”的特點(diǎn)，表明主梁斷面“軟顫振”的發(fā)生與攻角和斷面形狀有關(guān).

值得指出的是，由圖10和圖11可知，薄平板斷面通過數(shù)值模擬得到的不同攻角下的顫振性能均表現(xiàn)為硬顫振特性，而試驗(yàn)結(jié)果在+8°攻角下由硬顫振轉(zhuǎn)變?yōu)檐涱澱? 其次，兩種斷面通過數(shù)值模擬方法得到的位移響應(yīng)幅值和通過風(fēng)洞試驗(yàn)得到的位移響應(yīng)幅值之間具有較大的誤差，即兩種方法得到的后顫振性能不同，尤其是流線型箱梁斷面（兩種方法得到的位移響應(yīng)幅值范圍均可在文獻(xiàn)中找到相似的結(jié)果[13，15]）.以上現(xiàn)象的存在是由于在數(shù)值模擬方法中考慮了主梁斷面附加攻角及氣動(dòng)力非線性的影響，而風(fēng)洞試驗(yàn)方法在考慮主梁斷面附加攻角和氣動(dòng)力非線性的基礎(chǔ)上，還考慮了支承彈簧位置變化引起的非線性阻尼因素，使主梁斷面在后顫振過程中表現(xiàn)出強(qiáng)非線性自限幅極限環(huán)振動(dòng)，由此造成了主梁斷面顫振后振動(dòng)響應(yīng)的數(shù)值模擬結(jié)果較風(fēng)洞試驗(yàn)結(jié)果偏大.

4.3? ?分析與討論

為了進(jìn)一步分析攻角對(duì)主梁斷面顫振臨界風(fēng)速的影響機(jī)理，圖12和圖13分別給出了薄平板斷面和流線型箱梁斷面在0°和+8°攻角下以顫振臨界風(fēng)速振動(dòng)時(shí)一個(gè)周期內(nèi)不同時(shí)刻（T/4、T/2、3T/4、T）的流場(chǎng)渦量圖. 由圖12和圖13可知，0°攻角下，薄平板（Vcr = 14.5 m/s）和流線型箱梁斷面（Vcr =18.0 m/s）僅在尾部存在明顯的旋渦脫落，前沿上下表面沒有明顯的前沿渦，僅渦量大小發(fā)生輕微變化. 而當(dāng)攻角增大至+8°時(shí)，兩種斷面表現(xiàn)出“鈍體”特征，薄平板（Vcr =7.0 m/s）和流線型箱梁斷面（Vcr =10.5 m/s）不僅在尾部存在嚴(yán)重的旋渦脫落，其前沿也出現(xiàn)了邊界層分離，以致前沿存在明顯的周期性變化的大尺度旋渦. 這種前沿渦的出現(xiàn)會(huì)導(dǎo)致橋梁結(jié)構(gòu)表面壓力分布不均，產(chǎn)生較大的扭矩，進(jìn)一步激發(fā)結(jié)構(gòu)的氣動(dòng)不穩(wěn)定狀態(tài)，導(dǎo)致結(jié)構(gòu)在較低的風(fēng)速下提前進(jìn)入顫振.

5? ?結(jié)? ?論

本文基于計(jì)算流體動(dòng)力學(xué)軟件ANSYS FLUENT用戶自定義函數(shù)UDF和動(dòng)網(wǎng)格技術(shù)，結(jié)合Newmark-β法建立了橋梁主梁斷面二維流固耦合分析方法，采用該方法對(duì)不同攻角下（0°、±3°、±5°、±8°）薄平板和流線型箱梁斷面顫振臨界風(fēng)速進(jìn)行數(shù)值模擬研究，并與風(fēng)洞試驗(yàn)研究結(jié)果進(jìn)行比較，得到如下主要結(jié)論：

1）薄平板斷面和流線型箱梁斷面顫振臨界風(fēng)速數(shù)值模擬結(jié)果與風(fēng)洞試驗(yàn)結(jié)果及已有文獻(xiàn)結(jié)果吻合較好，表明本文所建立的橋梁斷面流固耦合分析方法的精度滿足要求.

2）在正攻角范圍內(nèi)，薄平板斷面和流線型箱梁斷面顫振臨界風(fēng)速隨著攻角的增加顯著降低;在負(fù)攻角范圍內(nèi)，流線型箱梁斷面的顫振臨界風(fēng)速隨攻角絕對(duì)值的增加先增加后降低. 薄平板在+8°攻角時(shí)的顫振性能由“硬顫振”變?yōu)椤败涱澱瘛?，而流線型箱梁斷面在不同攻角下的顫振性能均表現(xiàn)為“軟顫振”.

3）當(dāng)攻角較大時(shí)，薄平板斷面和流線型箱梁斷面均表現(xiàn)出“鈍體”特征，來流空氣繞過斷面前緣時(shí)發(fā)生分離，形成周期性變化的大尺度旋渦，造成主梁結(jié)構(gòu)表面壓力分布不均，產(chǎn)生較大的扭矩，進(jìn)一步激發(fā)結(jié)構(gòu)的氣動(dòng)不穩(wěn)定狀態(tài)，從而導(dǎo)致主梁斷面顫振臨界風(fēng)速降低.

參考文獻(xiàn)

[1]? ? YAYLA S，CANPOLAT C，SAHIN B，et al. The effect of angle of attack on the flow structure over the nonslender lambda wing[J].Aerospace Science and Technology，2013，28（1）：417—430.

[2]? ? 汪清，蔡金獅. 飛機(jī)大攻角氣動(dòng)力建模研究進(jìn)展[J]. 氣動(dòng)實(shí)驗(yàn)與測(cè)量控制，1994，8（3）：7—17.

WANG Q，CAI J S. Advances in aerodynamic modeling of airplane at high angles of attack[J]. Aerodynamic Experiment and Measurement Control，1994，8（3）：7—17. （In Chinese）

[3]? ? REINHOLD T A，BRINCH M，DAMSGAARD A. Wind tunnel tests for the Great Belt Link[C]// Proceedings of the First International Symposium on Aerodynamics of Large Bridges. Rotterdam：Balkema Publishers，1992，2：255—267.

[4]? ? LAROSE G L，LIVESEY F M. Performance of streamlined bridge decks in relation to the aerodynamics of a flat plate[J]. Journal of Wind Engineering and Industrial Aerodynamics，1997，69/70/71：851—860.

[5]? ? 許福友，陳艾榮. 蘇通大橋三維顫振分析[J]. 工程力學(xué)，2008，25（8）：139—144.

XU F Y，CHEN A R. 3-D flutter analysis of Sutong Bridge[J]. Engineering Mechanics，2008，25（8）：139—144. （In Chinese）

[6]? ? 朱青，朱樂東，郭震山. 節(jié)段模型靜風(fēng)附加攻角對(duì)顫振性能的影響[C]//第十四屆全國(guó)結(jié)構(gòu)風(fēng)工程學(xué)術(shù)會(huì)議論文集. 北京：中國(guó)土木工程學(xué)會(huì)，2010：719—724.

ZHU Q，ZHU L D，GUO Z S. Effect of additional angle of attack on flutter performance of segmental mode[C]// Papers Collection of the Fourth National Academic Conference on Structural Wind Engineering. Beijing：China Civil Engineering Society，2010：719—724. （In Chinese）

[7]? ? 朱樂東，朱青，郭震山. 風(fēng)致靜力扭角對(duì)橋梁顫振性能影響的節(jié)段模型試驗(yàn)研究[J]. 振動(dòng)與沖擊，2011，30（5）：23—26.

ZHU L D，ZHU Q，GUO Z S. Effect of wind-induced static torsional angle on flutter performance of bridges via sectional model test[J]. Journal of Vibration and Shock，2011，30（5）：23—26. （In Chinese）

[8]? ? 張宏杰，朱樂東. 附加風(fēng)攻角對(duì)1400 m斜拉橋顫振分析結(jié)果的影響[J]. 振動(dòng)與沖擊，2013，32（17）：95—99.

ZHANG H J，ZHU L D. Influence of additional attack angle on flutter performance of a 1 400 m long cable stayed bridge [J]. Journal of Vibration and Shock，2013，32（17）：95—99. （In Chinese）

[9]? ? 歐陽克儉，陳政清. 中央穩(wěn)定板提高顫振穩(wěn)定性能的細(xì)觀作用機(jī)理[J]. 振動(dòng)與沖擊，2016，35（1）：11—16.

OUYANG K J，CHEN Z Q. Micro-mechanism of a central stabilizer for improving a bridges flutter stability[J]. Journal of Vibration and Shock，2016，35（1）：11—16. （In Chinese）

[10]? 劉志文，謝普仁，陳政清，等. 大跨度流線型箱梁懸索橋顫振穩(wěn)定性氣動(dòng)優(yōu)化[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2019，46（3）：1—9.

LIU Z W，XIE P R，CHEN Z Q，et al. Aerodynamic optimization of flutter stability for a long-span streamlined box girder suspension bridge[J]. Journal of Hunan University（Natural Sciences），2019，46（3）：1—9. （In Chinese）

[11]? DAITO Y，MATSUMOTO M，ARAKI K. Torsional flutter mechanism of two-edge girders for long-span cable-stayed bridge[J]. Journal of Wind Engineering and Industrial Aerodynamics，2002，90（12/13/14/15）：2127—2141.

[12]? 朱樂東，高廣中. 典型橋梁斷面軟顫振現(xiàn)象及影響因素[J]. 同濟(jì)大學(xué)學(xué)報(bào)（自然科學(xué)版），2015，43（9）：1289—1294.

ZHU L D，GAO G Z. Influential factors of soft flutter phenomenon for typical bridge deck sections[J]. Journal of Tongji University（Natural Science），2015，43（9）：1289—1294. （In Chinese）

[13]? 朱青，陳文天，朱樂東，等. 大攻角下超大跨度斜拉橋顫振性能節(jié)段模型風(fēng)洞試驗(yàn)[J]. 中國(guó)公路學(xué)報(bào)，2019，32（10）：67—74.

ZHU Q，CHEN W T，ZHU L D，et al. Flutter performance of a super-long-span cable-stayed bridge under large attack angles via wind tunnel sectional model tests[J]. China Journal of Highway and Transport，2019，32（10）：67—74. （In Chinese）

[14]? 吳長(zhǎng)青，張志田，張偉峰. 考慮幾何非線性的橋梁后顫振極限環(huán)特性[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2018，45（5）：1—10.

WU C Q，ZHANG Z T，ZHANG W F. Limit cycle oscillation of bridge post-flutter with geometric nonlinearities[J]. Journal of Hunan University（Natural Sciences），2018，45（5）：1—10. （In Chinese）

[15]? 竇然.大跨橋梁顫振與流動(dòng)控制數(shù)值模擬[D].哈爾濱：哈爾濱工業(yè)大學(xué)，2015：45—55.

DOU R. Numerical simulation of flutter and flow control for long-span bridges[D]. Harbin：Harbin Institute of Technology，2015：45—55. （In Chinese）

[16]? 劉慕廣，陳政清. 典型鈍體斷面大攻角下的顫振自激力特性[J]. 振動(dòng)與沖擊，2013，32（10）：22—25.

LIU M G，CHEN Z Q. Characteristics of self-excited forces in flutter of typical blunt body under large attack angles[J]. Journal of Vibration and Shock，2013，32（10）：22—25. （In Chinese）

[17]? TANG H J，LI Y L. Flutter performance of central-soltted plate under large wind attack angles[J]. Journal of Wind Engineering and Industrial Aerodynamics，2014，186：192—203.

[18]? TANG H J，LI Y L，WANG Y F，et al. Aerodynamic optimization for flutter performance of steel truss stiffening girder at large angles of attack[J]. Journal of Wind Engineering and Industrial Aerodynamics，2017，168：260—270.

[19]? 王騎，李郁林，李志國(guó)，等. 不同風(fēng)攻角下薄平板的顫振導(dǎo)數(shù)[J]. 工程力學(xué)，2018，35（10）：10—16.

WANG Q，LI Y L，LI Z G，et al. Flutter derivatives of a thin plate model under different attack angles[J]. Engineering Mechanics，2018，35（10）：10—16. （In Chinese）

[20]? 楊陽，張亮亮，吳波，等. 大攻角及橋面粗糙度對(duì)扁平鋼箱梁渦振性能的影響[J]. 土木建筑與環(huán)境工程，2015，37（6）：32—38.

YANG Y，ZHANG L L，WU B，et al. Effects of large attack angle and deck roughness on flat steel box girder vortex-induced vibration[J]. Journal of Civil，Architectural & Environmental Engineering，2015，37（6）：32—38. （In Chinese）

[21]? 劉小兵，張海東，劉慶寬. 大攻角下分離雙鋼箱梁間距對(duì)渦振特性的影響[J]. 振動(dòng)與沖擊，2017，36（14）：202—207.

LIU X B，ZHANG H D，LIU Q K. Influence of the spacing on vortex-induced vibration characteristics of twin separate steel box girders under a large wind attacking angle[J]. Journal of Vibration and Shock，2017，36（14）：202—207. （In Chinese）

[22]? 陳艾榮，艾輝林. 計(jì)算橋梁空氣動(dòng)力學(xué)：大渦模擬[M]. 北京：人民交通出版社，2010：162—173.

CHEN A R，AI H L. Computational bridge aerodynamics：large eddy simulation[M]. Beijing：China Communications Press，2010：162—173. （In Chinese）

[23]? 詹昊，廖海黎. 橋梁主梁端部翼板顫振主動(dòng)控制流固耦合計(jì)算[J]. 振動(dòng)工程學(xué)報(bào)，2018，31（2）：276—282.

ZHAN H，LIAO H L. Flutter active control studies of bridge with flaps attached to deck edge by FSI calculation[J]. Journal of Vibration Engineering，2018，31（2）：276—282. （In Chinese）

[24]? 龔慧星，劉志文. 矩形斷面柱體渦激振動(dòng)數(shù)值模擬與機(jī)理分析[J]. 公路交通科技，2016，33（8）：76—85.

GONG H X，LIU Z W. Numerical simulation and mechanism analysis of vortex-induced vibration of rectangular cylinder[J]. Journal of Highway and Transportation Research and Development，2016，33（8）：76—85. （In Chinese）

[25]? 祝志文，陳政清. 數(shù)值模擬橋梁斷面氣動(dòng)導(dǎo)數(shù)和顫振臨界風(fēng)速[J]. 中國(guó)公路學(xué)報(bào)，2004，17（3）：41—50.

ZHU Z W，CHEN Z Q. Numerical simulations for aerodynamic derivatives and critical flutter of bridge deck[J]. China Journal of Highway and Transport，2004，17（3）：41—50. （In Chinese）

[26]? WALTHER J H. Discrete vortex method for two-dimensional flow past bodies of arbitrary shape under-going prescribed rotary and translational motion[D]. Lyngby：Technical University of Denmark，1994：131—166.

[27]? POULSEN N K，DAMSGAARD A，REINHOLD T A.Determination of flutter derivatives for the Great Belt Bridge[J]. Journal of Wind Engineering and Industrial Aerodynamics，1992，41（1/2/3）：153—164.