祝志文　顏爽　王欽華　李加武

摘要： 為研究兩自由度薄平板大幅運(yùn)動的氣動特性，評價其氣動穩(wěn)定性，基于任意拉格朗日?歐拉描述法的動網(wǎng)格技術(shù)，通過有限差分法求解描述任意流變區(qū)域不可壓流動的控制方程，開展了不同折減風(fēng)速下薄平板豎彎和扭轉(zhuǎn)運(yùn)動繞流場的CFD （Computational Fluid Dynamics）模擬。研究認(rèn)為，單自由度薄平板小幅運(yùn)動的氣動力系統(tǒng)是線性和穩(wěn)定的，即使單自由度大幅豎彎運(yùn)動也是線性和氣動穩(wěn)定的。但單自由度大幅扭轉(zhuǎn)運(yùn)動的平板氣動力系統(tǒng)出現(xiàn)非線性，并隨折減風(fēng)速的提高非線性變得顯著，且平板將進(jìn)入氣動不穩(wěn)定狀態(tài)。另外，大幅扭轉(zhuǎn)耦合不同豎彎振幅運(yùn)動的平板，氣動力系統(tǒng)均為非線性并隨折減風(fēng)速的提高越加顯著，而該非線性主要來自扭轉(zhuǎn)自由度的大幅運(yùn)動;對該兩自由度耦合系統(tǒng)，當(dāng)豎彎振幅較小和折減風(fēng)速較高時，氣動力系統(tǒng)是不穩(wěn)定的;但當(dāng)豎彎振幅較大時，氣動力系統(tǒng)將是穩(wěn)定的。

關(guān)鍵詞： 氣動穩(wěn)定性; 薄平板; 大幅運(yùn)動; CFD; 氣動力非線性

中圖分類號： U448.21+3; V211.3 ? ?文獻(xiàn)標(biāo)志碼： A ? ?文章編號： 1004-4523（2021）02-0271-12

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

引 ?言

橋梁跨徑的不斷增大使得結(jié)構(gòu)的頻率和阻尼不斷降低，導(dǎo)致橋梁對自然風(fēng)作用的敏感程度明顯增加，因而可能引起橋梁主梁的大幅運(yùn)動。這些大幅運(yùn)動，包括主梁的大幅渦激運(yùn)動[1]和顫振[2]，以及拉索的大幅振動[3?4]。對鈍體外形的橋梁主梁，或主梁因大幅運(yùn)動產(chǎn)生較大的相對攻角效應(yīng)，橋梁主梁的氣動力系統(tǒng)可能會因這種大幅運(yùn)動而呈現(xiàn)顯著的非線性[5?6]。開展橋梁主梁大幅豎彎和扭轉(zhuǎn)運(yùn)動下的氣動特性研究，有助于評價現(xiàn)有線性氣彈理論的應(yīng)用條件和不足，有利于揭示橋梁大幅運(yùn)動氣動力系統(tǒng)的非線性特征和影響因素，以及橋梁氣動穩(wěn)定性的影響因素[7]。

在風(fēng)洞中開展橋梁主梁大幅豎彎和扭轉(zhuǎn)運(yùn)動試驗(yàn)，存在試驗(yàn)驅(qū)動機(jī)構(gòu)設(shè)計和布置、測力和測壓、模型保護(hù)等諸多困難，而基于CFD開展這一研究，不僅能克服上述困難，還能方便地設(shè)置運(yùn)動參數(shù)和邊界條件，具有風(fēng)洞試驗(yàn)無法相比的優(yōu)勢。當(dāng)基于CFD開展橋梁主梁的大幅運(yùn)動模擬時，因主梁周圍空氣所占據(jù)的實(shí)際空間隨時間在發(fā)生顯著的變化，具有大幅運(yùn)動的邊界，對此類顯著流變區(qū)域流體的描述，需采用任意拉格朗日?歐拉描述法[8]。下面以橋梁氣彈研究常采用的薄平板為研究對象，研究其大幅運(yùn)動的氣動特征和氣動穩(wěn)定性。

1 薄平板強(qiáng)迫運(yùn)動繞流場控制方程

基于二維流體運(yùn)動的任意拉格朗日?歐拉描述方程[8]，通過建立計算域和物理域的隨時間變化的映射關(guān)系，可在不隨時間變化的二維計算域上，采用歐拉描述法實(shí)現(xiàn)對二維物理域上任意時變區(qū)域的流體流動的描述。根據(jù)橋梁主梁剛性斷面的事實(shí)，以及主梁豎彎和扭轉(zhuǎn)運(yùn)動的特點(diǎn)，可簡化得到薄平板強(qiáng)迫運(yùn)動繞流場控制方程。

如圖1所示為大幅運(yùn)動的薄平板，設(shè)OXY表示絕對參考系下的直角坐標(biāo)系，并用oxy表示隨平板一起做剛性運(yùn)動的動參考系下的直角坐標(biāo)系，平板具有繞平衡位置的豎向自由度和繞轉(zhuǎn)動軸的扭轉(zhuǎn)自由度。設(shè)平板豎向運(yùn)動速度為V0，繞通過其中點(diǎn)垂直于OXY平面軸的扭轉(zhuǎn)角速度為ω（t）=dψ（t）/dt。對薄平板上述給定的兩自由度運(yùn)動，其對應(yīng)的任意拉格朗日?歐拉描述方程可簡化為動參考系下關(guān)于原參變量u，v和p的N?S方程[9]，也即：

3.2 平板豎彎和扭轉(zhuǎn)運(yùn)動CFD模擬方法驗(yàn)證

為驗(yàn)證薄平板豎彎或扭轉(zhuǎn)運(yùn)動CFD模擬數(shù)值方法的合理性，本文以非定常氣動力模擬為依據(jù)，通過有解析解的薄平板顫振導(dǎo)數(shù)Theodorsen解析解與本文識別結(jié)果的對比，驗(yàn)證本文方法的合理性。對豎彎或扭轉(zhuǎn)運(yùn)動的二自由度運(yùn)動薄平板，其受到的非定常氣動力作用可用八個顫振導(dǎo)數(shù)來表達(dá)[11?12]，即：

式中 ?k=2πfb/U_0為折減頻率;H_i^*和A_i^* （i=1?4）是顫振導(dǎo)數(shù)。

在不同折減風(fēng)速下開展CFD計算，將分別得到作用在豎彎和扭轉(zhuǎn)運(yùn)動薄平板上的氣動力，利用式（18）和（19）并通過最小二乘法可識別平板對應(yīng)折減風(fēng)速下的八個氣動導(dǎo)數(shù)。本文CFD模擬首先開展了網(wǎng)格無關(guān)性檢查，檢查對象是顫振導(dǎo)數(shù)中相對重要的A_2^*導(dǎo)數(shù)在折減風(fēng)速為4，8和12上的值，物面第一個網(wǎng)格高度分別取0.0005b，0.001b和0.002b三種情況，將計算結(jié)果與Theodorsen解析解對比表明，在0.0005b和0.001b兩種情況下得到的結(jié)果與解析解一致性最好，也好于0.002b的情況。從非定常計算較大的計算量上考慮，本文CFD計算平板物面第一個網(wǎng)格高度取0.001b。

如圖4所示，將本文識別的顫振導(dǎo)數(shù)和平板Theodorsen解析解作了比較?？梢娫诘驼蹨p風(fēng)速下，二者吻合很好;隨著折減風(fēng)速的增大，二者出現(xiàn)較小偏差。偏差的原因可能是解析解是基于無黏、無流動分離的勢流理論，但本文是真實(shí)的有黏空氣流動，且平板運(yùn)動會導(dǎo)致流動分離，而黏性效應(yīng)和流動分離可能隨折減風(fēng)速的增大而影響顯著。本文識別結(jié)果與Theodorsen解析解[13]趨勢性非常一致，這說明本文數(shù)值模擬能有效地獲得運(yùn)動薄平板上作用的氣動力，因而能用于研究運(yùn)動薄平板的氣動力和能量特征。

3.3 單自由度小幅運(yùn)動

為便于對比平板大幅運(yùn)動的氣動力和能量特征，下面給出Vr=8下，平板分別作小幅豎彎（振幅為B/80）或小幅扭轉(zhuǎn)運(yùn)動（振幅0.2°）的氣動力時程。圖5（a）和（b）分別是從多個氣動力周期中取出一個完整周期（每個周期氣動力完全相同），對應(yīng)平板單自由度豎彎或單自由度扭轉(zhuǎn)運(yùn)動氣動力系數(shù)時程曲線，可見無論是升力還是扭矩系數(shù)，其曲線具有明顯的諧波特性，其時程頻率等于扭轉(zhuǎn)強(qiáng)迫運(yùn)動頻率。因此當(dāng)平板單自由度小幅運(yùn)動時，其氣動力系統(tǒng)具有良好的線性特性。

圖6為薄平板單自由度小幅扭轉(zhuǎn)運(yùn)動在一個周期內(nèi)4個關(guān)鍵時刻的風(fēng)壓拓?fù)鋱D，對應(yīng)從平衡位置出發(fā)后達(dá)到最大正攻角位置（圖6（a）），然后回到平衡位置（圖6（b）），再達(dá)到最大負(fù)攻角位置（圖6（c）），最后回到零攻角的平衡位置（圖6（d））?？梢姀囊粋€狀態(tài)到另外一個狀態(tài)，壓力呈現(xiàn)規(guī)律性變化，反映了平板扭轉(zhuǎn)運(yùn)動方向的改變。平板前緣為風(fēng)壓值大和梯度高的區(qū)域，后緣為低風(fēng)壓區(qū)且壓力變化遠(yuǎn)不及前緣劇烈。

圖7（a）和（b）是在不同折減風(fēng)速下豎彎運(yùn)動的平板，一個周期內(nèi)氣流分別在平板上表面各點(diǎn)、以及整個平板上所做的無量綱總功?？梢姎饬鲗ζ桨灞砻嫠悬c(diǎn)均做負(fù)功，也即平板表面均消耗氣流能量，前緣（L.E.）局部耗能最強(qiáng)，后緣（T.E.）耗能最弱。圖7（b）給出了Vr=8?60，豎彎運(yùn)動平板一個周期內(nèi)氣流作用總功隨折減風(fēng)速的變化。從能量的角度可見，單自由度小幅豎彎運(yùn)動的平板是氣動穩(wěn)定的。

圖8（a）和（b）分別為單自由度小幅扭轉(zhuǎn)運(yùn)動的平板，一個周期內(nèi)氣流在平板上表面各點(diǎn)、以及整個平板所做的無量綱總功。從圖8（a）可見平板前緣為穩(wěn)定區(qū)，其后到轉(zhuǎn)動軸（A.R）為不穩(wěn)定區(qū)，也即氣流對平板做正功，平板從氣流中吸收能量，因此為激勵區(qū);轉(zhuǎn)動軸到平板后緣氣流做負(fù)功，因而是耗散氣流能量的穩(wěn)定區(qū)，此時平板上激勵區(qū)與穩(wěn)定區(qū)共存。圖8（b）給出了Vr=8?60，一個扭轉(zhuǎn)周期內(nèi)氣流對平板所做的總功，可見在所計算的折減風(fēng)速范圍內(nèi)，一個周期內(nèi)氣流對扭轉(zhuǎn)運(yùn)動的平板均做負(fù)功，也即平板消耗氣流的能量，因此單自由度小幅扭轉(zhuǎn)運(yùn)動的平板也是氣動穩(wěn)定的。

3.4 單自由度大幅運(yùn)動

為研究薄平板單自由度大幅運(yùn)動的氣動特性，需對豎彎或扭轉(zhuǎn)自由度方向設(shè)定遠(yuǎn)大于上述單自由度小幅運(yùn)動的幅值。根據(jù)文獻(xiàn)[7]對橋梁斷面氣動力非線性特性的研究結(jié)果，本文研究中，取大幅豎彎振幅為B/4，大幅扭轉(zhuǎn)振幅為20°，以研究其氣動力的線性特征和氣動穩(wěn)定性。

3.4.1 單自由度大幅豎彎運(yùn)動（LAH）

圖9（a）和（b）分別給出Vr=12和Vr=40時，分別從多個氣動力周期時程中，取出一個完整周期（每個周期氣動力完全相同），對應(yīng)單自由度大幅豎彎運(yùn)動的氣動力時程?？梢婋m然折減風(fēng)速差別大，但升力和扭矩系數(shù)時程的相位差變化不明顯，特別是升力和扭矩系數(shù)時程曲線均具有較好的諧波特性，也即沒有觀察到明顯的氣動力非線性，因此可以認(rèn)為，大幅豎彎運(yùn)動的薄平板氣動力系統(tǒng)的非線性不明顯。

圖10（a）和（b）分別是在一個周期內(nèi)，氣流在大幅豎彎運(yùn)動平板上表面各點(diǎn)、以及對整個平板所做的無量綱總功?？梢姎饬鲗ωQ彎運(yùn)動的平板表面所有點(diǎn)均做負(fù)功，且隨著折減風(fēng)速的提高，功的絕對值雖越來越小，但仍為負(fù)。同時，在全部折減風(fēng)速上，氣流對整個平板所做的無量綱總功也均為負(fù)。因此，在本文所考慮的較大折減風(fēng)速范圍內(nèi)，一個周期內(nèi)氣流對大幅度豎彎運(yùn)動的平板均做負(fù)功，也即消耗氣流的能量，因此單自由度大幅豎彎運(yùn)動的平板是氣動穩(wěn)定的。

3.4.2 單自由度大幅扭轉(zhuǎn)運(yùn)動（LAP）

考察大幅扭轉(zhuǎn)運(yùn)動的情況，圖11（a）和（b）分別為折減風(fēng)速Vr=8和Vr=40時，分別從多個氣動力周期時程中，取出一個完整周期（每個周期氣動力完全相同），該周期內(nèi)的氣動力系數(shù)時程曲線。與圖9相比，可見隨著折減風(fēng)速的提高，升力和扭矩系數(shù)間相位逐漸變化，且時程曲線波形發(fā)生畸變，也即失去諧波特性，特別是高折減風(fēng)速下氣動力時程曲線的畸變更顯著。對低和高折減風(fēng)速下氣動扭矩系數(shù)時程的歸一化功率譜分析表明，在低折減風(fēng)速下扭矩系數(shù)時程的頻率等于扭轉(zhuǎn)運(yùn)動頻率，如圖12（a）所示;但高折減風(fēng)速下扭矩系數(shù)時程雖然峰值頻率為扭轉(zhuǎn)運(yùn)動頻率，但頻譜圖中明顯出現(xiàn)了倍頻成分，如圖12（b）所示，且倍頻成分的能量明顯小于扭轉(zhuǎn)運(yùn)動頻率對應(yīng)峰值能量，可見出現(xiàn)的非線性不是特別顯著。因此，當(dāng)薄平板大幅扭轉(zhuǎn)運(yùn)動時，其氣動力系統(tǒng)不再為線性，且隨著折減風(fēng)速的提高，其非線性將越來越顯著。這可能是折減風(fēng)速增大，繞薄平板流動的分離加劇，因此平板表面的流態(tài)隨著折減風(fēng)速的提高在不斷改變，使得氣動力的非線性變得明顯。

圖13給出了一個扭轉(zhuǎn)運(yùn)動周期內(nèi)氣流對平板所做無量綱總功隨折減風(fēng)速的變化，可見當(dāng)折減風(fēng)速小于16時，一個周期內(nèi)氣流對平板所做總功為負(fù)，也即平板消耗氣流的能量;但當(dāng)折減風(fēng)速大于20后，氣流對上游側(cè)半個平板估的正功將大于對下游側(cè)半個平板所做負(fù)功，氣流對平板所做的總功由負(fù)變?yōu)檎藭r一個周期內(nèi)扭轉(zhuǎn)運(yùn)動平板將從氣流中吸收能量，表明大幅扭轉(zhuǎn)運(yùn)動平板已進(jìn)入氣動不穩(wěn)定狀態(tài)。

3.5 兩自由度耦合運(yùn)動

從上面研究可知，當(dāng)扭轉(zhuǎn)運(yùn)動幅值角為20°且折減風(fēng)速較高后，單自由度大幅扭轉(zhuǎn)運(yùn)動的薄平板將進(jìn)入氣動不穩(wěn)定狀態(tài)。下面針對扭轉(zhuǎn)運(yùn)動幅值角20°的情況，考察不同豎彎運(yùn)動幅值平板兩自由度耦合系統(tǒng)的氣動特征，豎彎幅值分別采用小幅、中幅和大幅，對應(yīng)的幅值分別是B/50，B/16和B/4，并設(shè)定耦合運(yùn)動的豎彎和扭轉(zhuǎn)自由度相位角為0o（CFD可設(shè)定任意相位角）。分別開展不同折減風(fēng)速的CFD模擬，比較作用在平板上的氣動力和一周期內(nèi)氣流對平板做功的能量特征。

3.5.1 小幅豎彎大幅扭轉(zhuǎn)（SAH?LAP）

圖14（a）和（b）分別是折減風(fēng)速Vr=10和Vr=40時，從各自多個氣動力周期時程中，取出一個完整周期（每個周期氣動力完全相同），對應(yīng)小幅豎彎大幅扭轉(zhuǎn)運(yùn)動平板一個周期內(nèi)的氣動力系數(shù)時程曲線。可見耦合運(yùn)動平板的氣動力系統(tǒng)表現(xiàn)出非線性，且隨著折減風(fēng)速的提高，其非線性特性越來越顯著，其特征與圖11相似。

圖15給出了兩個折減風(fēng)速下耦合運(yùn)動薄平板，在一個周期內(nèi)氣流在運(yùn)動平板上表面各點(diǎn)所做的無量綱總功。可見每個自由度方向氣流作用的能量特征與單自由度類似，也即豎彎自由度方向平板表面點(diǎn)均消耗氣流能量，但扭轉(zhuǎn)自由度方向是迎風(fēng)側(cè)上半個平板吸收氣流能量，下游側(cè)半個平板消耗氣流能量。當(dāng)折減風(fēng)速較高時，上游側(cè)半個平板吸收的能量明顯大于下游側(cè)半個平板消耗的能量，導(dǎo)致扭轉(zhuǎn)自由度方向吸收氣流能量。如綜合兩個自由度方向的總能量，在低折減風(fēng)速下豎彎自由度消耗氣流的能量大于扭轉(zhuǎn)自由度吸收的能量，因此平板是氣動穩(wěn)定的;但當(dāng)折減風(fēng)速較高后，豎彎自由度消耗氣流的能量小于扭轉(zhuǎn)自由度方向吸收的能量，此時兩自由度耦合運(yùn)動系統(tǒng)是氣動不穩(wěn)定的，如圖16所示。從圖15也可見，由于豎彎自由度振幅小，其消耗的氣流能量小，因此氣流在扭轉(zhuǎn)自由度方向作用的能量主導(dǎo)了兩自由度耦合運(yùn)動系統(tǒng)氣流作用的總能量。

3.5.2 中幅豎彎大幅扭轉(zhuǎn)（MAH?LAP）

維持薄平板兩自由度運(yùn)動的扭轉(zhuǎn)幅值角為20°，并將豎彎運(yùn)動位移幅值增大到B/16。圖17給出了折減風(fēng)速分別是Vr=10和Vr=40時，分別從多個氣動力周期時程中，取出一個完整周期（每個周期氣動力完全相同），對應(yīng)一個周期內(nèi)薄平板的升力和扭矩系數(shù)時程?？梢娂词沟驼蹨p風(fēng)速，薄平板氣動力的諧波特性變差，而隨著折減風(fēng)速的提高，氣動力時程的非線性特性越加顯著。

圖18給出了高低兩個折減風(fēng)速下耦合運(yùn)動的薄平板，一個周期內(nèi)氣流在運(yùn)動平板上表面各點(diǎn)所做的無量綱總功。其曲線特征與圖15相似，也即豎彎自由度方向平板表面所有點(diǎn)均消耗氣流能量，但扭轉(zhuǎn)自由度方向是迎風(fēng)側(cè)上半個平板吸收氣流能量，下游側(cè)半個平板消耗氣流能量。但由于豎彎自由度方向的振幅增大，豎彎運(yùn)動消耗氣流的能量明顯增大，從兩個自由度方向的能量總和來看，在所有折減風(fēng)速下，平板均消耗氣流的能量，如圖19所示。因此，以中幅豎彎大幅扭轉(zhuǎn)運(yùn)動的薄平板是氣動穩(wěn)定的。與前述小幅豎彎大幅扭轉(zhuǎn)運(yùn)動情況相比，可見豎彎幅度的增大提高了平板的氣動穩(wěn)定性。

3.5.3 兩自由度耦合大幅運(yùn)動（LAH?LAP）

進(jìn)一步增大豎彎運(yùn)動幅值到B/4，使得兩個自由度方向均為大幅運(yùn)動，也即大幅耦合運(yùn)動。圖20給出了折減風(fēng)速分別是Vr=10和Vr=40時，分別從多個氣動力周期時程中，取出一個完整周期（每個周期氣動力完全相同），對應(yīng)薄平板一個周期內(nèi)大幅耦合運(yùn)動的升力和扭矩系數(shù)時程?？梢姎鈩恿Φ姆蔷€性明顯，也隨折減風(fēng)速的增大而顯著。與圖17相比，相同的折減風(fēng)速下氣動力的非線性并沒有明顯變化。

圖21為低折減風(fēng)速（Vr=10）和高折減風(fēng)速（Vr=40）下大幅耦合運(yùn)動的薄平板，一個運(yùn)動周期內(nèi)氣流在運(yùn)動平板上表面各點(diǎn)所做的無量綱總功。豎彎和扭轉(zhuǎn)運(yùn)動氣流作用能量沿平板表面的分布與圖18稍有不同，可能是兩自由度大幅耦合運(yùn)動導(dǎo)致。但因豎彎自由度方向的大幅運(yùn)動，導(dǎo)致豎彎自由度方向消耗氣流的能量大幅度增加，使得一個周期內(nèi)氣流對平板表面所有點(diǎn)均做負(fù)功，因此氣流作用在運(yùn)動平板上的總能量為負(fù)，并明顯大于圖16小幅豎彎大幅扭轉(zhuǎn)和圖19中幅豎彎大幅扭轉(zhuǎn)運(yùn)動的情況，如圖22所示，因此大幅耦合運(yùn)動的平板是氣動穩(wěn)定的，原因是豎彎自由度方向的振幅增大，顯著提高了其氣動穩(wěn)定性。

4 結(jié) ?論

本文基于任意拉格朗日?歐拉描述法，開展了薄平板豎彎和扭轉(zhuǎn)強(qiáng)迫運(yùn)動繞流場CFD模擬，研究了薄平板大幅運(yùn)動的氣動特征和穩(wěn)定性，得到下述結(jié)論：

1） 單自由度小幅運(yùn)動的薄平板，其氣動力系統(tǒng)是線性的，也是氣動穩(wěn)定的。

2） 單自由度大幅豎彎運(yùn)動薄平板是氣動穩(wěn)定的，且氣動力系統(tǒng)非線性不明顯;但單自由度大幅扭轉(zhuǎn)運(yùn)動薄平板氣動力系統(tǒng)是非線性的，其非線性隨折減風(fēng)速提高而顯著，且折減風(fēng)速較高后，其氣動力系統(tǒng)將變得不穩(wěn)定。

3） 大幅扭轉(zhuǎn)的耦合運(yùn)動平板，氣動力系統(tǒng)均為非線性，其非線性主要來自扭轉(zhuǎn)自由度，非線性隨折減風(fēng)速的提高更顯著。

4） 大幅扭轉(zhuǎn)的耦合運(yùn)動平板氣動穩(wěn)定性，取決于豎彎自由度消耗的氣流能量與扭轉(zhuǎn)自由度吸收氣流能量的相對大小。當(dāng)豎彎振幅較小時和折減風(fēng)速較高時，氣動力系統(tǒng)是不穩(wěn)定的;當(dāng)豎彎自由度方向振幅較大時，其氣動力系統(tǒng)將是穩(wěn)定的。

參考文獻(xiàn)：

[1] Zhu Q， Xu Y L， Zhu L D， et al. Vortex-induced vibration analysis of long-span bridges with twin-box decks under non-uniformly distributed turbulent winds[J]. Journal of Wind Engineering & Industrial Aerodynamics， 2018，172： 31-41.

[2] Malík Josef. Sudden lateral asymmetry and torsional oscillations in the riginal Tacoma suspension bridge[J]. Journal of Sound and Vibration， 2013，332（15）：3772-3789.

[3] 祝志文.超臨界雷諾數(shù)下拉索順風(fēng)向自激力特性研究[J].振動工程學(xué)報，2014， 27（3）：377-384.

Zhu Zhiwen. Characteristics of along wind self-excited forces of a cable under super-critical Reynolds number[J]. Journal of Vibration Engineering， 2014， 27（3）：377-384.

[4] 祝志文，陳 魏，李健朋，等.多塔斜拉橋加勁索渦激振動實(shí)測與時域解析模態(tài)分解研究[J].中國公路學(xué)報，2019，32（10）：247-256.

Zhu Zhiwen. CHEN Wei， LI Jianpeng，et al. Field observation of vortex-induced vibration of stiffening cables in a multi-tower cable-stayed bridge with application of analytical mode decomposition[J]. China Journal of Highway and Transport， 2019，32（10）：247-256.

[5] Wu T， Kareem A. A nonlinear analysis framework for bluff-body aerodynamics： A Volterra representation of the solution of Navier-Stokes equations[J]. Journal of Fluids and Structures， 2015，54：479-502.

[6] 祝志文，石亞光，李健朋.扁平箱梁基于Volterra理論的氣動力非線性特性研究[J].中國公路學(xué)報，2018， 31 （1）： 74-81.

Zhu Zhiwen， Shi Yaguang， Li Jianpeng. Research on nonlinear aerodynamics of flat steel box girder based on Volterra theory[J]. China Journal of Highway and Transport ， 2018， 31 （1）： 74-81.

[7] 祝志文，顧 明，陳政清.兩自由度系統(tǒng)氣彈響應(yīng)特性研究[J].工程力學(xué)，2008，25（8）：23-30.

Zhu Zhiwen， Gu Ming， Chen Zhengqing. Study on response characteristics of two-DOF aeroelastic system[J]. Engineering Mechanics， 2008， 25（8）：23-30.

[8] Filipovic Nenad， Mijailovic Srboljub， Tsuda Akira， et al. An implicit algorithm within the arbitrary Lagrangian-Eulerian formulation for solving incompressible fluid flow with large boundary oscillations[J]. Computer Methods in Applied Mechanics and Engineering， 2006，195（44-47）： 6347-6361.

[9] Zhu Zhiwen， Chen Zhengqing， Gu Ming. CFD based simulations of flutter characteristics of ideal thin plates with and without central slot[J]. Wind and Structures， 2009， 12（1）： 1-19.

[10] Tau E Y. A second-order Projection method for incompressible Navier-Stokes equations in arbitrary domains[J]. Journal of Computational Physics， 1994， 115：147-152.

[11] 祝志文，陳政清.數(shù)值模擬橋梁斷面的顫振導(dǎo)數(shù)和顫振臨界風(fēng)速[J].中國公路學(xué)報，2004，15（4）：41-50.

Zhu Zhiwen， Chen Zhengqing. Numerical simulations for aerodynamic derivatives and critical flutter velocity of bridge deck[J]. China Journal of Highway and Transport， 2004， 15（4）：41-50.

[12] 祝志文，顧 明，陳政清.基于3211多階躍激勵CFD模型的顫振導(dǎo)數(shù)識別研究[J].振動工程學(xué)報，2008， 21（1）：18-23.

Zhu Zhiwen， Gu Ming， Chen Zhengqing. Identification of flutter derivatives based on 3211 multi-step excitation of CFD model[J]. Journal of Vibration Engineering， 2008，21（1）：18-23.

[13] Scanlan R H， Tomko J J. Airfoil and bridge deck flutter derivatives[J]. Journal of Engineering Mechanics， ASCE， 1971，97（6）： 1171-1737.

CFD investigation on aerodynamics and stability of a two degree-of-freedom thin plate undergoing large-amplitude oscillation

ZHU Zhi-wen1， YAN Shuang1， WANG Qin-hua1， LI Jia-wu2

（1. Key Laboratory of Structure and Wind Tunnel of Guangdong Higher Education Institutes （Shantou University），

Shantou 515063， China; 2. Highway School， Chang'an University， Xi'an 710064， China）

Abstract： In order to investigate the aerodynamic characteristics and aerodynamic stability of the two degree-of-freedom （2-DOF） thin plate undergoing large-amplitude oscillation， the governing equations for arbitrary-deformed fluid domain in the two-dimensional incompressible form are numerical solved based on the Arbitrary-Lagrangian-Eulerian description and the finite difference method in moving-grid system. Then， computational fluid dynamics （CFD） is employed to simulate the flow field around the thin plate undergoing heave and/or pitch oscillation at various reduced wind speeds. The research finds that the aerodynamic system around the thin plate is linear and stability when the thin plate experiences small-amplitude vibration in heave or pitch， and it is also linear and stable even when the thin plate undergoes large-amplitude heave motion. However， when the thin plate undergoes large-amplitude pitch motion， the aerodynamic system of the thin plate presents nonlinearity which will become significant and even turn into instability with the increase of reduced wind speed. Meanwhile， when the thin plate experiences 2-DOF oscillation with large amplitude in pitch coupled with different amplitude in heave， nonlinearity will also be presented and be more significant with the increase of reduced wind speed. It is noted that such nonlinearity is generated from the large-amplitude oscillation in pitch. For the coupled vibration system with 2-DOF， when the vibration amplitude is small while the reduced wind speed is high， its aerodynamic system is instable. While when the vibration amplitude in heave is large， the aeroelastic system will be stable.

Key words： aerodynamics stability; thin plate; large-amplitude oscillation; CFD; aerodynamic nonlinearity

作者簡介： 祝志文（1968-），男，教授，博士生導(dǎo)師。電話： 13574876655; E-mail： zhuzw@shu.edu.cn