楊 奇， 劉紅軍，2， 王建收

(1.中國海洋大學(xué)， 山東 青島 266100； 2.山東省海洋環(huán)境地質(zhì)工程重點實驗室，山東 青島 266100; 3.青島地質(zhì)工程勘察院，山東 青島 266100)

70%以上的地球表面被海洋覆蓋。對于人類來講，海洋是一個巨大的財富，特別是隨著科技的不斷發(fā)展，人們對能源的需求量日益增長，陸地上的不可再生能源已經(jīng)不能滿足人們的需求，因此對海洋能源的開發(fā)成為了當(dāng)務(wù)之急[1]。海上風(fēng)能是一種新型的再生能源，相比陸地上的風(fēng)力發(fā)電，有著一定優(yōu)勢。隨著國外海上風(fēng)機的快速發(fā)展，中國的海上風(fēng)力發(fā)電亟待進一步的研究[2]。

海上風(fēng)機基礎(chǔ)將上部荷載傳遞給地基，它的作用尤為重要。目前風(fēng)機基礎(chǔ)的形式有許多，如群樁、導(dǎo)管架、三腳架、重力式、承臺式、錨吸式基礎(chǔ)以及單樁基礎(chǔ)。其中單樁基礎(chǔ)，因其結(jié)構(gòu)較為簡單，占用面積小，承載力高，沉降量小且均勻，適用于灘涂和淺海0～30 m較大水深范圍和工程地質(zhì)條件復(fù)雜的場所，被廣泛應(yīng)用于海上風(fēng)機基礎(chǔ)，使用比例高達60%以上[3]。

海上風(fēng)機所處的海洋環(huán)境導(dǎo)致了樁基礎(chǔ)長期承受波浪荷載、風(fēng)荷載、潮流荷載等，增加了風(fēng)機結(jié)構(gòu)分析的難度[4]。郭健等[5]選用了不同的單元類型開展了海上風(fēng)機支撐結(jié)構(gòu)模態(tài)分析；陳前等[6-7]對海上風(fēng)機支撐結(jié)構(gòu)的動力特性進行探討。目前，海上風(fēng)機支撐結(jié)構(gòu)的設(shè)計主要采用集中質(zhì)量模型和完整模型，但對兩者的差異研究很少。

本文建立了兩種風(fēng)機支撐系統(tǒng)計算模型，通過輸入P-y曲線，并模擬施加波浪荷載、潮流荷載、風(fēng)荷載等，分別對兩種模型進行靜力分析和動力分析，探討在風(fēng)浪流荷載作用下兩種支撐系統(tǒng)模型的響應(yīng)。

1 數(shù)學(xué)模型

1.1 樁-土相互作用

地基反力法應(yīng)用Winkler地基模型，把樁周土離散為一個個單獨作用的彈簧。而在地基反力法中應(yīng)用最多的是荷載傳遞曲線法，即P-y曲線法，基本方程如式(1)。P-y曲線法就是在水平力H的作用下，泥面以下深度X處的土反力P與該點樁的撓度y之間的關(guān)系曲線。它綜合反映了樁周土的非線性、樁的剛度和外載荷作用性質(zhì)等特點[8-9]。王騰等[10]通過水平荷載下的單樁模型試驗，以API砂土的P-y曲線表達式為基礎(chǔ)，得到了黃河粉土的P-y曲線表達式(2)，并加以驗證，該表達式適用于黃河粉土。本文采用地基反力法應(yīng)用Winkler地基模型，把樁周土離散為若干單獨作用的彈簧，P-y曲線計算基于式(2)，樁基及土體參數(shù)見下文。

(1)

P=1.26Putanh(kzy/Pu)。

(2)

式中：P(x,y)為地基反力；y值為撓度；x值為沿樁身的長度；EI為樁的抗彎剛度；b為樁截面的計算寬度；Pu為承載力極限；k為初始地基反力模量；z為深度。

1.2 模態(tài)分析基本方程

風(fēng)機支撐系統(tǒng)典型無阻尼自由振動的基本方程如下式：

(3)

令：

{U}={φ}sin(ωt+φ)，

(4)

則有：

代入(3)式可得：

([K]-ω2[M]){φ}={0}。

(5)

由(4)式可得出方程的特征值ωi和特征向量{φi}，從而求得風(fēng)機支撐系統(tǒng)的振動頻率[11]。

1.3 瞬態(tài)分析動力方程

瞬態(tài)分析就是時間歷程分析，用于計算風(fēng)機支撐結(jié)構(gòu)在循環(huán)荷載作用下的動力響應(yīng)[12]。其基本運動方程如下：

(6)

在進行循環(huán)荷載作用下的瞬態(tài)分析時，必須確定阻尼系數(shù)。本文采用比例阻尼法[13]，即：

C=αM+βK。

(7)

根據(jù)振型分解法，可分別求出阻尼系數(shù)α、β：

(8)

式中：ξ為振型阻尼比，對于海上風(fēng)機，建議取值0.02～0.05，本文取0.03[13]；ωi為風(fēng)機系統(tǒng)的i階固有圓頻率，根據(jù)ω=2πf得到。

2 風(fēng)機系統(tǒng)及其環(huán)境參數(shù)

黃河三角洲擬建海上風(fēng)機發(fā)電系統(tǒng)，風(fēng)機機型為W3600，風(fēng)輪和機艙總重260 t。海面以上為風(fēng)機塔架和三葉式風(fēng)輪，葉片長度65 m，葉輪直徑122 m。塔架與過渡段視為整體，長64 m，直徑從4 m過渡到3 m[14]?；A(chǔ)形式采用單樁基礎(chǔ)，設(shè)計樁徑4 m，壁厚55 mm，樁長30 m，插入海床面以下20 m，海床土體主要為粉土，粘聚力c=12 kPa，內(nèi)摩擦角φ=25°。50年一遇的最大風(fēng)速為31 m/s[15]，海水深度為10 m，流速在2.5～40 cm/s(本文取40 cm/s)，采用50年一遇的波浪荷載，有效波高6.7 m，有效周期8.6 s[16]。風(fēng)機葉片材料采用環(huán)氧玻璃鋼，其余結(jié)構(gòu)材料均采用同一種鋼材。

3 兩種模型的建立

3.1 集中質(zhì)量模型

通過ANSYS有限元軟件建模并劃分網(wǎng)格計算，模型采用PIPE16、PIPE59、COMBIN39和MASS21四種單元類型。由于海床面以下樁體結(jié)構(gòu)受土體約束可不考慮塑性，采用PIPE16單元即可。而海床面以上的自由樁體結(jié)構(gòu)在風(fēng)浪流作用下可能產(chǎn)生部分塑性變形，故采用PIPE59單元(考慮塑性)。COMBIN39用于模擬樁土相互作用的關(guān)系，輸入王騰等人經(jīng)驗證的P-y曲線，模型合理可靠。MASS21用于動力分析時的集中質(zhì)量結(jié)構(gòu)。邊界條件采用樁端固定約束以及土體水平方向的位移約束，塔頂風(fēng)荷載等效為水平集中力。波浪荷載以及潮流作用通過Water Table輸入，自動施加在結(jié)構(gòu)上，可以得出樁底最大反力，與其對應(yīng)的波浪潮流耦合作用最大，即波浪相位角為Ф=61°，樁底反力與相位角的關(guān)系如圖1，集中質(zhì)量模型如圖2(a)。

圖1 樁底反力與相位角的關(guān)系Fig.1 The relationships between phase angle and reaction force in bottom of pile

3.2 風(fēng)機完整模型

如圖2(b)所示，模型采用PIPE16、PIPE59、COMBIN39和BEAM188四種單元類型。PIPE16、PIPE59、COMBIN39單元除上述功能外，PIPE16還用于機艙和輪轂的建立，BEAM188用于模擬風(fēng)機葉片。波浪荷載以及潮流作用同上。

4 結(jié)果分析

4.1 靜力分析

數(shù)值計算輸入的波浪潮流參數(shù)詳見第二節(jié)，將各參數(shù)輸入Water table自動換算成波浪力，風(fēng)速通過風(fēng)壓理論轉(zhuǎn)化為塔頂水平集中荷載。靜力求解得到了節(jié)點水平方向的位移分布如圖3所示，集中質(zhì)量模型中支撐結(jié)構(gòu)的最大水平位移(0.378 2 m)出現(xiàn)在塔頂，完整風(fēng)機模型的最大水平位移出現(xiàn)在葉片上，支撐結(jié)構(gòu)塔頂?shù)乃轿灰茷?.364 6 m；節(jié)點的水平方向應(yīng)力分布如圖4所示，集中質(zhì)量模型在海床面附近的水平應(yīng)力值最大為0.661×108Pa，而完整模型在機艙附近的水平應(yīng)力值最大為0.443×109Pa；圖5為結(jié)構(gòu)彎矩圖，兩種模型支撐結(jié)構(gòu)在土層中部彎矩最大，分別為0.245×108和0.263×108N·m。兩種模型支撐結(jié)構(gòu)的水平位移、水平應(yīng)力和彎矩均存在一定差異，且最大水平應(yīng)力差異很大。從樁底至塔頂?shù)乃轿灰品植既鐖D6所示，兩種模型支撐結(jié)構(gòu)的水平位移均是從樁底至塔頂非線性增加，但集中質(zhì)量模型塔架上部的水平位移低于完整模型，下部及樁基水平位移高于完整模型。為了保證風(fēng)機支撐系統(tǒng)的承載能力滿足要求，選用完整風(fēng)機模型更為適合。

圖2 兩種有限元模型Fig.2 Finite elements of two models

圖3 水平位移分布Fig.3 Horizontal displacement distribution

圖4 水平應(yīng)力分布Fig.4 Horizontal stress distribution

圖5 彎矩分布Fig.5 Moment distribution

圖6 樁底至塔頂水平位移變化曲線Fig.6 Horizontal displacement curve from bottom in pile to top in tower

4.2 模態(tài)分析

基于式(3)～(5)，模態(tài)分析時采用Block Lanczos法(可以提取中型到大型，50 000～100 000個自由度的大量振型，常用在實體單元和殼單元模型中)，得到振動頻率隨階數(shù)的變化如圖7，兩種模型的振動頻率在前兩階相差較小，高階振動頻率相差很大，原因可能是兩種模型采用了不同類型的單元。完整風(fēng)機模型的振動頻率隨階數(shù)的變化相對平緩，而集中質(zhì)量模型的振動頻率隨著階數(shù)的變化急劇增大。兩種模型求得的振動頻率差異顯著。

提取前10階模態(tài)，前6階振型如圖8。完整風(fēng)機模型在低階振動時，主要為葉片的轉(zhuǎn)動和擺動以及塔架的擺動，在高階自振時，才會出現(xiàn)單樁振動；集中質(zhì)量模型在前兩階主要為塔架的擺動，三階和四階出現(xiàn)樁基的擺動以及塔架的轉(zhuǎn)動。說明了風(fēng)機在低階振動時對支撐結(jié)構(gòu)的影響不大，主要是葉片、機艙、輪轂和塔架的振動。

圖7 振動頻率隨階數(shù)的變化曲線Fig.7 The curve of natural frequencies change with order

4.3 瞬態(tài)分析

基于式(6)～(8)瞬態(tài)分析采用完全法，阻尼系數(shù)α、β由模態(tài)分析得出的前兩階振動頻率代入式(8)求得。前100 s響應(yīng)內(nèi)兩種模型計算結(jié)果見表1。兩種模型計算結(jié)果均表明支撐結(jié)構(gòu)塔頂處將產(chǎn)生很大的水平位移，在設(shè)計時應(yīng)重點考慮。完整風(fēng)機模型的水平位移較集中質(zhì)量模型的水平位移相差0.016 m，在海面處的水平位移相差最大為0.020 m。完整風(fēng)機模型相對于集中質(zhì)量模型的最大水平位移出現(xiàn)時間有一定的滯后性。各關(guān)鍵位置(見圖9中從上至下曲線對應(yīng)位置依次為塔頂、海面、海底和樁底)的水平位移時程曲線如圖9。支撐結(jié)構(gòu)塔頂?shù)乃轿灰葡鄬τ谄渌恢玫乃轿灰聘资苡绊憽?/p>

各荷載作用前期(前10 s)，位移隨時間無規(guī)律變化，而后期呈固定波形變化，但集中質(zhì)量模型相對完整風(fēng)機模型呈現(xiàn)更明顯的規(guī)律性，本文將風(fēng)荷載等效為靜力荷載，前期(前10 s)為風(fēng)浪流共同作用，后期是波浪潮流耦合作用。因此，進行設(shè)計時必須考慮前期風(fēng)浪流的耦合作用。

圖8 兩種模型的前六階振型Fig.8 The first six order modes of two models

圖9 水平位移時程曲線

/m

5 結(jié)論

本文通過ANSYS有限元分析軟件，建立海上風(fēng)機支撐系統(tǒng)的兩種有限元計算模型并進行靜力和動力響應(yīng)對比分析，得出以下結(jié)論：

(1)風(fēng)機受風(fēng)浪流共同作用，而波浪相位角是波流耦合作用的影響因素。進行風(fēng)機動力分析時，必須考慮波浪的相位角。文中波流耦合作用最大時的相位角為Фmax=61°，最小為Фmin=14°。

(2)兩種模型的最大水平位移、最大水平應(yīng)力和最大彎矩分別相差3.73%、85.1%、7.3%。為了保證風(fēng)機支撐系統(tǒng)的安全，選取完整風(fēng)機模型進行設(shè)計分析更能降低風(fēng)機支撐系統(tǒng)承載破壞的可能。

(3)兩種模型各階振動形式不同，振動頻率相差顯著，隨著階數(shù)的增加，頻率差異越大。

(4)兩種模型對風(fēng)機支撐系統(tǒng)的動力響應(yīng)有顯著影響，且塔頂水平位移相對于其他位置的水平位移更易受影響。

[1] 李俊峰. 2012中國風(fēng)電發(fā)展報告[M]. 北京： 中國環(huán)境科學(xué)出版社， 2012.

Li Jun-Feng. China Wind Power Development Report in 2012 [M]. BeiJing: China Environmental Science Press, 2012.

[2] 辛華龍. 中國海上風(fēng)能開發(fā)研究展望[J]. 中國海洋大學(xué)學(xué)報(自然科學(xué)版)， 2010(6): 147-152.

XIN Hua-Long. Aspect on the development of offshore wind energy in China [J]. Periodical of Ocean University of China, 2010, 40(6): 146-152.

[3] 王淼， 吳云青， 蘇萌， 等. 單樁式基礎(chǔ)應(yīng)用于我國海上風(fēng)電的可行性探討[J]. 電力建設(shè), 2013(4): 63-66.

WANG Miao, WU Yunqing, SU Meng, et al. Application of single pile foundation in offshore wind power in China [J]. Electric Power Construction, 2013, 34(4): 63-66.

[4] 王國粹， 王偉， 楊敏. 3. 6MW海上風(fēng)機單樁基礎(chǔ)設(shè)計與分析[J]. 巖土工程學(xué)報， 2011, 33(S2)： 95-100.

WANG Guocui,WANG Wei,YANG Min. Design and analysis of monopile foundation for 3. 6 MW offshore wind turbine [J]. Chinese Journal of Geotechnical Engineering, 2011, 33(S2): 95-100.

[5] 郭健， 楊敏， 王偉. 海上風(fēng)機高樁承臺基礎(chǔ)模態(tài)建模分析研究[J]. 巖土工程學(xué)報， 2013， 35(S2): 1172-1175.

GUO Jian, YANG Min, WANG Wei. Modeling analysis of modals of high-pile cap foundation of offshore wind turbines [J]. Chinese Journal of Geotechnical Engineering, 2013, 35(92): 1172-1175.

[6] 陳前， 付世曉， 鄒早建. 海上風(fēng)力發(fā)電機組支撐結(jié)構(gòu)動力特性分析[J]. 振動與沖擊， 2012， 31(2): 86-90.

CHEN Qian, FU Shi-Xiao, ZOU Zao-Jian. Dynamic characteristics analysis for support structure of an offshore wind turbine[J]. Journal of Vibration and Shock, 2012, 31(2): 86-90.

[7] 陳前， 付世曉， 鄒早建. 海上風(fēng)力發(fā)電機組支撐結(jié)構(gòu)動力響應(yīng)特性研究[J]. 船舶力學(xué)， 2012， 3(4): 408-415.

CHEN Qian, FU Shi-Xiao, ZOU Zao-Jian. Investigation on dynamic response characteristics of offshore wind turbine support structure [J]. Journal of Ship Mechanics, 2012, 31(4): 408-415.

[8] 蘇靜波， 邵國建， 劉寧. 基于P-Y曲線法的水平受荷樁非線性有限元分析[J]. 巖土力學(xué)， 2006， 27(10): 1781-1785.

SU Jingbo, SHAO Guojian, LIU Ning. Nonlinear finite element analysis of piles under lateral load based on P-Y curves [J]. Rock and Soil Mechanics, 2006, 27(10): 1781-1785.

[9] Jin W. Evaluation of Application of p-y curve method in offshore engineering foundation [J]. China Ocean Engineering, 1993(4): 451-466.

[10] 王騰， 王天霖. 粉土 p-y 曲線的試驗研究[J]. 巖土力學(xué)， 2009， 30(5): 1343-1346.

WANG Teng, WANG Tianlin. Experimental research on silt p-y curves[J]. Rock and Soil Mechanics, 2009, 30(5): 1343-1346.

[11] 靳軍偉， 楊敏， 王偉， 等. 海上風(fēng)電機組單樁基礎(chǔ)模態(tài)及參數(shù)敏感性分析[J]. 同濟大學(xué)學(xué)報(自然科學(xué)版)， 2014， 43(3): 386-392.

JIN Junwei, YANG Min, WANG Wei, et al. Offshore wind turbine monopile foundation modal and parameter sensitivity analysis[J]. Journal of Tongji University: Natural Science, 2014, 42(3): 386-392.

[12] Gao Ming, Chen Jin-Zhen, Zhen Guo-Fang. Soil pile interaction under static, dynamic and cycle lateral loads and a proposal of p-y Curve formula[J]. China Ocean Engineering, 1989(3): 259-270.

[13] 楊進， 劉書杰， 謝仁軍. ANSYS在海洋石油工程中的應(yīng)用[M]. 北京： 石油工業(yè)出版社， 2010.

YANG Jin, LIU Shu-Jie, XIE Ren-Jun. Application of ANSYS in Offshore oil Engineering[M]. BeiJing: Petroleum Industry Press, 2010.

[14] Ganeriwala N S, Yang J, Richardson M. Using modal analysis for detecting cracks in wind turbine blades[J]. Sound & Vibration, 2011(5): 10.

[15] 劉晉秀， 江崇波， 范學(xué)煒. 黃河三角洲近40年來氣候變化趨勢及異常特征[J]. 海洋預(yù)報， 2002, 19(2): 31-35.

Liu Jinxiu, Jiang Chongbo, Fan Xuewei. The climate change tendency and the abnormal characters of the Yellow River Delta in recent 40 years [J]. Marine Forecasts, 2002, 19(2): 31-35.

[16] 劉紅軍， 王虎， 張民生， 等. 波浪作用下黃河三角洲粉質(zhì)土海床動力響應(yīng)分析[J]. 巖土力學(xué)， 2013， 34(7): 2065-2071.

LIU Hong-Jun, WANG Hu, ZHANG Min-Sheng, et al. Analysis of wave-induced dynamic response of silty seabed in Yellow River Delta [J]. Rock and Soil Mechanics, 2013, 34(7): 2065-2071.