徐萬(wàn)海 ， 楊 猛 ，， 蘆 燕

（1.天津大學(xué) 水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津 300072；2.天津大學(xué) 建筑工程學(xué)院，天津300072）

0 引 言

渦激振動(dòng)（Vortex-induced vibration，VIV）是引起海洋立管、海上浮式風(fēng)機(jī)系泊纜等大長(zhǎng)徑比柔性圓柱結(jié)構(gòu)疲勞損傷的重要因素。如何抑制VIV是工程界與學(xué)術(shù)界長(zhǎng)期關(guān)注的焦點(diǎn)問(wèn)題。渦激振動(dòng)的抑制方式主要有兩種，即主動(dòng)控制和被動(dòng)抑制。主動(dòng)控制通過(guò)監(jiān)測(cè)流場(chǎng)和結(jié)構(gòu)受力，運(yùn)用計(jì)算機(jī)自動(dòng)控制技術(shù)，引入適當(dāng)?shù)臄_動(dòng)改變流場(chǎng)。但其操作復(fù)雜，成本較高，技術(shù)難度大，工程實(shí)際中的應(yīng)用受到較大限制。被動(dòng)抑制直接改變結(jié)構(gòu)表面形狀或安裝附加裝置，改變結(jié)構(gòu)周圍的流場(chǎng)，這種方式在實(shí)際工程中的應(yīng)用十分廣泛[1-6]。

根據(jù)機(jī)理的不同，被動(dòng)抑制裝置可分為3類：（1）表面突出結(jié)構(gòu)，如螺旋列板（helical strakes）、突起（bumps）等；（2）尾流控制裝置，如整流罩（fairings）、分隔板（splitter plate）等；（3）覆蓋結(jié)構(gòu)，如軸向板條（axial rod shroud）、控制桿（control rod）等[7]。 上述VIV抑制裝置各有優(yōu)劣。 Vandiver等[8]開(kāi)展了現(xiàn)場(chǎng)實(shí)驗(yàn)，觀測(cè)了螺旋列板的VIV抑制效果，發(fā)現(xiàn)螺旋列板在減小高階模態(tài)振動(dòng)幅值的同時(shí)還能降低三次諧波成分，但同時(shí)也增大了結(jié)構(gòu)順流向受到的拖曳力。Allen[9]通過(guò)模型實(shí)驗(yàn)觀測(cè)發(fā)現(xiàn)：整流罩對(duì)于渦激振動(dòng)和拖曳力的抑制效果明顯，但流速過(guò)大時(shí)，立管將發(fā)生很大的橫向振動(dòng)，存在動(dòng)態(tài)不穩(wěn)定性。Akilli等[10]實(shí)驗(yàn)發(fā)現(xiàn)分隔板對(duì)VIV具有較好的抑制效果，但只對(duì)單個(gè)來(lái)流方向起作用。Tsutsui和Igarashi[11]在實(shí)驗(yàn)中發(fā)現(xiàn)單根控制桿能夠降低主管的脈動(dòng)升力和拖曳力，并優(yōu)化了控制桿的直徑d=0.25D（D為主管直徑），但單根控制桿僅對(duì)某些來(lái)流方向效果明顯。宋吉寧等[12]的實(shí)驗(yàn)研究表明：三根控制桿抑制裝置對(duì)來(lái)流方向適應(yīng)性很強(qiáng)。Wu等[13]開(kāi)展附有四根控制桿的海洋立管渦激振動(dòng)實(shí)驗(yàn)，研究了覆蓋率和間距（控制桿外表面到立管模型外表面的距離）對(duì)抑制效果的影響，表明覆蓋率越高，間距越小，抑制效果越好。

1 實(shí)驗(yàn)方案

實(shí)驗(yàn)觀測(cè)在天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室的拖曳水池中完成，水池長(zhǎng)137.0 m，寬7.0 m，深3.3 m。開(kāi)展了三組實(shí)驗(yàn)：（1）柔性圓柱模型渦激振動(dòng)實(shí)驗(yàn)；（2）附有三根控制桿抑制裝置的柔性圓柱模型渦激振動(dòng)實(shí)驗(yàn)；（3）附有四根控制桿抑制裝置的柔性圓柱模型渦激振動(dòng)實(shí)驗(yàn)。三組實(shí)驗(yàn)中的柔性圓柱模型兩端通過(guò)萬(wàn)向節(jié)固定在拖車底部的剛架上，兩端施加的預(yù)張力均為450 N。三組實(shí)驗(yàn)分別進(jìn)行，使用相同的柔性圓柱模型。拖曳水池上方的拖車拖動(dòng)模型勻速前進(jìn)，可模擬均勻來(lái)流條件，待圓柱振動(dòng)達(dá)到穩(wěn)定狀態(tài)時(shí)，開(kāi)始數(shù)據(jù)采集，采樣頻率100 Hz，采集時(shí)間50 s。實(shí)驗(yàn)中拖車的速度范圍為0.05～1.00 m/s，間隔為 0.05 m/s，共計(jì)測(cè)量 60個(gè)實(shí)驗(yàn)工況，實(shí)驗(yàn)中雷諾數(shù)最大可達(dá)到16 000。

圖1 實(shí)驗(yàn)?zāi)Ｐ虵ig.1 Experimental models

柔性圓柱模型的內(nèi)芯為銅管，沿軸向均勻設(shè)置七個(gè)測(cè)點(diǎn)，每個(gè)測(cè)點(diǎn)的橫流向和順流向布置應(yīng)變片，用于采集振動(dòng)信息。模型外層為硅膠管，以保證外表面光滑，并防水絕緣。控制桿模型采用外徑為0.25D的圓形截面塑膠條，間距為0.5D，固定在塑料材質(zhì)的支撐裝置上。支撐裝置有兩種，分別可以等間距地固定三根控制桿和四根控制桿。圖1為上述三組實(shí)驗(yàn)中的模型?？刂茥U與柔性圓柱模型的排布方式如圖2所示。制作控制桿和支撐裝置的材料在水中的比重近似等于1，對(duì)柔性圓柱模型整體質(zhì)量比的影響很小。詳細(xì)的圓柱模型結(jié)構(gòu)參數(shù)如表1所示。

表1 光滑圓柱模型參數(shù)Tab.1 Physical property for the cylinder model

圖2 控制桿與圓柱模型的布置圖Fig.2 Sketch of cylinder model with control rods

采用模態(tài)法處理采集得到的應(yīng)變信息。模態(tài)法可以通過(guò)柔性圓柱上若干個(gè)離散點(diǎn)的應(yīng)變信息獲得整個(gè)結(jié)構(gòu)上任一點(diǎn)的位移信息[14-15]，以橫流向?yàn)槔榻B模態(tài)法。

建立如下坐標(biāo)系：x軸對(duì)應(yīng)順流向，y軸對(duì)應(yīng)橫流向，z軸對(duì)應(yīng)柔性圓柱模型軸向。對(duì)于長(zhǎng)度為L(zhǎng)的柔性圓柱振動(dòng)問(wèn)題，其橫流向振動(dòng)位移用模態(tài)疊加的方法可以表示為：

其中：y（ z,t）為柔性圓柱橫流向振動(dòng)位移，Wn（t）為權(quán)重函數(shù)，φn（z）為模態(tài)函數(shù)，n為柔性圓柱振動(dòng)模態(tài)階次。本實(shí)驗(yàn)中柔性圓柱模型兩端均為簡(jiǎn)支邊界條件，因此模態(tài)函數(shù)φn（z）可表示為：

若在柔性圓柱模型上布置M個(gè)應(yīng)變片，其振動(dòng)位移可以采用N階模態(tài)的線性疊加來(lái)表示，則橫流向的曲率矩陣Py可以用下式表達(dá)：

其中：

其中：ε（ z,t）為測(cè)量的應(yīng)變信號(hào)，R為內(nèi)芯銅管半徑。實(shí)驗(yàn)中，柔性圓柱模型橫方向激發(fā)的最高階模態(tài)為4，即N=4，應(yīng)變片的布置數(shù)M=7，所以N＜M，權(quán)重函數(shù)可以采用最小二乘法來(lái)求得：

圖4(a)、(b)分別為改進(jìn)型ADRC和雙閉環(huán)PID對(duì)小車位置的控制曲線，圖4(c)、(d)分別是改進(jìn)型ADRC和雙閉環(huán)PID對(duì)小車位置的控制曲線。通過(guò)對(duì)比可以看出改進(jìn)型ADRC在2s以后就可以穩(wěn)定控制擺桿角度和小車位置，而雙閉環(huán)PID在3s以后失去對(duì)擺桿和小車位置的控制；通過(guò)仿真可以看出基于改進(jìn)型的ADRC的一級(jí)直線倒立擺系統(tǒng)具有較好的魯棒性。

將（7）式代入（1）式即可求出柔性圓柱模型的橫流向振動(dòng)位移。順流向的振動(dòng)位移確定方式與橫流向類似，在此不再贅述。

2 抑制效果分析

控制頻率是反映柔性圓柱渦激振動(dòng)特性的重要信息。圖3和圖4給出了三組實(shí)驗(yàn)中圓柱模型渦激振動(dòng)橫流向（CF）和順流向（IL）的控制頻率隨流速變化情況?？梢钥闯?，三種模型的控制頻率均隨流速線性增加。附有控制桿的圓柱模型控制頻率顯著低于不帶控制桿的光滑圓柱，四根控制桿的控制頻率比三根控制桿低，順流向尤為顯著。對(duì)比圖3和圖4可發(fā)現(xiàn)，光滑圓柱和附有三根控制桿的抑制圓柱順流向控制頻率大約為橫流向的2倍，而附有四根控制桿的抑制圓柱順流向和橫流向的控制頻率幾乎一致。當(dāng)流速為0.95 m/s和1.0 m/s時(shí)，三根控制桿圓柱模型的橫流向渦激振動(dòng)控制頻率出現(xiàn)了突然下降，產(chǎn)生該現(xiàn)象的原因?yàn)榭刂颇B(tài)由高階轉(zhuǎn)化為低階，導(dǎo)致了主控頻率的降低?？偠灾?，多根控制桿能降低圓柱模型的控制頻率，四根控制桿的效果優(yōu)于三根控制桿，前者能極大地降低結(jié)構(gòu)順流向的控制頻率。

圖3 橫流向控制頻率隨來(lái)流速度變化Fig.3 CF dominate frequency versus flow velocity

圖4 順流向控制頻率隨來(lái)流速度變化Fig.4 IL dominate frequency versus flow velocity

圖5 橫流向控制模態(tài)隨來(lái)流速度變化Fig.5 CF dominate mode versus flow velocity

圖6 順流向控制模態(tài)隨來(lái)流速度變化Fig.6 IL dominate mode versus flow velocity

柔性圓柱模型渦激振動(dòng)特性的另一重要參數(shù)是控制模態(tài)。圖5和圖6反映了三種實(shí)驗(yàn)?zāi)Ｐ蜏u激振動(dòng)的控制模態(tài)。從圖中可以看出，所有模型橫流向和順流向最低控制模態(tài)均為1階。光滑圓柱橫流向和順流向的最高控制模態(tài)分別為4階和6階，附有三根控制桿的模型為3階和4階，而附有四根控制桿模型橫流向與順流向均為3階?？梢?jiàn)控制桿能降低圓柱模型渦激振動(dòng)激發(fā)的最高模數(shù)，對(duì)順流向的抑制效果更加顯著，四根控制桿對(duì)控制模態(tài)的抑制效果略優(yōu)于三根控制桿。

應(yīng)變是實(shí)驗(yàn)中直接測(cè)量得到的信號(hào)，與結(jié)構(gòu)受到的應(yīng)力直接相關(guān)，是反映結(jié)構(gòu)渦激振動(dòng)疲勞破壞的重要指標(biāo)。圖7和圖8為三組實(shí)驗(yàn)中圓柱模型的最大應(yīng)變均方根值。圖9和圖10給出了來(lái)流速度為0.5 m/s時(shí)，圓柱模型測(cè)點(diǎn)處的應(yīng)變時(shí)間歷程曲線和對(duì)應(yīng)的頻譜分析圖?？梢钥闯?，控制桿顯著減小了圓柱模型渦激振動(dòng)的應(yīng)變，對(duì)順流向的抑制效果更為明顯，四根控制桿的作用效果比三根控制桿略好。從頻譜圖中可以看出，光滑圓柱的應(yīng)變響應(yīng)頻率最高，三根控制桿和四根控制桿響應(yīng)頻率差異不明顯，三根控制桿略高。光滑圓柱的應(yīng)變響應(yīng)頻率出現(xiàn)了幅值較高的高頻成分，控制桿能降低圓柱模型的應(yīng)變響應(yīng)頻率，并抑制高頻成分的幅值。

圖7 橫流向最大應(yīng)變均方根隨來(lái)流速度變化Fig.7 Max RMS CF strain versus flow velocity

圖8 順流向最大應(yīng)變均方根隨來(lái)流速度變化Fig.8 Max RMS IL strain versus flow velocity

圖9 測(cè)點(diǎn)處橫流向應(yīng)變時(shí)間歷程曲線和頻譜分析結(jié)果（黑色為光滑圓柱，藍(lán)色為三根控制桿，紅色為四根控制桿，V=0.5 m/s）Fig.9 CF strain history and its spectra at measured positions(Black for bare,blue for 3 control rods,red for 4 control rods,V=0.5 m/s)

圖10 測(cè)點(diǎn)處順流向應(yīng)變時(shí)間歷程曲線和頻譜分析結(jié)果（黑色為光滑圓柱，藍(lán)色為三根控制桿，紅色為四根控制桿，V=0.5 m/s）Fig.10 IL strain history and its spectra at measured positions(Black for bare,blue for three control rods,red for four control rods,V=0.5 m/s)

圖11和圖12為三組實(shí)驗(yàn)中模型的最大位移均方根。圖13和圖14為來(lái)流速度為0.5m/s時(shí)，圓柱模型測(cè)點(diǎn)處的位移時(shí)間歷程曲線和對(duì)應(yīng)的頻譜分析圖。在流速低于0.35 m/s時(shí)，附有控制桿的圓柱模型順流向響應(yīng)位移比高流速范圍時(shí)顯著增大，這是不合理的，其原因在于模型的響應(yīng)位移是通過(guò)模態(tài)法求得的，低階模態(tài)振動(dòng)對(duì)應(yīng)的應(yīng)變信號(hào)非常小，對(duì)測(cè)量噪聲比較敏感，模態(tài)法很難辨別低階模態(tài)響應(yīng)的真實(shí)信號(hào)和噪聲誤差，來(lái)流速度較低時(shí)測(cè)量噪聲會(huì)將較大的誤差引入位移計(jì)算中，圓柱模型安裝控制桿后，低流速下的噪聲干擾比光滑圓柱時(shí)更加明顯[16]?；谏鲜龇治?，可以認(rèn)為來(lái)流速度低于0.35 m/s時(shí)，附有多根控制桿圓柱模型的計(jì)算位移主要是由噪聲誤差造成的，在分析其位移時(shí)僅考慮流速大于0.35 m/s的結(jié)果。

可以看出，控制桿能顯著減小柔性圓柱模型的響應(yīng)幅值，四根控制桿的響應(yīng)幅值比三根控制桿小。從頻率分析圖中可以看出，光滑圓柱的位移響應(yīng)頻率最高，四根控制桿最低，三根控制桿位于兩者之間。光滑圓柱的位移響應(yīng)頻率出現(xiàn)了高頻成分，而控制桿對(duì)高頻成分有抑制效果。

圖14 測(cè)點(diǎn)處順流向位移時(shí)間歷程曲線和頻譜分析結(jié)果（黑色為光滑圓柱，藍(lán)色為三根控制桿，紅色為四根控制桿，V=0.5 m/s）Fig.14 IL displacement history and its spectra at measured positions(Black for bare,blue for three control rods,red for four control rods,V=0.5 m/s)

3 結(jié) 論

本文開(kāi)展了三組高雷諾數(shù)條件下柔性圓柱模型渦激振動(dòng)實(shí)驗(yàn)，對(duì)柔性圓柱模型、附有三根控制桿和四根控制桿的圓柱模型的渦激振動(dòng)特性進(jìn)行了觀測(cè)，通過(guò)對(duì)控制頻率、控制模態(tài)、結(jié)構(gòu)應(yīng)變、響應(yīng)位移等方面對(duì)比分析，分析了控制桿對(duì)柔性圓柱渦激振動(dòng)的抑制效果，得到了如下結(jié)論：

（1）控制桿對(duì)柔性圓柱渦激振動(dòng)抑制效果明顯，對(duì)順流向的抑制作用更加突出；

（2）總體而言，四根控制桿的抑制效果優(yōu)于三根控制桿；

（3）控制桿能抑制柔性圓柱渦激振動(dòng)響應(yīng)頻率的高頻成分。

實(shí)際工程中，海洋立管的周圍會(huì)部署若干管徑和剛度相對(duì)于立管均很小的附屬管，如在設(shè)計(jì)時(shí)合理確定附屬管與立管的空間位置就相當(dāng)于安裝控制桿并優(yōu)化其排布。排布合理的多根控制桿對(duì)流向適應(yīng)性強(qiáng)，對(duì)抑制大長(zhǎng)徑比柔性圓柱渦激振動(dòng)具有很廣闊的前景。因此，有待進(jìn)一步深入研究控制桿的空間排布對(duì)柔性圓柱VIV抑制的影響規(guī)律。