文档介绍：数值分析论文
论文发表专家
v^
海上平台桶形基础沉贯过程的数值分析
摘要：探讨了桶形基础沉贯过程中负压引起的渗流场分布及其对桶壁与土体间的接触应力和摩擦阻力
的影响。计ne等[2，3](2002，
2004)利用离心机模拟了竖向动荷载作用下吸力式桶形基础的水平荷载与力矩荷载耦合作用
论文发表专家一
團■中国学木期刊洌
[4](2007)通过离心机模型实验和数值计算，研究了桶基在循环荷载作用下的抗拔能力和桶
裙径向应力的变化。国内一些单位也进行了相关的模型桶基的室内和野外试验研究[5-10]。
吸力式桶形基础能否安装就位、正常运作，需要明确沉贯阻力、负压等因素之间的关系。然而，目前国内外关于负压沉贯过程的研究多集中于渗流场的分析，并未对负压引起的桶土间相互作用(接触应力及摩擦阻力)的改变进行定性分析。本文针对粘土地基中吸力式桶形基础沉贯过程进行分析，分别按照是否施加负压两种不同工况进行数值计算。探讨了负压引起的渗流场分布及其对桶-土间接触应力和摩擦阻力的
影响。通过与求解沉贯过程中桶-土接触应力的解析公式进行对比，印证了本文计算结果的合理性，并提岀了考虑初始应力的计算沉贯过程中桶-土接触应力的简化公式。
1桶形基础沉贯阻力渗流有限元计算
由于桶形基础的下沉是借助桶内外的水压差来实现的，所以在其下沉过程中必然引起桶外土体中的孔隙水向桶内运动，从而在桶周围形成渗流场。该渗流场符合达西定律的三维轴对称稳定渗流，水头分布满足三维轴对称拉普拉斯方程。
取桶体及周围土体作研究对象，沿周向取一个单位弧度，单元的划分采用三维轴对称单元。通过分批删除桶壁单元，可以模拟桶的动态下沉过程，通过改变桶内泥面上的水头边界条件可以模拟不同的负压条件,而在一定的负压值下，假定渗流场是一个稳态的过程。由于桶形基础负压沉贯时伴生的渗流场的作用,
使其周围土体中的有效应力发生了变化。在桶外土体中，由于渗流方向是自上而下的，故引起土体中的有效应力增加。在桶内土体中，由于渗流方向是自下而上的，故引起土体中的有效应力减小。根据桶内、外土体和端部土体由于渗流的作用使有效应力变化，可求出桶体内外侧阻力和桶端阻力计算公式以及桶基下沉总阻力计算公式。桶基外壁侧摩阻力为：
刀f=nDktg3丫'z+yZdZ
桶基内壁侧摩阻力为：
刀f=nDktgSy-YZdZ(2)
桶端阻力为：刀R=D-Dt，L+yLN3)
桶基下沉总阻力为：刀f=刀f+刀f+E(4)
式中丫为水的容重；丫为土的浮容重；L为土体的厚度；t为筒壁的厚度；H0为筒外泥面处的水头值；h'为筒内的真空度；Z为计算点至泥面的深度；D为土与桶基的摩擦角；D为桶基外径；D为桶基内径；h为桶壁某结点在桶外土体中的水头值；h为桶内土体中的水
头值。
2数值计算模型
对于单个桶体结构，考虑到结构与地基耦合体系几何形状与加载条件的对称性，建立计
算模型时取桶体及其周围土体的一半结构作为研究对象，如图1所示。桶体的基本尺寸为：
直径4m，，壁厚2cm。为了便于有限元网格的剖分，地基区域取为圆柱形，水平