基于整體耦合模型的海上風(fēng)機結(jié)構(gòu)樁-土相互作用影響研究

柏延強， 宋雨果， 邵慶梧， 齊曉海， 王文華， 李 昕

（1. 中廣核研究院有限公司， 廣東 深圳 518000； 2. 大連理工大學(xué) 建設(shè)工程學(xué)部抗震研究所， 遼寧 大連116024）

0 引言

海上發(fā)電機組由于其本身所處的復(fù)雜的環(huán)境條件，長期受到風(fēng)、浪、流、海冰、地震等環(huán)境荷載的聯(lián)合作用，因此，多極限狀態(tài)下海上風(fēng)電機組基礎(chǔ)結(jié)構(gòu)的安全成為保障海上風(fēng)電機組平穩(wěn)、安全運行的重要影響因素。 現(xiàn)階段，主要采用解耦或半整體分析方法開展隨機荷載作用下固定式海上風(fēng)機結(jié)構(gòu)動力特性分析及安全評價。

夏露[1]將作用于風(fēng)機轉(zhuǎn)子系統(tǒng)的空氣動力荷載等效為作用于基礎(chǔ)結(jié)構(gòu)的風(fēng)機荷載，基于海上風(fēng)機基礎(chǔ)結(jié)構(gòu)有限元模型，采用半整體方法研究了隨機風(fēng)浪作用下海上風(fēng)機結(jié)構(gòu)的運動響應(yīng)變化規(guī)律。 文獻[2]～[4]采用半整體方法研究了p-y 曲線、m 法、接觸單元法和耦合彈簧單元等不同樁土相互作用模型對于固定式海上風(fēng)機基礎(chǔ)結(jié)構(gòu)響應(yīng)特性的影響。 文獻[5]通過將樁-土相互作用等效為線性樁基超單元矩陣，通過修正超單元矩陣剛度，探討了地基土壤參數(shù)對于海上風(fēng)機結(jié)構(gòu)承載力評價指標(biāo)的影響。

綜合上述研究，半整體方法的主要技術(shù)特點為風(fēng)電機組結(jié)構(gòu)與基礎(chǔ)結(jié)構(gòu)的分離式建模和計算， 該方法所采用的分離式建模直接忽略了風(fēng)電機組與基礎(chǔ)結(jié)構(gòu)之間的耦合效應(yīng)， 以及氣動阻尼等效應(yīng)對于波浪、海流、海冰等荷載作用下基礎(chǔ)結(jié)構(gòu)響應(yīng)的影響。 需要指出，上述荷載、荷載組合和樁基邊界條件同時作用于風(fēng)機結(jié)構(gòu)， 由此所產(chǎn)生的耦合效應(yīng)對于海上風(fēng)機結(jié)構(gòu)動力特性和動力響應(yīng)具有顯著影響， 而現(xiàn)階段廣泛采用的半整體方法在一定程度上簡化或忽略了上述耦合效應(yīng)及其影響。 綜合考慮國內(nèi)外研究進展， 本文通過對FAST v8[6]進行二次開發(fā)，建立了包含轉(zhuǎn)子機艙組件-塔筒-基礎(chǔ)結(jié)構(gòu)-樁-土相互作用的海上風(fēng)機整體耦合數(shù)值仿真模型， 開展典型設(shè)計工況下海上風(fēng)機整體結(jié)構(gòu)動力響應(yīng)計算。 依據(jù)海上風(fēng)機結(jié)構(gòu)時頻域運動響應(yīng)變化規(guī)律，揭示并探討樁-土相互作用對于海上風(fēng)機結(jié)構(gòu)響應(yīng)的影響機理。

1 整體耦合模型

FAST 軟件SubDyn 模塊[7]將基礎(chǔ)結(jié)構(gòu)節(jié)點分為3 類，分別為過渡節(jié)點（Interface node）、樁基節(jié)點（Pile node）和內(nèi)部節(jié)點（Internal node），其中前兩類節(jié)點統(tǒng)稱為邊界節(jié)點。 圖1 所示為FAST v8中海上風(fēng)機整體耦合模型圖。

圖1 FAST v8 中海上風(fēng)機整體耦合模型圖Fig.1 Fully coupled analysis model of offshore wind turbine in FAST v8

SubDyn 模塊默認的基礎(chǔ)結(jié)構(gòu)樁基節(jié)點邊界條件為固端約束[圖1（a）]，該類邊界條件為

本文基于式（4）所示修正后基礎(chǔ)結(jié)構(gòu)耦合運動方程， 通過對FAST v8 軟件SubDyn 模塊進行二次開發(fā)，建立了風(fēng)、浪聯(lián)合作用下包含氣彈-伺服-水動-彈性-PSI 的海上風(fēng)機整體耦合分析模型。該模型能夠更為準(zhǔn)確地揭示風(fēng)、浪作用下海上風(fēng)機結(jié)構(gòu)耦合響應(yīng)特性。

2 樣本風(fēng)機及設(shè)計工況組合

2.1 樣本風(fēng)機基本參數(shù)

選取NREL OC3 單樁基礎(chǔ)海上風(fēng)機作為本次研究的樣本風(fēng)機[8]，該風(fēng)機輪轂中心高度為90.0 m，塔筒高度為87.6 m，風(fēng)機部分基本參數(shù)[9]如表1所示。單樁基礎(chǔ)直徑為6.0 m，壁厚為0.06 m，設(shè)計水深為20.0 m，風(fēng)機樁長為36.0 m，樁徑為6.0 m，沿樁長方向土壤參數(shù)如圖2 所示。

表1 樣本風(fēng)機基本參數(shù)Table 1 Basic parameters of reference OWT

圖2 樣本風(fēng)機土壤參數(shù)Fig.2 Soil profile of reference OWT

2.2 設(shè)計工況及荷載參數(shù)

依據(jù)海上風(fēng)機設(shè)計規(guī)范DNV GL-ST-0437[10]和我國東南海域某風(fēng)電場實測海洋環(huán)境資料，選取海上風(fēng)機典型設(shè)計工況及海洋環(huán)境參數(shù) （表2）。

表2 基本荷載組合工況Table 2 Selected load cases

依據(jù)表2 所列風(fēng)速參數(shù)，采用IEC Kaimal 譜[11]擬合得到不同工況下隨機風(fēng)速時程[圖3（a）]。 同時，采用JONSWAP 譜[12]，依據(jù)所選取波浪參數(shù)生成隨機波高時程[圖3（b）]，所對應(yīng)傅里葉譜如圖3（c）所示。

圖3 擬合得到的隨機環(huán)境荷載Fig.3 Synthesized stochastic environmental load conditions

3 考慮樁-土相互作用的海上風(fēng)機結(jié)構(gòu)耦合響應(yīng)特性研究

3.1 樣本風(fēng)機整體結(jié)構(gòu)動力特性

通過將白噪聲激勵施加于所建立的耦合數(shù)值仿真模型，結(jié)合小阻尼體系假定，選取各節(jié)點運動響應(yīng)傅里葉譜得到了如表3 所示的前兩階頻率，所對應(yīng)的各階振型如圖4 所示。

由表3 可知， 相比于固定邊界條件， 在考慮樁-土相互作用后， 海上風(fēng)機結(jié)構(gòu)自振頻率明顯降低，尤其對于二階頻率。 由圖4 對比可知：不同邊界條件下， 海上風(fēng)機耦合模型的一階振型變化規(guī)律基本相同， 最大振型位移幅值均位于塔筒位置；對于二階振型，由于樁-土相互作用的影響，振型位于幅值所在高度發(fā)生明顯變化； 對于固定邊界和彈性邊界耦合模型， 二階振型位移幅值所在高度分別為48.8 m 和44.1 m。

圖4 海上風(fēng)機整體結(jié)構(gòu)前兩階振型Fig.4 The first two mode shapes of OWT

表3 海上風(fēng)機整體結(jié)構(gòu)前兩階頻率Table 3 First two natural frequencies of OWT

3.2 隨機風(fēng)浪作用下海上風(fēng)機結(jié)構(gòu)耦合響應(yīng)機理

3.2.1 塔筒結(jié)構(gòu)運動響應(yīng)耦合機理

基于建立的整體耦合計算模型開展隨機風(fēng)浪聯(lián)合作用下的海上風(fēng)機整體結(jié)構(gòu)動力響應(yīng)， 對比固定邊界、 彈性邊界整體耦合計算模型響應(yīng)特性差異。

在隨機風(fēng)浪作用下，固定邊界、彈性邊界整體耦合模型塔頂位移時頻域時程如圖5 所示。

圖5 各設(shè)計工況下塔頂位移時程響應(yīng)Fig.5 Histories of tower top displacement under load cases

由圖5 可知，在隨機風(fēng)浪作用下，固定邊界、彈性邊界整體耦合模型所得塔頂位移時程相應(yīng)差異明顯。 選取圖4 所示塔筒各節(jié)點展開塔筒截面位移、彎矩時程響應(yīng)標(biāo)準(zhǔn)差對比，對比結(jié)果如圖6所示。

圖6 塔筒結(jié)構(gòu)時程響應(yīng)標(biāo)準(zhǔn)差統(tǒng)計值對比Fig.6 Comparisons of standard deviations of the histories of tower responses

由圖6 可知：在所選取設(shè)計工況下，基于包含樁-土相互作用耦合數(shù)值仿真模型所得塔筒各截面位移及彎矩時程響應(yīng)標(biāo)準(zhǔn)差明顯大于忽略樁-土相互作用耦合模型的計算結(jié)果； 對于塔頂截面位移，隨著塔筒高度的增加，二者差異明顯增大；樁-土相互作用對于塔筒各截面彎矩時程標(biāo)準(zhǔn)差的影響，以工況LC3 下塔基截面彎矩為例，此時有無樁-土相互作用， 可造成該截面彎矩時程響應(yīng)標(biāo)準(zhǔn)差的相對誤差為8.68%。 在工況LC1～LC3作用下，依據(jù)式（7）計算得到的有無樁-土相互作用耦合模型的塔頂位置位移時程響應(yīng)標(biāo)準(zhǔn)差相對誤差Δ 分別為3.68%，6.57%和12.57%。式中：δ* 為基于包含樁-土相互作用耦合模型所得結(jié)構(gòu)時程響應(yīng)標(biāo)準(zhǔn)差；δ 為忽略樁-土相互作用耦合模型所得結(jié)構(gòu)時程響應(yīng)標(biāo)準(zhǔn)差。

圖7 所示為各設(shè)計工況下塔頂位移頻域響應(yīng)曲線。

圖7 各設(shè)計工況下塔頂位移頻域響應(yīng)Fig.7 Fourier amplitudes of tower top displacement under selected load cases

由圖7 可知：相比于固定邊界，對于彈性邊界耦合模型，隨機風(fēng)、浪作用下整體結(jié)構(gòu)基頻及風(fēng)、浪荷載頻率的影響顯著增加； 樁-土相互作用對于海上風(fēng)機支撐系統(tǒng)耦合運動響應(yīng)特性存在顯著影響。與固定邊界條件相比，在隨機風(fēng)、浪作用下，考慮樁-土相互作用后的海上風(fēng)機塔筒結(jié)構(gòu)運動響應(yīng)明顯增大。

3.2.2 基礎(chǔ)結(jié)構(gòu)耦合響應(yīng)機理

選取風(fēng)、 浪作用下不同邊界條件單樁基礎(chǔ)海上風(fēng)機整體耦合模型基底傾覆力矩時程響應(yīng)及相應(yīng)統(tǒng)計值進行對比，對比結(jié)果如表4 和圖8 所示。

表4 各設(shè)計工況下基底傾覆力矩時程響應(yīng)統(tǒng)計Table 4 Statistics of bending moments under load cases

圖8 各設(shè)計工況下基底傾覆力矩時程響應(yīng)Fig.8 Histories of bending moment under selected load cases

由時程響應(yīng)標(biāo)準(zhǔn)差對比可驗證邊界條件對于海上風(fēng)機整體結(jié)構(gòu)運動響應(yīng)的影響。 由統(tǒng)計值對比可知，隨著來流風(fēng)速和設(shè)計波高的增大，由邊界條件差異所造成的基底傾覆力矩時程響應(yīng)標(biāo)準(zhǔn)差相對誤差由3.48%增加至18.90%。

圖9 所示為各設(shè)計工況下基底傾覆力矩頻域響應(yīng)曲線。由基底傾覆力矩傅里葉譜對比可知，隨著將海上風(fēng)機基礎(chǔ)樁基點固定邊界條件修正為彈性邊界條件，整體結(jié)構(gòu)基頻和風(fēng)、浪頻率對于基底傾覆力矩的影響愈發(fā)明顯， 由此造成了上述環(huán)境荷載作用下海上風(fēng)機結(jié)構(gòu)運動時程響應(yīng)統(tǒng)計值的增大。

圖9 各設(shè)計工況下基底傾覆力矩頻域響應(yīng)Fig.9 Fourier amplitudes of bending moment under selected load cases

4 結(jié)論

本文依據(jù)線性化樁-土相互作用模型， 通過對時域耦合數(shù)值仿真軟件FAST v8 的基礎(chǔ)結(jié)構(gòu)模塊SubDyn 進行二次開發(fā)， 建立了包含轉(zhuǎn)子機艙組件-塔筒-基礎(chǔ)結(jié)構(gòu)-樁基礎(chǔ)的海上風(fēng)機整體耦合計算模型。以單樁基礎(chǔ)海上風(fēng)機為研究對象，開展了隨機荷載作用下固定邊界和彈性邊界海上風(fēng)機耦合響應(yīng)特性研究，得到以下結(jié)論。

①在考慮樁-土相互作用后， 相比于固定邊界，海上風(fēng)機結(jié)構(gòu)自振頻率明顯降低，其中結(jié)構(gòu)基頻相對誤差為13%， 二階頻率相對誤差可達到35%，對于二階頻率影響尤為顯著。 同時，整體結(jié)構(gòu)二階振型位移幅值分布亦存在明顯差異。

②隨機風(fēng)、浪作用下，由于有無樁-土相互作用影響， 整體耦合分析方法得到的海上風(fēng)機結(jié)構(gòu)響應(yīng)時程存在顯著差異， 以塔基截面彎矩和基底傾覆力矩為例， 此時時程響應(yīng)標(biāo)準(zhǔn)差最大相對誤差可達到31.59%和18.90%。

③除時程響應(yīng)差異外， 由于邊界條件差異的影響，相對于固定邊界，在風(fēng)、浪荷載作用下整體結(jié)構(gòu)頻率及環(huán)境荷載頻率對于彈性邊界耦合模型所得結(jié)構(gòu)運動響應(yīng)的影響顯著增加。