海上風(fēng)電單樁基礎(chǔ)軸向承載特性及軟弱土層影響分析

李超杰，朱洪澤，蘇浩然，蘇 凱，3，4

（1.浙江華東工程咨詢有限公司，浙江 杭州 310014；2.武漢大學(xué)水資源與水電工程科學(xué)國家重點(diǎn)實驗室，湖北 武漢 430072；3.武漢大學(xué)水工巖石力學(xué)教育部重點(diǎn)實驗室，湖北 武漢 430072；4.武漢大學(xué)海綿城市建設(shè)水系統(tǒng)湖北省重點(diǎn)實驗室，湖北 武漢 430072）

0 引 言

海上風(fēng)電基礎(chǔ)的主要作用是固定風(fēng)機(jī)整體結(jié)構(gòu)并維持其穩(wěn)定性，根據(jù)基礎(chǔ)結(jié)構(gòu)是否與海床直接接觸可分為固定式基礎(chǔ)與漂浮式基礎(chǔ)。當(dāng)前我國近海風(fēng)電基本采用固定式基礎(chǔ)，其中單樁基礎(chǔ)較其他固定式基礎(chǔ)型式有結(jié)構(gòu)簡單、施工簡便快捷、經(jīng)濟(jì)性好等優(yōu)點(diǎn)［1］，是目前應(yīng)用最為廣泛的基礎(chǔ)之一［2，3］。近年來，隨著我國風(fēng)電事業(yè)的快速發(fā)展，超長大直徑單樁基礎(chǔ)在海上風(fēng)電中應(yīng)用越來越多［4，5］。

對于單樁基礎(chǔ)而言，其水平承載力關(guān)系到風(fēng)電結(jié)構(gòu)運(yùn)行期位移變形情況，而軸向承載力則直接關(guān)系到施工過程的順利與否。目前，眾多學(xué)者對單樁基礎(chǔ)軸向承載力展開了相關(guān)研究。梁超［6］對比了多個樁基軸向承載力計算理論公式及設(shè)計規(guī)范，并指出理論計算難以精確評估；董淑海［7］對水下大直徑單樁基礎(chǔ)進(jìn)行了原位試驗以分析結(jié)構(gòu)受力情況及樁-土作用特性；Lehane［8］對松砂條件下單樁軸向承載能力進(jìn)行了試驗研究，并分析了樁基礎(chǔ)內(nèi)壁與土體應(yīng)力水平、內(nèi)摩擦角及剪脹角間關(guān)系；Paik［9］對比了開口樁與閉口樁軸向承載力的差異，研究并推導(dǎo)了土塞效應(yīng)的樁端承載力公式；劉冰雪［10］通過數(shù)值分析討論了樁周土體參數(shù)、樁間距等因素對軸向承載力的影響；徐國賓［11］基于數(shù)值模擬研究并確定了影響基礎(chǔ)軸向承載力的主要因素；張明遠(yuǎn)［12］基于FLAC3D 對超長大直徑鋼管樁進(jìn)行數(shù)值分析，模擬得到基礎(chǔ)軸向承載力與實測結(jié)果一致；劉文白［13］、黃周泉［14］對樁桶符合型海上風(fēng)機(jī)基礎(chǔ)的軸向承載力進(jìn)行了研究，并對該結(jié)構(gòu)提出了設(shè)計優(yōu)化方法。

目前大多數(shù)研究針對單一或特定海域土體基礎(chǔ)，并未討論地基中包含軟弱土層及其變化對于樁基軸向承載性的影響，然而我國部分海域砂土上存在不同厚度軟弱土層，例如上覆淤泥質(zhì)粉砂土的等［15］，勢必將對樁基軸向承載力產(chǎn)生影響，因此有必要進(jìn)行深入研究。本文借助大型通用有限元軟件平臺ABAQUS，首先討論了軸向承載特性判斷標(biāo)準(zhǔn)的差異，并進(jìn)行了土體參數(shù)敏感性分析，隨后討論了軟弱土層位置及其厚度對軸向承載力的影響，以期為后續(xù)結(jié)構(gòu)優(yōu)化設(shè)計提供一定參考。

1 軸向承載特性判據(jù)及模型

1.1 軸向承載特性判據(jù)

對于摩擦-端承樁而言，單樁極限軸向荷載指當(dāng)樁側(cè)摩阻力和樁端阻力導(dǎo)致樁周及樁端土體都達(dá)到塑性狀態(tài)，此時較小的軸向荷載增量將引起較大的沉降量。學(xué)術(shù)界與工程界以及不同行業(yè)、國家地區(qū)對于確定摩擦-端承樁極限荷載目前并沒有一個較為統(tǒng)一的標(biāo)準(zhǔn)［10］，主要通過荷載-位移曲線（Q-S曲線）類別進(jìn)行判別［16］：①對于陡變型，通過確定平緩段與陡降段的拐點(diǎn)來確定樁基礎(chǔ)軸向極限承載力。此類方法受Q-S曲線比例尺影響較大，須選用合理的比例尺，以求得較為合理的軸向極限荷載；②對于緩變型曲線，無明確拐點(diǎn)，則根據(jù)樁頂沉降量來確定極限荷載，根據(jù)承載能力極限狀態(tài)和正常使用極限狀態(tài)不同而確定。這種方法具有明確且簡單的特征，一般情況下采用與樁徑相關(guān)的標(biāo)準(zhǔn)來進(jìn)行判斷。Vesic［17］等人建議對于打入樁和灌注樁分別采用樁頂下沉量等于10%、25%樁徑時的荷載為極限荷載。表1給出了兩種軸向極限承載力控制標(biāo)準(zhǔn)。

表1 樁基礎(chǔ)軸向極限承載力控制標(biāo)準(zhǔn)Tab.1 Different control standards of axial ultimate bearing capacity

1.2 模型及驗證

為比較不同變形控制標(biāo)準(zhǔn)所確定軸向承載力大小，依托江蘇某800MW 海上風(fēng)電場項目為工程背景，建立相關(guān)三維有限元模型。參考工程統(tǒng)一采用單樁基礎(chǔ)式風(fēng)機(jī)，樁徑D=6.0 m，樁長L=74 m，其中泥面以下埋深MCD=52 m，樁基礎(chǔ)簡化為等壁厚，其值為d=70 mm。根據(jù)工程實測資料設(shè)鋼管樁的彈性模量E=210 GPa，泊松比ν=0.30，密度ρ=7.80 g∕cm3。參考風(fēng)場海域地貌屬于南黃海濱海相沉積地貌單元，海底地形起伏明顯，根據(jù)工程勘測資料，各土層厚度及材料參數(shù)如表2所示，其中土層⑤為持力層，土層②為流塑狀軟土，工程性能差，屬于軟弱土層。

表2 地基各層土體參數(shù)表Tab.2 Soil parameters of each layer of foundation

利用ABAQUS建立相關(guān)有限元模型，并使用Mohr-Coulomb彈塑性屈服準(zhǔn)則模擬地層土體的力學(xué)性質(zhì)，樁基礎(chǔ)則假定為線彈性材料；模型水平范圍取為10 倍樁徑，縱向范圍為2 倍基礎(chǔ)埋深；整個模型均采用C3D8R 三維8 節(jié)點(diǎn)六面體單元進(jìn)行劃分；模型邊界方面，地層土體底部約束其全部位移，同時約束地基土側(cè)面的法向位移，地基土體頂部為自由面；樁-土相互作用關(guān)系采用面-面接觸模擬，其接觸屬性法向為“硬接觸”，切向類型為“罰”摩擦函數(shù)，摩擦系數(shù)取0.3。同時，建立以土層⑤-中砂材料參數(shù)為基準(zhǔn)的均質(zhì)地基模型，并將其作為基礎(chǔ)方案，對應(yīng)有限元模型如圖1所示。

圖1 有限元模型（基礎(chǔ)方案）Fig.1 Finite element model （benchmark scheme）

為對比不同控制標(biāo)準(zhǔn)下單樁基礎(chǔ)軸向承載力差異，在樁頂處梯度施加0～150 MN 豎直荷載，并提取樁頂軸向荷載-沉降曲線，如圖2所示。由圖2可知，該曲線具有較為明顯的拐點(diǎn)和陡降段，其破壞特征點(diǎn)比較顯著，標(biāo)準(zhǔn)1所對應(yīng)的曲線轉(zhuǎn)折點(diǎn)處軸向荷載的大小為60.4 MN；若采用標(biāo)準(zhǔn)2，樁頂沉降位移達(dá)到0.1D時，單樁頂部所對應(yīng)的軸向荷載為141.3 MN。

圖2 基礎(chǔ)方案中軸向荷載-沉降曲線Fig.2 Axial load-displacement curve in benchmark scheme

圖3為樁基礎(chǔ)沉降剛好達(dá)到變形控制標(biāo)準(zhǔn)時其樁身軸力分布曲線。由圖3可知泥面以上范圍內(nèi)，由于單樁存在自重，樁身軸力沿樁頂向下緩慢增長；泥面以下由于樁內(nèi)外摩阻力的存在，樁身軸力隨埋深增加而顯著減少，且深度越大，其衰減速率越快；同時樁端處仍存在著一定軸力，說明在上述標(biāo)準(zhǔn)下樁內(nèi)外摩阻力均不能完全抵消樁頂荷載，仍有一部分荷載由樁端阻力承擔(dān)。摩擦型樁大部分軸向荷載通過樁側(cè)摩阻力以剪應(yīng)力的形式傳遞到樁周土體中，樁側(cè)摩阻力大小決定了樁身軸力的傳遞速度。該典型樁樁周上部側(cè)摩阻力較小，樁身軸力傳遞較慢。

圖3 不同控制標(biāo)準(zhǔn)下的單樁軸力分布Fig.3 Axial force distribution under different control standards

如圖4展示了不同控制標(biāo)準(zhǔn)下的土體塑性區(qū)分布。當(dāng)樁頂沉降達(dá)到變形控制標(biāo)準(zhǔn)1 時，地基土體幾乎沒有產(chǎn)生等效塑性應(yīng)變；而樁頂沉降位移達(dá)到標(biāo)準(zhǔn)2時，樁基底端及側(cè)面出現(xiàn)較大范圍塑性區(qū)，且由于單樁剛度遠(yuǎn)大于地基土體，單樁正下方處的塑性區(qū)出現(xiàn)明顯的延伸。

圖4 不同控制標(biāo)準(zhǔn)下地基土體等效塑性應(yīng)變Fig.4 Equivalent plastic strain distribution of soil under different control standards

由上可知變形控制標(biāo)準(zhǔn)1、2所對應(yīng)單樁軸向承載力相差較大。當(dāng)采用與樁徑相關(guān)的標(biāo)準(zhǔn)，即以打入樁樁頂沉降10%樁徑作為其軸向承載力控制標(biāo)準(zhǔn)時，樁端土體出現(xiàn)較大范圍塑性應(yīng)變，此種變形控制標(biāo)準(zhǔn)難以適用于海上風(fēng)機(jī)大直徑鋼管樁在土體參數(shù)較差的分層地基中的情況。當(dāng)樁頂沉降達(dá)到變形控制標(biāo)準(zhǔn)1 后，軸向荷載-沉降曲線進(jìn)入陡降段，這說明以Q-S曲線之中轉(zhuǎn)折點(diǎn)處所對應(yīng)的軸向荷載作為參考標(biāo)準(zhǔn)可以更好地反映摩擦-端承樁單樁極限軸向荷載的物理意義，且較為直觀。綜上所述，本文以荷載-沉降（Q-S）曲線中出現(xiàn)較為明顯陡降段時的轉(zhuǎn)折點(diǎn)所對應(yīng)的軸向荷載的大小來確定基礎(chǔ)的極限承載力。

2 土體參數(shù)影響分析

2.1 地基土體彈性模量對樁基軸向承載力影響分析

采用圖1所示均質(zhì)土體模型，參照表1中不同土層土體彈性模量范圍，在保證其他條件不變的前提下梯度設(shè)置10～35 MPa多個土體彈性模量方案進(jìn)行加載計算。

相同荷載作用下，土體彈模越大對應(yīng)單樁樁頂沉降量越小。土體彈性模量從10 MPa 增長到35 MPa 的過程中，在60 MN 同級荷載作用下，單樁樁頂沉降從大約150 mm 下降到60 mm，降幅達(dá)到60%；軸向承載力從58.53 MN 增長到60.79 MN，僅增長3.9%。由此說明地基土體彈性模量的增大對提升單樁軸向承載力作用較小，這與劉冰雪［10］等人研究結(jié)果一致。

2.2 地基土體摩擦角對樁基軸向承載力影響分析

圖5 地基土體彈性模量的影響Fig.5 Influence of elastic modulus of soil

與前文一致，參照實際工程中地層土體摩擦角范圍，設(shè)置25°～35°多個地基摩擦角方案進(jìn)行加載計算。

與彈性模量影響相似，樁頂沉降位移與土體摩擦角成反比，但總體而言變化幅度較小。同時，土體摩擦角從25°增長到35°的過程中，軸向承載力增長了56.1%，說明地基土體摩擦角的增大對提升單樁軸向承載力作用顯著，且基本呈線性增長。

3 軟弱土層的影響

在實際工程中，地基土層存在性質(zhì)較差的軟弱土層，同時上文也表明土體參數(shù)對基礎(chǔ)軸向承載力存在較大影響，因此為研究軟弱土層相對位置、厚度因素對單樁基礎(chǔ)軸向承載力影響，選取表1中強(qiáng)度較低的②-粉砂夾粉土作為軟弱土層，模型中其余土體部分仍采用⑤-中砂。

3.1 軟弱土層位置影響

為研究分層地基中軟弱土層相對位置對單樁軸向承載力的影響，取軟弱土層厚度為2D，并以軟弱土層頂面為基準(zhǔn)，設(shè)置軟弱土層埋深由0 m 至泥面以下64 m，共設(shè)置17 種分層地基方案進(jìn)行加載計算，并與前文基礎(chǔ)方案進(jìn)行對比。圖7分別為軟弱土層頂部高程距離泥面0、64 m時模型示意圖。

圖7 軟弱土層埋深示意圖Fig.7 Schematic diagram of buried depth of soft soil layer

圖8（a）～（c）展示了部分方案Q-S曲線。相較于無軟弱土層時，隨著軟土從海床泥面處逐漸向下移動，單樁軸向承載力逐漸下降，當(dāng)軟弱土層頂部深度為44 m時單樁軸向承載力降幅最大，此時單樁軸向承載力由無軟弱土層時60.40 MN 下降到50.24 MN，降幅為16.8%。隨著軟弱土層埋深進(jìn)一步增加，單樁軸向承載力有較大提升，逐漸逼近無軟弱土層時單樁軸向承載力且最終趨于穩(wěn)定。結(jié)合樁基礎(chǔ)埋深信息，上述結(jié)果說明樁端周圍土體對單樁軸向承載力影響較大，而上層土體對單樁軸向承載力影響相對較小。

圖8 軟弱土層位置的影響Fig.8 Influence of soft soil layer position

3.2 軟弱土層厚度影響

為進(jìn)一步研究不同軟土厚度對軸向承載力的影響大小，結(jié)合上一小節(jié)結(jié)論，將軟弱土層底部固定在樁端以下12 m 處，并設(shè)置不同厚度軟土地基方案進(jìn)行加載計算。圖9為軟弱土層厚度為4、64 m時的三維有限元模型示意圖。

圖9 軟弱土層厚度示意圖Fig.9 Schematic diagram of thickness of soft soil layer

如圖10所示為不同軟弱土層厚度條件下樁基礎(chǔ)軸向承載力變化情況。由圖可知軟土層厚度與軸向承載力成反比。同時，對于初期，增加軟弱土層厚度直到16 m時，單樁軸向承載力迅速降低；而后當(dāng)軟弱土層厚度為16～32 m 時，降速基本恒定；再增加軟弱土層厚度，單樁軸向承載力持續(xù)降低，但降速也有所下降，直至上部土體全部為軟弱土層時，軸向承載力為31.919 MN，相較于無軟弱土層時，總降幅達(dá)到47%。

圖6 土體摩擦角的影響Fig.6 Influence of friction angle of soil

圖10 不同軟弱土層厚度下單樁軸向承載力變化曲線Fig.10 Variation of axial ultimate bearing capacity under variable thicknesses of soft soil layer

在實際工程應(yīng)用中，通過增大樁端附近樁土之間摩阻力來提高單樁整體軸向承載力是一種性價比較高的工程措施，然而對于海上風(fēng)電單樁基礎(chǔ)而言，改變樁端處土層泥質(zhì)參數(shù)是較為困難的，因此可以通過對靠近樁端處樁身表面進(jìn)行磨砂處理等工藝提高單樁軸向承載力。

4 結(jié) 論

本文依托江蘇某實際工程，首先對單樁基礎(chǔ)軸向承載極限控制標(biāo)準(zhǔn)進(jìn)行了對比，并討論了不同土體參數(shù)對樁身沉降及軸向承載力的影響，最后考慮土體地基中軟弱土層位置及厚度所產(chǎn)生的影響，最終得到如下結(jié)論。

（1）相較于通過樁基礎(chǔ)沉降量為確定基礎(chǔ)軸向承載力的控制標(biāo)準(zhǔn)，根據(jù)荷載-位移曲線拐點(diǎn)確定均質(zhì)土層下單樁極限承載力為60.4 MN，僅為前者控制標(biāo)準(zhǔn)的42.7%，二者差距較大。同時，達(dá)到控制標(biāo)準(zhǔn)時，前者樁端附近土體已出現(xiàn)較大范圍塑性區(qū)，控制標(biāo)準(zhǔn)難以適用于土體參數(shù)較差情況，因此建議采用根據(jù)荷載-位移曲線拐點(diǎn)作為軸向極限承載力計算控制標(biāo)準(zhǔn)。

（2）單樁基礎(chǔ)屬于端承摩擦型樁，所受軸向荷載大多通過樁側(cè)摩阻力以剪應(yīng)力的形式傳遞到樁周土體中，因此提高土體摩擦角可顯著提高基礎(chǔ)軸向承載力；提高土體彈?？纱蠓档蜆俄敵两滴灰疲珜μ岣咻S向承載力效果并不顯著。

（3）當(dāng)土體中存在軟弱土層時，軟弱土層越接近樁端，軸向承載力降幅越大，假定軟弱土層厚度為2倍樁徑時，最大降幅達(dá)到16.8%；當(dāng)軟弱土層位于樁端時，其厚度的增加將持續(xù)降低樁基軸向承載力，并且呈現(xiàn)三段式分布，最終軸向極限承載力降幅為47%。因此建議通過表面磨砂等工藝對靠近樁端處樁身表面進(jìn)行處理以提高單樁軸向承載力。