正常固结黏土中扭矩对负压沉箱承载力的影响分析

周松望，张艳，王栋

(1.中海油田服务股份有限公司 物探事业部工程勘察作业公司，天津 300459；2.中国海洋大学 环境科学与工程学院，山东 青岛 266100)

负压沉箱常用作固定平台和海底管汇的支撑基础[1-4]。在海底管道的日常运行阶段，受管道内热应力和海底底流的影响，负压沉箱不仅要承受竖向荷载(V)、水平荷载(H)和弯矩(M)的作用，还会受到扭矩(T)的影响，扭矩可能引起沉箱其他承载力分量的降低。Finnie等[5]给出了考虑扭矩影响的浅基础和桩基础的水平与竖向承载力计算公式；针对扭矩和水平荷载联合作用下的浅埋矩形基础，Nouri等[6]提供了塑性极限分析和三维有限元解答。对于长径比L/D(L和D分别代表沉箱的入土长度和直径)大于1、作为锚泊基础的负压沉箱，扭矩会造成竖向抗拔承载力和水平承载力的降低，但当施加的扭矩不超过20%的扭转承载能力时，扭矩的影响可以忽略[7-10]。支撑管汇的负压沉箱承受的竖向荷载为压力，而不是拉力，长径比L/D又常介于1～2之间[4]，目前，还不清楚该长径比范围内沉箱承受的扭矩对其他承载力分量(即V、H和M)的影响程度。

笔者采用有限元方法，模拟不同扭矩条件下L/D=1～2的沉箱与正常固结黏土的相互作用，探索扭矩对竖向承载力、水平承载力和抗弯能力的削弱机制。在分析大量变动参数的基础上，提出考虑扭矩影响的承载力计算公式。

1 有限元模型与参数设置

采用Abaqus软件建立沉箱与土相互作用的三维有限元模型。与已有的研究[11-12]类似，取沉箱顶面中心为讨论水平荷载和弯矩的参考点。图1给出了沉箱尺寸、荷载与位移的符号和方向规定：z向向下为正，u和w分别为水平和竖向位移；转角θ对应弯矩M，代表沉箱围绕参考点的转动角度；β为对应扭矩T的沉箱扭转角度。取典型长径比L/D=1、1.5和2，沉箱壁厚为0.01D。各承载力结果将采用归一化表达，试算表明，沉箱直径D的取值不影响归一化公式，因此，除特殊声明外，均以D=10 m进行讨论。为避免边界对承载力的影响，模拟的土体范围为：径向由沉箱侧壁向外延伸3.5D；深度方向由沉箱底部向下延伸3L。土体侧边界径向位移为零，土体底部为固定边界。以L/D=2为例，有限元网格如图2所示，土体剖分采用线性六面体单元，完全积分，沉箱附近土体采用细网格，侧壁和刃角下的典型单元大小为0.01D。负压沉箱刚度远大于土体，因此，将其简化为刚体。

图1 荷载与位移的方向规定

图2 有限元网格

中国南海、墨西哥湾和东南亚海域广泛分布正常固结黏性土，其不排水抗剪强度su随深度线性增加。

su=sum+kz

(1)

其中：sum为表层土的不排水强度；k为不排水强度随深度的增长梯度；z为土体深度。有限元变动参数分析中取sum为0、5、10、20 kPa，k为0、1、1.25、1.5、2.5 kPa/m。采用总应力形式，将土体视为Tresca理想弹塑性材料。由于不排水条件下土的体积不变，泊松比理论上为0.5，为避免数值计算中不可压缩材料引起的收敛困难，取泊松比为0.49。假定很高的土体弹性模量E=10 000su，以提高计算效率，试算表明，E的取值对沉箱极限承载力影响不大。在沉箱内侧、外侧和刃角分别设置一层软弱带，以考虑沉箱贯入施工造成的土体软化[12-13]，软弱带不排水强度为αsu，α为弱化系数，参考Supachawarote[12]、Kay等[13]及DNV规范[14]，取典型值α=0.65。当软弱带厚度不超过0.02D时，对承载力的影响可以忽略[12]，此处取0.01D。沉箱和软弱带为tie连接。

2 有限元模型验证

2.1 沉箱顶盖底面粗糙度的影响

沉箱顶盖底面与土体之间的摩擦对沉箱承载力影响很小。考察两种极端情况：1)顶盖底面与土体之间为tie连接，代表顶盖完全粗糙；2)顶盖底面与土体之间设置为光滑接触，代表顶盖完全光滑。对于3组典型工况，无量纲的扭矩T与扭转角度β的关系如图3所示，其中sut=sum+kL，代表刃角深度处的土体不排水强度。由图3可见，顶盖底面的粗糙度几乎没有任何影响，水平荷载、竖向荷载和弯矩作用下的单向承载力分析也有类似结论。以下分析中取保守情况，假定顶盖底面为光滑。

图3 顶盖底面光滑程度对扭矩的影响

2.2 与以往单向竖向承载力的比较

为验证建立的有限元模型的可靠性，将竖向承载力结果与Hu等[15]、Hung等[16]的结果进行对比。他们的计算均未考虑沉箱周围土体的弱化，即弱化系数α=1.0。无量纲化的竖向承载力V0如图4所示，其中，A为沉箱底端截面面积。在L/D=0.5～2范围内，所得结果与以往成果很接近，初步表明建立的有限元模型可靠。

图4 竖向承载力结果对比

3 单向最大扭矩

当沉箱发生扭转时，与外侧壁接触的土体首先被激发；随着转动增大，沉箱内侧壁和刃角深度截面处的土体逐渐发生变形，抗扭转能力取决于二者中先达到破坏的部分。假定内侧壁附近土体能够提供的扭矩为Twi，外侧壁附近土体提供的扭矩为Two，沉箱底部截面(即刃角深度处)土体提供的扭矩为Tb。由于沉箱壁很薄，内、外侧壁面积基本相同，因而

采用极限平衡方法

1)如果Twi≥Tb此时，沉箱底部截面土体先被“扭断”，破坏面沿着沉箱外侧壁和底部截面，沉箱与其内部土体之间没有相对扭转。假定沉箱底部截面上土体的剪应力同时达到不排水强度，单向最大扭矩表示为

T0=Two+Tb=

απD2L(sum+kL/2)/2+πD3(sum+kL)/12

(4)

2)如果Twi

T0=Two+Twi=απD2L(sum+kL/2)

(5)

图5为有限元给出的归一化“扭矩-转角”曲线与极限平衡结果对比，两种方法给出的单向最大扭矩很接近，有限元结果略高，再次证明了建立的有限元模型的可靠性。当L/D=1～2、α=0.65时，Twi几乎总是大于Tb，这意味着扭转破坏基本发生在沉箱外侧壁和底部截面。

图5 有限元法与极限平衡法获得扭矩对比

4 扭矩的影响

4.1 扭矩对竖向承载力的影响

沉箱竖向承载力V0由端部承载力Vb和沿外侧壁的摩擦力Vw组成。

V0=Vw+Vb

(6)

Vb激发的土体破坏场主要在沉箱底部截面之下，与扭矩引起的破坏面基本不重合，因而，可以不考虑扭矩的影响。

Vb=πD2(sum+kL)NcV/4

(7)

式中：NcV为正常固结黏土中沉箱端部承载力系数。由变动长径比的有限元参数分析得到

NcV=9.73+0.4(L/D-1)

(8)

扭矩和竖向力的共同作用造成沉箱外侧壁土体破坏，所以，Vw的计算需要考虑扭矩影响，按照Finnie等[5]的建议，当扭矩T作用下沉箱转角为β时，沿外侧壁上造成的水平剪应力τH为

τH=βDGT/2

(9)

式中：GT为沉箱外壁扭转刚度，GT=4G/D，G为土体剪切模量。τH与T之间的关系为

(10)

通过式(9)和式(10)，可以获得对应T的转角β和τH。沿沉箱外侧壁土体所承受的竖向剪应力τV、水平剪应力τH和抗剪强度αsu的关系为

(11)

τV沿沉箱外侧壁积分，即可得到扭矩T削弱(即引起τH)后外侧壁能提供的竖向力Vw

(12)

(13)

当T较大时，式(13)根号下部分可能小于0，这意味着αsu全部用于抵抗扭转，不能够再贡献于竖向承载力，剩余的扭矩和全部竖向力将作用于沉箱底部土体，此时，式(13)不再成立。

图6 扭矩作用下的竖向荷载-位移曲线

为简化设计应用，定义扭矩影响系数λT考虑扭矩对竖向承载力的降低

(14)

由于实际应用中扭矩的安全系数一般大于1.5[17]，所以，只拟合T/T0<0.8范围的有限元结果。变动sum、k和L/D，考察了50个工况，取最不利情况，拟合得到

λT=1-0.07tan(1.5T/T0)

(15)

拟合公式效果如图7所示。

图7 扭矩对竖向承载力的影响

4.2 扭矩对水平承载力的影响

图8为通过有限元得到的扭矩和水平荷载同时施加时的位移云图。由图8可见，沉箱内部和下部土体发生绕某一中心的转动，同时，沉箱外侧上部楔形土体位移显著。随着T的增加，主要表现为沉箱底部附近土体的运动范围减小。图9为存在扭矩的情况下“水平力-水平位移”关系曲线，水平力很快达到定值。仍然采用4.1节中的算例，图10给出了T/T0对归一化的水平承载力的影响：当T/T0≤0.2时，扭矩对水平承载力影响极小，当T/T0介于0.2～0.8之间时，扭矩造成水平承载力最大降低20%。

图8 扭矩对水平荷载破坏模式的影响

图9 水平荷载-位移曲线

图10 扭矩对水平承载力的影响

(16)

经拟合，扭矩影响系数λT的取值基本与式(15)相同，所以，仍采用式(15)表示。

4.3 扭矩对抗弯能力的影响

(17)

式中：λT经验算仍可采用式(15)计算。

图11 扭矩对弯矩的影响

5 结 论

正常固结黏土中支撑管汇的负压沉箱长径比大多在1～2之间，采用理论分析和有限元方法探索了扭矩对沉箱竖向承载力、水平承载力和抗弯能力的影响。结果表明，扭矩会造成其他承载力分量的降低，降低程度与扭矩大小有关：对于长径比介于1～2的沉箱，在常规土体强度分布情况下(sum=0～20 kPa，k=0～2.5 kPa/m)，当扭矩不超过20%的单向最大扭矩时，扭矩对其他承载力分量的影响非常小，可以忽略；当扭矩介于20%～80%的单向最大扭矩时，其他承载力分量最大降低20%。总结大量的有限元变动参数分析结果，给出了扭矩不超过80%的单向最大扭矩时扭矩对其他荷载分量的影响系数表达式。