双层地基变形特性的有限元分析.docVIP

双层地基变形特性的有限元分析 　　摘 要：基于Mohr-Coulomb屈服准则，利用有限元软件Abaqus建立有软弱下卧层的双层体系地基模型，研究其应力和变形特性，以及不同压力作用下对地基沉降变形的影响。结果表明，由于上层土较软，地基中塑性区首先出现在交界面处的上层土中，且发展较快；随着荷载的增加，塑性区向交界面处上下土层均有扩展，最后达到地基的极限承载力。该模型为研究具有应力扩散作用的上硬下软的双层地基的应力应变及沉降变形具有一定的实用价值。 　　关键词：双层地基； 沉降变形；有限元分析； 　　中图分类号： TU433 文献标识码： A 文章编号： 　　引 言：地基土具有压缩性能，处理不合理将引起地基变形，导致建筑物上部结构产生附加应力。若地基变形量过大，将会影响建筑物的正常使用，甚至危及建筑物的安全。对于砂性土地基，其沉降主要是瞬时沉降。对于一般的粘性土地基，当荷载刚加上时，由于土骨架的弹塑性变形的结果，地基将产生较小的瞬时沉降，随后将是大量的随时间而发展主固结沉降和次固结沉降。而我国东部沿海地区，上硬下软的双层地基是极为常见的。国内外研究学者对此进行了大量的研究工作，并提出了基于弹性理论和弹塑性理论的线弹性理论范围和弹塑性理论范围内不同边界条件下的线性、非线性解析解和相关数值模型，对进一步研究具有应力扩散效应的软土硬壳层的应力应变和沉降问题起到了很好的促进作用。本文将通过工程实例，利用有限元软件ABAQUS进行弹塑性分析，研究了荷载作用下该类地基中塑性区的产生和发展规律。 　　1 地基变形算例分析 　　1.1 土的本构模型 　　Mohr-Coulomb破坏和强度准则在岩土工程中应用十分广泛。该模型是基于材料破坏时应力状态的莫尔圆提出的，破坏线是与这些莫尔圆相切的直线。Mohr-Coulomb强度准则为 　　τ=с-σtanφ（1） 　　式中，τ为剪切强度，σ为正应力，с为材料的粘聚力，φ为材料的内摩擦角。 　　其中莫尔圆中存在以下关系 　　τ=scosφ （2） 　　σ=σm+ssinφ（3） 　　把τ和σ代入公式（1），则Mohr-Coulomb准则可以写为 　　 s+σmsinφ-ccosφ=0（4） 　　式中，s=(σ1-σ3)/2；σm=(σ1+σ3)/2。 　　1.2弹塑性分析 　　双层地基的承载力将随软硬地层条件的变化发生重大变化。由于硬壳层比软卧下层土的压缩模量大，土体强度高，在一定的外荷载作用下可以承担部分压力，软卧下层土由于硬壳层的扩散作用而没有发生较大变形。在不同层土体中，土的竖向不均匀性要比各层之间的差异要小得多，因此我们可以简化计算，假定土体在每一层中各向同性。同时，土是松散的颗粒体，且一般情况下，土体在实际工程中是受压的，计算中认为土不能承受拉力。 　　结合某工程实例，对于基础作用下的上硬下软双层地基变形弹塑性分析，上层土厚度Z=0.27m，础长B=1.44 m，计算模型的宽度和高度为7.2 m。 　　 　　图1.1 地基模型图1.2 地基计算模型 　　模型中，土层材料均采用Mohr Coulomb弹塑性本构关系，上下土层的物理力学参数见表1.1所示。 　　表1.1 模型材料计算参数 　　 　　 　　由于土不承受拉力，假定各层之间完全连接，不发生分离现象。土为半无限体，模型采用二维实体单元。有限元模型中的边界条件与实际土体状态相同，底层固结，由于侧面周围土压力的挤压与磨差作用，只约束土的水平自由度。加载方式与计算模型相一致。 　　2. 有限元分析结果 　　采用Mohr-Coulomb屈服准则，利用大型有限元程序Abaqus进行建模，图2.1、2.2及2.3所示的为土层应力应变分布云图。 　　 　　 图 2.1 水平应力分布云图图 2.2 竖向应力分布云图 　　 　　 图 2.3 土层沉降变形分布云图 图2.4Pu=25kPa时塑性应变分布图 　　由应力分布云图可以看出：在沿基底中心点X向，水平应力分布在软硬土层交界面处发生显著变化，上硬壳层应力明显大于下软卧层，承担外部主要荷载；竖向应力在上硬壳层土中变化坡度较陡，而在下层土中变化坡度较为平缓。两应力云图可以看出在两层地基交界面处出现了应力扩散现象。土的沉降变形也主要是由于上层土压缩产生。采用Mohr-Coulomb屈服准则，经过有限元计算得到的塑性应变图。当基底压力接近25kPa时，地基土首先在交界面的上层土中开始出现塑性区，且随着基底压力的增大，塑性区逐渐增大，当基底压力大于56kPa时，计算结果开始不收敛。 　　图2.4、2.5分别为基底压力为25kPa和30kPa时地基的塑性区分布图。 　　 　　图2.5Pu=30kPa时塑性应变分布图图2.6Pu=50kPa时塑性应