基于FLAC3D二次开发的不连续变形模拟方法探讨.pdfVIP

198 岩土工程学报 2013年 小于黏结强度时，接触面单元为弹性：当剪应力大于 1 FLAC3D不连续变形分析的基础 剪切强度或法向拉应力(有效应力)大于抗拉强度时， 接触面问的连接断开。同时，当接触面在连接未断开 FLAC3D软件中提供了interface接触面单元，该 时，不产生分离变形，仅产生剪切滑移变形。 单元可用于模拟节理、断层、桩土接触面等物体间的 (1)弹性阶段的控制方程： 接触、滑移、张开等问题¨6‘17J。该接触面单元具有摩 擦角、黏聚力、膨胀角、法向刚度、切向刚度以及拉、 碍“∽=knu。A+o-A，(1) ， (2) 剪黏结强度等参数，能很好地模拟裂纹面间的滑移、 F。s，”山’=F。s，D+蜒△材：，o佗皿”A+gf彳 张开直至完全断开的过程。 刻的剪切力。 接触节点 接触单元 (2)库仑滑移阶段的控制方程 。 (3) 只。。=cA+tan《b(F．一pA) 从以上描述中，可以看出FLAC3D中的接触面单 ‘衬N． 元可以用于模拟实际工程中裂纹的张开扩展过程。但 是由于FLAC3D中接触面的建立较为复杂，且其动态 ‘心 的二次开发。因此，实现FLAC3D不连续变形分析模 魁． 节点代表区域 二次开发。 图1接触单元模型 Interfacemodel Fig．1 2不连续变形模拟程序的二次开发 1．1接触面单元的组成成分 2．1 FLAC3D不连续变形模拟思路 接触面由一系列三角形单元组合形成，其中每一 莫尔库仑模型可用于模拟岩土体的受力变形过 个三角形单元由3个节点组成。一般来说，接触面单 程，利用其塑性屈服准则可给出塑性区岩土体的分布 元与实体单元的表面黏结在一起；接触面的节点自动 范围，但是该模型不能考虑进入塑性区岩土体的强度 在每一个接触面单元的顶点创建而成。当其他网格的 折减和损伤效应。本文提出的FLAC3D不连续变形模 表面和接触面单元发生接触时，接触面的节点就会发 拟方法中，充分考虑了塑性区岩土体的强度折减及损 现这种接触，并且由法向、切向刚度及滑动的相关性 伤效应。 质所定义。每一个接触面