程友良，薛占璞，李占嶺，趙懷宇，趙洪嵩，賀 忠

(1.華北電力大學(xué) 能源動力與機械工程學(xué)院，河北 保定 071003;2.河北電力勘測設(shè)計研究院，石家莊 050000；3.石家莊鐵道大學(xué)，石家莊 050000)

?

海上半潛式風(fēng)力發(fā)電塔基礎(chǔ)的水動力學(xué)特性分析

程友良1，薛占璞1，李占嶺2，趙懷宇2，趙洪嵩1，賀 忠3

(1.華北電力大學(xué) 能源動力與機械工程學(xué)院，河北 保定 071003;2.河北電力勘測設(shè)計研究院，石家莊 050000；3.石家莊鐵道大學(xué)，石家莊 050000)

提出一種強度與流場分析相結(jié)合的方法,利用Solidworks軟件對半潛式基礎(chǔ)進行建模及應(yīng)力分析，并進行模態(tài)分析，比較不同頻率之間振型，得出振動特性，利用FLUENT軟件單向耦合分析法進行速度場和壓力場數(shù)值模擬及水動力學(xué)特性分析。結(jié)果表明，橫梁三角形各個角度的流場和阻力特性變化呈非線性關(guān)系，影響流速大小，進而對于側(cè)向流及漩渦的形成產(chǎn)生影響，并且側(cè)向流越大，越容易形成漩渦。

海上半潛式風(fēng)力發(fā)電塔；應(yīng)力分析；振動特性；水動力學(xué)特性；數(shù)值模擬

近年來海上風(fēng)電發(fā)展迅速，海上風(fēng)力發(fā)電塔架高度逐漸增加，高度超過50 m，甚至達到120 m或更高[1]。國外對于浮式基礎(chǔ)的研究較早，半潛式基礎(chǔ)概念和設(shè)計許多由荷蘭多家公司及倫敦大學(xué)等機構(gòu)提出并且實施。國內(nèi)相關(guān)的研究工作[2-12]僅從結(jié)構(gòu)動力學(xué)或者僅從流體動力學(xué)方面分析海上風(fēng)機基礎(chǔ)，而將兩者結(jié)合分析整體基礎(chǔ)強度較少。

為此，考慮從結(jié)構(gòu)動力學(xué)和流體動力學(xué)兩個方面分析半潛式浮式基礎(chǔ)，采用數(shù)值模擬和湍流理論結(jié)合方法，創(chuàng)新地將機械強度與流體分析結(jié)合應(yīng)用到模型中，分析其在不同載荷下應(yīng)力變化幅度情況，不同頻率下模態(tài)振型，定量分析半潛式浮式基礎(chǔ)整體周圍流場分布及相關(guān)參數(shù)。

1 數(shù)學(xué)模型及有限元分析

1.1 模型建立及網(wǎng)格劃分

半潛式基礎(chǔ)在運輸、控制方面有良好的優(yōu)越性，近年來在近海區(qū)域利用率逐漸提高[13-17]。以海上風(fēng)力發(fā)電塔半潛式Windfloat基礎(chǔ)發(fā)電系統(tǒng)實例，見圖1。浮式基礎(chǔ)材料為高強度鋼，屈服極限為355 MPa。為了較好地分析Windfloat基礎(chǔ)，對其進行簡化，見圖2。

圖1 半潛式Windfloat基礎(chǔ)風(fēng)力發(fā)電系統(tǒng)

圖2 Windfloat基礎(chǔ)模型

利用Solidworks軟件模塊simulation對模型進行網(wǎng)格劃分，設(shè)置網(wǎng)格的質(zhì)量為良好，網(wǎng)格類型為三角形網(wǎng)格，采用自動劃分的方式，網(wǎng)格化模型見圖3，節(jié)點數(shù)為117 615，單元數(shù)為99 485。

圖3 網(wǎng)格化模型

1.2 仿真分析

1.2.1 數(shù)學(xué)模型驗證

按照文獻[8]極限海況下的載荷，許用應(yīng)力284 MPa,最大等效應(yīng)力202.364 MPa，最大位移39.538 mm進行設(shè)置。對比文獻[8]中LC6浮式基礎(chǔ)整體應(yīng)力云圖4b)，可知圖4a)、b)應(yīng)力變化幅度相同，模擬是可行的。

在極限工況下，施加載荷，目的在于分析其承載的極限工況[11-13]。

圖4 Windfloat基礎(chǔ)應(yīng)力云圖

從應(yīng)力云圖可知，承受風(fēng)塔的樁受力在應(yīng)力范圍允許之內(nèi)，整個樁應(yīng)力變化不大，在安全范圍內(nèi)，符合強度要求。

從位移云圖可知(圖5)，承受風(fēng)塔的樁穩(wěn)定性尤為重要，直接關(guān)系風(fēng)塔與樁的相對位置。進而影響整個系統(tǒng)的穩(wěn)定。圖中位移在0.15 mm與0.02 mm之間浮動，符合穩(wěn)定性要求。

圖5 Windfloat基礎(chǔ)位移云圖

1.2.1 模態(tài)分析

在Windfloat基礎(chǔ)固有頻率情況下，提取前5階非零模態(tài)振型。與前5階非零模態(tài)振型對應(yīng)的頻率分別為58、102、238、307及348 Hz。

Windfloat基礎(chǔ)自由振動方程為

(1)

式中：ω——自由振動固有頻率，特征方程為

(2)

即可求出多項式的根(特征值)即為離散模型的固有頻率，求出特征向量，獲取相應(yīng)的振型。

由振型圖分析可知，Windfloat基礎(chǔ)第3階模態(tài)(f=238 Hz)變形量大，最大位移在22 mm(圖6)，應(yīng)在實際運行中監(jiān)測頻率，避免在此頻率下運行，其余振型位移量變化幅度不大，屬于正常運行。

圖6 Windfloat基礎(chǔ)3階振型

經(jīng)過應(yīng)力及模態(tài)分析可知，此模型能夠滿足強度及穩(wěn)定性要求。

2 水動力學(xué)特性分析

2.1 Windfloat基礎(chǔ)模型

將Solidworks建立的基礎(chǔ)模型導(dǎo)入ANSYS ICEM CFD中，進行模型設(shè)置，幾何檢查及修復(fù)，如圖7為ICEM CFD環(huán)境中的幾何模型。

圖7 ICEM CFD幾何模型

2.2 控制方程及湍流模型

Windfloat基礎(chǔ)附近海域為粘性不可壓縮流體，運動類型為湍流運動?？刂品匠虨?/p>

連續(xù)性方程：

(3)

動量方程：

(4)

湍動能k的方程：

(5)

湍動能耗散率ε的方程：

(6)

為了使湍流運動分布均勻，連續(xù)性強，這里采用k-ε湍流模型，魯棒性好，此類模型可以有效模擬海上水域流動特性[18-22]。

2.3 計算域及網(wǎng)格劃分

Windfloat基礎(chǔ)計算域如圖8所示，計算域大小為長40 m、寬30 m、高40 m。Windfloat基礎(chǔ)位置分布在長方體區(qū)域長寬高1/4處，在ICEM CFD環(huán)境中劃分網(wǎng)格，網(wǎng)格類型為Tetra/Mixed，節(jié)點數(shù)為101 699，單元數(shù)為98 488，結(jié)構(gòu)中立柱網(wǎng)格邊長及密度較大，邊長為0.09，斜撐和橫梁密度較小。

圖8 計算域位置圖

網(wǎng)格質(zhì)量檢查良好，能夠保證在有限元計算中準(zhǔn)確性。圖9為ICEM CFD環(huán)境中的網(wǎng)格模型。增加網(wǎng)格的密度，其模擬結(jié)果誤差在0.4%之內(nèi)，進行網(wǎng)格無關(guān)性驗證，證明此網(wǎng)格模擬可行[23-27]。

圖9為網(wǎng)格化的Windfloat基礎(chǔ)，通過網(wǎng)格收斂性驗證，有限元網(wǎng)格單元尺寸大小對計算結(jié)果的影響，單元尺寸減小而應(yīng)力變化不明顯，控制在1%之內(nèi)，因此計算是可行的。

圖9 網(wǎng)格化模型

2.4 邊界條件設(shè)置

計算域中左端面為速度進口，數(shù)值為1.8 m/s，這符合海上實際流速情況；右端面為壓力出口；2個側(cè)面、頂面為對稱邊界條件；底面以及Windfloat基礎(chǔ)各面為無滑移壁面。

2.5 數(shù)學(xué)模型的驗證

邵萬駿等[3]在人工魚礁水動力學(xué)特性實驗作出相關(guān)曲線，唐友剛[5]在海上風(fēng)機半潛式水動力性能得出相關(guān)實驗數(shù)據(jù)，這里基于FLUENT14.0單向耦合分析,數(shù)值模擬測得流速數(shù)據(jù)曲線，實驗測量點位置分布見圖10。

圖10 測量點位置

通過文獻[3]、文獻[5]流速的實驗值與本人模擬值的比較可知，圖11各個點實驗值與模擬值的誤差在允許范圍之內(nèi)，此數(shù)值模擬是可行的。

圖11 測量點實驗值與模擬值比較

2.6 數(shù)值模擬結(jié)果及分析

當(dāng)計算域中Windfloat基礎(chǔ)三角橫梁角度為60°,60°,60°，來流速度為1.2 m/s時，出口垂直剖面的壓力點云分布如圖12。

圖12 垂直剖面壓力點云分布圖(x=0)

壓力分布是衡量海水流動方向及速度的重要參數(shù)，從圖12可知，垂直剖面壓力分布，接近對稱分布，連續(xù)性頻率的波浪力作用，導(dǎo)致壓力點云分布對稱，這說明在Windfloat基礎(chǔ)周圍水流動呈連續(xù)性，在實際的水運動中有一定的能量循環(huán)，符合水流動特性。當(dāng)位置從0～250 m及-250～0 m進行對比，可知壓力呈現(xiàn)類似正態(tài)分布點云圖，最大變化量為1 200 Pa，點云對稱分布，這為測量數(shù)據(jù)提供參考。

為了準(zhǔn)確描述Windfloat基礎(chǔ)周圍流場分布，設(shè)置來流速度仍為1.2 m/s，如圖13為Windfloat基礎(chǔ)速度矢量圖。與此對應(yīng)的整體速度矢量見圖14。從圖13中可知，當(dāng)水流經(jīng)過Windfloat基礎(chǔ)樁時，流向改變，向兩邊排開，在樁的上部和下部形成小的旋渦，流速變小，沒有危險區(qū)域。在樁正前部，會產(chǎn)生側(cè)向流。從圖14可知，樁在經(jīng)過流水時，速度變化幅度不大，側(cè)向流在樁附近速度疊加，導(dǎo)致速度變化幅度加大。

圖13 垂直剖面速度矢量圖(x=0)

圖14 Windfloat基礎(chǔ)速度矢量圖

水經(jīng)過Windfloat基礎(chǔ)的阻力曲線是監(jiān)測流速平緩與否的重要參數(shù)，如圖15為阻力曲線，橫坐標(biāo)為迭代步數(shù)，從中可知，阻力從零開始有一段跳躍性變化，這是在剛開始監(jiān)測初始化。在200～500步時，阻力曲線趨于平滑，沒有大幅度增減。經(jīng)過三角橫梁角度改變，多次數(shù)值模擬阻力曲線，三角形橫梁的角度與阻力曲線為非線性關(guān)系。4條曲線由上到下對應(yīng)三角形橫梁的角度依次為：60°，60°，60°；30°,60°,90°；120°,30°,30°；45°,45°,90°；圖16為不同角度下的阻力，橫坐標(biāo)為迭代步數(shù)，從中可知，三角形橫梁角度應(yīng)根據(jù)海域流速大小及阻力變化幅度決定。

圖15 Windfloat基礎(chǔ)阻力

圖16 三角形橫梁不同角度阻力變化

2.7 機械結(jié)構(gòu)強度與流場結(jié)合分析

從以上結(jié)構(gòu)動力學(xué)與流體動力學(xué)分析可知，半潛式Windfloat基礎(chǔ)結(jié)構(gòu)強度應(yīng)力云圖與整體流場基礎(chǔ)矢量圖之間有一定聯(lián)系，只有在滿足強度要求下，周圍流場才能夠符合流體運動的一般規(guī)律。因此，在分析流場之前，應(yīng)進行結(jié)構(gòu)強度分析。此Windfloat基礎(chǔ)的應(yīng)力變化幅度，符合強度要求，剖面及整體速度矢量圖符合流體湍流運動的一般規(guī)律。

3 結(jié)論

1)通過對Windfloat基礎(chǔ)應(yīng)力及模態(tài)分析，得到不同應(yīng)力變化幅值及模態(tài)振型圖對結(jié)構(gòu)破壞至關(guān)重要，在實際應(yīng)用中，應(yīng)該在一定時間內(nèi)，監(jiān)測應(yīng)力變化及模態(tài)情況，避免發(fā)生結(jié)構(gòu)破壞。

2)Windfloat基礎(chǔ)周圍流場分布、阻力變化情況與流速、橫梁角度有關(guān)，當(dāng)橫梁角度全為60°時，阻力變化平緩，適合海水流速變化平穩(wěn)的情況；其余橫梁角度適合于波浪與海流耦合的情況。這些橫梁角度結(jié)合剖面的速度矢量分布，給實際制造及監(jiān)測提供參考。

3)橫梁三角形角度的不同將影響側(cè)向流大小，進而形成旋渦，側(cè)向流越大，形成速度疊加，越容易形成旋渦。利用結(jié)構(gòu)動力學(xué)與流體動力學(xué)結(jié)合的方法，將進一步對半潛式基礎(chǔ)的結(jié)構(gòu)設(shè)計，提供技術(shù)支持。但在實際阻力監(jiān)測中，數(shù)值模擬與實際監(jiān)測有誤差是必然的，在一定允許誤差范圍內(nèi)數(shù)據(jù)是正確的。海上半潛式基礎(chǔ)不僅受到波浪、海流作用，還受到風(fēng)、地震等載荷作用，因此，對半潛式基礎(chǔ)進行多載荷響應(yīng)分析是未來的研究重點。

[1] 李健.海上風(fēng)力發(fā)電塔架風(fēng)渦激振動及疲勞分析[D].哈爾濱：哈爾濱工程大學(xué)，2011.

[2] 唐友剛,王涵,陶海成,等.海上風(fēng)機半潛型浮式基礎(chǔ)結(jié)構(gòu)設(shè)計及整體強度分析[J].中國造船，2013,54(3)：85-93.

[3] 邵萬駿,劉長根,聶紅濤,等.人工魚礁的水動力學(xué)特性及流場效應(yīng)分析[J].水動力學(xué)研究與進展,2014,29(5):580-585.

[4] 陳志強.海上風(fēng)力發(fā)電機樁基礎(chǔ)波流荷載研究[D].大連：大連理工大學(xué)，2013.

[5] 唐友剛,桂龍,曹菡,等.海上風(fēng)機半潛式基礎(chǔ)概念設(shè)計與水動力性能分析[J].哈爾濱工程大學(xué)報，2014,35(11):1314-1319.

[6] 夏宗朝.半潛漂浮式風(fēng)力發(fā)電機組仿真技術(shù)研究[D].重慶：重慶大學(xué)，2014.

[7] 夏露.海上風(fēng)機基礎(chǔ)在風(fēng)浪作用下的動力分析與疲勞分析[D].大連：大連理工大學(xué)，2012.

[8] 陶海成.海上風(fēng)電浮式基礎(chǔ)結(jié)構(gòu)設(shè)計及整體強度分析[D].天津：天津大學(xué)，2012.

[9] 劉坤寧.海上浮式風(fēng)電結(jié)構(gòu)風(fēng)機載荷模擬試驗技術(shù)研究[D].青島：中國海洋大學(xué)，2013.

[10] 汪洋.海上浮式風(fēng)電基礎(chǔ)設(shè)計方案及垂直軸風(fēng)機性能分析[D].鎮(zhèn)江：江蘇科技大學(xué)，2014.

[11] 葉小嶸.海上浮式風(fēng)力機系統(tǒng)環(huán)境載荷及耦合運動性能研究[D].哈爾濱：哈爾濱工程大學(xué)，2014.

[12] 王涵.小水深半潛型風(fēng)電浮式基礎(chǔ)的耦合動力分析與試驗研究[D].天津：天津大學(xué)，2013.

[13] 唐衍力.人工魚礁水動力的實驗研究與流場的數(shù)值模擬[D].青島：中國海洋大學(xué)，2013.

[14] 陳小波,陳健云,李靜.海上風(fēng)力發(fā)電塔脈動風(fēng)速時程數(shù)值模擬[J].中國電機工程學(xué)報,2008,28(32):111-116.

[15] 王人鳳,顧曉波,蔣志勇.三浮體式海上風(fēng)電基礎(chǔ)結(jié)構(gòu)響應(yīng)研究[J].船舶工程, 2012(3):93-96.

[16] 姜昭陽.人工魚礁水動力學(xué)與數(shù)值模擬研究[D].青島：中國海洋大學(xué)，2009.

[17] 黃遠東，龍催，鄧濟通.三棱柱型人工魚礁繞流流場的CFD分析[J].水資源與水工程學(xué)報,2013,24(1):1-4.

[18] 李毓湘，逄勇.珠江三角洲地區(qū)河網(wǎng)水動力學(xué)模型研究[J].水動力學(xué)研究與進展,2001,16(2):143-154.

[19] 李東東，陳陳.風(fēng)力發(fā)電系統(tǒng)動態(tài)仿真的風(fēng)速模型[J].中國電機工程學(xué)報,2005,25(3):115-119.

[20] 李珺，林軍，章守宇.方形人工魚礁通透性及其對礁體周圍流場影響的數(shù)值實驗[J].上海海洋大學(xué)學(xué)報,2010,19(6):836-840.

[21] GUO Peng.Influence analysis of wind shear and tower shadow on load and power based on blade element theory[C]. Proceedings of the 2011 Chinese Control and Decision Conference(CCDC),Sichuan Mianyang,China.,2011.

[22] MATSUSHITA Daisuke, MATSUMIYA Hikaru, HARA Yoshinori, et al. Studies on numerical site calibration over complex terrain for wind turbines[J]. Science China,2010(1):8-12.

[23] 曹菡.海上風(fēng)機半潛型浮式基礎(chǔ)水動力性能研究[D].天津：天津大學(xué)，2012.

[24] 劉中柏.海上風(fēng)電浮式基礎(chǔ)運動特性試驗研究[D].天津：天津大學(xué)，2014.

[25] FEI Chaoyang, WANG Nan, ZHOU Bo, et al.Dynamic Performance Investigation for Large-scale Wind Turbine tower[C]. Proceedings of the Eighth International Conference on Electrical Machines and Systems(Volume Ⅱ), Nanjing, china, 2005.

[26] FENG Gao, PENG Guo. Individual Pitch Control of Doubly-fed Wind Turbine Generator Based on Blade Element Theory[C]. Twenty-fourth Chinese Control and Decision Making,Shanxi,Taiyuan,China,2012.

[27] 鄧濟通,黃遠征,姜劍偉,等.布設(shè)間距對三棱柱形人工魚礁繞流影響的數(shù)值模擬[J].水資源與水工程學(xué)報,2013,24(2):98-103.

On Hydrodynamic Characteristics of the Foundation of Offshore Semi-submersible Type of Wind Power Tower

CHENG You-liang1, XUE Zhan-pu1, LI Zhan-ling2, ZHAO Huai-yu2, ZHAO Hong-song1, HE Zhong3

(1.School of Energy Power and Mechanical Engineering, North China Electric Power University, Baoding Hebei 071003, China;2.Hebei Electric Power Survey and Design Institute, Shijiazhuang 050000, China;3.Shijiazhuang Railway University, Shijiazhuang 050000, China)

For the semi-submersible foundation of offshore wind power tower, the stress and hydrodynamic characteristics were studied numerically. In Solidworks, the structural model was built up to carry out the stress and modal analysis. The vibration modes of different frequencies were compared to study its vibration characteristics. In FLUENT, the velocity field and pressure field were simulated numerically to study the hydrodynamic characteristics. The results showed that the angle of the cross beam triangle has a certain nonlinear influence on the flow velocity and the resistance, as well as the formation of lateral flow and the vortex; and the larger the lateral flow is, the more easily a whirlpool forms.

offshore semi-submersible type of wind power tower; stress analysis; vibration characteristics; hydrodynamic characteristics; numerical simulation

10.3963/j.issn.1671-7953.2016.05.038

2016-04-12

國家自然科學(xué)基金重點基金(11232012)

程友良(1963—)，男，博士，教授

U661.1;TK 83

A

1671-7953(2016)05-0151-06

修回日期：2016-04-29

研究方向:流體動力學(xué)理論及其應(yīng)用、流體設(shè)備與節(jié)能、可再生能源理論及其應(yīng)用、清潔能源利用技術(shù)與設(shè)備

E-mail:shenghuo66@163.com