李爭霖　胡金鵬

摘? ? 要：海上風(fēng)電機組高承臺復(fù)合樁基礎(chǔ)為同時存在大中尺度承臺和小尺度桿件基樁的復(fù)合結(jié)構(gòu)形式，按現(xiàn)有設(shè)計規(guī)范對此類結(jié)構(gòu)進行波浪水動力載荷計算存在較大誤差。為準確計算復(fù)合樁基礎(chǔ)小尺度桿件所受的波浪載荷，不僅要考慮入射波浪場的作用， 還應(yīng)考慮臨近大中尺度承臺結(jié)構(gòu)引起的波浪繞射及反射對小尺度桿件波浪載荷的影響。本文基于VOF法在Fluent中構(gòu)建波浪對風(fēng)電復(fù)合樁基礎(chǔ)作用的三維波浪數(shù)值水槽，首先模擬規(guī)則波單獨沖擊下部樁柱的過程，驗證模型計算小尺度桿件波浪載荷的有效性;然后通過數(shù)值計算得到作用于復(fù)合樁基礎(chǔ)上的入射波浪力、繞射波浪力以及總力。結(jié)果顯示，在一定的波浪條件下，繞射場作用在樁基礎(chǔ)上的波浪力隨著上部承臺尺寸的增大而增大。本文還進一步研究了極端波浪對小尺度樁柱波浪載荷的影響，有關(guān)的研究結(jié)果可為具有復(fù)合結(jié)構(gòu)形式的海洋工程基礎(chǔ)中的小尺度桿件上極端波浪載荷的準確計算提供一定的參考。

關(guān)鍵詞：高承臺樁基礎(chǔ);波浪繞射;極端波浪;數(shù)值模擬;波浪載荷

中圖分類號：P731.22 ?? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標識碼：A

Abstract： Because of the composite structure of large and medium scale pile cap and small scale bar pile， there is a large error in the wave hydrodynamic load calculation of high-cap composite pile foundation of offshore wind turbine according to the current design code. In order to accurately calculate the wave load on the small-scale members of the composite pile foundation， we should not only consider the effect of the incident wave field， but also consider the influence of wave diffraction and reflection on the wave load of small-scale members caused by the structure near the large and medium scale cap. This paper constructs a three-dimensional wave numerical flume of wave action on the composite pile foundation in Fluent based on VOF method. Firstly， it simulates the process of regular wave impact on lower pile column alone， verifies the validity of the model to calculate the wave load of small scale members， and then obtains the incident wave force， diffraction wave force and total force acting on the composite pile foundation by numerical calculation. The results show that under certain wave conditions， the wave force acting on the pile foundation increases with the increase of the size of the upper cap. In this paper， the effect of extreme wave on wave load of small scale pile column is further studied.

Key words： High-cap composite pile foundation; Wave diffraction; Extreme wave; Numerical simulation ;Wave load

1? ? ?引言

海上風(fēng)電機組高承臺群樁基礎(chǔ)是東南沿海地區(qū)廣泛應(yīng)用的新型風(fēng)電基礎(chǔ)形式。其中，波浪水動力載荷在結(jié)構(gòu)設(shè)計中起到了決定性作用。樁基礎(chǔ)所受的波浪載荷一般都根據(jù)莫里森公式計算[2]，但在極端波浪條件下，承臺與波長的比值往往接近甚至超過臨界值，此時莫里森公式存在較大的誤差;另外，較大的承臺結(jié)構(gòu)使得波浪場發(fā)生改變，生成影響樁基礎(chǔ)的入射波浪場、繞射波浪場，這將大大提高波浪水動力載荷計算的難度[3]。

目前，國內(nèi)外學(xué)者對波浪繞射場的樁柱荷載進行了大量的研究。MacCamy 和Fuchs[4]利用Laplace方程求解大直徑直立圓柱的入射波浪力和繞射波浪力;Linton和Evans [5]提出了規(guī)則波作用下樁柱繞射波浪力的計算方法;耿寶磊[6]計算了繞射波浪場中小直徑圓柱的波浪載荷;姜勝超[7]在線性理論假設(shè)下建立了有限水深水下垂直圓柱波浪繞射的解析解?？紤]到高承臺風(fēng)電基礎(chǔ)的實際應(yīng)用廣泛，其相關(guān)水動力問題需要作進一步研究，尤其是受到承臺大中尺度結(jié)構(gòu)影響下小直徑結(jié)構(gòu)上的波浪力研究，對于工程設(shè)計和后期的安全生產(chǎn)具有理論和實際意義。

本文以大承臺復(fù)合小樁柱的組合結(jié)構(gòu)為研究對象，考慮繞射場帶來的影響，研究組合結(jié)構(gòu)樁柱上的波浪荷載。本文基于Fluent軟件，以N-S方程作為控制方程，以VOF方法捕捉流體自由表面，采用UDF二次開發(fā)功能實現(xiàn)仿物理造波，利用PISO速度耦合方式及附加源項法實現(xiàn)阻尼層消波，減少遠方墻壁的反射，從而建立了能夠產(chǎn)生穩(wěn)定線性規(guī)則波的三維數(shù)值波浪水槽，通過數(shù)值計算得到作用到小尺度桿件的波浪力;同時還分析了極端波浪畸形波對繞射作用的影響。本文的研究方法，可以作為高承臺群樁復(fù)合式海洋工程結(jié)構(gòu)物考慮極端波浪力設(shè)計的一個有益的探索性研究。

2? ? ?數(shù)學(xué)模型原理

2.1? ?數(shù)值波浪水槽的基本控制方程

本文所建立的數(shù)值波浪水槽中的流體假設(shè)為不可壓縮粘性流體。

2.2? ? 自由液面處理方法

波浪問題屬于具有自由表面的兩相流問題，其最大特點是存在隨時間不斷變化的水-氣交界面，所以對自由液面的捕捉是處理波浪問題的關(guān)鍵。FLUENT軟件中采用VOF法來追蹤流體的自由液面。

VOF法由于具有可以處理自由面重構(gòu)、容易擴展到三維、計算時間短、占用存儲小等優(yōu)點，成為一種廣泛應(yīng)用的計算自由表面高度的方式。VOF法的核心思想是計算第q相流體體積在單元中的占據(jù)情況，即占比函數(shù)（F函數(shù)）。若q=1，表示整個單元都是q相流體;若aq=0，表示單元內(nèi)為空;若 0

由于用一般形式的差分法計算流體體積函數(shù)有一定困難，因此VOF法常用各單元之間的兩相關(guān)系表征函數(shù)的輸運情況：

2.3? ?造波方法

本文在二維及三維數(shù)值波浪水槽中進行模擬造波是根據(jù)物理水池里的推板造波方法，將水槽起動側(cè)設(shè)置為墻邊界，邊界的運動通過UDF編程定義，即相當于造波板。運動邊界附近網(wǎng)格點信息的動態(tài)傳遞，則通過動態(tài)層（Layering）方法，從而實現(xiàn)水槽的波浪啟動。運動邊界的速度信號需要通過目標波浪的波面方程推導(dǎo)得來。

水力傳遞函數(shù) Tr是一個多種周期都滿足的造波邊界附近波浪振幅 a和邊界振幅 e 的關(guān)系式，對于推板式造波方法，傳遞函數(shù)Tr 的表達式為：

3? ? ?數(shù)值模型建立

3.1? ?幾何模型

本文的研究主要集中在大尺度結(jié)構(gòu)對小尺度結(jié)構(gòu)波浪力的影響，考慮到三維模型的計算復(fù)雜性，為簡化起見本文的結(jié)構(gòu)物模型僅考慮了單個基樁，而更為完整的結(jié)構(gòu)物模型將在下一步研究中開展。

參考廣東省電力設(shè)計研究院目前設(shè)計的粵東深水風(fēng)電場（水深50 m，此水深條件為近海極端波浪多發(fā)區(qū)域）中海上風(fēng)電樁基礎(chǔ)的設(shè)計條件，本文的三維數(shù)值波浪水槽模型的幾何示意圖見圖1。水槽設(shè)計尺寸為：長度20 m、高度4 m、寬度4 m;水槽水深2m;其余部分為空氣，消波區(qū)位于水槽末端16～20 m處;復(fù)合結(jié)構(gòu)的上部大尺度圓柱的半徑為a=0.4 m、圓柱吃水T=1 m、樁柱半徑b=0.1m，樁的長度L=1 m，樁處于上部圓柱靠入射波的左側(cè)。最終劃分后的網(wǎng)格如圖2所示。在ICEM中劃分好的結(jié)構(gòu)網(wǎng)格導(dǎo)入Fluent，頂部為壓力入口邊界，模擬水槽上部與外界大氣相連通。

3.2? ?參數(shù)設(shè)置

通過UDF的UDEFINE_CG_MOTION編寫動邊界的運動規(guī)律，以及DEFINE_SOURCE控制區(qū)域動量實現(xiàn)源項消波，并導(dǎo)入三維數(shù)值波浪水槽，水槽左壁面相當于造波板，通過動網(wǎng)格模型實現(xiàn)壁面的運動。采用VOF方法獲取自由水面，壓力速度項用Piso算法計算，使用動態(tài)層方法（Layering）配合網(wǎng)格重構(gòu)來實現(xiàn)網(wǎng)格的運動，時間間為0.01 s、模擬時長為20 s。

3.3? ?模型驗證

模擬規(guī)則波單獨沖擊下部樁柱的過程來確認模型的可靠性，將數(shù)值模擬受力與莫里森公式結(jié)果對比分析。初始入射波浪為線性規(guī)則波，入射波高設(shè)定為0.15m、周期為1.6 s。從圖3中可以看出，數(shù)值模擬得到的波浪力與莫里森公式結(jié)果大致吻合。所以可以初步認定FLUENT在數(shù)值模擬波浪與小尺度垂直立柱式結(jié)構(gòu)物相互作用方面具有相當?shù)木取?/p>

4? ? 數(shù)值模擬結(jié)果分析

4.1? 上部承臺尺寸對樁基波浪力的影響

當上部圓柱的半徑取a的不同倍數(shù)，設(shè)不同圓柱的半徑為R=na，且n分別為n=1、2、3、4，考察不同圓柱半徑情況下承臺產(chǎn)生的繞射場對下部樁基礎(chǔ)波浪荷載的影響。入射波浪特征為波幅A/a=0.75、波數(shù)ka=0.8的規(guī)則波。圖4是上部圓柱半徑R不同時，樁所受的入射波浪作用力、繞射波浪作用力以及合力幅值的比較。

從圖4可以看出，上部承臺結(jié)構(gòu)物半徑的增加導(dǎo)致流場的繞射作用越加明顯，入射和繞射波浪場影響下的樁基礎(chǔ)波浪載荷也逐漸增加，但波浪力增加的幅度隨著圓柱半徑增加逐漸減緩;同時，由于繞射力的相位與入射力有一定差異，使得合力幅值在大多數(shù)時間小于入射波浪力。

4.2? ?極端波浪對樁基波浪力的影響

畸形波是一種特殊的非線性極端波浪，無論是良好還是惡劣的海況均有出現(xiàn)。海洋畸形波具有能量高度集中和出現(xiàn)隨機性的特點，對于海上結(jié)構(gòu)物具有極大的威脅，因此有必要進行相關(guān)研究。

在Longuet-Higgins隨機波浪模型基礎(chǔ)上，劉贊強等[8]提出根據(jù)區(qū)間調(diào)整初始相位的方法，從而實現(xiàn)在特定時間、特定點位生成畸形波;林炅增[9]等改變不同規(guī)則波疊加后波面方程的初始相位，使得疊加波浪在預(yù)定聚焦點位置和聚焦時刻疊加形成大波。本文采用的波譜為JONWAP譜，有效波高為0.15 m、譜峰周期為1.6 s，按等分頻率法將波譜分為100份，波頻率區(qū)間為 [0.5～3.5]。本文模型所生成的畸形波波面時程曲線圖，見圖5。

圖6為畸形波入射時，承臺復(fù)合樁基礎(chǔ)受到的總力、入射波浪作用力和繞射波浪作用力的對比圖?？梢钥闯觯喝肷洳ɡ藞鲇绊懴碌臉痘A(chǔ)受力最大值為26N;繞射波浪場影響下的樁基礎(chǔ)受力最大值為12N，為入射波浪力的46.2%;入射和繞射波浪場共同影響下的樁基礎(chǔ)受力最大值為32 N，是入射波浪場影響下的1.23倍。由此可知，繞射波浪力的大小幾乎達到了正常規(guī)則波條件下的波浪力大小，說明極端海況下考慮繞射波浪力的必要性。雖然繞射產(chǎn)生的波浪場與入射波浪場的相位有一定差異，但繞射力對合力的貢獻為正，總力達到了正常規(guī)則波條件下波浪力的3倍，其中入射波浪力與繞射波浪力的相位差異與承臺下的樁柱位置和入射波方向有關(guān)。

5? ? ?結(jié)論

本文基于VOF法在Fluent中構(gòu)建波浪對承臺復(fù)合樁結(jié)構(gòu)作用的數(shù)值波浪環(huán)境，模擬了在規(guī)則波以及極端波浪畸形波入射，考慮繞射場影響下樁柱所受的波浪荷載情況。對于本文給定的算例，可得出以下結(jié)論：（1）承臺尺寸的不同改變了波浪場中繞射場的形態(tài)，繞射場作用在樁基礎(chǔ)上的波浪力隨著上部承臺尺寸的增大而增大;（2）在畸形波作用下，入射和繞射波浪場共同影響下的樁基礎(chǔ)受力最大值，是入射波浪場影響下的1.23倍，繞射波浪力大小幾乎達到了正常規(guī)則波條件下的波浪力大小，說明了畸形波影響下繞射作用的影響顯著。

參考文獻

[1]林毅峰， 陸忠民， 李彬. 海上風(fēng)電場風(fēng)機基礎(chǔ)結(jié)構(gòu).

[2] TS 145-2015. 港口與航道水文規(guī)范[S].中華人民共和國交通運輸部，2015.

[3]林毅峰， 陸忠民， 黃俊， et al. 海上風(fēng)電機組高承臺群樁基礎(chǔ)設(shè)計特點及關(guān)鍵力學(xué)問題[J]. 海洋技術(shù)學(xué)報， 2016， 35（5）：29-36.

[4]MacCamy R C， Fuchs R A 1954 US Army Coastal Engineering Research?Center， 69.

[5]Linton C M， Evans D V 1990 J.Fluid Mech. （215）：549-569.

[6]Geng B L， Teng B， Ning D Z and Gou Y 2009 Journal of Dalian Maritime?University. 35（3）：5-8.

[7]Jiang S C， Teng B， and Ning D Z 2010 Ocean Engineering. 28（3）：68-75.

[8]劉贊強， 張寧川， 俞聿修， et al. 改進的相位調(diào)制法模擬畸形波：I—理論模型與驗證[J]. 水動力學(xué)研究與進展A輯， 2010（03）：111-118.

[9]林炅增， 胡金鵬. 畸形波作用下海上風(fēng)機單樁基礎(chǔ)所受波浪力研究[J].?廣東造船， 2019， 38（01）：23-26.