桶形基礎二次加載的大變形有限元模擬

陳霄 鄭敬賓 王棟

摘 要：桶形基礎用于支撐海上風電升壓站平臺時，平臺與電氣設備分兩階段安裝，中間存在數(shù)周到數(shù)月的間隔期。為探究桶形基礎在二次加載階段的峰值貫入阻力，發(fā)展了基于網(wǎng)格周期更新的有效應力大變形有限元方法，避免基礎貫入引起周圍土體單元扭曲，再現(xiàn)了黏土中基礎“貫入固結(jié)二次加載”的全過程。對比已有的離心模型試驗，驗證了發(fā)展的大變形數(shù)值方法的可靠性。通過比較初始貫入阻力與固結(jié)后的貫入阻力，發(fā)現(xiàn)桶基與平臺的安裝間隔期能夠有效增加桶基承載力，固結(jié)后的貫入阻力峰值僅需要很小的基礎位移即可達到。開展變動參數(shù)分析，考察影響峰值貫入阻力比的主要因素，提出考慮基礎尺寸、土體滲透性與固結(jié)時間的歸一化時間參數(shù)，并給出峰值貫入阻力比隨歸一化時間增長的趨勢線。

關鍵詞：桶形基礎;有限元分析;固結(jié);黏土;貫入阻力

中圖分類號：TU470 ? 文獻標志碼：A ? 文章編號：2096-6717（2022）01-0068-07

收稿日期：2021-03-25

基金項目：國家自然科學基金（41772294、U1806230、42025702、51809247）

作者簡介：陳霄（1996- ），男，主要從事海洋巖土工程研究，E-mail：chenxiao5185@stu.ouc.edu.cn。

王棟（通信作者），男，教授，博士生導師，E-mail：dongwang@ouc.edu.cn。

Abstract：

When caissons are used to support offshore wind substation platforms， the platform and electrical equipment are installed in two phases， with a waiting period of several weeks to months between the two events. To evaluate the peak resistance during secondary loading phase，this 'penetration-consolidation-penetration' process in clay is reproduced by the coupled effective stress large deformation finite element analysis to avoid the severe distortion of elements during installation. The numerical model is first verified by simulating a centrifuge test.A comparison of the initial penetration resistance with the penetration resistance after consolidation shows that the waiting period can cause a dramatically increase in soil strength. And the depth of peak penetration resistance during secondary loading phase is almost equal to the initial installation， implying that only a small displacement is required to achieve the compressive ultimate limit state of the caisson after the waiting period. Parametric study is subsequently carried out to investigate the effect of consolidation. A normalized time parameter considering foundation size， soil permeability and consolidation duration is proposed， and a trend line describing the increase of peak caisson resistance with normalized time is provided.

Keywords：caissons; finite element analysis; consolidation; clay; penetration resistance

近年來，負壓桶形基礎（以下簡稱為桶基）作為單樁基礎和導管架樁基的替代形式，已用于支撐多個海上風電場的風機和升壓站，如英國的Hornsea 3項目[1]。與風機的一次性吊裝不同，桶基支撐的升壓站平臺分兩個階段安裝：第1階段，桶基與平臺貫入土中預定深度，施工船返回港口，桶基靜置數(shù)周至數(shù)月;第2階段，各種大型電氣模塊運輸至場地，并吊放到平臺上，桶基實質(zhì)上承受二次加載。第1階段的桶基貫入在黏土中引起超靜孔壓，之后孔壓部分或完全消散[2-3]，導致第2階段開始時桶基周圍黏土的不排水強度高于第1階段結(jié)束時。

類似于一般的淺基礎承載力問題，常規(guī)有限元方法被廣泛用于分析不排水條件下桶基的豎向承載力及復合承載力[4-6]，但這些研究均假定桶基預埋就位，不考慮桶基安裝過程對承載力的影響。這是因為桶基長距離貫入導致側(cè)壁端部周圍土體單元的嚴重扭曲，常規(guī)有限元方法無法捕捉這一大變形過程。

研究者借助離心機試驗[7-11]或大變形有限元模擬[12-13]再現(xiàn)了桶基的貫入，探究了土層分布、安裝方式、基礎尺寸等一系列影響桶基貫入阻力的因素，提出了確定第1階段貫入阻力的方法。少量的數(shù)值模擬還進一步給出了桶基安裝后的孔壓消散過程[3]，但這些結(jié)果僅評估了地基土的固結(jié)程度，不能預測升壓站平臺第2階段的抗壓承載力。Jeanjean[14]在墨西哥灣海域進行了現(xiàn)場試驗，并總結(jié)了歷史上該區(qū)域的施工記錄，得到桶基安裝就位時的貫入阻力以及固結(jié)一段時間后的抗拔力，發(fā)現(xiàn)桶基在安裝后存在顯著的固結(jié)效應，最終抗拔力最高可達初始貫入阻力的1.6倍，但其試驗針對的是抗拔力，而不是抗壓承載力。

鑒于與桶基二次加載有關的研究較為有限，筆者發(fā)展了有效應力形式的大變形有限元方法，模擬黏性土中桶基的“貫入固結(jié)二次加載”完整過程，探索桶基初始貫入和部分排水固結(jié)后的超靜孔壓分布，進而考察固結(jié)效應對二次加載的影響，合理評估二次加載的豎向承載力。

1 研究方法

1.1 大變形有限元方法

桶形基礎的安裝是典型的巖土工程大變形問題，因此使用一種名為RITSS的大變形有限元方法（Remeshing and interpolation technique by small strain）。這種方法最初由Hu等[15]提出。Wang等[16-17]將其與有限元商業(yè)軟件Abaqus結(jié)合并擴展至有效應力分析。RITSS方法中，整個計算過程被劃分為一系列小變形計算步，通過網(wǎng)格重劃分及場變量映射，避免了單元的嚴重扭曲，具體細節(jié)可參見文獻[17-18]。

采用基于Abaqus的RITSS方法，將桶基的“貫入固結(jié)二次加載”過程簡化為二維軸對稱問題進行分析，建立如圖1所示的數(shù)值模型。桶基的長度、直徑和厚度分別為L、D和t，支撐電氣平臺的桶基長徑比L/D一般為1左右。為避免邊界效應，土體的徑向與豎向邊界分別取5D與8D。約束底邊徑向邊界的水平與豎向位移，豎向邊界的水平位移。桶基周圍附近土體的網(wǎng)格尺寸為1.5t。為提高計算效率，并避免貫入初期土體表面單元的劇烈變形，將桶基預埋一定深度。后文將證明該設置對后續(xù)貫入阻力的影響可以忽略。采用二階軸對稱縮減積分單元劃分土體，以保證映射精度[17]。

采用罰函數(shù)法定義“桶土”界面的摩擦力，在有效應力分析中，根據(jù)庫倫摩擦定律，摩擦力取決于法向有效應力及界面的庫倫摩擦系數(shù)。安裝造成桶基周圍土體的應變軟化，設計時一般假定周圍土體完全擾動。Koh等[3]采用大變形方法模擬桶基貫入時取摩擦系數(shù)為0.3tan ′/St，其中′為有效內(nèi)摩擦角，St為靈敏度，得到了與離心機試驗結(jié)果相近的貫入阻力。筆者的數(shù)值分析也采用相同的摩擦系數(shù)。

以往離心機試驗[7-11]比較了吸力式安裝（即通過氣泵從桶蓋與土體之間抽水以施加吸力）與純壓入式安裝（即采用位移控制使桶基貫入）兩種桶基貫入過程，發(fā)現(xiàn)安裝方式對黏土中桶基貫入阻力的影響可以忽略不計。因此，通過位移控制模式實現(xiàn)桶基的初始貫入和二次加載。在固結(jié)階段，桶基上的豎向力保持為定值，豎向力大小等于當前深度的初始貫入阻力。在整個“貫入固結(jié)二次加載”過程，僅允許土體表面邊界排水，其余各邊均為不排水邊界。

1.2 土體本構(gòu)模型

采用經(jīng)典的修正劍橋模型描述黏性土的有效應力與應變之間的關系，該模型已被成功用于分析靜力觸探儀、樁靴基礎的大變形貫入[18-19]。考慮兩種廣泛用于離心機試驗中的黏性土，即UWA高嶺土和Malaysian高嶺土，相關計算參數(shù)總結(jié)于表1。側(cè)壓力系數(shù)K0=1

2 數(shù)值模擬驗證

目前，還沒有公開發(fā)表的桶形基礎二次加載模型試驗或現(xiàn)場監(jiān)測數(shù)據(jù)。通過對比桶基首次安裝時的貫入阻力，驗證數(shù)值模擬的可靠性。Westgate等[9]完成了6個桶基在UWA高嶺土中貫入的離心機試驗。使用RITSS方法模擬其中名為Type1的試驗。桶基的原型尺寸為D=11.3 m，L=6 m，t=0.05 m。試驗前和試驗后分別進行兩次T-bar試驗，得到土體的平均抗剪強度su=10+2.8z kPa。利用兩次循環(huán)插拔T-bar試驗獲得土體的靈敏度St為2.5和3。大變形有限元分析中使用St=3，確定桶土摩擦系數(shù)。試驗中測得γ′=5.9 kN/m3，根據(jù)式（1）推導得到修正劍橋模型中的su=1.44z kPa。為便于比較，將貫入阻力F歸一化為F/suA，A為桶形基礎豎向投影面積，A=πD2/4。圖2比較了離心模型試驗[9]、有效應力形式的RITSS法、總應力形式的耦合歐拉拉格朗日有限元法（Coupled Eulerian-Lagrangian，CEL）結(jié)果[13]、DNV規(guī)范的預測[23]。圖中的縱坐標為歸一化貫入深度d/D，即貫入深度與桶基直徑之比。

DNV規(guī)范與CEL的結(jié)果均基于總應力分析，假定安裝時地基土完全不排水。其中DNV規(guī)范建議周圍土體的強度折減系數(shù)α=1/St=0.33，表示桶的最大側(cè)壁摩阻力為αsu。文獻[13]中CEL分析使用的系數(shù)為α=0.2，故圖2也給出了α=0.2時的DNV方法預測結(jié)果。由圖2可見，有效應力形式的RITSS方法和總應力形式的CEL方法給出的結(jié)果都與試驗數(shù)據(jù)接近，且均在DNV公式預測的α=0.2和α=0.33兩條曲線范圍內(nèi)。采用RITSS方法時，預埋使得計算曲線從0.2D開始，因此，無法比較0～0.2D范圍內(nèi)的結(jié)果。然而，當桶基的貫入深度超過0.25D后，貫入阻力曲線與試驗非常接近，表明一旦土體破壞模式被完全激發(fā)，預埋對貫入阻力不會產(chǎn)生明顯影響，這與已有結(jié)論一致[18-19]。圖2的比較也表明大變形有限元分析中的摩擦系數(shù)取值合理，可以在很大程度上考慮安裝擾動對貫入阻力的影響。

3 結(jié)果分析與討論

實際工程中鋼制桶基的壁厚一般小于0.05 m。試算發(fā)現(xiàn)，當桶基側(cè)壁厚度較小時，有效應力形式的大變形分析不易收斂，后續(xù)分析中統(tǒng)一取壁厚為0.19 m。Xiao等[13]比較了側(cè)壁較厚（t=0.2～0.6 m）與較薄（t=0.05 m）的桶基端部承載力系數(shù)，發(fā)現(xiàn)其范圍在10～12之間，表明一定范圍內(nèi)壁厚的增加沒有改變桶基貫入時端部以下土體的破壞模式，而側(cè)壁附近的土體破壞模式不依賴壁厚。因此，增加壁厚不會顯著影響大變形結(jié)果。

表2總結(jié)了變參數(shù)分析中開展的算例。為探究固結(jié)系數(shù)cv的影響，算例A～算例D分別模擬了桶基在UWA高嶺土和Malaysian高嶺土中的“貫入固結(jié)二次加載”過程。這兩種土的本構(gòu)模型參數(shù)接近，但滲透系數(shù)相差很大，相同應力水平下Malaysian高嶺土的cv約為UWA高嶺土的10倍。算例A～算例D還將用于考察首次安裝深度zp對二次加載的影響。算例E和算例F將考察桶基直徑D對二次加載的影響。表2中的固結(jié)時間tc在7～360 d之間變化，其下限是實際工程中最短的施工間隔，而上限人為取得較大，以便展示充分固結(jié)造成的二次加載效應。

修正劍橋模型并不能考慮土體結(jié)構(gòu)性的影響，分析中摩擦系數(shù)始終保持為0.3tan ′/St，這也意味著固結(jié)前后的土體靈敏度未改變，即不考慮擾動的土體結(jié)構(gòu)在固結(jié)過程中部分恢復，所以，計算結(jié)果偏于安全。以下分析中統(tǒng)一取典型值St=2.5。

3.1 固結(jié)時間對二次加載的影響

算例A、C（D=9.5 m，首次貫入深度zp=0.9D），固結(jié)7、30、60、120、270 d后，二次加載貫入阻力隨深度的變化如圖3所示。為了直觀顯示固結(jié)時間的影響，圖中包含了初始連續(xù)貫入阻力曲線作為參考。參考Wang等[19]對樁靴基礎貫入阻力的處理方法，利用峰值貫入阻力比Ncp/Ncr表征貫入阻力的改變，如圖3（a）所示。Ncp為二次加載時歸一化的峰值貫入阻力，Ncr為Ncp出現(xiàn)深度處歸一化的初始貫入阻力。

從圖3可以看出：

1）桶基的二次加載阻力在很短位移內(nèi)即可達到峰值，Ncp出現(xiàn)時的深度僅略大于初始貫入結(jié)束時的深度zp。

2）UWA高嶺土中即使固結(jié)時間僅有7 d，峰值貫入阻力也會增加約20%，見圖3（a）。如果固結(jié)時間為7～120 d，二次加載的Ncp隨固結(jié)時間增大，說明固結(jié)導致的桶基豎向承載力顯著提高。但120 d后的Ncp基本不再增加。

3）二次加載的貫入阻力曲線達到峰值后迅速降低，并向初始連續(xù)貫入阻力曲線靠近。這表明，盡管施工間隔期間桶基周圍土體固結(jié)，強度增加，但強度增加土體的范圍有限。這一特征與樁靴基礎的二次貫入類似[19]。

4）Malaysian高嶺土的結(jié)果表現(xiàn)出與UWA高嶺土類似的趨勢，如圖3（b）所示，Ncp/Ncr最大約為1.35。然而，固結(jié)60 d后Ncp就不再有顯著增長。產(chǎn)生這種差異的原因是兩種土的固結(jié)系數(shù)cv不同。Malaysian高嶺土的固結(jié)系數(shù)遠高于UWA高嶺土，初始貫入造成的累積孔壓在固結(jié)階段消散得更快，桶基周圍土體更快達到強度增長極限。

圖4展示了固結(jié)過程中桶基附近3個位置的土體超靜孔壓消散過程。3個位置分別為桶壁端部中心、端部上方0.02zp處與桶壁緊鄰的內(nèi)側(cè)與外側(cè)。為避免固結(jié)初期孔壓重分布對消散曲線的影響，利用Sully等[24]提出的方法確定每一個位置的名義初始孔壓Δumax，得到歸一化孔壓Δu/Δumax隨固結(jié)時間的變化曲線。圖4表明，3個位置處的超靜孔壓消散規(guī)律幾乎一致。對于UWA高嶺土，孔壓在最初的120 d迅速消散，120 d后的消散速度逐漸減緩，這與峰值貫入阻力的增長規(guī)律一致（圖3）。同樣，Malaysian高嶺土中超靜孔壓在前60 d急劇下降，且消散速度比UWA高嶺土更快。

由圖6可以看出，zp/D在0.7～1.5范圍內(nèi)時，兩種高嶺土參數(shù)計算得到的峰值貫入阻力比隨歸一化時間的變化趨勢一致。曲線整體上可以分成兩段：當T小于0.015時，Ncp/Ncr隨T大致呈線性增長關系;當T大于等于0.015時，Ncp/Ncr趨近于一個穩(wěn)定值。為了建立Ncp/Ncr與T的關系，在圖6中添加了趨勢線。地基土充分固結(jié)完成時，算例A的Ncp/Ncr達到1.36，算例F的Ncp/Ncr約為1.26，其他算例結(jié)果介于二者之間。也就是說，完全固結(jié)時二次加載得到的峰值貫入阻力最少增長約26%。

實際條件下，峰值貫入阻力的增長可能超出圖6的有限元結(jié)果，這是因為在固結(jié)過程中孔壓消散以及土體結(jié)構(gòu)性的恢復都可以造成后續(xù)承載力的上升，且二者的貢獻程度相當[14]。受限于本構(gòu)模型，本文的有限元分析還不能反映土體結(jié)構(gòu)性的恢復，所得峰值貫入阻力增長程度偏于保守。

4 結(jié)論

采用有效應力形式的大變形有限元方法模擬海上風電升壓站平臺的安裝過程，再現(xiàn)了桶基“貫入固結(jié)二次加載”的全過程。數(shù)值模擬中考慮了兩種力學性質(zhì)接近，但固結(jié)系數(shù)相差10倍的高嶺土，考察了固結(jié)時間、首次安裝深度、桶基直徑等關鍵變量對二次加載時貫入阻力的影響。主要結(jié)論如下：

1）當固結(jié)時間相同時，高固結(jié)系數(shù)土中二次加載的貫入阻力增長更迅速，但固結(jié)系數(shù)對峰值貫入阻力比極限值的影響不大，完全固結(jié)時的峰值貫入阻力比在1.26～1.36之間。

2）式（4）定義的歸一化固結(jié)時間較好地反映了峰值貫入阻力比的變化趨勢。采用歸一化固結(jié)時間表征固結(jié)程度時，峰值貫入阻力比隨歸一化固結(jié)時間約為線性增加，直至達到極限值。

3）實際施工中，應充分利用二次加載前的施工間隔時間，發(fā)揮地基土固結(jié)后的強度，從而提高二次加載時的極限承載力，達到優(yōu)化基礎尺寸，降低工程造價的目的。參考文獻：

[1] Orsted Hornsea Project Three （UK）. Hornsea Project Three Offshore Wind Farm. Examining authoritys report of findings and conclusions and recommendation to the Secretary of State for Business， Energy and Industrial [EB/OL]. https：//infrastructure.planninginspectorate.gov.uk/projects/eastern/hornsea-project-three-offshore-wind-farm/

[2] ANDERSEN K H， JOSTAD H P. Shear strength along outside wall of suction anchors in clay after installation [C]// Proceedings of the 12th International Offshore and Polar Engineering Conference. Kitakyushu： International Society of Offshore and Polar Engineers， 2002： 785-794.

[3] ?KOH K X， WANG D， HOSSAIN M S. Numerical simulation of caisson installation and dissipation in Kaolin clay and calcareous silt [J]. Bulletin of Engineering Geology and the Environment， 2018， 77（3）： 953-962.

[4] ?周松望， 張艷， 王棟. 正常固結(jié)黏土中扭矩對負壓沉箱承載力的影響分析[J]. 土木與環(huán)境工程學報（中英文）， 2020， 42（1）： 18-23.

ZHOU S W， ZHANG Y， WANG D. Influence of torsion on bearing capacity of suction caisson in normally consolidated clay [J]. Journal of Civil and Environmental Engineering， 2020， 42（1）： 18-23. （in Chinese）

[5] 陳佳瑩，滕竟成，吳則祥. 黏土中吸力式桶形破壞包絡面數(shù)值模擬研究[J/OL]. 土木與環(huán)境工程學報（中英文）. https：//kns.cnki.net/kcms/detail/50.1218.TU.20210125.1115.002.html.

CHEN J Y， TENG J C， WU Z X. Numerical modelling of failure envelope for suction bucket foundation in clay [J/OL]. Journal of Civil and Environmental Engineering. https：//kns.cnki.net/kcms/detail/50.1218.TU.20210125.1115.002.html. （in Chinese）

[6] 程健，姜君，王棟.考慮側(cè)壁土體脫離影響的桶基承載力[J]. 土木與環(huán)境工程學報（中英文）， 2021， 43（4）： 52-57.

CHENG J， JIANG J， WANG D. Capacities of bucket foundation considering the influence of sidewall-soil separation [J]. Journal of Civil and Environmental Engineering， 2021， 43（4）： 52-57. （in Chinese）.

[7] ?CHEN W， RANDOLPH M F. External radial stress changes and axial capacity for suction caissons in soft clay [J]. Géotechnique， 2007， 57（6）： 499-511.

[8] ?CHEN W， ZHOU H， RANDOLPH M F. Effect of installation method on external shaft friction of caissons in soft clay [J]. Journal of Geotechnical and Geoenvironmental Engineering， 2009， 135（5）： 605-615.

[9] WESTGATE Z J， TAPPER L， LEHANE B M， et al. Modelling the installation of stiffened caissons in overconsolidated clay [C]//Proceedings of the 28th International Conference on Offshore Mechanics and Arctic Engineering. Honolulu： American Society of Mechanical Engineers， 2009. 119-126.

[10] KIM J H， KIM D S. Soil displacement near a bucket foundation installed in sand by suction and jacking in a centrifuge [J]. Journal of Geotechnical and Geoenvironmental Engineering， 2019， 145（11）： 06019015.

[11] RAGNI R， BIENEN B， OLOUGHLIN C D， et al. Observations of the effects of a clay layer on suction bucket installation in sand [J]. Journal of Geotechnical and Geoenvironmental Engineering， 2020， 146（5）： 04020020.

[12] ?ZHOU H J， RANDOLPH M F. Large deformation analysis of suction caisson installation in clay [J]. Canadian Geotechnical Journal， 2006， 43（12）： 1344-1357.

[13] ?XIAO Z， FU D F， ZHOU Z F， et al. Effects of strain softening on the penetration resistance of offshore bucket foundation in nonhomogeneous clay [J]. Ocean Engineering， 2019， 193： 106594.

[14] ?JEANJEAN P. Setup characteristics of suction anchors for soft gulf of Mexico clays： Experience from field installation and retrieval [C]//Offshore Technology Conference. Houston， Texas， USA. Offshore Technology Conference， 2006： OTC18005.

[15] ?HU Y， RANDOLPH M F. A practical numerical approach for large deformation problems in soil [J]. International Journal for Numerical and Analytical Methods in Geomechanics， 1998， 22（5）： 327-350.

[16] WANG D， HU Y， RANDOLPH M F. Effect of loading rate on the uplift capacity of plate anchors [C]// Proceedings of the 18th International Offshore and Polar Engineering Conference. California： International Society of Offshore and Polar Engineers， 2008. 727-731

[17] ?WANG D， BIENEN B， NAZEM M， et al. Large deformation finite element analyses in geotechnical engineering [J]. Computers and Geotechnics， 2015， 65： 104-114.

[18] ?MAHMOODZADEH H， WANG D， RANDOLPH M F. Interpretation of piezo ball dissipation testing in clay [J]. Géotechnique， 2015， 65（10）： 831-842.

[19] ?WANG D， BIENEN B. Numerical investigation of penetration of a large-diameter footing into normally consolidated Kaolin clay with a consolidation phase [J]. Géotechnique， 2016， 66（11）： 947-952.

[20] STEWART， D. P. Lateral loading of piled bridge abutments due to embankment construction [D]. University of Western Australia， 1992.

[21] ?PURWANA O A， LEUNG C F， CHOW Y K， et al. Influence of base suction on extraction of jack-up spudcans [J]. Géotechnique， 2005， 55（10）： 741-753.

[22] ?WROTH C P. The interpretation of in situ soil tests [J]. Géotechnique， 1984， 34（4）： 449-489.

[23] DNV （Det Norske Veritas）. Geotechnical design and installation of suction anchors in clay： DNV Recommended Practice RP-E303 [S]. Hvik， 2005.

[24] ?SULLY J P， ROBERTSON P K， CAMPANELLA R G， et al. An approach to evaluation of field CPTU dissipation data in overconsolidated fine-grained soils [J]. Canadian Geotechnical Journal， 1999， 36（2）： 369-381.

（編輯 胡玲）