趙志新，李 昕，施 偉，王文華

(1.大連理工大學(xué)海岸與近海工程國家重點實驗室，遼寧大連116024；2.大連理工大學(xué)水利工程學(xué)院，遼寧大連116024；3. 大連理工大學(xué)深海工程研究中心，遼寧大連116024；4. 大連理工大學(xué)海洋可再生能源研究中心，遼寧大連116024)

0 引 言

我國海上風(fēng)能資源儲量為陸上的3倍，且海上風(fēng)能具有風(fēng)速高、風(fēng)況穩(wěn)定、選址不受空間限制等優(yōu)勢，使得浮式海上風(fēng)力機(jī)成為開發(fā)利用風(fēng)能的必然選擇。然而，浮式海上風(fēng)力機(jī)是剛?cè)峄旌蠌?qiáng)非線性的多體系統(tǒng)，其動力特性分析異常復(fù)雜。Jonkman等[1]開發(fā)了嵌入FAST的水動力學(xué)計算模塊HydroDyn，建立了浮式風(fēng)力機(jī)氣動-水動-伺服-彈性全耦合計算模型，對5 MW駁船式浮式風(fēng)力機(jī)風(fēng)浪組合工況下的動力特性進(jìn)行了分析。Bachynski等[2]利用非線性氣動-水動耦合程序Simo-Riflex-AeroDyn，分析了故障工況下TLP、Spar和半潛三種典型浮式風(fēng)力機(jī)的動力特性。Karimirad等[3]設(shè)計了張力腿與Spar相結(jié)合的新型浮式風(fēng)力機(jī)支撐平臺，并使用HAWC2和USFOS/vpOne對其動力特性進(jìn)行了對比分析。馬鈺等[4]使用FAST對OC3-Hywind Spar式浮式風(fēng)力機(jī)在不同海況下的動力特性進(jìn)行分析，并研究了風(fēng)力機(jī)系統(tǒng)各運動模態(tài)之間的耦合效應(yīng)。

考慮到深海風(fēng)能開發(fā)的成本和效益，浮式海上風(fēng)力機(jī)大型化是必然趨勢。目前丹麥科技大學(xué)和維斯塔斯風(fēng)力技術(shù)公司聯(lián)合提出的DTU 10 MW風(fēng)力機(jī)[5]成為國內(nèi)外學(xué)者廣泛接受的樣本風(fēng)力機(jī)。Islam[6]、Xue[7]和Tian[8]采用DTU 10 MW風(fēng)力機(jī)，分別針對半潛、Spar式和TLP 3種典型浮式平臺結(jié)構(gòu)，運用簡化耦合時域分析方法，對比研究了3種典型浮式風(fēng)力機(jī)在工作和極端海況下的動力特性。然而，浮式風(fēng)力機(jī)大型化后，除面臨更為復(fù)雜的海洋環(huán)境荷載外，其巨型柔性葉片的氣彈性效應(yīng)、下部支撐平臺的大幅、強(qiáng)非線性運動響應(yīng)以及與系泊系統(tǒng)之間的強(qiáng)非線性耦合作用還有待深入研究。

以DTU 10 MW風(fēng)力機(jī)為研究對象，設(shè)計了無撐桿的半潛浮式風(fēng)力機(jī)支撐平臺，基于“氣動-水動-控制-彈性”全耦合計算模型，使用FAST軟件對超大型浮式風(fēng)力機(jī)系統(tǒng)在典型海況下的動力特性進(jìn)行了分析，并與NREL 5 MW無撐桿的半潛浮式風(fēng)力機(jī)的動力特性進(jìn)行了對比分析。

1 浮式海上風(fēng)力機(jī)結(jié)構(gòu)系統(tǒng)設(shè)計

以DTU 10 MW風(fēng)力機(jī)為研究對象，參照Luan等[9]提出的NREL 5 MW無撐桿的半潛風(fēng)力機(jī)支撐平臺，基于海洋工程普遍采用的Froude數(shù)相似準(zhǔn)則進(jìn)行放大[10]，完成了DTU 10 MW浮式風(fēng)力機(jī)支撐平臺的設(shè)計。半潛浮式風(fēng)力機(jī)整體設(shè)計方案如圖1所示，系泊系統(tǒng)布置如圖2所示，其中DTU 10 MW風(fēng)力機(jī)轉(zhuǎn)子為逆風(fēng)向、3葉片；控制方式為變速變槳調(diào)節(jié)；切入、額定、切出風(fēng)速分別為4、11.4、25 m·s-1；切入、額定轉(zhuǎn)速分別為6、9.6 r·min-1；轉(zhuǎn)子、輪轂直徑分別為178.3、5.6 m；輪轂中心高度為119 m；額定葉尖速比為90 r·min-1；轉(zhuǎn)子、機(jī)艙、塔架的質(zhì)量分別為227.962、446.036、527.362 t，更為詳細(xì)的參數(shù)請參考文獻(xiàn)[11]。平臺設(shè)計水深為100 m；平臺吃水為36.9 m；中心柱、邊柱直徑為8 m；中心柱和邊柱高度分別為41.82、54.12 m；底部浮筒的長、寬、高分別為55.965、11.07、7.38 m；平臺重心位置為(0 m，0 m，-30.092 m)； 平臺正浮時排水體積為19 257.13 m3； 平臺質(zhì)量(包括壓艙水)為17 942.21 t。此外，系泊系統(tǒng)選用懸鏈線式錨鏈；錨鏈數(shù)目為3 根，且相鄰錨鏈之間的夾角為120°。

圖1 浮式風(fēng)力機(jī)整體輪廓

圖2 系泊系統(tǒng)布置

2 浮式風(fēng)力機(jī)系統(tǒng)運動分析理論

2.1 系統(tǒng)運動方程

基于風(fēng)力機(jī)空氣動力學(xué)、水動力學(xué)及多體動力學(xué)基本理論，浮式風(fēng)力機(jī)系統(tǒng)時域耦合運動方程為：

(1)

2.2 波浪荷載

使用ANSYS-AQWA軟件，基于三維勢流理論計算浮式風(fēng)力機(jī)支撐平臺的水動力系數(shù)和波浪激勵力，以生成FAST-HydroDyn模塊的水動力輸入文件；并在浮式海上風(fēng)力機(jī)的水動力計算過程中，依據(jù)Morison公式的二次阻尼項考慮水體粘性，其中ANSYS-AQWA所建立平臺的面源模型如圖3所示。

圖3 ANSYS-AQWA水動力計算模型

圖4 縱蕩衰減運動時程曲線

2.3 風(fēng)荷載

使用FAST-Aerodyn模塊計算風(fēng)力機(jī)的氣動荷載，其計算方法基于經(jīng)典的葉素動量理論(BEM)。依據(jù)葉素理論，單個葉素上的氣動荷載可由下式得出：

dQ=dLcosφ+dDsinφ=

(2)

dT=dLsinφ-dDcosφ=

(3)

式中，dL和dT為葉素上的氣動升力和氣動阻力；CL和CD為葉素翼型的升力系數(shù)和阻力系數(shù)；ρa(bǔ)為空氣密度；c為葉素弦長；W為相對入流風(fēng)速；r為葉素到輪轂中心的距離；φ為入流角。

2.4 錨鏈張力

使用FAST-MoorDyn模塊來模擬系泊線，以考慮錨鏈慣性力和阻尼力等動態(tài)效應(yīng)，其主要基于梁理論，將錨鏈抽象為離散的動力學(xué)模型，并根據(jù)動力學(xué)平衡條件在時域內(nèi)對錨鏈的運動進(jìn)行模型[12]。

3 環(huán)境條件

參考IEC 61400—3規(guī)范設(shè)定的典型工況如表1所示。主要針對海上浮式風(fēng)力機(jī)所受的風(fēng)和波浪荷載來研究10 MW浮式海上風(fēng)力機(jī)在不同海洋環(huán)境條件下的動力特性，其中穩(wěn)態(tài)風(fēng)速為風(fēng)力機(jī)額定風(fēng)速11.4 m·s-1，波浪譜選取JONSWAP譜，且風(fēng)浪作用同向，如圖2所示。

表1 環(huán)境參數(shù)

4 計算結(jié)果分析

4.1 單自由度自由衰減運動結(jié)果

對支撐平臺剛體運動6自由度以及塔架順風(fēng)向和側(cè)向兩個方向的自由衰減運動分別進(jìn)行了數(shù)值模擬，得到縱蕩衰減運動時程曲線如圖4所示，并通過傅里葉變換計算得到支撐平臺縱蕩(橫蕩)、垂蕩、橫搖(縱搖)及艏搖固有頻率分別為0.013、0.036、0.039、0.011 Hz；塔架順風(fēng)向和側(cè)向固有頻率分別為0.368 Hz和0.379 Hz。

4.2 風(fēng)力機(jī)運動響應(yīng)分析

基于浮式風(fēng)力機(jī)“氣動-水動-控制-彈性”全耦合計算模型，依據(jù)表2中的環(huán)境參數(shù)，使用FAST對DTU10 MW浮式風(fēng)力機(jī)系統(tǒng)進(jìn)行全耦合分析，計算總時長為4 200 s。為消除在數(shù)值計算初始階段瞬態(tài)反應(yīng)對計算結(jié)果的影響，僅截取3 200 s到4 200 s時段的穩(wěn)態(tài)反應(yīng)進(jìn)行分析。因縱蕩、垂蕩、縱搖運動對半潛浮式風(fēng)力機(jī)而言較為關(guān)鍵，故本文僅研究縱蕩、垂蕩、縱搖運動響應(yīng)。

表2 不同工況下計算結(jié)果統(tǒng)計

比較表2中浮式風(fēng)力機(jī)在3種工況下的縱蕩運動統(tǒng)計值：在單獨風(fēng)作用時，風(fēng)力機(jī)系統(tǒng)運動達(dá)到準(zhǔn)平衡狀態(tài)后，縱蕩運動沿縱蕩方向平衡在7.75 m處；在單獨波浪作用時，風(fēng)力機(jī)系統(tǒng)在初始位置以小幅度做振蕩運動；在風(fēng)浪聯(lián)合作用下，風(fēng)力機(jī)系統(tǒng)在縱蕩運動平均位置7.76 m處振蕩。故對于浮式風(fēng)力機(jī)而言，風(fēng)荷載決定縱蕩運動的平衡位置，而波浪荷載主要引起縱蕩運動的振蕩，決定縱蕩運動的幅值。同樣比較風(fēng)力機(jī)系統(tǒng)垂蕩和縱搖運動統(tǒng)計值也具有與縱蕩運動相同的規(guī)律。

由圖5縱蕩、垂蕩和縱搖運動響應(yīng)頻譜圖可知：在風(fēng)浪聯(lián)合作用下，平臺縱蕩和縱搖運動響在縱蕩固有頻率0.013 Hz處和縱搖固有頻率0.039 Hz處均出現(xiàn)較大峰值，而平臺垂蕩運動在垂蕩固有頻率0.036 Hz處也出現(xiàn)較大峰值，這表明風(fēng)力機(jī)系統(tǒng)的運動包括風(fēng)荷載激勵產(chǎn)生的低頻共振響應(yīng)以及波浪荷載激勵產(chǎn)生的波頻運動響應(yīng)，且縱蕩和縱搖運動之間存在耦合效應(yīng)。此外，單獨波浪荷載激勵產(chǎn)生的平臺波頻運動反應(yīng)幅值顯著高于風(fēng)浪聯(lián)合作用下平臺波頻運動反應(yīng)的幅值，這表明風(fēng)力機(jī)氣動荷載產(chǎn)生的氣動阻尼能夠削弱平臺波頻運動反應(yīng)。

圖5 響應(yīng)頻譜

4.3 風(fēng)力機(jī)錨鏈張力反應(yīng)分析

在工況M3下，由圖6錨鏈張力反應(yīng)頻譜圖可知：浮式風(fēng)力機(jī)系泊系統(tǒng)在海洋環(huán)境中主要受到低頻、波頻和高頻成分的激勵作用。低頻成分與半潛浮式風(fēng)力機(jī)支撐平臺的低頻運動相關(guān)，主要包括平臺的縱蕩和縱搖運動，同時這也說明錨鏈張力反應(yīng)與支撐平臺的低頻運動之間存在耦合效應(yīng)；波頻成分主要由于波浪荷載的波頻激勵作用；而高頻成分主要是由塔架彈性反應(yīng)、轉(zhuǎn)子旋轉(zhuǎn)(3P)引起。

圖6 錨鏈張力反應(yīng)頻譜

圖7 5 MW與10 MW風(fēng)力機(jī)縱蕩運動響應(yīng)頻譜對比

5 10 MW與5 MW浮式風(fēng)力機(jī)動力特性對比

在工況M3下，使用FAST對NREL 5 MW無撐桿的半潛海上浮式風(fēng)力機(jī)[10]進(jìn)行全耦合分析，并與同工況下DTU 10 MW風(fēng)力機(jī)支撐平臺運動反應(yīng)和結(jié)構(gòu)受力對比如表3和表4所示。通過對比可知：10 MW浮式風(fēng)力機(jī)系統(tǒng)的受力與5 MW浮式風(fēng)力機(jī)相比均增加2倍以上，但10 MW浮式風(fēng)力機(jī)縱蕩運動響應(yīng)的平均值比5 MW浮式風(fēng)力機(jī)的僅增大了2%，這表明10 MW浮式風(fēng)力機(jī)需要錨鏈提供較大的縱蕩回復(fù)力剛度來降低縱蕩運動響應(yīng)的均值。此外，盡管10 MW浮式風(fēng)力機(jī)因其大的結(jié)構(gòu)自重和浮力可為縱搖運動和垂蕩運動提供較大的回復(fù)力和力矩，但同時也遭受更大的風(fēng)和波浪荷載，導(dǎo)致10 MW浮式風(fēng)力機(jī)縱搖和縱蕩運動響應(yīng)的平均值有所增大，但10 MW浮式風(fēng)力機(jī)整體運動性能良好，適合未來超大型浮式風(fēng)力機(jī)的應(yīng)用。

比較圖7中5 MW和10 MW浮式風(fēng)力機(jī)縱蕩運動響應(yīng)頻譜可知：10 MW浮式風(fēng)力機(jī)由風(fēng)荷載所激勵的低頻縱蕩運動的幅值明顯高于5 MW浮式風(fēng)力機(jī)低頻縱蕩運動幅值；而5 MW浮式風(fēng)力機(jī)由波浪荷載所激勵的波頻運動的幅值高于10 MW浮式風(fēng)力機(jī)波頻運動的幅值。故10 MW浮式風(fēng)力機(jī)的縱蕩運動以風(fēng)荷載激勵所產(chǎn)生的低頻運動為主，而5 MW浮式風(fēng)力機(jī)的縱蕩運動以波浪荷載激勵所產(chǎn)生的波頻運動為主。因此可預(yù)測：浮式海上風(fēng)力機(jī)大型化后，風(fēng)荷載對風(fēng)力機(jī)系統(tǒng)運動的低頻激勵作用更突出。

表4 5 MW與10 MW風(fēng)力機(jī)結(jié)構(gòu)受力對比

表3 5 MW與10 MW風(fēng)力機(jī)運動響應(yīng)對比

6 結(jié) 語

通過對DTU 10 MW浮式風(fēng)力機(jī)動力特性分析，并與NREL 5 MW浮式風(fēng)力機(jī)動力特性對比可知：

(1)在風(fēng)浪聯(lián)合作用下，半潛浮式風(fēng)力機(jī)平臺運動主要受到風(fēng)荷載的低頻激勵作用以及波浪荷載的波頻激勵作用；而錨鏈在海洋環(huán)境中主要受到波浪荷載的波頻激勵作用，以及由支撐平臺低頻運動引起的低頻激勵作用和塔架彈性反應(yīng)、轉(zhuǎn)子旋轉(zhuǎn)引起的高頻激勵作用。

(2)風(fēng)荷載對轉(zhuǎn)子作用所產(chǎn)生的氣動阻尼效應(yīng)能夠減小風(fēng)力機(jī)平臺波頻運動反應(yīng)的幅值；風(fēng)力機(jī)支撐平臺的縱蕩運動與縱搖運動，錨鏈張力與浮式風(fēng)力機(jī)支撐平臺的運動之間均存在耦合效應(yīng)。

(3)浮式海上風(fēng)力機(jī)大型化后，風(fēng)力機(jī)平臺運動反應(yīng)和結(jié)構(gòu)受力大幅增加，且風(fēng)力機(jī)氣動荷載激勵效應(yīng)更為突出。但風(fēng)力機(jī)系統(tǒng)整體運動性能良好，這為我國未來超大型浮式海上風(fēng)力機(jī)的設(shè)計與應(yīng)用提供參考。