ALE方法评估黏土场地的自升式平台桩靴承载力.pdf

第38 卷第5期2023年10 月文章编号：10 0 1-450 0(2 0 2 3)0 5-0 0 0 7-0 5中国海洋平台CHINA OFFSHORE PLATFORMVol.38 No.5Oct.,2023DOI:10.12226/j.issn.1001-4500.2023.05.20230502ALE方法评估黏土场地的自升式平台桩靴承载力冯士伦，陈晓桐，王进（1.天津大学建筑工程学院，天津30 0 0 7 2；2.中国土木工程集团有限公司，北京10 0 0 38)摘要：为研究自升式平台桩靴在黏土场地中的极限承载力，建立Abaqus桩靴-土模型，采用任意拉格朗日-欧拉（ArbitraryLagrangian-Eulerian，A LE）方法对桩靴的贯入过程进行模拟，提出大值法、中值法、小值法等3种数据处理方法，并拟合桩靴承载力估算公式。将模拟计算结果与SNAME规范公式计算结果进行对比分析，对比结果证明ALE方法模拟桩靴贯入的可行性。为提高工程的安全性，提出对SNAME公式计算结果进行折减的具体公式，以使其与ALE小值法处理结果更吻合。关键词：桩靴贯入；任意拉格朗日-欧拉（ArbitraryLagrangian Eulerian，A LE）方法；极限承载力中图分类号：TE951文献标志码：AEvaluation of Bearing Capacity of Jack-up Platform Pile Bootin Clay Using ALE MethodFENG Shilun,CHEN Xiaotong,WANG Jin?(1.Schoolof Civil Engineering,Tianjin University,Tianjin 300072,China;2.China Civil Engineering Construction Corporation,Beijing 100038,China)Abstract:In order to study the ultimate bearing capacity of jack-up platform pile boot in the clay,theAbaqus pile boot-soil model is established,and the penetration process of the pile boot is simulated by ArbitraryLagrangian-Eulerian(ALE)method.Three data processing methods of large-value method,middle-valuemethod and small-value method are put forward,and the formulae for estimating the bearing capacity of pileboot are fitted.The simulation results are compared with those of the SNAME formula,and the feasibility ofALE method to simulate the penetration of the pile boot is proved.Considering the safety of the practicalengineering,in order to modify the calculation result of SNAME formula and make it more consistent with thatof ALE small-value method,a reduction method is proposed for the calculation formula of SNAME.Key words:pile boot penetration;Arbitrary Lagrangian-Eulerian(ALE)method;ultimate bearing capacity基金项目：国家自然科学基金（编号：52 0 7 8 336）作者简介：冯士伦（197 6 一），男，副教授，主要研究方向为海洋工程80 引 言自升式钻井平台是海洋油气资源开发的重要装备1。自升式平台就位作业频繁，插桩安装过程安全风险大，插桩事故时有发生。比如2 0 2 1年，马来西亚钻井平台NAGA7号在预压插桩时桩腿沉陷，导致平台倾覆，所幸无人员伤亡，但也造成该项目巨大损失2 。因此，有必要对桩靴贯入过程进行更全面、更准确的研究，以避免插桩事故发生，保证安全生产。已有学者对桩靴贯入问题进行模拟和讨论。HOSSAIN等3-4结合离心模型试验与大变形有限元分析方法研究桩靴贯入土体破坏模式及承载特性。SHIAU等5 采用有限元方法模拟上砂土下黏土的土层与桩靴的相互作用，调查桩靴几何参数的影响。BURD等6 采用有限元方法研究上砂土下黏土的土层中桩靴贯入的极限承载力及其影响因素。张浦阳等采用饱和黏土条件下土-桩靴系统的有限元模型对桩靴贯入均质与非均质黏土时的桩靴承载力系数进行数值计算，提出较传统方法更合理的计算方法。陈洋彬等8 通过有限元分析研究桩靴贯人形成孔洞对承载力的影响。范庆来等9-10 1采用欧拉-拉格朗日耦合（Coupled Eulerian-Lagrangian,CEL）算法模拟上砂土下黏土双层地基中桩靴的贯入过程，研究土体破坏模式与各参数对贯入阻力的影响。陈大江等11以实际工程平台为对象，采用数值模拟方法研究安全插桩深度等桩靴稳定性问题。王东石等12 设计插桩模型试验装置，研究桩靴贯人的穿刺现象。冯加伟等13使用CEL方法模拟桩靴贯人，研究土体响应及插桩深度等问题。张湖宾等14使用CEL方法研究开孔、裙摆等改进对桩靴穿刺现象的影响。综上所述，自升式平台桩靴贯入问题的研究已有很大进展。研究人员较多采用有限元方法进行研究，尤其是采用CEL方法。不过，采用任意拉格朗日-欧拉（Arbitrary Lagrangian-Eulerian，A LE)方法进行自升式平台桩靴贯入的研究成果较少，该方法可避免网格畸变问题，有利于大变形分析，对自升式平台桩靴贯人问题有较好的适用性。希望本研究能为该课题的进展提供一个不同的视角，为自升式平台桩靴的安全贯人、稳定工作提供参考。中国海洋平台1ALE方法简介ALE方法是纯拉格朗日分析法与纯欧拉分析法的结合，一方面引进拉格朗日分析法对边界运动的处理思路，另一方面借鉴欧拉分析法使网格单元独立于实体物质存在。这种结合克服了拉格朗日分析法不利于大变形计算的特点，也克服了欧拉分析法在物质边界处理上的困难，是一种十分利于大变形分析的有限元分析方法。在进行模型分析时，在每一步计算中先进行拉格朗日运动计算，然后重新划分单元网格，再进行下一步运动计算，使网格的大小和位置随着分析步不断适应变形，调整后再进行下一步计算，以避免网格发生畸变导致计算失败或中断，从而适应大变形分析。如图1所示：图1（a）中圆圈区域尖端的网格是未使用ALE方法的网格变形情况，显然四边形单元压缩成三角形单元，网格严重畸变；图1（b）中圆圈区域尖端的网格是使用ALE方法后的网格变形情况，显然通过变形后网格的重新划分，网格形态更好。ALE方法除了大变形后网格形态好以外，可更精准地定义边界条件和接触条件，避免大变形引起的网格纠缠和局部误差，从而提高计算精度。（a）未使用ALE方法网格变形（b）使用ALE方法网格变形图1网格变形情况对比2应用算例为考察ALE方法对自升式平台桩靴贯人问题的有效性，针对工程中常见的饱和软黏土层和桩靴，建立桩靴-土模型，进行自升式平台桩靴贯入分析。桩靴采用如图2（a）所示典型的纺锤状桩靴，具体尺寸如图2（b）所示。基于模型与边界条件的对称性，为减少计算量、缩短计算时间，建立桩靴与土的1/4模型。在插桩过程中，桩靴变形很小，定义为刚性体。为消除边界效应，土体模型水平向尺寸取10 0 m、竖直向尺寸设50 m，整体有限元模型图如图2（c）所示。土体侧向边界约束水平位移，土体底边界约束竖向位移。在桩靴刚体参考点上施加恒定速度模拟桩靴贯入。土体中心靠近第38 卷第5期第5期桩靴的部分变形较大，将此区域土体作为自适应网格，远离桩靴的部分变形较小，在此处采用较粗大冯士伦，等ALE方法评估黏土场地的自升式平台桩靴承载力的网格，网格划分如图2（d）所示。这样既可实现大变形分析，又可在一定程度上减小计算量。18(a）桩靴模型为考察黏土强度指标与桩靴竖向承载力之间的关系，选取黏土不排水抗剪强度Cu如表1所示，黏土浮容重取8 kN。表1黏土强度指标编号黏土1黏土2 黏土3黏土4黏土5不排水抗剪强度/kPa10152025303计算结果处理3.1数据处理方法ALE方法是一种动力分析方法，因此计算的输出结果曲线呈现波动。为方便处理和应用波动的曲线，提出小值法、中值法和大值法等3种数据处理方法。3种方法对应取原始数据的最小值、中间值和最大值分别拟合该曲线，如图3所示。小值法具体做法为找出波动曲线各阶段的最小值点，将其连线，得到相对平滑的包络线，可体现曲线走势，大值法以相同方法取最大值，中值法则取小值法与大值法的平均值连线。由于在不同C土体内桩靴贯人的数据处理结果相似，此处作为举例，图3仅列出黏土3中的处理结果。10204-6一8-10L图3波动曲线3种处理方式（b）桩靴尺寸示例图2 桩靴-土模型示例由图3可知：经过3种方法处理后的曲线光滑明了，方便使用；在贯入初期即贯入深度为2 m以内的阶段，采用小值法与大值法承载力相差3MN左右；在贯人中后期，采用小值法与大值法承载力相差12 13MN。在实际工程应用中，偏于安全考虑，建议选用小值法进行处理，若结合实际考虑工程安全裕量较大，也可使用中值法和大值进行处理。3.2ALE方法结果偏于安全考虑，对不同Cu土体内桩靴贯入的ALE方法计算结果统一采用小值法进行处理，处理结果如图4所示。承载力/MN1020304050607080C=10C-15一2.C20-C25.C.30/一4一6一8-10L承载力/MN3040原始数据小值涛中值涛大值法（c）整体模型5060(d）网格划分图4不同C黏土ALE结果对比对于黏土1黏土 5，Cu由小到大的增量均为5kPa，桩靴承载力也随之由小到大增长，并且不同曲线间的增量也大致相同。在贯入初期，当贯入深度约3m时，相邻编号的2 种土体中桩靴承载力相差约5MN；在贯中期，即当深度约6 m时，差值约8 MN；在贯人后期，当深度为9 10m时，差值约10 MN。在黏聚力不同的土体中，桩靴贯人深度越大，桩靴的承载力提高得越10多。由图4可知，曲线起始段呈弧线状，后段形似直线，直线起点约在贯人深度为4m的位置（桩靴顶面入土的位置）。这是由于桩靴在此位置的横截面面积最大，与土体接触的侧向面积也最大。超过此位置后，桩靴承载力开始随贯人深度呈线性增大。图4直线段数据在数值上可拟合为以下基于Cu和贯人深度的方程：Fv=(0.205Cu+0.45)D+0.01C+0.15Cu+5.5)式中：Fv为桩靴承载力；D为贯入深度。根据式（1）可估计桩靴在其他不同Cu黏土中的承载力，为预估桩靴承载力提供参考。4与SNAME结果的对比分析国内桩靴承载力分析大多借鉴SNAME（T h eSociety of Naval Architects and Marine Engineers)规范。SNAME出版的 Recommended Practice forSite Specific Assessment of Mobile Jack-up Units提出预测桩靴插深的方法及桩靴地基承载力的计算公式15。根据该规范，在黏土环境中地基极限垂直承载力的计算公式为Fv=(CuNesede+po)A式中：N。为承载力因子，取5.14；s。=1十（Ng/N。)(B/L)为承载力形状系数，其中，N。=e a n g tan（45十/2)为承载力系数（其中为内摩擦角），B为桩靴直径，L为基础长度，对于圆形截面基础B=L；d。=1十0.4(D/B)为承载力深度系数；p为深度D处的有效覆土压力；A为桩靴与土接触的最大截面。般齐麟等16 1采用RITSS方法探讨插桩过程土体变形规律且计算插桩各深度承载力，并与SNAME规范结果进行对比分析以验证该方法的有效性和合理性。凌有临等17 采用Abaqus 的 CEL计算方法模拟某桩靴在单层土中的插桩过程，并将其结果与土工试验结果和 SNAME规范计算结果进行对比，指出其中的区别与联系。包彩虹等18 7对比SNAME与其他规范的修正公式，并基于此建立插桩计算软件。显然，SNAME规范结果是个很好的参照物。以黏土3计算结果为例，将ALE方法的计算结果曲线与SNAME公式的计算结果曲线进行对比分析，如图5所示。中国海洋平台-2-4一6一8(1)-10L图5ALE数值模拟与SNAME公式结果对比从整体上来看，ALE方法的计算结果与SNAME公式的计算结果相似，都是随着贯人深度的增大，桩靴承载力提高；不过，在大部分情况下SNAME公式的计算结果较ALE方法的计算结果大，即在贯入同一深度时，SNAME估计承载力大于ALE计算承载力，仅在桩靴贯人深度较大（大于7 m）时，ALE大值法的结果大于SNAME公式结果。鉴于在实际工程中自升式平台插桩事故时有发生，事故原因往往是桩靴承载力估计出现偏差。基于上述计算结果，建议将SNAME桩靴承载力公式的结果进行折减，使其接近ALE小值法的结果，以提高桩靴承载力计算结果的安全性，即(2)建议把SNAME桩靴承载力计算公式调整为Fv=(CuNesede+po)Ao.3 VD式中：D的最大值为9 m。式（3）即在原 SNAME桩靴承载力计算公式上乘以一个不大于1的折减系数。折减后的SNAME桩靴承载力估计值与ALE小值法结果更接近，对自升式平台插桩工程的安全更具有实际指导意义，可作为一种安全措施参考。黏土1黏土5各种土体中的比较结果相似，为简洁起见，仅在图6 中给出黏土3中SNAME桩靴承载力折减后的结果与ALE小值法结果的对比。承载力/MN1020一2/-4一6一8-10L图6 折减后SNAME公式结果与ALE方法结果对比第38 卷第5期承载力/MN10203040一小值法中值法大值法.-SNAME公式3040一ALE方法-SNAME公式折减5060(3)50第5期5结论ALE方法是纯拉格朗日分析法与纯欧拉分析法的结合。使用ALE方法模拟不同黏土中自升式平台桩靴的贯入过程，同时与SNAME规范公式计算结果进行相应的对比分析，得出主要结论如下：（1）使用ALE方法模拟桩靴贯入可反映土体属性对桩靴承载力的影响，桩靴承载力随着Cu的增大而增大，说明采用ALE方法模拟桩靴贯人的可行性。（2）对于ALE方法结果出现的波动问题，提出小值法、中值法和大值法等3种处理方法，并建议根据不同的工程安全裕量采用不同处理方法。（3）观察到桩靴顶面入土后，桩靴承载力-贯人深度曲线接近直线，拟合出基于黏土Cu与贯人深度的承载力公式，采用该公式可预估不同C黏土中的桩靴承载力，方便在可行性研究阶段对桩靴承载力进行评估。（4）采用SNAME公式计算的桩靴承载力在大部分情况下都较ALE模拟结果大，为安全起见，拟合出SNAME公式的折减系数，使折减后的结果与ALE小值法结果吻合，为桩靴承载力的安全评估提供参考。参考文献1汪张棠，赵建亭.我国自升式钻井平台的发展与前景：中国海洋平台，2 0 0 8（4）：8-13.2突发!马来西亚一座海上钻井平台倒塌沉没EB/OL.2 0 2 2-7-6 .h t t p s:/w w w.s o h u.c o m/a/464583790_726570.31HOSSAIN M S,RANDOLPH M F.New mechanism-based design approach for spudcan foundations onsingle layer clay J.Journal of Geotechnical andGeoenvironmental Engineering,2009,135(9):1264-1274.冯士伦，等ALE方法评估黏土场地的自升式平台桩靴承载力11.4HOSSAIN M S,RANDOLPH M F.Deep-penetratingspudcan foundations on layered clays:numericalanalysisJ.Gotechnique,2010,60(3):171-184.5SHIAU J S,LYAMIN A V,SLOAN S W.Bearingcapacity of a sand layer on clay by finite element limitanalysisJ.Canadian Geotechnical Journal,2003,40(5):900-915.6 1BURD H J,FRYDMAN S.Bearing capacity of plane-strain footings on layered soils J.C a n a d ia nGeotechnical Journal,1997,34(2):241-253.7张浦阳，于晓洋，丁红岩.海上自升式钻井平台插桩阶段桩靴承载力计算J.石油勘探与开发，2 0 11,38(5):613-619.8陈洋彬，郑敬宾，王栋.弱超固结黏土中桩靴贯人形成孔洞对承载力影响J.海洋工程，2 0 2 139（1）：112-120.9 范庆来，苏煜茹，张炳磊.双层地基中纺锤形桩靴基础贯人阻力及机理分析.应用基础与工程科学学报，2020,28(2):452-459.10范庆来，张炳磊，贾小伟.双层黏土地基方形桩靴贯入过程与穿透机制.土木建筑与环境工程，2 0 18,40(1):62-67.11陈大江，陈洪飞，沈佳轶，等.舟山某海域海上自升式平台桩基稳定性研究.科技通报，2 0 2 1,37（9）：6 9-7 4.12王冬石，袁烨，徐文祥，等.非连续压载下黏土层中桩靴承载力变化规律.中国海上油气，2 0 18,30（5）：145-150.13冯加伟，周立臣，梁思颖.钻井船插桩对海洋平台桩基影响的数值研究J.中国海洋平台，2 0 2 2,37（1）：43-47.14张湖宾，陈晨.基于CEL法的海洋自升式平台插桩研究J.舰船电子工程，2 0 2 2,42(4)：18 6-19 0.15The Society of Naval Architects and Marine Engineers.Recommended Practice for Site Specific Assessment ofMobile Jack-up UnitsS.2nd ed.2002.16殷齐麟，董胜，樊敦秋.复杂地层中自升式平台插桩的数值模拟.工程力学,2 0 16,33（9）：2 0 4-2 11.17凌有临，杨栋，张建.自升式平台插桩过程CEL分析及验证门.科学技术与工程，2 0 15,15（7)：7 0-7 3.18包彩虹，王俊杰：自升式风电安装船在不同土质插桩深度研究.中国港湾建设，2 0 19,39（8）：46-52.(上接第6 页）1oBREBBIA C A,WALKER S.Response of offshorestructures J.D y n a m ic A n a ly s is o f O ffs h o r eStructures,1979,3(4):292-320.11王卫远，何倩倩，周鹏飞，等.福建南日群岛海域波浪数值模拟研究.海洋预报，2 0 13,30（5）：2 6-30.12袁方超，吴向荣，卢君峰.福建中南部沿海风暴潮统计特征分析.海洋预报，2 0 18，35（3）：6 8-7 5.13张军，许金电，郭小钢.福建沿海海域波浪能资源分析与评价门.台湾海峡，2 0 12,31（1）：130-135.