桶形基础二次加载的大变形有限元模拟

陈霄 郑敬宾 王栋

摘 要：桶形基础用于支撑海上风电升压站平台时，平台与电气设备分两阶段安装，中间存在数周到数月的间隔期。为探究桶形基础在二次加载阶段的峰值贯入阻力，发展了基于网格周期更新的有效应力大变形有限元方法，避免基础贯入引起周围土体单元扭曲，再现了黏土中基础“贯入固结二次加载”的全过程。对比已有的离心模型试验，验证了发展的大变形数值方法的可靠性。通过比较初始贯入阻力与固结后的贯入阻力，发现桶基与平台的安装间隔期能够有效增加桶基承载力，固结后的贯入阻力峰值仅需要很小的基础位移即可达到。开展变动参数分析，考察影响峰值贯入阻力比的主要因素，提出考虑基础尺寸、土体渗透性与固结时间的歸一化时间参数，并给出峰值贯入阻力比随归一化时间增长的趋势线。

关键词：桶形基础;有限元分析;固结;黏土;贯入阻力

中图分类号：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阶段，桶基与平台贯入土中预定深度，施工船返回港口，桶基静置数周至数月;第2阶段，各种大型电气模块运输至场地，并吊放到平台上，桶基实质上承受二次加载。第1阶段的桶基贯入在黏土中引起超静孔压，之后孔压部分或完全消散[2-3]，导致第2阶段开始时桶基周围黏土的不排水强度高于第1阶段结束时。

类似于一般的浅基础承载力问题，常规有限元方法被广泛用于分析不排水条件下桶基的竖向承载力及复合承载力[4-6]，但这些研究均假定桶基预埋就位，不考虑桶基安装过程对承载力的影响。这是因为桶基长距离贯入导致侧壁端部周围土体单元的严重扭曲，常规有限元方法无法捕捉这一大变形过程。

研究者借助离心机试验[7-11]或大变形有限元模拟[12-13]再现了桶基的贯入，探究了土层分布、安装方式、基础尺寸等一系列影响桶基贯入阻力的因素，提出了确定第1阶段贯入阻力的方法。少量的数值模拟还进一步给出了桶基安装后的孔压消散过程[3]，但这些结果仅评估了地基土的固结程度，不能预测升压站平台第2阶段的抗压承载力。Jeanjean[14]在墨西哥湾海域进行了现场试验，并总结了历史上该区域的施工记录，得到桶基安装就位时的贯入阻力以及固结一段时间后的抗拔力，发现桶基在安装后存在显著的固结效应，最终抗拔力最高可达初始贯入阻力的1.6倍，但其试验针对的是抗拔力，而不是抗压承载力。

鉴于与桶基二次加载有关的研究较为有限，笔者发展了有效应力形式的大变形有限元方法，模拟黏性土中桶基的“贯入固结二次加载”完整过程，探索桶基初始贯入和部分排水固结后的超静孔压分布，进而考察固结效应对二次加载的影响，合理评估二次加载的竖向承载力。

1 研究方法

1.1 大变形有限元方法

桶形基础的安装是典型的岩土工程大变形问题，因此使用一种名为RITSS的大变形有限元方法（Remeshing and interpolation technique by small strain）。这种方法最初由Hu等[15]提出。Wang等[16-17]将其与有限元商业软件Abaqus结合并扩展至有效应力分析。RITSS方法中，整个计算过程被划分为一系列小变形计算步，通过网格重划分及场变量映射，避免了单元的严重扭曲，具体细节可参见文献[17-18]。

采用基于Abaqus的RITSS方法，将桶基的“贯入固结二次加载”过程简化为二维轴对称问题进行分析，建立如图1所示的数值模型。桶基的长度、直径和厚度分别为L、D和t，支撑电气平台的桶基长径比L/D一般为1左右。为避免边界效应，土体的径向与竖向边界分别取5D与8D。约束底边径向边界的水平与竖向位移，竖向边界的水平位移。桶基周围附近土体的网格尺寸为1.5t。为提高计算效率，并避免贯入初期土体表面单元的剧烈变形，将桶基预埋一定深度。后文将证明该设置对后续贯入阻力的影响可以忽略。采用二阶轴对称缩减积分单元划分土体，以保证映射精度[17]。

采用罚函数法定义“桶土”界面的摩擦力，在有效应力分析中，根据库伦摩擦定律，摩擦力取决于法向有效应力及界面的库伦摩擦系数。安装造成桶基周围土体的应变软化，设计时一般假定周围土体完全扰动。Koh等[3]采用大变形方法模拟桶基贯入时取摩擦系数为0.3tan ′/St，其中′为有效内摩擦角，St为灵敏度，得到了与离心机试验结果相近的贯入阻力。笔者的数值分析也采用相同的摩擦系数。

以往离心机试验[7-11]比较了吸力式安装（即通过气泵从桶盖与土体之间抽水以施加吸力）与纯压入式安装（即采用位移控制使桶基贯入）两种桶基贯入过程，发现安装方式对黏土中桶基贯入阻力的影响可以忽略不计。因此，通过位移控制模式实现桶基的初始贯入和二次加载。在固结阶段，桶基上的竖向力保持为定值，竖向力大小等于当前深度的初始贯入阻力。在整个“贯入固结二次加载”过程，仅允许土体表面边界排水，其余各边均为不排水边界。

1.2 土体本构模型

采用经典的修正剑桥模型描述黏性土的有效应力与应变之间的关系，该模型已被成功用于分析静力触探仪、桩靴基础的大变形贯入[18-19]。考虑两种广泛用于离心机试验中的黏性土，即UWA高岭土和Malaysian高岭土，相关计算参数总结于表1。側压力系数K0=1

2 数值模拟验证

目前，还没有公开发表的桶形基础二次加载模型试验或现场监测数据。通过对比桶基首次安装时的贯入阻力，验证数值模拟的可靠性。Westgate等[9]完成了6个桶基在UWA高岭土中贯入的离心机试验。使用RITSS方法模拟其中名为Type1的试验。桶基的原型尺寸为D=11.3 m，L=6 m，t=0.05 m。试验前和试验后分别进行两次T-bar试验，得到土体的平均抗剪强度su=10+2.8z kPa。利用两次循环插拔T-bar试验获得土体的灵敏度St为2.5和3。大变形有限元分析中使用St=3，确定桶土摩擦系数。试验中测得γ′=5.9 kN/m3，根据式（1）推导得到修正剑桥模型中的su=1.44z kPa。为便于比较，将贯入阻力F归一化为F/suA，A为桶形基础竖向投影面积，A=πD2/4。图2比较了离心模型试验[9]、有效应力形式的RITSS法、总应力形式的耦合欧拉拉格朗日有限元法（Coupled Eulerian-Lagrangian，CEL）结果[13]、DNV规范的预测[23]。图中的纵坐标为归一化贯入深度d/D，即贯入深度与桶基直径之比。

DNV规范与CEL的结果均基于总应力分析，假定安装时地基土完全不排水。其中DNV规范建议周围土体的强度折减系数α=1/St=0.33，表示桶的最大侧壁摩阻力为αsu。文献[13]中CEL分析使用的系数为α=0.2，故图2也给出了α=0.2时的DNV方法预测结果。由图2可见，有效应力形式的RITSS方法和总应力形式的CEL方法给出的结果都与试验数据接近，且均在DNV公式预测的α=0.2和α=0.33两条曲线范围内。采用RITSS方法时，预埋使得计算曲线从0.2D开始，因此，无法比较0～0.2D范围内的结果。然而，当桶基的贯入深度超过0.25D后，贯入阻力曲线与试验非常接近，表明一旦土体破坏模式被完全激发，预埋对贯入阻力不会产生明显影响，这与已有结论一致[18-19]。图2的比较也表明大变形有限元分析中的摩擦系数取值合理，可以在很大程度上考虑安装扰动对贯入阻力的影响。

3 结果分析与讨论

实际工程中钢制桶基的壁厚一般小于0.05 m。试算发现，当桶基侧壁厚度较小时，有效应力形式的大变形分析不易收敛，后续分析中统一取壁厚为0.19 m。Xiao等[13]比较了侧壁较厚（t=0.2～0.6 m）与较薄（t=0.05 m）的桶基端部承载力系数，发现其范围在10～12之间，表明一定范围内壁厚的增加没有改变桶基贯入时端部以下土体的破坏模式，而侧壁附近的土体破坏模式不依赖壁厚。因此，增加壁厚不会显著影响大变形结果。

表2总结了变参数分析中开展的算例。为探究固结系数cv的影响，算例A～算例D分别模拟了桶基在UWA高岭土和Malaysian高岭土中的“贯入固结二次加载”过程。这两种土的本构模型参数接近，但渗透系数相差很大，相同应力水平下Malaysian高岭土的cv约为UWA高岭土的10倍。算例A～算例D还将用于考察首次安装深度zp对二次加载的影响。算例E和算例F将考察桶基直径D对二次加载的影响。表2中的固结时间tc在7～360 d之间变化，其下限是实际工程中最短的施工间隔，而上限人为取得较大，以便展示充分固结造成的二次加载效应。

修正剑桥模型并不能考虑土体结构性的影响，分析中摩擦系数始终保持为0.3tan ′/St，这也意味着固结前后的土体灵敏度未改变，即不考虑扰动的土体结构在固结过程中部分恢复，所以，计算结果偏于安全。以下分析中统一取典型值St=2.5。

3.1 固结时间对二次加载的影响

算例A、C（D=9.5 m，首次贯入深度zp=0.9D），固结7、30、60、120、270 d后，二次加载贯入阻力随深度的变化如图3所示。为了直观显示固结时间的影响，图中包含了初始连续贯入阻力曲线作为参考。参考Wang等[19]对桩靴基础贯入阻力的处理方法，利用峰值贯入阻力比Ncp/Ncr表征贯入阻力的改变，如图3（a）所示。Ncp为二次加载时归一化的峰值贯入阻力，Ncr为Ncp出现深度处归一化的初始贯入阻力。

从图3可以看出：

1）桶基的二次加载阻力在很短位移内即可达到峰值，Ncp出现时的深度仅略大于初始贯入结束时的深度zp。

2）UWA高岭土中即使固结时间仅有7 d，峰值贯入阻力也会增加约20%，见图3（a）。如果固结时间为7～120 d，二次加载的Ncp随固结时间增大，说明固结导致的桶基竖向承载力显著提高。但120 d后的Ncp基本不再增加。

3）二次加载的贯入阻力曲线达到峰值后迅速降低，并向初始连续贯入阻力曲线靠近。这表明，尽管施工间隔期间桶基周围土体固结，强度增加，但强度增加土体的范围有限。这一特征与桩靴基础的二次贯入类似[19]。

4）Malaysian高岭土的结果表现出与UWA高岭土类似的趋势，如图3（b）所示，Ncp/Ncr最大约为1.35。然而，固结60 d后Ncp就不再有显著增长。产生这种差异的原因是两种土的固结系数cv不同。Malaysian高岭土的固结系数远高于UWA高岭土，初始贯入造成的累积孔压在固结阶段消散得更快，桶基周围土体更快达到强度增长极限。

图4展示了固结过程中桶基附近3个位置的土体超静孔压消散过程。3个位置分别为桶壁端部中心、端部上方0.02zp处与桶壁紧邻的内侧与外侧。为避免固结初期孔压重分布对消散曲线的影响，利用Sully等[24]提出的方法确定每一个位置的名义初始孔压Δumax，得到归一化孔压Δu/Δumax随固结时间的变化曲线。图4表明，3个位置处的超静孔压消散规律几乎一致。对于UWA高岭土，孔压在最初的120 d迅速消散，120 d后的消散速度逐渐减缓，这与峰值贯入阻力的增长规律一致（图3）。同样，Malaysian高岭土中超静孔压在前60 d急剧下降，且消散速度比UWA高岭土更快。

由图6可以看出，zp/D在0.7～1.5范围内时，两种高岭土参数计算得到的峰值贯入阻力比随归一化时间的变化趋勢一致。曲线整体上可以分成两段：当T小于0.015时，Ncp/Ncr随T大致呈线性增长关系;当T大于等于0.015时，Ncp/Ncr趋近于一个稳定值。为了建立Ncp/Ncr与T的关系，在图6中添加了趋势线。地基土充分固结完成时，算例A的Ncp/Ncr达到1.36，算例F的Ncp/Ncr约为1.26，其他算例结果介于二者之间。也就是说，完全固结时二次加载得到的峰值贯入阻力最少增长约26%。

实际条件下，峰值贯入阻力的增长可能超出图6的有限元结果，这是因为在固结过程中孔压消散以及土体结构性的恢复都可以造成后续承载力的上升，且二者的贡献程度相当[14]。受限于本构模型，本文的有限元分析还不能反映土体结构性的恢复，所得峰值贯入阻力增长程度偏于保守。

4 结论

采用有效应力形式的大变形有限元方法模拟海上风电升压站平台的安装过程，再现了桶基“贯入固结二次加载”的全过程。数值模拟中考虑了两种力学性质接近，但固结系数相差10倍的高岭土，考察了固结时间、首次安装深度、桶基直径等关键变量对二次加载时贯入阻力的影响。主要结论如下：

1）当固结时间相同时，高固结系数土中二次加载的贯入阻力增长更迅速，但固结系数对峰值贯入阻力比极限值的影响不大，完全固结时的峰值贯入阻力比在1.26～1.36之间。

2）式（4）定义的归一化固结时间较好地反映了峰值贯入阻力比的变化趋势。采用归一化固结时间表征固结程度时，峰值贯入阻力比随归一化固结时间约为线性增加，直至达到极限值。

3）实际施工中，应充分利用二次加载前的施工间隔时间，发挥地基土固结后的强度，从而提高二次加载时的极限承载力，达到优化基础尺寸，降低工程造价的目的。参考文献：

[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]  周松望， 张艳， 王栋. 正常固结黏土中扭矩对负压沉箱承载力的影响分析[J]. 土木与环境工程学报（中英文）， 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] 陈佳莹，滕竟成，吴则祥. 黏土中吸力式桶形破坏包络面数值模拟研究[J/OL]. 土木与环境工程学报（中英文）. 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] 程健，姜君，王棟.考虑侧壁土体脱离影响的桶基承载力[J]. 土木与环境工程学报（中英文）， 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.

（編辑 胡玲）