静压桩沉桩挤土效应的有限元分析.docxVIP

静压桩沉桩挤土效应的有限元分析

0静压桩沉桩分析的应用

静压桩为压缩桩。虽然较之打入桩,静压桩的贯入过程要平稳得多,但它的贯入将使周围土体发生侧向移动,地表发生隆起,产生挤土效应。静压桩沉桩挤土效应的研究已成为各国学者和工程界广泛关注的焦点。

圆孔扩张理论被广泛应用于静压桩沉桩挤土的分析。该方法最初被使用于金属加工方面,后被推广用于岩土工程中的压力测试分析、深基础的承载力分析以及由于沉桩而产生的土体扰动分析。其基本思想就是将静压桩的压桩过程看作是在半无限体中扩张出一个与桩同直径R的圆形小孔,扩孔完成后,在孔的周围将形成塑性区。然而其求解过程需要假设沉桩问题是轴对称或球对称的,严格限定各变量场如位移场、应变场、孔压场等仅仅依赖于径向坐标,并且依赖于土体的弹塑性模型。但是可以从该理论的思想得到启发,用有限元方法来实现这一过程,并且由于有限元方法的广泛适用性,可以推广于分析方桩扩孔问题的研究,使问题的研究不再局限于圆孔扩张问题,扩大了应用范围。

1使用fem模拟需要解决的问题

在进行FEM模拟的时候,需要解决以下2个问题。

1.1初始孔径a0的影响

实际沉桩情况下,小孔的初始半径为零,但用有限元分析小孔扩张时,小孔初始半径R0必须为非零值,以避免在r=0处产生无穷大的应变。根据文献的研究,这种限制不会造成结果的矛盾。解析解表明,在同样的物质中,一小孔从r=a0开始扩张,在经历了较大的变形后才达到极限扩张压力,有限元最后计算值和解析解相差在6%以内。因此用初始半径a0的小孔扩张到2a0的过程来模拟0～R的小孔扩张半径为R桩的沉入过程,在精度上是能接受的。

R与a0的关系可以从体积关系推出,即排开的土体积相等,那么π[(2a0)2-a2002]=πR2,所以a0=R/3√a0=R/3.对于截面为a×a的方形桩的扩张问题,同样可以转为化a0×a0的小孔扩张为2a0×2a0的扩孔问题,且由(2a0)2-a2002=a2,得到a0=a/3√.a0=a/3.

1.2沉桩过程分析

沉桩时,桩周土体发生较大的位移,即涉及大变形问题,桩周土体处于大变形、大应变状态。研究表明,大变形结果和小变形结果差别不大,约6.4%。因此,采用小变形来分析沉桩过程是可以接受的。文中分析计算时采用小位移、小应变。为了反映大变形,在每一个增量计算结束后,将有限元网格的坐标重新计算,重新生成劲度矩阵,代入下一次增量计算。这种方法近似于修正的Lagrangian法,但是应变的表示式仍然是一次的,不同于大变形的二次应变计算式。

2为沉桩和膨胀孔辩护的fem模拟

2.1应力应力应变张量和应力偏量的计算

在模拟过程中,采用弹塑性模型描述桩周土体。文中采用熟知的Druker-Prager准则进行分析。其函数形式为f=αI1+J2??√?k=0.(1)f=αΙ1+J2-k=0.(1)

式中α,k为材料常数;I1,J2分别为应力张量第一不变量和应力偏量第二不变量。对于平面应变问题

α=sin?3(3+sin2?)√,k=3√cos??c3+sin2?√.(2)α=sin?3(3+sin2?),k=3cos??c3+sin2?.(2)

式中c,φ分别为材料的粘聚力和内摩擦角。

弹塑性本构关系的增量形式为

d{σ}=[Dep]d(ε).(3)

式中d{σ}为应力增量列阵;[Dep]为弹塑性矩阵;d{ε}为应变增量列阵。

在进行模拟时采用增量形式加载,基本方程为

[Kˉˉˉˉ]{Δu}i={Δf}i+{Δfp}i.(4)[Κˉ]{Δu}i={Δf}i+{Δfp}i.(4)

式中{Δf}i为初始边界荷载或第i次开挖边界释放等效节点力列阵;{Δfp}i为由塑性应变引起的附加等效节点列阵。

2.2小方孔扩孔模拟

陕西渭北某煤矿于1999年进行了静压桩地基加固工程。文中对76号桩到达地面下10m时的扩孔进行模拟。在该项目中,采用的是微型静压桩,桩体截面为20cm×20cm,扩孔模拟时将11.55cm×11.55cm的等效小方孔扩张为23.1cm×23.1cm.初始小孔采用空单元进行模拟。该土层的计算参数见表1。计算区域取为22倍的桩体边长,即4400cm×4400cm的方形区域。计算网格如图1所示,共划分了3600个单元。计算时,左右两侧采用限制水平位移的滑动支承约束,上下2侧采用限制竖向位移的滑动支承约束。

2.3fem模拟的结果和分析

2.3.1fem模拟输出结果

由图2可以大致看出,桩体周围的土体变形主要发生在桩体周围,距离桩体较远处变形比较微弱。

如图3所示,根据FEM的模拟输出结果分别做出了x和y方向的位移等值线。其形状是完全相同的,只是方向不同。由等值线可以看出,土体位移主要发生在桩体5倍半径(此处为方桩,半径认为是桩体边长的一半