陳 偉, 宗 智
(1.大連理工大學,大連 116023;2.工業(yè)裝備結構分析國家重點實驗室,大連 116023)
立管渦激振動的離散渦數(shù)值模擬
陳 偉1,2, 宗 智1,2
(1.大連理工大學,大連 116023;2.工業(yè)裝備結構分析國家重點實驗室,大連 116023)
通過離散渦方法求解立管所受渦激升力,并與立管動力響應方程相耦合來對渦激振動問題進行時域內(nèi)的數(shù)值模擬。通過在均勻來流和剪切來流下的模擬,所得的結論與其他文獻中的結論比較吻合,較好地揭示了立管渦激振動的特征。
立管;渦激振動;離散渦方法;CFD
Abstract:Discrete vortex method is used to calculate the lifting force induced by the vortex shedding from the structure.The results will couple with the riser dynamic response equation to realize a numerical simulation of VIV in time domain.Some conclusions are got after simulating the VIV phenomenon of riser under uniform flow and shear flow,then they compared with other literatures’conclusions show a good agreement,and the VIV characteristics of riser are revealed better.
Key words:riser;Votex-Induced Vibration;discrete vortex method;CFD
立管是海洋結構物不可或缺但又十分薄弱的子系統(tǒng)。在海洋洋流環(huán)境下,立管結構由于發(fā)生泄渦而產(chǎn)生的渦激振動會縮短結構的疲勞壽命,甚至會引起結構的破壞。渦激振動問題是流體力學中復雜的問題之一,一直困擾著學術、工程界[1]。對立管的渦激振動研究可分為實驗和數(shù)值模擬等兩方面研究。前者主要集中在模型試驗,現(xiàn)場實驗因花費巨大而很少進行;后者因借助于計算能力的飛速發(fā)展而多有進行。但目前全三維的數(shù)值模擬仍不可行。關于渦激振動的諸多計算方法,可見文獻[2]。
本文采用離散渦方法(DVM)求解 N-S方程,DVM方法將速度與壓力表示的控制方程轉(zhuǎn)變?yōu)闇u量場來求解,減小了計算區(qū)域,具有不依賴于計算網(wǎng)格等優(yōu)點,較適合高雷諾數(shù)條件下的流場計算[3]。離散渦方法至今在提高精度和計算效率方面已有很大的改進,在工程領域得以廣泛應用[4]。
立管結構尺度比很大,對其進行三維數(shù)值模擬目前還不現(xiàn)實。本文借助船舶水動力計算的切片思路,沿立管長度方向取若干個水動力剖面,在這些剖面上進行二維流場求解,同時與立管的動力響應方程相互耦合來實現(xiàn)立管結構的渦激振動模擬。
以渦量ω和流函數(shù)ψ表示的二維不可壓縮非定常粘性流動的控制方程為
式中v為流體的運動粘性系數(shù)。離散渦方法將連續(xù)的渦量場離散,把離散的渦元視為粒子,在拉格朗日框架內(nèi)離散求解上述控制方程。
Chorin提出算子分裂法[5]來處理考慮粘性時的(1)式,將其分為對流和粘性擴散2個部分:
對流部分:
對于(3)式的求解即為不考慮粘性時渦元粒子之間的相互誘導,而(4)式的解與隨機運動過程的解相同,可以采用隨機走位法來求解。離散渦方法主要分為物面渦的生成、渦元的對流和擴散、渦元的消除和融合、物體受力計算等幾個方面,其具體列式和步驟可參考文獻[6,7]。
考慮軸向變形時立管的控制方程可寫為[8]
式(5)中:EI是立管的彎曲剛度;T是作用在管壁的軸向張力;w是單位長度立管的重量(包括內(nèi)部液體);f是單位長度的橫向載荷。取y軸垂直向上,x為水平方向,坐標原點在海底。式(6)中:A為立管的實體截面面積;u為軸向位移;q為軸向載荷。采用迦遼金法對上兩式處理可得到立管彈性剛度矩陣 KE和幾何剛度矩陣KG。由此可求得一定來流下立管的靜平衡位置:
式中:d為廣義位移向量;F為載荷向量。本文采用直接迭代法求得立管的靜平衡位置,此處所得到的總剛度矩陣在動力求解過程中保持不變。立管的阻尼采用瑞利阻尼,振型阻尼取5%。動力響應方程為
式中:M、B、K分別為系統(tǒng)的質(zhì)量矩陣、阻尼矩陣和剛度矩陣;F(t)為外載荷矩陣。本文模型采用的是桿梁的組合模型。所選取的立管和海水參數(shù)如表1所示,立管為 TTR型,兩端簡支處理。
表1 立管和海水的相關參數(shù)[9,10]
圖1 不同來流速度下的橫向振動位移
為簡化計算模型,在動力求解中假定結構剛度矩陣保持不變,立管在一定來流下的靜平衡位置附近振動。因此,首先計算得到立管的靜平衡位置。為后文分析方便還需要對立管進行模態(tài)求解。進入動力求解后,在同一時刻,對所有的水動力剖面,即集中質(zhì)量所在位置進行DVM求解,得到結構的外載荷矩陣;然后用Newmark方法求解結構在下一時刻的動力響應,得到各節(jié)點處的運動物理量和新位置后,再利用DVM求外載荷。幾經(jīng)如此反復后,實現(xiàn)流固耦合模擬。
實際環(huán)境下的海流沿垂向分布形式往往較復雜,本文采用FORTRAN編程計算了均勻來流和剪切來流情況下的立管渦激響應。立管沿長度方向劃分為60個長度相同的單元,需要計算的流體剖面有50個,各水動力剖面每次新生120個點渦。采用TECPLOT軟件示意計算結果。
4.1 均勻來流下立管的渦激振動
在均勻來流情況下,所取來流速度大小范圍為0.16 m/s~1.5 m/s,其對應的雷諾數(shù)范圍為4.0×104~3.75×105,此時的斯托哈爾數(shù)在0.21左右。來流速度分別為0.2、0.45、0.76和1.15 m/s時立管的橫向振動位移如圖1所示。圖中粗實線是出現(xiàn)最大節(jié)點位移時的主振型,細實線描述了出現(xiàn)最大節(jié)點位移時的次要振型,而細虛線則表示其他時刻的振型。由圖1可知,隨著來流速度的增加,立管振動的模態(tài)也隨之增加。對于同一立管模型,文獻[10]中取80個單元,計算出各模態(tài)下的最大振幅在0.6左右,本文為0.42,誤差的原因還需進一步分析。
圖2所示在上述四個來流速度下,計算500個時間步長后立管各水動力剖面內(nèi)的尾流渦元分布。對比尾流渦元分布圖及其相應的振動位移響應圖可知,尾流的渦形式沿立管長度方向發(fā)生變化:在振動劇烈的部位,尾渦呈現(xiàn)出2P模式;在振動的“節(jié)點”處,尾流為2S模式。這與 Techet等人進行的試驗[10]結果一致。各來流下的渦模式變化都已在圖中標出。對于上述四個來流速度,取 Z=24 m這一水動力剖面,其位移和升力時程曲線如圖3所示。
根據(jù)折合速度Ur的定義,其倒數(shù)稱為折合頻率,即
在均勻來流情況下,fn取振動的模態(tài)頻率;U為來流速度;D為圓柱直徑。表2列出了各個來流速度下計算的時間步長、主要的響應模態(tài)、對應的折合速度以及折合頻率。
由表2可見,立管的渦激振動的折合速度大都保持在4~7之間,即在這一區(qū)間內(nèi)立管從外界流體獲取能量,對應的折合頻率為0.25~0.14,正好包含斯托哈爾數(shù) St所在的區(qū)間。由于 St表示的是泄渦頻率特征,這說明渦激振動中泄渦頻率有較寬的頻帶,不同來流下鎖定發(fā)生的頻率亦有不同。所有這些都表明立管的振動模態(tài)是由泄渦和鎖定發(fā)生的區(qū)間來控制的。
圖3 不同速度下的振動位移和升力系數(shù)時程曲線,y=24 m
4.2 剪切來流下立管的渦激振動
本文考慮梯形和三角形兩種形式下的剪切來流,最大來流速度取1.2 m/s,最小速度取0.4 m/s。計算所得的立管振動響應如圖4所示??梢?在剪切來流情況下,立管的振動位移較均勻來流情況下要小很多,這是由于沿立管長度方向雷諾數(shù)與泄渦頻率不同所造成的。但主要的振動頻率還是由平均雷諾數(shù)下的泄渦所控制。此時的平均來流為0.8 m/s,參考表2可知,正是第三階頻率占主要地位,這與模擬結果一致。
表2 不同來流速度下立管響應特性
通過對立管在均勻來流和剪切來流下的渦激振動數(shù)值模擬,得到如下結論:
(1)在均勻來流情況下,隨著來流速度的增加,立管振動的模態(tài)會隨之增加。立管最大的位移并不一定發(fā)生在主要模態(tài)上。尾渦的形式沿立管長度方向發(fā)生變化,在振動劇烈的區(qū)域為2P模式,振動的“節(jié)點”區(qū)域為2S模式。(2)剪切來流情況下由于沿立管長度方向泄渦頻率各不一樣,導致立管的振動位移比均勻來流情況下的要小。泄渦頻率不同,導致立管受激發(fā)獲取能量的區(qū)域同樣比均勻來流情況下的要小。(3)立管的振動模態(tài)是由泄渦頻率和鎖定的區(qū)間所控制的。
圖4 剪切流下立管的振動位移
[1] 潘志遠.海洋立管渦激振動機理與預報方法研究[D].上海:上海交通大學,2006.
[2] 秦延龍,王世澎.海洋立管渦激振動計算方法進展[J].中國海洋平臺,2008,23(4):14-17.
[3] 祝寶山.非定常流動的快速拉格朗日渦方法數(shù)值模擬[J].力學學報,2008,40(1):9-18.
[4] Kamemoto K.On contribution of advanced vortex element methods toward virtual reality of unsteady vortical flows in the new generation of CFD[J].ABCM,2004,XXVI(4):368-378.
[5] Chorin A J.Numerical study of slightly viscous flow[J].Journal of Fluid Mechanics,1973,57:785-796.
[6] Mustto A A,Bodstein G C R.Improved vortex method for the simulation of the flow around circular cylinders[C].AIAA,2001,NO:A01-31110,pp:1-11.
[7] 潘巖松.高層建筑二維流場的離散渦方法數(shù)值模擬[D].武漢:華中科技大學,2004.
[8] Patel M H,Witz J A.Complaint Offshore structures[M].Butterworth-Heinemann,England,1991.
[9] Silvio B C Martins,Sergio H Sphaier,Severino F Silva Neto.Integrated fluid-structure model(RVM-FEM)for riser analysis[C].International Offshore and Polar Engineering Conference,2004,2:105-112.
[10] Yamamoto C T,Meneghini J R,Saltara F,Fregonesi R A,Ferrari J A Jr.Numerical simulation of vortex-induced vibration on flexible cylinders[J].Journal of Fluids and Structures,2004,19:467-489.
Numerical Simulation of Vortex-Induced Vibration for Risers by Discrete Vortex Method
CHEN Wei1,2, ZONG Zhi1,2
(1.Dalian University of Technology,Dalian 116023,China;2.State Key Laboratory of Structural Analysis for Industrial Equipment,Dalian 116023,China)
O327,P751
A
1001-4500(2010)01-0026-06
2009-08-19; 修改稿收到日期:2009-10-30
陳 偉(1985-),男,碩士研究生,從事船舶與海洋結構物設計與制造工作。