橋梁主梁端部翼板顫振主動(dòng)控制流固耦合計(jì)算

詹昊 廖海黎

摘要： 安裝主動(dòng)控制翼板是提高大跨度橋梁顫振穩(wěn)定性的一種有效方法。運(yùn)用流固耦合技術(shù)對(duì)橋梁顫振主動(dòng)控制進(jìn)行計(jì)算分析，可以考慮氣動(dòng)翼板的大幅度扭轉(zhuǎn)以及氣動(dòng)翼板和主梁端部的氣動(dòng)干擾效應(yīng)。通過(guò)對(duì)商用軟件FLUENT二次開(kāi)發(fā)，建立了豎彎和扭轉(zhuǎn)流固耦合數(shù)值仿真計(jì)算模型，并對(duì)主梁端部安裝了主動(dòng)控制翼板的大海帶橋的顫振穩(wěn)定性進(jìn)行了數(shù)值仿真計(jì)算分析。系統(tǒng)地研究了前后翼板相對(duì)于主梁的角速度對(duì)顫振性能的影響。數(shù)值仿真計(jì)算結(jié)果表明：沒(méi)有采用主動(dòng)控制翼板時(shí)，顫振臨界風(fēng)速計(jì)算值和風(fēng)洞實(shí)驗(yàn)值吻合良好。采用主動(dòng)控制翼板后，當(dāng)前翼板角速度與主梁反向，后翼板角速度與主梁同向時(shí)控制效果良好。且隨著氣動(dòng)翼板角速度增大，主梁扭轉(zhuǎn)位移減小。旋渦脫落圖表明：作用在翼板的流場(chǎng)和作用在主梁的流場(chǎng)相互干擾，因此作用在整個(gè)系統(tǒng)上的力矩變化不僅來(lái)源于氣動(dòng)翼板的力矩，而且來(lái)源于流場(chǎng)形態(tài)的改變。計(jì)算表明在上述良好的控制時(shí)，氣動(dòng)翼板提供反向力矩，且與作用在主梁上的力矩相位相反，最大限度地平衡了作用在主梁上的力矩，使作用在主梁系統(tǒng)上的力矩均值減小，這是主動(dòng)控制翼板提高橋梁顫振穩(wěn)定性的原因之一。最后研究了氣動(dòng)翼板合適長(zhǎng)度，計(jì)算表明：當(dāng)氣動(dòng)翼板長(zhǎng)度為主梁寬度的10%～15%時(shí)，顫振主動(dòng)控制效果較好。

關(guān)鍵詞： 橋梁工程； 流固耦合； 顫振臨界風(fēng)速； 顫振主動(dòng)控制； 主梁端部翼板

中圖分類(lèi)號(hào)： U441+.3文獻(xiàn)標(biāo)志碼： A文章編號(hào)： 1004-4523（2018）02-0276-07

DOI：10.16385/j.cnki.issn.1004-4523.2018.02.010

引言

顫振通常是一種具有發(fā)散性質(zhì)的自激振動(dòng)，會(huì)對(duì)橋梁結(jié)構(gòu)造成災(zāi)難性破壞。安裝主動(dòng)控制翼板是提高橋梁顫振穩(wěn)定性的有效方法。自1992年Ostenfeld和Larsen提出利用主動(dòng)控制翼板抑制橋梁顫振的概念以來(lái)[1]，關(guān)于主動(dòng)控制翼板的實(shí)驗(yàn)與理論研究正逐漸引起橋梁工程界的關(guān)注。1992年Kobayashi 和Nagaoka 發(fā)表了第一篇關(guān)于橋梁主動(dòng)控制的文章，通過(guò)節(jié)段模型風(fēng)洞試驗(yàn)表明主動(dòng)控制翼板能將顫振臨界風(fēng)速提高一倍[2]。1999年Wilde 和 Fujino通過(guò)理論分析和風(fēng)洞試驗(yàn)得出運(yùn)用主動(dòng)控制翼板能夠使橋梁顫振臨界風(fēng)速提高57%，并發(fā)現(xiàn)當(dāng)翼板小幅度運(yùn)動(dòng)時(shí)理論分析和風(fēng)洞試驗(yàn)吻合，當(dāng)翼板大幅度運(yùn)動(dòng)時(shí)理論分析結(jié)果遠(yuǎn)小于風(fēng)洞試驗(yàn)結(jié)果[3]。Del Aroc等基于控制面和加勁梁運(yùn)動(dòng)不會(huì)干擾對(duì)方流場(chǎng)的假定，通過(guò)理論分析發(fā)現(xiàn)當(dāng)迎風(fēng)側(cè)控制控制面反向于加勁梁扭轉(zhuǎn)運(yùn)動(dòng)，而背風(fēng)側(cè)控制面同向于加勁梁扭轉(zhuǎn)運(yùn)動(dòng)時(shí)，控制效果最佳，可使懸索橋的顫振臨界風(fēng)速提高2倍以上[4]。論文通過(guò)風(fēng)洞試驗(yàn)研究表明，在如文[4]最佳控制效果時(shí)，隨著翼板相對(duì)于主梁的扭轉(zhuǎn)速度增加，主梁的最大扭轉(zhuǎn)角度有減小的趨勢(shì)[5]。P Omenzetter 提出了一種裝置，利用主纜和加勁梁之間的相對(duì)運(yùn)動(dòng)來(lái)驅(qū)動(dòng)風(fēng)嘴運(yùn)動(dòng)實(shí)現(xiàn)顫振主動(dòng)控制，通過(guò)時(shí)域計(jì)算表明迎風(fēng)風(fēng)嘴相對(duì)主梁反相扭轉(zhuǎn)，背風(fēng)風(fēng)嘴相對(duì)主梁同相扭轉(zhuǎn)提高顫振穩(wěn)定性效果良好[6-7]。2002年劉高基于非定常氣動(dòng)力理論，推算安裝主動(dòng)控制翼板后作用在整個(gè)橋梁的氣動(dòng)力表達(dá)式，從增加系統(tǒng)扭轉(zhuǎn)阻尼的角度，研究了翼板扭轉(zhuǎn)參數(shù)的選取[8]。2013年郭增偉以一座3000 m的懸索橋?yàn)槔治隽丝刂泼媾まD(zhuǎn)運(yùn)動(dòng)相位、運(yùn)動(dòng)方向和運(yùn)動(dòng)振幅對(duì)顫振的影響[9]。以上主動(dòng)控制結(jié)果表明：合理選取主動(dòng)控制翼板扭轉(zhuǎn)運(yùn)動(dòng)參數(shù)，能夠有效地提高橋梁的顫振穩(wěn)定性。

上述文獻(xiàn)中所采用的計(jì)算方法為頻域或時(shí)域顫振計(jì)算方法，其中氣動(dòng)力由基于顫振導(dǎo)數(shù)的Scanlan公式表示，以上方法存在局限性：1）只能假設(shè)氣動(dòng)翼板作小位移運(yùn)動(dòng)，不能考慮流動(dòng)分離現(xiàn)象。2）難以考慮主梁和氣動(dòng)翼板之間的氣動(dòng)力干擾效應(yīng)。不同以往的研究方法。本文采用了新的方法——流固耦合技術(shù)對(duì)此問(wèn)題進(jìn)行研究。對(duì)商用軟件FLUENT進(jìn)行二次開(kāi)發(fā)，建立主梁豎彎和扭轉(zhuǎn)流固耦合數(shù)值仿真計(jì)算模型，設(shè)定氣動(dòng)翼板相對(duì)于主梁運(yùn)動(dòng)的運(yùn)動(dòng)規(guī)律，研究主動(dòng)控制翼板對(duì)大跨度橋梁顫振穩(wěn)定性的影響。該方法允許氣動(dòng)翼板的大幅度運(yùn)動(dòng)，可以考慮主梁和翼板之間的氣動(dòng)干擾效應(yīng)。并通過(guò)結(jié)構(gòu)受力分析初步揭示了顫振主動(dòng)控制機(jī)理。

1數(shù)值仿真計(jì)算原理

將主梁作為質(zhì)量、彈簧和阻尼系統(tǒng)，這里沒(méi)有考慮氣動(dòng)翼板的質(zhì)量。數(shù)值仿真計(jì)算原理示意圖如圖1所示。

圖1數(shù)值仿真計(jì)算原理示意圖

Fig.1Schematic diagram of numerical simulation calculation

主梁豎彎運(yùn)動(dòng)和扭轉(zhuǎn)運(yùn)動(dòng)方程如下所示：m+ch+khy=Ltotal（1）

Iθ+cθ+kθθ=Mtotal（2）式中m和Iθ分別表示主梁的質(zhì)量和扭轉(zhuǎn)慣性矩，kh和kθ分別表示主梁的豎向剛度和扭轉(zhuǎn)剛度，ch和cθ表示主梁的豎彎和扭轉(zhuǎn)阻尼，y和θ分別表示主梁的豎向位移和扭轉(zhuǎn)位移。Ltotal和Mtotal表示作用在主梁系統(tǒng)上的總氣動(dòng)力。

圖2主梁系統(tǒng)受力示意圖

Fig.2Schematic diagram of force acting on the main beam system

主梁系統(tǒng)受力如圖2所示，其中Ltotal=Ldeck+Lle+Ltr，Mtotal=Mdeck+Mle+Mtr，Mdeck，Mle和Mtr 分別表示主梁、前翼板和后翼板對(duì)于系統(tǒng)彈性中心的氣動(dòng)扭矩。氣動(dòng)翼板以靠近主梁的一端為軸進(jìn)行轉(zhuǎn)動(dòng)，氣動(dòng)翼板豎向位移與主梁保持一致，氣動(dòng)翼板的角速度根據(jù)主梁運(yùn)動(dòng)的角速度進(jìn)行調(diào)整。ωle=Aleωdeck，ωtr=Atrωdeck，其中ωdeck為主梁的角速度，ωle和ωtr分別為前翼板和后翼板的角速度。Ale和Atr是氣動(dòng)翼板相對(duì)于主梁角速度的系數(shù)。

Ale=Atr=1表示ωle=ωtr=ωdeck，即氣動(dòng)翼板的角速度與主梁的角速度相同，相當(dāng)于固定翼板。對(duì)于不可壓縮流體的連續(xù)方程和納維-斯脫克斯方程如下所示：·V=0（3）

Vt+（V·）V=-1ρp+ν2V（4）式中ρ為流體密度；V，p，t分別表示速度向量，壓強(qiáng)和時(shí)間，ν為運(yùn)動(dòng)黏度。通過(guò)FLUENT求解方程（3），（4），求出作用在物體上的力。并提取升力和力矩帶入主梁振動(dòng)方程（1），（2），將Newmark方法的C語(yǔ)言代碼同F(xiàn)LUENT軟件連接，通過(guò)Newmark方法求解主梁的動(dòng)力響應(yīng)。然后運(yùn)用FLUENT的動(dòng)網(wǎng)格技術(shù)，將主梁的速度傳遞給網(wǎng)格得到網(wǎng)格位置的更新，開(kāi)始下一個(gè)時(shí)間步的計(jì)算，如此循環(huán)得到各時(shí)間步振動(dòng)位移，最后設(shè)定氣動(dòng)翼板相對(duì)于主梁運(yùn)動(dòng)的數(shù)學(xué)函數(shù)關(guān)系式，實(shí)現(xiàn)顫振主動(dòng)控制。

2數(shù)值仿真計(jì)算模型

數(shù)值仿真計(jì)算模型按照實(shí)際橋梁尺寸建模，橋梁計(jì)算參數(shù)如表1所示。主梁截面如圖3和4所示，計(jì)算區(qū)域和計(jì)算網(wǎng)格如圖5所示。

表1計(jì)算參數(shù)

Tab.1Calculation parameters

質(zhì)量/

（kg·m-1）質(zhì)量慣性矩/

（kg·m2·m-1）豎彎頻率/

Hz扭轉(zhuǎn)頻率/

Hz236872.501×1060.0970.270

圖3主梁截面尺寸（單位：m）

Fig.3Section size of main beam（Unit：m）圖4主梁計(jì)算模型

Fig.4Calcultion model of main beam

圖5計(jì)算區(qū)域和計(jì)算網(wǎng)格（設(shè)置氣動(dòng)翼板）

Fig.5Calculation region and calculation mesh（with active wing plate）

由圖5可知，計(jì)算區(qū)域長(zhǎng)、寬分別為170 m和160 m。風(fēng)向從左至右，左側(cè)設(shè)定為速度入口，右側(cè)設(shè)定為自由出流。上下邊界為無(wú)滑移固壁邊界，靠近物體網(wǎng)格加密，遠(yuǎn)離物體網(wǎng)格逐漸稀疏。由于壁面是能量交換的主要地方，主梁和氣動(dòng)翼板設(shè)定了邊界層網(wǎng)格。邊界層網(wǎng)格不變形，隨主梁和翼板一起運(yùn)動(dòng)，保證了主梁和氣動(dòng)翼板壁面周?chē)W(wǎng)格的質(zhì)量。網(wǎng)格總數(shù)4萬(wàn)左右，第一層網(wǎng)格距離主梁壁面0.04 m，計(jì)算時(shí)間步長(zhǎng)取0.002 s，湍流強(qiáng)度取6.5%，采用有限體積法求解，其中對(duì)流項(xiàng)采用二階中心差分格式，壓力和速度的耦合采用SMPLEC算法，計(jì)算采用LES湍流模型。沒(méi)有設(shè)置氣動(dòng)翼板的計(jì)算區(qū)域和計(jì)算網(wǎng)格設(shè)置同圖5。

3數(shù)值仿真計(jì)算結(jié)果〖*2〗3.1沒(méi)有安裝氣動(dòng)翼板豎彎阻尼比和扭轉(zhuǎn)阻尼比取0.002，0°風(fēng)攻角時(shí)豎向位移時(shí)程和扭轉(zhuǎn)位移時(shí)程曲線如圖6和7所示。

圖6豎向位移時(shí)程曲線（未安裝主動(dòng)控制翼板）

Fig.6Time histories of vertical displacements（without active wing plate）

圖7扭轉(zhuǎn)位移時(shí)程曲線（未安裝主動(dòng)控制翼板）

Fig.7Time histories of torsional displacements（without active wing plate）如圖6和7所示，當(dāng)風(fēng)速V=72 m/s時(shí)，主梁豎向位移幅值隨時(shí)間基本保持不變，主梁扭轉(zhuǎn)位移幅值隨時(shí)間呈現(xiàn)衰減的趨勢(shì)；當(dāng)風(fēng)速V=74 m/s時(shí)，主梁豎向位移和扭轉(zhuǎn)位移幅值隨時(shí)間逐漸增大，發(fā)生了彎扭耦合顫振，此時(shí)豎向振動(dòng)頻率和扭轉(zhuǎn)振動(dòng)頻率大約是0.2 Hz，介于豎向自振頻率和扭轉(zhuǎn)自振頻率之間。顫振臨界風(fēng)速如表2所示。

表2顫振臨界風(fēng)速[10]

Tab.2Critical flutter wind speed[10]

風(fēng)攻

角/（°）顫振臨界風(fēng)速/（m·s-1）節(jié)段模型風(fēng)洞實(shí)驗(yàn)數(shù)值計(jì)算070.0～72.972.0～74.0由表2可知，數(shù)值仿真計(jì)算值和節(jié)段模型風(fēng)洞試驗(yàn)值吻合良好。

3.2安裝氣動(dòng)翼板

豎彎阻尼比和扭轉(zhuǎn)阻尼比取0.005，翼板運(yùn)動(dòng)方式主梁最大扭轉(zhuǎn)位移的影響如表3所示。

由表3可知，以下工況顫振控制效果不好：1）當(dāng)前翼板和后翼板都與主梁的扭轉(zhuǎn)方向相同；2）前翼板和后翼板都與主梁的扭轉(zhuǎn)方向相反；3）前翼板扭轉(zhuǎn)方向與主梁的扭轉(zhuǎn)方向相同，后翼板扭轉(zhuǎn)方向與主梁的扭轉(zhuǎn)方向相反。只有當(dāng)前翼板扭轉(zhuǎn)方向與主梁的扭轉(zhuǎn)方向相反，后翼板扭轉(zhuǎn)方向與主梁的扭轉(zhuǎn)方向相同時(shí)，顫振控制效果良好。在這種情況下，隨著翼板扭轉(zhuǎn)速度的增加，主梁扭轉(zhuǎn)角度減小。這種變化規(guī)律和文獻(xiàn)[5]中風(fēng)洞實(shí)驗(yàn)的結(jié)論一致。這時(shí)橋梁顫振臨界風(fēng)速達(dá)到120 m/s，提高到1.6倍以上。旋渦脫落如圖8所示。

由圖8可見(jiàn)，氣動(dòng)翼板發(fā)生了流動(dòng)分離現(xiàn)象。由于氣動(dòng)翼板和主梁的運(yùn)動(dòng)，作用在翼板的流場(chǎng)和作用在主梁的流場(chǎng)相互干擾。作用在主梁系統(tǒng)上的力矩變化不僅來(lái)源于氣動(dòng)翼板的力矩，而且來(lái)源于流場(chǎng)形態(tài)的改變。當(dāng)V=90 m/s，ωle=-12ωdeck，ωtr=12ωdeck時(shí)，扭轉(zhuǎn)位移時(shí)程曲線和扭轉(zhuǎn)速度時(shí)程曲線如圖9和10所示。

由圖9和10可知，主梁扭轉(zhuǎn)位移很小，前翼板與主梁扭轉(zhuǎn)運(yùn)動(dòng)方向相反，后翼板與主梁扭轉(zhuǎn)運(yùn)動(dòng)

表3翼板角速度對(duì)主梁最大扭轉(zhuǎn)位移的影響

Tab.3The influence of wing plate rotational velocity on maximum torsional displacement amplitude of the main beam

（a） V=75 m/s翼板扭轉(zhuǎn)角速度/[（°）·s-1]主梁最大扭

轉(zhuǎn)角度/（°）ωle=ωdeck，ωtr=ωdeck發(fā)散ωle=6ωdeck，ωtr=6ωdeck發(fā)散ωle=6ωdeck，ωtr=-6ωdeck發(fā)散ωle=-6ωdeck，ωtr=-6ωdeck發(fā)散ωle=-6ωdeck，ωtr=6ωdeck-2.2（b） V=90 m/s翼板扭轉(zhuǎn)角速度/[（°）·s-1]主梁最大扭

轉(zhuǎn)角度/（°）ωle=-4ωdeck，ωtr=4ωdeck發(fā)散ωle=-6ωdeck，ωtr=6ωdeck-3.528ωle=-12ωdeck，ωtr=12ωdeck-2.002ωle=-18ωdeck，ωtr=18ωdeck-1.177（c） V=120 m/s翼板扭轉(zhuǎn)角速度/[（°）·s-1]主梁最大扭

轉(zhuǎn)角度/（°）ωle=-6ωdeck，ωtr=6ωdeck發(fā)散ωle=-12ωdeck，ωtr=12ωdeck-2.656ωle=-18ωdeck，ωtr=18ωdeck-1.396圖8旋渦脫落圖

Fig.8Vortex shedding diagram

圖9扭轉(zhuǎn)位移時(shí)程曲線（安裝主動(dòng)控制翼板，V=90 m/s）

Fig.9Time histories of torsional displacements（with active wing plate，V=90 m/s ）

方向相同，前后翼板扭轉(zhuǎn)速度和位移是主梁扭轉(zhuǎn)速度和位移的12倍，運(yùn)動(dòng)規(guī)律與所設(shè)定的運(yùn)動(dòng)規(guī)律吻合。當(dāng)V=120 m/s， ωle=-12ωdeck，ωtr=12ωdeck，時(shí)，對(duì)系統(tǒng)彈性中心取矩，作用在單位長(zhǎng)度上的力矩位移時(shí)程曲線如圖11所示。

由圖11知，不考慮氣動(dòng)翼板的作用時(shí)，作用于主梁的氣動(dòng)力矩Mdeck為負(fù)值，均值約為-600000 N·m；前翼板和后翼板對(duì)主梁提供的氣動(dòng)力矩均值為正值，且與作用在主梁的氣動(dòng)力矩Mdeck相位相反。考慮氣動(dòng)翼板的作用后，作用在主梁系統(tǒng)的力矩Mtotal均值約為-300000 N·m，減小了一倍。同理分析得到當(dāng)V=120 m/s， ωle=-18ωdeck，ωtr=18ωdeck ，前翼板和后翼板提供了更大的反向力矩，作用于整個(gè)主梁系統(tǒng)的力矩均值減小了，約為-200000 N·m。

Wilde等對(duì)設(shè)置不同氣動(dòng)翼板寬度（1.5，3.0，4.5，6 m）主梁的橋梁進(jìn)行了顫振主動(dòng)控制研究，Wilde認(rèn)為氣動(dòng)翼寬度取為橋?qū)挼?.1倍時(shí)，最為經(jīng)濟(jì)合理[11]。本文運(yùn)用數(shù)值計(jì)算方法對(duì)此進(jìn)行了研究，當(dāng)V=120 m/s， ωle=-12ωdeck，ωtr=12ωdeck ，時(shí)，主梁最大扭轉(zhuǎn)角度如表4所示。

圖10扭轉(zhuǎn)速度時(shí)程曲線（安裝主動(dòng)控制翼板，V=90 m/s）

Fig.10Time histories of torsional velocity（with active wing plate，V=90 m/s ）

由表4可知，氣動(dòng)翼板長(zhǎng)度為2 m時(shí)，顫振主動(dòng)控制效果差。隨著氣動(dòng)翼板長(zhǎng)度增加，主梁最大扭轉(zhuǎn)角度減小。隨著氣動(dòng)翼板長(zhǎng)度繼續(xù)增加，主梁最大扭轉(zhuǎn)角度又開(kāi)始增加。氣動(dòng)翼板長(zhǎng)度為主梁寬度的10%～15%，顫振主動(dòng)控制效果較好，這與Wilde的研究結(jié)論基本相同。

表4不同翼板長(zhǎng)度時(shí)主梁最大扭轉(zhuǎn)角度（V=120 m/s）

Tab.4The maximum torsion angle of main beam under different length of wing plate （V=120 m/s）

翼板長(zhǎng)度/m扭轉(zhuǎn)角度/（°）2發(fā)散振動(dòng)3-2.8344-2.0025-2.1156-2.3334結(jié)論

本文通過(guò)對(duì)FLUENT軟件進(jìn)行二次開(kāi)發(fā)，運(yùn)

圖11作用在單位長(zhǎng)度上的力矩 （V=120 m/s）

Fig.11Moment force acting on the structure of per unit length （V=120 m/s）

用流固耦合方法研究了主梁端部主動(dòng)控制翼板對(duì)橋梁顫振穩(wěn)定性的影響，主要結(jié)論如下：

（1）大海帶橋顫振臨界風(fēng)速數(shù)值仿真計(jì)算值與風(fēng)洞實(shí)驗(yàn)結(jié)果吻合良好，顫振形態(tài)為彎扭耦合顫振，符合流線型箱梁的顫振形態(tài)，顫振頻率介于豎彎頻率和扭轉(zhuǎn)頻率之間。

（2） 在主梁風(fēng)嘴前后端安裝主動(dòng)控制翼板，當(dāng)前翼板扭轉(zhuǎn)運(yùn)動(dòng)方向與主梁的扭轉(zhuǎn)方向相反，后翼板扭轉(zhuǎn)運(yùn)動(dòng)方向與主梁扭轉(zhuǎn)方向相同時(shí)，顫振控制效果良好。在此情況下，隨著翼板相對(duì)于主梁的扭轉(zhuǎn)速度增加，主梁的最大扭轉(zhuǎn)角度有減小的趨勢(shì)。 比沒(méi)有安裝氣動(dòng)翼板，顫振臨界風(fēng)速提高到1.6倍以上。

（3） 由旋渦脫落圖可知作用在翼板的流場(chǎng)和作用在主梁的流場(chǎng)相互干擾，因此作用在整個(gè)系統(tǒng)上的力矩變化不僅來(lái)源于氣動(dòng)翼板的力矩，而且來(lái)源于流場(chǎng)形態(tài)的改變。在如（2）所述的良好控制時(shí)，由于氣動(dòng)翼板提供反向力矩，且與單獨(dú)作用在主梁的氣動(dòng)力矩相位相反，最大限度地平衡了作用在主梁上的力矩，使作用在主梁系統(tǒng)上的力矩均值減小，這可能是主動(dòng)控制翼板提高顫振穩(wěn)定性的原因之一。

（4）氣動(dòng)翼板的長(zhǎng)度過(guò)小或過(guò)大，顫振主動(dòng)控制效果不好。當(dāng)氣動(dòng)翼板長(zhǎng)度為主梁寬度的10%～15%時(shí)，顫振主動(dòng)控制效果較好。

流固耦合計(jì)算方法為橋梁顫振主動(dòng)控制提供了一種新的研究手段，但主動(dòng)控制翼板提高橋梁顫振穩(wěn)定性的機(jī)理還需進(jìn)一步深入研究。

參考文獻(xiàn)：

[1]Ostenfeid K，Larsen A.Bridge engineering and aerodynamics[C] Proceedings of the First International Symposium on Aerodynamics of Large Bridges. A.A.Balkema， Rotterdam，1992：3—22.

[2]Kobayashi K，Nagaoka H.Active control of flutter of a suspension bridge[J]. Journal of Wind Engineering and Industrial Aerodynamics， 1992，41&44：143—151.

[3]Wilde K，F(xiàn)ujino Y，Kawakami T. Analytical and experimental study on passive aerodynamic control of flutter of a bridge deck[J]. Journal of Wind Engineering and Industrial Aerodynamics， 1999，80（1/2）：105—119.

[4]DEL ARCO D C，APARICIO A C.Improving suspension bridge wind stability with aerodynamic appendages[J].Journal of Structural Engineering，1999， 125（12）： 1367—1375.

[5]具有主動(dòng)振動(dòng)控制橋梁模型的風(fēng)洞試驗(yàn)[J].楊永賢，譯.公路交通技術(shù)，2002，（4）：78—81.

Wind-tunnel test of a bridge model with active vibration control[J].Translated by YANG Yong-xian.Technology of Highway and Transport，2002，（4）：78—81.

[6]Omenzetter P， Wilde K， Fujino Y， Suppression of wind-induced instabilities of a long span bridge by a passive deck-flaps control system Part I： Numerical simulations[J]. Journal of Wind Engineering and Industrial Aerodynamics， 2000：61—79.

[7]Omenzetter P，Wilde K， Fujino Y. Suppression of wind-induced instabilities of a long span bridge by a passive deck-flaps control system Part II： Numerical simulations[J]. Journal of Wind Engineering and Industrial Aerodynamics， 2000：81—91.

[8]劉高，強(qiáng)士中，王秀偉.主動(dòng)控制翼板抑制懸索橋顫振研究[J].應(yīng)用力學(xué)學(xué)報(bào)，2002，19（4）：52—57.

LIU Gao， QIANG Shi-zhong ，WANG Xiu-wei，Study of active control wing plate for flutter suppression of long span suspension bridge[J]. Chinese Journal of Applied Mechanics， 2002，19（4）：52—57.

[9]郭增偉，楊詠昕，葛耀君.大跨懸索橋顫振主動(dòng)控制面理論研究[J].中國(guó)公路學(xué)報(bào)，2013，26（2）：119—126.

GUO Zeng-wei，YANG Yong-xin，GE Yao-jun.Theoretical study of active control surface for flutter suppression of long span suspension bridge[J].China Journal of Highway and Transport，2013，26（2）： 119—126.

[10]Niels J Gimsing.大貝耳特海峽：東橋[M].西南交通大學(xué)土木工程學(xué)院橋梁工程系，中鐵大橋局武漢橋梁科學(xué)研究院譯.成都：西南交通大學(xué)出版社，2008：136—137.

Niels J Gimsing. Great Belt Strait： East Bridge[M]. Translated by School of Civil Engineering， Southwest Jiaotong University， Bridge Science Research Institute Ltd.，China Railway Major Bridge Engineering Group. Chengdu：Southwest Jiaotong University Press，2008：136—137.

[11]Wilde K，Omenze P，F(xiàn)ujino Y.Suppression of bridge flutter by active deck-flaps control system[J].Journal of Engineering Mechanics，2001，127（1）：80—90.

Flutter active control studies of bridge with flaps attached to

deck edge by FSI calculation

ZHAN Hao1，2， LIAO Hai-li1

（1.School of Civil Engineering， Southwest Jiaotong University， Chengdu 610031， China；

2.China Railway Major Bridge Reconnaissance & Design Institute Co.， Ltd.， Wuhan 430056， China）

Abstract： Long span suspension bridges， due to their flexibility and lightness， are much prone to flutter instability. A relatively new research area on aerodynamic stability for long span bridges is based on actively controlled flaps attached along the girders. Different from previous research， this paper proposes a new method to study the question. By secondary development of commercial computational fluid dynamics software FLUENT， the paper establishes a two-dimensional bending and torsional fluid-structure interaction numerical model to study flutter stability of a long span suspension bridges with active control wing plate which is attached directly to the edges of the deck. This method not only allows the large amplitude motion of the wing plate， but also considers the interaction between the aerodynamic force of the main beam and wing plate. The flutter stability is studied by controlling the rotational velocity of the wing plates with respect to that of the bridge deck. Numerical results show that the flutter critical wind speed of girder without wing plate is in good agreement with wind tunnel test. It is an optimal control law that the leading surface rotates in the opposite direction and the trailing surface rotates in the same direction with respect to the deck motion. The maximum torsional displacement amplitude decreases when increasing the rotational velocity of the winglet with respect to that of the deck. The rule is consistent with the conclusion of the wind tunnel experiment of relevant literature. The vortex shedding flow shows that in this system the flow patterns around the deck and flap are interacted by each other. So the stabilizing moment not only comes from the aerodynamic forces generated on control wing plates but can also be achieved through modification of the aerodynamic forces induced on the bridge deck. The calculated results show that in the optimal control， the aerodynamic wing provides the reverse moment and the moment phase is opposite to that of the main beam. It will maximize the moment balance of the main beam and decreases the average moment on the main beam system. This may be one of the reasons to improve the flutter stability of bridge with active wing plate. Finally， the proper length of the aerodynamic wing plate is studied. The calculation results show that the flutter active control is better when the wing length is 10%-15% of the main beam width.

Key words： bridge engineering；fluid-structure interaction； critical flutter wind speed； flutter active control； wing plate attached to edges of deck