王金梅, 张迎宾, 赵兴权,2, 余鹏程, 王 潘

(1. 西南交通大学土木工程学院, 四川 成都 610031; 2. 山东建筑大学土木工程学院, 山东 济南 250101; 3. 四川省地质工程勘察院, 四川 成都 610072)

DDA (discontinuous deformation analysis)方法[1]是以研究非连续块体系统不连续位移和变形为目的的一种数值方法,能够很好地模拟块体间的滑动、张开和闭合,自提出以来,已广泛地应用于滑坡[2-6]、落石[7-9]、隧洞坍塌[10-11]、开挖爆破[12-13]和倾倒破坏[14]等许多工程领域.

许多学者对DDA进行了更加深入的研究.作为一种非连续介质力学数值分析方法,DDA方法的最大特点是可以反映岩石块体之间接触面的滑移、分离与倾覆旋转等大变形大位移问题,因而,块体间的接触相关问题成为DDA方法的核心问题之一.针对DDA方法在接触上存在的一些缺陷,许多学者做了进一步的修正和改进.焦玉勇等[15]在DDA方法中加入块体自适应剖分算法,并引入虚实节理转换来模拟裂纹的扩展、贯通到碎裂、坍塌的全过程.夏才初等[16]利用表征虚节理向实节理转化程度的虚节理连通率对虚节理力学性质参数的弱化规律进行研究,基于Jennings强度准则和最大抗拉强度准则导出了节理力学性质参数弱化函数,以模拟断续节理的开裂.初始的DDA程序由于其内在的接触判断、接触处理方式导致材料粘聚力不能准确地被考虑,使DDA的稳定性分析功能被大打折扣,对于此问题,Wang等[17]采用接触界面相对位移量作为界面粘聚力的失效判定,提出了与接触界面相对位移相关的剪切强度……