断层失稳滑动极易引发断层型冲击地压和矿震等灾害。鉴于利用引入库仑摩擦定律的连续-非连续方法不能很好地模拟断层黏滑过程,在自主开发的拉格朗日元与离散元耦合的连续-非连续方法的基础上,引入了速率-状态依赖摩擦定律。并通过模拟直接剪切条件下滑动-保持-滑动实验,验证了该方法的正确性。开展了双轴压缩条件下平直走滑断层黏滑过程的数值模拟研究。在一次黏滑中,岩样上端面平均差应力的绝对值-时步数目曲线可被分为近似线性阶段、偏离线性阶段和软化阶段,其中,软化阶段又可被分为亚失稳Ⅰ阶段、亚失稳Ⅱ阶段和失稳滑动阶段。在平均差应力的绝对值-时步数目曲线的偏离线性阶段,位于断层上部的监测单元的剪应力最先开始释放。在亚失稳Ⅰ阶段,所有监测单元均完成了由应力累积到应力释放的转化;各测点的位移值有协同化的趋势;在亚失稳Ⅱ阶段,所有监测单元均完成了由应力缓慢释放到应力加速释放的转化;各测点的位移值协同化程度加强。开展了不同倾角时平直走滑断层黏滑过程的数值模拟研究,考察了亚失稳阶段摩擦系数的变化规律。随着断层倾角的增加,应力降由规则变得不规则,初次黏滑出现时的时步数目减小,平均应力降减小,黏滑周期减小,断层容易发生黏滑,每次黏滑释放的能量越小。在黏滑前近似线线增长阶段,断层面上的剪应力和摩擦系数分布逐渐趋于均匀,在亚失稳阶段,断层面上节点摩擦系数的平均值快速下降,这可作为断层进入亚失稳阶段的特征之一。以济三煤矿6303工作面地质条件为背景建立了数值模型,开展了采动诱发断层滑动过程数值模拟研究,并讨论了基于速率-状态依赖摩擦定律和基于库仑摩擦定律得到的位移值不同的原因。在位移值的缓慢增长阶段,基于两种摩擦定律得到的位移值相差不大;在快速增长阶段,基于速率-状态依赖摩擦定律得到的平均位移值较小,且平均位移值有2次突增（当工作面距断层40m时,断层共出现第1次失稳滑动;当工作面距断层10m时,断层共出现第2次失稳滑动）。这说明,利用速率-状态依赖摩擦定律可以较好地模拟出工作面推进过程中,断层突然的失稳滑动。当工作面距断层40m时,应加强顶板控制,以预防冲击地压等灾害的发生。在采动条件下,基于速率-状态依赖摩擦定律得到的位移值小于基于库仑摩擦定律的,这是由于,引入速率-状态依赖摩擦定律后,断层长时间处于速率强化状态和黏着阶段。这与直剪条件下得到的结果不同,引入速率-状态依赖摩擦定律后,由于滑块处于速率弱化状态,断层面上的摩擦系数减小,导致基于速率-状态依赖摩擦定律得到的位移值大于库仑摩擦定律的。该论文有图37幅,表2个,参考文献89篇。