深基坑复合土钉墙支护三维数值模拟分析.docxVIP

深基坑复合土钉墙支护三维数值模拟分析 0 基于数值模拟方法的基坑支护问题 复合土钉的保护是一个三维问题，现有的理论和方法大多是将复合土钉简化为平面问题。在采用平面模型时，土钉在垂直于计算平面方向被展开成板状，从而明显夸大了钉–土界面的实际面积，而且这样展开的土钉将其上下土层完全分开，在这种情况下无论是否加入界面单元，也无论界面单元的力学参数如何取值，在理论上都不能符合实际。 同时，基坑的变形或应力应变状态是一个支护结构与土体的共同作用问题，在基坑开挖的不同工况、支护结构与土体是一个相互作用的渐变过程。 因此，一般采用数值模拟方法来研究基坑支护问题。在目前众多的数值模拟分析程序中，基于有限元差分法的FLAC-3D程序不仅能很好的反映基坑的三维问题，考虑支护结构与土体的相互作用问题，而且FLAC-3D程序能够很好的模拟地质材料受力作用下屈服、塑性流动、软化直至大变形等力学行为，因此FLAC-3D程序模拟的结果更加符合工程实际情况。 有鉴于此,本文采用美国Itasca公司开发的FLAC(Fast Lagrangian Analysis of Continua)分析软件,对深圳假日广场基坑工程中的预应力锚索复合土钉支护结构进行开挖支护施工过程的三维动态模拟分析,并将模拟模型与现场原位测试数据的拟合及对比分析,为深基坑复合土钉墙支护技术的设计和施工提供有益的指导,同时对深入了解复合土钉墙支护作用的机理具有较大意义。 1 材料本构模型或后整理后的三维力学行为 FLAC-3D是美国Itasca Consulting Group Inc开发的，属于三维快速拉格朗日分析程序。三维快速拉格朗日分析是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。FLAC-3D包含了10种弹塑性材料本构模型，有静力、动力、蠕变、渗流、温度5种计算模式，各种模式间可以互相耦合，并且FLAC-3D还可以模拟多种结构形式，如岩体、土体或其他材料实体如梁、桩、壳、衬砌、锚索以及土工织物等等。 同时，FLAC-3D程序具有强大的前后处理功能，该程序以时间步长的形式来推动程序的计算，可以记录每一时间步长下结构的受力状态，从而可以对结构从刚刚受力直至结构发生屈服、破坏的整个过程进行分析，并可以形成一个flac.mov动画文件，通过动画的方式显示整个过程，使得FLAC-3D成为一个求解三维岩土问题的最理想工具之一。 2 三维数值模拟分析模型的构建 2.1 喷射混凝土社会面 土钉、锚索采用FLAC-3D中的索结构(cable structure)。索结构仅能够承受拉应力和压应力，不能承受弯矩的作用。Buhan根据屈服设计理论认为，除了很大直径的土钉外，土钉由于抗弯作用而产生的剪应力作用是较小的，故采用索结构模拟土钉能够符合工程实际。 另外，采用索结构的最大优点在于它可以模拟土钉、锚索与周围土体或水泥注浆体的相互作用，能够准确反映土钉、锚索的轴向拉应力和剪应力的大小和分布形式。 喷射混凝土面层采用FLAC-3D中的壳结构(shell-type structure)。壳结构由多个具有三节点18个自由度的三角形单元组成(图1)。用户可以根据计算的精确程度定义三角形单元的大小，还可以规定壳结构单元的厚度、材料属性以及壳结构与土体之间的作用方式。通过对不同位置上的节点设置监测历史记录，可以得到作用于喷射混凝土面层的土压力变化过程。 腰梁采用FLAC-3D中的梁结构(beam structure)。梁结构由多个具有12个自由度的和截面双轴对称的节单元(segment)组成（图2），各节单元之间通过节点连接。梁截面为矩形，由梁截面的y和z方向的惯性矩来定义。 2.2 mb弹性–完全塑性模型 目前在土钉支护数值模拟分析中，采用的本构模型主要有：非线性E-v模型、Mohr-Coulomb弹性–完全塑性模型、非线性E-B模型、渐进单屈服面模型和Drucker-Prager模型。 本文土体本构模型选用Drucker-Prager模型，Drucker-Prager模型考虑了静水压力对屈服强度的影响，模型参数少且易于确定，计算简单，而且考虑了岩土材料的剪胀性。 2.3 模型尺寸和约束 Duncan和Goodman认为，边坡的计算范围可取边坡坡脚以下H深、坡脚左右水平2～3H的范围，H为边坡高度。据此，建立的模型尺寸如图3所示。 由于模型尺寸足够大，可以认为开挖对模型边界的应力和应变影响较小，故对模型的前面、后面和底面采用固定约束，在侧面仅对y方向施加约束，限制土体y方向的位移。至于模型侧面的应力边界条件，可认为近似服从静止土压力的分布形式（图4所示）。 2.4 土钉的flac-3d模型 由于建立整个基坑南侧边坡模型太大以至于网格和节点太多，且土钉、锚索数量众多，软件计算速