层状地基下高强混凝土管桩压桩过程数值分析.docxVIP

层状地基下高强混凝土管桩压桩过程数值分析

静压预制混凝土（phc）管桩具有噪声小、无泥浆污染、施工周期长、桩重、单位承担费低等技术优势，在软土地利用工程中得到广泛应用。但是,压桩过程中由于挤土效应导致的诸多问题,尤其是对周边环境的不利影响,始终限制着该技术的进一步推广应用。

静力压桩过程是桩端土在压桩力作用下连续发生塑性破坏,产生压密和侧向挤出,并进一步在桩侧摩擦力作用下产生剪切变形的过程。静力压桩过程的力学机理十分复杂,力学行为分析涉及到的问题有大变形、土体本构关系、桩土接触和固结等。目前基于圆孔扩张理论的解析解尚难以全面考虑上述诸多问题。

本文在已有有限单元法模拟压桩过程的基础上[6,7,8,9,10,11,12,13,14,15],选用某实际工程的层状地基条件,考虑桩在穿越不同土层时桩土接触界面性质的变化,进行层状地基静力压桩过程的数值模拟分析。

1数值模拟的重要力学问题

1.1桩土接触界面性质的模拟

静力压桩过程是一个桩土相互作用的动态连续过程,其本质是桩土在接触界面互相挤压和相对滑移的过程。为了模拟层状地基条件下的整个压桩过程,需要预先在桩土所有可能接触的位置设置接触面,并考虑桩在穿越不同土层时桩土接触界面性质的变化。接触面算法有拉格朗日乘子法和罚函数法,对于存在明显相对滑动的桩土界面应采用拉格朗日乘子法。

1.2确定变力的模拟

压桩过程是在压桩力作用下完成的,由于压桩过程中压桩力是个不确定的变力,直接施加力荷载不能有效地运用有限元程序模拟整个压桩过程。由于实际压桩过程是通过控制桩体下沉速率调整压桩力的,因此采用位移荷载法,即在桩顶施加定位移,模拟压桩过程。

1.3小变形假定适用条件

静力压桩过程中桩周土体应力、应变很大,在拉应力和剪应力作用下进入塑性破坏状态,发生大位移和大转动。应变与位移梯度不再符合线性关系,因此小变形假定不再适用。应采用连续介质力学大变形理论中的物质描述方法,以质点位移为基本位移场变量,实现对土体介质质点运动过程的跟踪模拟。

1.4土的性质与土种

层状地基静力压桩过程数值模拟选用的土体本构关系既要合理模拟桩周土塑性区的应力、位移特征,还要能反映不同土层强度性质的差异。压桩过程中,桩周土的性质发生变化,结构性较强的土呈现应变软化性状,而密实度不高的砂性土则呈现应变硬化性状。以往文献中有采用边界面模型和HiSS模型的,而更多的是采用满足Mohr-Coulomb屈服方程的理想弹塑性模型和能够较好反映软黏土应力—应变关系的修正剑桥模型。

2有限模型

2.1管桩的入土深度

本文选取直径为40cm,壁厚为8cm的PHC管桩,其入土深度为36m,容重22.0kN·m-3,弹性模量25GPa,泊松比0.17。层状地基土层的分布及相关物理、力学性质指标见表1。

2.2管桩桩长丝网应力分析的准则

(1)各层土体均为均质连续各向同性理想弹塑性体,满足Mohr-Coulomb屈服准则,采用总应力法进行计算分析。

(2)根据压桩结束后PHC管桩中土塞实际高度,将压桩过程简化为初始孔径10cm的柱形孔扩张至40cm。

(3)忽略土的结构性和扰动后土体力学性质的变化。

(4)桩土接触后不再分离,即桩土接触后处于接触并滑动状态。

2.3界面接触参数

有限元模型采用二维轴对称模型,连续压桩过程通过施加桩顶位移模拟。桩体材料采用线弹性本构关系。考虑几何非线性,桩周土体变形采用大变形理论描述。

在桩土所有可能接触处设置接触属性,桩土界面采用面—面摩擦接触单元,接触面的摩擦类型为库仑摩擦,接触算法采用拉格朗日乘子法,法向接触设置为硬接触。不同桩土界面摩擦接触计算参数的选取参考勘察的桩侧极限摩阻力、静止侧压力系数及压桩力变化曲线综合确定,结果见表2。

模型的地基体深度取60m,宽度取55m,上表面为自由表面,下表面和右侧面为固定边界。计算采用通用有限元分析软件ABAQUS进行。

3结果分析

3.1水平位移沿深度方向的变化

桩周土体水平位移量是挤土效应强弱的直接反映,也是静力压桩施工影响场地周边环境的主要因素。有限单元法计算得到的桩周弹性区径向位移衰减曲线如图1所示。

由图1可见,层状地基条件下,不同深度处弹性区水平位移沿径向的衰减规律一致。由于浅层地基上覆土重较小,压桩过程中桩周土体以隆起为主,水平位移相应减小,故1层(深度2m)曲线与其他深层土体水平位移曲线差异明显。

图2为距桩中心不同径向距离处水平位移沿深度的变化曲线。由图2可知,距桩中心相同距离处,不同深度、不同土层的径向位移比较接近,但近地表处水平位移均明显偏小;距桩中心近处水平位移较大,并随径向距离增加而迅速衰减。

3.2桩周压力沿深度增长的速率

数值分析得到的不同土层桩周压力沿深度的变化曲线如图3所示。比较图中曲线可知,压桩过程