漂浮式風機時域耦合分析及結(jié)構(gòu)計算研究

摘要：近年來，漂浮式海上風機的發(fā)展受到了廣泛的關(guān)注。由于其結(jié)構(gòu)的特殊性，傳統(tǒng)的海上浮式結(jié)構(gòu)物設(shè)計方法無法完全滿足其整體系統(tǒng)的分析和支撐結(jié)構(gòu)的強度評估。本文對浮式風機設(shè)計中的關(guān)鍵技術(shù)進行了研究，首先，提出了一套包含風機、浮式基礎(chǔ)及系泊系統(tǒng)的全系統(tǒng)的耦合分析流程;在此基礎(chǔ)上，采取時域勢流理論及載荷傳遞方法，研究了支撐基礎(chǔ)結(jié)構(gòu)的強度計算和疲勞損傷的評估方法。

本文之后完成了某一漂浮式風機概念的算例研究，采用相應(yīng)的方法建立了水動力及結(jié)構(gòu)有限元模型，完成了頻域、時域水動力分析，系統(tǒng)耦合分析及結(jié)構(gòu)強度計算和疲勞校核。對這一方法的使用，本文總結(jié)了相關(guān)的經(jīng)驗，并對后續(xù)進一步的研究提出了展望。

關(guān)鍵詞：漂浮式風機，時域耦合分析，結(jié)構(gòu)設(shè)計

介紹

風能的開發(fā)在可再生能源中的應(yīng)用十分重要，世界各國都在大力發(fā)展包含風電技術(shù)在內(nèi)的各種新能源技術(shù)。相較于陸上風電開發(fā)，海洋風能資源較陸地更為豐富穩(wěn)定，適合大兆瓦風機的開發(fā)和利用。目前新開發(fā)的海上固定式風機已普遍采用6兆瓦以上機型，而10兆瓦甚至12兆瓦的大型風機也已完成開發(fā)和生產(chǎn)組裝。大型風機的使用，將進一步降低風電開發(fā)和運營成本。經(jīng)過30年的發(fā)展與探索，隨著海洋風電技術(shù)日益進步，海洋風電開發(fā)已經(jīng)成為當下新能源技術(shù)開發(fā)和利用的主要形式之一，具有廣闊的應(yīng)用前景。

當水深超過50-60米后，傳統(tǒng)的固定式基礎(chǔ)在環(huán)境載荷及結(jié)構(gòu)設(shè)計、安裝施工等方面都會面臨較大的挑戰(zhàn)。而另一方面，深遠海離岸較遠，風場較淺海海域，風速更強且更穩(wěn)定，加之深遠海海域空間廣闊，風能儲藏量極大，特別適合安裝大型風機機組。如能長期、大規(guī)?；乩么笮惋L機進行海洋風力發(fā)電，將能有效降低海洋風力發(fā)電的度電成本。因此，適用于深遠海并能夠支撐大型化風機的浮式基礎(chǔ)結(jié)構(gòu)被廣泛的研究。從安裝的角度講，浮式風機可以在近岸碼頭進行組裝，從而免去海洋固定式風機耗時、昂貴的海上起重、吊裝作業(yè)。浮式風機對水深條件的適用性良好，因此能夠?qū)崿F(xiàn)單一浮式風機概念的批量化設(shè)計與生產(chǎn)，進一步降低浮式風機的開發(fā)成本。

在過去的十幾年間，多座試驗性的非全尺寸樣機在歐美和亞洲的部分國家安裝運行，取得了大量有價值的設(shè)計和運行數(shù)據(jù)。目前國際上建成或者在建的主要示范商業(yè)項目包括蘇格蘭的Hywind項目和葡萄牙的WindFloat項目，如圖1和圖2所示。

我國國內(nèi)目前有一個示范項目（三峽陽江示范項目）已完成浮式基礎(chǔ)和風機建造下線，正等待運輸安裝，后續(xù)還有2-3個示范項目正在開展。

從技術(shù)層面上，浮式風機的基礎(chǔ)平臺主要參考浮式油氣平臺的設(shè)計，但是兩者的安全級別和用途具有差異性，因此浮式風機基礎(chǔ)仍具有較大的優(yōu)化空間。浮式風機的控制系統(tǒng)由于平臺運動與氣動載荷耦合導致的“氣動負阻尼”問題，也將挑戰(zhàn)相關(guān)控制系統(tǒng)的開發(fā)。由于水平氣動推力的作用，浮式風機將不可避免產(chǎn)生平均位移偏置和傾斜，使得系泊系統(tǒng)產(chǎn)生不均衡的張力狀態(tài)，這也有別于傳統(tǒng)油氣工業(yè)的浮式系統(tǒng)。

浮式風機設(shè)計關(guān)鍵技術(shù)研究

本文主要研究浮式風機設(shè)計中的兩個關(guān)鍵技術(shù)，一是包含風機（葉片、控制系統(tǒng)等）、浮式支撐結(jié)構(gòu)和系泊系統(tǒng)的全系統(tǒng)耦合計算;二是采用時域載荷傳遞的結(jié)構(gòu)分析。下面對這些關(guān)鍵技術(shù)及最終的全系統(tǒng)分析進行說明。

1. 時域耦合計算

對于傳統(tǒng)的陸上或者固定式海上風機，支撐結(jié)構(gòu)能提供相對剛性的支撐，“超單元”（Superelement）分析方法可以將基礎(chǔ)結(jié)構(gòu)以一組矩陣的形式傳遞給交接面/點（Interface Point），能夠在非全耦合的情況下對風機控制系統(tǒng)和基礎(chǔ)變形進行較好的全局考慮。而對于漂浮式風機，初步研究表明，用于固定式基礎(chǔ)的風機機型在漂浮式風機概念上的應(yīng)用中，主要應(yīng)考慮塔筒、尤其是控制系統(tǒng)的優(yōu)化。然而，與固定式基礎(chǔ)相比，浮式基礎(chǔ)具有較大的剛體位移（與之相比其自身變形成為小量），而風機葉片對于這種運動，例如縱傾十分敏感，葉片轉(zhuǎn)動面5度的縱傾或能帶來約20%的受力變化。另一方面，浮式基礎(chǔ)平臺的系泊定位系統(tǒng)本身也是一個非線性系統(tǒng)，如果考慮相對淺水或者組合錨鏈形式、非金屬材料等的使用等，系統(tǒng)非線性進一步增強。因而整個浮式基礎(chǔ)和系泊系統(tǒng)無法像固定式基礎(chǔ)一樣作為一組矩陣來和風機一起計算。

基于上述挑戰(zhàn)，本文認為，從系統(tǒng)分析的角度講，可以考慮將空氣動力學、浮式結(jié)構(gòu)物勢流理論和柔性桿件模型統(tǒng)一放在一個系統(tǒng)下，分別模擬風機、浮式基礎(chǔ)平臺和系泊系統(tǒng)，采用全時域耦合分析的方法 [3]。

其中風機葉片的載荷采用的是葉片元素動量（Blade Element Momentum，BEM）理論。該理論結(jié)合了動量理論和葉片元理論，基于經(jīng)驗升力和阻力系數(shù)，在機翼上局部產(chǎn)生的力與流過轉(zhuǎn)子盤的空氣動量的變化相平衡[4]。

對浮式支撐結(jié)構(gòu)而言，首先通過頻域勢流理論得到其水動力系數(shù)，如質(zhì)量/剛度矩陣、附加質(zhì)量和勢流阻尼系數(shù)、一階及二階波浪力等。在時域求解內(nèi)，波浪力可以表達為以下形式[5]

其中 和 分別為一階和二階波浪力， 為三階（或更高階）波浪力， 和 則為一階和二階脈沖響應(yīng)函數(shù)。

在本文的耦合系統(tǒng)分析中，系泊系統(tǒng)（及其它柔性結(jié)構(gòu)物）采用空間離散有限元系統(tǒng)模型，其動態(tài)平衡一般可以表示為[6]

其中 和 分別為慣性力向量、阻尼力向量、結(jié)構(gòu)內(nèi)響應(yīng)力向量和外部載荷向量，而 和 分別代表結(jié)構(gòu)運動、速度和加速度。

在分析過程中，上述的風機模型、浮體水動力模型和系泊系統(tǒng)模型通過超節(jié)點（Supernode）進行連接，從而實現(xiàn)在每一時間步長的實時耦合求解。

2. 時域載荷傳遞及結(jié)構(gòu)分析

浮式平臺基礎(chǔ)的結(jié)構(gòu)設(shè)計是浮式風機設(shè)計中的另一關(guān)鍵技術(shù)。其相關(guān)規(guī)范通常以滿足風機運營的設(shè)計標準為前提，針對不同結(jié)構(gòu)物形式，以傳統(tǒng)海工結(jié)構(gòu)物設(shè)計標準為基礎(chǔ)制定相應(yīng)的參考和設(shè)計準則，同時兼顧漂浮式風機的特點，提出了很多新的有針對性的要求。

從結(jié)構(gòu)分析的方法來看，傳統(tǒng)的船舶和海洋平臺結(jié)構(gòu)設(shè)計遵循的通常是頻域設(shè)計波法，即首先在頻域水動力載荷計算及長期統(tǒng)計預報的基礎(chǔ)上，根據(jù)不同平臺形式設(shè)計的結(jié)構(gòu)特點和規(guī)范要求，確定極限環(huán)境載荷工況（設(shè)計波），之后針對確定的設(shè)計工況進行結(jié)構(gòu)設(shè)計。疲勞分析也通常是基于頻域載荷在長期波浪分布條件下的累計損傷來確定。

然而，對于浮式風機設(shè)計而言，首先風機設(shè)計規(guī)范對設(shè)計工況有十分明確的定義，其計算通常要求在時域內(nèi)完成。其次，即使是“借鑒”油氣平臺設(shè)計概念，浮式風機基礎(chǔ)形式的結(jié)構(gòu)特點并無充分的歷史經(jīng)驗和參考，因而無法明確地找到所謂“典型截面”進行統(tǒng)計學預報。進一步的，在疲勞計算中，因為疲勞損傷同載荷及結(jié)構(gòu)應(yīng)力并不是線性關(guān)系，不同載荷過程，例如風機載荷、風浪流載荷以及錨鏈載荷等，產(chǎn)生的疲勞損傷并不能簡單的分別計算再加以疊加。因而風機基礎(chǔ)結(jié)構(gòu)設(shè)計必須遵從時域耦合計算和載荷傳遞的方法，綜合考慮不同載荷過程對平臺結(jié)構(gòu)產(chǎn)生的影響。

本文在結(jié)構(gòu)計算中采用了時域勢流理論和載荷傳遞[7]，在每一時間步長下滿足如下邊界條件：

其中， 為定常流（Basic Flow）速度勢， 和 分別為Local Flow，Memory Flow和入射波（Incoming wave）速度勢。 和 分別是Memory flow和入射波的自由液面升高。

結(jié)構(gòu)的時域疲勞評估采用雨流法[8]。

3. 全系統(tǒng)的耦合分析和結(jié)構(gòu)評估

對于整個系統(tǒng)的設(shè)計而言，圖3顯示了一種時域耦合分析、載荷傳遞及結(jié)構(gòu)疲勞和強度校核的流程[3]。其主要步驟簡述如下：

· 建立浮式基礎(chǔ)的水動力和結(jié)構(gòu)模型;

· 采用頻域勢流理論計算浮式基礎(chǔ)的水動力系數(shù)，如附加質(zhì)量和阻尼系數(shù)、一階和二階波浪力傳遞函數(shù)等;

· 將上述系數(shù)導入耦合分析系統(tǒng)，建立風機和系泊系統(tǒng)模型，進行全系統(tǒng)的時域耦合分析;

· 將風機運動、系泊系統(tǒng)及塔筒受力等時歷輸出，在時域勢流理論求解器中重現(xiàn)相關(guān)波浪及運動，進行時域載荷傳遞;

· 利用上述載荷，在結(jié)構(gòu)求解器中進行時域分析，之后進行強度和疲勞評估。

漂浮式風機分析實例

本文分析的對象基于工業(yè)聯(lián)合項目OC5（Offshore Code Comparison Collaboration，Continuation）中提供的三柱半潛式平臺DeepCWind概念，其主要尺寸如圖4和表1所示[9]。其結(jié)構(gòu)設(shè)計主要參考[10]。

圖5和6顯示的是浮式基礎(chǔ)的概念模型和艙室模型。這兩個模型用于生成頻域水動力計算的面元模型和載荷傳遞及分析的結(jié)構(gòu)模型，分別如圖7和8所示。在頻域分析中，基于勢流理論的計算僅需要面元單元，因而圖7中無Morison桿件模型。

圖9和10顯示的是葉片的某一單元的截面幾何翼型參數(shù)及其升阻力系數(shù)。

完整的風機和塔筒模型如圖11所示。本次計算中，采用3根懸鏈線形式的系泊纜，圖12顯示的是整體耦合分析模型。

圖13顯示的是時域勢流理論分析模型，注意到這里包含自由液面的有限元模型，用于求解滿足方程（3）中的邊界條件。

在結(jié)構(gòu)疲勞分析的過程中，對于關(guān)鍵結(jié)構(gòu)位置，例如支撐桿件同立柱的連接處，桿單元模型無法得到詳細的應(yīng)力分布，因而通常要使用局部細化子模型。圖14顯示的經(jīng)細化后的連接處的有限元模型。

頻域水動力分析及時域耦合計算結(jié)果

圖15-17顯示的浮式基礎(chǔ)的垂蕩、縱搖和橫搖在0°、15°、30°和45°浪向下的運動傳遞函數(shù)（Response Amplitude Operator，RAO），其中三個方向上分別添加了5%和3%的臨界阻尼系數(shù)。在不考慮風機及系泊系統(tǒng)的情況下，三個方向上的運動都避開了常見的波浪周期，幅值也在合理的范圍內(nèi)。由于前述的耦合因素的影響，這些RAO結(jié)果僅能作為初步設(shè)計的參考，接下來需進行時域耦合分析。

圖18-20顯示的一階波浪力在x和z方向上的RAO，以及x方向上的平均波浪力RAO。這些參數(shù)將用于時域耦合計算。

在本文的時域耦合分析中，波浪和風的方向設(shè)定為0°，采用3參數(shù)JONSWAP譜，其有義波高和譜峰周期分別為6米和10秒。定常風速為8米/秒。

圖21-23顯示的是縱蕩、垂蕩和縱搖三個方向上運動的（部分）時歷曲線?？梢钥吹剑v蕩方向上主要為低頻主導，垂蕩和縱搖方向則主要為波頻運動。

圖24顯示的是三根錨鏈張力的時歷曲線，其中迎浪方向上的錨鏈受力較大。三根錨鏈的受力均由低頻力占主導，符合平臺水平運動特性。

進一步的，通過傅里葉變換，得到縱蕩、垂蕩和系泊纜受力的能量譜密度，如圖26-28所示。圖中清楚地顯示了波頻響應(yīng)和低頻響應(yīng)的范圍和大小。這些分析進一步驗證了時域耦合分析的必要性。

浮式風機結(jié)構(gòu)分析

通過時域勢流理論計算，得到圖29顯示的（某一時刻）的載荷分布，可以看到計算中同時考慮了面元載荷和支撐桿件上受到的Morison力。

通過時域結(jié)構(gòu)分析，圖30和31顯示的是某一時刻面單元上的VonMises應(yīng)力分布和桿件單元上軸向力分布圖。

圖32和圖33顯示的是局部細化模型上的應(yīng)力分布，以及（部分）疲勞分析的結(jié)果。

結(jié)論和展望

本文首先總結(jié)了漂浮式風機系統(tǒng)分析的相關(guān)理論，包含葉片元素動量理論、頻域和時域勢流理論及耦合系統(tǒng)分析理論等，研究了整體系統(tǒng)及支撐結(jié)構(gòu)設(shè)計中的關(guān)鍵技術(shù)和挑戰(zhàn)，提出了全耦合時域分析及載荷傳遞和結(jié)構(gòu)分析的流程。之后基于DeepCwind概念，完成了上述流程的算例計算研究。首先建立了頻域分析面元模型，評估了該結(jié)構(gòu)的水動力特性，并得到了一階和二階波浪力等水動力系數(shù)。在時域分析系統(tǒng)中，建立了風機、浮式基礎(chǔ)和系泊系統(tǒng)的耦合模型，初步分析了平臺運動響應(yīng)、系泊系統(tǒng)受力等。在時域載荷傳遞分析的基礎(chǔ)上，針對某一海況計算了支撐浮體的結(jié)構(gòu)響應(yīng)，建立了局部細化模型并校核了疲勞損傷。本文的研究與嘗試，對于漂浮式風機這一新興浮式結(jié)構(gòu)物的設(shè)計具有重要的意義。

通過本文的研究與嘗試發(fā)現(xiàn)，采用新的全時域耦合分析和載荷傳遞，對傳統(tǒng)的分析方法和分析工具都帶來了巨大的挑戰(zhàn)。例如，由于風載荷的頻率較高，時域分析的時間步長通常較小，而平臺設(shè)計的不規(guī)則波浪時長較長，從而造成巨大的計算量和數(shù)據(jù)儲存量，加之風機設(shè)計規(guī)范中的工況組合要求較多，使得整個分析和結(jié)果評估過程變得非常耗時。而結(jié)構(gòu)設(shè)計通常是個循環(huán)的過程，如果每一輪的設(shè)計都實現(xiàn)這一完整的時歷耦合計算分析的話，將大大延長整個設(shè)計周期。后續(xù)將進一步研究在不同設(shè)計階段，采用不同設(shè)計方法相結(jié)合的可能性，同時考慮引入云計算，希望通過優(yōu)化設(shè)計過程和不同設(shè)計階段的數(shù)據(jù)接口，大幅提高設(shè)計效率。

參考文獻：

[1]https：//www.offshorewind.biz/2017/08/31/siem-moxie-starts-hywind-scotland-charter/

[2]https：//www.windcrete.com/windfloat-atlantic-has-begun-producing-power/

[3]張帆 2020. Recent developments & considerations of FOWT 大連理工大學 海岸和近海工程國家重點實驗室 學術(shù)講堂 http：//slcoe.dlut.edu.cn/info/1091/3472.htm

[4]Martin O.L. Hansen. 2008. Aerodynamics of Wind Turbines. Second Edition.

[5]Sintef Ocean. 2021. SIMO Theory Manual.

[6]Sintef Ocean. 2021. Riflex Theory Manual.

[7]DNV. 2021. Wasim User Manual

[8]DNV. 2021. Stofat User Manual.

[9]2017. OC5 Project Phase II：Validation of global loads of the DeepCwind floating Semisubmersible Wind Turbine.

[10]DNV. 2021. Sesam for Floating Offshore Wind design.

挪威船級社（中國）有限公司