基于MARC软件的岩体节理裂缝数值模拟.pdfVIP

维普资讯 第31卷 第22期 山 西 建 筑 VoI．31No．22 2OO5年 11月 SHANXI ARCHITECTURE Nov． 2005 ·91 · 文章编号：1009—6825(2005 基于MARC软件的岩体节理裂缝数值模拟 彭海军 蔡汝铭 摘 要：在地下洞室的数值模拟计算分析中，研究的对象主要是岩体，研究方法很多，而有限元方法在实践中是使用最为 广泛的，简析了裂缝模拟的基本原理，使用有限元分析程序MARC模拟裂缝，分析裂缝存在时对洞室结构的各种影响。 关键词：有限元法，数值模拟，岩体裂缝，MARC 中图分类号：TU457 文献标识码：A 引言 本构方程：io-} [D][B]i}， 地下洞室～般是指由人工开挖或生产活动在地层中形成的 由虚功原理得到的平衡方程：iF} [B [D][B]以i}。 地下空间，其研究的对象主要是岩体 ，岩体是一种天然形成的复 其中，}}和{sf分别为单元的应力与应变分量；[D]为平 杂的地质介质而非性质单一的材料。而含有软弱夹层及节理裂 面应变弹性矩阵；[B]为位移一应变关系矩阵；t为单元厚度；A 缝的岩体又极为常见，因而用有限元分析岩体时模拟裂缝方面显 为单元面积。 得十分重要，完全模拟这种地质介质难度相当大，在实际的数值 则单元劲度矩阵[忌e=[B]T[D][B]以。根据单元劲度矩 模拟过程中，都要进行适当的合理的简化，以方便模拟。文中用 阵总装整体劲度矩阵得出结构的平衡方程[K]{}=[R]。 MARC软件对数值模拟岩体裂缝进行了简化，进而得到合适的计 2 节理裂缝的模拟简化 算模型，在此基础上进行计算得出相应的结论。 岩体是一种地质介质，在它形成的过程中，大多经历了许多 1 将空间问题简化为平面应变问题 的地质构造活动，这些地质构造运动在地层中形成了一系列的构 由于地下洞室相对于洞室断面来讲往往是很长的，可以认为 造形迹，大型的有断层、褶曲等，小型的有节理、断裂、裂缝等，而 变形分量和位移分量都不会沿 z轴方向(洞室长度方向)发生变 小型的节理裂缝的宽度几乎接近于零，这时如仍沿用普通的平面 化 ，所以各点都只会有x和Y方向的位移而不会有z方向的位移。 单元模拟就会出现困难，目前普遍采用的方法是用 goodman单元 因此可以将空间的数值模拟问题简化为平面应变问题 J。如果 模拟节理裂缝。按平面应力问题考虑，其基本原理如下 J： 采用三角形单元划分所需模拟截面，则有限元的计算三大方程可 假设一种 “零宽度矩形节理单元”，如图1所示节理单元图。 以简单用节点位移 参数表达出来。 底面任一点的位移为：{}B=[N]{U}B， 几何方程：{￡} [B]{}， 顶面任一点的位移为：{}．r=[N]{U}．r， 基坑深度和长度的加大，将大大地增加现有结构的维护难 临界流速的任何地质条件下使用。6)对于软土深基坑 ，冻胀和冻 度。但对冻结法而言，冻结深度的增加 ，施工的难度并不增大，却 融沉陷不利于冻结法在基坑工程中应用，但可采取文中所述措施 相对节省大量的辅助工作量，总造价的增加比例也较小。由此可 抑制冻胀，消除下沉。 见 ，基坑越深，冻结法施工越具有优越性。 参考文献： 3 结语 [1]李永盛、深基坑维护结构的变形、稳定与相邻构筑物的影响 [J]．建筑施工，1994，16(4)：29—38． 1)淤泥质粘土冻结