水動力相互作用對MOB 動力響應的影響研究

2022-10-19 10:43:10張波陳林董明海王化明

艦船科學技術 2022年17期

張波，陳林，董明海，王化明

(1.浙江國際海運職業技術學院，浙江舟山 316021；2.浙江海洋大學東?？茖W技術學院，浙江舟山 316022；3.浙江海洋大學船舶與海運學院，浙江舟山 316022)

0 引言

MOB 是一種能夠漂浮在海上的多浮體系統，是一種海洋環境、浮體與連接器相互耦合的動力學系統。連接器主要分為柔性連接器和剛性連接器。由于連接器限制了浮體的自由運動，所以連接器承受著極大的載荷。準確預報連接器載荷是目前MOB 研究領域的重點內容之一。和單浮體水動力計算不同，當多個浮體存在于同一流場時，由于物面條件變化將引起流場輻射勢和繞射勢的改變，這就導致多浮體水動力系數計算比單個浮體復雜得多。

多模塊MOB 動力響應的計算方法主要采用頻域法和時域法，所采用的計算模型主要有RMFC，RMRC，FMFC 和FMRC 模型。余瀾對超大型海上移動基地進行了較為深入的研究，涉及MOB 的動力響應，考慮了浮體水彈性的影響，并進行了實驗驗證。王永恒等提出由3 個相同模塊組成的超大型浮體，研究了不同環境不同連接器剛度下連接器載荷并與試驗結果進行了比較，通過比較提出了較為合理的連接器剛度值。宋斌等提出了一種由多個浮式平臺合圍成的人工港灣式超大型海上基地VLFOB，設計了一種可分離式連接器裝置將8 個模塊連接成整體。研究了不同連接器軸向剛度對VLFOB 固有頻率的影響，得出了連接器剛度的合理范圍避免與海洋環境條件產生共振。魯達以鉸接多浮體為研究對象，研究了浪向、波頻、阻尼等因素對波浪能發電裝置發電功率的影響，并提出了一種用于多浮體水彈性分析的簡化計算方法。李志偉以島礁附近的超大型浮式結構物為研究對象，研究了海底變化對浮體運動響應和連接器載荷的影響，并將數值計算結果與試驗結果進行了比較。L.Sun計算了兩鉸接方形浮式結構物的水動力響應并與Newman 的算例作比較以驗證程序正確性，又計算了2 艘大型液化天然氣船并排?？繒r在波浪中的運動響應和剛性連接器載荷響應。J.N.Newman計算了2 個鉸接連接的駁船在波浪中的運動響應，2 條駁船尺寸相同，長40 m，寬10 m，吃水5 m。J.N.Newman 的論文中又提供了另外一個算例，以1，3，5，9，17 個直線排列的圓柱形浮體陣列為研究對象，計算每個圓柱體的波浪力和水動力系數。

1 計算原理

本文研究的MOB 系統由7 個模塊組成，模塊編號為～，共有12 個連接器，連接器編號為～。整體坐標系和局部坐標系的定義如圖1 所示。

圖1 MOB（7MS）坐標系定義Fig.1 Coordinate system of MOB(7MS)

當多個模塊位于同一流場中時，總速度勢表達為：

忽略多體間相互影響的剛性浮體柔性連接器模型計算方法與考慮多體間相互影響的計算方法主要區別在于水動力系數矩陣的組成不同。當考慮浮體間相互影響時，第模塊第階輻射力，可表示為：

考慮浮體間相互影響時，認為第模塊第方向的運動會影響第模塊第方向上的力（矩）。如果不考慮浮體間相互影響，則在計算水動力系數時將多浮體系統視為單個剛體。

將個模塊的質量矩陣合成一個總質量矩陣，總質量矩陣M和靜水恢復力系數矩陣K可以表達為如下：

式中：M為第模塊的質量矩陣；K為第模塊的靜水恢復力矩陣。

連接器載荷為連接器總剛度矩陣乘以連接器位移矩陣，將連接器視為線性彈簧，表達式如下：

Δu為兩相鄰模塊和+1的位移差。

位于第和+1兩個模塊間的連接器剛度矩陣可表示為：

個模塊組成的多浮體系統的總剛度矩陣可表示為：

頻域下多浮體運動方程的矩陣表達式：

式中：M為系統質量矩陣,6×6；為附加質量矩陣，6×6；為阻尼系數矩陣，6×6；K為連接器總剛度矩陣，6×6；K為靜水恢復力矩陣，6×6；和矩陣大小 6×1。

短期預報所采用的波浪譜為Bretshneider 譜，北太平洋開敞海域7 級海況（SS7）的特征波高為7.5 m，譜峰周期為13.8 s。

2 計算模型

MOB 浮體模塊采用半潛式結構，下浮體長度270 m，設計吃水28.75 m，設計排水量2.975×10kg，方向慣性半徑50.412 m，方向慣量半徑80.494 m，方向慣量半徑90.671 m。網格劃分后的MOB 水動力模型如圖2 所示。

圖2 模塊水動力計算網格Fig.2 Mesh of MOB module

3 結果與分析

本文計算8 個浪向角（0°，15°，30°，45°，60°，75°，85°，90°），46 個波浪頻率（0.1～1.5 Hz），2 種連接器剛度下浮體的運動響應和連接器載荷，并在北太平洋7 級海況下進行短期預報。計算過程中，連接器剛度選擇（K=1 0N/m，K=1 0N/m，K=1 0N/m）和（K=1 0N/m，K=1 0N/m，K=1 0N/m）。對采用2 種計算方法得出的結果進行比較時，用“U-”表示不考慮浮體間相互作用時的數值計算結果，“C-”表示考慮浮體間相互作用時的數值計算結果。

本文通過自主開發程序Multifloat 計算MOB 系統的波浪誘導響應。為了驗證本程序的有效性，計算了5 模塊MOB 系統在連接器剛度 45波浪入射角下連接器的垂向載荷，并與文獻[11]中的結果進行了比較（見圖3），兩者吻合良好。5MS 表示5 個模塊組成的MOB 系統。

圖3 C3 垂向載荷（5MS，K6）Fig.3 RAO of Fz of C3 (5MS,K6)

3.1 波浪激勵力和水動力系數矩陣計算

利用水動力計算軟件SESAM 計算出考慮浮體間水動力相互作用的波浪激勵力矩陣和水動力系數矩陣（主要包括附加質量矩陣、阻尼系數矩陣、靜水恢復力矩陣）。圖4 為波浪入射角為 0時7 個模塊的波浪激勵力（矩）RAO。由圖可知，不同模塊的波浪力（矩）存在一定差別，在特定頻率下差別較為顯著。由圖4(c)和圖4 (d)可知，當波浪入射角為 0時，與各浮體的縱搖波浪力矩相比，橫搖波浪力矩明顯較小但不為零，這可能是由于模塊間水動力相互作用導致的。

圖4 各模塊波浪力(矩)幅值Fig.4 RAO of force and moment

圖5 附加質量系數Fig.5 Added mass coefficients

3.2 連接器載荷響應

圖6為和剛度下不同波頻和波浪入射角度下連接器縱向和垂向載荷幅值圖，計算過程中均考慮了模塊間的相互作用。由圖可知，連接器剛度下，波頻在0.3～0.5 之間時連接器載荷值較大，其余頻率下連接器載荷較小。連接器剛度下，波頻位于0.3～0.5 時，連接器縱向和垂向載荷值較大，當波頻為0.475，波浪入射角為 7 5時，載荷達到極值。

圖6 K5和K11 剛度下C1 連接器載荷幅值Fig.6 Load amplitude of C1 under K5 and K11

3.3 模塊間相互作用對浮體運動的影響

考慮浮體間相互作用時的計算結果是最為準確的，但是由于求解模塊間相互作用的水動力系數矩陣過程較為復雜。為了簡化計算，采用不考慮浮體間相互影響的方法進行水動力計算目前也是比較流行的。

不考慮模塊間相互作用時得到的計算結果與考慮模塊間相互作用時的結果存在一定差別。圖7 和圖8分別為連接器剛度下采用2 種計算方法得到的模塊縱蕩和垂蕩響應幅值。從圖中可以看出2 種方法得出的模塊運動值存在一定差別，考慮水動力相互作用時模塊的縱蕩和垂蕩響應極值略大于不考慮相互作用時的值，2 種方法計算得出的模塊縱蕩響應極值最多相差27.1%。

圖7 M1 模塊縱蕩響應幅值（K5）Fig.7 Surge RAO of M1(K5）

圖8 M1 模塊垂蕩響應幅值（K5）Fig.8 Heave RAO of M1(K5）

3.4 模塊間相互作用對連接器載荷的影響

圖9 和圖10 分別為剛度下連接器縱向和垂向載荷。由圖可知，不同波浪入射角度下連接器載荷極值差異較大。在同一個波頻下使用2 種計算方法得出的連接器縱向和垂向載荷值存在一定差別，但差別并不顯著。

圖9 C1 縱向載荷RAO（K5）Fig.9 RAO of Fx of C1(K5)

圖10 C1 垂向載荷RAO（K5）Fig.10 RAO of Fz of C1(K5)

圖11為和剛度時北太平洋7 級海況下和連接器縱向和垂向載荷的短期預報值。由圖可知，剛度下連接器縱向載荷極值比剛度時小，但剛度下連接器縱向載荷比剛度時大。當浪向角在 6 0以下時，連接器載荷較小。隨著浪向角增大，連接器各方向載荷有增大的趨勢，0～ 6 0之間相對平緩，60之后顯著增大，并在 75～ 8 5時達到了峰值。采用2 種計算方法得到的連接器載荷短期預報值雖然存在一定偏差，但是變化趨勢是一致的。