成欣　葉舟　高月文　李春

摘 要：分別針對(duì)MIT/NREL TLP和UmaineHywind Spar兩種海上風(fēng)力機(jī)浮式平臺(tái)主體的水動(dòng)力特性進(jìn)行研究，旨在分析兩者的系泊穩(wěn)定性.基于海洋水動(dòng)力學(xué)和結(jié)構(gòu)動(dòng)力學(xué)理論，建立了平臺(tái)/纜索系統(tǒng)耦合模型，在水深和外界載荷激勵(lì)相同的情況下，利用有限元分析ANSYS軟件中的水動(dòng)力學(xué)計(jì)算模塊進(jìn)行時(shí)域、頻域響應(yīng)分析，研究了兩種平臺(tái)在海風(fēng)、海流和隨機(jī)波聯(lián)合作用下的動(dòng)態(tài)響應(yīng)，并分析了兩種平臺(tái)隨波浪頻率的響應(yīng)變化.結(jié)果表明：MIT/NREL TLP平臺(tái)的動(dòng)態(tài)響應(yīng)較大，而UmaineHywind Spar平臺(tái)動(dòng)態(tài)響應(yīng)較??；兩平臺(tái)均在低頻波浪作用下產(chǎn)生響應(yīng)峰值.

關(guān)鍵詞：海上風(fēng)力機(jī)； 漂浮式平臺(tái)； 時(shí)域； 頻域； 動(dòng)態(tài)特性

中圖分類(lèi)號(hào)： TK 83 文獻(xiàn)標(biāo)志碼： A

海上有豐富的風(fēng)能資源和廣闊平坦的區(qū)域，使得近海風(fēng)力發(fā)電技術(shù)成為近年來(lái)研究和應(yīng)用的熱點(diǎn).我國(guó)已建成了東海大橋海上風(fēng)電場(chǎng)，計(jì)劃建設(shè)的近海淺水區(qū)域風(fēng)電場(chǎng)有10多個(gè)[1].除渤海外，我國(guó)黃海、東海和南海平均水深均在40 m以上，淺水區(qū)域有限，非常適合漂浮式風(fēng)力機(jī)系統(tǒng)的開(kāi)發(fā).海上漂浮式風(fēng)力機(jī)的最大特點(diǎn)是克服了在海床底部安裝基礎(chǔ)結(jié)構(gòu)時(shí)受水深影響的缺點(diǎn)，使海上風(fēng)電場(chǎng)的建設(shè)可向深水區(qū)發(fā)展[2].在海上風(fēng)力機(jī)平臺(tái)和開(kāi)發(fā)深海油氣的平臺(tái)中，張力腿平臺(tái)（TLP）和浮柱式平臺(tái)（Spar）是公認(rèn)的優(yōu)秀平臺(tái)形式[3].TLP是一種典型的深水平臺(tái)，以其半固定、半順應(yīng)的運(yùn)動(dòng)特征在深水海洋工程中廣為應(yīng)用，其最重要的特點(diǎn)是平臺(tái)的豎向運(yùn)動(dòng)很小[4-5].Spar平臺(tái)屬于順應(yīng)式平臺(tái)的范疇，憑借其優(yōu)良的性能和相對(duì)較低的造價(jià)，成為世界深海風(fēng)力機(jī)和油氣開(kāi)采的主力平臺(tái)類(lèi)型之一[6].近年來(lái)，我國(guó)風(fēng)電產(chǎn)業(yè)發(fā)展迅猛，海上石油開(kāi)發(fā)業(yè)從近海走向深海[7]，為適應(yīng)深海平臺(tái)開(kāi)發(fā)的需要，有關(guān)深海平臺(tái)的相關(guān)理論和技術(shù)的研究得到了蓬勃發(fā)展.

浮式平臺(tái)的動(dòng)態(tài)特性能夠反映平臺(tái)在外界載荷作用下的響應(yīng)特性和穩(wěn)定性，對(duì)整個(gè)漂浮式風(fēng)力機(jī)系統(tǒng)的穩(wěn)定性和安全性有重大影響.本文分別選取設(shè)計(jì)水深相同的MIT/NREL TLP平臺(tái)和UmaineHywind Spar平臺(tái)，借助有限元分析ANSYS軟件中的水動(dòng)力計(jì)算模塊進(jìn)行時(shí)域、頻域響應(yīng)分析計(jì)算，分析兩種平臺(tái)在相同風(fēng)、浪、流作用下的位移動(dòng)態(tài)響應(yīng)以及平臺(tái)響應(yīng)隨波浪頻率的變化.

1 基本環(huán)境載荷

環(huán)境載荷指直接或間接由環(huán)境作用引起的載荷，包括由環(huán)境載荷引起的所有外力，如系泊力、運(yùn)動(dòng)慣性力、液艙晃蕩力等.漂浮式風(fēng)力機(jī)所受到的環(huán)境載荷主要來(lái)自風(fēng)、浪、流，還有內(nèi)波、地震和海冰（大塊浮冰或冰山）.本文僅考慮前三種載荷.

1.1 風(fēng)載荷

海風(fēng)對(duì)海洋結(jié)構(gòu)物的工作影響很大，通常作用在結(jié)構(gòu)上的風(fēng)載荷是動(dòng)態(tài)的，但一些結(jié)構(gòu)對(duì)風(fēng)載荷的反應(yīng)幾乎為靜態(tài)形式.在實(shí)際工程中常用設(shè)計(jì)風(fēng)速法得到結(jié)構(gòu)所受到的海風(fēng)載荷.根據(jù)我國(guó)海上移動(dòng)平臺(tái)入級(jí)規(guī)范可計(jì)算作用于構(gòu)件上的風(fēng)力F，并應(yīng)確定合力作用點(diǎn)的垂直高度[8]，即

式中：S為平臺(tái)在正浮或傾斜狀態(tài)時(shí)受風(fēng)構(gòu)件的正投影面積；Ch、Cs分別為受風(fēng)構(gòu)件的高度系數(shù)和形狀系數(shù)；P為風(fēng)壓；V為設(shè)計(jì)風(fēng)速.

1.2 波浪載荷

波浪載荷是漂浮式風(fēng)力機(jī)平臺(tái)所受到的主要環(huán)境載荷，且較風(fēng)載荷更為復(fù)雜，這也是與陸上風(fēng)力機(jī)最大不同之處.作用在漂浮式平臺(tái)上的水動(dòng)力包括波浪激振力、由結(jié)構(gòu)運(yùn)動(dòng)產(chǎn)生的輻射力和考慮波浪高階效應(yīng)的漂移力.其中：波浪激振力包括由入射波形成的不穩(wěn)定壓力所引起的F-K力（佛汝德-克雷洛夫力）和由于結(jié)構(gòu)存在影響波浪密度分布產(chǎn)生壓差而引起的繞射力；輻射力包括附加質(zhì)量力和輻射阻尼力.對(duì)于大尺度結(jié)構(gòu)，一般采用繞射理論對(duì)作用于物體表面上的整個(gè)水動(dòng)壓力進(jìn)行積分計(jì)算波浪載荷.在繞射理論中，流場(chǎng)用速度勢(shì)函數(shù)描述，該速度勢(shì)在流體各處應(yīng)滿足拉普拉斯方程，并滿足物體表面、自由表面、海底及無(wú)窮遠(yuǎn)處的邊界條件.通?？偹俣葎?shì)由入射勢(shì)、繞射勢(shì)（假定物體固定）和運(yùn)動(dòng)著的物體在靜水中產(chǎn)生的輻射勢(shì)組成，從而得到結(jié)構(gòu)物的運(yùn)動(dòng)與作用力.描述波浪運(yùn)動(dòng)的理論有許多，根據(jù)不同要素可分為線性和非線性、有旋和無(wú)旋、規(guī)則和不規(guī)則、單向和多向、淺水和深水等.漂浮式風(fēng)力機(jī)平臺(tái)一般遠(yuǎn)離海岸，局部水深幾乎不變；與波長(zhǎng)相比，水深相對(duì)較大.通常用無(wú)旋模型描述一個(gè)波浪在其生成區(qū)域外的傳播或在水池內(nèi)造波機(jī)產(chǎn)生的波浪傳播.理想流體和無(wú)旋流動(dòng)的假設(shè)可使問(wèn)題大大簡(jiǎn)化，其方法是引入速度勢(shì)函數(shù)Φ（x，y，z，t），利用速度勢(shì)給出速度場(chǎng)V（x，y，z，t）=

Φ（x，y，z，t），于是流動(dòng)的無(wú)旋性自動(dòng)得到保證.在流域內(nèi)質(zhì)量守恒（對(duì)于不可壓縮流體，divV=0）可表示為速度勢(shì)函數(shù)的拉普拉斯方程[9]，即

通過(guò)伯努利拉格朗日關(guān)系式可得到流域內(nèi)的壓力

式中： p0為參考?jí)毫Γ创髿鈮毫?；ρw為海水密度；g為重力加速度.

波浪理論的分類(lèi)與3個(gè)幾何參數(shù)（波高H、波長(zhǎng)L、水深h）有關(guān).淺水理論與深水（或中等水深）理論之間的區(qū)別取決于厄塞爾參數(shù)Ur，有

式中：A為波幅；k=2π/L；非線性參數(shù)ε=A/h；色散參數(shù)μ=kh.

當(dāng)Ur<1，意味著色散效應(yīng)較非線性效應(yīng)占優(yōu)勢(shì)，用參數(shù)kA≡ε對(duì)勢(shì)函數(shù)Φ進(jìn)行展開(kāi)，可實(shí)現(xiàn)對(duì)問(wèn)題的求解，即

式中，（i）為勢(shì)函數(shù)的i階導(dǎo)數(shù).

這就是所謂的斯托克斯方法，由此可得到斯托克斯規(guī)則波模型.相反，當(dāng)Ur>1，h/L為小參數(shù)，從而得到橢圓余弦波模型和孤立波模型.在大部分海洋工程模型中，水深足夠大，故可應(yīng)用斯托克斯模型[10].

小尺度海洋結(jié)構(gòu)上的波浪力Fw通常采用Morison方程計(jì)算，由拖曳力和慣性力組成，即

式中：CD、CM分別為曳力系數(shù)和慣性力系數(shù)；u為垂直于構(gòu)件軸線水質(zhì)點(diǎn)的速度分量；D為直立圓柱的直徑.

1.3 海流載荷

海流載荷是由海流作用在海洋工程結(jié)構(gòu)物上所產(chǎn)生的載荷.海流力是作用在海洋結(jié)構(gòu)物上的一種流動(dòng)阻力.根據(jù)水下結(jié)構(gòu)物上的阻力是流體動(dòng)能函數(shù)的原理，可按照穩(wěn)定流動(dòng)條件下阻力的數(shù)學(xué)表達(dá)式得到海流力，其計(jì)算需合理確定海平面以下某深度的海流速度、阻力系數(shù)和慣性力系數(shù).其中，海流設(shè)計(jì)流速應(yīng)取為在平臺(tái)作業(yè)海區(qū)范圍內(nèi)可能出現(xiàn)的最大流速值，包括潮流流速、風(fēng)暴涌流速和風(fēng)成流流速，還應(yīng)考慮作業(yè)海區(qū)流速的垂向分布.當(dāng)波浪存在時(shí)，還應(yīng)對(duì)無(wú)波浪時(shí)的流速垂向分布進(jìn)行修正，以使瞬時(shí)波面處的流速保持不變.

當(dāng)只考慮海流作用時(shí)，作用在平臺(tái)水下部分構(gòu)件的海流載荷為式中，Vw為設(shè)計(jì)海流流速.

應(yīng)注意海流與波浪的相互作用.當(dāng)采取Morison公式計(jì)算波浪載荷時(shí)，應(yīng)將波浪水質(zhì)點(diǎn)速度與海流速度矢量相加；當(dāng)采用繞射理論計(jì)算波浪載荷時(shí)，海流載荷應(yīng)按式（8）計(jì)算，并與波浪載荷矢量相加.

2 平臺(tái)建模與計(jì)算

2.1 平臺(tái)模型建立

本文研究對(duì)象為MIT/NREL TLP和UmaineHywind Spar平臺(tái).TLP平臺(tái)由麻省理工學(xué)院設(shè)計(jì)的張力腿平臺(tái)改進(jìn)得到，以混凝土為壓艙物，8根纜索分成4組構(gòu)成4根筋腱，分別連接4根從底部水平延伸而出的輻條和位于海底的樁基，纜索由于預(yù)張力作用始終處于完全張緊狀態(tài)；Spar平臺(tái)是緬因大學(xué)根據(jù)國(guó)際能源署（IEA）研發(fā)的OC3Hywind Spar改變水深以便于試驗(yàn)比較而得來(lái)，主體為一深吃水的細(xì)長(zhǎng)浮筒，周身附連3根懸鏈線纜索并通過(guò)三腳架連接，以增加平臺(tái)系泊的抗偏剛度，相鄰纜索夾角為120°，且在一定預(yù)張力作用下處于半張緊半松弛狀態(tài)[11].

本文采用ANSYS軟件中水動(dòng)力計(jì)算模塊進(jìn)行模型導(dǎo)入和表面網(wǎng)格劃分，分析中對(duì)系泊系統(tǒng)進(jìn)行了適當(dāng)簡(jiǎn)化，忽略了系泊系統(tǒng)的阻力.對(duì)UmaineHywind Spar作了去除三角形連接、纜索直接連在浮筒上的簡(jiǎn)化，這在靜態(tài)分析中是合適的，但不適用于所有的動(dòng)態(tài)條件[12].計(jì)算中沒(méi)有對(duì)纜索進(jìn)行線性簡(jiǎn)化.兩平臺(tái)參數(shù)如表1所示.圖1分別為MIT/NREL TLP和UmaineHywind Spar平臺(tái)幾何模型.

2.2 響應(yīng)算子與自由度

浮體波頻運(yùn)動(dòng)6個(gè)自由度響應(yīng)的穩(wěn)態(tài)幅值稱(chēng)為幅值運(yùn)動(dòng)響應(yīng)算子（response amplitude operator）.該算子是由波浪激勵(lì)到船體或浮體運(yùn)動(dòng)的傳遞函數(shù)，為船體或浮體運(yùn)動(dòng)譜與波浪譜的比值，表征單位波幅的特征響應(yīng).對(duì)于平動(dòng)，其幅值響應(yīng)算子Y2yζ（ω）=Syζ（ω）/Sζ（ω）；對(duì)于轉(zhuǎn)動(dòng)，其幅值響應(yīng)算子Y2θζ（ω）=Sθζ（ω）/Sζ（ω），其中：Syζ（ω）和Sθζ（ω）分別為海上結(jié)構(gòu)物的平動(dòng)和轉(zhuǎn)動(dòng)運(yùn)動(dòng)譜；Sζ（ω）為波浪譜；ω為入射波浪的頻率.

隨機(jī)波浪可視為無(wú)數(shù)個(gè)振幅、頻率不等，初相位隨機(jī)并沿與x軸成不同角度的方向傳播的簡(jiǎn)單余弦波的疊加.通常把風(fēng)浪和由此引起的搖蕩運(yùn)動(dòng)都看成是具備各態(tài)歷經(jīng)性的平穩(wěn)隨機(jī)過(guò)程，平臺(tái)對(duì)任一波浪成分的響應(yīng)是該成分波波幅的線性函數(shù)并與它對(duì)其它波浪成分的響應(yīng)無(wú)關(guān).利用平臺(tái)各自由度的運(yùn)動(dòng)幅值響應(yīng)算子給出在每一個(gè)波浪頻率下的平臺(tái)響應(yīng)并疊加求和，可得到在多個(gè)波浪作用下的平臺(tái)運(yùn)動(dòng)方程[13].可用平臺(tái)位置的幅值響應(yīng)算子反映在外界環(huán)境載荷激勵(lì)下平臺(tái)在6個(gè)自由度上的運(yùn)動(dòng)情況，從而體現(xiàn)平臺(tái)穩(wěn)定性.6個(gè)自由度上的運(yùn)動(dòng)如圖2所示.由于風(fēng)、浪、流均為-180°入射，因此主要分析沿x軸的平動(dòng)（縱蕩）和繞y軸的轉(zhuǎn)動(dòng)（縱搖）位置幅值響應(yīng)算子.

2.3 計(jì)算工況

在外界激勵(lì)條件相同時(shí)對(duì)兩種平臺(tái)進(jìn)行時(shí)域響應(yīng)分析，風(fēng)速恒定取為10 m·s-1，設(shè)定海流速度從海平面至海底由0.95 m·s-1線性減小至0 m·s-1，隨機(jī)波譜選擇P-M譜，外界海況條件參數(shù)如表2所示.利用表2的數(shù)據(jù)建立平臺(tái)與纜索在風(fēng)、浪、流聯(lián)合作用下的耦合計(jì)算模型，分析時(shí)間為300 s，時(shí)間步長(zhǎng)為0.01 s.

3 結(jié)果分析

3.1 MIT/NREL TLP平臺(tái)

3.1.1 時(shí)域響應(yīng)分析

圖3（a）、（b）分別為MIT/NREL TLP平臺(tái)縱蕩和縱搖時(shí)域位移幅值響應(yīng)算子，其中縱坐標(biāo)分別表示對(duì)應(yīng)單位波幅的位移和偏轉(zhuǎn)角.從圖3中可看出，平臺(tái)在風(fēng)、浪、流載荷作用下繞起始位置作往復(fù)運(yùn)動(dòng)，縱蕩和縱搖運(yùn)動(dòng)均在約180 s后開(kāi)始進(jìn)入峰值，其中縱蕩最大位置幅值響應(yīng)算子約為5 m·m-1，縱搖最大位置幅值響應(yīng)算子最大值約為17°·m-1.

3.1.2 頻域響應(yīng)分析

頻域響應(yīng)分析是研究平臺(tái)隨波浪頻率變化的運(yùn)動(dòng)響應(yīng)，考慮到載荷的作用方向與平臺(tái)自身特點(diǎn)，此處分析縱蕩、垂蕩和縱搖的頻域響應(yīng).圖4（a）、（b）、（c）分別為MIT/NREL TLP平臺(tái)縱

蕩、垂蕩和縱搖運(yùn)動(dòng)的頻域響應(yīng)算子.從圖中可看

出，平臺(tái)縱蕩運(yùn)動(dòng)在0.07 rad·s-1左右的波浪頻率下產(chǎn)生響應(yīng)峰值，之后隨著波浪頻率增大響應(yīng)急劇降低，約在0.4 ～0.7 rad·s-1之間上升，隨后便呈下降趨勢(shì)；垂蕩運(yùn)動(dòng)約在波浪頻率0.5 rad·s-1時(shí)

出現(xiàn)響應(yīng)峰值，該值前后均為急升急降趨勢(shì)；而縱搖運(yùn)動(dòng)的響應(yīng)峰值出現(xiàn)在波浪頻率0.8 rad·s-1左右；此外，TLP平臺(tái)的垂蕩和縱搖頻域響應(yīng)的峰值接近.

3.2 UmaineHywind Spar平臺(tái)

3.2.1 時(shí)域響應(yīng)分析

圖5為UmaineHywind Spar平臺(tái)在相同載荷作用下的縱蕩和縱搖位置幅值響應(yīng)算子，其中縱坐標(biāo)分別表示對(duì)應(yīng)單位波幅的位移和偏轉(zhuǎn)角.從圖中可看出，平臺(tái)在風(fēng)、浪、流載荷的作用下繞起始位置作往復(fù)運(yùn)動(dòng)，且比MIT/NREL TLP的運(yùn)動(dòng)更為規(guī)律.相比于TLP平臺(tái)，其縱蕩位移響應(yīng)算子均為cm級(jí)，平動(dòng)位移變化很?。挥捎谄脚_(tái)自身結(jié)構(gòu)特性，縱搖位置響應(yīng)亦很小，最大值不到3°·m-1.

3.2.2 頻域響應(yīng)分析

圖6為UmaineHywind Spar平臺(tái)的縱蕩、垂蕩和縱搖運(yùn)動(dòng)的頻域響應(yīng)算子.由圖6可知，

Spar平臺(tái)縱蕩運(yùn)動(dòng)在0.07 rad·s-1的波浪頻率下即出現(xiàn)響應(yīng)峰值，且除0.20～0.34 rad·s-1之間有小幅升高外，響應(yīng)幅值呈下降趨勢(shì)，且峰值略小于TLP平臺(tái)的峰值；垂蕩運(yùn)動(dòng)的響應(yīng)峰值出現(xiàn)在頻率小于0.5 rad·s-1處，峰值較TLP平臺(tái)的?。豢v搖運(yùn)動(dòng)的響應(yīng)峰值出現(xiàn)在頻率0.4 ～0.5 rad·s-1之間，

且與TLP平臺(tái)相比峰值較??；Spar平臺(tái)的垂蕩響應(yīng)峰值大于縱蕩和縱搖的響應(yīng)峰值.

以上分析說(shuō)明，在外界風(fēng)、浪、流載荷作用下：

（1） MIT/NREL TLP平臺(tái)產(chǎn)生較大動(dòng)態(tài)位移，并引發(fā)較大幅度的往復(fù)運(yùn)動(dòng)和搖擺運(yùn)動(dòng)；

（2） UmaineHywind Spar平臺(tái)動(dòng)態(tài)位移較小，考慮其重心遠(yuǎn)遠(yuǎn)低于浮心，穩(wěn)定性較好；

（3） 兩平臺(tái)的縱蕩、垂蕩和縱搖運(yùn)動(dòng)響應(yīng)對(duì)頻率變化均較敏感，并在低頻率時(shí)出現(xiàn)響應(yīng)峰值.

4 結(jié) 語(yǔ)

本文在給定海況條件下，考慮風(fēng)、浪、流聯(lián)合作用，借助有限元分析ANSYS軟件中的水動(dòng)力模塊計(jì)算分析了兩種經(jīng)典風(fēng)力機(jī)平臺(tái)的時(shí)域和頻域響應(yīng).當(dāng)然，由于在計(jì)算過(guò)程中尚未考慮平臺(tái)主體與風(fēng)力機(jī)正常運(yùn)行下的耦合作用，因此實(shí)際耦合情況下的穩(wěn)定性情況尚有待繼續(xù)分析.若能在后續(xù)研究工作中加以解決，并能同時(shí)考慮平臺(tái)與波浪之間的共振效應(yīng)，則可對(duì)這兩種經(jīng)典平臺(tái)的性能和實(shí)用性作出更具現(xiàn)實(shí)意義的對(duì)比，以期得到在不同海況條件下的最佳平臺(tái)選擇，或者通過(guò)對(duì)平臺(tái)的改進(jìn)以降低危險(xiǎn)和建造成本.這將具有更大的實(shí)用價(jià)值，且將為我國(guó)深海漂浮式風(fēng)力機(jī)的發(fā)展提供更多有益的參考.

參考文獻(xiàn)：

[1] 黃維平，劉建軍，趙戰(zhàn)華.海上風(fēng)電基礎(chǔ)結(jié)構(gòu)研究現(xiàn)狀及發(fā)展趨勢(shì)[J].海洋工程，2009，27（2）：130-134.

[2] 高偉，李春，劉全.深海漂浮式風(fēng)力機(jī)的概念設(shè)計(jì)與氣動(dòng)水動(dòng)力耦合特性評(píng)述[J].能源研究與信息，2011，27（3）：168-173.

[3] 唐友剛，張素俠，張若瑜，等.深海系泊系統(tǒng)動(dòng)力特性研究進(jìn)展[J].海洋工程，2008，26（1）：120-126.

[4] 胡志敏，董艷秋，張建民.張力腿平臺(tái)波浪載荷計(jì)算[J].中國(guó)海洋平臺(tái)，2002，17（3）：6-11.

[5] 董艷秋，胡志敏，張翼.張力腿平臺(tái)及基礎(chǔ)設(shè)計(jì)[J].海洋工程，2000，18（4）：63-68.

[6] 王興剛，孫昭晨，梁書(shū)秀，等.水深和螺旋板對(duì)Classic Spar水動(dòng)力特性的影響[J].中國(guó)海洋平臺(tái)，2011，26（1）：23-30.

[7] 曾曉輝，沈曉鵬，吳應(yīng)湘，等.深海平臺(tái)分析和設(shè)計(jì)中的關(guān)鍵力學(xué)問(wèn)題[J].船舶工程，2005，27（5）：18-21.

[8] 中國(guó)船級(jí)社.海上移動(dòng)平臺(tái)入級(jí)規(guī)范[M].北京：人民交通出版社，2012.

[9] 唐友剛.海洋工程結(jié)構(gòu)動(dòng)力學(xué)[M].天津：天津大學(xué)出版社，2008.

[10] MORLIN B.海洋工程水動(dòng)力學(xué)[M].劉水庚，譯.北京：國(guó)防工業(yè)出版社，2010.

[11] ROBERTSON A N，JONKMAN J M.Loads analysis of several offshore floating wind turbine concepts[C]∥International Society of Offshore and Polar Engineers 2011 Conference，Hawaii，2011.

[12] JONKMAN J.Definition of the floating system for phase IV of OC3[R].NREL/TP-500-47535，2010.

[13] 吳秀恒.船舶操縱性與耐波性[M].2版.北京：人民交通出版社，1999.