深海立管渦激振動預(yù)報模型及影響因素

唐友剛, 青兆熹, 張杰, 王賓

(1. 天津大學(xué) 水利工程仿真與安全國家重點實驗室，天津 300072； 2.上海海事大學(xué) 海洋科學(xué)與工程學(xué)院，上海 201306)

深海立管渦激振動預(yù)報模型及影響因素

唐友剛, 青兆熹, 張杰, 王賓

(1. 天津大學(xué) 水利工程仿真與安全國家重點實驗室，天津 300072； 2.上海海事大學(xué) 海洋科學(xué)與工程學(xué)院，上海 201306)

渦激振動是立管發(fā)生破壞的主要原因之一，深海立管自重大、柔度高、頂部張力集中等會導(dǎo)致出現(xiàn)新的動力特性。為預(yù)報深海立管渦激振動并揭示其動力特性，考慮立管自重影響及線性剪切流，本文建立立管渦激振動方程，基于Van der Pol尾流振子模型，采用有限差分法計算立管的振動響應(yīng)，并設(shè)計渦激振動試驗進(jìn)行驗證，最后研究流速及頂張力對立管渦激振動的影響。結(jié)果表明：流速越大，立管渦激振動頻率越高，振動應(yīng)力越大；同等流速下，頂張力越大，立管渦激振動主頻率變化不大，但振動位移增大，振動應(yīng)力減小。

深海立管；渦激振動；預(yù)報模型；試驗；振動響應(yīng)；尾流振子模型

海流經(jīng)過立管時，在立管下游產(chǎn)生尾流和漩渦，周期發(fā)放的漩渦對立管產(chǎn)生垂直于流向的渦激升力，引起立管的渦激振動，渦激振動是立管發(fā)生破壞的關(guān)鍵因素之一[1]。由于海流速度沿水深變化的非均勻性及流固耦合的復(fù)雜性，渦激振動的準(zhǔn)確預(yù)報是一個巨大的難題[2]。

目前，研究渦激振動的方法主要有兩種，即模型試驗法和數(shù)值模擬法[3]。數(shù)值模擬法主要有3類方法：1)基于切片理論的CFD方法，這類模型大多將二維Navier-Stokes方程的求解與梁模型結(jié)合在一起，通過求解方程得到作用在圓柱上的力，然后將力反饋到圓柱上求得結(jié)構(gòu)的動力響應(yīng)。2)基于圓柱受迫振動實驗數(shù)據(jù)的經(jīng)驗?zāi)Ｐ?，常用的有DNV模型、LIC模型、MARINTE模型等。這些模型的基本假定是渦激振動發(fā)生在一個或有限個離散頻率上，基于這些經(jīng)驗?zāi)Ｐ瓦M(jìn)行渦激振動預(yù)報的軟件有SHEAR7、VIVA、VIVANA等。3)介于上述兩者之間的尾流振子模型，也稱半經(jīng)驗?zāi)Ｐ?。通過相互獨(dú)立的方式分別建立圓柱振子運(yùn)動方程和流體振子運(yùn)動方程，然后利用它們共同預(yù)報流體-彈性系統(tǒng)的動力響應(yīng)[4-6]。

立管下端一般通過萬向節(jié)與海底井口相連，上端通過升沉補(bǔ)償裝置(又稱張緊器)與平臺相連。張緊器給立管提供較大的頂張力，支持立管的重量，使立管保持張緊垂直狀態(tài)，避免立管長度過大致使底部發(fā)生屈曲。對于淺海立管(水深一般小于300 m)，由于其剛度較大，自重對立管軸向張力分布影響不大，多忽略自重影響做恒張力處理[7]；淺海海流流速沿水深變化不如深海明顯，柱體受波流聯(lián)合作用明顯，一般渦激振動分析時可作均勻流處理[8]；此外淺海立管固有頻率主要受自身剛度控制，渦激振動“鎖頻”現(xiàn)象明顯，最大響應(yīng)多出現(xiàn)在立管中上部[9]。

對于深海立管，長徑比迅速增大，立管固有頻率主要受張力控制，而且深海海底流速低，渦泄頻率與立管低階固有頻率接近，此外立管自重大、柔度高、頂部張力集中等因素也增加了立管渦激振動問題的復(fù)雜性。因此對深海立管的渦激振動現(xiàn)象亟需更深入的研究。

為預(yù)報深海立管渦激振動并揭示其動力特性，本文考慮立管自重影響及線性剪切流，建立立管渦激振動方程，基于Van der Pol尾流振子模型，采用有限差分法計算立管的振動響應(yīng)，并通過與模型試驗對比驗證該數(shù)值模型的可靠性，最后研究流速及頂張力對立管渦激振動的影響。

1 立管渦激振動分析模型

針對長細(xì)比很大的深海立管，考慮立管自重引起的軸向變張力和剪切流引起的渦激升力，簡化后的立管模型示意如圖1所示。

基于歐拉梁的彎曲振動理論，建立深海立管渦激振動方程：

(1)

圖1 立管模型示意圖Fig.1 Simply supported riser

fx(z,t)是單位長度外激勵，由漩渦泄放產(chǎn)生的渦激升力fL(z,t)和立管振動產(chǎn)生的流體阻尼力ft(z,t)組成：

fy(z,t)=fL(z,t)-ft(z,t)

(2)

立管兩端邊界條件視為鉸接：

(3)

考慮立管與流體之間的流-固耦合影響，基于VanderPol方程描述漩渦的尾流特性[10]：

(4)

2.1 兩組紅細(xì)胞計數(shù)、血紅蛋白、血清鐵蛋白比較 兩組干預(yù)后紅細(xì)胞計數(shù)、血紅蛋白、血清鐵蛋白較干預(yù)前顯著提高(P<0.05)。干預(yù)前組間紅細(xì)胞計數(shù)、血紅蛋白、血清鐵蛋白比較，差異無統(tǒng)計學(xué)意義(P>0.05)。觀察組干預(yù)后紅細(xì)胞計數(shù)、血紅蛋白、血清鐵蛋白顯著高于對照組(P<0.05)。見表1。

聯(lián)立方程(1)、(4)，深海立管渦激振動預(yù)報模型為

(5)

2 數(shù)值算法

由式(5)可知，深海立管渦激振動預(yù)報模型體現(xiàn)了流固耦合的特性，拖曳力阻尼是非線性的，所以模態(tài)疊加法不再適用，本文應(yīng)用有限差分法對方程(5)在時間和空間上離散，以求解立管的振動響應(yīng)。

應(yīng)用中心差分格式，將立管長度平均分為N份，得到N個計算單元和N+1個節(jié)點，單元長度為ΔN=L/N，立管頂端的節(jié)點編號為0，第二個為1，依次分布直到N。根據(jù)邊界條件(3)，頂端節(jié)點0和底端節(jié)點N的運(yùn)動形式是已知的，于是得到關(guān)于時間t的2(N-1)個二階常微分方程組，應(yīng)用四階Runge-Kutta法數(shù)值求解[3]。

3 渦激振動試驗

為驗證深海立管渦激振動數(shù)值計算的可靠性，在天津大學(xué)水利工程仿真與安全國家重點實驗室船舶拖曳水池中(137 m×7 m×3 m)進(jìn)行了深海立管渦激振動試驗。

受水池深度限制，立管試驗?zāi)Ｐ蜋M向布置，由拖車勻速帶動模擬均勻流。采用應(yīng)變片監(jiān)測立管應(yīng)力變化。試驗裝置設(shè)備主要包括立管模型(Teflon管)、萬向節(jié)、鋼絲繩、彈簧、滑輪、測速儀和張力計等。試驗測量設(shè)備包括電阻應(yīng)變片、電阻應(yīng)變儀、信號采集儀和振動信號采集分析軟件等。立管實體及模型參數(shù)如表1所示。

表1 立管模型參數(shù)

4 試驗與計算結(jié)果對比

理論計算得到振動應(yīng)力、振動頻率和振動位移與試驗中測量的數(shù)據(jù)對比如圖2～4所示。

經(jīng)過大量的試驗對比驗證，說明本文采用的渦激振動預(yù)報模型是合理的，Matlab編程計算的結(jié)果是可靠的。

圖2 中點應(yīng)力時程曲線Fig.2 Bending stress time-history of riser at midpoint

圖3 振動頻譜分析(U=0.5 m/s, Tt=100 N)Fig.3 Frequency spectrum of riser (U=0.5 m/s, Tt=100 N)

圖4 振動最大位移[11]Fig.4 The maximum displacement of riser[11]

5 影響因素分析

不考慮波浪的影響，海面流速一般為0.1～1 m/s[12]。當(dāng)頂張力系數(shù)ftop=1.3時，計算不同流速下立管的振動響應(yīng)，根據(jù)立管每一點的位移時程曲線和彎曲應(yīng)力時程曲線，得到立管振動的最大位移和最大彎曲應(yīng)力，圖5、6是選取高、中、低三種流速：0.2、0.5和0.8m/s計算得到的最大位移和最大彎曲應(yīng)力對比圖。

由圖5看出，流速越大，立管渦激振動激發(fā)模態(tài)越高，但最大振動位移略有減小，說明立管渦激振動位移主要受低階模態(tài)控制；低流速時振動位移最大值在立管底部，高流速時振動位移最大值在立管中上部。

由圖6看出，流速越大，立管渦激振動彎曲應(yīng)力越大，說明立管渦激振動彎曲應(yīng)力主要受高階模態(tài)控制；不管何種流速，彎曲應(yīng)力最大值始終在立管底部。

不同流速下立管振動最大位移、最大彎曲應(yīng)力及振動主頻率的統(tǒng)計如表2所示?？梢钥闯觯?/p>

1)振動位移：流速越大，最大位移越小，最大值出現(xiàn)的位置越靠近立管上部。原因分析：小流速產(chǎn)生的渦激升力小，主要激起低階模態(tài)，低階模態(tài)產(chǎn)生的位移較大，特別是立管底部由于內(nèi)部張力較小，因此有較大的位移；大流速激起高階模態(tài)，高階模態(tài)產(chǎn)生的位移較小，立管中上部由于渦激升力大，導(dǎo)致有較大的振動位移。

圖5 立管渦激振動位移最大值Fig.5 The maximum displacement of riser in VIV

圖6 立管渦激振動彎曲應(yīng)力最大值Fig.6 The maximum bending stress of riser in VIV

2)振動彎曲應(yīng)力：流速越大，激發(fā)模態(tài)越高，彎曲應(yīng)力越大，最大彎曲應(yīng)力均集中在立管底部。原因分析：彎曲應(yīng)力大小由振動位移的二次導(dǎo)數(shù)決定，振動位移二次導(dǎo)數(shù)隨模態(tài)及水深的變化趨勢同立管模態(tài)曲率隨模態(tài)和水深的變化趨勢是一致的，即深海立管渦激振動位移的二次導(dǎo)數(shù)最大幅值在立管底部，在振動位移相差不大的情況下，位移二次導(dǎo)數(shù)的最大幅值隨模態(tài)階次升高而增大[3]。

表2 不同流速下立管渦激振動響應(yīng)統(tǒng)計表

5.2 頂張力對立管渦激振動的影響

當(dāng)U0=0.5 m/s時，計算頂張力系數(shù)ftop分別為1.1、1.3和1.5時立管振動的最大位移和最大彎曲應(yīng)力，如圖7、8所示。

可以看出，頂張力增大，立管振動位移略有增大，但立管彎曲應(yīng)力明顯減小。這是因為頂張力系數(shù)增大，則立管內(nèi)部張力增大，相當(dāng)于增大了立管的彎曲剛度，立管模態(tài)固有頻率升高，在同等流速下由于渦激頻率不變，導(dǎo)致立管激起模態(tài)降低，從而使立管振動位移增大，彎曲應(yīng)力減小。這是深海立管所特有的性質(zhì)，與普通梁或淺海立管橫向振動隨軸向力變化的規(guī)律是不一致的。

立管振動響應(yīng)頻譜分析如圖9所示。由圖9可以發(fā)現(xiàn)，改變頂張力系數(shù)，立管渦激振動主頻率變化不大。說明立管振動主頻率主要受渦激升力頻率的控制，也就是受流速的控制。

圖7 立管渦激振動位移最大值Fig.7 The maximum displacement of riser in VIV

圖8 立管渦激振動彎曲應(yīng)力最大值Fig.8 The maximum bending stress of riser in VIV

圖9 立管渦激振動響應(yīng)頻譜對比圖Fig.9 Frequency spectrum of riser in VIV

6 結(jié)論

1)流速越大，渦泄頻率越高，立管渦激振動激發(fā)模態(tài)越高，立管振動主頻率越大，振動應(yīng)力也越大；

2)頂張力增大，立管固有頻率升高，在同等流速下雖然立管振動響應(yīng)主頻率變化不大，但導(dǎo)致立管激起振動的模態(tài)降低，從而使立管振動應(yīng)力減小。

[1]XUE Hongxiang, WANG Kunpeng, TANG Wenyong. A practical approach to predicting cross-flow and in-line VIV response for deepwater risers[J]. Applied ocean research, 2015, 52: 92-101.

[2]SRINIL N. Analysis and prediction of vortex-induced vibrations of variable-tension vertical risers in linearly sheared currents[J]. Applied ocean research, 2011, 33(1): 41-53.

[3]張杰. 深海立管參激-渦激聯(lián)合振動與疲勞特性研究[D]. 天津: 天津大學(xué), 2014. ZHANG Jie. Vortex-induced vibration and fatigue analysis of deepwater risers considering parametric excitations[D]. Tianjin: Tianjin University, 2014.

[4]魏東澤, 白興蘭, 顧恩凱. 考慮流固耦合的立管渦激振動研究[J]. 中國造船, 2014, 55(3): 108-116. WEI Dongze, BAI Xinglan, GU Enkai. Study on vortex-induced vibration of riser considering fluid-structure interaction[J]. Shipbuilding of China, 2014, 55(3): 108-116.

[5]唐世振. 考慮順流向振動的深水頂張力立管渦激振動分析[D]. 青島: 中國海洋大學(xué), 2010. TANG Shizhen. Study of the vortex induced vibration for the deepwater top tensioned risers considering the in-line vibration[D]. Qingdao: Ocean University of China, 2010.

[6]YAMAMOTO C T, MENEGHINI J R, SALTARA F, et al.

Numerical simulations of vortex-induced vibration on flexible cylinders[J]. Journal of fluids and structures, 2004, 19(4): 467-489.

[7]張杰, 唐友剛. 深海立管固有振動特性的進(jìn)一步分析[J]. 船舶力學(xué), 2014, 18(1/2): 165-171. ZHANG Jie, TANG Yougang. Further analysis on natural vibration of deep-water risers[J]. Journal of ship mechanics, 2014, 18(1/2): 165-171.

[8]董艷秋. 波、流聯(lián)合作用下海洋平臺張力腿的渦激非線性振動[J]. 海洋學(xué)報, 1994, 16(3): 121-129. DONG Yanqiu. Nonlinear vibration analysis of tension leg under combined wave-current[J]. Acta oceanologica sinica, 1994, 16(3): 121-129.

[9]郭海燕, 傅強(qiáng), 婁敏. 海洋輸液立管渦激振動響應(yīng)及其疲勞壽命研究[J]. 工程力學(xué), 2005, 22(4): 220-224. GUO Haiyan, FU Qiang, LOU Min. Vortex-induced vibrations and fatigue life of marine risers conveying flowing fluid[J]. Engineering mechanics, 2005, 22(4): 220-224.

[10]FACCHINETTI M L, DE LANGRE E, BIOLLEY F. Coupling of structure and wake oscillators in vortex-induced vibrations[J]. Journal of fluids and structures, 2004, 19(2): 123-140.

[11]FUJIWARA T, UTO S, KANADA S. An experimental study of the effects that change the vibration mode of riser VIV[C]//ASME 2011 30th International Conference on Ocean, Offshore and Arctic Engineering. Rotterdam, The Netherlands: ASME, 2011: 487-492.

[12]GAO Yun, ZONG Zhi, SUN Lei. Numerical prediction of fatigue damage in steel catenary riser due to vortex-induced vibration[J]. Journal of hydrodynamics, ser. B, 2011, 23(2): 154-163.

Prediction model and influence factors on vortex-induced vibration of deepwater risers

TANG Yougang1, QING Zhaoxi1, ZHANG Jie1,2, WANG Bin1

(1. National Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072, China; 2. College of Ocean Science and Engineering, Shanghai Maritime University, Shanghai 201306, China)

Vortex-induced vibration (VIV) is one of the main causes of the destruction of marine risers. The dynamic VIV characteristics of deepwater risers may include their large self-weight, high flexibility, and the stress concentration caused by top tension, among others. In order to predict the VIV of deepwater risers and determine their dynamic characteristics, we established a governing equation of VIV in risers that takes into account the self-weight and linear shear flow. We used the finite difference method to calculate the dynamic response of a riser, based on Van der Pol wake-oscillator model. We then conducted a VIV experiment to verify the prediction model. In this study, we discuss the effects of flow velocity and top tension on the VIV of deepwater risers. Our study results show that the vibration frequency and vibration stress of a riser increases, as the flow velocity increases. Under the same flow conditions, the vibration displacement of a riser increases and the vibration stress decreases as the top tension increases. The dominant frequency of a riser is hardly influenced at all by the top tension.

deepwater risers; vortex-induced vibration; prediction model; experiment; vibration response； wake-oscillator model

2016-03-28.

日期：2017-01-11.

國家自然科學(xué)基金項目(51239008)；天津大學(xué)水利工程仿真與安全國家重點實驗室基金項目(HESS-1510).

唐友剛(1952-), 男, 教授, 博士生導(dǎo)師； 張杰(1982-), 男,講師，博士.

張杰，E-mail: zhangjie@shmtu.edu.cn.

10.11990/jheu.201603096

TE58;TU311.3

A

1006-7043(2017)03-0338-06

唐友剛, 青兆熹, 張杰，等.深海立管渦激振動預(yù)報模型及影響因素[J]. 哈爾濱工程大學(xué)學(xué)報, 2017, 38(3):338-343.

TANG Yougang, QING Zhaoxi, ZHANG Jie，et al.Prediction model and influence factors on vortex-induced vibration of deepwater risers[J]. Journal of Harbin Engineering University, 2017, 38(3):338-343.

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/23.1390.u.20170111.1509.008.html