循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗粒流模拟

工程科学学报 Chinese Journal of Engineering 循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗粒流棋拟 张杰郭奇峰蔡美峰张英汪炳锋吴星辉 Particle flow simulation of the crack propagation characteristics of granite under cyclic load ZHANG Jie.GUO Qi-feng.CAI Mei-feng.ZHANG Ying,WANG Bing-feng.WU Xing-hui 引用本文： 张杰，郭奇峰，蔡美峰，张英，汪炳锋，吴星辉.循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗粒流模拟).工程科学 学报，2021,43(5)：636-646.doi:10.13374.issn2095-9389.2020.03.15.003 ZHANG Jie,GUO Qi-feng,CAI Mei-feng.ZHANG Ying,WANG Bing-feng,WU Xing-hui.Particle flow simulation of the crack propagation characteristics of granite under cyclic load[J].Chinese Journal of Engineering,2021,43(5):636-646.doi: 10.13374-issn2095-9389.2020.03.15.003 在线阅读View online:https::/doi.org10.13374.issn2095-9389.2020.03.15.003 您可能感兴趣的其他文章 Articles you may be interested in 双孔聚能爆破煤层裂隙扩展贯通机理 Crack propagation and coalescence mechanism of double-hole cumulative blasting in coal seam 工程科学学报.2020,42(12：1613 https:ldoi.org/10.13374.issn2095-9389.2020.05.19.001 循环动力扰动下花岗岩细观损伤特性试验研究 Experimental study on microscopic damage characteristics of granite under cyclic dynamic disturbances 工程科学学报.2017,3911)：1634htps:1doi.org10.13374.issn2095-9389.2017.11.004 含孔洞大理岩破坏特性的颗粒流分析 Particle flow analysis of fracture characteristics of marble with a single hole 工程科学学报.2017,3912：1791 https:/1doi.org10.13374j.issn2095-9389.2017.12.003 循环荷载下花岗岩应力门槛值的细观能量演化及岩爆倾向性 Meso-energy evolution and rock burst proneness of the stress thresholds of granite under triaxial cyclic loading and unloading test 工程科学学报.2019,41(7)：864 https::/1oi.org/10.13374.issn2095-9389.2019.07.004 GH4169合金高温疲劳裂纹扩展的微观损伤机制 Microscopic damage mechanisms during fatigue crack propagation at high temperature in GH4169 superalloy 工程科学学报.2018,40(7)：822 https:ldoi.org10.13374j.issn2095-9389.2018.07.008 煤层深孔聚能爆破有效致裂范围探讨 Effective fracture zone under deep-hole cumulative blasting in coal seam 工程科学学报.2019,41（⑤）：582htps:/1doi.org/10.13374.issn2095-9389.2019.05.004

循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗粒流模拟 张杰 郭奇峰 蔡美峰 张英 汪炳锋 吴星辉 Particle flow simulation of the crack propagation characteristics of granite under cyclic load ZHANG Jie, GUO Qi-feng, CAI Mei-feng, ZHANG Ying, WANG Bing-feng, WU Xing-hui 引用本文: 张杰, 郭奇峰, 蔡美峰, 张英, 汪炳锋, 吴星辉. 循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗粒流模拟[J]. 工程科学 学报, 2021, 43(5): 636-646. doi: 10.13374/j.issn2095-9389.2020.03.15.003 ZHANG Jie, GUO Qi-feng, CAI Mei-feng, ZHANG Ying, WANG Bing-feng, WU Xing-hui. Particle flow simulation of the crack propagation characteristics of granite under cyclic load[J]. Chinese Journal of Engineering, 2021, 43(5): 636-646. doi: 10.13374/j.issn2095-9389.2020.03.15.003 在线阅读 View online: https://doi.org/10.13374/j.issn2095-9389.2020.03.15.003 您可能感兴趣的其他文章 Articles you may be interested in 双孔聚能爆破煤层裂隙扩展贯通机理 Crack propagation and coalescence mechanism of double-hole cumulative blasting in coal seam 工程科学学报. 2020, 42(12): 1613 https://doi.org/10.13374/j.issn2095-9389.2020.05.19.001 循环动力扰动下花岗岩细观损伤特性试验研究 Experimental study on microscopic damage characteristics of granite under cyclic dynamic disturbances 工程科学学报. 2017, 39(11): 1634 https://doi.org/10.13374/j.issn2095-9389.2017.11.004 含孔洞大理岩破坏特性的颗粒流分析 Particle flow analysis of fracture characteristics of marble with a single hole 工程科学学报. 2017, 39(12): 1791 https://doi.org/10.13374/j.issn2095-9389.2017.12.003 循环荷载下花岗岩应力门槛值的细观能量演化及岩爆倾向性 Meso-energy evolution and rock burst proneness of the stress thresholds of granite under triaxial cyclic loading and unloading test 工程科学学报. 2019, 41(7): 864 https://doi.org/10.13374/j.issn2095-9389.2019.07.004 GH4169合金高温疲劳裂纹扩展的微观损伤机制 Microscopic damage mechanisms during fatigue crack propagation at high temperature in GH4169 superalloy 工程科学学报. 2018, 40(7): 822 https://doi.org/10.13374/j.issn2095-9389.2018.07.008 煤层深孔聚能爆破有效致裂范围探讨 Effective fracture zone under deep-hole cumulative blasting in coal seam 工程科学学报. 2019, 41(5): 582 https://doi.org/10.13374/j.issn2095-9389.2019.05.004

工程科学学报.第43卷.第5期：636-646.2021年5月 Chinese Journal of Engineering,Vol.43,No.5:636-646,May 2021 https://doi.org/10.13374/j.issn2095-9389.2020.03.15.003;http://cje.ustb.edu.cn 循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗 粒流模拟 张杰2，郭奇峰1,2)区，蔡美峰，2，张英12)，汪炳锋1,2，吴星辉2 1)北京科技大学土木与资源工程学院.北京1000832)北京科技大学金属矿山高效开采与安全教育部重点实验室，北京100083 ☒通信作者，E-mail:guoqifeng@ustb.edu.cn 摘要从细观角度、采用颗粒离散元法开展了预制裂隙花岗岩循环加卸载的数值模拟试验.首先，使用图像处理技术识别 花岗岩中的不同细观组分、结合室内单轴压缩试验结果对细观力学参数进行了标定.然后，通过编制颗粒流代码追踪裂隙的 类型和扩展过程，分析岩石破坏过程中裂隙发展的阶段性特征.结果表明：不同倾角裂隙岩石的新生裂隙走向与预制裂隙贯 通方向基本一致；根据新生裂隙的优势倾向分组得到裂隙起裂角与预制裂隙倾角的关系：倾角B≤45时剪切和张拉裂隙的起 裂角单调递减，倾角B≥60时剪切和张拉裂隙的起裂角单调递增：循环扰动荷载增加了裂隙岩体的轴向变形，轴向累积残余 应变曲线呈反S形、提高扰动荷载应力上限促使曲线进入加速阶段：试件峰值强度随裂隙倾角增大表现出先减小后增大的 趋势，峰值强度为实验室完整岩石单轴抗压强度的63%~89%，反映了较为明显的劣化现象：在循环荷载作用下，剪切裂隙和 张拉裂隙增长曲线表现出明显的变化特点，在裂隙不稳定扩展阶段中张拉裂隙数目增长速率显著大于剪切裂隙，对分析岩石 变形破坏过程具有一定的参考意义 关键词细观结构：循环扰动：裂隙产状；裂隙扩展：颗粒流 分类号TD315.3 Particle flow simulation of the crack propagation characteristics of granite under cyclic load ZHANG Jie2)GUO Qi-feng CAl Mei-feng2,ZHANG Ying2 WANG Bing-feng2,WU Xing-hui2 1)School of Civil and Resource Engineering,University of Science and Technology Beijing,Beijing 100083,China 2)Key Laboratory of Ministry of Education of China for High-Efficient Mining and Safety of Metal Mines,University of Science and Technology Beijing,Beijing 100083,China Corresponding author,E-mail:guogifeng @ustb.edu.cn ABSTRACT The microcracks in natural rock masses considerably impact the stability of the underground engineering structures.The mechanical properties of the cracked rock masses contribute considerably to the strength of the rock masses and their compression failure mechanism.The instability and failure of the surrounding rocks are often induced by the propagation and penetration of these internal cracks.In practical engineering,rock mass excavation is a process involving dynamic disturbance.The mechanical properties of the rocks under cyclic load are considerably different from those of the rocks under static load.The characteristics and development of microcracks are the main factors influencing rock fatigue failure.From the microscopic viewpoint,the particle-based discrete element method is used to conduct the cyclic loading and unloading tests of the preexisting cracked granite.First,the microcompositions of granite are determined using image processing techniques,and the micromechanical parameters are calibrated based on the indoor 收稿日期：2020-03-15 基金项目：中央高校基本科研业务费资助项目(FRF-TP.18-015A3):国家自然科学基金资助项目(51974014，U2034206)

循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗 粒流模拟 张    杰1,2)，郭奇峰1,2) 苣，蔡美峰1,2)，张    英1,2)，汪炳锋1,2)，吴星辉1,2) 1) 北京科技大学土木与资源工程学院，北京 100083    2) 北京科技大学金属矿山高效开采与安全教育部重点实验室，北京 100083 苣通信作者，E-mail：guoqifeng@ustb.edu.cn 摘    要    从细观角度、采用颗粒离散元法开展了预制裂隙花岗岩循环加卸载的数值模拟试验. 首先，使用图像处理技术识别 花岗岩中的不同细观组分、结合室内单轴压缩试验结果对细观力学参数进行了标定. 然后，通过编制颗粒流代码追踪裂隙的 类型和扩展过程，分析岩石破坏过程中裂隙发展的阶段性特征. 结果表明：不同倾角裂隙岩石的新生裂隙走向与预制裂隙贯 通方向基本一致；根据新生裂隙的优势倾向分组得到裂隙起裂角与预制裂隙倾角的关系：倾角 β≤45°时剪切和张拉裂隙的起 裂角单调递减，倾角 β≥60°时剪切和张拉裂隙的起裂角单调递增；循环扰动荷载增加了裂隙岩体的轴向变形，轴向累积残余 应变曲线呈反 S 形、提高扰动荷载应力上限促使曲线进入加速阶段；试件峰值强度随裂隙倾角增大表现出先减小后增大的 趋势，峰值强度为实验室完整岩石单轴抗压强度的 63% ~ 89%，反映了较为明显的劣化现象；在循环荷载作用下，剪切裂隙和 张拉裂隙增长曲线表现出明显的变化特点，在裂隙不稳定扩展阶段中张拉裂隙数目增长速率显著大于剪切裂隙，对分析岩石 变形破坏过程具有一定的参考意义. 关键词    细观结构；循环扰动；裂隙产状；裂隙扩展；颗粒流 分类号    TD315.3 Particle flow simulation of the crack propagation characteristics of granite under cyclic load ZHANG Jie1,2) ，GUO Qi-feng1,2) 苣 ，CAI Mei-feng1,2) ，ZHANG Ying1,2) ，WANG Bing-feng1,2) ，WU Xing-hui1,2) 1) School of Civil and Resource Engineering, University of Science and Technology Beijing, Beijing 100083, China 2) Key Laboratory of Ministry of Education of China for High-Efficient Mining and Safety of Metal Mines, University of Science and Technology Beijing, Beijing 100083, China 苣 Corresponding author, E-mail: guoqifeng@ustb.edu.cn ABSTRACT    The microcracks in natural rock masses considerably impact the stability of the underground engineering structures. The mechanical properties of the cracked rock masses contribute considerably to the strength of the rock masses and their compression failure mechanism. The instability and failure of the surrounding rocks are often induced by the propagation and penetration of these internal cracks.  In  practical  engineering,  rock  mass  excavation  is  a  process  involving  dynamic  disturbance.  The  mechanical  properties  of  the rocks  under  cyclic  load  are  considerably  different  from  those  of  the  rocks  under  static  load.  The  characteristics  and  development  of microcracks are the main factors influencing rock fatigue failure. From the microscopic viewpoint, the particle-based discrete element method  is  used  to  conduct  the  cyclic  loading  and  unloading  tests  of  the  preexisting  cracked  granite.  First,  the  microcompositions  of granite  are  determined  using  image  processing  techniques,  and  the  micromechanical  parameters  are  calibrated  based  on  the  indoor 收稿日期: 2020−03−15 基金项目: 中央高校基本科研业务费资助项目（FRF-TP-18-015A3）；国家自然科学基金资助项目（51974014，U2034206） 工程科学学报，第 43 卷，第 5 期：636−646，2021 年 5 月 Chinese Journal of Engineering, Vol. 43, No. 5: 636−646, May 2021 https://doi.org/10.13374/j.issn2095-9389.2020.03.15.003; http://cje.ustb.edu.cn

张杰等：循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗粒流模拟 637. uniaxial compression test results.The stage of crack development during rock failure is analyzed by compiling particle flow code to track the type and propagation process of cracks.Results indicate that the orientations of new cracks in fractured rocks with different dip angles are similar to those of the preexisting cracks.Further,the relation between the crack initiation angle and the inclination angle of the preexisting cracks is obtained according to the tendency of new cracks.The crack initiation angle of shear and tension cracks decreases and increases monotonically,respectively,when the inclination angle B s45 and B>60.The cyclic disturbance load increases the axial deformation of the fractured rock mass,and the axial cumulative residual strain curve exhibits an inverse S-shape when entering the acceleration stage faster with the increasing upper stress limit.The peak strength of the model specimen shows a decreasing trend followed by an increasing trend with the increasing fracture inclination.The peak strengths of the laboratory-intact rock are 63%to 89%,indicating an obvious deterioration phenomenon in the rock materials.The growth of shear and tension cracks show different characteristics under cyclic load;the growth rate of tension cracks is considerably higher than that of shear cracks during the unstable crack development stage.The results presented in this study may be used as reference to investigate the deformation and failure mechanisms of rock materials. KEY WORDS microstructure;cyclic load;crack orientation;crack propagation;particle flow 天然岩体中存在的微裂隙和宏观裂纹对地下 的裂纹比单调加载下的裂纹更加明显.Ghamgosar 工程结构的稳定性有着显著的影响.裂隙岩体的 等进行了循环荷载作用下岩石破裂过程的实验 力学性质决定了岩体的强度和岩体受压破坏的机 研究，认为岩石中裂纹扩展由微裂纹的萌生与合 制-)，围岩的失稳和破坏往往是由这些内部裂隙 并导致.祝艳波等I使用MTS815力学试验系统 的扩展和贯通诱发引起为此，许多学者采用 研究了石膏质岩在不同循环应力加载条件下的微 天然岩石、类岩石材料模型和数值分析模型开展 观结构特征和疲劳寿命.何明明等研究了不同 了相关工作6，从不同角度研究了单裂隙和多裂 应力振幅条件下的砂岩变形特性，以屈服应力作 隙岩石在特定试验条件下裂隙萌生和扩展问题， 为循环加载过程中硬化和软化特性变化的分界 对裂隙岩体翼裂纹和次生裂纹的扩展形态、破坏 点.任松等20采用声发射与核磁共振分析技术研 模式和应力场的变化规律有了较深入的认识 究，发现饱水砂岩细碎裂隙较大裂隙体积之比与 岩石起裂应力研究以裂纹应变模型为基础， 循环过程相关 通常以完整岩体起裂应力定义裂隙岩体的起裂应 上述研究成果有助于正确认识裂隙岩体的破 力，Cai等山根据试验结果计算出多节理岩体起裂 坏机理，但现有循环扰动荷载作用下岩体新生裂 应力为其峰值应力的0.4~0.6倍，王桂林等21计算裂 隙的扩展特征研究还不够充分.目前，扰动荷载作 隙砂岩的起裂应力水平在0.6~0.7之间.Afolagboye 用下裂隙岩石的研究仍以室内试验为主，裂隙岩 等研究了裂隙倾角和岩桥长度对裂纹起裂位 体取样、预制裂隙精度控制和矿物颗粒黏结等问 置、扩展路径的影响.Bobet和Einstein!4探究了 题均对研究结果产生影响.采用数值分析模型进 岩体在受压过程中，翼裂纹和次生裂纹的扩展形 行研究，能够考虑岩石中细观组分的类型和分布， 态和破坏模式，赵海军等)分析了非连续结构面 便于探究微细观破裂的具体演化过程.因此，本文 对岩体力学特性的影响，认为裂隙产状与最大主 以三山岛金矿花岗岩试块为研究对象，采用图像 压应力取向夹角及围压大小决定其破裂特征.唐 处理技术建立基于花岗岩试件细观组分的颗粒流 礼忠和宋徉霖基于颗粒流离散元法，探究了裂 模型，开展预制单裂隙花岗岩循环加卸载试验，以 隙倾角、岩体抗压强度、裂隙局部应力场的变化 便揭示循环扰动应力下裂隙花岗岩中裂隙发展的 规律，认为缓倾斜裂隙更利于岩体的塑性变形、陡 阶段性特征，为地下工程岩体稳定性研究提供理 倾斜裂隙则表现出明显的弹脆性性质 论基础 在实际工程中，岩体开挖是一个动态扰动的 1颗粒流模型建立 过程.岩石在循环荷载作用与静荷载条件下的力 学性能差异显著，岩石本身特性和微裂隙发展是 1.1花岗岩组分识别 岩石疲劳破坏的主要因素)近些年，裂隙岩体的 岩石模型试样在宏观压缩和拉伸条件下的力 疲劳损伤研究逐渐增多.Chen等I6比较了循环加 学性能与颗粒的细观力学参数有关，同时也受矿 载和单调加载下的裂纹扩展特征，发现循环加载下 物组分的影响.花岗岩由石英、长石和云母等不

uniaxial compression test results. The stage of crack development during rock failure is analyzed by compiling particle flow code to track the  type  and  propagation  process  of  cracks.  Results  indicate  that  the  orientations  of  new  cracks  in  fractured  rocks  with  different  dip angles are similar to those of the preexisting cracks. Further, the relation between the crack initiation angle and the inclination angle of the  preexisting  cracks  is  obtained  according  to  the  tendency  of  new  cracks.  The  crack  initiation  angle  of  shear  and  tension  cracks decreases  and  increases  monotonically,  respectively,  when  the  inclination  angle β ≤ 45°  and β ≥ 60°.  The  cyclic  disturbance  load increases the axial deformation of the fractured rock mass, and the axial cumulative residual strain curve exhibits an inverse S-shape when  entering  the  acceleration  stage  faster  with  the  increasing  upper  stress  limit.  The  peak  strength  of  the  model  specimen  shows  a decreasing trend followed by an increasing trend with the increasing fracture inclination. The peak strengths of the laboratory-intact rock are 63% to 89%, indicating an obvious deterioration phenomenon in the rock materials. The growth of shear and tension cracks show different characteristics under cyclic load; the growth rate of tension cracks is considerably higher than that of shear cracks during the unstable crack development stage. The results presented in this study may be used as reference to investigate the deformation and failure mechanisms of rock materials. KEY WORDS    microstructure；cyclic load；crack orientation；crack propagation；particle flow 天然岩体中存在的微裂隙和宏观裂纹对地下 工程结构的稳定性有着显著的影响. 裂隙岩体的 力学性质决定了岩体的强度和岩体受压破坏的机 制[1−3] ，围岩的失稳和破坏往往是由这些内部裂隙 的扩展和贯通诱发引起[4−5] . 为此，许多学者采用 天然岩石、类岩石材料模型和数值分析模型开展 了相关工作[6−9] ，从不同角度研究了单裂隙和多裂 隙岩石在特定试验条件下裂隙萌生和扩展问题， 对裂隙岩体翼裂纹和次生裂纹的扩展形态、破坏 模式和应力场的变化规律有了较深入的认识. 岩石起裂应力研究以裂纹应变模型为基础[10] ， 通常以完整岩体起裂应力定义裂隙岩体的起裂应 力，Cai 等[11] 根据试验结果计算出多节理岩体起裂 应力为其峰值应力的 0.4～0.6 倍，王桂林等[12] 计算裂 隙砂岩的起裂应力水平在 0.6～0.7 之间. Afolagboye 等[13] 研究了裂隙倾角和岩桥长度对裂纹起裂位 置、扩展路径的影响. Bobet 和 Einstein[14] 探究了 岩体在受压过程中，翼裂纹和次生裂纹的扩展形 态和破坏模式. 赵海军等[5] 分析了非连续结构面 对岩体力学特性的影响，认为裂隙产状与最大主 压应力取向夹角及围压大小决定其破裂特征. 唐 礼忠和宋徉霖[9] 基于颗粒流离散元法，探究了裂 隙倾角、岩体抗压强度、裂隙局部应力场的变化 规律，认为缓倾斜裂隙更利于岩体的塑性变形、陡 倾斜裂隙则表现出明显的弹脆性性质. 在实际工程中，岩体开挖是一个动态扰动的 过程. 岩石在循环荷载作用与静荷载条件下的力 学性能差异显著，岩石本身特性和微裂隙发展是 岩石疲劳破坏的主要因素[15] . 近些年，裂隙岩体的 疲劳损伤研究逐渐增多. Chen 等[16] 比较了循环加 载和单调加载下的裂纹扩展特征，发现循环加载下 的裂纹比单调加载下的裂纹更加明显. Ghamgosar 等[17] 进行了循环荷载作用下岩石破裂过程的实验 研究，认为岩石中裂纹扩展由微裂纹的萌生与合 并导致. 祝艳波等[18] 使用 MTS815 力学试验系统 研究了石膏质岩在不同循环应力加载条件下的微 观结构特征和疲劳寿命. 何明明等[19] 研究了不同 应力振幅条件下的砂岩变形特性，以屈服应力作 为循环加载过程中硬化和软化特性变化的分界 点. 任松等[20] 采用声发射与核磁共振分析技术研 究，发现饱水砂岩细碎裂隙-较大裂隙体积之比与 循环过程相关. 上述研究成果有助于正确认识裂隙岩体的破 坏机理，但现有循环扰动荷载作用下岩体新生裂 隙的扩展特征研究还不够充分. 目前，扰动荷载作 用下裂隙岩石的研究仍以室内试验为主，裂隙岩 体取样、预制裂隙精度控制和矿物颗粒黏结等问 题均对研究结果产生影响. 采用数值分析模型进 行研究，能够考虑岩石中细观组分的类型和分布， 便于探究微细观破裂的具体演化过程. 因此，本文 以三山岛金矿花岗岩试块为研究对象，采用图像 处理技术建立基于花岗岩试件细观组分的颗粒流 模型，开展预制单裂隙花岗岩循环加卸载试验，以 便揭示循环扰动应力下裂隙花岗岩中裂隙发展的 阶段性特征，为地下工程岩体稳定性研究提供理 论基础. 1    颗粒流模型建立 1.1    花岗岩组分识别 岩石模型试样在宏观压缩和拉伸条件下的力 学性能与颗粒的细观力学参数有关，同时也受矿 物组分的影响. 花岗岩由石英、长石和云母等不 张    杰等： 循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗粒流模拟 · 637 ·

638 工程科学学报，第43卷，第5期 同的矿物胶结而成，不同矿物若以均质材料进行 中长石、石英和云母所占图像面积百分比，分别 赋值校准很可能导致不准确的模拟结果.可以基 为49.85%、29.74%和20.41%，其中云母含量过高， 于真实岩石材料数字图像、以矿物颗粒表面不同 这与阈值选取有关.实际花岗岩矿物组分中除上 颜色或亮度变化准确识别花岗岩表面细观结构， 述三种主要矿物外，还有角闪石、少量辉石等次要 进而将所得结构作为建立离散元模型的基本条件 矿物，可将长石和石英以外的矿物百分比视为黑 试验岩块尺寸为：50mm×50mm×100m,使 云母矿物含量 用数码相机获得花岗岩试块的数字图像，图像中 的每个像素点都可以用像素点矩阵表示： Gray value 0 f(1,1)f1,2) f1,N) Feldspar f2,1) f2,2) f2,N) F(i,)= (1) Quartz Mica fM,1)f(M,2) f(M,N) 255 式中，f亿，)表示在0,255)范围内的像素点(，) 的灰度值，M和N代表像素点的行数和列数.采 图1花岗岩图像矿物识别.(a)标准试件：(b)局部放大图 Fig.I Mineral recognition from a granite image:(a)standard test 用MATLAB中图像处理方法，对原始图像进行平 specimen;(b)partial enlarged detail 滑处理以消除噪声影响.然后采用阈值分割算法， 采用图像处理技术得到的花岗岩细观组分， 识别花岗岩矿物组分特征. 并不能完全展现其表面形态特征.在二维数值计 改进的最大类间方差法适用于图像的多阈值 算模型中，通常以圆形颗粒为基础来构建矿物组 分割，以类间特性差最大或类内特性差最小为判 分的分布区域，从而建立反映试件组分分布的细 据，求得图像分割的最佳阈值.设灰度范围O, 观模型：然而，矿物组分的三维形态特征难以直接 L]的图像中，第i级像素点个数为N,像素点总个 获取，PFC3D中的球形颗粒很难用于直接建立表征 数为N=方N,则第1级像素点的出现概率为 矿物组分分布的多面体.为简化计算模型，本文根 =0 据试样二维图像，通过编程识别不同矿物组分，计 N 选择灰度阈值组T={,2,…,tm-1小，将图像 算矿物组分占比，进而建立相应的三维分析模型. 分割为m个区域，灰度区域的灰度值分别为 1.2细观参数标定 C1=[0,t,C2=[t1,2],,Cm=[tm-1,L,则方差最 在岩石微观结构研究中21-2四，通常认为矿物成 大值的最优阈值T为： 分的强度控制着裂纹的萌生，长石、石英、黑云母 (T")=arg ma ,- 矿物强度关系为1：1.17：0.42，压痕试验下长石、 (2) 石英、黑云母的断裂韧度比值为1：2.08：0.76.根 据不同矿物破坏应力阈值的裂纹演化特征，发现 =∑Pd,=∑iP:lwx:X= 3) 微裂纹在黑云母中首先产生、在石英中产生最晚) icC. ieC, 在岩石宏细观参数关系研究中以长石组分颗 其中，o(s=1,2,3,…,m)为每一灰度区域像素点占 粒的黏结强度为平行黏结力学参数基本值，石英 图像总比的比例，C,为区域灰度值的集合，入为平 和云母组分颗粒的黏结强度取基本值的1.6倍和 均灰度值，为图像的总体灰度.采用该算法进行 0.3倍，不同组分的黏结强度取相应组分的平均 双阈值求解，得到灰度图像的最佳分割阈值为 值，颗粒的最小半径Rmn取1.2mm,最大半径和最 124和168，从而对花岗岩中的长石、石英、云母图 小半径比为1.66 像进行分割.本文尚未考虑图像中不同矿物组分 模型细观力学性质参数需要根据室内试验结 之间存在灰度交叉的情况，由此导致的误差予以 果进行标定，以保证模型试样与实验试样的宏观 忽略 力学响应基本一致.本文使用AW2000型全数字 花岗岩标准试样的灰度图如图1(a)所示，岩 电液控制刚性压力试验机进行完整岩石单轴压缩 石结构特征和阈值分割如图1(b)所示.通过对灰 试验，位移控制加载速率为0.03 mm:min、加载至 度图像中各像素点进行数量统计，得到矿物成分 试样破坏.室内单轴压缩试验得到的宏观力学性

同的矿物胶结而成，不同矿物若以均质材料进行 赋值校准很可能导致不准确的模拟结果. 可以基 于真实岩石材料数字图像、以矿物颗粒表面不同 颜色或亮度变化准确识别花岗岩表面细观结构， 进而将所得结构作为建立离散元模型的基本条件. 试验岩块尺寸为：50 mm × 50 mm × 100 m，使 用数码相机获得花岗岩试块的数字图像，图像中 的每个像素点都可以用像素点矩阵表示： F(i, j) =   f(1,1) f(1,2) ··· f(1,N) f(2,1) f(2,2) ··· f(2,N) . . . . . . . . . f(M,1) f(M,2) ··· f(M,N)   （1） 式中，f (i，j) 表示在 [0, 255] 范围内的像素点 (i，j) 的灰度值，M 和 N 代表像素点的行数和列数. 采 用 MATLAB 中图像处理方法，对原始图像进行平 滑处理以消除噪声影响. 然后采用阈值分割算法， 识别花岗岩矿物组分特征. N = ∑ L i=0 Ni Pi = Ni N T = {t1,t2,··· ,tm−1} C1 = [0,t1] C2 = [t1,t2], ··· Cm = [tm−1,L] δ 2 B T ∗ 改进的最大类间方差法适用于图像的多阈值 分割，以类间特性差最大或类内特性差最小为判 据，求得图像分割的最佳阈值. 设灰度范围 [0， L] 的图像中，第 i 级像素点个数为 Ni，像素点总个 数 为 ， 则 第 i 级像素点的出现概率为 . 选择灰度阈值组 ，将图像 分 割 为 m 个区域 ，灰度区域的灰度值分别为 ， ， ，则方差 最 大值的最优阈值 为： δ 2 B (T ∗ ) = arg max   ∑m s=1 ωs(λs −λ) 2   （2） ωs = ∑ i∈Cs Pi ;λs = ∑ i∈Cs iPi/ωs ;λ = ∑ L i=0 iPi （3） ωs(s = 1,2,3,··· ,m) Cs λs λ 其中， 为每一灰度区域像素点占 图像总比的比例， 为区域灰度值的集合， 为平 均灰度值， 为图像的总体灰度. 采用该算法进行 双阈值求解 ，得到灰度图像的最佳分割阈值为 124 和 168，从而对花岗岩中的长石、石英、云母图 像进行分割. 本文尚未考虑图像中不同矿物组分 之间存在灰度交叉的情况，由此导致的误差予以 忽略. 花岗岩标准试样的灰度图如图 1（a）所示，岩 石结构特征和阈值分割如图 1（b）所示. 通过对灰 度图像中各像素点进行数量统计，得到矿物成分 中长石、石英和云母所占图像面积百分比，分别 为 49.85%、29.74% 和 20.41%，其中云母含量过高， 这与阈值选取有关. 实际花岗岩矿物组分中除上 述三种主要矿物外，还有角闪石、少量辉石等次要 矿物，可将长石和石英以外的矿物百分比视为黑 云母矿物含量. (a) (b) Feldspar Mica Quartz 255 0 Gray value 图 1    花岗岩图像矿物识别. （a）标准试件；（b）局部放大图 Fig.1     Mineral  recognition  from  a  granite  image:  (a)  standard  test specimen; (b) partial enlarged detail 采用图像处理技术得到的花岗岩细观组分， 并不能完全展现其表面形态特征. 在二维数值计 算模型中，通常以圆形颗粒为基础来构建矿物组 分的分布区域，从而建立反映试件组分分布的细 观模型；然而，矿物组分的三维形态特征难以直接 获取，PFC3D 中的球形颗粒很难用于直接建立表征 矿物组分分布的多面体. 为简化计算模型，本文根 据试样二维图像，通过编程识别不同矿物组分，计 算矿物组分占比，进而建立相应的三维分析模型. 1.2    细观参数标定 在岩石微观结构研究中[21−22] ，通常认为矿物成 分的强度控制着裂纹的萌生，长石、石英、黑云母 矿物强度关系为 1∶1.17∶0.42，压痕试验下长石、 石英、黑云母的断裂韧度比值为 1∶2.08∶0.76. 根 据不同矿物破坏应力阈值的裂纹演化特征，发现 微裂纹在黑云母中首先产生、在石英中产生最晚[23] . 在岩石宏细观参数关系研究中[24] ，以长石组分颗 粒的黏结强度为平行黏结力学参数基本值，石英 和云母组分颗粒的黏结强度取基本值的 1.6 倍和 0.3 倍，不同组分的黏结强度取相应组分的平均 值，颗粒的最小半径 Rmin 取 1.2 mm，最大半径和最 小半径比为 1.66. 模型细观力学性质参数需要根据室内试验结 果进行标定，以保证模型试样与实验试样的宏观 力学响应基本一致. 本文使用 AW2000 型全数字 电液控制刚性压力试验机进行完整岩石单轴压缩 试验，位移控制加载速率为 0.03 mm·min−1、加载至 试样破坏. 室内单轴压缩试验得到的宏观力学性 · 638 · 工程科学学报，第 43 卷，第 5 期

张杰等：循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗粒流模拟 639. 质参数为：弹性模量E=43.12GPa,泊松比=0.163， 形颗粒刚度比值以匹配泊松比；再设侧限压力为 峰值应力o=126.83MPa. 0,逐次降低平行黏接强度，得到岩石试件的峰值 以邓树新等、Si等P的成果作为细观力学 应力；最后改变球形颗粒接触模量和平行黏结的 性质参数取值的基础，平行黏结强度比值取为 弹性模量比值，微调平行黏结的弹性模量，得到与 1.0,将试样强度作为峰值应力.先将球形颗粒和 室内试验更为匹配的应力-应变关系.得到的细观 平行黏结的接触模量比值设为1.0，通过调节颗粒 力学性质参数见表1，室内试验与模拟试验应力- 接触模量得到岩石试样的弹性模量，然后改变球 应变全过程曲线如图2所示 表1花岗岩细观力学性质参数 Table 1 Microscale mechanical parameters of granite Particles forming grains Linear parallel bond model Mineral Ratio of normal Ratio of normal Minimum particle Maximum to Young's Young's to shear Shear bond Tensile-shear component minimum modulus, to shear Friction radius forming grain,Rmin/mm radius ratio, stiffness of the modulus, stiffness of the strength, bond strength ratio, RmayRmin EGPa particle,k/k E/GPa parallel bond, Te/MPa ratio,dc/斤e kn/ks Feldspar 1.2 1.66 45.5 1.15 28.0 1.6 51.0 0.5 1.0 Quartz 1.2 1.66 33.0 1.15 22.6 1.6 81.6 0.8 1.0 Mica 1.2 1.66 11.2 1.15 5.9 1.6 15.3 0.15 1.0 140 120 -Laboratory test curve ☐Feldspar Simulation test curve Mica 100 Quart edW/ssan 80 Newly-generated cracks 员 60 40 20 0 Pre-existing cracks 0.1 0.2 0.3 0.4 Strain/% 图3裂隙花岗岩试件模型 图2试件应力-应变曲线 Fig.2 Stress-strain curves of a specimen Fig.3 Numerical model of a granite specimen with a single crack 1.3试验方案 石强度降低、模量减小，但变化规律基本一致.因 在数值模型中预制单裂隙位置删除相应的球 此，模拟试验中计算模型采用边界墙体进行位移 体颗粒.裂隙位于试样中心，长度20mm,宽度0.3mm, 加载控制，通过调整上下墙的位移给模型施加轴 与水平方向夹角B,试样模型如图3所示.预制单 向荷载，将循环加卸载速率与静力加载阶段速率 裂隙夹角B取值0°、30°、45°、60°、90°，新生裂隙 取为一致 的起裂角为0，新生裂隙倾角为α.以模型试件顶 岩体开挖产生的应力重分布和应力集中，导 部为正北方向，左右两侧方向分别为270°和90°， 致很大范围内的岩体处于静力屈服阶段.为分析 则预制裂隙的倾角为B、倾向均为90° 扰动荷载上限对花岗岩力学性质产生的影响，荷 矿山开采爆破和岩体破坏产生的低频波频率 载上限至少应大于屈服应力.本文完整花岗岩模 为几赫兹至几十赫兹2可，这种低频扰动荷载是工 型试件单轴峰值荷载为125.30MPa,屈服应力为 程岩体的一种重要受力形式.室内试验研究能够 峰值荷载的86%.如图2所示：考虑本次试验主要 以试验机进行的循环加卸载试验来模拟低频波对 研究扰动应力对岩石裂隙萌生及扩展的影响，扰 岩石的扰动作用.不同循环频率的动力加载试验 动荷载下限设置为起裂应力范围.因此，模型中加 表明，随着施加动力载荷时循环频率的增高，岩 载轴向应力采用分级增加荷载形式进行循环加卸

质参数为：弹性模量 E =43.12 GPa，泊松比 μ= 0.163， 峰值应力 σ = 126.83 MPa. 以邓树新等[25]、Shi 等[26] 的成果作为细观力学 性质参数取值的基础 ，平行黏结强度比值取为 1.0，将试样强度作为峰值应力. 先将球形颗粒和 平行黏结的接触模量比值设为 1.0，通过调节颗粒 接触模量得到岩石试样的弹性模量，然后改变球 形颗粒刚度比值以匹配泊松比；再设侧限压力为 0，逐次降低平行黏接强度，得到岩石试件的峰值 应力；最后改变球形颗粒接触模量和平行黏结的 弹性模量比值，微调平行黏结的弹性模量，得到与 室内试验更为匹配的应力−应变关系. 得到的细观 力学性质参数见表 1，室内试验与模拟试验应力− 应变全过程曲线如图 2 所示. 表 1 花岗岩细观力学性质参数 Table 1   Microscale mechanical parameters of granite Mineral component Particles forming grains Linear parallel bond model Minimum particle radius forming grain, Rmin/mm Maximum to minimum radius ratio, Rmax/Rmin Young’s modulus, Ec /GPa Ratio of normal to shear stiffness of the particle, kn /ks E¯ c Young’s modulus, /GPa ¯kn ¯ks Ratio of normal to shear stiffness of the parallel bond, / τc Shear bond strength, /MPa φ¯ Friction ratio, σ¯ c τ¯c Tensile–shear bond strength ratio, / Feldspar 1.2 1.66 45.5 1.15 28.0 1.6 51.0 0.5 1.0 Quartz 1.2 1.66 33.0 1.15 22.6 1.6 81.6 0.8 1.0 Mica 1.2 1.66 11.2 1.15 5.9 1.6 15.3 0.15 1.0 0 0.1 0.2 0.3 0.4 0 20 40 60 80 100 120 140 Laboratory test curve Simulation test curve Stress/MPa Strain/% 图 2    试件应力–应变曲线 Fig.2    Stress–strain curves of a specimen 1.3    试验方案 在数值模型中预制单裂隙位置删除相应的球 体颗粒. 裂隙位于试样中心，长度20 mm，宽度0.3 mm， 与水平方向夹角 β，试样模型如图 3 所示. 预制单 裂隙夹角 β 取值 0°、30°、45°、60°、90°，新生裂隙 的起裂角为 θ，新生裂隙倾角为 α. 以模型试件顶 部为正北方向，左右两侧方向分别为 270°和 90°， 则预制裂隙的倾角为 β、倾向均为 90°. 矿山开采爆破和岩体破坏产生的低频波频率 为几赫兹至几十赫兹[27] ，这种低频扰动荷载是工 程岩体的一种重要受力形式. 室内试验研究能够 以试验机进行的循环加卸载试验来模拟低频波对 岩石的扰动作用. 不同循环频率的动力加载试验 表明[28] ，随着施加动力载荷时循环频率的增高，岩 石强度降低、模量减小，但变化规律基本一致. 因 此，模拟试验中计算模型采用边界墙体进行位移 加载控制，通过调整上下墙的位移给模型施加轴 向荷载，将循环加卸载速率与静力加载阶段速率 取为一致. 岩体开挖产生的应力重分布和应力集中，导 致很大范围内的岩体处于静力屈服阶段. 为分析 扰动荷载上限对花岗岩力学性质产生的影响，荷 载上限至少应大于屈服应力. 本文完整花岗岩模 型试件单轴峰值荷载为 125.30 MPa，屈服应力为 峰值荷载的 86%，如图 2 所示；考虑本次试验主要 研究扰动应力对岩石裂隙萌生及扩展的影响，扰 动荷载下限设置为起裂应力范围. 因此，模型中加 载轴向应力采用分级增加荷载形式进行循环加卸 Mica β θ Feldspar Quart α Newly-generated cracks 20 mm Pre-existing cracks 图 3    裂隙花岗岩试件模型 Fig.3    Numerical model of a granite specimen with a single crack 张    杰等： 循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗粒流模拟 · 639 ·

640 工程科学学报，第43卷，第5期 载试验.通过单轴抗压强度试验，获得5组不同裂 裂隙分组统计采用相似归并法.这一方法原 隙倾角花岗岩的峰值荷载；循环荷载上限分别取 理是：归并产状相近的两条裂隙，以平均产状为新 峰值荷载的60%、70%、80%和90%，循环荷载下 的裂隙组，总裂隙数量相应减少一组，如此往复进 限为峰值荷载的40%：待到循环荷载结束后，继续 行.裂隙所在空间平面以其法线矢量表示，根据产 加载到试样破坏.在研究扰动荷载对岩石裂隙萌 状可换算出其平面法线矢量余弦（即单位矢量）： 生的影响分析中-，该循环荷载下限值大于试 {sinBsin叱，sinBcosx,cosB引，X为裂隙倾向，B为裂隙倾 件压密阶段的应力，可以避免压密阶段对循环过 角.视法线矢量夹角最小的两个平面裂隙为相近 程中花岗岩变形特性的影响.每级荷载循环六次， 产状（采用法线矢量夹角的余弦值进行比较），计 以模拟扰动压应力的环境 算如下： 2裂隙岩体破裂过程分析 cosΦ=sinp1sinX1sinf2sinK2+sinβ1 cosxI sinB2cos2+cosB cosB2 (4) 2.1新生裂隙分布规律 则平均产状对应的法线矢量为： 2.11裂隙萌生走向分布 (KYZ={∑singsin此am-∑sinB.cos,∑cosB 岩石中的局部应力和集中应力促进了新生裂 (5) 隙的萌生与扩展，裂隙主要包括剪切裂隙和拉伸 式中，X轴表示平均产状对应正北方向的法线矢量， 裂隙两种类型.通过内置FISH语言编制跟踪颗粒 Y轴为正东方向的法线矢量，Z轴垂直水平面的法 接触点破坏情况的程序，确定模型中新生剪切和 线矢量，n为平面裂隙组数.平均产状可按下式计算： 拉伸裂隙的位置、类型和裂隙面产状.裂隙花岗 X=arctan(lX/IYD） (6) 岩在实验条件下，起裂强度为峰值强度的0.4~ 0.8倍2本次研究中，峰值荷载前应力水平达到 B=arccos(IZI/X+2) (7) 0.4时，花岗岩模型萌生少量的裂隙：当第一级循 由此得到裂隙岩石试件的走向玫瑰花图如图4 环荷载（应力水平0.4~0.6）结束后，对模型中预制 所示.裂隙倾向和倾角统计结果见表2，其中裂隙 裂隙周边新生裂隙情况进行统计分析 的优势倾向分组数取2 (a) Shear crack (b) Shear crack Tension crack 270° 90°270° 90°270° 90°270° Sheaerack Tension crack (d) Shear crack Tension crack 09 00 08 270° 90°270° 90°270° 90270° (e) Shear crack Tension crack 0° 270° 90°270° 90° 图4裂隙岩石试件走向政瑰花图.(a)B=0°：(b)B=30:(cB=45°：(d)B=60:(e)B=90° Fig.4 Strike rose diagrams of a cracked rock specimen:(a)=(b)B=30(c)B=45(d)60:(e)B=9 根据统计结果，当裂隙倾角=0时张拉裂隙 和剪切裂隙优势走向在0°左右，与预制裂隙走向 的优势走向在331°~350°；当裂隙倾角=30°和 一致；当裂隙倾角=90时剪切裂隙多在351°~ =45时，新生剪切裂隙的优势走向分布于351°~ 360° 360°，张拉裂隙的优势走向在341°~350°，两类裂 由于试样模型中预制裂隙贯通，在单轴压缩 隙分布较为集中；当裂隙倾角=60时，张拉裂隙 荷载作用下新生裂隙的走向大体与预制裂隙一

载试验. 通过单轴抗压强度试验，获得 5 组不同裂 隙倾角花岗岩的峰值荷载；循环荷载上限分别取 峰值荷载的 60%、70%、80% 和 90%，循环荷载下 限为峰值荷载的 40%；待到循环荷载结束后，继续 加载到试样破坏. 在研究扰动荷载对岩石裂隙萌 生的影响分析中[11−12] ，该循环荷载下限值大于试 件压密阶段的应力，可以避免压密阶段对循环过 程中花岗岩变形特性的影响. 每级荷载循环六次， 以模拟扰动压应力的环境. 2    裂隙岩体破裂过程分析 2.1    新生裂隙分布规律 2.1.1    裂隙萌生走向分布 岩石中的局部应力和集中应力促进了新生裂 隙的萌生与扩展，裂隙主要包括剪切裂隙和拉伸 裂隙两种类型. 通过内置 FISH 语言编制跟踪颗粒 接触点破坏情况的程序，确定模型中新生剪切和 拉伸裂隙的位置、类型和裂隙面产状. 裂隙花岗 岩在实验条件下，起裂强度为峰值强度的 0.4～ 0.8 倍[29] . 本次研究中，峰值荷载前应力水平达到 0.4 时，花岗岩模型萌生少量的裂隙；当第一级循 环荷载（应力水平 0.4～0.6）结束后，对模型中预制 裂隙周边新生裂隙情况进行统计分析. {sinβsinχ,sinβcosχ, cosβ} χ β 裂隙分组统计采用相似归并法. 这一方法原 理是：归并产状相近的两条裂隙，以平均产状为新 的裂隙组，总裂隙数量相应减少一组，如此往复进 行. 裂隙所在空间平面以其法线矢量表示，根据产 状可换算出其平面法线矢量余弦（即单位矢量）： ， 为裂隙倾向， 为裂隙倾 角. 视法线矢量夹角最小的两个平面裂隙为相近 产状（采用法线矢量夹角的余弦值进行比较），计 算如下： cosΦ = sinβ1 sinχ1 sinβ2 sinχ2+sinβ1 cosχ1 sinβ2 cosχ2+cosβ1 cosβ2 （4） 则平均产状对应的法线矢量为： {XYZ} = {∑ sinβnsinχn, ∑ sinβncosχn, ∑ cosβn } （5） 式中，X 轴表示平均产状对应正北方向的法线矢量， Y 轴为正东方向的法线矢量，Z 轴垂直水平面的法 线矢量，n 为平面裂隙组数. 平均产状可按下式计算： χ = arctan(|X| / |Y|) （6） β = arccos(|Z| / √ |X| 2 +|Y| 2 ) （7） 由此得到裂隙岩石试件的走向玫瑰花图如图 4 所示. 裂隙倾向和倾角统计结果见表 2，其中裂隙 的优势倾向分组数取 2. (a) Shear crack Tension crack 0° 270° 0° 270° 90° (b) Shear crack Tension crack 0° 270° 90° 0° 270° 90° (c) Shear crack Tension crack 0° 270° 90° 0° 270° 90° (d) Shear crack Tension crack 0° 270° 90° 0° 270° 90° (e) Shear crack Tension crack 0° 90° 270° 90° 0° 90° 270° 图 4    裂隙岩石试件走向玫瑰花图. （a）β = 0°；（b）β = 30°；（c）β = 45°；（d）β = 60°；（e）β = 90° Fig.4    Strike rose diagrams of a cracked rock specimen: (a) β = 0°; (b) β = 30°; (c) β = 45°; (d) β = 60°; (e) β = 90° 根据统计结果，当裂隙倾角 β=0°时张拉裂隙 的优势走向 在 331°～ 350°；当裂隙倾 角 β=30°和 β=45°时，新生剪切裂隙的优势走向分布于 351°～ 360°，张拉裂隙的优势走向在 341°～350°，两类裂 隙分布较为集中；当裂隙倾角 β=60°时，张拉裂隙 和剪切裂隙优势走向在 0°左右，与预制裂隙走向 一致；当裂隙倾角 β=90°时剪切裂隙多在 351°～ 360°. 由于试样模型中预制裂隙贯通，在单轴压缩 荷载作用下新生裂隙的走向大体与预制裂隙一 · 640 · 工程科学学报，第 43 卷，第 5 期

张杰等：循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗粒流模拟 641 表2新生裂隙倾向和倾角分布统计 Table 2 Statistics of the distribution of tendencies and inclinations for newly generated cracks Tendencies and Shear cracks Tension crack inclinations for preexisting cracks Tendency grouping Average inclination/()Percentage/% Tendency grouping Average inclination/()Percentage/% 151°-160° 6 63 61°-70° 77 8.0 90°∠0° 211°-220 60 6.3 251°-260° 62 8.0 261°-270° 62 11.9 71°-80° 61 10.5 90°∠30° 241°-250° 好 9.5 141°-150 72 10.5 261°-270° 4 17.4 71°-80° 52 18.5 90°∠45° 251°-260° 45 15.2 91°-100° 64 11.1 261°-270° 51 9.5 81°-90° 64 14.3 90°∠60° 251°-260° 7.9 61°-70° 69 10.7 261°-270° 45 8.5 131°-140° 88 11.1 90°∠90° 181°-190° 6.4 21°-30° 48 5.6 致，仅裂隙倾角=90时张拉裂隙优势走向在 70k 41°~50°，而裂隙倾角=0时的剪切裂隙走向分布 -Shear crack 60 -。Tension crack 较为分散.模型试件在裂隙萌生阶段首先出现剪 50 切裂隙；张拉裂隙出现较晚，出现于峰前峰值强度 40 的0.6倍处，仅占裂隙总数的6%~10%. 2.1.2裂隙起裂角 20 在裂隙萌生阶段，可根据新生裂隙的倾角分 10 布进行裂隙花岗岩的起裂角研究.朱维申等0提 0 出，起裂角度与裂隙的空间位置和加载方向有关， 0 102030405060708090 双轴压缩荷载下雁形裂纹的起裂角为70°~80°： Crack dip,/( 李银平等发现，单轴压缩荷载作用下大理岩的 图5预制裂隙倾角B与起裂角0的关系 起裂角为52°~68°.裂隙的起裂、扩展和贯通是压 Fig.5 Relation between the crack initial angle 0 and crack dip B 剪应力场中岩石的主要破坏形式，翼裂纹通常产 在单轴加载条件下更倾向于最大主应力方向，因 生于原生裂隙端部、以弧形向前扩展、偏向主应 此张拉裂隙发展最快.同理，裂隙倾角B≥60的 力方向发展. 岩石，剪切裂隙更倾向于最大主应力方向，发展最快 根据表2中岩石试样新生裂隙的优势倾向分 2.2裂隙扩展及破裂的阶段性特征 组，取两组倾角平均值，得到岩石试样的新生裂隙 岩石变形破坏过程中，裂隙经历了压密、起 倾角α，则起裂角0为新生裂隙倾角a与预制裂隙 裂、扩展、贯通四个阶段.通过监测新生裂隙出 倾角B的差值.由此得到剪切裂隙和张拉裂隙的 现时间及其数目增长曲线拐点，可以进行岩石变 起裂角度与裂隙倾角的变化情况，如图5所示. 形破坏阶段特征分析.试件模型新生裂隙数目与 从图5可以看出，起裂角随预制裂隙倾角的增 轴向应变的变化情况如图6所示. 大呈现非线性变化：倾角B≤45时，剪切和张拉裂 岩石裂隙在压密阶段，应变随应力的增加呈 隙的起裂角单调递减，且张拉裂隙比剪切裂隙的 线性增长，未出现新生裂隙.随着裂隙进一步发 起裂角大；倾角B≥60时，剪切和张拉裂隙的起裂 展，岩石压缩进入AB阶段、岩石内部开始出现剪 角单调递增，剪切裂隙的起裂角较大.新生裂隙在 切裂隙，裂隙数目极少并保持稳定；在BC阶段，随 倾角B=45°和B=60°的岩石试样中起裂角最小， 荷载增加，剪切裂隙少量增长，开始出现张拉裂 在倾角B=0°的岩石试样中起裂角最大. 隙.微裂隙在外荷载作用下压密闭合，剪切裂隙和 由于裂隙的扩展方向总是与最大拉应力方向 张拉裂隙萌生，新生裂隙出现的时间与预制裂隙 保持垂直，当裂隙扩展角保持不变时，裂隙倾角 倾角无关.由于应力水平不足以使得原有裂纹扩 B≤45的岩石，张拉裂隙比剪切裂隙的起裂角大， 展，新萌生的张拉裂隙和剪切裂隙数目较为稳定

致 ，仅裂隙倾 角 β=90°时张拉裂隙优势走向 在 41°～50°，而裂隙倾角 β=0°时的剪切裂隙走向分布 较为分散. 模型试件在裂隙萌生阶段首先出现剪 切裂隙；张拉裂隙出现较晚，出现于峰前峰值强度 的 0.6 倍处，仅占裂隙总数的 6% ～10%. 2.1.2    裂隙起裂角 在裂隙萌生阶段，可根据新生裂隙的倾角分 布进行裂隙花岗岩的起裂角研究. 朱维申等[30] 提 出，起裂角度与裂隙的空间位置和加载方向有关， 双轴压缩荷载下雁形裂纹的起裂角为 70°～80°； 李银平等[31] 发现，单轴压缩荷载作用下大理岩的 起裂角为 52°～68°. 裂隙的起裂、扩展和贯通是压 剪应力场中岩石的主要破坏形式，翼裂纹通常产 生于原生裂隙端部、以弧形向前扩展、偏向主应 力方向发展. 根据表 2 中岩石试样新生裂隙的优势倾向分 组，取两组倾角平均值，得到岩石试样的新生裂隙 倾角ɑ，则起裂角 θ 为新生裂隙倾角ɑ与预制裂隙 倾角 β 的差值. 由此得到剪切裂隙和张拉裂隙的 起裂角度与裂隙倾角的变化情况，如图 5 所示. 从图 5 可以看出，起裂角随预制裂隙倾角的增 大呈现非线性变化：倾角 β ≤ 45°时，剪切和张拉裂 隙的起裂角单调递减，且张拉裂隙比剪切裂隙的 起裂角大；倾角 β ≥ 60°时，剪切和张拉裂隙的起裂 角单调递增，剪切裂隙的起裂角较大. 新生裂隙在 倾角 β = 45°和 β = 60°的岩石试样中起裂角最小， 在倾角 β = 0°的岩石试样中起裂角最大. 由于裂隙的扩展方向总是与最大拉应力方向 保持垂直，当裂隙扩展角保持不变时，裂隙倾角 β ≤ 45°的岩石，张拉裂隙比剪切裂隙的起裂角大， 在单轴加载条件下更倾向于最大主应力方向，因 此张拉裂隙发展最快. 同理，裂隙倾角 β ≥ 60°的 岩石，剪切裂隙更倾向于最大主应力方向，发展最快. 2.2    裂隙扩展及破裂的阶段性特征 岩石变形破坏过程中，裂隙经历了压密、起 裂、扩展、贯通四个阶段[32] . 通过监测新生裂隙出 现时间及其数目增长曲线拐点，可以进行岩石变 形破坏阶段特征分析. 试件模型新生裂隙数目与 轴向应变的变化情况如图 6 所示. 岩石裂隙在压密阶段，应变随应力的增加呈 线性增长，未出现新生裂隙. 随着裂隙进一步发 展，岩石压缩进入 AB 阶段、岩石内部开始出现剪 切裂隙，裂隙数目极少并保持稳定；在 BC 阶段，随 荷载增加，剪切裂隙少量增长，开始出现张拉裂 隙. 微裂隙在外荷载作用下压密闭合，剪切裂隙和 张拉裂隙萌生，新生裂隙出现的时间与预制裂隙 倾角无关. 由于应力水平不足以使得原有裂纹扩 展，新萌生的张拉裂隙和剪切裂隙数目较为稳定. 表 2 新生裂隙倾向和倾角分布统计 Table 2 Statistics of the distribution of tendencies and inclinations for newly generated cracks Tendencies and inclinations for preexisting cracks Shear cracks Tension crack Tendency grouping Average inclination/(°) Percentage/% Tendency grouping Average inclination/(°) Percentage/% 90°∠0° 151°–160° 65 6.3 61°–70° 77 8.0 211°–220° 60 6.3 251°–260° 62 8.0 90°∠30° 261°–270° 62 11.9 71°–80° 61 10.5 241°–250° 48 9.5 141°–150° 72 10.5 90°∠45° 261°–270° 54 17.4 71°–80° 52 18.5 251°–260° 45 15.2 91°–100° 64 11.1 90°∠60° 261°–270° 51 9.5 81°–90° 64 14.3 251°–260° 51 7.9 61°–70° 69 10.7 90°∠90° 261°–270° 45 8.5 131°–140° 88 11.1 181°–190° 43 6.4 21°–30° 48 5.6 0 10 20 30 40 50 60 70 80 90 0 10 20 30 40 50 60 70 Shear crack Tension crack Crack initial angle, θ/(°) Crack dip, β/(°) 图 5    预制裂隙倾角 β 与起裂角 θ 的关系 Fig.5    Relation between the crack initial angle θ and crack dip β 张    杰等： 循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗粒流模拟 · 641 ·

642 工程科学学报，第43卷，第5期 Stress-strain curve Stress-strain curve 100 (a) Stress-strain curve under cyclic load 3.0 100 Stress-stran cuve under cyclic load (b)3.0 --Tension crack Shear crack D Tension crack 2.5 ---Shear crack 80 0 2.5 nsion crack u nder cyclic load Sension crack under cyclic load Shear crack under cyclic load 2.0 -Shear crack under cyclic load 2.0 60 60 15 1.5 0 0 10 1.0 20 05 0.5 0 0.1 0.2 0 0.3 0.4 01 0.2 0.3 0.4 Strain/% Strain/% -Stress-strain curve (c) 120 Stress-strain curve (d4.0 3.0 100 Stress-strain curve under cyclic load Tension crack 3.5 -.Shear crack 2.5 100 Shear crack 3.0 80 Sension crack under cvclic load Sension crack under cyclic load Shear crack under cyelic load 80 Shear crack under cyclic load 2.0 60 1.5 60 2.0 40 1.0 20 0.5 0.5 0.2 0.3 0.4 0.2 0.3 0.4 Strain/% Strain/% 140 -Stress-strain curve (e) 120 Stress-strain curve under cyclic load Tension crack Shear crack 100 Sension crack under cyclic o Shear crack under cyclic load 80 60 % 20 0.1 0.2 0.3 0.4 Strain/% 图6新生裂隙数目与轴向应变的变化情况.(a)B=0°：(b)B=30°：(c)B=45°；(d)B=60°：(e)B=90° Fig.6 Number of newly generated cracks and the change in axial strain:(a)B=0°，(b)B=30°：(c)B=45°，(d)B=60°，(e)B=90° 在裂隙稳定扩展阶段(CD阶段)，岩石应力水 30°、45°的试件，张拉裂隙数目增长速率大于剪切 平达到起裂应力，促使内部微裂纹开始扩展.随着 裂隙，张拉裂隙数目超过剪切裂隙时开始出现应 荷载的逐渐增加，张拉裂隙和剪切裂隙共同增长， 力峰值：对于=60°、90°的试件，张拉裂隙和剪切 剪切裂隙数目增长速率大于张拉裂隙，以剪切裂 裂隙数目增长速率基本相同，裂隙扩展以剪切裂 隙扩展为主.在预制裂隙倾角B=45的试件中，新 隙为主.在该阶段，岩石应力-应变曲线的滞回环 生裂隙开始扩展的时间最早，B=30的试件次之 逐渐稀疏，循环后期裂隙倾角=30°、45°、60的岩 在该阶段中，B=45的岩石试件更利于裂隙的扩展. 石试件峰值强度弱化，应力峰值小于加载应力水 在裂隙不稳定扩展阶段(DE阶段)，岩石内部 平上限 裂隙进一步扩展，形成裂隙网络，出现宏观裂纹 裂隙贯通过后，试件完整性并没有完全丧失、 根据图6（新生裂隙增长数目变化图）可以看出，即 岩石仍具有一定的承载力.随着荷载的持续增加， 使在卸荷初期，破裂仍会持续发展，扰动荷载作用 应力快速下降，张拉裂隙和剪切裂隙数目快速增 下微破裂变化比较显著.对于预制裂隙倾角=0°、 长，新生裂隙数目与预制裂隙倾角呈正相关性

在裂隙稳定扩展阶段（CD 阶段），岩石应力水 平达到起裂应力，促使内部微裂纹开始扩展. 随着 荷载的逐渐增加，张拉裂隙和剪切裂隙共同增长， 剪切裂隙数目增长速率大于张拉裂隙，以剪切裂 隙扩展为主. 在预制裂隙倾角 β =45°的试件中，新 生裂隙开始扩展的时间最早，β =30°的试件次之. 在该阶段中，β =45°的岩石试件更利于裂隙的扩展. 在裂隙不稳定扩展阶段（DE 阶段），岩石内部 裂隙进一步扩展，形成裂隙网络，出现宏观裂纹. 根据图 6（新生裂隙增长数目变化图）可以看出，即 使在卸荷初期，破裂仍会持续发展，扰动荷载作用 下微破裂变化比较显著. 对于预制裂隙倾角 β=0°、 30°、45°的试件，张拉裂隙数目增长速率大于剪切 裂隙，张拉裂隙数目超过剪切裂隙时开始出现应 力峰值；对于 β=60°、90°的试件，张拉裂隙和剪切 裂隙数目增长速率基本相同，裂隙扩展以剪切裂 隙为主. 在该阶段，岩石应力−应变曲线的滞回环 逐渐稀疏，循环后期裂隙倾角 β=30°、45°、60°的岩 石试件峰值强度弱化，应力峰值小于加载应力水 平上限. 裂隙贯通过后，试件完整性并没有完全丧失、 岩石仍具有一定的承载力. 随着荷载的持续增加， 应力快速下降，张拉裂隙和剪切裂隙数目快速增 长，新生裂隙数目与预制裂隙倾角呈正相关性. 0 0.1 0.2 0.3 0.4 0 20 40 60 80 100 . C E. . D . B A. 0 0.5 1.0 1.5 2.0 2.5 3.0 Number of cracks/10 3 Number of cracks/10 3 Number of cracks/10 3 Number of cracks/10 3 0 0.5 1.0 1.5 2.0 2.5 3.0 Number of cracks/10 3 . E B . . D . C . A 0 0.1 0.2 0.3 0.4 . E D . . C A. B. 0 0.5 1.0 1.5 2.0 2.5 3.0 0 20 40 60 80 100 120 A. . E . D C. . B 0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 0 0.1 0.2 0.3 0.4 0 20 40 60 80 100 120 140 A. Strain/% E. D. C. . B 0 1 2 3 4 5 (a) (b) (c) (e) (d) Stress/MPa 0 20 40 60 80 100 Stress/MPa Stress/MPa Stress/MPa 0 20 40 60 80 100 Stress/MPa Stress-strain curve Stress-strain curve under cyclic load Tension crack Shear crack Sension crack under cyclic load Shear crack under cyclic load Stress-strain curve Stress-strain curve under cyclic load Tension crack Shear crack Sension crack under cyclic load Shear crack under cyclic load Stress-strain curve Stress-strain curve under cyclic load Tension crack Shear crack Sension crack under cyclic load Shear crack under cyclic load Stress-strain curve Stress-strain curve under cyclic load Tension crack Shear crack Sension crack under cyclic load Shear crack under cyclic load Stress-strain curve Stress-strain curve under cyclic load Tension crack Shear crack Sension crack under cyclic load Shear crack under cyclic load Strain/% 0 0.1 0.2 0.3 0.4 Strain/% Strain/% 0 0.1 0.2 0.3 0.4 Strain/% 图 6    新生裂隙数目与轴向应变的变化情况. （a）β = 0°；（b）β = 30°；（c）β = 45°；（d）β = 60°；（e）β = 90° Fig.6    Number of newly generated cracks and the change in axial strain: (a) β = 0°; (b) β = 30°; (c) β = 45°; (d) β = 60°; (e) β = 90° · 642 · 工程科学学报，第 43 卷，第 5 期

张杰等：循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗粒流模拟 643 3裂隙岩体力学特征 =45的岩石初始应变较大，-0°和30的岩石塑性 变形最小.当循环荷载为峰值强度的40%70% 3.1轴向残余应变特征 和40%~80%时处于稳定阶段（Ⅱ），残余应变随着 在循环荷载作用下，岩体塑性变形随着循环次 循环次数的增加，应变缓慢增长.当循环应力上限 数N的增加而逐渐累积.假设轴向残余应变ε为循 提高到峰值应力0.9倍时，应变值开始迅速增长， 环荷载水平下限两个连续循环之间轴向应变的差值， 进入加速阶段（Ⅲ）：该阶段中，=0的岩石经历 则循环次数与累计残余应变关系如图7所示，图7 3次循环扰动应力首先出现破坏、总循环次数累 中，曲线呈反S形，加载应力初次达到峰值荷载的 计21次，一=30°、45°和60的岩石均在第6次循环 0.4倍为循环次数的起，点.随着循环次数的增加，裂 中残余应变发生突变、产生破环，=90°的岩石经 隙岩体试件在循环应力限值之间表现出不同的特点 过24次循环扰动荷载后残余应变依旧稳定增长、 进入裂隙不稳定扩展阶段(DE阶段)、峰值应力仅 0.28 次于单轴抗压强度.由此可见，扰动荷载应力上限 026 的提高能够加快岩体轴向应变进入加速阶段，促 'ure.ns 0.24 进裂隙的扩展和贯通. 0.22 岩石的轴向应变不仅受应力水平、循环次数 影响，还与岩石预制裂隙的倾角有关.其中，倾角 0.20 =30°和=45的岩石，轴向残余应变最大，受循环 0.18 =0°，-·=30°，+-=45° =60°：◆=90° 荷载影响也最大；大倾角裂隙岩体(B≥60)，轴向 0.16 0 6 12 18 残余应变受循环荷载影响相对较小 Cycles,N 3.2裂隙岩体的强度特征 图7循环次数与轴向应变关系 岩石裂纹从裂隙尖端开始，在扰动应力作用 Fig.7 Relation between the number of cycles and axial strain 下逐渐扩展、破坏模式由剪切破坏为主转变为张 根据模拟结果，当循环荷载为峰值强度的 拉破坏占主导.图8为不同预制裂隙倾角岩石试 40%~60%时岩石处于初始阶段(I),裂隙倾角 件在循环加卸载结束后的破裂模式. Quartz Shear crack Tension crack Pre-existing crack 特 因8不同预制裂隙倾角岩石试件的破裂模式.(a)B=0°：(b)B=30°：(c)B=45°：(d)B=60°：(e)B=90° Fig.8 Fracture modes of a rock specimen with different crack angles::(a)B=0°，(b)B=30°；(c)B=45o;(d)B=60°：(e)B=90 根据图8中裂隙的分布可以看出，裂隙的扩展 90的岩石试件破裂模式与完整试件最为接近、呈 方向与加载方向一致.预制裂隙两端部形成剪切 现劈裂破环.岩石试件模型在单轴压缩和循环扰 破裂带，进一步向上侧和下侧中部发展为破裂面、 动荷载作用下的强度统计结果见表3，表中峰值强 分别与加载两端的局部剪切带连接，导致岩石整 度弱化岩石试件取最大加载应力水平上限 体失稳破坏.在不同预制裂隙倾角（β）中的岩石试 循环扰动应力和裂隙倾角对岩体强度影响很 件中，当B=0时（图8(a)),裂隙中部应力集中区 大：模型试件的峰值强度随裂隙倾角的增大表现 域产生破裂，与端部裂隙共同沿外荷载方向扩展， 出先减小后增大的趋势：扰动应力加剧了新生裂 发生横向破裂面贯通破坏；当B=30°~60时，岩石 隙的扩展，在循环次数一定条件下扰动应力与循 多为横向破裂与竖向劈裂组合贯通破坏；当B= 环荷载的应力上限成正比例增加.预制裂隙倾角

3    裂隙岩体力学特征 3.1    轴向残余应变特征 在循环荷载作用下，岩体塑性变形随着循环次 数 N 的增加而逐渐累积. 假设轴向残余应变 ε 为循 环荷载水平下限两个连续循环之间轴向应变的差值， 则循环次数与累计残余应变关系如图 7 所示，图 7 中，曲线呈反 S 形，加载应力初次达到峰值荷载的 0.4 倍为循环次数的起点. 随着循环次数的增加，裂 隙岩体试件在循环应力限值之间表现出不同的特点. 0 6 12 18 24 0.16 0.18 0.20 0.22 0.24 0.26 0.28 Ⅰ Ⅱ Ⅲ Axial residual strain, ε/% β=0°; β=30°; β=60°; β=90° β=45° Cycles, N 图 7    循环次数与轴向应变关系 Fig.7    Relation between the number of cycles and axial strain 根据模拟结果 ，当循环荷载为峰值强度 的 40%～60% 时岩石处于初始阶段 (Ⅰ)，裂隙倾角 β=45°的岩石初始应变较大，β=0°和 30°的岩石塑性 变形最小. 当循环荷载为峰值强度的 40%～70% 和 40%～80% 时处于稳定阶段 (Ⅱ)，残余应变随着 循环次数的增加，应变缓慢增长. 当循环应力上限 提高到峰值应力 0.9 倍时，应变值开始迅速增长， 进入加速阶段 (Ⅲ)；该阶段中 ， β=0°的岩石经历 3 次循环扰动应力首先出现破坏、总循环次数累 计 21 次，β=30°、45°和 60°的岩石均在第 6 次循环 中残余应变发生突变、产生破坏，β=90°的岩石经 过 24 次循环扰动荷载后残余应变依旧稳定增长、 进入裂隙不稳定扩展阶段（DE 阶段）、峰值应力仅 次于单轴抗压强度. 由此可见，扰动荷载应力上限 的提高能够加快岩体轴向应变进入加速阶段，促 进裂隙的扩展和贯通. 岩石的轴向应变不仅受应力水平、循环次数 影响，还与岩石预制裂隙的倾角有关. 其中，倾角 β=30°和 β=45°的岩石，轴向残余应变最大，受循环 荷载影响也最大；大倾角裂隙岩体（β≥60°），轴向 残余应变受循环荷载影响相对较小. 3.2    裂隙岩体的强度特征 岩石裂纹从裂隙尖端开始，在扰动应力作用 下逐渐扩展、破坏模式由剪切破坏为主转变为张 拉破坏占主导. 图 8 为不同预制裂隙倾角岩石试 件在循环加卸载结束后的破裂模式. Tension crack Mica Feldspar Quartz Shear crack (a) (b) (c) (d) (e) Pre-existing crack 图 8    不同预制裂隙倾角岩石试件的破裂模式. （a）β =0°；（b）β =30°；（c）β =45°；（d）β =60°；（e）β =90° Fig.8    Fracture modes of a rock specimen with different crack angles: (a) β = 0°; (b) β = 30°; (c) β = 45°; (d) β = 60°; (e) β = 90° 根据图 8 中裂隙的分布可以看出，裂隙的扩展 方向与加载方向一致. 预制裂隙两端部形成剪切 破裂带，进一步向上侧和下侧中部发展为破裂面、 分别与加载两端的局部剪切带连接，导致岩石整 体失稳破坏. 在不同预制裂隙倾角（β）中的岩石试 件中，当 β =0°时（图 8（a）），裂隙中部应力集中区 域产生破裂，与端部裂隙共同沿外荷载方向扩展， 发生横向破裂面贯通破坏；当 β =30°～60°时，岩石 多为横向破裂与竖向劈裂组合贯通破坏；当 β = 90°的岩石试件破裂模式与完整试件最为接近、呈 现劈裂破坏. 岩石试件模型在单轴压缩和循环扰 动荷载作用下的强度统计结果见表 3，表中峰值强 度弱化岩石试件取最大加载应力水平上限. 循环扰动应力和裂隙倾角对岩体强度影响很 大：模型试件的峰值强度随裂隙倾角的增大表现 出先减小后增大的趋势；扰动应力加剧了新生裂 隙的扩展，在循环次数一定条件下扰动应力与循 环荷载的应力上限成正比例增加. 预制裂隙倾角 张    杰等： 循环扰动荷载作用下花岗岩中裂隙萌生扩展过程的颗粒流模拟 · 643 ·

644 工程科学学报，第43卷，第5期 表3蜂值强度统计 160 (a)Feldspar 60%,Quartz 20%,Mica 20% Table 3 Statistics of peak strengths 140 (b)Feldspar 40%Quartz40%,Mica20% 120 (c)Feldspar 50%,Quartz 40%,Mica 10% Inclination angle of Peak strength Peak strength under cyclic rock specimen,B/() under cyclic load to the uniaxial strength Cycles load/MPa of the intact rock ratio 三100 0 84.4 0.67 号 30 79.3 0.63 60 (a) b)(c) 45 91.3 0.73 24 9 60 96.7 0.7 24 90 112 0.89 以 0.1 0.2 0.30.40.5 0.6 Strain/% B=30时，岩石破裂过程更为迅速，峰值强度仅为 图9不同矿物比例的岩石应力-应变曲线和破坏模式 实验室完整岩石单轴抗压强度的63%；B=90时， Fig.9 Stress-strain curves and failure modes of rocks with different 石强度最高、为完整岩石单轴抗压强度的89%，循 mineral ratios 环扰动荷载对垂直裂隙影响最小. 矿物组分，按物特性、结合室内单轴抗压强度试 3.3矿物组分比例的影响 验结果进行了细观力学性质参数的赋值，进而模 岩石由一种或多种矿物组合而成，其矿物组 拟了岩石试样循环加卸载试验，重现了不同倾角 分所占比例各不相同.花岗岩矿物组分中，长石质 裂隙的扩展演化过程，得出以下结论： 量分数为40%~60%，石英质量分数为20%~ (1)岩石裂隙端部的集中应力促进了裂隙的 40%.为研究矿物组分比例对裂隙分布及力学性 萌生，不同倾角裂隙对应岩石中的新生裂隙走向 质的影响，设定三种不同矿物配比，长石、石英、 与预制贯通裂隙走向基本一致：裂隙萌生阶段首 云母的矿物颗粒体积比依次为6：2：2、4：4：2 先出现剪切裂隙，张拉裂隙出现较晚 和5：4：1.由于裂隙倾角为30时岩石试件强度 (2)根据新生裂隙的优势倾向分组得到裂隙 受到影响最为显著，故以此试件模型进行单轴压 起裂角与预制裂隙倾角的关系：倾角B≤45时，剪 缩试验 切和张拉裂隙的起裂角单调递减，张拉裂隙比剪 裂隙倾角为30的岩石试件，在不同矿物组分 切裂隙的起裂角大：倾角B≥60时，剪切和张拉裂 比例下的应力-应变曲线和破坏模式如图9所示 隙的起裂角单调递增，剪切裂隙的起裂角较大.由 花岗岩中云母矿物强度最低，随其含量的减少，岩 于裂隙扩展方向与最大拉应力方向保持垂直，当 石试件的峰值应力显著增加；当长石含量最高时， 裂隙扩展角保持不变时，起裂角大的裂隙更倾向 岩石初始压密阶段显著缩短；当云母矿物含量一 于最大主应力方向、裂隙发展最快 定、石英含量逐渐增大时，岩石初始压密阶段特征 (3)循环扰动荷载增加了岩石裂隙的发育程 逐渐突出.预制裂隙岩石试件中，端部裂隙逐渐发 度，这种增加主要体现在裂隙不稳定扩展阶段，张 展和贯通，最终导致失稳破坏.石英含量最高的模 拉裂隙数目增长速率显著提高.相比于静态载荷， 型中，岩石破坏主要是剪切破裂面上裂隙贯通所 循环扰动荷载增加了裂隙岩体的轴向变形；提高 致，破坏模式呈典型的破裂面横向破坏，此时裂隙 循环扰动荷载应力上限，促使岩体轴向残余应变 发展数目最少：随长石含量的增高，裂隙面端部裂 进入加速阶段.岩石轴向残余变形发展情况在不 隙逐渐向着外荷载方向扩展，试件裂隙数目最大、 同裂隙倾角下存在差异，在倾角B=45°和B=90° 更倾向于竖向劈裂破坏 时分别达到最大和最小 岩石中矿物组分比例是影响其宏观力学性质 (4)岩石在裂隙尖端起裂并发展为剪切破裂 的因素之一.因此，在采用数值方法分析岩石材料 带、这一破裂带沿加载方向与加载端部的局部破 变形破坏过程时，应考虑岩石中不同矿物比例的 裂带连接，导致岩石整体失稳破坏.循环扰动应力 影响 下，模型试件的峰值强度为实验室完整岩石单轴 抗压强度的63%~89%，随裂隙倾角的增大表现出 4结论 先减小后增大的趋势.本文所建离散元模型，考虑 建立了反映花岗岩非均质结构特征的颗粒流 了岩石中不同矿物组分的影响，对岩石变形破坏 模型，采用灰度分界阈值分割法识别了岩石中的 机理研究具有参考意义

β =30°时，岩石破裂过程更为迅速，峰值强度仅为 实验室完整岩石单轴抗压强度的 63%；β =90°时 ， 石强度最高、为完整岩石单轴抗压强度的 89%，循 环扰动荷载对垂直裂隙影响最小. 3.3    矿物组分比例的影响 岩石由一种或多种矿物组合而成，其矿物组 分所占比例各不相同. 花岗岩矿物组分中，长石质 量分数 为 40% ～ 60%，石英质量分数 为 20%～ 40%. 为研究矿物组分比例对裂隙分布及力学性 质的影响，设定三种不同矿物配比，长石、石英、 云母的矿物颗粒体积比依次为 6∶2∶2、4∶4∶2 和 5∶4∶1. 由于裂隙倾角为 30°时岩石试件强度 受到影响最为显著，故以此试件模型进行单轴压 缩试验. 裂隙倾角为 30°的岩石试件，在不同矿物组分 比例下的应力−应变曲线和破坏模式如图 9 所示. 花岗岩中云母矿物强度最低，随其含量的减少，岩 石试件的峰值应力显著增加；当长石含量最高时， 岩石初始压密阶段显著缩短；当云母矿物含量一 定、石英含量逐渐增大时，岩石初始压密阶段特征 逐渐突出. 预制裂隙岩石试件中，端部裂隙逐渐发 展和贯通，最终导致失稳破坏. 石英含量最高的模 型中，岩石破坏主要是剪切破裂面上裂隙贯通所 致，破坏模式呈典型的破裂面横向破坏，此时裂隙 发展数目最少；随长石含量的增高，裂隙面端部裂 隙逐渐向着外荷载方向扩展，试件裂隙数目最大、 更倾向于竖向劈裂破坏. 岩石中矿物组分比例是影响其宏观力学性质 的因素之一. 因此，在采用数值方法分析岩石材料 变形破坏过程时，应考虑岩石中不同矿物比例的 影响. 4    结论 建立了反映花岗岩非均质结构特征的颗粒流 模型，采用灰度分界阈值分割法识别了岩石中的 矿物组分，按矿物特性、结合室内单轴抗压强度试 验结果进行了细观力学性质参数的赋值，进而模 拟了岩石试样循环加卸载试验，重现了不同倾角 裂隙的扩展演化过程，得出以下结论： （1）岩石裂隙端部的集中应力促进了裂隙的 萌生，不同倾角裂隙对应岩石中的新生裂隙走向 与预制贯通裂隙走向基本一致；裂隙萌生阶段首 先出现剪切裂隙，张拉裂隙出现较晚. （2）根据新生裂隙的优势倾向分组得到裂隙 起裂角与预制裂隙倾角的关系：倾角 β ≤ 45°时，剪 切和张拉裂隙的起裂角单调递减，张拉裂隙比剪 切裂隙的起裂角大；倾角 β ≥ 60°时，剪切和张拉裂 隙的起裂角单调递增，剪切裂隙的起裂角较大. 由 于裂隙扩展方向与最大拉应力方向保持垂直，当 裂隙扩展角保持不变时，起裂角大的裂隙更倾向 于最大主应力方向、裂隙发展最快. （3）循环扰动荷载增加了岩石裂隙的发育程 度，这种增加主要体现在裂隙不稳定扩展阶段，张 拉裂隙数目增长速率显著提高. 相比于静态载荷， 循环扰动荷载增加了裂隙岩体的轴向变形；提高 循环扰动荷载应力上限，促使岩体轴向残余应变 进入加速阶段. 岩石轴向残余变形发展情况在不 同裂隙倾角下存在差异，在倾角 β = 45°和 β = 90° 时分别达到最大和最小. （4）岩石在裂隙尖端起裂并发展为剪切破裂 带. 这一破裂带沿加载方向与加载端部的局部破 裂带连接，导致岩石整体失稳破坏. 循环扰动应力 下，模型试件的峰值强度为实验室完整岩石单轴 抗压强度的 63%～89%，随裂隙倾角的增大表现出 先减小后增大的趋势. 本文所建离散元模型，考虑 了岩石中不同矿物组分的影响，对岩石变形破坏 机理研究具有参考意义. 表 3    峰值强度统计 Table 3    Statistics of peak strengths Inclination angle of rock specimen, β/(°) Peak strength under cyclic load/MPa Peak strength under cyclic load to the uniaxial strength of the intact rock ratio Cycles 0 84.4 0.67 21 30 79.3 0.63 24 45 91.3 0.73 24 60 96.7 0.77 24 90 112 0.89 24 0 0.1 0.2 0.3 0.4 0.5 0.6 0 20 40 60 80 100 120 140 160 Strain/% Stress/MPa (a) Feldspar 60%, Quartz 20%, Mica 20% (b) Feldspar 40%, Quartz 40%, Mica 20% (c) Feldspar 50%, Quartz 40%, Mica 10% (a) (b) (c) 图 9    不同矿物比例的岩石应力–应变曲线和破坏模式 Fig.9     Stress –strain  curves  and  failure  modes  of  rocks  with  different mineral ratios · 644 · 工程科学学报，第 43 卷，第 5 期