基于AQWA的沉船打撈作業安全性分析*

張　明　代　路　甘　進　李佳衛　高曉圓

(武漢九通汽車廠裝備研究院1)　武漢　430034)　(武漢理工大學交通學院2)　武漢　430063)

基于AQWA的沉船打撈作業安全性分析*

張明1)代路1)甘進2)李佳衛2)高曉圓2)

(武漢九通汽車廠裝備研究院1)武漢430034)(武漢理工大學交通學院2)武漢430063)

摘要：打撈工程中，起重船舶在水上安全作業十分重要，尤其是在起吊沉船過程中.兩船與錨鏈，吊纜等在風、浪、流的耦合作用下的運動響應，以及吊纜、生根位置的合理配置是打撈工程的關鍵影響因素.基于多體水動力學軟件AQWA對沉船起吊過程進行了數值仿真，分析了1 000 t全回轉起重船在打撈沉船的過程中起重船與沉船的運動響應，提出了一種基于AQWA的沉船打撈作業安全評估方法.

關鍵詞：沉船打撈；起重船；AQWA；運動響應

0引言

隨著長江航運的發展和標準化船型進程的推進，通行船舶呈現出噸位大型化，船型結構復雜化的趨勢.加上長江干線流速大、易淤積，尤其是三峽庫區深水大流速的等復雜打撈環境條件，給沉船打撈工作帶來較大難度和挑戰[1].長久以來，由于打撈設備，打撈技術，作業環境，打撈人才缺乏等多方面的原因，長江干線沉船打撈一直是一個高風險、粗放式的作業模式，尤其是沉船結構受力計算，船舶姿態監控，復雜打撈環境對沉船打撈的影響分析計算等都是嚴重依靠打撈人員和潛水員的經驗判斷，或者邊撈邊試，或者以幾倍的打撈代價參與打撈，難以達到精準計算、快速高效地進行沉船打撈的目的.

AQWA軟件在工程中的應用十分廣泛，從單一浮體的運動響應到多體耦合運動響應計算都有所研究.它的求解方法基于勢流理論[2]，首先求解在單位簡諧波中，船舶六自由度運動的幅值響應算子RAO，再通過線性系統的分析方法以及頻率特性到脈沖響應函數進行傅里葉變換，可以得到船舶在規則或者不規則波下的時域運動響應預報[3].李宇辰等[4]以方形系數為0.7的系列60船為例，利用AQWA 計算得到的船體運動響應合理地反映了系列60船的實際耐波性能，可用于船舶波浪載荷預報.胡毅等[5]應用AQWA-DRIFT和AQWA-LINE 2個模塊，研究了大型LNG船碼頭系泊時，在風、浪、流聯合載荷作用下船舶的運動響應，得到了LNG 船碼頭系泊時整體運動響應及系纜繩所受張力.

本文著眼于對起重船起吊沉船過程中的運動安全性評估，基于AQWA軟件對整個過程中船舶的運動響應進行數值分析，分析數據可為吊船打撈沉船技術方案設計提供參考.

1計算流程

本文采用多體水動力學軟件AQWA建立1 000 t起重船與沉船的三維模型，模擬打撈工程中起吊沉船的作業過程.主要應用到AQWA-LINE與AQWA-DRIFT 2個模塊計算起重船和沉船的運動響應.

在AQWA-LINE模塊中，計算頻域下的浮體(或潛體)在單位簡諧波下的波浪載荷以及運動響應.在該模塊采用有限元加邊界元的方法，在船殼上劃分有限元，在有限元上布置源匯來模擬船殼的物面邊界條件，用邊界元法求解格林函數得到浮體結構的波浪力，并且可以得到一階波浪力、二階漂移力以及物體的附加質量、輻射阻尼等，為后續的計算提供輸入參數.

在AQWA-DRIFT模塊中，模擬浮體(或潛體)在不規則波中實時的運動響應，可以考慮波頻運動以及低周震蕩漂移運動，同時在此模塊下也能考慮作用在物體上的風和流的作用以及通過懸鏈線或鉸鏈等方式考慮多個浮體(或潛體)之間的耦合作用，風和流的模擬通過在重心上施加力與力矩來實現，從而實現吊船打撈起吊作業模擬，計算流程見圖1.

圖1　基于AQWA的吊船打撈計算流程

2理論基礎

線性勢流理論中，對于不可壓縮的理想流體，無旋即有勢，速度勢需要滿足以下定解條件.

線性自由面條件：

在浮體的濕表面s(t)上有物面條件：

底部有邊界條件：

將總的速度勢中分出入射波部分和擾動勢φp即

擾動勢φp也是簡諧的，滿足疊加原理，有理由把擾動勢分解成7個組成部分，即

式中：φj(j=1,2,…,6)為稱為輻射勢；φ7為繞射勢.

得到了上面所述的三維浮體諸輻射勢φj(j=1,2,…,6)和繞射勢φ7，用伯努利方程(只保留一階項)可寫出作用于浮體的一階動壓力.即

對其在物面上積分可得流體動力(F1,F2,F3)，kN，和動力矩(F4,F5,F6)，kN·m.

其中

式中：f0k為入射波的力(矩)，f7k為繞射力(矩)，這二者之和f0k+f7k是浮體在規則波種受到的波浪激勵力(矩)，f0k占主要成分，又稱為弗勞德-克雷洛夫力(矩)，Tkj為浮體做單位速度j態運動時受到的k方向輻射力(矩).

通過速度勢的計算不難得到浮體的運動響應微分方程：

式中：Ms為結構質量，Ma為附加質量，B為阻尼系數，K為水靜力剛度，F為波浪力.

利用線性系統的分析方法，頻率特性和脈沖響應函數互為傅里葉變換，輻射問題的時域解表達形式為：

式中：vj為浮體的j態運動速度.

同樣得到繞射勢的是與表達形式

式中規定，x=y=0處波面起浮為η(t).

3計算模型

3.1尺度參數

表1為1 000t起重船主尺度，表2為沉船主尺度參數，分別建立2艘船的三維幾何模型.起重船吊機選取固定式尾吊最大起重能力，吊高75m.圖2為在劃分網格后的起重船和沉船運動響應計算模型.

船體上的質量分布對船體在波浪中的運動影響很大，將船體重量的統計值以質量點的形式加在重心位置，也可以分布多個質量點來模擬載荷的分布情況，本文選擇第一種加載方式.

表1　起重船主尺度參數

表2　沉船主尺度參數

圖2　運動響應計算模型

3.2系泊及吊纜布置

1 000 t起重船同時配備有錨泊定位與動力定位系統，以應對三峽庫區中大水深的特殊情況，如無法拋錨定位，則采取動力定位.在較淺水深區域作業時可以以錨泊定位為主，輔助以動力定位.本文中計算水深取60 m，不考慮動力定位的作用.圖3為起重船錨泊定位布置圖.表3為所設置的錨鏈參數.

圖3　起重機系泊布置

參數數值材料鋼纜長度/m215等效截面積/m20.01自重/(kg·m-1)150剛度EA/kN600000

在實際的打撈作業中，沉船的重量重心往往是比較難以獲取到的信息，主要由于沉船在水中的情況不明，水下探測的作業難度大，精度差，加上一些船舶不規范的運營或改裝，使得可搜集的資料不完全、不準確，使得打撈作業難以開展.沉船重量重心往往也是通過估算獲得，精度不高.在進行數值模擬時，可以輸入估算的重量重心位置，在計算機中模擬起吊的過程，從而對沉船上的生根位置進行調整.本文中將起重船的起吊掛鉤位置垂直置于沉船重心上方.本例中起吊前沉船的狀態為正坐于水底，4根吊纜對稱布置在重心位置的兩側，見圖4.沉船生根位置到吊鉤的直線距離為140 m，即吊纜的長度，選擇剛度為500 kN·m-1的線彈性繩進行模擬.

圖4　沉船起吊吊纜布置

3.3計算工況及環境參數

本計算所去的風、浪、流等環境載荷參數見表4.其中，選擇不規則PM波譜模擬波浪環境，選擇的波浪參數偏于危險情況風浪流的方向為180°表示沿船首指向船尾，90°則表示垂直于船長方向作用在船側，以下工況分別計算風浪流沿縱向作用與橫向作用兩種情況.

通常情況下被打撈的船只、物品都是沉坐在水底，打撈力除了需要克服其自身水中重量以外，還需要克服泥沙對沉船的吸附力，吸附力的大小和作用點與沉船水中重量、船體與水底泥沙的接觸面積、形狀、坐沉時間、底質性質等因素有關，十分復雜，無法進行精確計算.在實際工程中往往簡化認為吸附力與沉船水中重量成正比，以吸附力系數的形式添加到沉船水中重量中，以此進行打撈力的計算，吸附力系數同時也與水底底質有關.根據調研資料，取吸附力系數10%.

整個打撈起吊過程的模擬分為3個階段，沉底階段、起吊階段、出水階段.沉底階段中沉船的水中重量為另兩個階段重量的110%.

表4　環境載荷參數

4計算結果

4.1頻域計算

研究表明，不規則的海浪可以簡化為無數多個頻率、方向和波幅不同的簡諧波的疊加，可視為均值為零的正太隨機過程.頻域分析就是為了計算出浮體在不同頻率、方向的單位波高簡諧波種的運動響應，根據線性疊加原理，即可計算出在不規則波、各波高中的運動響應.對于本文中的起重船，起吊位置在船尾部，起吊過程的縱搖響應以及在橫向外載荷作用下的橫搖響應值得著重關注.

圖5為各浪向下起重船橫搖幅值響應算子，圖6為縱搖幅值響應算子.圖中，對橫搖響應最大的浪向為90°，即波浪垂直于船長方向，在波頻為0.137 Hz(圓頻率即為0.863 rad/s)時出現峰值，此時根據微幅波理論，此時的水深處于深水情況，易求得此危險工況的波長為82.8m，波浪周期約為7.28 s.

圖5　橫搖幅值響應算子

圖6　縱搖幅值響應算子

而當起重船與來浪呈180°或135°，即迎浪作業與迎浪45°時，縱搖響應最大，根據上式，135°時波頻為0.124 Hz(圓頻率0.778 rad/s)，此時波長為101.8 m，周期為8.07 s.180°時波頻為0.097 Hz(圓頻率0.609 rad/s)，此時波長為166 m，周期為10.32 s.

4.2時域計算

在時域計算中考慮系泊、吊纜、風浪流以及兩艘船之間的耦合作用，得到兩艘船的六自由度的運動，通過統計運動響應的時程曲線，提取運動響應幅值、頻率等信息，對起重船的安全性加以評估.圖7為耦合作用下吊船打撈作業時域計算模型.

圖7　吊船打撈時域計算模型

通過對打撈過程的時域分析，分別得出起重船和沉船的運動響應時程曲線.圖8為浪向角180°時沉船在起吊狀態下的橫搖響應時間歷程曲線(此時對沉船為橫向載荷).由于起重船尾吊作業，起吊過程中將產生較大幅度的縱搖與縱蕩，而當起重船受橫向外載荷作用時，橫搖運動響應較大，而橫蕩運動雖受低頻二階慢漂力影響較大[6]，但有錨鏈提供恢復力，不作為考察重點.

圖8　沉船起吊工況下受縱向外載荷的橫搖響應時程曲線

將沉船與起重船在三種起吊狀態、兩種外載荷方向下的垂蕩、橫搖以及縱搖的時間歷程曲線進行匯總，提取運動響應的幅值整理于表5，表6.

通過分析比較沉船與起重船在三種起吊狀態、兩種外載荷方向下的垂蕩、橫搖以及縱搖的時間歷程曲線及運動響應幅值可知：

1) 在每一個工況下，起重船的橫搖響應的極值的絕對值十分相近，也就是說起重船在正浮位置左右橫搖震蕩，可見掛鉤受力位置在起重船中縱線上，沒有產生橫傾的力矩.而垂蕩與縱搖的幅值相差較大，縱搖幅值最為明顯，可見尾吊作業導致了起重船產生了較為嚴重的尾傾，在此尾傾角度附近產生縱搖震蕩，尾傾在1°～3°，沉底與出水階段的尾傾最大，這2個工況為最不利工況，應在打撈過程中予以重視.

表5　起重船運動響應幅值匯總

表6　沉船運動響應幅值匯總

2) 起重船與沉船的運動響應均會受外載荷方向的不同而變化，起重船的運動響應明顯小于沉船運動響應，說明在實際吊船打撈施工中為保證整個沉船起吊過程的安全性重點之一是采取相應措施控制沉船在打撈過程中的運動響應.

5結 束 語

本文基于AQWA軟件對1 000 t全回轉起重船尾吊打撈某沉船進行仿真計算，考慮兩船與風、浪、流、系纜以及吊纜之間的耦合作用，選取了沉底、起吊、出水3種工況，進行了頻域與時域的分析.驗證了應用AQWA軟件開展沉船吊船打撈計算和分析的可行性，并提出了一套基于AQWA的吊船打撈起重船和沉船運動響應計算方法，可為打撈工程中吊船打撈方案設計和制定提供有利的依據.

參 考 文 獻

[1]王兵．長江航道局三峽庫區應急搶險打撈起重船建設方案研究[J]．中國水運(下半月)，2014(1)：51-56.

[2]戴仰山,沈進威,宋競正,船舶波浪載荷[M].北京：國防工業出版社，2007.

[3]戴遺山,段文洋.船舶在波浪中運動的勢流理論[M].北京：國防工業出版社，2008.

[4]李宇辰，吳晞，袁小龍，等．耐波性預報的數值方法[J]．艦船科學技術，2013(10)：67-72.

[5]胡毅，胡紫劍，劉元丹，等．AQWA的大型LNG船碼頭系泊分析[J]．艦船科學技術，2012(3):84-89.

[6]KIM Y. Dynamic analysis of multiple-body floating platforms coupled with mooring lines and risers[D]．College Station，TX：Texas A&M University，2003．

Safety Analysis of Shipwreck Salvage Operations Based on AQWA

ZHANG Ming1)DAI Lu1)GAN Jin2)LI Jiawei2)GAO Xiaoyuan2)

(WuhanJiutongCarFactory,EquipmentResearchInstitute,Wuhan430034,China)1)(SchoolofTransportation,WuhanUniversityofTechnology,Wuhan430063，China)2)

Abstract：The security of crane ship is very important n the process of salvage, especially in lifting shipwrecks. The motion response of crane ship, shipwrecks, anchor chains and lifting ropes associated with wind, wave and current and the arrangement of ropes and roots have significant contribution to wreck salvage. The process of a shipwreck salvaged by a 1000 t fully swing crane ship is simulated by using hydrodynamics software AQWA. The motion responses of crane ship and shipwreck are obtained, and a security evaluation method of ship-sunken salvage process is put forward.

Key words：ship-sunken salvage; crane ship; AQWA; motion response

收稿日期：2016-03-12

中圖法分類號：U674.35

doi:10.3963/j.issn.2095-3844.2016.03.015

張明(1977- )：男，助理工程師，主要研究領域為工程結構安全可靠性

*中央高校基本科研業務費專項資金資助(2013-IV-069)