LNG储罐桩土非线性有限元模拟.pdfVIP

LNG储罐桩土非线性有限元模拟

赵泽钦;陈团海;张超

【摘要】在对LNG储罐进行全罐结构分析和安全校核时,通常通过在承台底部或

桩身施加等效非线性弹簧等方法来简化模拟土体对上部结构的承载作用,但这种简

化的近似方法无法准确地反映桩土之间复杂的相互作用关系.采用三维非线性有限

元软件ANSYS对LNG大型储罐的桩基进行桩土作用分析,分析时通过三种方法估

算了土体的变形模量,同时考虑了初始地应力和土体宽度的边界效应,并把计算得到

的结果与现场试桩实测结果进行比较.研究表明,三维非线性有限元模拟结果与现场

实测的荷载-沉降数据拟合度较好;三种方法估算土体的变形模量,方法一(基于多孔

介质理论,由压缩模量推算变形模量)的偏差相对最大,方法二(通过标贯击数估算变

形模量)的准确性依赖于所选取经验公式的适用性,方法三(通过参数反分析法拟合变

形模量)在缺乏原位测试数据时能得到较好的拟合结果;桩侧摩阻力总体随着深度的

增加而增大,桩端处附近达到极限值,桩侧土体的总侧摩阻力占总荷载的69.7％,桩端

土体的端阻力占总荷载的30.3％.

【期刊名称】《石油工程建设》

【年(卷),期】2017(043)003

【总页数】7页(P15-20,27)

【关键词】桩土作用;有限元模拟;非线性接触;变形模量;荷载-沉降曲线

【作者】赵泽钦;陈团海;张超

【作者单位】中海石油气电集团有限责任公司,北京100028;中海石油气电集团有

限责任公司,北京100028;中海石油气电集团有限责任公司,北京100028

【正文语种】中文

LNG大型储罐的基本结构主要包括桩基、承台、预应力混凝土外罐、X7Ni9钢内

罐、吊顶及顶梁框架等，上部罐体和罐内液体的重力荷载通过承台和桩基传递到下

方的土体上，因而桩土在各种工况条件下的荷载传递和承载能力是LNG大型储罐

设计和建造过程中主要考虑的问题之一。在对LNG储罐进行全罐结构分析和安全

校核时，通常通过在承台底部或桩身施加等效非线性弹簧等方法来简化模拟土体对

上部结构的承载作用。这种简化的近似方法无法准确地反映土体复杂的本构关系以

及桩土之间的非线性作用关系。文章分析了桩土接触作用机理及适用的本构模型，

使用国际通用有限元软件ANSYS建立了桩土非线性接触的三维实体模型，通过三

种方法估算土体的变形模量，并与现场试桩实测结果进行比较和分析，为LNG大

型储罐全罐应力分析中的桩基分析部分提供了一种新的思路和方法。

1.1桩土作用研究方法

桩土交互作用是一种高度非线性接触问题，可以使用有限元法进行数值求解，分析

时需要选择桩和土体适用的本构模型和屈服准则。对于桩体，一般软件中都有成熟

的钢筋混凝土模型。土体的本构模型则比较复杂：一种为非线性弹性本构模型，如

邓肯-张模型[1]、K-G模型；一种为弹塑性本构模型，主要为Dracker-Prager模

型[2]（以下简称为DP模型）。LNG储罐桩基所处的土体符合弹塑性本构特征，

使用DP模型进行模拟。

1.2DP模型

ANSYS软件中目前能用于近似描述岩土材料的模型只有DP模型。DP模型是理

想弹塑性模型，是对Mohr-Coulomb准则的修正。它主要包含三个参数：内聚力

c（Pa或kPa）、摩擦角准（°）和剪胀角Φf（°），其屈服准则为：

3βσm+[0.5{s}T[M]{s}]0.5-σy=0（1）

式中：β为与摩擦角准有关的无量纲材料常量；σm为平均应力，即土力学中的静

水压力，Pa或kPa；{s}为偏应力，Pa或kPa；T为矩阵转置符；[M]为Mises屈

服准则中的系数项，无量纲；σy为与内聚力c和摩擦角准相关的材料屈服参数，

Pa或kPa。

1.3桩土荷载传递

桩顶在受到向下压的荷载时，桩身压缩产生向下的位移，桩土之间发生相对滑移。

当施加荷载超过桩侧土体的极限侧摩阻力与桩端土体的极限端阻力之和时，桩体发

生持续的大位移沉降。

桩土接触是一个典型的刚性面-柔性面非线性接触问题，可以引入Goodman无厚

度四节点接触单元[3]。在ANSYS分析中，选用TARGE170描述桩体的刚性目标

面，CONTAC173描述桩侧土的柔性接触面，组成接触单元，接触