徐萬海,馬燁璇, 羅浩, 欒英森, 吳夢寧

(天津大學 水利工程仿真與安全國家重點實驗室, 天津 300072)

小傾角傾斜柔性圓柱渦激振動實驗研究

徐萬海,馬燁璇, 羅浩, 欒英森, 吳夢寧

(天津大學 水利工程仿真與安全國家重點實驗室, 天津 300072)

海洋工程中柔性圓柱結構疲勞破壞的主要誘因是海流引起的渦激振動。垂直來流作用下柔性圓柱渦激振動研究已取得諸多成果，然而在實際工程中，圓柱軸向與海流并不完全垂直，存在一定傾斜角度。針對這種普遍現(xiàn)象，學術界提出了傾斜圓柱渦激振動的不相關原則，但其正確與否至今存在爭議。本文設計了室內(nèi)模型實驗，開展了長徑比350，雷諾數(shù)800～16 000，質量比1.9，傾斜角度為15°的柔性圓柱渦激振動實驗，通過對比垂直圓柱與小傾角傾斜圓柱的實驗結果發(fā)現(xiàn)，在模態(tài)轉化區(qū)域外不相關原則適用于小傾角傾斜柔性圓柱渦激振動分析，同時觀測到在模態(tài)轉化區(qū)傾斜圓柱更容易激發(fā)高階模態(tài)。

傾斜柔性圓柱; 渦激振動; 不相關原則; 流-固耦合

海洋油氣開采及海洋可再生能源開發(fā)利用過程中，海底懸空管道、海洋立管以及海洋浮式風機系泊纜等柔性圓柱結構十分常見。在海流的作用下，圓柱結構兩側會出現(xiàn)交替的漩渦脫落，產(chǎn)生周期變化的流體力，促使結構振動，結構振動又反作用于流場，這種典型的流-固耦合現(xiàn)象稱為“渦激振動(vortex-induced vibrations, VIV)”。近年來，人們針對垂直來流作用下柔性圓柱渦激振動開展了廣泛的理論和實驗研究，并取得了大量成果[1-4]。

實驗是探究柔性圓柱渦激振動現(xiàn)象最有效的手段之一。Trim等[5]開展了垂直均勻流作用的柔性立管渦激振動實驗，發(fā)現(xiàn)順流向渦激振動對結構疲勞損傷的貢獻與橫流向同等重要。Song等[6]設計了大長徑比柔性立管渦激振動實驗，研究表明柔性圓柱立管振動存在多個模態(tài)，并存在3倍、4倍和5倍頻率等高頻現(xiàn)象。Huera-Huarte等[7]開展了垂直階梯流作用下低質量比(1.1、2.7)柔性圓柱渦激振動實驗，觀測發(fā)現(xiàn)低質量比情況下柔性圓柱響應幅值顯著增大，橫流向與順流向最大幅值可分別達到3倍和1.5倍結構直徑。

上述研究的來流條件均與圓柱軸向垂直。實際工程應用中，圓柱結構與來流方向并不完全垂直，常存在不同的傾斜角度。為了簡化問題，學術界提出了不相關原則[8-9]，認為傾斜圓柱渦激振動問題等效為與圓柱軸向垂直的來流速度分量引起的渦激振動。現(xiàn)階段，傾斜圓柱渦激振動研究才剛剛起步，不相關原則的合理性仍存在爭議。Franzini等[10]開展了傾斜剛性圓柱的渦激振動實驗，結果表明：傾角小于20°時，不相關原則成立，傾角大于20°時則不完全成立。Bourguet等[11]采用直接數(shù)值模擬研究了傾角為60°，長徑比為50柔性圓柱渦激振動，發(fā)現(xiàn)在低軸向力條件時傾斜圓柱渦激振動特性與垂直圓柱差異較大，不相關原則不合理，高軸向力時不相關原則合理。

由于軸向二次流的存在，柔性圓柱在傾斜來流條件下的渦激振動十分復雜，已開展的研究不能完全反映柔性圓柱在傾斜來流作用下的渦激振動行為，缺乏大長徑比、小質量比、高雷諾數(shù)條件下的傾斜圓柱渦激振動實驗研究。基于此，本文設計開展了小傾角傾斜柔性圓柱渦激振動實驗，目的是探究不相關原則的適用性，并觀測小傾角傾斜圓柱渦激振動特性。

1 小傾角柔性圓柱VIV實驗

實驗在天津大學水利工程仿真與安全國家重點實驗室的拖曳水池中完成。圖1為實驗水池及裝置示意圖。拖曳水池長137.0 m，寬7.0 m，深3.3 m，柔性圓柱模型通過萬向節(jié)安裝在與拖車相連鋼架上，圓柱結構軸向力為500 N。調整鋼架和拖車的相對位置可實現(xiàn)垂直來流條件和15°傾角來流條件。實驗中拖車速度范圍為0.05～1.00 m/s，間隔為0.05 m/s，雷諾數(shù)范圍近似800～16 000。拖車采用模擬控制和數(shù)字控制兩套系統(tǒng)，可保證速度精度為0.1%，如圖1所示。

柔性圓柱模型內(nèi)芯材料為銅管，銅管沿軸向被7個截面等分8段，應變片粘貼于上述7個截面的橫流向和順流向，用于采集振動信息，應變片采用半橋接法，測量導線沿銅管端部引出。銅管外部包裹硅膠管，目的是使柔性圓柱模型外表面規(guī)則光滑，同時保護內(nèi)部的應變片和測量導線，并起到防水和絕緣的作用。

拖車勻速拖動柔性圓柱模型，實現(xiàn)均勻來流條件，待圓柱振動達到穩(wěn)定開始數(shù)據(jù)采集，采樣頻率為100 Hz，穩(wěn)定段采集時間為50 s。本文實驗完成了兩個工況：一為柔性圓柱模型軸向與來流方向垂直(即傾角為0)，二為柔性圓柱模型傾角為15°。共計開展40組觀測與測量。詳細的圓柱模型結構參數(shù)如表1所示。

表1 柔性圓柱模型參數(shù)

圖1 實驗示意圖Fig.1 Schematic of the experimental set-up

2 數(shù)據(jù)處理-模態(tài)法

運用模態(tài)法處理實驗采集的應變數(shù)據(jù)，模態(tài)法可以根據(jù)柔性圓柱模型上離散點應變信息獲得整個圓柱模型上任一點的位移信息[12-13]。下面以橫流向(CF)為例詳細說明模態(tài)法的運用過程，順流向(IL)與之類似不再贅述。

建立坐標系：x軸對應IL方向，y軸對應CF方向，z軸對應圓柱模型軸向(如圖1(b)所示)。圓柱模型CF方向的曲率Cy(z,t)可表示為

(1)

式中：t為時間，ε(z,t)為應變信號，R為圓柱模型的半徑。對于長度為L的圓柱振動問題，根據(jù)模態(tài)法可以將其CF振動位移表示為

(2)

式中：U(z,t)為圓柱CF方向振動位移，n為模態(tài)階次，L為圓柱模型總長度，An(t)為權重系數(shù)，φn(z)為模態(tài)函數(shù)。本實驗中柔性圓柱模型兩端為簡支條件，因此φn(z)為

(3)

根據(jù)力學關系，曲率Cy(z,t)還可表示為

(4)

若在圓柱模型上布置M個應變片，其振動位移采用N階模態(tài)的線性疊加來表示，式(4)可寫成矩陣形式：

(5)

其中，

(6)

(7)

(8)

數(shù)據(jù)處理時權重系數(shù)采用最小二乘法確定：

(9)

將式(9)代入式(2)，可得到CF方向振動位移。

3 小傾角柔性圓柱渦激振動特性

控制模態(tài)是描述柔性圓柱渦激振動特性的一個重要的指標。圖2和圖3分別繪制了圓柱橫流向和順流向振動控制模態(tài)隨約化速度變化情況。其中橫坐標為約化速度Vr=Vcos(α)/(f1D)，V是來流速度，a表示傾斜角度，垂直時為0，傾斜時為15°，f1為結構一階固有頻率，本實驗測量了圓柱結構靜水中的固有頻率，發(fā)現(xiàn)測量的結果與式(10)得到的理論結果誤差非常小?？紤]到柔性圓柱垂直和傾斜時測量結果存在微小差異，為了方便分析，統(tǒng)一采用式(10)得到的理論值計算約化速度：

(10)

從圖2和圖3中可發(fā)現(xiàn)無論是橫流向還是順流向渦激振動，垂直與傾斜圓柱兩者的控制模態(tài)變化規(guī)律基本一致。橫流向渦激振動最高控制模態(tài)為3階，順流向為5階。在不同控制模態(tài)轉化區(qū)域，如約化速度Vr在8.0、18.0附近范圍，相比于垂直柔性圓柱，傾斜柔性圓柱更容易被激發(fā)更高階模態(tài)。產(chǎn)生上述現(xiàn)象的原因是軸向二次流引起的，這也充分證明了傾斜柔性圓柱渦激振動的復雜性。

圖2 橫流向控制模態(tài)隨流速變化Fig.2 CF dominate mode number versus reduced velocity

圖3 順流向控制模態(tài)隨流速變化Fig.3 IL dominate mode number versus reduced velocity

圖4和圖5為最大位移均方根隨約化速度的變化圖像。

圖4 橫流向最大均方根位移幅值隨約化速度變化Fig.4 Max RMS CF displacement versus reduced velocity

圖5 順流向最大均方根位移幅值隨約化速度變化Fig.5 Max RMS IL displacement versus reduced velocity

由圖4、5可以看出在模態(tài)轉化區(qū)域外，垂直圓柱與傾斜圓柱橫流向和順流向的最大位移均方根基本接近。但在模態(tài)轉化區(qū)域，傾斜圓柱得到的響應幅值，無論是橫流向還是順流向均會更小。由于實驗所用的圓柱模型質量比為1.9，垂直圓柱和小傾角圓柱均出現(xiàn)了較大響應幅值現(xiàn)象，這與近期Huera-Huarte等[7]的研究成果一致。

柔性圓柱渦激振動的控制頻率是評估結構疲勞壽命的重要參數(shù)。橫流向與順流向控制頻率隨約化速度變化情況如圖6和圖7所示。觀測發(fā)現(xiàn)控制頻率隨約化速度變化呈線性增長，0與15°傾角工況下，圓柱渦激振動控制頻率變化規(guī)律基本相同。但在低階模態(tài)向高階模態(tài)轉化的區(qū)域，傾斜圓柱的控制頻率高于垂直圓柱，再次說明傾斜柔性圓柱更容易激發(fā)高階模態(tài)。

圖6 橫流向控制頻率隨約化速度變化Fig.6 CF dominate frequency versus reduced velocity

圖7 順流向控制頻率隨約化速度變Fig.7 IL dominate frequency versus reduced velocity

圖8給出了中等來流速度(V=0.5 m/s)時七個測點處橫流向位移時間歷程曲線和對應的頻譜圖。

注：實線為垂直圓柱，虛線為傾斜圓柱圖8 測點處橫流向位移時間歷程曲線和頻譜分析結果(V=0.5 m/s)Fig.8 CF displacement history and its spectra at measured positions (V=0.5 m/s)

V=0.5 m/s位于模態(tài)轉化區(qū)域外，垂直圓柱和傾斜圓柱的位移歷程曲線兩者之間幾乎無差異，證明在模態(tài)轉化區(qū)域范圍外，不相關原則對于小傾角傾斜柔性圓柱是適用的。從橫流向位移的頻譜分析圖中可以看出傾斜圓柱控制頻率為4.7 Hz，而垂直圓柱為5.0 Hz，均出現(xiàn)了2倍頻和3倍頻等高頻諧波成分。

圖9給出了較大來流速度(V=0.8 m/s)時垂直圓柱和傾斜圓柱在的測點處的橫流向位移時間歷程曲線和對應的頻譜圖。V=0.8 m/s時正好處于模態(tài)轉化區(qū)域，從位移歷程曲線可看出，垂直圓柱和傾斜圓柱測點處的位移差異較大，因此在模態(tài)轉化區(qū)域內(nèi)不相關原則并不能很好的適用于小傾角傾斜柔性圓柱。從頻譜分析圖中可以看出，傾斜圓柱橫流向的控制頻率為8.13 Hz，垂直圓柱橫流向的控制頻率為6.90 Hz，均出現(xiàn)了高頻的成分。

注：實線為垂直圓柱，虛線為傾斜圓柱圖9 測點處橫流向位移時間歷程曲線和頻譜分析結果(V=0.8 m/s)Fig.9 CF displacement history and its spectra at measured positions(V=0.8 m/s)

4 結論

本文對小傾角傾斜柔性圓柱渦激振動進行了實驗研究，采用了大長徑比、小質量比的圓柱實驗模型，雷諾數(shù)最大可達到16 000，通過與垂直柔性圓柱的渦激振動實驗結果對比得到了如下結論：

1)小傾角傾斜圓柱與垂直圓柱測點處的位移在模態(tài)轉化區(qū)域外差異很小，而在模態(tài)轉化區(qū)域內(nèi)差異較大，并且均發(fā)現(xiàn)了高頻的成分。

2)通過分析小傾角傾斜圓柱和垂直圓柱的響應幅值，發(fā)現(xiàn)在模態(tài)轉化區(qū)域外，兩者的最大位移均方根差異很小。由于柔性圓柱模型質量比較小，均出現(xiàn)了大響應幅值現(xiàn)象。

3)根據(jù)小傾角傾斜圓柱和垂直圓柱的控制模態(tài)、控制頻率、響應幅值、測點處的位移觀測結果，發(fā)現(xiàn)在模態(tài)轉化區(qū)域內(nèi)傾斜圓柱比垂直圓柱更容易激發(fā)高階模態(tài)，不相關原則在模態(tài)轉化區(qū)域范圍外適用于小傾角柔性圓柱的渦激振動。

[1]SARPKAYA T. A critical review of the intrinsic nature of vortex-induced vibrations[J]. Journal of fluids and structures, 2004, 19(4): 389-447.

[2]GABBAI R D, BENAROYA H. An overview of modeling and experiments of vortex-induced vibration of circular cylinders[J]. Journal of sound and vibration, 2005, 282(3/4/5): 575-616.

[3]WILLIAMSON C H K, GOVARDHAN R. A brief review of recent results in vortex-induced vibrations[J]. Journal of wind engineering and industrial aerodynamics, 2008, 96(6/7): 713-735.

[4]WU Xiaodong, GE Fei, HONG Youshi. A review of recent studies on vortex-induced vibrations of long slender cylinders[J]. Journal of fluids and structures, 2012, 28: 292-308.

[5]TRIM A D, BRAATEN H, LIE H, et al. Experimental investigation of vortex-induced vibration of long marine risers[J]. Journal of fluids and structures, 2005, 21(3): 335-361.

[6]SONG Jining, LU Lin, TENG Bin, et al. Laboratory tests of vortex-induced vibrations of a long flexible riser pipe subjected to uniform flow[J]. Ocean engineering, 2011, 38(11/12): 1308-1322.

[7]HUERA-HUARTE F J, BANGASH Z A, GONZLEZ L M. Towing tank experiments on the vortex-induced vibrations of low mass ratio long flexible cylinders[J]. Journal of fluids and structures, 2014, 48: 81-92.

[8]HANSON A R. Vortex shedding from yawed cylinders[J]. AIAA journal, 1966, 4(4): 738-740.

[9]VAN ATTA C W. Experiments on vortex shedding from yawed circular cylinders[J]. AIAA journal, 1968, 6(5): 931-933.

[10]FRANZINI G R, GON?ALVES R T, MENEGHINI J R, et al. One and two degrees-of-freedom vortex-induced vibration experiments with yawed cylinders[J]. Journal of fluids and structures, 2013, 42: 401-420.

[11]BOURGUET R, KARNIADAKIS G E, TRIANTAFYLLOU M S. On the validity of the independence principle applied to the vortex-induced vibrations of a flexible cylinder inclined at 60°[J]. Journal of fluids and structures, 2015, 53: 58-69.

[12]LIE H, KAASEN K E. Modal analysis of measurements from a large-scale VIV model test of a riser in linearly sheared flow[J]. Journal of fluids and structures, 2006, 22(4): 557-575.

[13]CHAPLIN J R, BEARMAN P W, HUERA HUARTE F J, et al. Laboratory measurements of vortex-induced vibrations of a vertical tension riser in a stepped current[J]. Journal of fluids and structures, 2005, 21(1): 3-24.

Vortex-induced vibration of an inclined flexible cylinder with a small yaw angle

XU Wanhai，MA Yexuan, LUO Hao，LUAN Yingsen，WU Mengning

(State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072，China)

Vortex-induced vibrations (VIV) could cause fatigue damage in flexible cylinder structures for ocean engineering. VIV has been extensively studied through the canonical problem of a circular cylinder free to oscillate within a flow perpendicular to the body axis. In practical applications, the cylinder structures are frequently inclined with respect to the direction of the oncoming flow. Some researchers proposed the independence principle (IP) to simplify this complicated problem,but it remains controversial up till now. In this study, laboratory tests were conducted on VIV of an inclined flexible cylinder with a small yaw angle (15°), an aspect ratio of 350 and a mass ratio of 1.9. The Reynolds number ranged from about 800 to 16 000. The experimental results from a vertical cylinder and an inclined cylinder were compared. It was found that the IP is appropriate for VIV analysis of an inclined flexible cylinder with a small yaw angle outside the modal transformation region. In addition, a high mode is more easily excited by an inclined cylinder in the modal transformation region.

inclined flexible cylinder; vortex-induced vibration (VIV); independence principle (IP); fluid-solid interaction

2015-11-17.

日期：2016-11-14.

國家重點基礎研究計劃 (2014CB046801)； 國家自然科學基金項目(51379144，51479135，51679167)；國家自然科學基金創(chuàng)新研究群體科學基金項目(51621092)；天津市自然科學基金項目(15JCQNJC07700).

徐萬海(1981-)，男，副教授.

徐萬海， E-mail：xuwanhai@tju.edu.cn.

10.11990/jheu.201511035

http://www.cnki.net/kcms/detail/23.1390.u.20161114.1036.021.html

TB123;TV312

A

1006-7043(2017)02-0195-06

徐萬海，馬燁璇, 羅浩, 等. 小傾角傾斜柔性圓柱渦激振動實驗研究[J]. 哈爾濱工程大學學報， 2017, 38(2): 195-200. XU Wanhai，MA Yexuan, LUO Hao， et al. Vortex-induced vibration of an inclined flexible cylinder with a small inclined angle[J]. Journal of Harbin Engineering University, 2017, 38(2): 195-200.