許 超，孫麗萍，趙志堅(jiān)，張利軍

（1.哈爾濱工程大學(xué) 船舶工程學(xué)院，哈爾濱 150001；2.中遠(yuǎn)船務(wù)工程集團(tuán)技術(shù)中心，遼寧大連 116600）

風(fēng)電安裝船運(yùn)動(dòng)響應(yīng)分析及模型試驗(yàn)研究

許 超1，孫麗萍1，趙志堅(jiān)2，張利軍2

（1.哈爾濱工程大學(xué) 船舶工程學(xué)院，哈爾濱 150001；2.中遠(yuǎn)船務(wù)工程集團(tuán)技術(shù)中心，遼寧大連 116600）

為了分析風(fēng)電安裝船在實(shí)際海況下的運(yùn)動(dòng)響應(yīng)及阻力性能，基于三維勢(shì)流理論，本文對(duì)風(fēng)電安裝船進(jìn)行頻域和時(shí)域運(yùn)動(dòng)響應(yīng)數(shù)值分析，并開(kāi)展風(fēng)電安裝船模型試驗(yàn)，對(duì)航行阻力、幅頻響應(yīng)函數(shù)和在位運(yùn)動(dòng)響應(yīng)進(jìn)行試驗(yàn)研究。通過(guò)試驗(yàn)得到風(fēng)電安裝船在航速13 kn時(shí)阻力為899 kN，RAO響應(yīng)曲線和不規(guī)則波中的運(yùn)動(dòng)幅值驗(yàn)證了數(shù)值計(jì)算的準(zhǔn)確性。研究結(jié)果表明風(fēng)電安裝船在十年一遇的北大西洋海況中運(yùn)動(dòng)幅值較大，不宜進(jìn)行安裝作業(yè)。

風(fēng)電安裝船；模型試驗(yàn)；耐波性；阻力

0 引言

海上風(fēng)電場(chǎng)的開(kāi)發(fā)已受到越來(lái)越多的重視和關(guān)注，其中的一個(gè)難點(diǎn)在于風(fēng)機(jī)整套設(shè)備的海上安裝。風(fēng)電安裝船是近些年來(lái)新發(fā)展的一種特種工程船舶，它將運(yùn)輸船、海上作業(yè)平臺(tái)和起重船的各項(xiàng)功能融為一體，可以獨(dú)立完成整套海上風(fēng)機(jī)設(shè)備的運(yùn)輸、安裝和維護(hù)。中遠(yuǎn)船務(wù)技術(shù)中心作為設(shè)計(jì)單位，對(duì)某型號(hào)風(fēng)電安裝船母型船進(jìn)行了優(yōu)化設(shè)計(jì)，以滿(mǎn)足航行性能和作業(yè)性能的要求。本文以中遠(yuǎn)船務(wù)提供的優(yōu)化船型為研究對(duì)象，對(duì)新船型進(jìn)行數(shù)值仿真和模型試驗(yàn)研究，以獲得新船型在不同航速下的阻力性能和在實(shí)際海況中的運(yùn)動(dòng)響應(yīng)。為風(fēng)電安裝船的設(shè)計(jì)開(kāi)發(fā)和在實(shí)際海況中的作業(yè)提供參考［1］。

1 水動(dòng)力性能數(shù)值分析

首先根據(jù)三維勢(shì)流理論方法，在整個(gè)頻域范圍內(nèi)對(duì)風(fēng)電安裝船進(jìn)行輻射和繞射分析，獲得在不同頻率單位波幅規(guī)則波作用下的水動(dòng)力參數(shù)，包括一階波浪力傳遞函數(shù) RAO、二階波浪力傳遞函數(shù)QTF、附加質(zhì)量以及輻射阻尼等，頻域分析的水動(dòng)力結(jié)果為后續(xù)風(fēng)電安裝船的時(shí)域運(yùn)動(dòng)響應(yīng)分析打下基礎(chǔ)［2］。

1.1 頻域水動(dòng)力分析

應(yīng)用輻射和繞射分析計(jì)算理論對(duì)無(wú)航速的風(fēng)電安裝船在規(guī)則波下的運(yùn)動(dòng)進(jìn)行分析，對(duì)于流場(chǎng)中的浮體，浮體在流場(chǎng)中的壓力根據(jù)線性伯努利方程，可得：

對(duì)浮體整個(gè)濕表面單元進(jìn)行壓力積分，則作用在風(fēng)電安裝船上的波浪激勵(lì)力為：

式中，fi為i方向的波浪激勵(lì)力；ni為面元法向量；S為浮體濕表面。

浮體在受到波浪激勵(lì)力的同時(shí)，浮體自身的運(yùn)動(dòng)也帶動(dòng)附近流場(chǎng)的變化，即產(chǎn)生了輻射勢(shì)，則浮體對(duì)流體的反作用力可以寫(xiě)成：

式中，μij為浮體附加質(zhì)量；λij為輻射阻尼［3］。

將波浪激勵(lì)力和輻射力線性疊加，就得到了入射頻率為ω時(shí)，作用在風(fēng)電安裝船的單位波幅規(guī)則波下的一階波浪力：

浮體在頻域范圍內(nèi)的一階運(yùn)動(dòng)方程為：

式中，M為風(fēng)電安裝船總的質(zhì)量矩陣；λ為附加質(zhì)量矩陣；C為系統(tǒng)線性阻尼矩陣；K為系統(tǒng)總剛度矩陣；F（w）為波浪載荷（單位波幅），w為波浪頻率；x為幅值響應(yīng)算子（RAO）。

基于上述理論，用水動(dòng)力分析程序AQWA對(duì)風(fēng)電安裝船進(jìn)行頻域分析，二階波浪力采用近場(chǎng)方法求得。其中，風(fēng)電安裝船橫搖和縱搖的一階波浪力傳遞函數(shù)RAO，如圖1所示，縱蕩和橫蕩的二階波浪力傳遞函數(shù)QTF如圖2所示。

圖1 滿(mǎn)載狀態(tài)橫搖和縱搖RAO結(jié)果

圖2 滿(mǎn)載狀態(tài)橫蕩和縱蕩QTF

1.2 時(shí)域運(yùn)動(dòng)響應(yīng)分析

頻域分析通常只適用于穩(wěn)態(tài)問(wèn)題，無(wú)法解決瞬態(tài)問(wèn)題和強(qiáng)非線性問(wèn)題。對(duì)于風(fēng)電安裝船在惡劣海況下的大幅值運(yùn)動(dòng)，須在時(shí)域內(nèi)模擬船體在波浪環(huán)境中的運(yùn)動(dòng)響應(yīng)。在時(shí)域內(nèi)建立風(fēng)電安裝船和系泊系統(tǒng)的運(yùn)動(dòng)控制方程，所受的外載荷主要包括波浪力和系泊力，運(yùn)動(dòng)耦合方程為：

式中，M為風(fēng)電安裝船質(zhì)量矩陣；λ為附加質(zhì)量矩陣；C為阻尼矩陣；K為剛度矩陣；Fwave為波浪載荷（包括高頻波浪力和低頻慢漂力）；Fmooring為系泊系統(tǒng)對(duì)風(fēng)電安裝船提供的系泊力［4］。

基于上述理論和海洋環(huán)境條件，用AQWA軟件的DRIFT模塊對(duì)風(fēng)電安裝船進(jìn)行時(shí)域運(yùn)動(dòng)響應(yīng)分析，計(jì)算出結(jié)構(gòu)總力，求得結(jié)構(gòu)加速度，進(jìn)而求得結(jié)構(gòu)最新位置，重復(fù)以上步驟求得模型的運(yùn)動(dòng)響應(yīng)。模型采用四根水平系泊線布置，系泊線剛度為2.78E6（N/m），1號(hào)和3號(hào)系泊線長(zhǎng)度為1 023 m，2號(hào)和4號(hào)系泊線長(zhǎng)度為1 062 m，計(jì)算模型如圖3所示。波浪為JONSWAP譜波浪，有義波高8.4 m，譜峰周期12.16 s，浪向角為135°和180°，時(shí)間步長(zhǎng)為0.12，步數(shù)為90 000步，每個(gè)工況模擬時(shí)間為3小時(shí)［5，6］。其中圖4為風(fēng)電安裝船135°浪向角時(shí)橫搖、縱搖和垂蕩時(shí)域運(yùn)動(dòng)響應(yīng)模擬結(jié)果。

圖3 時(shí)域運(yùn)動(dòng)響應(yīng)計(jì)算模型

2 模型試驗(yàn)研究

2.1 試驗(yàn)條件介紹

模型試驗(yàn)在哈爾濱工程大學(xué)多功能深水池和船模拖曳水池進(jìn)行。多功能深水池長(zhǎng)50 m，寬30 m，深10 m，水池配備75單位的搖板式多單元造波系統(tǒng)，在本次RAO試驗(yàn)和不規(guī)則波試驗(yàn)中需要用到。船模拖曳水池長(zhǎng)108 m，池寬7 m，水深3.5 m。在本次靜水拖航阻力試驗(yàn)中需要用到日本制造的四自由度適航儀。

圖4 時(shí)域運(yùn)動(dòng)響應(yīng)模擬結(jié)果

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

本次試驗(yàn)主要是為了得到風(fēng)電安裝船的相關(guān)水動(dòng)力特性，為后續(xù)設(shè)計(jì)與分析提供依據(jù)。需要開(kāi)展靜水拖航阻力試驗(yàn)、RAO試驗(yàn)和不規(guī)則波試驗(yàn)。分別得到風(fēng)電安裝船的阻力性能、RAO特征和實(shí)際海況中的運(yùn)動(dòng)特征，并與數(shù)值計(jì)算結(jié)果相比較。根據(jù)國(guó)際海洋工程模型試驗(yàn)相關(guān)規(guī)范，本次試驗(yàn)按照縮尺比1:50制作模型。模型用玻璃鋼材料由倒模技術(shù)制成，實(shí)船型線如圖5所示，模型如圖6所示。為了方便調(diào)整模型屬性參數(shù)，將模型沿甲板邊線增高350 mm，不影響模型原船體型線，原模型具體參數(shù)如表1所示。

圖5 風(fēng)電安裝船橫剖線圖

圖6 風(fēng)電安裝船模型

表1 模型屬性值

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

3.1 阻力試驗(yàn)結(jié)果

在哈爾濱工程大學(xué)船模拖曳水池開(kāi)展了風(fēng)電安裝船的拖航阻力試驗(yàn)［7］，圖7為進(jìn)行拖航阻力試驗(yàn)的過(guò)程。拖航航速由實(shí)船航速6 kn～13 kn換算得到，船模吃水分別為壓載和滿(mǎn)載狀態(tài)，航向角為180°，船模阻力由測(cè)力天平測(cè)得。采用傅汝德?lián)Q算法求得實(shí)船阻力，摩擦阻力由1957ITTC公式求得，其中滿(mǎn)載狀態(tài)阻力結(jié)果如表2所示。試驗(yàn)過(guò)程中發(fā)現(xiàn)，船舶首部與平行中體的過(guò)渡段附近有較大的興波，原因在于此處型線變化梯度大。實(shí)船阻力曲線如圖8所示，阻力曲線為二次拋物線形態(tài)，符合理論基礎(chǔ)［8］。

表2 風(fēng)電船滿(mǎn)載狀態(tài)阻力試驗(yàn)結(jié)果

圖7 拖航阻力試驗(yàn)

圖8 風(fēng)電安裝船阻力試驗(yàn)結(jié)果

3.2 RAO試驗(yàn)結(jié)果

模型在規(guī)則波中的試驗(yàn)，根據(jù)線性系統(tǒng)假設(shè)，輸入的單變量是波浪，記做x，一般由波浪譜密度函數(shù)Sx（ω）表示輸入波浪的大小。輸出的單變量是模型響應(yīng)（運(yùn)動(dòng)），記做y，一般由響應(yīng)譜密度函數(shù)Sy（ω）表示模型運(yùn)動(dòng)響應(yīng)大小，模型在單位波幅下的運(yùn)動(dòng)響應(yīng)幅值（RAO）記作Hj（ω）［9］，則：

圖9 風(fēng)電安裝船RAO試驗(yàn)過(guò)程

RAO試驗(yàn)過(guò)程如圖9所示。根據(jù)模型試驗(yàn)得到風(fēng)電安裝船在白噪聲波譜中的運(yùn)動(dòng)響應(yīng)，將時(shí)域信號(hào)進(jìn)行傅立葉變換，得到頻域信號(hào)，再進(jìn)行平滑處理得到船舶在不同浪向角下各自由度的RAO。限于篇幅有限，選擇兩種典型工況下的運(yùn)動(dòng)狀態(tài)對(duì)比結(jié)果進(jìn)行分析，其中風(fēng)電安裝船 90°浪向角橫蕩和橫搖RAO對(duì)比分析結(jié)果如圖10和圖11所示，180°浪向角縱蕩和縱搖RAO對(duì)比分析結(jié)果如圖12和圖13所示。從圖中可以看出，數(shù)值計(jì)算結(jié)果和試驗(yàn)結(jié)果比較吻合，驗(yàn)證了數(shù)值算法的準(zhǔn)確性，風(fēng)電安裝船滿(mǎn)載狀態(tài)下橫搖固有周期為13.5 s，縱搖固有周期為11 s。

圖10 滿(mǎn)載狀態(tài)90°橫蕩RAO對(duì)比分析

圖11 滿(mǎn)載狀態(tài)90°橫搖RAO對(duì)比分析

圖12 滿(mǎn)載狀態(tài)180°縱蕩RAO對(duì)比分析

3.3 不規(guī)則波試驗(yàn)結(jié)果

針對(duì)風(fēng)電安裝船滿(mǎn)載狀態(tài)，180°和135°浪向角兩種工況開(kāi)展不規(guī)則波試驗(yàn)。圖10為風(fēng)電安裝船滿(mǎn)載狀態(tài)135°浪向角時(shí)不規(guī)則波試驗(yàn)過(guò)程，系泊布置方式與數(shù)值計(jì)算的布置方式一致［10］。波浪為JONSWAP譜波浪，有義波高8.4 m，譜峰周期12.16 s。其中對(duì)風(fēng)電安裝船滿(mǎn)載狀態(tài) 135°浪向角工況下的運(yùn)動(dòng)響應(yīng)幅值進(jìn)行統(tǒng)計(jì)分析，換算為實(shí)船運(yùn)動(dòng)響應(yīng)，并與時(shí)域運(yùn)動(dòng)響應(yīng)數(shù)值計(jì)算結(jié)果進(jìn)行對(duì)比分析，所得結(jié)果如表3所示。由于數(shù)值計(jì)算不能準(zhǔn)確模擬船舶在惡劣海況下的強(qiáng)非線性運(yùn)動(dòng)，所以時(shí)域響應(yīng)計(jì)算與試驗(yàn)結(jié)果有一定的誤差。

圖13 滿(mǎn)載狀態(tài)180°縱搖RAO對(duì)比分析

圖14 風(fēng)電安裝船不規(guī)則波試驗(yàn)過(guò)程

表3 時(shí)域運(yùn)動(dòng)響應(yīng)對(duì)比分析結(jié)果

4 結(jié)論

在阻力試驗(yàn)中，船首與平行中體的過(guò)渡區(qū)域有較為明顯的型波，型波逐漸向船尾及兩側(cè)擴(kuò)散，說(shuō)明此處型線變化梯度大。通過(guò)RAO試驗(yàn)得到風(fēng)電安裝船滿(mǎn)載狀態(tài)下的橫搖和縱搖固有周期，并且驗(yàn)證了AQWA程序?qū)AO計(jì)算的可行性，風(fēng)電安裝船在波浪周期為 11 s～13 s的海況下作業(yè)非常危險(xiǎn)。隨著波浪周期的增大，風(fēng)電安裝船在單位波幅下的橫蕩和縱蕩運(yùn)動(dòng)趨近于波高。在不規(guī)則波試驗(yàn)中，風(fēng)電安裝船的縱橫搖角達(dá)到10°，船首波浪抨擊現(xiàn)象嚴(yán)重，并且船體有較大范圍的縱蕩運(yùn)動(dòng)。時(shí)域運(yùn)動(dòng)響應(yīng)計(jì)算結(jié)果與試驗(yàn)結(jié)果誤差較大，原因在于AQWA程序不能模擬自由液面的動(dòng)態(tài)變化，模型縱橫搖運(yùn)動(dòng)角度越大，模擬結(jié)果越不準(zhǔn)確，這是需要進(jìn)一步研究的問(wèn)題。本文的研究結(jié)果對(duì)風(fēng)電安裝船的設(shè)計(jì)開(kāi)發(fā)和實(shí)際作業(yè)具有指導(dǎo)意義。

［1］郝金鳳，強(qiáng)兆新，丁君海，等.風(fēng)電安裝船線型設(shè)計(jì)及性能分析［J］.船舶工程，2014（S1）: 226-230.

［2］張鳳偉.風(fēng)機(jī)安裝船海上作業(yè)運(yùn)動(dòng)及其流體動(dòng)力分析［D］.哈爾濱: 哈爾濱工程大學(xué)，2012.

［3］馬汝建.任意結(jié)構(gòu)形狀的大型海洋結(jié)構(gòu)物的附加質(zhì)量［J］.中國(guó)海洋平臺(tái)，1995，102（2）: 65-71.

［4］Zhang X，Zou Z.Estimation of the Hydrodynamic Coefficients from Captive Model Test Results by Using Support Vector machines［J］.Ocean Engineering，2013，73（4）: 25-31.

［5］趙亞楠.海上自升式安裝船環(huán)境載荷分析與樁腿驅(qū)動(dòng)控制研究［D］.哈爾濱工程大學(xué)，2012.

［6］張婷，賀捷，黃錦林.海洋平臺(tái)波浪荷載數(shù)值模擬研究［J］.船海工程，2013（5）: 150-154.

［7］周廣利，黃德波，李鳳來(lái).船模拖曳阻力試驗(yàn)的不確定度分析［J］.哈爾濱工程大學(xué)學(xué)報(bào)，2006（3）: 377-381，390.

［8］閆蕾，宗智，倪少玲，等.半滑行前三體船型模型阻力對(duì)比試驗(yàn)研究［J］.船海工程，2013（1）: 41-44.

［9］楊建民，肖龍飛，盛振邦.海洋工程水動(dòng)力學(xué)試驗(yàn)研究［M］.上海: 上海交通大學(xué)出版社，2008.

［10］劉建峰.自升式海上風(fēng)電安裝船的模型試驗(yàn)分析研究［J］.造船技術(shù)，2014（1）: 21-25.

中國(guó)船級(jí)社（CCS）將發(fā)布全球首部LNG運(yùn)輸船兼作浮式儲(chǔ)存裝置實(shí)施指南

中國(guó)船級(jí)社（CCS）制定的《液化天然氣運(yùn)輸船兼作浮式儲(chǔ)存裝置實(shí)施指南》順利通過(guò)了來(lái)自海事管理機(jī)構(gòu)、航運(yùn)企業(yè)、能源公司、設(shè)計(jì)單位、科研院所、造船廠、設(shè)備廠家等單位的30余位專(zhuān)家和代表的評(píng)審，于2016年7月正式發(fā)布實(shí)施。屆時(shí)，該指南將成為全球首部LNG運(yùn)輸船兼作浮式儲(chǔ)存裝置（FSU）實(shí)施指南（以下簡(jiǎn)稱(chēng)LNG-FSU指南）。

LNG-FSU指南主要包括附加標(biāo)志授予、適應(yīng)性條件評(píng)估、連接系統(tǒng)、通訊系統(tǒng)、應(yīng)急系統(tǒng)、系泊系統(tǒng)、兼容性評(píng)估、貨物安全操作及風(fēng)險(xiǎn)評(píng)估等方面的技術(shù)要求。

從技術(shù)要求的安全可靠性、可實(shí)現(xiàn)性、可操作性的角度出發(fā)，基于中國(guó)自主設(shè)計(jì)建造的首艘LNG運(yùn)輸船“海洋石油301”在印尼作為FSU作業(yè)的實(shí)例專(zhuān)題研究，CCS對(duì)LNG運(yùn)輸船兼作FSU作業(yè)進(jìn)行了全面分析，充分考慮了LNG運(yùn)輸船兼作FSU作業(yè)時(shí)的風(fēng)險(xiǎn)、作業(yè)的兼容性和貨物安全操作等問(wèn)題，以及船舶水動(dòng)力特性、FSU與補(bǔ)給船之間相對(duì)運(yùn)動(dòng)對(duì)連接系統(tǒng)和系泊系統(tǒng)帶來(lái)的影響，提出了相應(yīng)的風(fēng)險(xiǎn)評(píng)估內(nèi)容、評(píng)估結(jié)果、風(fēng)險(xiǎn)控制措施和技術(shù)要求。

該指南緊密聯(lián)系工業(yè)界的最新需求，緊跟全球發(fā)展熱點(diǎn)，理論與實(shí)踐相結(jié)合，對(duì)促進(jìn)和引導(dǎo)LNG浮式儲(chǔ)存裝置的發(fā)展具有重要意義。

來(lái)源：中國(guó)船級(jí)社

Motion Response Analysis and Model Test Study of Wind Turbine Installation Vessel

Xu Chao1，Sun Li-ping1，Zhao Zhi-jian2，Zhang Li-jun2

（1.Harbin Engineering University，Harbin 150001，China； 2.COSCO Shipyard Group Co.，Ltd.，Liaoning Dalian 116600，China）

In order to study the motion response and resistance performance for the wind turbine installation vessel （WTIV），based on the three-dimensional potential flow theory，the numerical analysis of the vessel’s motion response is conducted in frequency and time domain.The model test of the WTIV is also carried out，including resistance test，RAO test and irregular test.Through the test，the resistance of the WTIV at the speed of 13 kn is 899 kN.The RAO response curve and the motion value of the irregular test verify the accuracy of the numerical calculation.The results show that the WTIV has large motion value in the serious sea condition of ten years on the North Atlantic Ocean，which is not suitable for installation.

wind turbine installation vessel （WTIV）； model test； sea-keeping performance； resistance

U661.43

A

10.14141/j.31-1981.2016.04.007

工信部高技術(shù)船舶科研項(xiàng)目（工信部聯(lián)裝2011536號(hào)）。

許超（1990—），男，碩士研究生。研究方向：海洋工程結(jié)構(gòu)水動(dòng)力學(xué)。