含孔洞大理岩破坏特性的颗粒流分析

工程科学学报，第39卷，第12期：1791-1801,2017年12月 Chinese Journal of Engineering,Vol.39,No.12:1791-1801,December 2017 D0l:10.13374/j.issn2095-9389.2017.12.003;http://journals..ustb.edu.cn 含孔洞大理岩破坏特性的颗粒流分析 黎崇金四，李夕兵，李地元 中南大学资源与安全工程学院，长沙410083 ☒通讯作者，E-mail:lcj2015@csu.edu.cm 摘要基于室内单轴压缩试验结果，利用颗粒流程序P℉C2D,模拟含预制孔洞大理岩在单轴和双轴压缩条件下的破坏过 程，分析预制孔洞形状、围压大小以及岩石非均质性对大理岩力学特性和裂纹扩展的影响.数值结果表明：与完整大理岩试 样相比，含孔洞试样的峰值强度显著降低，降低程度与孔洞形状有关：围压对含孔洞大理岩试样的力学特性和裂纹扩展有显 著影响，含孔洞试样的峰值强度随围压的增加而增加，但偏应力峰值随围压的增加呈先增大后减小的变化趋势：试样的破坏 模式与孔洞形状相关，含圆形孔洞试样为类X型剪切破坏，含矩形孔洞或马蹄形孔洞试样为对角剪切破坏；岩石内部的矿物 结核影响了裂纹的扩展路径，从而改变试样的宏观破坏模式.微观机理分析表明：孔洞周边裂纹的萌生与扩展过程伴随着应 力集中区的释放与转移：含孔洞试样的宏观裂纹有3种模式：孔壁剥落、拉伸裂纹和压剪裂纹 关键词颗粒流程序：孔洞试样：裂纹扩展：岩石非均质性：双轴压缩 分类号TU452 Particle flow analysis of fracture characteristics of marble with a single hole LI Chongjin,LI Xi-bing,LI Di-yuan School of Resources and Safety Engineering.Center South University,Changsha 410083,China Corresponding author,E-mail:lcj2015@csu.edu.cn ABSTRACT Based on uniaxial compression laboratory results,the failure process of marble with a single hole was simulated using the particle flow code PFC2D under uniaxial and biaxial compressions.The influence of the pre-existing hole shape,confining pressure and the rock heterogeneity on the mechanical properties of marble and the coalescence of cracks was analyzed.Numerical results show that compared to an intact marble specimen,the peak strength of a specimen containing a single hole reduces significantly,and the ex- tent of reduction is related to the hole shape.The confining pressure has a significant effect on the mechanical properties of marble and the coalescence of cracks.The peak strength increases with the increase of confining pressure,while the peak deviator stress increases first and then decreases with the confining pressure increasing.Specimens containing circular holes exhibit X-type shear failure,while specimens containing rectangular or U-shaped holes exhibit diagonal shear failure.The mineral nodules in the rock specimen signifi- cantly affect the propagation of cracks leading to a change in the failure mode of marble specimen.Analysis results of the microscopic mechanisms show that the initiation and propagation of cracks around the hole is always accompanied by the release and transfer of stress from a concentrated zone.Three types of macroscopic cracks can be classified in specimen containing a single hole as hole sur- face spalling,tensile crack,and compressive shear crack. KEY WORDS particle flow code (PFC);single hole specimen:crack propagation:rock heterogeneity:biaxial compression 岩石作为一种天然的非均质材料，其内部含有许开挖等地下工程也会在岩体中形成孔洞.外荷载作用 多宏观孔洞和微观裂隙、节理等缺陷：地下采矿与隧道下，裂纹沿孔洞周边起裂、扩展、贯通进而导致岩石整 收稿日期：201701-16 基金项目：国家自然科学基金资助项目(11472311,41630642)：中南大学中央高校基本科研业务费专项资金资助项目(2017zzt538)

工程科学学报，第 39 卷，第 12 期: 1791--1801，2017 年 12 月 Chinese Journal of Engineering，Vol． 39，No． 12: 1791--1801，December 2017 DOI: 10． 13374 /j． issn2095--9389． 2017． 12． 003; http: / /journals． ustb． edu． cn 含孔洞大理岩破坏特性的颗粒流分析 黎崇金，李夕兵，李地元 中南大学资源与安全工程学院，长沙 410083 通讯作者，E-mail: lcj2015@ csu． edu． cn 摘 要 基于室内单轴压缩试验结果，利用颗粒流程序 PFC2D，模拟含预制孔洞大理岩在单轴和双轴压缩条件下的破坏过 程，分析预制孔洞形状、围压大小以及岩石非均质性对大理岩力学特性和裂纹扩展的影响． 数值结果表明: 与完整大理岩试 样相比，含孔洞试样的峰值强度显著降低，降低程度与孔洞形状有关; 围压对含孔洞大理岩试样的力学特性和裂纹扩展有显 著影响，含孔洞试样的峰值强度随围压的增加而增加，但偏应力峰值随围压的增加呈先增大后减小的变化趋势; 试样的破坏 模式与孔洞形状相关，含圆形孔洞试样为类 X 型剪切破坏，含矩形孔洞或马蹄形孔洞试样为对角剪切破坏; 岩石内部的矿物 结核影响了裂纹的扩展路径，从而改变试样的宏观破坏模式． 微观机理分析表明: 孔洞周边裂纹的萌生与扩展过程伴随着应 力集中区的释放与转移; 含孔洞试样的宏观裂纹有 3 种模式: 孔壁剥落、拉伸裂纹和压剪裂纹． 关键词 颗粒流程序; 孔洞试样; 裂纹扩展; 岩石非均质性; 双轴压缩 分类号 TU452 Particle flow analysis of fracture characteristics of marble with a single hole LI Chong-jin ，LI Xi-bing，LI Di-yuan School of Ｒesources and Safety Engineering，Center South University，Changsha 410083，China Corresponding author，E-mail: lcj2015@ csu． edu． cn ABSTＲACT Based on uniaxial compression laboratory results，the failure process of marble with a single hole was simulated using the particle flow code PFC2D under uniaxial and biaxial compressions． The influence of the pre-existing hole shape，confining pressure and the rock heterogeneity on the mechanical properties of marble and the coalescence of cracks was analyzed． Numerical results show that compared to an intact marble specimen，the peak strength of a specimen containing a single hole reduces significantly，and the ex￾tent of reduction is related to the hole shape． The confining pressure has a significant effect on the mechanical properties of marble and the coalescence of cracks． The peak strength increases with the increase of confining pressure，while the peak deviator stress increases first and then decreases with the confining pressure increasing． Specimens containing circular holes exhibit X-type shear failure，while specimens containing rectangular or U-shaped holes exhibit diagonal shear failure． The mineral nodules in the rock specimen signifi￾cantly affect the propagation of cracks leading to a change in the failure mode of marble specimen． Analysis results of the microscopic mechanisms show that the initiation and propagation of cracks around the hole is always accompanied by the release and transfer of stress from a concentrated zone． Three types of macroscopic cracks can be classified in specimen containing a single hole as hole sur￾face spalling，tensile crack，and compressive shear crack． KEY WOＲDS particle flow code ( PFC) ; single hole specimen; crack propagation; rock heterogeneity; biaxial compression 收稿日期: 2017--01--16 基金项目: 国家自然科学基金资助项目( 11472311，41630642) ; 中南大学中央高校基本科研业务费专项资金资助项目( 2017zzts538) 岩石作为一种天然的非均质材料，其内部含有许 多宏观孔洞和微观裂隙、节理等缺陷; 地下采矿与隧道 开挖等地下工程也会在岩体中形成孔洞． 外荷载作用 下，裂纹沿孔洞周边起裂、扩展、贯通进而导致岩石整

·1792 工程科学学报，第39卷，第12期 体破坏是岩石类材料破坏的重要表现形式.因此，开 它不需要预先定义破坏准则，仅由颗粒间的微观参数 展含孔洞岩石破坏特性的研究对岩土工程的安全稳定 来表现材料的宏观力学行为阁 具有重要的工程意义.近年来，国内外学者对含预制 基于此，本文以室内物理试验结果为基础，利用颗 孔洞岩石试样的破坏特性进行了大量的理论分析、室 粒流程序PC2D模拟含不同孔洞形状的大理岩在单 内试验和数值模拟，取得了许多成果.Germanovich与 轴和双轴压缩条件下的破坏过程，进一步探索孔洞形 Dyskin四基于断裂力学理论分析了含孔洞缺陷岩石的 状、围压条件以及岩石非均质性对裂纹扩展演化的 断裂破坏机理.Fakhimi等网对含预制圆形孔洞的砂 影响. 岩进行了双轴压缩试验，研究了含圆形孔洞试样破坏 过程的声发射特性；Zhu等因利用有限元数值模拟软 1试验方案与数值模型建立 件RFPA2D研究了不同侧压系数条件下，含孔洞试样 1.1室内试验概况 的破坏过程，结果表明侧压系数对含孔洞岩石的破坏 试验所用大理岩取自意大利中北部的卡拉拉市， 模式有显著地影响：Wog等研究了含预制圆形孔洞 其矿物成分主要为方解石和白云石，平均密度为2700 花岗岩在单轴压缩下劈裂裂纹的演化过程，结果表明 kgm3.经实验室精密加工，制备成完整的大理岩试 劈裂裂纹起裂于孔洞内表面的上下两端，并沿平行于 样，尺寸为120mm×150mm×20mm.在此基础上，利 加载方向扩展，其演化过程受到试样尺寸、孔洞尺寸和 用水刀切割在完整试样的几何中心预制了圆形、矩形 岩石非均质性的影响：Wang等利用含预制孔洞的花 和马蹄形三种孔洞形状，三种孔洞的面积近似相等，模 岗岩进行室内试验，并利用RFPA2D进行数值模拟，研 型及孔洞尺寸如图1所示.采用NSTRON-1346型液压 究了开挖损伤区(EDZ)的形成机理：杨圣奇等-对 同服控制试验机对大理岩试样进行单轴压缩试验，试验 含孔洞缺陷的大理岩、砂岩进行了试验研究和数值模 过程采用位移加载控制，加载速率为0.l8 mm*min. 拟，研究了含缺陷岩石在单轴压缩条件下的破坏过程： a 刘招伟与李元海图利用数字散斑相关量测方法 (DSCM)和数字照相量测软件系统Photolnfor研究了 单轴压缩下含圆形孔洞岩石试样的变形破坏规律： Zhao等和Wang等0分别采用声发射技术和RF- PA3D数值模拟软件，对含孔洞岩石的变形破坏规律 进行了研究，结果表明，圆孔的破坏主要为圆孔上下孔 壁的拉伸破坏和圆孔两帮裂纹扩展：张社荣等四利用 Marble I P℉C3D对含圆形孔洞试样进行了单轴、双轴、和三轴 试验，研究了不同应力路径下岩石试样的破坏特征；李 6 y 地元等对两侧预制了方形孔洞的板状花岗岩进行 2r=25 mm 上a=22mm大 单轴压缩试验，并利用FLAC3D进行数值模拟，研究了 8 试样的声发射和变形破坏特征：李地元等国和李夕兵 等分别利用高速摄影仪和核磁共振(NMR)技术研 a=20 mm 究了冲击荷载作用下含孔洞岩石试样的动态变形破坏 特性. 图1大理岩试样及孔洞尺寸.(a)大理岩试样：(b)P℉C模型： 虽然前人对含预制孔洞岩石试样的破坏特性进行 (c)圆形孔洞：(d)矩形孔洞：（©）马蹄形孔洞 比较深入的研究，但由于实验条件的限制，室内试验多 Fig.1 Marble specimens and hole sizes:(a)marble specimen: (b)PFC model:(c)circular hole:(d)rectangular hole:(e)U- 以单轴压缩试验为主，这与工程中岩体的实际受力情 shaped hole 况存在较大差异的.近年来，随着计算机科学的快速 发展，使利用数值方法研究岩石裂纹扩展演化机理成 1.2数值模型建立及微观参数标定 为可能6一切.但以往的数值模拟大都使用基于有限单 以室内物理模型为基础，利用PFC2D建立120mm 元法或有限差分法的计算机软件来模拟裂纹的扩展， ×150mm的二维离散元模型，如图1(b)所示.在PFC 这些软件必须预先定义岩石的破坏准则，而且为了能 中，颗粒之间的黏结有接触黏结和平行黏结两种基本 够模拟岩石的非连续性而不得不采用大量的假设，进 模型，由于平行黏结模型更适用于模拟硬岩类材 而导致采用不同破坏准则和不同假设得到不同试验结 料，故本文选用平行黏结模型.利用P℉C模型进行 果的现象.而基于离散单元法的颗粒流程序P℉C则是 数值模拟之前，需要先标定试样的微观参数.以岩石 能够很好模拟岩石非均匀性和非连续性的数值软件， 试样的单轴压缩试验数据为基础，通过“试错法”标

工程科学学报，第 39 卷，第 12 期 体破坏是岩石类材料破坏的重要表现形式． 因此，开 展含孔洞岩石破坏特性的研究对岩土工程的安全稳定 具有重要的工程意义． 近年来，国内外学者对含预制 孔洞岩石试样的破坏特性进行了大量的理论分析、室 内试验和数值模拟，取得了许多成果． Germanovich 与 Dyskin［1］基于断裂力学理论分析了含孔洞缺陷岩石的 断裂破坏机理． Fakhimi 等［2］对含预制圆形孔洞的砂 岩进行了双轴压缩试验，研究了含圆形孔洞试样破坏 过程的声发射特性; Zhu 等［3］利用有限元数值模拟软 件 ＲFPA2D 研究了不同侧压系数条件下，含孔洞试样 的破坏过程，结果表明侧压系数对含孔洞岩石的破坏 模式有显著地影响; Wong 等［4］研究了含预制圆形孔洞 花岗岩在单轴压缩下劈裂裂纹的演化过程，结果表明 劈裂裂纹起裂于孔洞内表面的上下两端，并沿平行于 加载方向扩展，其演化过程受到试样尺寸、孔洞尺寸和 岩石非均质性的影响; Wang 等［5］利用含预制孔洞的花 岗岩进行室内试验，并利用 ＲFPA2D 进行数值模拟，研 究了开挖损伤区( EDZ) 的形成机理; 杨圣奇等［6--7］对 含孔洞缺陷的大理岩、砂岩进行了试验研究和数值模 拟，研究了含缺陷岩石在单轴压缩条件下的破坏过程; 刘招 伟 与 李 元 海［8］ 利 用 数 字 散 斑 相 关 量 测 方 法 ( DSCM) 和数字照相量测软件系统 PhotoInfor 研究了 单轴压缩下含圆形孔洞岩石试样的变形破坏规律; Zhao 等［9］和 Wang 等［10］ 分别采用声发射技术和 ＲF￾PA3D 数值模拟软件，对含孔洞岩石的变形破坏规律 进行了研究，结果表明，圆孔的破坏主要为圆孔上下孔 壁的拉伸破坏和圆孔两帮裂纹扩展; 张社荣等［11］利用 PFC3D 对含圆形孔洞试样进行了单轴、双轴、和三轴 试验，研究了不同应力路径下岩石试样的破坏特征; 李 地元等［12］对两侧预制了方形孔洞的板状花岗岩进行 单轴压缩试验，并利用 FLAC3D 进行数值模拟，研究了 试样的声发射和变形破坏特征; 李地元等［13］和李夕兵 等［14］分别利用高速摄影仪和核磁共振( NMＲ) 技术研 究了冲击荷载作用下含孔洞岩石试样的动态变形破坏 特性． 虽然前人对含预制孔洞岩石试样的破坏特性进行 比较深入的研究，但由于实验条件的限制，室内试验多 以单轴压缩试验为主，这与工程中岩体的实际受力情 况存在较大差异［15］． 近年来，随着计算机科学的快速 发展，使利用数值方法研究岩石裂纹扩展演化机理成 为可能［16--17］． 但以往的数值模拟大都使用基于有限单 元法或有限差分法的计算机软件来模拟裂纹的扩展， 这些软件必须预先定义岩石的破坏准则，而且为了能 够模拟岩石的非连续性而不得不采用大量的假设，进 而导致采用不同破坏准则和不同假设得到不同试验结 果的现象． 而基于离散单元法的颗粒流程序 PFC 则是 能够很好模拟岩石非均匀性和非连续性的数值软件， 它不需要预先定义破坏准则，仅由颗粒间的微观参数 来表现材料的宏观力学行为［18］． 基于此，本文以室内物理试验结果为基础，利用颗 粒流程序 PFC2D 模拟含不同孔洞形状的大理岩在单 轴和双轴压缩条件下的破坏过程，进一步探索孔洞形 状、围压条件以及岩石非均质性对裂纹扩展演化的 影响． 1 试验方案与数值模型建立 1. 1 室内试验概况 试验所用大理岩取自意大利中北部的卡拉拉市， 其矿物成分主要为方解石和白云石，平均密度为 2700 kg·m － 3 ． 经实验室精密加工，制备成完整的大理岩试 样，尺寸为 120 mm × 150 mm × 20 mm． 在此基础上，利 用水刀切割在完整试样的几何中心预制了圆形、矩形 和马蹄形三种孔洞形状，三种孔洞的面积近似相等，模 型及孔洞尺寸如图 1 所示． 采用 INSTＲON--1346 型液压 伺服控制试验机对大理岩试样进行单轴压缩试验，试验 过程采用位移加载控制，加载速率为0. 18 mm·min － 1 ． 图 1 大理岩试样及孔洞尺寸． ( a) 大理岩试样; ( b) PFC 模型; ( c) 圆形孔洞; ( d) 矩形孔洞; ( e) 马蹄形孔洞 Fig． 1 Marble specimens and hole sizes: ( a) marble specimen; ( b) PFC model; ( c) circular hole; ( d) rectangular hole; ( e) U￾shaped hole 1. 2 数值模型建立及微观参数标定 以室内物理模型为基础，利用 PFC2D 建立 120 mm × 150 mm 的二维离散元模型，如图 1( b) 所示． 在 PFC 中，颗粒之间的黏结有接触黏结和平行黏结两种基本 模型，由 于 平 行 黏 结 模 型 更 适 用 于 模 拟 硬 岩 类 材 料［19］，故本文选用平行黏结模型． 利用 PFC 模型进行 数值模拟之前，需要先标定试样的微观参数． 以岩石 试样的单轴压缩试验数据为基础，通过“试错法”标 · 2971 ·

黎崇金等：含孔洞大理岩破坏特性的颗粒流分析 ·1793· 定，最终得到的微观参数如表1所示，k和k,分别为正 行单轴和双轴压缩试验.单轴压缩试验首先删除左右 向和切向接触刚度，k。和k。分别为正向和切向平行黏 两侧的墙体，然后以0.05m·s的速度移动上下墙体 结刚度.表2为室内试验与数值模拟得到的试样宏观 进行加载：双轴压缩试验首先利用servo程序控制边界 力学参数，从表中可知，数值模型的宏观物理力学参数 墙体移动，使模型达到指定围压（本次试验所取围压 与室内试验结果基本一致，说明所标定的微观参数比 大小为2.5、5.0、7.5和10MPa),然后保持围压不变， 较合理，能够反映岩石的宏观力学特性. 同样以0.05ms的速度进行竖向加载.试验过程中， 通过删除颗粒的方法在模型内预制孔洞，然后进 当轴向的残余应力为峰值应力的60%时停止加载. 表1试样的微观力学参数 Table 1 Microscopic parameters of the rock specimens 摩擦系数， 接触模量， 接触刚度比， 平行黏结模量，平行黏结刚度平行黏结法向平行黏结切向 粒径/mm E./GPa knk。 E./GPa 比，kn店 强度，G。/MPa 强度，F。MPa 0.4-0.64 0.25 14.32 2.5 14.32 2.5 42.15±4.22 42.15±4.22 表2 单轴压缩条件下岩石试样的宏观力学参数 应力-变曲线及峰值应力.从图2(a)可以看出，单轴 Table 2 Macroscopic mechanical parameters of the rock specimens un- 压缩条件下大理岩试样的应力一应变曲线有明显的压 der uniaxial compression 密阶段、弹性阶段和屈服阶段.试样峰后应力快速跌 密度，p/ 单轴抗压强弹性模量， 泊松比， 试验条件 (kg'm-3) 度，d。/MPa E/GPa 落，表现为显著的脆性破坏特征；从峰后曲线来看，完 室内试验 2700 51.26 17.04 0.25 整试样的跌落较快，表明完整试样比含孔洞试样表现 数值模拟 3000 51.95 16.97 0.25 出更明显的脆性特征.值得注意的是，由于P℉C2D通 过压密的二维圆盘来模拟类岩石材料，难以完全准确 2试验结果与分析 地再现结构复杂的三维大理岩试样，所以由PC2D获 得的应力一应变曲线不能体现初始压密段：而且由模 2.1单轴压缩条件下的破坏特性 拟得到的峰值强度要略高于室内试验结果，但两者表 图2为单轴压缩条件下，室内试验与数值模拟的 60o 50 马蹄形孔试样 50 马蹄形孔试样 40 美 完整试样 完整试样 304 圆形孔试样 20 圆形孔 试样 20 矩形孔试样 10 矩形孔试样 0 34 5■ 6 78 2 轴向应变，8/103 轴向应变，e,103 ■室内试验结果 50 ☐数值模拟结果 40 0 20 圆形孔 矩形孔马蹄形孔 试样 试样 试样 试样 试样类型 图2单轴压缩条件下室内试验与数值模拟对比.(a)室内试验应力-应变曲线：(b)数值模拟应力-应变曲线：()峰值应力 Fig.2 Comparison of experiment and numerical simulation results under uniaxial compression:(a)stress-strain curves in laboratory tests:(b) stress-strain curves in numerical simulations:(c)peak stress

黎崇金等: 含孔洞大理岩破坏特性的颗粒流分析 定，最终得到的微观参数如表 1 所示，kn和 ks分别为正 向和切向接触刚度，kn 和ks 分别为正向和切向平行黏 结刚度． 表 2 为室内试验与数值模拟得到的试样宏观 力学参数，从表中可知，数值模型的宏观物理力学参数 与室内试验结果基本一致，说明所标定的微观参数比 较合理，能够反映岩石的宏观力学特性． 通过删除颗粒的方法在模型内预制孔洞，然后进 行单轴和双轴压缩试验． 单轴压缩试验首先删除左右 两侧的墙体，然后以 0. 05 m·s － 1的速度移动上下墙体 进行加载; 双轴压缩试验首先利用 servo 程序控制边界 墙体移动，使模型达到指定围压 ( 本次试验所取围压 大小为 2. 5、5. 0、7. 5 和 10 MPa) ，然后保持围压不变， 同样以 0. 05 m·s － 1的速度进行竖向加载． 试验过程中， 当轴向的残余应力为峰值应力的60%时停止加载． 表 1 试样的微观力学参数 Table 1 Microscopic parameters of the rock specimens 粒径/mm 摩擦系数， υ 接触模量， Ec /GPa 接触刚度比， kn /ks 平行黏结模量， Ec /GPa 平行黏结刚度 比，kn /ks 平行黏结法向 强度，σc /MPa 平行黏结切向 强度，τc /MPa 0. 4 ～ 0. 64 0. 25 14. 32 2. 5 14. 32 2. 5 42. 15 ± 4. 22 42. 15 ± 4. 22 表 2 单轴压缩条件下岩石试样的宏观力学参数 Table 2 Macroscopic mechanical parameters of the rock specimens un￾der uniaxial compression 试验条件 密度，ρ / ( kg·m － 3 ) 单轴抗压强 度，σc /MPa 弹性模量， E /GPa 泊松比， υ 室内试验 2700 51. 26 17. 04 0. 25 数值模拟 3000 51. 95 16. 97 0. 25 图 2 单轴压缩条件下室内试验与数值模拟对比． ( a) 室内试验应力--应变曲线; ( b) 数值模拟应力--应变曲线; ( c) 峰值应力 Fig． 2 Comparison of experiment and numerical simulation results under uniaxial compression: ( a) stress--strain curves in laboratory tests; ( b) stress--strain curves in numerical simulations; ( c) peak stress 2 试验结果与分析 2. 1 单轴压缩条件下的破坏特性 图 2 为单轴压缩条件下，室内试验与数值模拟的 应力--变曲线及峰值应力． 从图 2( a) 可以看出，单轴 压缩条件下大理岩试样的应力--应变曲线有明显的压 密阶段、弹性阶段和屈服阶段． 试样峰后应力快速跌 落，表现为显著的脆性破坏特征; 从峰后曲线来看，完 整试样的跌落较快，表明完整试样比含孔洞试样表现 出更明显的脆性特征． 值得注意的是，由于 PFC2D 通 过压密的二维圆盘来模拟类岩石材料，难以完全准确 地再现结构复杂的三维大理岩试样，所以由 PFC2D 获 得的应力--应变曲线不能体现初始压密段; 而且由模 拟得到的峰值强度要略高于室内试验结果，但两者表 · 3971 ·

·1794· 工程科学学报，第39卷，第12期 现出相同的变化趋势，所以不影响本文对含孔洞试样 同样先在孔洞上下端出现：孔洞下方的转角处由于压 力学变化规律的研究.从图2(c)可知，各试样的峰值 应力集中出现了初始剪切裂纹：随着轴向荷载增加，剪 应力大小关系为：完整试样>马蹄形孔试样>圆形孔 切裂纹继续扩展形成宏观破坏裂纹导致试样发生剪切 试样>矩形孔试样 破坏，孔洞左侧同样出现了明显的剥落现象 在单轴压缩试验过程中，采用高速摄影仪记录了 作为对比，图4给出了单轴压缩条件下含孔洞试 岩石的破坏过程，图3给出了含孔洞大理岩试样的裂 样裂纹扩展的数值模拟结果，图中黑色和红色分别代 纹扩展过程.如图3(a)所示，单轴压缩条件下，含圆 表张拉型微裂纹和剪切型微裂纹.如图4(a)所示，单 形孔洞的大理岩试样首先在孔洞上下边缘出现初始拉 轴压缩条件下，含圆形孔洞试样首先在孔洞上边缘出 伸裂纹白斑，在孔洞左侧出现了初始剪切裂纹，剪切裂 现了拉伸裂纹：当轴压增加到38.5MPa时，孔洞下边 纹向左上方延伸形成宏观破坏裂纹，同时孔洞左下方 缘出现了另一条拉伸裂纹，拉伸裂纹沿平行于加载方 出现了远场拉伸裂纹；随后孔洞右侧的剪切裂纹向右 向延伸，同时孔洞左下方形成了一条远场拉伸裂纹，试 下方延伸，导致试样发生宏观剪切破坏.对于矩形孔 样左下边角被压碎，形成了一条宏观裂纹：轴压增加到 洞试样，如图3(b)所示，初始拉伸裂纹白斑首先在孔 峰值应力后，微观裂纹快速增加，剪切裂纹从孔洞左右 洞上下端出现：随后孔洞左下角出现了初始剪切裂纹， 两侧沿试样对角方向延伸，导致试样发生剪切破坏，但 初始剪切裂纹向左上方延伸形成宏观破坏裂纹，孔洞 孔洞上下两端的拉伸裂纹没有延伸至试样外边缘.由 左侧出现了明显的剥落现象：最后另一条宏观剪切裂 图4()可知，矩形孔洞试样首先在孔洞上边缘出现了 纹在孔洞右下方形成，导致试件发生失稳破坏.对于 初始拉伸裂纹：随后，另一条拉伸裂纹出现在孔洞下边 马蹄形孔洞试样，如图3（©）所示，初始拉伸裂纹白斑 缘，孔洞周边转角处由于压应力集中出现了少量微观 一宏观破坏裂纹 初始拉 神裂纹 初始剪切裂纹 远场拉伸裂纹 孔壁剥落 宏观破坏裂纹 初始拉 伸裂纹 T2 T2 T2 初始剪切裂纹 孔壁剥落 ←一宏观破坏裂纹 (c 初始拉 伸裂纹 初始剪切裂纹 图3单轴压缩条件下含孔洞大理岩试样的裂纹扩展过程（室内试验结果）.()圆形孔洞：(b)矩形孔洞：()马蹄形孔洞 Fig.3 Crack coalescence process of marble specimens containing a single hole under uniaxial compression (laboratory results):(a)circular hole: (b)rectangular hole:(c)U-shaped hole

工程科学学报，第 39 卷，第 12 期 现出相同的变化趋势，所以不影响本文对含孔洞试样 力学变化规律的研究． 从图 2( c) 可知，各试样的峰值 应力大小关系为: 完整试样 ＞ 马蹄形孔试样 ＞ 圆形孔 试样 ＞ 矩形孔试样． 在单轴压缩试验过程中，采用高速摄影仪记录了 岩石的破坏过程，图 3 给出了含孔洞大理岩试样的裂 纹扩展过程． 如图 3( a) 所示，单轴压缩条件下，含圆 图 3 单轴压缩条件下含孔洞大理岩试样的裂纹扩展过程 ( 室内试验结果) ． ( a) 圆形孔洞; ( b) 矩形孔洞; ( c) 马蹄形孔洞 Fig． 3 Crack coalescence process of marble specimens containing a single hole under uniaxial compression ( laboratory results) : ( a) circular hole; ( b) rectangular hole; ( c) U-shaped hole 形孔洞的大理岩试样首先在孔洞上下边缘出现初始拉 伸裂纹白斑，在孔洞左侧出现了初始剪切裂纹，剪切裂 纹向左上方延伸形成宏观破坏裂纹，同时孔洞左下方 出现了远场拉伸裂纹; 随后孔洞右侧的剪切裂纹向右 下方延伸，导致试样发生宏观剪切破坏． 对于矩形孔 洞试样，如图 3( b) 所示，初始拉伸裂纹白斑首先在孔 洞上下端出现; 随后孔洞左下角出现了初始剪切裂纹， 初始剪切裂纹向左上方延伸形成宏观破坏裂纹，孔洞 左侧出现了明显的剥落现象; 最后另一条宏观剪切裂 纹在孔洞右下方形成，导致试件发生失稳破坏． 对于 马蹄形孔洞试样，如图 3( c) 所示，初始拉伸裂纹白斑 同样先在孔洞上下端出现; 孔洞下方的转角处由于压 应力集中出现了初始剪切裂纹; 随着轴向荷载增加，剪 切裂纹继续扩展形成宏观破坏裂纹导致试样发生剪切 破坏，孔洞左侧同样出现了明显的剥落现象． 作为对比，图 4 给出了单轴压缩条件下含孔洞试 样裂纹扩展的数值模拟结果，图中黑色和红色分别代 表张拉型微裂纹和剪切型微裂纹． 如图 4( a) 所示，单 轴压缩条件下，含圆形孔洞试样首先在孔洞上边缘出 现了拉伸裂纹; 当轴压增加到 38. 5 MPa 时，孔洞下边 缘出现了另一条拉伸裂纹，拉伸裂纹沿平行于加载方 向延伸，同时孔洞左下方形成了一条远场拉伸裂纹，试 样左下边角被压碎，形成了一条宏观裂纹; 轴压增加到 峰值应力后，微观裂纹快速增加，剪切裂纹从孔洞左右 两侧沿试样对角方向延伸，导致试样发生剪切破坏，但 孔洞上下两端的拉伸裂纹没有延伸至试样外边缘． 由 图 4( b) 可知，矩形孔洞试样首先在孔洞上边缘出现了 初始拉伸裂纹; 随后，另一条拉伸裂纹出现在孔洞下边 缘，孔洞周边转角处由于压应力集中出现了少量微观 · 4971 ·

黎崇金等：含孔洞大理岩破坏特性的颗粒流分析 ·1795· 35.7 MPa 38.5 MPa 40.7MPa.蜂值 24.4 MPa 初始拉 伸塑纹 (a) 远场拉伸裂纹 34.5 MPa 37.2 MPa 39.5MPa.峰值 23.6 MPa V型破碎区 b 初始拉 伸裂纹 远场拉伸裂纹 38.5 MPa 40.2 MPa 42.1MPa,峰值 25.3 MPa 初始拉 (c) 伸裂纹 V型破碎区 远场拉伸裂纹 图4单轴压缩条件下含孔洞试样的裂纹扩展过程（数值模拟结果）.(a)圆形孔洞：(b)矩形孔洞：()马蹄形孔洞 Fig.4 Crack coalescence process of a specimen containing a single hole under uniaxial compression (numerical results):(a)circular hole:(b) rectangular hole;(c)U-shaped hole 裂纹：当轴压增加到39.5MPa时，转角处的微观裂纹 纹首先从孔洞周边产生并向试样边缘延伸，但初始拉 向试样中部扩展、搭接形成V型破碎区，孔洞左下方 伸裂纹一般不扩展至试样边缘，试样的宏观破坏主要 出现了一组沿加载方向发育的远场拉伸裂纹：随加载 由剪切裂纹引起。对比不同孔洞形状试样的破坏过程 继续进行，剪切裂纹从V型破碎区的尖端开始向试样 可以发现，在矩形孔洞和马蹄形孔洞的侧面出现了V 的对角方向延伸，并与远场拉伸裂纹组成拉剪混合破 型破碎区（室内试验结果表现为孔壁剥落），剪切裂纹 碎带，最终导致试样发生宏观剪切破坏.对于马蹄形 从V型破碎区开始向试样对角延伸：而圆形孔洞边缘 孔洞试样，如图4(c)所示，当轴压为38.5MPa时，初 未出现V型破碎区，剪切裂纹直接从孔洞左右两侧开 始拉伸裂纹首先出现在孔洞下边缘，同时孔洞下方转 始向试样对角延伸.值得注意的是，由于岩石类材料 角处出现了少量微观裂纹；随轴压增加到40.2MPa,另 的非均质性，数值模拟的裂纹扩展过程难以得到与室 一条拉伸裂纹出现在孔洞上边缘，同时孔洞左下方出 内试验完全一致的结果. 现了一组沿加载方向发育的远场拉伸裂纹：当轴压增 2.2双轴压缩条件下的破坏特性 加到峰值时，左侧的剪切裂纹从孔洞左转角处开始向 为研究含孔洞试样破坏特性的围压效应，利用 左下方延伸并与远场拉伸裂纹搭接形成拉剪混合破碎 P℉C2D模拟含孔洞试样在不同围压条件下的双轴压 带，同时孔洞右侧同样出现了V型破碎区：试样最终 缩试验，结果如图5~7所示. 发生宏观剪切破坏 图5为双轴压缩条件下偏应力、微观裂纹数与轴 对比室内试验与数值模拟结果可以发现，数值模 向应变的关系曲线.从图中可以看出，虽然试样预制 拟的裂纹扩展过程与室内试验所得的裂纹扩展过程非 孔洞形状不同，但偏应力和微观裂纹数的变化规律基 常接近，说明P℉C2D能够较好地模拟含孔洞岩石试样 本相同.以圆形孔洞试样为例，加载初期试样的微观 的渐进破坏过程.单轴压缩条件下，含孔洞试样的裂 裂纹累积数为零，偏应力呈直线增加：当轴向应变增加

黎崇金等: 含孔洞大理岩破坏特性的颗粒流分析 图 4 单轴压缩条件下含孔洞试样的裂纹扩展过程 ( 数值模拟结果) ． ( a) 圆形孔洞; ( b) 矩形孔洞; ( c) 马蹄形孔洞 Fig． 4 Crack coalescence process of a specimen containing a single hole under uniaxial compression ( numerical results) : ( a) circular hole; ( b) rectangular hole; ( c) U-shaped hole 裂纹; 当轴压增加到 39. 5 MPa 时，转角处的微观裂纹 向试样中部扩展、搭接形成 V 型破碎区，孔洞左下方 出现了一组沿加载方向发育的远场拉伸裂纹; 随加载 继续进行，剪切裂纹从 V 型破碎区的尖端开始向试样 的对角方向延伸，并与远场拉伸裂纹组成拉剪混合破 碎带，最终导致试样发生宏观剪切破坏． 对于马蹄形 孔洞试样，如图 4( c) 所示，当轴压为 38. 5 MPa 时，初 始拉伸裂纹首先出现在孔洞下边缘，同时孔洞下方转 角处出现了少量微观裂纹; 随轴压增加到 40. 2 MPa，另 一条拉伸裂纹出现在孔洞上边缘，同时孔洞左下方出 现了一组沿加载方向发育的远场拉伸裂纹; 当轴压增 加到峰值时，左侧的剪切裂纹从孔洞左转角处开始向 左下方延伸并与远场拉伸裂纹搭接形成拉剪混合破碎 带，同时孔洞右侧同样出现了 V 型破碎区; 试样最终 发生宏观剪切破坏． 对比室内试验与数值模拟结果可以发现，数值模 拟的裂纹扩展过程与室内试验所得的裂纹扩展过程非 常接近，说明 PFC2D 能够较好地模拟含孔洞岩石试样 的渐进破坏过程． 单轴压缩条件下，含孔洞试样的裂 纹首先从孔洞周边产生并向试样边缘延伸，但初始拉 伸裂纹一般不扩展至试样边缘，试样的宏观破坏主要 由剪切裂纹引起． 对比不同孔洞形状试样的破坏过程 可以发现，在矩形孔洞和马蹄形孔洞的侧面出现了 V 型破碎区( 室内试验结果表现为孔壁剥落) ，剪切裂纹 从 V 型破碎区开始向试样对角延伸; 而圆形孔洞边缘 未出现 V 型破碎区，剪切裂纹直接从孔洞左右两侧开 始向试样对角延伸． 值得注意的是，由于岩石类材料 的非均质性，数值模拟的裂纹扩展过程难以得到与室 内试验完全一致的结果． 2. 2 双轴压缩条件下的破坏特性 为研究含孔洞试样破坏特性的围压效应，利 用 PFC2D 模拟含孔洞试样在不同围压条件下的双轴压 缩试验，结果如图 5 ～ 7 所示． 图 5 为双轴压缩条件下偏应力、微观裂纹数与轴 向应变的关系曲线． 从图中可以看出，虽然试样预制 孔洞形状不同，但偏应力和微观裂纹数的变化规律基 本相同． 以圆形孔洞试样为例，加载初期试样的微观 裂纹累积数为零，偏应力呈直线增加; 当轴向应变增加 · 5971 ·