考慮樁-彈性地基相互作用的單樁風機自振頻率

牛文杰

(遼寧工程技術(shù)大學力學與工程學院，遼寧 阜新 123000)

?

考慮樁-彈性地基相互作用的單樁風機自振頻率

牛文杰

(遼寧工程技術(shù)大學力學與工程學院，遼寧 阜新 123000)

考慮樁-彈性地基相互作用，采用集中質(zhì)量法和柔度法對泥面線以上的單樁風機結(jié)構(gòu)進行多自由度動力分析，確定單樁風機的自振頻率。通過算例給出單樁風機的自振頻率值，并與不考慮樁-彈性地基相互作用的單樁風機的自振頻率值進行比較。

單樁風機； 自振頻率； 樁基礎和彈性地基的相互作用； 柔度法； 集中質(zhì)量法

0 引言

對于水平軸風力發(fā)電機[1]，大型風力機塔筒通常都采用變截面的錐形筒體形式，可將其視為頂端受橫向力作用的變截面懸臂梁結(jié)構(gòu)[2]。橫向力由風力機葉片、輪轂和機艙的迎風阻力組成[3]。風力發(fā)電機組運行時，其葉輪上的風荷載和風機偏航引起的荷載通過結(jié)構(gòu)和傳動機構(gòu)作用在塔架頂端[4-5]。已建成的大部分海上風電場都采用單立柱基礎[4-6]。風機塔上作用有自重荷載、波浪荷載、流體流過立柱的流體壓力荷載[5]以及可能的近海地震作用與浮冰沖擊作用。此外，風機葉片的轉(zhuǎn)速一般是常數(shù)或在某一指定范圍，可通過閉環(huán)控制系統(tǒng)實現(xiàn)[7]。因此，對于葉片轉(zhuǎn)速一定的風機，當以塔筒為研究對象時，其上的動力荷載還主要有：輪轂等轉(zhuǎn)子系統(tǒng)施加的動荷載，其頻率為一定值，記為f；當風機葉片在轉(zhuǎn)動過程中遮擋住塔筒時，改變了環(huán)境流場，會施加給塔筒荷載，如果是三葉片的風機，則葉片轉(zhuǎn)動引起的作用在塔筒上的振動荷載頻率是3f[8]。

在上述動力荷載作用下，需要通過工程設計，使現(xiàn)場近海單樁風機的自振頻率避開上述動力荷載的激勵頻率，以避免近海單樁風機在動力荷載作用下發(fā)生共振。另外一個工程問題就是：對于近海風機這種高而柔的結(jié)構(gòu)物需要承受各種動力和循環(huán)荷載，其在長期的水平荷載作用下，由于樁周土在動應力作用下的材料弱化，風機的自振頻率會發(fā)生變化。因此需要研究確定近海單樁自振頻率的計算方法或?qū)嶒灧椒ā?/p>

確定近海單樁自振頻率的實驗方法通常采用Bhattacharya等(2011)建議的模型試驗，以保證模型與現(xiàn)場風機的相似準數(shù)一致[9]。

確定近海單樁自振頻率的理論解有以下兩方面的研究：

(1) 計算近海單樁風機的自振頻率需要考慮樁土相互作用。而關(guān)于樁土相互作用方面的研究有：Zaaijer[10]研究了單樁-土體的動力特性，土體被簡化為一系列彈簧；Padrn等[11]通過邊界元-有限元耦合模型對埋在半空間彈性體中的樁和群樁進行了動力分析，其中樁按照伯努利假設的梁單元進行模擬，而土按照連續(xù)的、半無限體、各向同性的線彈性介質(zhì)通過邊界元來進行模擬；Dash等[12]使用數(shù)值方法研究了土體中樁發(fā)生壓-彎耦合時的性狀。

(2) 牛文杰等[13]使用集中質(zhì)量法，采用剛度法對近海風力發(fā)電機的塔和泥面線上基礎按照2個自由度的動力系統(tǒng)進行了自由振動分析,并得到自振頻率。但考慮地基-風機相互作用時的計算模型更加符合實際。Bhattacharya[14]按照無限自由度體系，把地基對塔筒的作用簡化為轉(zhuǎn)動彈簧和平動彈簧，彈簧剛度可使用試驗測得，并使用解析解計算了風機的自振頻率，模型試驗和有限元解的結(jié)果表明計算結(jié)果是合理的。另外，可使用p-y曲線法對海洋平臺的樁土結(jié)構(gòu)系統(tǒng)進行模態(tài)分析[15]。

上述計算自振頻率的方法，或過于簡略，或過于復雜，或需要大量的現(xiàn)場測試數(shù)據(jù)，或需要借助于Ansys軟件進行分析計算，因此需要一種快速估算近海單樁風機自振頻率的計算方法。

風力發(fā)電機的上部塔筒和樁基礎都假定為彈性模量相同的鋼材制成的套筒。通過張氏法[16-17]得到半無限長樁基礎的控制方程；然后通過歐拉-伯努利經(jīng)典梁理論方程得到泥面線上以上風機上部塔筒的變形控制方程；再通過計算變形的疊加原理與套筒在泥面線處的變形連續(xù)性條件，得到泥面線上風機立柱(即套筒)的撓度函數(shù)；最后使用柔度法，給出按照多自由度分析的考慮樁-彈性地基相互作用時風機的自振頻率值。

1 考慮樁-彈性地基相互作用的自由振動理論

1.1 考慮樁土相互作用時單樁結(jié)構(gòu)物的撓度

圖1所示為單樁風力發(fā)電機。以下分析中，把風力機簡化為截面和彈性模量都相同的單樁結(jié)構(gòu)物(圖2)。風機高度為h。

當水平荷載作用于桿件頂端時，對于圖2中的單樁結(jié)構(gòu)物即為等截面桿件在泥面線以上的水平撓度y(x)F。這可按照兩部分來計算：一部分是樁底端固定時的位移曲線函數(shù)y(x)F0；另一部分是假定無外荷載作用，泥面線位置處單樁結(jié)構(gòu)物轉(zhuǎn)角θ0和位移y0引起的樁的水平撓度y(x)F1。最后通過疊加原理可確定y(x)F。

(1)

圖1 單樁式風力發(fā)電機Fig.1 A mono-pile wind turbine

圖2 單樁結(jié)構(gòu)物頂部(桿件)作用水平荷載FFig.2 Mono-pile structure top under horizontal load F

在圖2中，以泥面線與桿件的交點為坐標原點，記x是泥面線上桿件某點與原點的距離。根據(jù)材料力學，可確定y(x)F0：

(2)

在圖2中，泥面線位置處單樁結(jié)構(gòu)物的轉(zhuǎn)角和位移引起的單樁結(jié)構(gòu)物在泥面線以上的水平撓度

(3)

式中：θ0和y0分別是泥面線位置處單樁結(jié)構(gòu)物立柱的轉(zhuǎn)角和撓度(圖2)。

綜合式(1)、(2)和(3)，可得圖2中單樁結(jié)構(gòu)物在頂部作用有水平集中荷載時的水平撓度：

(4)

(1) 泥面線位置處單樁結(jié)構(gòu)物的轉(zhuǎn)角和位移

圖2中以泥面線與單樁結(jié)構(gòu)物的交點為坐標原點，記z是泥面線下樁基礎某點與原點的距離， 即樁基礎的水平撓度。當單樁結(jié)構(gòu)物頂部作用有水平集中力F時，泥面線以下半無限長樁在坐標z處的水平撓度y的控制方程為[16-17]：

(5)

式(5)中物理量的含義可查找文獻[16-17]。

泥面線位置處：

(6)

將式(6)代入式(5)，得出泥面線位置處單樁結(jié)構(gòu)物的撓度：

(7)

在單樁結(jié)構(gòu)物發(fā)生小撓度的假定下，單樁結(jié)構(gòu)物軸線的轉(zhuǎn)角θ在任一點滿足以下關(guān)系：

(8)

觀察圖2可發(fā)現(xiàn)：

(9)

由式(5)可得單樁結(jié)構(gòu)物軸線在泥面線以下任一點的轉(zhuǎn)角的正切值為：

(10)

泥面線位置處單樁結(jié)構(gòu)物軸線的縱坐標滿足式(6)，將式(6)代入式(10)，得出泥面線位置處單樁結(jié)構(gòu)物軸線轉(zhuǎn)角的正切值：

(11)

由式(9)和(11)得出在圖2所示的水平力F作用下，泥面線位置處單樁結(jié)構(gòu)物軸線轉(zhuǎn)角的正切值為：

(12)

根據(jù)式(8)和(12)得泥面線位置處單樁結(jié)構(gòu)物軸線轉(zhuǎn)角:

(13)

(2) 單樁結(jié)構(gòu)物的水平撓度y(x)F

結(jié)合式(4)、(7)和(13)，得出在圖2所示的單樁結(jié)構(gòu)物在頂部作用有水平集中荷載F時，泥面線以上的水平撓度：

(14)

單樁結(jié)構(gòu)物頂端的位移為：

(15)

因此，式(14)代入式(15),得出單樁結(jié)構(gòu)物頂端的水平撓度為：

(16)

1.2 多自由度分析

按照柔度法對圖1的結(jié)構(gòu)進行自振頻率計算。記機艙、輪轂和葉片的總重量為m0，泥面線以上套筒質(zhì)量為M0。按照集中質(zhì)量法[20]，在圖3中m1=0.25*M0+m0,m2=0.5*M0。

泥面線以上塔筒高度為h。根據(jù)式(16)和(14),圖3中在水平集中力F作用下，集中質(zhì)量位置處的水平撓度分別為：

圖3 計算m1位置處的柔度系數(shù)Fig.3 Computation of flexibility coefficient at the position of m1

因此在圖3中,柔度系數(shù)[21]

根據(jù)式(16),圖4中在水平集中力F的作用下，集中質(zhì)量位置處的水平撓度為：

(17)

根據(jù)式(14), 圖4中m2以下位置處的撓度函數(shù)為

(18)

圖4 計算m2 位置處的柔度系數(shù)Fig.4 Computation of flexibility coefficient at the position of m2

(19)

將式(17)和(18)代入式(19)，得：

(20)

因此在圖4中，柔度系數(shù)[21]

(21)

(22)

上述柔度系數(shù)確定后，即可建立集中質(zhì)點在慣性力作用下的位移方程，進而根據(jù)柔度法計算出圖3和圖4中2個自由度動力系統(tǒng)的自振頻率。

2 算例

圖1中，風力發(fā)電機的塔和基礎都假定為彈性模量[18]E=2.1×1011N/m2的鋼材制成的套筒。

套筒直徑d=4 000 mm,泥面線以上高度h=90 m,樁入土110 m。參考文獻[19]中水下樁的樁徑和壁厚的比值，圖1中套筒的壁厚采用52 mm。根據(jù)環(huán)形截面慣性矩計算公式得I=1.256 m4。因此鋼套筒的環(huán)形截面抗彎剛度EI=2.1×1011N/m2×1.256 m4=2.6×1011N·m2。

通過計算，按照鋼材密度7.85×103kg/m3,計算出泥面線以上套筒質(zhì)量M0=4.5×105kg。機艙重量1.0×105kg,輪轂重量0.2×105kg，葉片重量3×0.1×105kg[19]，則機艙、輪轂和葉片的總重量m0=1.5×105kg。

上述設計參數(shù)都是根據(jù)相關(guān)參考文獻的設計值組合得到的，與文獻[13]的數(shù)據(jù)相同，但其自振頻率分析沒有考慮樁基礎以及樁土相互作用。這些設計參數(shù)是否滿足初步設計要求，需根據(jù)現(xiàn)場記錄的荷載統(tǒng)計數(shù)據(jù)進行套筒的彎曲正應力和剪應力計算，然后根據(jù)材料力學的強度理論進行判斷，最后逐步進行精確有限元分析和優(yōu)化設計。

假設地基土的水平反力系數(shù)[16]k=200 kPa，即地基土的水平反力系數(shù)比極軟淤泥的水平反力系數(shù)還小，根據(jù)集中質(zhì)量法及算例參數(shù)可知：

m1=0.25*M0+m0=2.625*105kg，

m2=0.5*M0=2.25*105kg,

h=90m，EI=2.6×1011N·m2，β=0.030m-1

根據(jù)柔度法[21]，即可計算出圖1泥面線以上結(jié)構(gòu)的自振頻率(計算過程見附錄的fortran源程序)：

第一頻率為1.141 Hz,第二頻率為8.421 Hz。

3 討論

3.1 與其他算法計算的自振頻率的比較

對于參數(shù)完全相同的結(jié)構(gòu)(圖1), 當不考慮樁土相互作用時，牛文杰等[13]使用集中質(zhì)量法，采用剛度法對近海風力發(fā)電機的塔和泥面線上基礎按照2個自由度的動力系統(tǒng)進行自由振動分析, 得到自振頻率，經(jīng)過驗算，文中M應為(0.25*M0+m0),m應為0.5*M0。重新根據(jù)剛度法計算，第一頻率是2.050 Hz,第二頻率是13.559 Hz；根據(jù)柔度法計算，第一頻率是1.928 Hz,第二頻率是13.597 Hz。在地面按照9.8sin(4.626t)的規(guī)律運動時(單位m)，質(zhì)點M的位移幅值AM=-15.084 m，質(zhì)點m的位移幅值Am=-4.505 m。

而本文考慮樁土相互作用，計算出理論上更加精確的自振頻率：第一頻率1.141 Hz,第二頻率8.421 Hz。

3.2 影響自振頻率的因素

在自振頻率的計算過程中，通過觀察知泥面線以上風機結(jié)構(gòu)的自振頻率取決于樁的慣性矩EI、樁徑d、地基土的水平反力系數(shù)k、樁基礎入土深度、風機在泥面線以上的高度h及泥面線以上風機的質(zhì)量分布。

3.3 改進方向

在自振頻率的計算中，文獻[9,22]考慮了葉片、機頭等部件的重力作用，這使得計算結(jié)果更加符合實際。在考慮重力作用的情況下，根據(jù)柔度法，單樁結(jié)構(gòu)物的水平撓度會增大。實際上也應該考慮塔筒、樁基礎的重力作用，但其計算過程會更加復雜。

另外，對風機整體，當葉片的回轉(zhuǎn)直徑達到91.6m時，若風力發(fā)電機的葉片、機頭等部件簡化為一集中質(zhì)量塊進行自振頻率的計算，其結(jié)果是否合理還需實驗檢驗。風機的支撐結(jié)構(gòu)截面一般都是從下到上逐漸變化的(變小)，還需在今后計算時進行改進。

樁實際上是有一定長度的，樁端的支撐條件對風機體系的自振頻率存在一定影響。本文將樁作為半無限長考慮是一種工程假定，實際近海風機樁基礎的入土深度需要根據(jù)近海地質(zhì)資料、海洋環(huán)境荷載和正常使用極限等因素確定。實際工程中地基大多數(shù)是由不同巖土層組成的(例如海砂-軟土互層)，因此地基對樁的橫向支撐剛度也必然是隨深度變化的，本文沒有考慮這一變化，與實際情況不符，需在今后加以考慮，例如使用p-y曲線而非張氏法考慮樁土相互作用。樁周土在復雜動應力[23]作用下會發(fā)生應力-應變關(guān)系的變化，進而影響單樁風機的自振頻率，也需在今后計算中加以考慮[24]。

4 結(jié)論

以泥面線以上的風機上部結(jié)構(gòu)為研究對象，給出快速估算近海單樁風機自振頻率的計算方法。按照集中質(zhì)量法，把泥面線以上的風機結(jié)構(gòu)按照多自由度動力系統(tǒng)進行分析，計算出更加精確的自振頻率值。

觀察自振頻率的計算過程，可知泥面線以上風機結(jié)構(gòu)的自振頻率取決于樁的慣性矩EI、樁徑d、地基土的水平反力系數(shù)k、樁基礎入土深度、風機在泥面線以上的高度h及泥面線以上風機結(jié)構(gòu)的質(zhì)量分布。

)

[1] 趙先民.水平軸風力機動力特性和流場的數(shù)值模擬[D].上海:同濟大學,2007.

ZHAOXian-min.DynamicCharacteristicsofaHorizontalAxisWindTurbineandNumericalSimiluationofFlowField[D].Shanghai:TongjiUniversity,2007.(inChinese)

[2] Wang F,Chen Q,Yu G-C.Research on Tower Stiffness of Large Scale Wind Turbine Unit[J].New Energy and Technology,1997,20(6):38-39.

[3] 湯煒梁，袁奇，韓中合.風力機塔筒抗臺風設計[J].太陽能學報,2008,29(4):422-427.

TANG Wei-liang,YUAN Qi,HAN Zhong-he.Withstanding Typhoon Design of Wind Turbine Tower[J].Acta Energiae Solaris Sinica,2008,29(4):422-427.(in Chinese)

[4] 王懿，段夢蘭，尚景宏，等.海上風機基礎結(jié)構(gòu)力學分析[J].中國海洋平臺，2009,24(4):14-20.

WANG Yi，DUAN Meng-lan，SHANG Jing-hong，et al. Analysis of Mechanics for Offshore Wind Turbines Founadtion[J].China Offshore P1atform,2009,24(4):14-20.(in Chinese)

[5] Merle C Potter,David C Wiggert.Mechanics of Fluids[M].Beijing:China Machine Press,2003.

[6] 尤漢強,楊敏.循環(huán)荷載下海洋風機單樁基礎簡化分析模型[J].巖土工程學報,2010,32(增刊2):13-16.

YOU Han-qiang,YANG Min.Simplified Model for Single-pile of Offshore Wind Turbine Under Cyclic Loading[J].Chinese Journal of Geotechnical Engineering,2010,32(Supp.2):13-16.(in Chinese)

[7] Jan van der Tempel,David-pieter Molenaar.Wind Turbine Structural Dynamics——A Review of the Principles for Modern Power Generation,Onshore and Offshore[J].Wind Engineering,2002,26(4):211-220.

[8] Subhamoy Bhattacharya,James A Cox,Domenico Lombardi,et al.Dynamics of Offshore Wind Turbines Supported on Two Foundations[J].Geotechnical Engineering,2013,166(2):159-169.

[9] Subhamoy Bhattacharya，Domenico Lombardi，David Muir Wood.Similitude Relationships for Physical Modelling of Monopile-supported Offshore Wind Turbines[J].International Journal of Physical Modelling in Geotechnics,2011,11(2)：58-68.

[10] Zaaijerm M B,Foundation Modeling to Assess Dynamic Behaviour of Offshore Wind Turbines[J].Applied Ocean Research,2006,28(1):45-57.

[12] Suresh R Dash,Subhamoy Bhattacharya,Anthony Blakeborough.Bending-buckling Interaction as a Failure Mechanism of Piles in Liquefiable Soils[J].Soil Dynamic and Earthquake Engineering,2010,30(1/2):32-39.

[13] 牛文杰, 王振宇, 李洪然.海上風機-基礎水平振動的多自由度分析[J].土木建筑與環(huán)境工程, 2011, 33 (增刊1):121-124.

NIU Wen-jie,WANG Zhen-yu,LI Hong-ran.The Multiple-degree-of-freedom Analysis of Offshore Wind Turbine and Foundation in Horizontal Vibration[J].Journal of Civil, Architectural & Environmental Engineering,2011,33(Supp.1):121-124.(in Chinese))

[14] Bhattacharya S,Adhikari S.Experimental Validation of Soil-structure Interaction of Offshore Wind Turbines[J].Soil Dynamics and Earthquake Engineering,2011,31(5):805-816.

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

YANG Jin,LIU Shu-jie.The Application of ANSYS in Offshore Oil Engineering[M].Beijing:Petroleum Industry Press,2010.

[16] 劉金礪.樁基礎設計與計算[M].北京：中國建筑工業(yè)出版社，1990.

LIU Jin-Li.Design and Calculation of Pile Foundation[M].Beijing:China Architecture & Building Press,1990.(in Chinese)

[17] 盧世深,林亞超.樁基礎的計算和分析[M].北京:人民交通出版,1987.

LU Shi-shen,LIN Ya-chao.Calculation and Analysis of Pile Foundation[M].Beijing:China Communications Press,1987.

[18] Abdel-rahman K,Achmus M.Finite Element Modelling of Horizontally Loaded Monopile Foundations for Offshore Wind Energy Converters in Germany[C]// International Symposium on Frontiers in Offshore Geotechnics.Perth,Australia,2005:309-396.

[19] 郇彩云.海上樁式風機基礎結(jié)構(gòu)設計與研究[D].大連：大連理工大學,2009.

XUAN Cai-yun.Design,Research of Pile-style Foundation Structures for Offshore Wind Turbine[D].Dalian:Dalian University of Technology,2009.(in Chinese)

[20] 龍馭球,包世華.結(jié)構(gòu)力學[M].北京：高等教育出版社,1994.

LONG Yu-qiu,BAO Shi-hua.Structural Mechanics[M].Beijing:Higher Education Press,1994.(in Chinese)

[21] 周競歐,朱伯欽,許哲明.結(jié)構(gòu)力學[M].上海:同濟大學出版社,2004.

ZHOU Jing-ou,ZHU Bo-qin,XU Zhe-ming.Structural Mechanics[M].Shanghai:Tongji University Press,2004.(in Chinese)

[22] 余璐慶.海上風機桶形基礎安裝與支撐結(jié)構(gòu)動力特性研究[D].杭州：浙江大學,2014.

YU Lu-qing.Study on the Installation Behavior of Suction Caisson and the Dynamic Properties of Offshore Wind Turbine Structures[D].Hangzhou:Zhejiang University,2014.(in Chinese)

[23] 左照坤,童朝霞.近海風機單樁基礎樁周土應力特征分析[J].地震工程學報,2014,36(3):549-554.

ZUO Zhao-kun,TONG Zhao-xia.Stress Characteristics of Soil around the Pile of a Monopile Foundation in Offshore Wind Turbines[J].China Earthquake Engineering Journal,2014,36(3):549-554.(in Chinese)

[24] 余璐慶，王立忠，Bhattacharya.海上風機支撐結(jié)構(gòu)動力特性模型試驗研究[J].地震工程學報,2014,36(4):797-803.

YU Lu-qing,WANG Li-zhong,Bhattacharya.Scaled Model Test Study of the Dynamic Behavior of an Offshore Wind Turbine Support Structure[J].China Earthquake Engineering Journal,2014,36(4):797-803.(in Chinese)

附錄 計算圖1泥面線以上結(jié)構(gòu)自振頻率的FORTRAN源程序

REAL M1,M2,H,EI,BETA,F11,F21,F22,F12,LMDA1,LMDA2,OMIGA1,OMIGA2

M1=2.625*100000

M2=2.25*100000

H=90

EI=2.6*100000*100000*10

BETA=0.030

F11=-H*H*H/(6*EI)+H*H*H/(2*EI)+((1+BETA*H)/(2*EI*BETA*BETA)+H/(2*EI*BETA))*H+(1+BETA*H)/(2*EI*BETA*BETA*BETA)

F21=-H*H*H*0.5*0.5*0.5/(6*EI)+H*0.5*H*0.5*H/(2*EI)

F21=F21+((1+BETA*H)/(2*EI*BETA*BETA)+H/(2*EI*BETA))*0.5*H+(1+BETA*H)/(2*EI*BETA*BETA*BETA)F22=H*H*H*0.5*0.5*0.5/(3*EI)+((1+BETA*0.5*H)/(2*EI*BETA*BETA)+0.5*H/(2*EI*BETA))*0.5*H+(1+BETA*0.5*H)/(2*EI*BETA*BETA*BETA)

F12=H*H*H*0.5*0.5*0.5/(3*EI)+H*0.5*H*0.5*0.5*H/(2*EI)

F12=F12+((1+BETA*0.5*H)/(2*EI*BETA*BETA)+0.5*H/(2*EI*BETA))*H+(1+BETA*0.5*H)/(2*EI*BETA*BETA*BETA)

LMDA1=F11*M1+F22*M2

LMDA1=LMDA1+SQRT((F11*M1+F22*M2)*(F11*M1+F22*M2)-4*(F11*F22-F12*F12)*M1*M2)

LMDA1=LMDA1/2

LMDA2=F11*M1+F22*M2

LMDA2=LMDA2-SQRT((F11*M1+F22*M2)*(F11*M1+F22*M2)-4*(F11*F22-F12*F12)*M1*M2)

LMDA2=LMDA2/2

OMIGA1=SQRT(1/LMDA1)

OMIGA2=SQRT(1/LMDA2)

WRITE(*,*)F11

WRITE(*,*)F21

WRITE(*,*)F22

WRITE(*,*)F12

WRITE(*,*)LMDA1

WRITE(*,*)LMDA2

WRITE(*,*)OMIGA1

WRITE(*,*)OMIGA2

end program

Natural Frequency of Mono-pile Wind Turbine Considering Interaction between Pile and Elastic Foundation

NIU Wen-jie

(SchoolofMechanicsandEngineering,LiaoningTechnicalUniversity,Fuxin123000,Liaoning,China)

In this paper, we analyze the mono-pile wind turbine, comprising blade bodies, a hub, a cabin, a conical shell tower, a transition piece, and a mono-pile foundation. Our study objective was to determine the natural frequency of the mono-pile wind turbine above the mudline. We assumed the conical shell tower, transition piece, and mono-pile foundation to be a vertical steel tube (VST) with a uniform cross section and with the same elastic modulus. We used the lumped-mass method to simplify the mono-pile wind turbine above the mudline as two concentrated masses. One concentrated mass at the VST top is the summation of the mass of the blade bodies, hub, cabin and one-fourth the mass of the VST above the mudline. The other concentrated mass at the middle of the VST is one-half the mass of the VST above the mudline. Taking the mono-pile foundation below the mudline as a semi-infinite long pile, we used Zhang’s method to obtain the horizontal stiffness and rotational stiffness of the pile top. In this way, we can determine the deflection of the VST under horizontal load at any point of the wind turbine above the mudline. The flexibility coefficient of the unit horizontal load at any position of the VST can be determined accordingly. Then, using both the flexibility and lumped-mass methods, we can determine the natural frequency of the mono-pile wind turbine above the mudline. Lastly, with an example, we retrieved the natural frequency value of the mono-pile wind turbine above the mudline. We then compared the analysis result with the mono-pile wind turbine natural frequency obtained without considering the interaction between the pile and elastic groundwork, and found little difference between them. The new theory we propose in this paper is more accurate theoretically and should be validated experimentally.

mono-pile wind turbine; natural frequency; interaction between pile and elastic foundation; flexibility method; lumped-mass method

2016-04-26

遼寧工程技術(shù)大學博士科研啟動基金(11-415)

牛文杰(1982-)，男，河南漯河人，副教授，博士，從事近海單樁風機工程中巖土力學、結(jié)構(gòu)力學和流體力學問題的研究。E-mail: nwj1982@sohu.com。

TU352.2; TU411.3; TU473.1

A

1000-0844(2016)05-0713-07

10.3969/j.issn.1000-0844.2016.05.0713