联合载荷作用下不排水桩靴的承载能力

王立忠等人，中国海洋工程，第25卷，第一期，15-301中国海洋工程，第25卷，第一期，15-30页2011年中国海洋工程协会和柏林海德堡斯普林格出版社文件编号：10.1007/s13344-001-0002-0联合载荷作用下不排水桩靴的承载能力王立忠，舒恒，李玲玲，国振浙江大学土木工程与建筑学院，中国，杭州，301158中国城市建设控股集团第二高速公路顾问有限公司，中国，武汉430056(12月23日收到2009年12月23日收到初版；2010年8月17日收到修改版本，2010年9月27日接受)摘要采用一系列的三维有限元分析方法研究仅在垂向、水平、弯矩力分别单独作用、三种力两两联合作用、三种力联合作用下桩靴的极限承载能力。详细研究了埋置率和土壤非均质程度对承载能力的影响。桩靴在不同载荷单独作用时承载能力用无量刚的承载能力系数表示，并与已发表的论文结果做对比。组合载荷作用下的桩靴极限承载能力由无量纲化的和标准化后的破坏包络线呈现。论文所得的破坏包络线和已发表的论文结果之间的比较说明了土壤的埋置率和非均质情况决定了破坏包络线的大小和形状。关键词：软质粘土；桩靴；承载能力；联合载荷；有限元第一章介绍自升式海洋平台由于其在油田开发的低成本、良好的操纵性、生产的高效性和作业区域的灵活性而广泛应用于离岸油气的开发，尤其是边缘油田。典型的自升式海洋平台包括上部主体、桩腿、桩靴，上部的主体由三根相互独立的桁架式桩腿支撑，每一根桩腿又支撑在倒锥形的结构上，成为桩靴。图1b所示为一桩靴的典型结构，近似为圆形，直径达20米。挪威船级社提出，对于固结软粘土和非固结软粘土等土壤，由于土壤失效造成的位移很大，所以桩基础的承载能力的预测难度将会很大。在这些土壤情况下，最初的土壤失效都应当认定为地基稳定性的失效。在生命周期中，桩靴基础将经受垂直、水平、和倾覆力矩载荷的作用，这些载荷主要由平台的自重和环境载荷造成。为了预测自升式海洋平台的状态，必须充分研究在各种联合载荷作用下桩靴基础的响应情况。在海洋工业界，海森（1970）、威斯克（1975）提出的传统承载能力理论在海洋岩土工程中应用较为广泛。近几十年来数值分析已被越来越多的研究者认同，尤其是不排水离岸浅滩基础的极限王立忠等人，中国海洋工程，第25卷，第一期，15-302承载能力。已知的关于数值分析的文章，大多局限于垂直承载力或安装过程的载荷(胡锦涛和兰多夫，1998；侯赛因等人，2004；侯赛因，2004)，组合载荷作用下的承载能力研究较少。此外，海洋沉积物通常是非均质的，抗剪强度会随深度增加。非均匀粘土基础在垂直载荷作用下的承载能力研究(Davis和Booker，1973；Houlsby和Wroth，1983；Kusakabe等人)较为广泛。然而直到现在，组合载况下的承载能力研究依然很少，并且也仅仅限于长条状和圆形的桩基础(Bransby和兰多夫，1998；Taiebat和卡特，2000；Gourvenec和兰多夫，2003)。即使对于长条形的桩靴基础，倾斜和偏心载荷作用下的承载能力理论依然不可靠，抗剪强度随深度增加的地基土壤承载能力(Ukritchon等人，1998；Gourvenec和Randolph，2003)结果偏小。在联合载荷作用下的桩靴基础极限承载能力可以由位于V-H-M载荷空间破坏包络线很好的呈现(Roscoe和Schofield，1957；Bransby和Randolph，1998；Taiebat和Carter，2000；GourvenecandRandolph，2003；Gourvenec，2007、2008)。在不考虑埋置率的情况下，破坏包络线将会有相似的形状的猜想被广泛提出(Martin1994；Bransby和Randolph，1994)。在埋置支撑点的设计参考规范中指出，地基表面的破坏包络线会因为深度因子的不同扩大或缩小。同时，由于水平和弯曲方向载荷的相互作用，实际工程中破坏包络线也是不同的。为了研究上述猜想的合理性以及深度因子的具体影响机理，研究人员利用有限元的方法探索桩靴基础破坏包络线的形状尺寸与一系列埋置率的关系，这项研究已经利用商业有限元软件ABAQUS开展。桩靴基础仅受垂向力或水平力或弯矩时的承载能力研究已经很充分；在V-H、V-M、H-M和V-H-M联合荷载作用下，利用无量纲化和规范化的破坏包络线进行的研究也已经成熟。此外，研究对不同的埋置率和地基土壤的匀质性造成的不同的破坏包络线也做了充分的比较，所得的结论有利用桩靴的设计和安装。王立忠等人，中国海洋工程，第25卷，第一期，15-303第二章有限元模型和分析方法2.1几何和材料参数如图2所示，研究中采用直径为D的埋置式圆形桩靴。在软海床上，桩靴基础可能刺入土壤中长达1-2倍的桩靴直径，因此假设有充分的土壤回流。有限元模型中忽略了桩腿，这会带来一些简化并且在设计中会留有一定的安全余量。环境载荷持续时间的很短，所以应该假定土壤不排水。不排水土壤的响应可由线性理想弹塑性理论分析，不排水土壤的属性可由杨氏模量E和泊松比υ定义。失效由特雷斯卡准则定义，最大剪应力仅限于不排水剪切强度（SU）。杨氏模量和不排水抗剪强度之比为一常数，E/SU=500，泊松比υ取为0.49.对于均质土，不排水抗剪强度SU取为10千帕。对于正常固结粘土（非匀质粘土），假定剪切强度与深度呈线性变化，SU=1.1Z千帕。土壤的有效单位重量γ'承力为6.8Kn/m3。图3所示为论文中使用的土壤和桩靴基础简图。其中d是桩靴基础的埋置深度，SU0是参考剪切强度，k是剪切强度随深度z的变化梯度。正常压密粘土，取图3b中D/4深度处的剪切强度为SU0，SU0代表最宽埋置截面以下且在半径以内的平均剪切压强。2.2有限元网格因为桩靴基础是轴对称的，所以选择半圆柱形三维有限元网格。海床的半径和深度分别为7.5D和5.0D。外板边界在垂直方向没有位移，底部的网格在所有三个方向是固定的。假设地基土与桩靴充分结合，防止桩靴和土壤滑动或分离。有限元网格密度等参数取决于埋置率，典型的有限元网格图4所示，其中埋置率d=D。为了较好的确定破坏信息，桩靴附近使用较为密集的网格。因为不排水的土壤具有不可压缩性，所以在模型中采用一阶、完全集成、混合连续的C3D8H网格来模拟土壤。有限元加载和控制边界条件时，混合使用位移和压力变量(而不是仅仅位移)去近似模拟平衡方程和协调条件。把桩靴表示为离王立忠等人，中国海洋工程，第25卷，第一期，15-304散的刚体，荷载参考点的运动定义为桩靴的运动。在联合载荷作用下，参考点处的垂直、水平和弯矩载荷分别表示为V、H、M，位移响应分别表示为w、u、θ(图5)。2.3分析方法采用数值分析中广泛采用的小位移有限元方法。为实现单一载荷，利用位移控制加载的方法，通过垂直、水平位移和旋转的方式将载荷作用到载荷参考点上，直到达到破坏载荷。对于联合加载下的分析已经在二维加载情况V-H、V-M、H-M和三维加载V-H-M中展开。为了确定桩靴基础在联合载荷作用下的失效模式，在极限状态下应当按照一定的程序施加载荷并且获得加载区域的失效点。当桩靴基础即将达到极限状态并且将要失效时，可以确定破坏包络线上的一个点（V、H、M）。冲击实验方法和恒位移比探测可用于确定在V-H-M负荷空间中的失效包络线。冲击试验采用Tan在离心试验中的加载方式。Tan(1990)，Martin(1994)和Martin、Houlsby(2000)都认为，在实验期中加载路径非常接近破坏包络线。或者说(Gottardi等人，1999)，弹性刚度与塑料刚度的比值必须很大，否则需要进行适当的纠正。该方法已成功应用在很多研究中(Martin，1994；Bransby和Randolph，1998、1999；Martin和Houlsby，2000；Byrne等人，2002；Cassidy等人，2004；Taiebat和Carter，2000；Gourvenec，2008；Yun和Bransby，2007)。撞击试验包含二步。第一，对于i：j加载区域的破坏包络线，施加在i坐标方向的位移由零一直到极限加载情况；第二，当在i（Ui）方向的位移增量一直保持为零时，增加j(Uj)方向的位移，直到j(Uj)方向的极限载况。冲击载荷的一个例子如下图6所示。王立忠等人，中国海洋工程，第25卷，第一期，15-305在小张力有限元分析中，一旦达到极限荷载Fiult，Fi不会随着同方向位移的增加继续增加。因此，塑性刚度在破坏时为零，随着冲击试验的进行破坏轨迹不会继续扩散。然而，由于弹塑性屈服发生在破坏的轨迹线内，在实际破坏包络线内侧冲击试验会留下加载的轨迹(Gourvenec和Randolph，2003)。为了进一步检验由冲击试验获得的破坏包络线的准确性，本文后面采用了Bransby和Randolph于1997年提出的探测试验。加载路径从原点开始，变化梯度是由土壤起初的弹性刚度确定；随着塑性变形的发生，梯度将会逐渐的变化，直到位于破坏包络线上，加载路径沿着破坏包络线的方向直到法向与采用的位移比例一致。因为冲击试验会留下一个由内侧靠近破坏包络线的加载路径，论文中采用冲击试验结合探测试验的方法来确定破坏包络线。2.4SignConvention符号规则论文中的力和载荷遵从右手定则和顺时针正对流原则，参考1997年Butterfield等人的论文。力与载荷的符号如下表1所示。纯载荷作用下的极限载荷由下标“ult”表示。王立忠等人，中国海洋工程，第25卷，第一期，15-306第三章结果和分析3.1纯载荷作用下的极限承载能力图7-12表示分别在垂向载荷、水平载荷、弯矩单独作用下的破坏机理。表2总结了有限元网格在垂向、水平和弯矩的极限承载能力与Skempton公示推导的结果数据。桩靴在垂直、水平或弯矩载荷作用下的破坏，可能位于很浅，也可能位于很深的位置，这取决于埋深（图7-12）。从图形和表格2中可以看出，对于匀质粘土，极限的埋置深度大概为1.0D。当埋置深度小于临界埋置深度，失效位置很浅，可以到达土壤表面。否则，王立忠等人，中国海洋工程，第25卷，第一期，15-307破坏位置会较深，位于桩靴附近的局部区域，并且极限承载能力将随着埋深增加。在非均质粘土中，对于垂直和弯矩载荷作用，临界埋深为1.5D，但是在水平载荷作用下没有明显的临界埋深存在。王立忠等人，中国海洋工程，第25卷，第一期，15-308Skempton于1951年提出的经验公式已经广泛应用于预测桩靴在粘土中的极限承载能力。qu=6SU(1+0.2dD)+γ′≤9SU+γ′d(1)其中d是埋深；D是直径；SU是桩靴以下D/2处的平均粘土剪切强度；γ'是有效的单位重量。公示1中超出的承载力部分γ'd仅仅在桩靴上部空洞是开时才是适用的。如果假定土壤回流是充分的，那么方程可以修改为：qu=6SU(1+0.2dD)+γ′VSA≤9SU+γ′VSA(2)其中vs是桩靴的体积，表2表明，埋深高达2D的桩靴，由Skempton经验公式预测的极限承载能力要比有限元方法得到的小40%，这主要是由于不同的地基形状经验公式是对于长条状地基提出的，并且修改以适用于圆形的基础，但是实际上桩靴是纺锤形的。此外，关于埋置桩靴基础在水平力和弯矩作用下的极限承载能力的文章很少，桩靴基础受到垂直、水平和弯矩载荷作用下的有限元计算结果可以应用到实际的工程中。王立忠等人，中国海洋工程，第25卷，第一期，15-3093.2V-H承载平面的破坏包络线图13、14表示在充分考虑每一种埋置率和土壤强度分布的情况下，垂直和水平载荷联合作用下的桩靴极限承载能力。图13a、14a分别表示在均质和非均质粘土中的桩靴受到V-H联合载荷作用下，冲击试验和探测试验所得结果的对比。如图所示，不同的加载方法的载荷路径是不同的，他们在不同点具有变化。联合作用下的破坏包络线由这些点的连线决定。图像也表示出冲击试验得到的破坏包络线比探测方法要小，这也与Gourvenec和Randolph（2003）所得的结论一致。图13b、14b分别表示出匀质