MODELING OF ROCK DISCONTINUITY WITH MESHLESS METHOD
无网格法对岩体不连续面的模拟

鉴于以Goodman单元为代表的界面单元在传统有限元法中取得了广泛而成功的应用,因而采用无网格法进行岩土工程数值分析时,首先考虑引入Goodman单元以模拟不连续面.详细分析将Goodman单元引入无网格法模拟不连续面的可行性和当前研究所存在的问题.研究结果表明,由于Goodman单元最初的应用对象是有限元法,其位移模式实质是为了实现与其相连接的有限单元相协调;而无网格的位移模式是基于离散节点的,对于不同的待插值点,影响域覆盖该插值点的节点及其数量也不相同;因而必然存在着无网格法的位移模式与Goodman单元边界上假定的位移模式不相协调的问题.解决这个问题的关键是在考虑节理单元的刚度对总体刚度矩阵的贡献时,应采用数值积分计算得到节理单元的刚度矩阵而不是简单的将传统解析表达式累加到求解系统的总体刚度矩阵中去;考虑到上述特点,界面单元可以隐式或显式出现在计算模型中.以自然单元法为例详细介绍显式和隐式Goodman单元的具体实现方案,并给出相应的算例进行数值验证和对比;所提出的思路对于一般的无网格方法都是适用的.