尾流干涉下立管渦激振動試驗(yàn)研究*

趙德光 婁敏 王曉凱 張睿

(1.中國石油大學(xué)(華東)石油工程學(xué)院 2.海洋石油工程(青島)有限公司)

0 引 言

渦激振動是一種廣泛存在于自然界中的流固耦合現(xiàn)象，海水流經(jīng)立管時會在尾流區(qū)內(nèi)產(chǎn)生交替的漩渦，當(dāng)漩渦脫落頻率接近立管自振頻率時會發(fā)生“鎖定”現(xiàn)象，從而導(dǎo)致結(jié)構(gòu)發(fā)生渦激振動現(xiàn)象[1-2]，進(jìn)而縮短立管壽命[3]。

目前，國內(nèi)外學(xué)者針對渦激振動問題進(jìn)行了廣泛研究。在數(shù)值求解方面，CUI Z.D.等[4]采用有限元法對并聯(lián)排列的彈性耦合圓柱體渦激振動開展研究，發(fā)現(xiàn)兩個圓柱的主頻率在“鎖定”狀態(tài)內(nèi)、外始終保持一致。M.RAHMANIAN等[5]對兩個不同直徑并聯(lián)排列圓柱體的渦激振動進(jìn)行研究，兩圓柱僅考慮橫流向振動，結(jié)果表明，小圓柱體的存在對大圓柱體的漩渦脫落方式和振動幅值產(chǎn)生了較大影響。王曉凱等[6]對串聯(lián)排列圓柱的干涉效應(yīng)進(jìn)行研究，發(fā)現(xiàn)間距比和約化速度對下游圓柱的振動響應(yīng)有顯著影響。R.C.MYSA等[7]在低雷諾數(shù)下對串聯(lián)排列圓柱的流動誘導(dǎo)振動進(jìn)行研究，分析了升力組成中壓力和黏滯力的意義以及與流速的關(guān)系。FAN X.T.等[8]研究了不同間距比下，串聯(lián)排列三圓柱在亞臨界雷諾數(shù)下的渦激振動模式，發(fā)現(xiàn)位移峰值是渦激振動與尾流激振共同作用的結(jié)果。武磊等[9]基于Open FOAM求解器，對階梯流中不同浸沒長度下的串聯(lián)排列雙立管渦激振動展開研究，發(fā)現(xiàn)隨著浸沒長度的增加，下游立管的漩渦泄放受上游立管的尾流作用得到進(jìn)一步加強(qiáng)。在試驗(yàn)研究方面，G.R.S.ASSI[10-11]對兩自由度串聯(lián)、交錯排列的圓柱體尾流誘導(dǎo)振動進(jìn)行了試驗(yàn)研究，結(jié)果表明，當(dāng)下游圓柱的初始位置從尾跡中心線向兩側(cè)偏離時，由誘導(dǎo)振動引起的激振效果會減弱。S.W.KIM等[12]在亞臨界雷諾數(shù)范圍內(nèi)，研究了柔性圓柱雙自由度的耦合運(yùn)動，采用三維位移計(jì)和張力計(jì)測量模型的動態(tài)響應(yīng)，并將得到的響應(yīng)幅值和漩渦脫落頻率與他人的試驗(yàn)結(jié)果進(jìn)行對比分析。HU Z.M.等[13]利用風(fēng)洞試驗(yàn)研究了不同質(zhì)量-阻尼比下單自由度串聯(lián)排列圓柱的流動誘導(dǎo)振動特性，結(jié)果表明，隨著質(zhì)量-阻尼比的增大，兩圓柱的無量綱振幅、振動區(qū)域與滯后現(xiàn)象等所受影響將逐漸減小。梁文洲[14]系統(tǒng)地對串列排列、并列排列和交錯排列的雙圓柱體開展試驗(yàn)研究，得到了有別于單圓柱體和固定不動雙圓柱體或雙圓柱體單自由度渦激振動的變化規(guī)律。在渦激振動抑制方面，婁敏等[15]采用有限元軟件分析了不同約化速度下，仙人掌形截面對柔性圓柱體的抑振效果。劉志慧等[16]對具有盤球附體的立管模型進(jìn)行抑渦試驗(yàn)，發(fā)現(xiàn)螺距是影響盤球附體抑渦效果的重要因素。婁敏等[17]對帶三控制桿的串聯(lián)立管進(jìn)行試驗(yàn)，分析了不同立管間距對三控制桿渦激振動抑制效果的影響。宋振華等[18]研究了3根附屬桿抑制裝置對立管的抑制效果及抑制機(jī)理，發(fā)現(xiàn)間隙比和來流角度對立管的抑制效果有顯著影響。

目前在求解立管的渦激振動問題上，二維數(shù)值模擬采用的質(zhì)量-彈簧-阻尼系統(tǒng)忽略了立管非線性運(yùn)動響應(yīng)，難以適用于剛度較小且構(gòu)型復(fù)雜的柔性立管[19]。受限于計(jì)算成本，采用有限元求解三維頂張力立管和懸鏈線立管的渦激振動問題尚不成熟。在此背景下，本文在室內(nèi)波流聯(lián)合水槽中分別開展剛性立管和柔性立管渦激振動試驗(yàn)，選取不同材質(zhì)模型用于模擬實(shí)際工程中的剛性立管與柔性立管(頂張力立管和懸鏈線立管)，并通過試驗(yàn)得到立管模型在不同排列方式下的響應(yīng)特性隨來流速度的變化關(guān)系，對比分析不同工況下流動干涉效應(yīng)，以期為進(jìn)一步研究尾流干涉下海洋立管渦激振動提供參考。

1 試驗(yàn)?zāi)康募胺椒?/h2>

為深入探究剛性立管和柔性立管在不同排列方式下的響應(yīng)特性隨來流速度的變化關(guān)系，以及上游立管對下游立管尾流干涉的影響規(guī)律，本文于2021年3月采用模型試驗(yàn)的方法對立管進(jìn)行動力響應(yīng)分析。

1.1 試驗(yàn)器材

試驗(yàn)水槽為中國石油大學(xué)(華東)室內(nèi)波流聯(lián)合水槽，水槽尺寸為16.0 m(長)×0.8 m(寬)×1.4 m(高)。該水槽能夠造出指定波高和周期的規(guī)則波與不規(guī)則波，同時可以雙向造流，可產(chǎn)生0.05～0.70 m/s的來流速度。

試驗(yàn)用管材料分別為硅膠管和有機(jī)玻璃管(Polymethyl methacrylate，PMMA)。硅膠管質(zhì)地較為柔軟，在室內(nèi)水槽的小尺度條件下能夠隨水流自由運(yùn)動，呈現(xiàn)應(yīng)有的構(gòu)型，主要用于柔性立管的試驗(yàn)；PMMA管又叫亞克力管，綜合力學(xué)性能較好、機(jī)械強(qiáng)度較高，可用于剛性圓柱的模型試驗(yàn)。本次試驗(yàn)中選取的立管模型有效長度(1.5 m)、外徑(10 mm)和壁厚(2 mm)等參數(shù)均保持一致，其他物理特性參數(shù)存在一定差異，具體如表1所示。

表1 材質(zhì)參數(shù)Table 1 Material parameters

根據(jù)水槽尺寸定制鋁合金支架、插槽、插塊結(jié)構(gòu)，以實(shí)現(xiàn)水槽在無需放水的情況下快速更換立管模型，提高試驗(yàn)效率。由于鋁合金支架的剛度整體較大，在來流作用下支架與立管模型不會發(fā)生共振作用。為與實(shí)際工程中立管的運(yùn)動狀況相符，本次試驗(yàn)按照立管模型參數(shù)定做了柔性接頭，模型兩端采用萬向節(jié)進(jìn)行連接，保證其在任意方向上的運(yùn)動不受阻礙。

為捕獲立管的振動響應(yīng)，采用非接觸位移測量儀(3D-PDT精密位移測量儀)和光纖光柵測試系統(tǒng)兩種不同方式進(jìn)行測量。其中，非接觸位移測量憑借兩個分辨率為500萬像素的高清攝像頭對物體表面的標(biāo)記測點(diǎn)進(jìn)行實(shí)時捕捉(不會對流場造成影響)，可在立管模型上自由布置測點(diǎn)，方便對任意位置處的響應(yīng)進(jìn)行監(jiān)測。光纖光柵測試系統(tǒng)由布拉格光纖光柵傳感器(Fiber Bragg Grating sensors，F(xiàn)BG)和光纖光柵解調(diào)儀兩部分組成。FBG體積微小，相比于其他測試儀器對外部流場和結(jié)構(gòu)本身都不會產(chǎn)生明顯影響，當(dāng)傳感器數(shù)目較多時還可以采用焊接的方式將光纖串聯(lián)在一起，有效避免傳統(tǒng)電阻應(yīng)變片單個單線導(dǎo)致的線路冗雜。

1.2 試驗(yàn)步驟

試驗(yàn)開始前，使用白色有機(jī)硅密封膠對立管模型的下端進(jìn)行密封處理，目的是能夠向管內(nèi)注入流體以增大附加質(zhì)量，適當(dāng)降低自振頻率。使用PMMA管進(jìn)行剛性立管的模型試驗(yàn)，對立管模型布置光柵應(yīng)變測點(diǎn)和非接觸測點(diǎn)，每種測量儀器所對應(yīng)的測點(diǎn)個數(shù)和測點(diǎn)位置如圖1所示。使用硅膠管進(jìn)行柔性立管的模型試驗(yàn)，頂張力立管模型與懸鏈線立管模型的構(gòu)型如圖2所示。與剛性立管模型的非接觸測點(diǎn)布置情況一致，從模型的中部開始布置非接觸測點(diǎn)，從上到下依次為測點(diǎn)1、測點(diǎn)2和測點(diǎn)3。將立管模型安裝在鋁合金支架后垂直放置于水槽中，浸沒深度為750 mm。

圖1 立管模型測點(diǎn)布置示意圖Fig.1 Layout of measuring points on riser model

圖2 立管模型Fig.2 Riser model

試驗(yàn)期間對水槽造流機(jī)產(chǎn)生的實(shí)際流速進(jìn)行測試校準(zhǔn)。在波流水槽內(nèi)安裝聲學(xué)多普勒流速儀，對試驗(yàn)過程中的來流流速進(jìn)行實(shí)時監(jiān)測。為保證試驗(yàn)結(jié)果的準(zhǔn)確性，在后續(xù)試驗(yàn)中所采用的流速皆以校核后的實(shí)際流速為準(zhǔn)。

將材質(zhì)為PMMA的剛性立管模型內(nèi)部注滿流體后放置于水槽中，通過敲擊試驗(yàn)測定立管模型的自振頻率，通過非接觸攝像機(jī)測得的自振頻率fn1=2.52 Hz，通過光纖光柵系統(tǒng)測得的自振頻率fn2=2.44 Hz，兩者在數(shù)值上較為接近，確保了測試結(jié)果的準(zhǔn)確性及可靠性。考慮柔性立管的特殊構(gòu)型及運(yùn)動情況，采用非接觸位移測量儀實(shí)時記錄立管模型的振動特性，通過敲擊試驗(yàn)測量柔性立管的自振頻率為1.85 Hz。

2 試驗(yàn)內(nèi)容

為探究剛性立管和柔性立管尾流干涉下的渦激振動特性，本文分別進(jìn)行剛性孤立立管試驗(yàn)、剛性立管尾流干涉試驗(yàn)、柔性孤立立管試驗(yàn)和柔性串聯(lián)立管尾流干涉試驗(yàn)。

在剛性立管模型尾流干涉試驗(yàn)中，兩根立管保持固定間距比，即S/D=5，并采用串聯(lián)排列、交錯排列與并聯(lián)排列3種排列方式。圖3為不同排列方式下的雙立管模型。圖3中α表示來流方向與兩立管軸向中心連線之間的夾角，即來流角度的大小取決于兩立管模型的相對位置，3種不同的排列方式對應(yīng)來角度分別為0°、45°和90°。本文重點(diǎn)關(guān)注來流角度的改變對兩立管模型之間流動干涉效應(yīng)的影響。圖3中D為立管直徑，z軸垂于xy平面，表示立管模型上某點(diǎn)的高度坐標(biāo)，u為來流流速。

圖3 雙立管排布示意圖Fig.3 Schematic diagram of double riser arrangement

在柔性立管模型尾流干涉試驗(yàn)中，兩根立管保持間距比S/D=5并采用串聯(lián)排列方式，其中頂張力立管處于上游位置，懸鏈線立管模型布置在下游位置，測點(diǎn)設(shè)置與前文保持一致。

3 試驗(yàn)結(jié)果及分析

3.1 剛性孤立立管試驗(yàn)分析

當(dāng)漩渦脫落頻率(fs)與固有頻率(fn)的比值趨于相近時立管模型便處于“鎖定”狀態(tài)。圖4表示不同約化速度(Vr=u/(fnD))下立管模型的頻率比。由圖4可得，非接觸攝像機(jī)與光纖光柵系統(tǒng)的測試結(jié)果非常接近，頻率比隨約化速度的增大在數(shù)值上相應(yīng)增加。根據(jù)定義可判斷立管模型的“鎖定”區(qū)間為Vr=5～9。

圖4 頻率比Fig.4 Frequency ratio

重點(diǎn)關(guān)注立管模型中間位置,即非接觸測點(diǎn)1處響應(yīng)特性，立管模型在兩個流向上的幅值(以y/D或x/D來表示)如圖5所示。橫流向振幅(y/D)在“鎖定”區(qū)間內(nèi)顯著增加，當(dāng)Vr=6時取得最大值，約為0.46。離開“鎖定”區(qū)間后幅值迅速下降。順流向幅值隨來流速度的增加呈上升趨勢。觀察發(fā)現(xiàn)，兩種不同采集儀器得到的結(jié)果較為接近，再次印證了試驗(yàn)結(jié)果的準(zhǔn)確性。

圖5 立管模型振動幅值曲線Fig.5 Vibration amplitude curve of riser model

圖6為處理后得到的立管模型位移(以相對位置來表示，即z/L)包絡(luò)圖。圖6中不同顏色的曲線代表立管模型在來流作用時間內(nèi)不同時刻所處的振動位置，L表示立管模型的長度。

圖6 立管模型位移包絡(luò)圖Fig.6 Displacement envelope diagram of riser model

觀察發(fā)現(xiàn)振動幅值沿管長方向呈現(xiàn)出兩端小、中間大的分布特征，最大值出現(xiàn)在立管模型中部位置，同時在“鎖定”區(qū)間內(nèi)運(yùn)動幅值明顯增加，與圖5中呈現(xiàn)的結(jié)果相符。此外，隨著來流速度的不斷增大，立管模型承受的激勵作用不斷提升，在圖6f中可以觀察到2階振型的出現(xiàn)，標(biāo)志著振型由此發(fā)生轉(zhuǎn)變。

3.2 剛性立管尾流干涉試驗(yàn)分析

不同來流角度時下游立管受干涉作用后的響應(yīng)特性變化情況如圖7與圖8所示。圖7表示頻率比的變化情況。受上游立管模型流動干涉效應(yīng)的影響，下游立管模型的頻率比與孤立立管相比發(fā)生了顯著變化。觀察發(fā)現(xiàn)，當(dāng)α=0°時下游立管模型的“鎖定”區(qū)間為Vr=6～10，與孤立立管Vr=5～9相比存在一定滯后，原因在于上游立管的流動干涉效應(yīng)導(dǎo)致流經(jīng)下游立管的實(shí)際來流速度小于設(shè)定流速，因此“鎖定”區(qū)間將向后推移發(fā)生在較大Vr處，同時頻率比在數(shù)值上將小于孤立立管；而當(dāng)α=45°和90°時，下游立管模型的頻率比在數(shù)值上與α=0°相比逐漸增大，同時“鎖定”區(qū)間的范圍與孤立立管保持一致，由此判定α=0°即串聯(lián)排列時下游立管所受流動干涉效應(yīng)最為顯著，隨著來流角度的增大干涉作用逐步減弱。

圖7 不同來流角度下頻率比Fig.7 Frequency ratio in different inflow angles

圖8表示不同來流角度時下游立管模型無量綱振幅隨Vr的變化關(guān)系。由圖8可得，當(dāng)α=0°時，高約化速度時下游立管模型的橫流向振幅(y/D)相比孤立立管明顯提升，原因在于上游立管模型的流動干涉效應(yīng)引發(fā)了流體與下游立管模型間更劇烈的非線性相互作用，致使其振動幅值明顯大于孤立立管，隨來流角度的增大橫流向振幅逐漸降低，愈發(fā)接近于孤立立管；受上游立管模型流動干涉的影響，α=0°時下游立管模型的順流向振幅(x/D)小于孤立立管，隨著來流角度的增大干涉效應(yīng)逐漸減弱，順流向振幅相應(yīng)增加并在數(shù)值上逐漸超過孤立立管。

圖8 不同來流角度下的振幅曲線Fig.8 Amplitude curve in different inflow angles

3.3 柔性孤立立管試驗(yàn)研究

圖9表示柔性孤立立管模型在不同Vr下的橫流向振動幅值。觀察發(fā)現(xiàn)：柔性頂張力立管模型的響應(yīng)幅值變化與剛性立管較為相似，開始階段振動幅值較小、增長緩慢，進(jìn)入“鎖定”區(qū)間后迅速增大并趨于平穩(wěn)，離開“鎖定”區(qū)間后振幅隨之下降，同時立管中部的幅值最大；懸鏈線立管模型則有所不同，主要原因在于其獨(dú)特構(gòu)型導(dǎo)致多階振動模態(tài)的同時存在，使得在同一來流速度下不同測點(diǎn)處的振動幅值差異顯著。

圖9 柔性孤立立管模型振動幅值Fig.9 Vibration amplitude of flexible riser model

3.4 柔性串聯(lián)立管尾流干涉試驗(yàn)研究

在柔性立管尾流干涉試驗(yàn)中，下游懸鏈線立管模型受流動干涉作用后的振動幅值與孤立懸鏈線立管的對比情況如圖10所示。

圖10 懸鏈線立管模型振幅對比Fig.10 Amplitude comparison of catenary riser model

本文引入振幅折減效率η，以此表征干涉作用前、后懸鏈線立管模型振幅的變化情況，如式(1)所示。

(1)

式中：η表示振幅折減效率，yRMS表示孤立懸鏈線立管模型的橫流向振動幅值，yDRMS表示串聯(lián)排列下游懸鏈線立管模型的橫流向振動幅值。

η越接近于100%，表示干涉前、后懸鏈線立管模型的振幅差異越小，即說明下游懸鏈線立管模型受到的干涉效應(yīng)影響相對越弱。由圖10d綜合分析可得，下游懸鏈線立管模型非接觸測點(diǎn)1處所受干涉影響最弱，非接觸測點(diǎn)3處所受干涉影響最強(qiáng)。分析原因主要在于懸鏈線立管構(gòu)型獨(dú)特，不同截面所處的垂向位置存在一定偏差，下游懸鏈線立管模型從上端到下端與上游頂張力立管模型的水平間距逐漸減小，間距越小即受到的流動干涉效應(yīng)隨之越強(qiáng)，因此從非接觸測點(diǎn)1到非接觸測點(diǎn)3所受的干涉作用逐漸增強(qiáng)。

4 結(jié) 論

(1)孤立剛性立管模型的“鎖定”區(qū)間為Vr=5～9，在此區(qū)間內(nèi)流體與立管間產(chǎn)生劇烈的非線性相互作用，立管在橫流向持續(xù)以較大幅度振動，順流向振幅隨著約化速度的增大持續(xù)增大；立管橫流向振動幅值沿管長方向呈現(xiàn)出兩端小、中間大的分布特征，約化速度為12時可觀察到2階振型的出現(xiàn)。

(2)在剛性立管尾流干涉試驗(yàn)中，來流角度對下游立管模型的流動干涉作用有顯著影響。當(dāng)α=0°時，下游立管模型的“鎖定”區(qū)間相比孤立立管有一定滯后，在高約化速度下，下游立管橫流向振幅有所增大，受干涉作用影響順流向振幅小于孤立立管；隨著來流角度的不斷增大，下游立管所受流動干涉作用逐步減弱，橫流向振幅逐漸降低并接近于孤立立管，順流向振幅相應(yīng)提升并在數(shù)值上逐漸大于孤立立管。

(3)柔性孤立頂張力立管的振動特性與剛性立管相似，而孤立懸鏈線立管模型由于構(gòu)型獨(dú)特，導(dǎo)致在來流作用下多階振動模態(tài)共存，不同測點(diǎn)處的橫流向振動幅值差異相對較大。

(4)在柔性串聯(lián)立管尾流干涉試驗(yàn)中，受上游頂張力立管流動干涉影響，下游懸鏈線立管的響應(yīng)幅值與柔性孤立立管相比存在較大區(qū)別。由振幅折減效率分析可知，下游懸鏈線立管從非接觸測點(diǎn)1到非接觸測點(diǎn)3所受的干涉影響逐漸加強(qiáng)，即與上游立管水平間距越小的位置受到的流動干涉效應(yīng)越強(qiáng)。