基于接触问题的桩土相互作用机理研究.docxVIP

基于接触问题的桩土相互作用机理研究 结构或结构的摩擦接触通常发生在技术结构或结构的结构或结构的侧壁上，例如，在机械工程中，设备和弹性基础设施之间的压缩接触，以及在桩基工程中，桩侧土壤界面的摩擦接触。 1 节点力矢量f 一般的非线性有限元增量平衡方程可写成: 式中:K为结构切线刚度矩阵;Δu为位移增量矢量;R为参考荷载矢量 (包括接触力矢量) ;F为内部节点力矢量。右上角标 (m, j) 表示增量迭代次数为第m增量步第j迭代步。 假定荷载增量比例系数为λ, 采用子结构的思想, 把式 (1) 中的位移增量矢量Δu分解为接触面 (潜在接触面) 位移增量矢量Δu 利用自由度凝聚的方法, 并将Δu 式中符号的具体意义见文献[5]。 从式 (3) ~ (4) 可以看出, 平衡方程中位移增量 (Δu 2 接触处理的第j次迭代 定义可能发生接触面上的法向与切向拉格朗日乘子分别为λ 法向接触力的增量增广拉格朗日乘子迭代格式可写成: 式中:j表示当前加载步内进行接触处理的第j次迭代;符号为Macauley括号, 切向摩擦滑移法条件取如下形式: 采用向前差分法, 可以将式 (8) 表示为: 其中, k 式 (11) 中:Φ 当 将式 (9) 代入式 (12) , 则得到: 于是关于切向拉格朗日乘子的相应增量迭代格式可写成: 从以上推导过程可以看出, 上述增量迭代格式充分考虑了接触面的物理约束与几何约束条件, 将其与有限元方法结合起来, 便于数值计算。 3 桩身轴力及侧摩阻力分布特性 结合实际桩基工程, 取3根试桩进行现场静载检测试验。桩的直径为0.6 m, 总长为24 m。为了解每根桩在桩顶荷载作用下的荷载传递情况及各土层的侧摩阻力和端阻力, 在桩顶、桩端及土层分界面处埋设土压力盒与钢筋应变计。基于本文建立的数值计算方法进行桩土相互作用的计算。实测土层部分计算参数见表1。 数值计算结果与试验实测结果的比较见图1和图2。试验实测结果均取3根试桩的测试结果平均值, 简称实测均值。数值计算结果取计算截面处所有节点计算结果的平均值, 以下简称为计算值。 图1为测得桩顶荷载为1 000 kN时桩身轴力随深度的变化曲线。可以看出桩身轴力沿深度方向逐渐减小, 呈“倒三角锥形”分布, 桩身轴力上部大于下部。在桩顶荷载作用下, 桩身上部0~19.6 m的范围内, 淤泥质黏土层桩身轴力变化平缓, 而在19.6~24 m范围内的粉质黏土层, 桩身轴力传递变快。数值计算结果与实测结果吻合较好, 且均能反映出桩身轴力变化的特点。 图2为桩侧摩阻力随深度变化的情况, 可以看出, 桩身上部首先发生弹性压缩, 产生土对桩向上的侧摩阻力。随着荷载沿桩身自上向下传递, 桩身侧摩阻力逐渐发挥作用, 沿桩身呈非线性分布。淤泥质黏土层桩侧摩阻力沿桩长呈上大下小倒三角形分布, 淤泥质黏土层上半段桩侧摩阻力较大, 到了下半段时桩侧摩阻力显著下降, 在桩端附近的粉质黏土层, 侧摩阻力又显著增加。 通过上述分析可以知道:轴力随深度的变化和计算结果规律一致, 轴力计算值比实测值偏大;桩侧摩阻力的计算值和实测值有一定的偏差, 但总体趋势吻合较好。 上述分析表明, 计算值和实测值有一定的偏差, 但曲线变化规律基本一致, 能够真实反映桩土相互作用时荷载的传递规律。 4 竖向荷载作用下桩身轴力分布特性 1) 桩土相互作用时, 桩身侧摩阻力自上而下逐渐发挥作用, 沿桩身呈非线性分布。 2) 在竖向荷载作用下, 桩身轴力沿深度方向呈“倒三角锥形”分布, 桩身轴力在粉质黏土层衰减速度比在淤泥质黏土层中要快。 3) 通过对桩顶沉降、轴力和侧摩阻力分布对比可以发现, 计算值和实测值有一定的差别, 但是曲线变化趋势基本一致。