朱一鳴， 王 磊， 張 濤， 劉翰林

(上海交通大學(xué) 海洋工程國家重點實驗室, 高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心, 上海 200240)

?

·實驗技術(shù)·

8 000 m3耙吸式挖泥船動力定位系統(tǒng)時域模擬與模型試驗

朱一鳴， 王 磊， 張 濤， 劉翰林

(上海交通大學(xué) 海洋工程國家重點實驗室, 高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心, 上海 200240)

針對耙吸式挖泥船動力定位系統(tǒng)，以某8 000 m3耙吸式挖泥船為研究對象，分析其動力定位系統(tǒng)在模擬海況下的定位能力與特點,構(gòu)建出動力定位時域模型。對挖泥船在各種海況下的實時運動進行時域模擬分析，得到該挖泥船的偏移半徑時歷和功率消耗時歷，并進行相關(guān)模型試驗。時域模擬與模型試驗結(jié)果表明,該挖泥船動力定位系統(tǒng)僅可在較小艏向角范圍內(nèi)完成有效定位，由于舵轉(zhuǎn)角的限制，螺旋槳無法在較大艏向角的情況下提供足夠的轉(zhuǎn)船力矩完成定位。通過模型試驗與時域模擬結(jié)果的分析與比較，為動力定位時域模擬程序進一步優(yōu)化提供參考。

耙吸式挖泥船； 動力定位系統(tǒng)； 時域模擬； 模型試驗

0 引 言

耙吸式挖泥船通常裝備有耙頭挖掘機具以及水力吸泥裝置，是一種裝倉式自航挖泥船，具有自挖、自載、自卸等特點，航行性能優(yōu)良[1]。其優(yōu)點主要表現(xiàn)在機動靈活、效率高、操作簡單、經(jīng)濟性好等方面。鑒于耙吸式挖泥船的諸多優(yōu)點，現(xiàn)代疏浚與圍海造田等大型工程已廣泛采用耙吸式挖泥船，具有廣闊的應(yīng)用前景與市場需求[2]。

隨著國際疏浚市場的不斷發(fā)展，動力定位系統(tǒng)逐漸地應(yīng)用在工程船舶上，對配備有動力定位系統(tǒng)的耙吸式挖泥船的需求也越來越大。動力定位系統(tǒng)是通過推進器產(chǎn)生的推力，使得海洋結(jié)構(gòu)物在海上保持一定的位置和艏向角的一種先進的定位系統(tǒng)。其基本原理是通過獲取船舶的位置誤差數(shù)據(jù)，由控制系統(tǒng)計算所需的推力，并傳遞到推力系統(tǒng)以產(chǎn)生推力，以此達到抵抗海洋外部環(huán)境載荷，從而減少、消除位置誤差的目的[3]。

目前一些主要的生產(chǎn)和使用廠家都一致認可動力定位能力曲線是評價船舶動力定位能力的有效工具[4]。動力定位能力曲線是一種靜態(tài)的動力定位能力分析方法，該方法將推力系統(tǒng)理論上產(chǎn)生的推力與靜態(tài)環(huán)境載荷相平衡，精度不高，可在初步設(shè)計階段中使用。更加精確地分析定位能力需采用動態(tài)的時域模擬分析以及模型試驗等方法。動態(tài)時域模擬與模型試驗方法能夠考慮到動態(tài)載荷，獲得實時位置偏移數(shù)據(jù)，能較為可靠地反映動力定位狀態(tài)下的海洋結(jié)構(gòu)物在各種海況下的定位狀態(tài)，從而驗證動力定位系統(tǒng)的有效性，并可為動力定位控制系統(tǒng)進一步優(yōu)化提供試驗依據(jù)和指導(dǎo)，相比于靜態(tài)的動力定位能力曲線能更加直觀地評價系統(tǒng)的定位能力，是更為理想的研究方法[5]。

目前，挖泥船動力定位系統(tǒng)研究較少，大多處于對挖泥船動力定位系統(tǒng)的理論研究或模型試驗的研究階段[6-7]，尚未將兩者很好的結(jié)合。本文針對耙吸式挖泥船這一特殊船舶，采用時域模擬以及模型試驗兩種方法進行動力定位系統(tǒng)研究，最終確定該耙吸式挖泥船的最優(yōu)作業(yè)方案，為工程實踐提供借鑒。

1 理論基礎(chǔ)

動力定位時域模擬是建立時域運動模型，分析海洋結(jié)構(gòu)物在動力定位系統(tǒng)控制下的運動與受力，時域模擬能較好地模擬真實海況，模擬結(jié)果更加精確，并獲得定位精度、功率消耗等定位信息，為動力定位控制系統(tǒng)的進一步優(yōu)化提供借鑒。耙吸式挖泥船動力定位時域模擬基本原理及各個模塊如圖1所示，主要包括環(huán)境載荷計算模塊、低頻運動方程模塊、控制模塊以及推力分配模塊[8]。

1.1 風載荷計算

風載荷是海洋結(jié)構(gòu)物所承受的外載荷中很重要的一部分，其常見的估算方法為模塊法，也是本文采用的估算方法，即把整個海洋結(jié)構(gòu)物離散分類，得到若干個標準模塊，分別估算載荷，將各個模塊載荷疊加獲得總載荷。風載荷的計算表達式為[9-10]：

圖1 挖泥船動力定位時域模擬原理

(1)

式中：Fxw,Fyw,Mxyw分別為縱向、橫向與艏搖風力；ρw為空氣密度；Cxw(φwR),Cyw(φwR)和Cxyw(φwR)分別為縱向、橫向與艏搖風力系數(shù)；AT，AL分別為艏向與側(cè)向受風面積；LPP為兩柱間長；vwR為海平面上方10 m相對風速。

1.2 流載荷計算

海洋洋流的存在，流載荷也是環(huán)境載荷的一部分，由于挖泥船的實時運動，故需考慮船與海流相對速度的影響，與風載荷類似，流載荷計算表達式為[11-12]：

(2)

式中：Fxc,Fyc,Mxyc即為縱向、橫向和艏向流力；ρc為海水密度；Cxc(φcR),Cyc(φcR)和Cxyc(φcR)分別為縱向、橫向和艏向流力系數(shù)；vcR為船體與海浪相對速度；d為船舶平均吃水。

1.3 波浪載荷計算

波浪載荷可以分為一階波浪載荷和二階波浪載荷。一階波浪載荷主要使海洋結(jié)構(gòu)物在平衡位置產(chǎn)生波頻運動，不會發(fā)生漂移；而二階波浪載荷則會使海洋結(jié)構(gòu)物發(fā)生緩慢的漂移。二階波浪載荷主要由平均、差頻和合頻三部分組成，在常規(guī)的計算中一般只考慮差頻力部分和平均力部分，而忽略合頻力部分。平均波浪力和差頻波浪力合稱波浪漂移力，其大小與入射波波高的平方成正比，頻域下其比例關(guān)系用二次傳遞函數(shù)表示。由二次傳遞函數(shù)P(ω1,ω2)與Q(ω1,ω2)經(jīng)傅里葉變換，可得到時域脈沖響應(yīng)函數(shù)[13-14]：

e(iω1τ1-iω2τ2)dω1dω2

(3)

(4)

若給定波浪時歷ζ(t)，則有波浪漂移力時歷：

ζ(t-τ2)dτ1dτ2

(5)

式中，i=1,2，…，6。

1.4 低頻運動方程

挖泥船在風、浪、流以及推力等共同作用下產(chǎn)生的運動可分解為波頻運動和低頻運動。波頻運動由一階波浪力引起，其幅值較大，推力器產(chǎn)生的推力不可能完全抵消。同時高頻波浪載荷引起的高頻運動僅表現(xiàn)為在平均位置處的周期性振蕩運動，所以低頻運動是動力定位系統(tǒng)所要抵抗的主要運動，運動時域模擬中也僅考慮低頻波浪力的作用。

挖泥船自由度低頻運動方程為[15-16]：

(6)

式中:u,v,w為線位移速度；p，q，r為角位移速度；m為挖泥船質(zhì)量，包含附加質(zhì)量；Ixx，Iyy，Izz，Izx為挖泥船慣性矩，包括附加質(zhì)量矩；X，Y，Z分別為x，y和z方向的外力；K，M，N分別為x，y和z方向的外力矩。

2 挖泥船動力定位時域模擬

2.1 挖泥船模型

本船艙容約8 100 m3，雙機雙槳復(fù)合驅(qū)動，用于沿海疏浚和吹填作業(yè)。主要船型參數(shù)如下：總長LOA=134 m,型寬B=22 m,型深D=9.2 m,設(shè)計吃水d1=6.8 m,作業(yè)吃水d2=8.2 m,排水量(設(shè)計吃水)T=15 700 t,排水量(作業(yè)吃水)T=19 379.4 t。

該船采用2個主推力器與1個電驅(qū)動艏側(cè)推來控制挖泥船的定位，主推力器直徑3.4 m，最大功率4 500 kW，艏側(cè)推直徑2 m，最大功率600 kW。同時該船還配備有2個半平衡半懸掛舵，最大舵角35°。推力器具體安裝位置如圖2所示。

圖2 推力器安裝位置示意

2.2 動力定位時域模擬

本動力定位時域模擬主要利用上海交通大學(xué)海洋工程國家重點實驗室的軟件DPSIM進行計算。根據(jù)該船用途以及作業(yè)環(huán)境要求，模擬風速15 m/s，最大流速1.029 m/s，波浪模擬采用JONSWAP波浪譜，其有義波高3 m，譜峰周期6 s，取風浪流同向聯(lián)合作用。由于動力定位更關(guān)注挖泥船在水平面的運動偏移，故時域模擬僅對挖泥船水平面的3個自由度運動進行控制，即：橫蕩、縱蕩和艏搖。根據(jù)船體推力系統(tǒng)配置的特點，無法在船體艏向與外界環(huán)境力方向夾角過大的情況下進行有效定位。故在時域模擬中艏向角自0°～20°每隔5°進行計算，包含不同流速下工況的對比，其中流速分別為1和2 kn。通過不同工況的時域模擬，最終得到挖泥船動力定位狀態(tài)下的水平偏移和總功率消耗的時歷曲線及其統(tǒng)計值。

2.2.1 偏移半徑計算結(jié)果

定位精度的評估方式是通過測量船體上控制點與目標位置水平方向偏差實現(xiàn)的，同時要兼顧艏向角的偏差。水平方向的偏移半徑公式如下：

(7)

偏移半徑統(tǒng)計值時域模擬結(jié)果如表1所示。挖泥船在不同艏向角、不同流速下的水平偏移時歷曲線如圖3所示。

根據(jù)以上時域模擬的結(jié)果可以發(fā)現(xiàn)，在相同偏移角度情況下，隨著流速的增加，定位精度下降。當流速增加到2 kn以上時，無法進行有效定位。在相同流速情況下，隨著艏向角的增加，定位精度下降。由于舵轉(zhuǎn)角的限制，無法提供更大的Y向力及偏轉(zhuǎn)扭矩，當艏向角增加到20°以上時，無法進行有效定位(見圖4)。

對于動力定位系統(tǒng)，普遍允許的偏移量為水深的2%～4%，對于本船工作狀態(tài)下的水深40 m而言，允許平均偏移量為0.8～1.6 m，從結(jié)果上看，僅當艏向角為0°流速為1 kn的工況是符合要求的。

表1 挖泥船偏移半徑統(tǒng)計值 m

2.2.2 功率消耗計算結(jié)果

時域模擬計算結(jié)果的消耗功率統(tǒng)計值如表2所示。挖泥船在不同艏向角不同流速下的功率消耗時歷曲線如圖5所示。

表2 挖泥船功率消耗統(tǒng)計值 kW

根據(jù)以上數(shù)值模擬的結(jié)果可以發(fā)現(xiàn)，在相同艏向角情況下，隨著流速增加，動力定位消耗功率增加，包括平均功率、最大功率及最小功率。在相同流速情況下，隨著船體偏轉(zhuǎn)角度增加，消耗總功率總體呈增加的趨勢。當艏向角為20°，流速為2 kn時，兩個主推進器及艏側(cè)推并沒有接近于滿負荷工作，而再增加偏轉(zhuǎn)角度或者流速時卻無法進行有效定位。說明無法定位的原因不是推進器無法提供足夠的推力，因為在保證偏移半徑盡可能小的前提下，同時要兼顧艏向角盡可能小的擾動，由于舵轉(zhuǎn)角的限制，無法提供足夠的偏轉(zhuǎn)扭矩使艏向角保持在目標角度，使得船體無法在目標位置平衡，從而偏離目標位置，無法進行有效定位。

3 動力定位模型試驗

模型試驗是研究動力定位系統(tǒng)的重要方法之一，通過在上海交通大學(xué)海洋工程國家重點實驗室海洋工程水池開展挖泥船動力定位模型試驗，獲得該船相關(guān)定位信息，包括運動、受力以及消耗功率等。根據(jù)試驗的條件與要求，模型縮尺比λ=50。模擬的環(huán)境載荷如表3所示。

有義波高/m譜峰周期/s風速/(m·s-1)流速/(m·s-1)風浪流方向?qū)嶋H值模型值實際值模型值實際值模型值實際值模型值同向30.0660.849152.1211.0290.145

動力定位更關(guān)注挖泥船在水平面內(nèi)的3自由度運動，故本試驗主要測量分析挖泥船縱蕩、橫蕩及艏搖的運動時歷。模型分別在0°、10°艏向角下進行實驗，共有2個實驗工況。最終，艏向角為0°時，挖泥船具有一定定位能力。各個螺旋槳功率消耗在一定范圍之內(nèi)，10°艏向角時，挖泥船定位并沒有成功，推力器功率并未達到峰值，說明挖泥船所布置的3個推進器和2個舵不足以合理分配抵抗外環(huán)境力。

將模型試驗數(shù)據(jù)換算成實船數(shù)據(jù)，試驗結(jié)果見表4和表5。圖6、7給出了已換算為實船的挖泥船模型試驗數(shù)據(jù)的時歷曲線。

表4 0°艏向角下的位置統(tǒng)計分析

表5 0°艏向角下的功率消耗統(tǒng)計分析

對于動力定位系統(tǒng)，普遍允許的偏移距離為水深的2%～4%，對于本船工作狀態(tài)下的水深40 m而言，允許平均偏移量為0.8～1.6 m，從結(jié)果上看，偏移半徑的統(tǒng)計平均值是符合要求的。本試驗設(shè)計海況對挖泥船比較惡劣，因此當遭遇大風大浪時，出現(xiàn)最大漂移半徑較大并且平均漂移半徑靠近允許范圍的極限值的情況。

圖6 偏移半徑時歷曲線

圖7 推進器消耗功率時歷曲線

圖8 艏搖角時歷曲線

當艏向角為0°時，大部分風浪流力可以由主推進器直接抵抗，當艏向角有一定偏移時，可通過艏側(cè)推進器進行調(diào)整來保持一定的定位位置。本試驗中，當艏向角偏移一定角度時，無法進行有效定位，螺旋槳并未達到最大功率，說明推力不足不是無法定位的原因，由于舵轉(zhuǎn)角的限制，無法提供足夠的偏轉(zhuǎn)扭矩使艏向角保持在目標角度，使得船體無法在目標位置平衡，從而偏離目標位置。由試驗結(jié)果不難發(fā)現(xiàn)，在試驗設(shè)計的海況下，挖泥船總體定位較好，若出現(xiàn)惡劣海況，仍會出現(xiàn)定位失敗的情況，因此，挖泥船在實施作業(yè)過程中最好避免上述模擬以及試驗海況。

4 結(jié)果分析

本文通過數(shù)值時域模擬和模型實驗的方法從定位精度和功率消耗兩個角度對比分析挖泥船的動力定位系統(tǒng)，得到以下結(jié)論：

(1) 在相同艏向角度情況下，隨著流速的增加，動力定位的定位精度下降，且消耗功率增加。當流速增加到2 kn以上時，無法進行有效定位。

(2) 在相同流速情況下，隨著艏向角的增加，動力定位的定位精度下降，且消耗功率有增加趨勢。由于舵轉(zhuǎn)角的限制，無法提供更大的偏轉(zhuǎn)扭矩使艏向角保持在目標角度，使得船體無法在目標位置平衡，從而偏離目標位置。

(3) 模型試驗的結(jié)果與時域模擬計算得到的結(jié)果共同表明，該挖泥船僅能在艏向角為0°，流速較小(小于2 kn)的情況下，具有一定的定位能力。

(4) 模型試驗的結(jié)果與時域模擬計算得到的結(jié)果存在一定的偏差，與試驗過程中測量裝置的誤差有一定的聯(lián)系，但是總體趨勢吻合較好。

5 結(jié) 語

本文針對研究較少的耙吸式挖泥船動力定位系統(tǒng)，以某8 000 m3耙吸式挖泥船為研究對象，從定位精度和功率消耗兩個角度出發(fā)，采用模型試驗與時域模擬兩種方法對動力定位系統(tǒng)進行研究。通過構(gòu)建動力定位時域模型，對挖泥船在各種海況下的實時運動進行時域模擬分析，得到該挖泥船的偏移半徑時歷和功率消耗時歷，并進行模型試驗，將模型試驗與時域模擬結(jié)果的分析與比較，得到該挖泥船定位能力與相關(guān)信息，為實船動力定位系統(tǒng)提供參考與優(yōu)化空間，并為將來工程實踐應(yīng)用提供必要的理論與實踐指導(dǎo)。

[1] 林吉良. 耙吸式挖泥船疏浚仿真[D]. 南京: 河海大學(xué), 2004.

[2] 張躍東. 耙吸式挖泥船動力定位系統(tǒng)工程化研究[D]. 南京： 江蘇科技大學(xué) 2010

[3] 申 輝. 動力定位推力系統(tǒng)效率研究[D]. 上海: 上海交通大學(xué), 2012.

[4] 陳 恒. 深海半潛式平臺動力定位推力系統(tǒng)設(shè)計研究 [D]. 上海 上海交通大學(xué), 2008.

[5] 李 博. 動力定位系統(tǒng)的環(huán)境力前饋研究[D]. 上海 上海交通大學(xué), 2013.

[6] 袁 偉, 俞孟蕻, 丁樹友, 等. 耙吸挖泥船動力定位系統(tǒng)模型試驗研究[J]. 船舶工程, 2014(3): 23.

[7] 李德珊. 耙吸式挖泥船運動建模及動力定位控制研究[D]. 大連: 大連海事大學(xué), 2013.

[8] 楊 歡, 王 磊, 申 輝. 深水半潛式鉆井平臺動力定位實時功率模擬研究[J]. 海洋工程, 2012(2): 4.

[9] Sahin I, Aybar A. A survey of research on wind loads on semisubmersibles [J]. OCEANS, 1984, 12(3): 548-555.

[10] Ir F van Walree. Wind and Current Loads on Semi-subersibles [R]. MARIN:No.49817- 5- SE,1991．

[11] Wichers J E W. A Simulation Model for a Single Point Moored Tanker [M]. Publ. 797. MARIN, Wageningen, the Netherlands, 1988.

[12] Faltinsen O M. Sea Loads on Ships and Offshore Structures [M]. Cambridge University, 1990.

[13] Pinkster J A. Low Frequency Second Order Wave Exciting Forces on Floating Structures [M]. MARIN Publication, 1980.

[14] 劉應(yīng)中,廖國平. 船舶在波浪上運動理論[M]. 上海 上海交通大學(xué)出版社,1987

[15] Skejic R, Faltinsen O M. A unified seakeeping and maneuvering analysis of ships in regular waves [J]. Journal of Marine Science and Technology, 2008, 13(4): 371-394.

[16] 王 磊, 孫 攀, 王 亮. 深水半潛平臺動力定位時域模擬研究[J]. 中國科學(xué): 物理學(xué) 力學(xué) 天文學(xué), 2011(2): 3.

Numerical Research and Model Test on the Dynamic Positioning System of a 8 000 m3Drag Suction Dredger

ZHUYi-ming,WANGLei,ZHANTao,LIUHan-lin

(State Key Laboratory of Ocean Engineering, Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration, Shanghai Jiaotong University, Shanghai 200240, China)

The research focuses on the dynamic positioning system on a 8 000 m3drag suction dredger and builds the model of low frequency motion equation and model of thruster system which leads to a time domain simulation for analyzing the dynamic position system. The environmental loads include wind, wave and current. By these, the horizontal offset and power consumed as well as their statistic values can be analyzed. At the same time, model test is also conducted to verify the results The computed results agree well with test data. The results of simulation and model test show that the dredger can keep its state in certain environment. The distribution of thrusters may not provide enough moment to help the vessel keep its state. The result can be easily introduced into engineering practice and provide designers with useful information when one is designing a dynamic positioning system on a drag suction dredger.

drag suction dredger; dynamic position system; time domain simulation; model test

2016-01-26

國家自然科學(xué)基金項目(51179103)； 國家重點基礎(chǔ)研究發(fā)展計劃項目(2013CB036103)

朱一鳴(1991-)，男，江蘇淮安人，碩士生，從事動力定位系統(tǒng)研究。Tel.:18221372875； E-mail:zym262513@sjtu.edu.cn

王 磊(1971-)，男，安徽蚌埠人，副教授，從事動力定位系統(tǒng)研究。Tel.:021-62932025； E-mail:wanglei@sjtu.edu.cn

P 751

A

1006-7167(2016)09-0004-06