裂纹几何配置对裂隙砂岩力学特性及破坏模式的影响_田茂祥.pdf

收稿日期2020-02-15 基金项目重庆市研究生科研创新项目 （编号 CYB19015） 。 作者简介田茂祥 （1987） ， 男， 讲师， 硕士。 总第 528 期 2020 年第 6 期 金属矿山 METAL MINE 裂纹几何配置对裂隙砂岩力学特性及 破坏模式的影响 田茂祥 1 雷瑞德 2 张亮 3 陈小平 1 王丰 11 （1. 重庆市市政设计研究院， 重庆 400020； 2. 重庆大学资源与安全学院， 重庆 400044； 3. 重庆大学土木工程学院， 重庆 400045） 摘要天然岩体常常包含不连续体， 比如孔隙、 裂隙、 断裂以及其他缺陷等， 这些缺陷的存在将导致岩体力 学强度产生较大的劣化作用。为了研究不同裂纹几何配置对裂隙砂岩力学特性、 裂纹扩展演化及断裂机制的影 响， 基于室内砂岩力学参数测试结果， 借助颗粒流平行黏结模型对预制裂隙试样进行一系列的模拟研究。结果表 明 当裂隙倾角恒定时， 随着岩桥角度的增加， 砂岩峰值应力与峰值应变呈现出先降低后增加的趋势； 当岩桥角度 恒定时， 峰值应力与峰值应变随着裂隙倾角的增加而增加。初始裂纹从预制裂隙尖端萌生， 起裂应力水平为 0.68σc； 累积微裂纹数呈指数函数形式增加； 试样破坏形式主要有拉伸破坏、 剪切破坏及拉剪混合破坏； 岩桥贯通 模式主要有间接贯通和直接贯通两类。 关键词预制裂隙裂纹扩展PFC2D平行黏结模型应变能 中图分类号TD313文献标志码A文章编号1001-1250 （2020） -06-142-08 DOI10.19614/ki.jsks.202006022 Effect of Flaw Geometric Configurations on Mechanical Properties and Failure Modes of Cracked Sandstone Tian Maoxiang1Lei Ruide2Zhang Liang3Chen Xiaoping1Wang Feng12 （1. Chongqing Municipal Research Institute of Design， Chongqing 400020， China； 2. School of Resources and Safety Engineering， Chongqing University， Chongqing 400044， China； 3. School of Civil Engineering， Chongqing University， Chongqing 400045， China） AbstractNatural rock masses often contain discontinuities，such as pores，fissures，fractures and other defects． The existence of these defects can greatly degrade the mechanical strength of rock． In order to study the effects of different crack geometries on the mechanical properties，crack propagation and fracture evolution mechanisms of fractured sandstone，based on the test results of mechanical parameters of sandstone in laboratory，a series of simulation studies on the fissured samples were carried out by means of the particle flow parallel bonding model firstly． The results show that when the flaw inclination angle is constant，the peak stress and peak strain of sandstone decrease first and then increase with the increase of ligament angle． When the ligament angle is constant，the peak stress and peak strain increase with the increase of flaw inclination an⁃ gle． The initial crack originates from the pre-existing crack tip and the cracking stress level is 0．68σc． The cumulative micro⁃ cracks of sandstone increase exponentially with the increase of axial strain． The failure mode of the specimen mainly includes tensile failure，shear failure and tensile-shear mixture failure． The coalescence mode of the ligament mainly includes direct coalescence and indirect coalescence． KeywordsPre-existing flaw， Crack propagation， PFC2D， Parallel bonded model， Strain energy Series No． 528 June 2020 岩体中常常含有大量的缺陷， 比如孔隙、 裂隙及 节理等。这些不连续体通常是由地质构造、 长期风 化及人类开挖活动形成的 ［1-3］。岩体中裂隙结构的存 在很大程度上降低了岩体的整体力学强度， 此外， 地 下工程开挖、 隧道建设和矿山巷道支护等常常涉及 裂纹断裂萌生、 扩展及贯通过程 ［4-5］。因此， 对裂隙岩 采矿工程 142 ChaoXing 石断裂演化过程进行表征和预测是非常必要的。 国内外学者对裂隙岩石的力学特性及断裂机制 进行了大量的研究， 并从中取得了诸多的研究成 果 ［6-10］。在室内试验方面， Wong和 Einstein［6］采用石 膏预制的单裂纹试样进行单轴压缩试验， 得到翼型 裂纹和次生裂纹的扩展演化特征。Yang 和Jing ［7］实 施了单裂缝砂岩的单轴压缩试验， 发现其力学参数 与裂隙倾角具有一定的相关性。对于室内物理试验 来说， 试样在切割钻取过程中产生的初始损伤是很 难定量表征， 尤其预制裂隙岩样大多数是通过高压 水刀加工切割， 造成试样局部的初始损伤更大， 从而 得到的试验结果离散性较大。再加上物理试验很难 定量地捕捉和表征裂纹的细观演化过程。除了室内 试验外， 众多学者借助数值模拟的方法对不同材料 的裂纹扩展演化过程和贯通机制进行研究 ［11-13］。比 如， Zhang和 Wong ［11］采用离散元数值模拟方法研究 单裂纹类岩石材料的强度、 变形及断裂扩展过程， 并 总结了多种裂纹类型和岩桥贯通模式。然而， 该研 究未考虑多裂纹工况下裂隙岩石的断裂失稳机制。 王桂林等 ［12］基于PFC2D软件对Z型裂隙砂岩的强度 和断裂行为进行了详细的分析。得到裂纹的起裂应 力水平在0.6～0.7之间， 此外， 还分析了不同裂纹的破 坏模式。然而， 该研究未考虑非贯通岩桥对裂隙岩 石断裂行为的影响。截止目前， 大多数研究者只考 虑了单一裂隙因素， 例如裂隙倾角或岩桥倾角对断 裂机制的影响 ［14-16］。然而， 实际工程中裂隙岩体的失 稳破坏通常是由多个裂隙因素耦合作用诱发的。鉴 于此， 本项目对不同裂隙倾角和岩桥倾角裂隙砂岩 的断裂破坏过程进行了一系列的模拟研究。 本研究以郑万高铁专线雷家坡隧道局部掌子 面裂纹萌生扩展为背景， 采用数值模拟反演方法分 析不同裂隙几何结构下裂纹的扩展演化过程及强 度特性， 为隧道施工过程中遇到类似情况提供一定 的指导借鉴意义， 从而采取相应的止裂和支护加固 措施。 1工程背景 郑州至万州客运专线铁路 （重庆段） 站前工程 ZWCQZQ - 8 标 段 起 讫 里 程 为 DK806 949.06～ DK819278.697， 标段正线全长12.329 km， 其中隧道 2座， 总延长12.024 km， 占正线长97.5。标段内2座 隧道均为长隧道， 属于中度风险， 其中雷家坡隧道为 标段重难点工程， 全长8 033 m， 隧道最大埋深376 m， 位于重庆市万州区大周镇与熊家镇境内。 雷家坡隧道于2016年7月21日开工， 截止目前 进口上台阶开挖1 114 m、 二衬浇筑948 m， 斜井上台 阶开挖946 m、 二衬浇筑767 m， 出口上台阶开挖781 m、 二衬浇筑 660 m； 天城隧道于 2016 年 6 月 7 日开 工， 目前进口上台阶开挖875 m、 二衬浇筑756 m， 出 口上台阶开挖1 192 m、 二衬浇筑1 044 m。从现场实 测记录得知， 隧道施工工程中伴随有不同宽度和深 度的裂隙产生， 其中 12 裂纹宽度在 0.2 mm 以下， 35 裂纹宽度 0.2～0.5 mm， 53 裂纹宽度 0.5～1.5 mm。隧道的某个掌子面见图1。 2数值试验方案 2. 1平行黏结模型原理 本研究采用PFC2D数值模拟软件， 结合平行黏 结模型 （BPM） 模拟颗粒间运动与变形的力学行为， 该 软件的计算原理基于牛顿第二定律。由于平行黏结 键不仅能够在颗粒之间传递力和向量， 而且能够传 递接触点处的作用力。因此， 平行黏结模型用于本 研究模拟裂隙砂岩的力学强度和断裂行为。平行键 模型示意图如图2所示 ［13］。 PFC2D模型中， 应力是通过作用在每个颗粒上 的平行黏结力与接触的方法获得， 平均应力向量的 计算公式［17］ σij 1 - n ∑NpV P ∑ Np∑Nc ||xi C - xi P n C,P i F C j ,（1） 式中，Np为球的质心；Nc为球的接触；V P 为颗粒体积；n 为孔隙度；xi C 和xj P 分别为颗粒质心和接触的位置； n C,P i 为接触的单位法向量；F C j 为接触作用力。 2020年第6期田茂祥等 裂纹几何配置对裂隙砂岩力学特性及破坏模式的影响 143 ChaoXing 2. 2参数标定 基于PFC2D数值模拟软件建立75 mm150 mm （宽高） 二维离散元数值模型， 颗粒直径在0.2 ～ 0.3 mm之间， 颗粒总数为50 705个， 颗粒之间的接触个 数为132 439。首先， 基于物理试验结果确定模型的 细观参数， 然后， 采用反复调试的方法， 得到最终模 拟所需的细观参数。数值模型细观参数如表1所示。 此外， 为确保整个加载过程为准静态加载， 墙体加载 速率为0.05 m/s ［12］。 2. 3数值模拟方案设计 典型的数值计算模型几何结构如图3所示。预 制裂纹长度2a为14 mm， 岩桥长度2b为16 mm， 裂纹 宽度为1.5 mm。详细模拟方案为 ①预制裂隙倾角α 固定不变， 岩桥倾角β依次为0、 30、 60、 90、 120和 150； ②岩桥倾角 β 固定不变， 裂隙倾角α依次为 15、 45和75。 3裂隙砂岩强度及变形特征 3. 1应力应变曲线 图4 （a） 、（b） 、（c） 分别为裂隙倾角15、 45、 75时 轴向应力应变曲线。 从图 4 得知， 与完整试样的应力应变曲线相 比， 含预制裂隙试样的应力应变曲线均位于其下 方。此外， 试样在应力峰前和峰后均出现不同程度 的波动现象， 而完整试样的应力应变曲线未出现 该现象。该现象的主要原因是由于含预制裂隙试样 在加载过程中裂纹尖端易形成高应力积聚区， 而且 矿岩体属于非均质材料。当加载作用力超过其最大 拉伸应力时， 导致试样局部破断失稳。随着加载的 继续， 新的承载体出现， 从而使试样的承载能力再次 增加。从图中还可得知， 完整试样的峰值应力和峰 值应变均高于裂隙试样。 3. 2强度和变形特征 通过获得不同裂隙工况下应力应变曲线中应 力的最大值， 得到不同裂隙几何结构组合下峰值应 力变化规律如图5所示。 由图5可知， 砂岩试样的峰值应力与裂隙倾角和 岩桥倾角密切相关。总体来说， 随着岩桥倾角的增 加， 峰值应力呈现出先降低后增加的变化趋势。此 外， 当岩桥倾角为60时， 峰值应力达到最小值。该 现象的主要原因为常规岩石单轴压缩时剪切破裂角 金属矿山2020年第6期总第528期 144 ChaoXing 为45φ/2， 并且该岩样的内摩擦角为34， 因此， 该 试样的剪切破裂角为62。再加上试样预制裂纹的 存在， 裂纹尖端形成局部高应力积聚区。该区域的 裂纹扩展速度要高于其他区域， 导致裂隙尖端裂纹 最先发育。 同时， 不同裂隙倾角组合下峰值应力对应的峰 值应变如图6所示。 由图6可知， 峰值应变的变化趋势与峰值应力相 似。随着岩桥倾角的增加， 峰值应变呈现出先降低 后增加的变化趋势。同一裂隙倾角下， 岩桥倾角为 60时， 峰值应变取得最小值。对应的最小峰值应变 分别为0.004 81， 0.004 63和0.006 12。 3. 3应变能演化规律 通过调用PFC2D软件中应变能计算命令， 得到 整个加载过程中试样应变能的演化规律， 不同裂隙 几何结构组合下应变能演化规律如图7所示。 由图7可知， 初始加载阶段， 不同裂隙几何结构 组合下应变能呈现出向下凹的非线性变化趋势。该 现象的主要原因为试样内含有大量的初始孔隙、 微 裂隙等， 由于孔隙、 微裂隙的闭合， 加载初期较小的 作用力会产生较大的变形量。此外， 完整试样的应 变能均大于不同裂隙几何结构下的应变能。由于应 变能的计算原理是基于应力应变曲线所围成面积 的积分， 因此， 从应力应变曲线图中可以定性地获 得应变能的大小。与试样的应力应变曲线相对 应， 3种不同裂隙倾角下应变能最小值均在岩桥倾角 为 60时获得。另外， 从图中还可得知裂隙倾角为 45时， 不同岩桥倾角的应变能演化规律接近。 3. 4阻尼耗散能演化规律 众所周知， 当物体的弹性变形超过极限变形无法 回到原始状态时， 会导致机械能的损失， 尤其当物体 产生裂纹时， 会有大部分能量消耗。因此， 分析试样 的阻尼耗散能有助于理解其断裂破坏机制。不同裂 隙几何结构组合下阻尼耗散能演化规律如图8所示。 2020年第6期田茂祥等 裂纹几何配置对裂隙砂岩力学特性及破坏模式的影响 145 ChaoXing 从图8可以看出， 整个加载过程中， 阻尼耗散能 仅当试样趋近破坏时才开始出现， 尤其在试样破裂 的瞬间， 阻尼耗散能急剧增加。从图中还可得知， 随 着预制裂隙倾角的增加， 不同岩桥倾角对应的阻尼 耗散能逐渐增大。 3. 5滑移摩擦能演化规律 滑移摩擦能是表征试样加载过程中产生裂纹时 所消耗的能量， 该参数能够间接地反应加载过程中 裂纹数量的大小程度。不同裂隙几何结构组合下阻 尼耗散能演化规律如图9所示。 从裂纹滑移摩擦能与轴向应变演化曲线得知， 当试样进入屈服阶段后， 其裂纹滑移摩擦能逐渐出 现， 随着变形的增加， 裂纹滑移摩擦能逐渐增大。当 曲线接近峰值应力时， 裂纹滑移摩擦能呈直线趋势 上升， 并且裂纹滑移摩擦能随着裂隙倾角的增加而 增加。 4裂隙砂岩裂纹扩展特征分析 由于室内试验只能借助高速相机摄像技术捕捉某 个时刻试样的宏观裂纹特征， 而对试样内部细观微裂纹 无法定量获取。从细观机理上分析岩石的损伤演化有 助于更全面地理解岩石的宏观断裂失稳机制。而 PFC2D数值模拟能够定量地表征整个加载过程中宏细 观裂纹产生的位置及数量， 因此， 采用数值模拟方法对 裂隙试样的损伤断裂演化过程进行表征显得非常必要。 4. 1裂纹扩展演化过程分析 为了分析岩体内部微裂纹与宏观应力应变曲线 之间的对应关系， 限于篇幅， 本小节选取一组典型试样 的裂纹演化过程进行分析。轴向应力、 累积总裂纹、 拉 伸裂纹和剪切裂纹应变的演化规律如图10所示。不 同应力时刻对应的裂纹起裂、 扩展和贯通演化过程如图 11所示。 从图10可以看出， 拉伸裂纹先于剪切裂纹出现， 加载前期， 未出现任何类型的裂纹， 当荷载增至峰值 应力的0.68σc时， 拉伸裂纹开始缓慢地增加， 直到接 近峰值应力， 剪切裂纹才出现。此外， 整个加载过程 中， 拉伸裂纹起到了非常重要的角色， 拉伸裂纹占总 裂纹的89.08。 从图11可以看出， 裂纹的萌生位置出现在预制 裂隙尖端， 并以拉伸裂纹的形式出现。该现象主要 是由于试样内晶粒位错和断裂等微损伤产生的。通 过对比图10和图11， 发现试样的宏观断裂过程与其 微裂纹应变曲线一一对应。随着荷载的增加， 宏 观裂纹的扩展长度逐渐增大， 新的剪切裂纹沿着预 制裂隙的方向扩展。当加载至峰值应力时， 上预制 裂隙右端剪切裂纹的扩展程度相对于拉伸裂纹较 大， 当应力降至峰后0.98σc时， 上预制裂隙左端剪切 裂纹开始扩展。随着变形继续增加， 岩桥区域被拉 伸裂纹和剪切裂纹连接贯通。当应力降至峰后 0.40σc时， 裂纹的数量及扩展程度进一步增加， 宏观 裂纹贯穿整个试样。 4. 2岩桥贯通模式分析 图12 （a） 、（b） 、（c） 分别为裂隙倾角15、 45、 75 试样破坏模式示意图。基于Wong等 ［6］对裂纹类型的 分类， 主要有拉伸裂纹、 剪切裂纹和拉剪混合裂纹 等。 从图12 （a） 、（b） 可知， 试样的破坏模式由拉剪复 合向剪切过渡再向拉剪复合模式转换。该破坏模式 能够间接地解释图 5中试样峰值应力的变化趋势。 裂纹的起裂位置发生在预制裂隙尖端， 并且试样的 金属矿山2020年第6期总第528期 146 ChaoXing ［1］ ［2］ ［3］ ［4］ 断裂区域主要集中在岩桥区域， 试样的破坏模式主 要为拉剪复合形式。随着岩桥角度的增加， 岩桥贯 通模式由间接贯通逐渐转化为直接贯通。 由图12 （c） 可知， 当岩桥倾角小于120时， 试样 的破坏模式为沿着预制裂隙形成的一条剪切断裂 带。由此可推断， 当裂隙倾角增至一定程度时， 诱发 试样断裂失稳的裂纹类型中剪切裂纹占比例逐渐增 加。该现象也能够解释裂隙倾角与其对应的峰值应 力之间的联系。此外， 岩桥贯通模式不同于 15和 45工况， 贯通模式由直接贯通变为间接贯通。 5裂隙砂岩断裂锁固体理论讨论 众所周知， 矿岩体属于非均质性较高的材料， 其 生成条件、 矿物组分、 胶结程度的不同会造成岩样内 部强度不均， 对于不同尺度的结构体， 强度较大部分 控制着整个试样的稳定性 ［18-19］。因此， 从锁固体理论 分析裂隙砂岩的断裂失稳过程即为其内部锁固体不 断失稳破坏的过程。基于阿累尼乌斯方程得到砂岩 的微破裂速度公式 ［20］ v t A0exp ■ ■ ■ ■ ■■ ■ ■ ■ ■ ■■ A1E∫ 0 t ε dt - U0 kT ，（2） 式中，v为裂纹断裂速率；A0、A1为常数；T为岩体温 度；k为玻尔兹曼常数；U0为断裂活化能；ε 为应变率； E为弹性模量。 为了计算方便， 假设岩体加载过程为恒温过程， 且加载速率为定值， 则不同时刻t的累积微裂数 ［20］为 v t ∫ 0 t vdt kTA0 A1E ε exp - U0 kT ■ ■ ■ ■ ■ ■ ■ ■ exp A1E εt kT - 1 ．（3） 由式 （3） 可知， 累积微裂数呈指数函数增加， 结 合图10分析发现， 模拟加载过程中， 试样累积微裂数 应变曲线的演化也呈现出指数函数增加。该公式 从理论上也能够解释图11中的现象， 由于预制裂隙 附近锁固体较其他区域先发育， 因此， 岩桥锁固体区 域最先失稳破坏， 导致试样断裂区域裂纹扩展速度 大于其他区域。 6结论 （1） 对比完整试样的应力应变曲线特征， 含裂 隙试样的应力应变曲线在峰值附近出现不同程度 的波动， 另外， 裂隙试样的峰值应力和峰值应变均出 现不同程度的降低。峰值应力和峰值应变的变化规 律一致， 当裂隙倾角不变时， 二者随着岩桥角度的增 加呈现出先降低后增加的趋势。当岩桥角度不变 时， 二者随着裂隙倾角的增加而增加。 （2） 裂隙砂岩的应变能、 阻尼耗散能和滑移摩擦 能与裂隙的几何结构配置有关， 均在岩桥倾角为60 时取得最小值。此外， 阻尼耗散能从应力接近峰值 时开始出现， 裂纹滑移摩擦能从弹性阶段逐渐增加， 尤其在接近峰值处急剧增加。 （3） 当试样预制裂隙倾角为15和45时， 随着岩 桥角度的增加， 岩桥贯通模式由间接贯通向直接贯 通转换。但当裂隙倾角为75时， 随着岩桥角度的增 加， 岩桥贯通模式由直接贯通变为间接贯通。 （4） 通过理论分析和模拟结果对比发现， 加载过 程中试样累积微裂纹随时间演化均呈现指数函数的 形式。 参 考 文 献 Hoek E，Martin C D． Fracture initiation and propagation in intact rock-A review［J］ ． Journal of Rock Mechanics and Geotechnical Engineering， 2014， 6287-300． Horii H，Nemat Nasser S． Compression-induced microcrack growth in brittle solids Axial splitting and shear failure［J］ ． Journal of Geophysical Research， 1985， 903105-3125． Guo Weiyao， Tan Yunliang， Yu Fenghai， et al． Mechanical behav⁃ ior of rock-coal-rock specimens with different coal thicknesses［J］ ． Geomechanics and Engineering， 2018， 15 （4） 1017-1027． Zhao Zhihong，Zhou Dong． Mechanical properties and failure modes of rock samples with grout-infilled flaws A particle mechan⁃ 2020年第6期田茂祥等 裂纹几何配置对裂隙砂岩力学特性及破坏模式的影响 147 ChaoXing ［5］ ［6］ ［7］ ［8］ ［9］ ［10］ ［11］ ［12］ ics modeling［J］ ． Journal of Natural Gas Science and Engineering， 2016， 34702-715． Guo Weiyao，Gu Qingheng，Tan Yunliang，et al． Case studies of rock bursts in tectonic areas with facies change［J］ ． Energies， 2019， 121330-1340． Wong L N Y，Einstein H H． Systematic uation of cracking be⁃ havior in specimens containing single flaws under uniaxial com⁃ pression［J］ ． International Journal of Rock Mechanics and Mining Sciences， 2009， 46239-249． Yang Shengqi，Jing Hongwei． Strength failure and crack coales⁃ cence behavior of brittle sandstone samples containing a single fis⁃ sure under uniaxial compression［J］ ． International Journal of Frac⁃ ture， 2011， 168227-250． 侯恩科，从通，谢晓深，等． 基于颗粒流的浅埋双煤层斜交 开采地表裂缝发育特征 ［J］ ． 采矿与岩层控制工程学报， 2020， 2 （1） 1-9． Hou Enke，Cong Tong，Xie Xiaoshen，et al． Ground surface frac⁃ ture development characteristics of shallow double coal seam stag⁃ gered mining based on particle flow［J］ ． Journal of Mining and Strata Control Engineering， 2020， 2 （1） 1-9． 庄贤鹏， 王述红， 王凯毅， 等． 断续双裂隙砂岩的强度与起裂、 贯通模式的研究 ［J］ ． 科学技术与工程， 2019， 19 （20） 122-127． Zhuang Xianpeng，Wang Shuhong，Wang Kaiyi， et al．Experimen⁃ tal studies on peak strength， initiation and coalescence mode of sandstone specimen containing two pre-existing fissures［J］ ． Sci⁃ ence Technology and Engineering， 2019， 19 （20） 122-127． 王桂林， 张亮， 许明， 等． 单轴压缩下非贯通岩体损伤破坏 能量演化机制研究［J］ ． 岩土工程学报， 2019， 41 （4） 640-647． Wang Guilin，Zhnag Liang，Xu Ming， et al． Energy damage evolu⁃ tion mechanism of non-across jointed rock mass under uniaxial compression［J］ ． Chinese Journal of Geotechnical Engineering， 2019， 41 （4） 640-647． Zhang Xiaopeng，Wong Louisngaiyuen． Cracking processes in rock-like material containing a single flaw under uniaxial compres⁃ sion a numerical study based on parallel bonded-particle model ap⁃ proach［J］ ． Rock Mechanics and Rock Engineering，2012，45 711-737． 王桂林， 梁再勇， 张亮， 等． Z型裂隙对砂岩强度和破裂行为 影响机制研究［J］ ． 岩土力学， 2018， 39 （S2） 390-397． Wang Guilin，Liang Zaiyong，Zhang Liang，et al． Study of influ⁃ 金属矿山2020年第6期总第528期 148 ChaoXing ［13］ ［14］ ［15］ ［16］ ［17］ ［18］ ［19］ ［20］ ence mechanism of Z-type fissure on sandstone strength and frac⁃ ture behavior［J］ ． Rock Soil Mechanics，2018，39（S2） 390- 397． Cho N，Martin C D，Sego D C． A clumped particle model for rock ［J］ ． International Journal of Rock Mechanics and Mining Scienc⁃ es， 2007， 44 （7） 997-1010． 单仁亮，白瑶，孙鹏飞，等． 裂隙红砂岩冻胀力特性试验研 究［J］ ． 煤炭学报， 2019， 44 （6） 1742-1752． Shan Renliang，Bai Yao，Sun Pengfei，et al． Experimental study on frost heaving pressure properties in fractured red sandstone［J］ ． Journal of China Coal Society， 2019， 44 （6） 1742-1752． 赵国彦，李振阳，吴浩，等． 含非贯通裂隙砂岩的动力破坏 特性研究［J］ ． 岩土力学， 2019， 40 （S1） 74-81． Zhao Guoyan，Li Zhenyang，Wu Hao，et al． Dynamic failure char⁃ acteristics of sandstone with non-penetrating cracks［J］ ． Rock Soil Mechanics， 2019， 40 （S1） 74-81． 王笑然，王恩元，刘晓斐，等． 裂隙砂岩裂纹扩展声发射响应 及速率效应研究［J］ ． 岩石力学与工程学报， 2018， 37 （6） 1447- 1458． Wang Xiaoran，Wang Enyuan，Liu Xiaofei，et al． Macro-crack propagation process and corresponding AE behaviors of fractured sandstone under different loading rates ［J］ ． Chinese Journal of Rock Mechanics Engineering， 2018， 37 （6） 1447-1458． Potyonbdy D O，Cundall P A． A bonded-particle model for rock ［J］ ． International Journal of Rock Mechanics and Mining Scienc⁃ es， 2004， 411329-1364． 黄润秋， 陈国庆， 唐鹏． 基于动态演化特征的锁固段型岩质滑 坡前兆信息研究［J］ ． 岩石力学与工程学报，2017， 36 （3） 521- 533． Huang Runqiu， Chen Guoqing，Tang Peng． Precursor ination of locking segment landslides based on transient characteristics ［J］ ． Chinese Journal of Rock Mechanics and Engineering，2017， 36 （3） 521-533． 秦四清， 徐锡伟， 胡平， 等． 孕震断层的多锁固段脆性破裂机 制与地震预测新方法的探索［J］ ． 地球物理学报， 2010， 53 （4） 1001-1014． Qin Siqing， Xu Xiwei， Hu Ping， et al． Brittle failure mechanism of multiple locked paths in a seismogenic fault system and exploration on a new way for earthquake prediction［J］ ． Chinese Journal of Geophysics， 2010， 53 （4） 1001-1014． Zhurkov S N，Kuksenko V S，Petrov V A． Principles of the kinetic approach of fracture prediction［J］ ． Theoretical and Applied Frac⁃ ture Mechanics， 1984， 1 （3） 271-274． 2020年第6期田茂祥等 裂纹几何配置对裂隙砂岩力学特性及破坏模式的影响 149 ChaoXing