基于矿物晶体模型非均质岩石单轴压缩力学特性研究

刘黎旺; 李海波; 李晓锋; 张国凯; 武仁杰

doi:10.11779/CJGE202003016

基于矿物晶体模型非均质岩石单轴压缩力学特性研究 English Version

刘黎旺^{1, 2,},
李海波^{1, 2},
李晓锋¹,
张国凯³,
武仁杰^{1, 2}

1.
中国科学院武汉岩土力学研究所岩土力学与工程国家重点实验室,湖北　武汉　430071
2.
中国科学院大学,北京　100049
3.
南京理工大学机械工程学院,江苏　南京　210094

基金项目:

国家自然科学基金项目 51439008

国家自然科学基金青年基金项目 51809137

岩土力学与工程国家重点实验室课题 Z017015

详细信息

作者简介:
刘黎旺（1998— ）,男,博士研究生,主要从事岩石动力学及数值计算方法研究。E-mail：lwliu@mail.whrsm.ac.cn

中图分类号: TU452
计量
- 文章访问数: 452
- HTML全文浏览量: 44
- PDF下载量: 230
出版历程
- 收稿日期: 2019-04-03
- 网络出版日期: 2022-12-07
- 刊出日期: 2020-02-29

Research on mechanical properties of heterogeneous rocks using grain-based model under uniaxial compression

LIU Li-wang^{1, 2,},
LI Hai-bo^{1, 2},
LI Xiao-feng¹,
ZHANG Guo-kai³,
WU Ren-jie^{1, 2}

1.
State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China
2.
University of Chinese Academy of Science, Beijing 100049, China
3.
School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China

摘要

摘要: 为研究岩石微观结构非均质性对宏观力学特性及裂纹扩展规律的影响,结合花岗岩室内单轴试验及矿物晶体模型（Grain-based model, GBM）对岩石的微观矿物结构进行离散元建模,研究了非均质岩石宏观应力-应变曲线、声发射及晶体尺度裂纹扩展规律。通过改变岩石微观结构进一步研究了不同矿物组成条件下岩石宏观力学特性及其变化原因。结果表明：GBM可在考虑岩石微结构及强度非均质的条件下,有效模拟岩石加载过程中宏观及微观力学特性；加载初期,岩石内部晶体尺度裂纹首先以晶间裂纹为主,随后转变为以晶内裂纹为主,且裂纹破坏主要为拉伸破坏,岩石破坏时拉伸裂纹及晶内裂纹与总裂纹占比分别约为93.87%,60.95%；加载过程中,岩石微裂纹首先随机分布在岩石内部,后开始聚集并造成宏观破裂面的出现,且宏观破裂面的形成主要与矿物晶体内裂纹的扩展及聚合有关；随岩石内长石矿物含量的增加,岩石峰值应力与损伤应力整体呈增加的变化趋势,且该变化规律可能与长石矿物内裂纹数增加及云母矿物内裂纹数减小有关；通过矿物晶体模型对岩石微观结构进行模拟并对非均质性所造成的力学特性变化进行微观尺度解释,有利于更好理解岩石非均质性与宏观力学特性的关系。
- 非均质岩石 /
- 矿物晶体模型 /
- 晶体尺度裂纹 /
- 声发射 /
- 离散元方法
Abstract: In order to investigate the influences of structural heterogeneity on the mechanical properties and crack growth of rocks, the microstructure of granite is modelled by using grain-based model (GBM). In comparison with the results of laboratory experiments, the properties of stress-strain curves, acoustic emission and grain-scale crack growth of rocks with different mineralogies are investigated. The results show that GBM can be used to efficiently study the macro- and micro-mechanical properties of rocks when the mineral components have different structures and strengths. The grain-scale cracks in rocks are mainly intergranular cracks at the onset of loading, then the intragranular cracks become predominant inversely, and the failure pattern of cracks is always dominated by tensile cracks. When the model fails, the ratios of intergranular and intragranular cracks to total cracks are about 93.87% and 60.95%, respectively. During the whole loading, the locations of microcracks are first located randomly, then the clustering of microcracks leads to the appearance of macroscopic failure surface, and the formation of failure surface is related to the propagation and coalescence of intragranular cracks. When the feldspar mineral in rocks increases, the corresponding peak stress and damage stress increase due to the increase of intragranular cracks formed in feldspar mineral and the decrease of intragranular cracks formed in biotite mineral. In this study, the modelling of micro-structure based on GBM and the reason why heterogeneous rocks behave different mechanical properties may promote our understanding of the influences of rock heterogeneity on the mechanical properties.
- heterogeneous rock /
- grain-based model /
- grain-scale crack /
- acoustic emission /
- discrete element method

HTML全文

0. 引言

岩石材料在地质构造及矿物成形过程中,矿物晶体结构一般具有较明显的非均质性特点,在成分组成、几何尺寸及力学强度^[1]等方面具有较大差异。在外力作用下,岩石矿物可能发生晶体尺度局部应力集中,并导致破坏的发生。岩石破坏后,其宏观裂纹在微观尺度上主要分为两大类：①沿矿物晶体边界发展的晶间裂纹；②矿物晶体内部发展的晶内裂纹（图1）。

图 1 矿物晶体微观结构图^[2-6]

Figure 1. Micro-structures of mineral grains^[2-6]

下载: 全尺寸图片幻灯片

近年来,针对岩石微观结构非均质性对宏观力学特性的影响,国内外学者对非均质性岩石进行了较多的试验研究^[5-7]。Tugrul等^[8]通过试验及数据相关性分析发现,矿物平均晶体尺寸及矿物组成对岩石单轴压缩强度（UCS）有较大的影响；Ündül^[9]通过研究矿物组分与岩石强度之间的关系指出,斜长石矿物含量的增加将导致岩石UCS增加。此外,岩石特征应力与微观矿物结构也具有一定的关系。Ündül等^[10]研究表明岩石起裂应力 $σ_{ci}$ $σ_{ci}$ 受矿物成分影响较大,随斜长石矿物含量增加, $σ_{ci}$ $σ_{ci}$ 呈增加的变化趋势。Cowie等^[2]对花岗岩特征应力进行统计并拟合发现,云母矿物含量与 $σ_{ci}$ $σ_{ci}$ 及 $σ_{cd}$ $σ_{cd}$ 具有负相关的关系,且随岩石中最小矿物尺寸的增加, $σ_{cd}$ $σ_{cd}$ 与 $σ_{ci}$ $σ_{ci}$ 的比值呈增加的变化趋势。然而,Eberhardt等^[11]对比3种不同粒径花岗岩单轴压缩试验结果指出,岩石中矿物晶体的尺寸对 $σ_{ci}$ $σ_{ci}$ 影响不大,但随矿物晶体尺寸增加, $σ_{cd}$ $σ_{cd}$ 呈减小的变化趋势。目前,尽管国内外学者对非均质性岩石宏观力学特性进行了较多试验研究,但天然岩石材料微观矿物结构的随机性使得试验结果表现出较大的离散性,故通过试验手段定量化研究非均质岩石宏观力学仍有一定的困难。

Cundall等^[12]提出的离散单元法（discrete element method, DEM）可在有效研究岩石宏观力学特性的前提下,从微观尺度对岩石内部裂纹的起裂、扩展及聚合规律进行研究。尽管颗粒流软件（particle flow code, PFC）先后开发的两种团簇模型cluster和clump模型^[13-15]可有效克服传统圆盘颗粒所带来的摩擦角及强度值较小问题,但上述模型在研究矿物晶体结构对宏观力学特性影响方面仍具有不足,如cluster模型不能有效控制颗粒簇尺寸,clump模型不能模拟晶体内破裂等；经过多年发展,PFC中矿物晶体模型（PFC-GBM）的提出使得通过岩石内部矿物晶体组成、尺寸等信息进行定量化晶体建模成为了可能,且通过对矿物晶体赋予不同的微观参数即可定量化模拟非均质性岩石宏观力学特性^[16-18]。周喻等^[19]利用上述模型从细观角度揭示了岩石加载条件下的破裂机制及强度特性。Bewick等^[20]研究了非均质岩石直剪条件下宏观力学特性及微裂纹扩展规律。Hofmann等^[21]和Peng等^[22]研究了花岗岩单轴及三轴压缩条件下岩石内部微裂纹的变化规律。综上,采用GBM对非均质性岩石进行矿物晶体尺度模拟是一种十分有效的方法。

本文结合花岗岩室内试验,采用颗粒流矿物晶体模型（PFC-GBM）对花岗岩矿物晶体结构进行建模,通过对比试验与模拟声发射结果并分析晶体尺度微裂纹发展规律研究了单轴压缩条件下非均质岩石微观力学特性变化规律。同时,改变岩石矿物组成进行了一系列单轴压缩试验模拟,分析并讨论了矿物结构非均质性对岩石宏观力学特性的影响及变化原因。

1. PFC-GBM简介及多尺度断裂机理

1.1 PFC-GBM简介

Potyondy^[23]将Voronoi多边形网格引入PFC软件,通过在矿物晶体内赋予平行黏结模型(Parallel bonded model, PBM),矿物晶体间赋予光滑节理模型(Smooth joint model,SJM)可实现非均质岩石内微观矿物晶体建模,其建模流程主要包括以下3个步骤（图2）：①生成大尺寸颗粒以控制矿物晶体尺寸及含量（图2（a））；②连接接触多边形质心形成矿物晶体边界（图2（b）,2（c））；③在晶体边界内填充小尺寸颗粒并进行分组形成矿物晶体模型（图2（d））。

图 2 Potyondy方法建模流程^[23]

Figure 2. Modelling process developed by Potyondy^[23]

下载: 全尺寸图片幻灯片

尽管上述方法的提出使得GBM在PFC软件中实现流程化,但并不能定量化控制岩石的矿物结构,尤其晶体组成、尺寸大小。因此,本文采用其他方法对岩石微观结构进行定量化建模,并将模型数据与PFC软件相结合,成功实现定量化研究岩石非均质性对宏观力学特性的影响。

1.2 多尺度断裂机理

根据断裂机理不同,GBM颗粒间接触破坏所产生的晶体尺度裂纹（晶内/晶间裂纹）可进一步分为拉伸裂纹和剪切裂纹。当矿物晶体内颗粒法向拉应力大于拉伸强度时,颗粒间接触将发生拉伸破坏并产生晶内拉伸裂纹；当矿物晶体内颗粒切向剪切应力大于剪切强度时,颗粒间接触将发生剪切破坏并产生晶内剪切裂纹。同理,矿物晶体间拉伸及剪切裂纹断裂原理与上述相同。

2. 花岗岩模型建立及微观参数确定

2.1 花岗岩模型建立

图3为花岗岩室内试验示意图,试样采用尺寸为50 mm×100 mm的标准圆柱试样。单轴压缩试验在中国科学院武汉岩土力学研究所MTS815.04型电液伺服岩石试验系统上进行,轴向加载速率设置为0.001 mm/s。

图 3 花岗岩室内试验示意图

Figure 3. Schematic of experiments on granite

下载: 全尺寸图片幻灯片

试验过程中,采用轴向引伸计和环向应变规实时测量岩石的轴向应变及环向应变,同时采用声发射装置对岩石破裂过程中声发射特征参数进行监测与采集。

结合室内试验,本文通过Neper软件对花岗岩非均质结构进行建模。建模流程如图4,主要分为3个步骤,首先在50 mm×100 mm的模型区域内生成直径大小服从均匀分布的圆盘颗粒,颗粒半径最小为0.15 mm,且最大与最小半径的比值为1.66（步骤1）；随后在Neper软件中生成矿物晶体多边形,各矿物晶体尺寸使用软件内相关命令设置为服从正态分布（步骤2）；最后将矿物晶体信息以相应的格式导入PFC软件,根据步骤2所生成的矿物晶体种类将步骤1所生成的颗粒分组以形成GBM模型（步骤3）。

图 4 花岗岩矿物晶体模型建模流程图

Figure 4. Grain-based modelling processes of granite

下载: 全尺寸图片幻灯片

2.2 微观参数确定

结合文献[21, 22, 24, 25]中GBM参数校正方法及模拟效果,模型中矿物晶体力学特性（弹性模量、微观拉伸强度及黏结强度等）设置满足以下关系：石英矿物>长石矿物>云母矿物>其他矿物。另外,由于光滑节理模型并不能较好地模拟实际矿物晶体间凹凸不平的接触面及矿物界面间的内锁问题^[26-27],本文采用Li等^[16]的方法对矿物晶体内及晶体间接触赋予平行黏结模型（PBM）,且晶体间微观参数采用表2公式进行折减以简化GBM的参数校正过程。使用“试错法”将微观参数进行校正后,得模拟与试验力学特性比对表见表1。

表 2 模型微观参数表

Table 2. Micro-parameters of grain-based model

参数	石英矿物	长石矿物	云母矿物	其他矿物	矿物晶体间参数
体积含量/%	24	49	24	3
矿物晶体尺寸d/mm	2.43±0.25	3.42±0.25	1.98±0.25	1.17±0.25
颗粒最小半径R^min/mm	0.15
颗粒半径比R^max/R^min	1.66
密度/(kg·m^-3)	2650	2600	3050	1650
线性弹性模量E/GPa	51	42	33	24	0.2×2E¹E²/(E¹+E²)
摩擦系数μ	0.35	0.41	0.33	0.44	0.8
线性刚度比k_n/k_s	1.1	1.3	1.7	3.7	1.3
平行弹性模量pb_emod/GPa	51	42	33	24	0.2×2E¹E²/(E¹+E²)
黏结刚度比k_n/k_s	1.1	1.3	1.7	3.7	1.3
拉伸强度pb_ten/MPa	126±16	105±16	98±13	77±9	21±0
黏结强度pb_coh/MPa	196±42	162±28	146±22	105±0	50±0
摩擦角φ/(°)	19.5	22.4	17.3	23.7	38.6
注：表中E¹,E²分别为晶间两侧矿物颗粒的接触模量值。

下载: 导出CSV

| 显示表格

表 1 宏观力学特性比对表

Table 1. Comparison of mechanical properties

宏观力学特性	试验结果	模拟结果	误差/%
UCSσ_c/MPa	172	171	0.6
拉伸强度σ_t/MPa	13.6	14.3	5.1
弹性模量E/GPa	52.4	55.4	5.7
泊松比 $ν$ $ν$	0.229	0.228	0.44
峰值应变ε_p/%	0.36	0.34	5.6

下载: 导出CSV

| 显示表格

由表可知,使用矿物晶体模型所得宏观力学特性与试验所得结果误差较小,故该模型可在保证有效模拟宏观力学特性的前提下,进一步研究岩石非均质性对宏观力学特性的影响,表2为该模型微观参数表。

3. 矿物晶体模型模拟结果

3.1 应力-应变曲线

图5为花岗岩试验及模拟所得应力-应变曲线。由图5（a）可知,岩石单轴压缩试验力学响应包括以下5个阶段：（Ⅰ）裂纹闭合阶段；（Ⅱ）弹性变形阶段；（Ⅲ）裂纹稳定扩展阶段；（Ⅳ）裂纹不稳定扩展阶段；（Ⅴ）峰后阶段。通过分析裂纹体积应变及试样体积应变的变化关系,可得花岗岩试样特征应力,即闭合应力 $σ_{cc}$ $σ_{cc}$ 、起裂应力 $σ_{ci}$ $σ_{ci}$ 、损伤应力 $σ_{cd}$ $σ_{cd}$ 分别为44.8,71及133.2 MPa。

图 5 试验及模拟应力-应变曲线

Figure 5. Stress-strain curves of experiment and simulation

下载: 全尺寸图片幻灯片

由于没有考虑材料内部的初始缺陷,模拟结果并不能有效地模拟岩石裂纹闭合阶段,故本文模拟所得力学响应存在4个阶段：（Ⅱ）弹性变形阶段；（Ⅲ）裂纹稳定扩展阶段；（Ⅳ）裂纹不稳定扩展阶段；（Ⅴ）峰后阶段。结合文献[28～30]判断起裂应力及损伤应力的方法,当AE事件率开始变化时,此时所对应的轴向应力可认为是起裂应力 $σ_{ci}$ $σ_{ci}$ ；随着轴向加载的进行,AE事件率开始进入稳定变化段,此时岩石内部微裂纹开始稳定起裂并扩展；当AE事件率变化速率快速变化且在后一段加载过程中波动性较大时,可认为此临界应力为损伤应力 $σ_{cd}$ $σ_{cd}$ ,此时岩石内部微裂纹开始不稳定扩展及聚合。本文通过上述方法对特征应力进行判断,可得模拟过程中岩石 $σ_{ci}$ $σ_{ci}$ 及 $σ_{cd}$ $σ_{cd}$ 分别为79.2,143.3 MPa。该结果与试验所得特征应力相差在10 MPa以内,差异性较小。

3.2 声发射演化规律

试验过程中,通过对声发射特征参数进行空间定位处理可反演岩石真实破裂过程。模拟过程中,本文借鉴文献[31, 32]中声发射事件的确定方法将一个微裂纹的产生对应于1次声发射事件,通过对岩石内微裂纹空间位置进行监测并记录,可得到岩石的破裂过程。结果表明,模拟与试验所得声发射事件具有相似的变化规律（图6）。当轴向应力达到岩石起裂应力时,试样进入裂纹稳定扩展阶段,此时微裂纹开始起裂且在空间上呈随机分布状态；当轴向应力达到岩石损伤应力时,试样内微裂纹进一步起裂并发生非稳定扩展,此时微裂纹发生聚集成簇的现象；一旦轴向应力超过岩石峰值应力,微裂纹聚集成簇的程度增加进而导致岩石宏观破裂面的产生。从模拟结果可知,岩石发生破坏时其理想声发射空间定位结果应主要聚集在岩石宏观破坏面附近；相比而言,尽管声发射试验受多方面因素影响,如试验环境、岩石非均质结构等,定位的精确性及规律性并不是很完美,但其结果仍可提高对岩石破裂机理的了解。

图 6 模拟与试验累计声发射规律

Figure 6. Accumulative acoustic emission properties of simulation and experiment

下载: 全尺寸图片幻灯片

3.3 晶体尺度裂纹演化规律

图7为岩石加载过程中晶体尺度裂纹演化规律。当轴向应力超过起裂应力后,岩石内部裂纹起裂主要以矿物晶间裂纹为主,且破坏模式大多为拉伸型破坏,这与Kranz^[33]发现材料拉伸强度普遍小于剪切强度的结论所一致。随轴向加载的进行,晶内拉伸裂纹数量将超过晶间拉伸裂纹并占主导地位。当试样发生破坏时,晶内裂纹与总裂纹数的占比约为60.95%,拉伸裂纹与总裂纹数的占比约为93.87%。

图 7 晶体尺度裂纹变化曲线

Figure 7. Variation curves of growth of grain-scale cracks

下载: 全尺寸图片幻灯片

为验证上述模型模拟所得晶体尺度裂纹发展规律,将试验后花岗岩试样分别进行偏光显微镜及电子显微镜观察。通过在试样表面多处进行观察发现,花岗岩试样受力破坏后所产生的宏观裂纹在矿物晶体尺度主要以不同矿物内部裂纹为主,少部分为矿物晶间裂纹,且该试验结果与文献[34]所得结论一致,故GBM模型可较真实地模拟岩石外力作用下晶体尺度裂纹演化规律。花岗岩代表性微裂纹分析结果如图8所示,图中Qz为石英矿物,Bt为云母矿物,Pl为斜长石,Kf为钾长石。

图 8 花岗岩微裂纹分析结果

Figure 8. Micro-cracks of granite

下载: 全尺寸图片幻灯片

图9为岩石加载过程中矿物晶体内裂纹演化规律。加载初期,长石矿物晶体内裂纹首先发生起裂,随后石英、云母等矿物内裂纹依次开始起裂并扩展；当试样发生破坏时,矿物内裂纹数量关系如下：长石内裂纹>云母内裂纹>石英内裂纹>其他矿物内裂纹。由于长石矿物在此模型中的含量较高且力学特性相对石英矿物较差,故其内部裂纹最先发生起裂且数量始终占主导地位；相比云母矿物与石英矿物,尽管两者含量相近,但云母矿物相比石英矿物具有较差的力学性质,故岩石发生破坏时,云母内裂纹数量多于石英内裂纹。

图 9 矿物晶体内裂纹演化规律

Figure 9. Variation of intragranular cracks

下载: 全尺寸图片幻灯片

在此,为进一步研究矿物晶体裂纹空间演化特征,图9中各加载点（A,B,C,D点）处晶体尺度裂纹分布图如图10所示。由图可知,当轴向应力达到A点时,模型中晶体尺度裂纹主要为矿物晶间裂纹,且裂纹随机分布在模型中；当轴向应力达到B点时,晶体尺度裂纹开始发生聚集现象,此时宏观破裂面开始形成。随着微裂纹的聚集程度加深,岩石宏观破裂面的发展及贯通将导致模型发生破坏。另外,观察宏观破裂面处晶体内裂纹分布发现,长石及云母矿物内裂纹占主导地位,石英矿物内裂纹次之,故在此可推断岩石发生破坏时微裂纹扩展及聚合主要以矿物晶体内裂纹为主,尤其长石及云母矿物内裂纹。

图 10 晶体尺度裂纹分布演化图

Figure 10. Evolution of crack distributions in grain-scale

下载: 全尺寸图片幻灯片

3.4 不同长石矿物含量岩石力学特性

在上述模型基础上,为研究不同矿物组成岩石宏观力学特性,本文在控制石英及其他矿物成分含量不变的条件下,通过Neper软件改变岩石中长石含量w以10%的梯度从20%变化至60%,生成5组不同的GBM模型,且每组模型进行3次共计15次数值模拟。模型如图11所示,黑色多边形为云母矿物,灰色多边形为长石矿物,白色多边形为石英矿物,蓝色多边形为其他矿物。

图 11 不同长石含量矿物晶体模型

Figure 11. Grain-based models with different feldspar contents

下载: 全尺寸图片幻灯片

图12为不同长石含量岩石峰值应力及特征应力变化规律。由图12（a）可知,模拟所得峰值应力整体随长石矿物含量增加呈增加的变化趋势；当长石含量从20%增加至60%时,模拟所得峰值应力从168 MPa变化至178 MPa,但其变化幅度与文献[2]中统计拟合结果相比较小,这可能与模拟及文献[2]中岩石微观结构差异较大具有一定的联系。另外,随长石矿物含量增加,岩石损伤应力呈轻微增加的变化趋势,而起裂应力受长石矿物含量影响较小（图12（b））。

图 12 不同长石含量模型峰值应力及特征应力变化曲线

Figure 12. Curves of peak and characteristic stresses of models with different feldspar contents

下载: 全尺寸图片幻灯片

Ündül^[9]通过试验指出,斜长石矿物与矿物基质含量占比的增加将增大单轴压缩强度,而矿物含量与起裂应力之间并没有较明确的关系,故此模拟所得结果与试验结果相吻合。

现从微观角度对不同长石含量岩石单轴压缩条件下晶体尺度裂纹变化规律进行分析,以进一步研究并讨论峰值应力变化的原因。将上述模型破坏时矿物晶体内裂纹数量取平均值并与总裂纹数进行正则化处理后,得不同长石含量GBM模型矿物晶体内裂纹占比变化规律如图13。

图 13 不同长石含量模型矿物内裂纹变化规律

Figure 13. Variations of intragranular cracks of models with different feldspar contents

下载: 全尺寸图片幻灯片

由图13可知,随长石矿物含量增加,岩石云母矿物内裂纹占比呈递减的变化趋势,而长石矿物内裂纹占比呈递增的变化趋势,石英及其他矿物内裂纹占比则变化不大,故长石矿物含量的升高将导致岩石内部裂纹形式从以云母内裂纹为主转变为以长石内裂纹为主,模型破坏时相应的晶体尺度裂纹分布见图14。随长石矿物含量增加,岩石宏观破裂面处云母内裂纹数量减少,而长石内裂纹数量增加；相比云母矿物而言,长石矿物由于具有更好的力学性质且内部裂纹发生扩展所需要的能量越多,故随着长石含量的增加,岩石内部裂纹形式从以云母内裂纹为主转变为以长石内裂纹为主可能导致岩石峰值应力整体呈增加变化趋势。

图 14 不同长石含量模型晶体尺度裂纹分布图

Figure 14. Distributions of grain-scale cracks of models with different feldspar contents非均质岩石力学特性的数值方法。

下载: 全尺寸图片幻灯片

4. 结论

本文结合花岗岩室内单轴试验并采用矿物晶体模型（GBM）对岩石矿物非均质结构进行了离散元模拟,改变岩石矿物组成研究了岩石微观非均质性对宏观力学特性的影响,并从矿物晶体尺度对其变化原因进行了较为合理的解释。主要结论如下：

（1）矿物晶体模型（GBM）可在考虑岩石矿物晶体结构及强度非均质性的前提下,对其应力-应变曲线、声发射特性进行有效模拟,并可研究矿物晶体尺度微裂纹起裂、扩展及聚合规律,是一种有效研究

（2）岩石单轴压缩过程中,微裂纹起裂首先以晶间裂纹为主,随后以晶内裂纹为主,且主要破坏机理为拉伸型破坏。岩石破坏时,拉伸裂纹及晶内裂纹与总裂纹的占比分别约为93.87%,60.95%。

（3）岩石单轴压缩过程中,长石矿物内裂纹首先起裂,随后云母、石英等矿物内裂纹依次发生起裂并扩展。岩石破坏时,微裂纹聚集将导致宏观破裂面的形成,且宏观破裂面的形成主要与晶体内裂纹扩展及聚合有关。

（4）随长石矿物含量增加,岩石峰值应力及损伤应力整体呈增加的变化趋势,而起裂应力变化并不明显。从矿物晶体尺度裂纹变化规律可推断,岩石峰值应力增大的原因可能与长石内裂纹占比增加、云母内裂纹占比减少有关。

图 1 矿物晶体微观结构图^[2-6]

Figure 1. Micro-structures of mineral grains^[2-6]

下载: 全尺寸图片幻灯片

图 2 Potyondy方法建模流程^[23]

Figure 2. Modelling process developed by Potyondy^[23]

下载: 全尺寸图片幻灯片

图 3 花岗岩室内试验示意图

Figure 3. Schematic of experiments on granite

下载: 全尺寸图片幻灯片

图 4 花岗岩矿物晶体模型建模流程图

Figure 4. Grain-based modelling processes of granite

下载: 全尺寸图片幻灯片

图 5 试验及模拟应力-应变曲线

Figure 5. Stress-strain curves of experiment and simulation

下载: 全尺寸图片幻灯片

图 6 模拟与试验累计声发射规律

Figure 6. Accumulative acoustic emission properties of simulation and experiment

下载: 全尺寸图片幻灯片

图 7 晶体尺度裂纹变化曲线

Figure 7. Variation curves of growth of grain-scale cracks

下载: 全尺寸图片幻灯片

图 8 花岗岩微裂纹分析结果

Figure 8. Micro-cracks of granite

下载: 全尺寸图片幻灯片

图 9 矿物晶体内裂纹演化规律

Figure 9. Variation of intragranular cracks

下载: 全尺寸图片幻灯片

图 10 晶体尺度裂纹分布演化图

Figure 10. Evolution of crack distributions in grain-scale

下载: 全尺寸图片幻灯片

图 11 不同长石含量矿物晶体模型

Figure 11. Grain-based models with different feldspar contents

下载: 全尺寸图片幻灯片

图 12 不同长石含量模型峰值应力及特征应力变化曲线

Figure 12. Curves of peak and characteristic stresses of models with different feldspar contents

下载: 全尺寸图片幻灯片

图 13 不同长石含量模型矿物内裂纹变化规律

Figure 13. Variations of intragranular cracks of models with different feldspar contents

下载: 全尺寸图片幻灯片

图 14 不同长石含量模型晶体尺度裂纹分布图

Figure 14. Distributions of grain-scale cracks of models with different feldspar contents非均质岩石力学特性的数值方法。

下载: 全尺寸图片幻灯片

表 2 模型微观参数表

Table 2 Micro-parameters of grain-based model

参数	石英矿物	长石矿物	云母矿物	其他矿物	矿物晶体间参数
体积含量/%	24	49	24	3
矿物晶体尺寸d/mm	2.43±0.25	3.42±0.25	1.98±0.25	1.17±0.25
颗粒最小半径R^min/mm	0.15
颗粒半径比R^max/R^min	1.66
密度/(kg·m^-3)	2650	2600	3050	1650
线性弹性模量E/GPa	51	42	33	24	0.2×2E¹E²/(E¹+E²)
摩擦系数μ	0.35	0.41	0.33	0.44	0.8
线性刚度比k_n/k_s	1.1	1.3	1.7	3.7	1.3
平行弹性模量pb_emod/GPa	51	42	33	24	0.2×2E¹E²/(E¹+E²)
黏结刚度比k_n/k_s	1.1	1.3	1.7	3.7	1.3
拉伸强度pb_ten/MPa	126±16	105±16	98±13	77±9	21±0
黏结强度pb_coh/MPa	196±42	162±28	146±22	105±0	50±0
摩擦角φ/(°)	19.5	22.4	17.3	23.7	38.6
注：表中E¹,E²分别为晶间两侧矿物颗粒的接触模量值。

下载: 导出CSV

表 1 宏观力学特性比对表

Table 1 Comparison of mechanical properties

宏观力学特性	试验结果	模拟结果	误差/%
UCSσ_c/MPa	172	171	0.6
拉伸强度σ_t/MPa	13.6	14.3	5.1
弹性模量E/GPa	52.4	55.4	5.7
泊松比 $ν$ $ν$	0.229	0.228	0.44
峰值应变ε_p/%	0.36	0.34	5.6

下载: 导出CSV

参考文献(34)

[1]	BASS J D. Elasticity of minerals, glasses, and melts[J]. Mineral Physics & Crystallography: a Handbook of Physical Constants, 1995, 2: 45-63.
[2]	COWIE S, WALTON G. The effect of mineralogical parameters on the mechanical properties of granitic rocks[J]. Engineering Geology, 2018, 240: 204-225. doi: 10.1016/j.enggeo.2018.04.021
[3]	YILMAZ N G, GOKTAN R M, KIBICI Y. Relations between some quantitative petrographic characteristics and mechanical strength properties of granitic building stones[J]. International Journal of Rock Mechanics and Mining Sciences, 2011, 3(48): 506-513.
[4]	LI X F, LI X, LI H B, et al. Dynamic tensile behaviours of heterogeneous rocks: the grain scale fracturing characteristics on strength and fragmentation[J]. International Journal of Impact Engineering, 2018, 118: 98-118. doi: 10.1016/j.ijimpeng.2018.04.006
[5]	SAJID M, COGGAN J, ARIF M, et al. Petrographic features as an effective indicator for the variation in strength of granites[J]. Engineering Geology, 2016, 202: 44-54. doi: 10.1016/j.enggeo.2016.01.001
[6]	TANDON R S, GUPTA V. The control of mineral constituents and textural characteristics on the petrophysical & mechanical (PM) properties of different rocks of the Himalaya[J]. Engineering Geology, 2013, 153: 125-143. doi: 10.1016/j.enggeo.2012.11.005
[7]	PŘIKRYL R. Some microstructural aspects of strength variation in rocks[J]. International Journal of Rock Mechanics and Mining Sciences, 2001, 38(5): 671-682. doi: 10.1016/S1365-1609(01)00031-4
[8]	TUĞRUL A, ZARIF I H. Correlation of mineralogical and textural characteristics with engineering properties of selected granitic rocks from Turkey[J]. Engineering Geology, 1999, 51(4): 303-317. doi: 10.1016/S0013-7952(98)00071-4
[9]	ÜNdül Ö. Assessment of mineralogical and petrographic factors affecting petro-physical properties, strength and cracking processes of volcanic rocks[J]. Engineering Geology, 2016, 210: 10-22. doi: 10.1016/j.enggeo.2016.06.001
[10]	ÜNDÜL Ö, AMANN F, AYSAL N, et al. Micro-textural effects on crack initiation and crack propagation of andesitic rocks[J]. Engineering Geology, 2015, 193: 267-275. doi: 10.1016/j.enggeo.2015.04.024
[11]	EBERHARDT E, STIMPSON B, STEAD D. Effects of grain size on the initiation and propagation thresholds of stress-induced brittle fractures[J]. Rock Mechanics and Rock Engineering, 1999, 32(2): 81-99. doi: 10.1007/s006030050026
[12]	CUNDALL P A, STRACK O D L. A discrete numerical model for granular assemblies[J]. Géotechnique, 1979, 29(1): 47-65. doi: 10.1680/geot.1979.29.1.47
[13]	LI X F, LI H B, LIU Y Q, et al. Numerical simulation of rock fragmentation mechanisms subject to wedge penetration for TBMs[J]. Tunnelling and Underground Space Technology, 2016, 53: 96-108. doi: 10.1016/j.tust.2015.12.010
[14]	CHO N, MARTIN C D, SEGO D C. A clumped particle model for rock[J]. International Journal of Rock Mechanics and Mining Sciences, 2007, 44(7): 997-1010. doi: 10.1016/j.ijrmms.2007.02.002
[15]	刘广, 荣冠, 彭俊, 等. 矿物颗粒形状的岩石力学特性效应分析[J]. 岩土工程学报, 2013, 35(3): 540-550. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201303020.htm LIU Guang, RONG Guan, PENG Jun, et al. Mechanical behaviors of rock affected by mineral particle shapes[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(3): 540-550. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201303020.htm
[16]	LI X F, ZHANG Q B, LI H B, et al. Grain-based discrete element method (GB-DEM) modelling of multi-scale fracturing in rocks under dynamic loading[J]. Rock Mechanics and Rock Engineering, 2018, 51(12): 3785-3817. doi: 10.1007/s00603-018-1566-2
[17]	LAN H, MARTIN C D, HU B. Effect of heterogeneity of brittle rock on micromechanical extensile behavior during compression loading[J]. Journal of Geophysical Research: Solid Earth, 2010, 115(B01202): 1-14.
[18]	LI X F, LI H B, ZHAO J. 3D polycrystalline discrete element method (3PDEM) for simulation of crack initiation and propagation in granular rock[J]. Computers and Geotechnics, 2017, 90: 96-112. doi: 10.1016/j.compgeo.2017.05.023
[19]	周喻, 高永涛, 吴顺川, 等. 效晶质模型及岩石力学特征细观研究[J]. 岩石力学与工程学报, 2015, 34(3): 511-519. ZHOU Yu, GAO Yong-tao, WU Shun-chuan, et al. An equivalent crystal model for mesoscopic behaviour of rock[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(3): 511-519. (in Chinese)
[20]	BEWICK R P, KAISER P K, BAWDEN W F, et al. DEM simulation of direct shear: 1 Rupture under constant normal stress boundary conditions[J]. Rock Mechanics and Rock Engineering, 2014, 47(5): 1647-1671. doi: 10.1007/s00603-013-0490-8
[21]	HOFMANN H, BABADAGLI T, ZIMMERMANN G. A grain based modeling study of fracture branching during compression tests in granites[J]. International Journal of Rock Mechanics and Mining Sciences, 2015, 77: 152-162. doi: 10.1016/j.ijrmms.2015.04.008
[22]	PENG J, WONG L N Y, TEH C I, et al. Modeling micro-cracking behavior of Bukit Timah granite using grain-based model[J]. Rock Mechanics and Rock Engineering, 2018, 51(1): 135-154. doi: 10.1007/s00603-017-1316-x
[23]	POTYONDY D O. A grain-based model for rock: approaching the true microstructure[C]//Proceedings of the Rock Mechanics in the Nordic Countries, 2010, Kongsberg: 225-234.
[24]	ZHANG Y, WONG L N Y, CHAN K K. An Extended Grain‐Based Model Accounting for Microstructures in Rock Deformation[J]. Journal of Geophysical Research: Solid Earth, 2019, 124(1): 125-148. doi: 10.1029/2018JB016165
[25]	PENG J, WONG L N Y, TEH C I. Influence of grain size heterogeneity on strength and microcracking behavior of crystalline rocks[J]. Journal of Geophysical Research: Solid Earth, 2017, 122(2): 1054-1073. doi: 10.1002/2016JB013469
[26]	LI X F, LI H B, ZHAO J. The role of transgranular capability in grain-based modelling of crystalline rocks[J]. Computers and Geotechnics, 2019, 110: 161-183. doi: 10.1016/j.compgeo.2019.02.018
[27]	MEHRANPOUR M H, KULATILAKE P H S W. Improvements for the smooth joint contact model of the particle flow code and its applications[J]. Computers and Geotechnics, 2017, 87: 163-177. doi: 10.1016/j.compgeo.2017.02.012
[28]	DIEDERICHS M S, KAISER P K, EBERHARDT E. Damage initiation and propagation in hard rock during tunnelling and the influence of near-face stress rotation[J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(5): 785-812. doi: 10.1016/j.ijrmms.2004.02.003
[29]	张国凯, 李海波, 王明洋, 等. 岩石单轴压缩下损伤表征及演化规律对比研究[J]. 岩土工程学报, 2019, 41(6): 1074-1082. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201906013.htm ZHANG Guo-kai, LI Hai-bo, WANG Ming-yang, et al. Comparative study on damage characterization and damage evolution of rock under uniaxial compression[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(6): 1074-1082. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201906013.htm
[30]	朱俊, 邓建辉, 黄弈茗, 等. 饱和大理岩特征强度试验研究[J]. 岩石力学与工程学报, 2019, 38(6): 1129-1138. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201906005.htm ZHU Jun, DENG Jian-hui, HUANG Yi-min, et al. Experimental study on the characteristic strength of saturated marble[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(6): 1129-1138. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201906005.htm
[31]	HAZZARD J F, YOUNG R P, MAXWELL S C. Micromechanical modeling of cracking and failure in brittle rocks[J]. Journal of Geophysical Research: Solid Earth, 2000, 105(B7): 16683-16697. doi: 10.1029/2000JB900085
[32]	穆康, 李天斌, 俞缙, 等. 围压效应下砂岩声发射与压缩变形关系的细观模拟[J]. 岩石力学与工程学报, 2014, 33(增刊1): 2786-2793. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2014S1028.htm MU Kang, LI Tian-bin, YU Jin, et al. Mesoscopic simulation of relationship of acoustic emission and compressive deformation behavior in sandstone under confining pressure effect[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(S1): 2786-2793. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2014S1028.htm
[33]	KRANZ R L. Crack-crack and crack-pore interactions in stressed granite[C]//International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1979, Pergamon: 37-47.
[34]	MOORE D E, LOCKNER D A. The role of microcracking in shear-fracture propagation in granite[J]. Journal of Structural Geology, 1995, 17(1): 95-114. doi: 10.1016/0191-8141(94)E0018-T

施引文献(42)

期刊类型引用(20)

1.	孟珂，王笑梅，杜晓冉，张晓曼，罗娟. 罕见多矿物晶体共生标本的综合鉴定. 矿产综合利用. 2025(01): 200-205 . 百度学术
2.	李红丽. 非均质岩石单轴压缩下损伤演化规律数值模拟研究. 有色矿冶. 2024(01): 43-48 . 百度学术
3.	鞠明和，陶泽军，蔚立元，姜礼杰，郑彦龙，邹春江. 钢粒子迟滞重复冲击破岩硬岩损伤破裂特征研究. 岩土力学. 2024(04): 1242-1255 . 百度学术
4.	吴泽兵，袁若飞，张文溪，王刚，胡诗尧. PDC混合布齿钻头破碎非均质花岗岩数值模拟. 天然气工业. 2024(05): 105-117 . 百度学术
5.	周元，吕威帆，王颖轶. 基于块体离散元法的盾构掘进围岩与管片变形模拟研究. 都市快轨交通. 2024(03): 125-134 . 百度学术
6.	裴书锋，郝文锋，王营利，王一汀，曾凤娟. 双江口水电站花岗岩单轴压缩微观破坏机制研究. 西北水电. 2024(04): 62-68 . 百度学术
7.	张国桥，孙鹏，吴祥业，王婧雅，郭文斌，田宇航. 基于PFC-GBM非均质模型的砂岩裂纹演化细观规律研究. 中国矿业. 2024(09): 158-169 . 百度学术
8.	冯龙飞，王双明，王晓东，解嘉豪，窦林名. 煤单轴峰后动态冲击破坏特征及差异机制模拟研究. 煤炭学报. 2024(S2): 714-730 . 百度学术
9.	马文强，王酒婷. 花岗岩受压宏-细观破坏特征及能量演化规律. 信阳师范学院学报(自然科学版). 2023(02): 314-320 . 百度学术
10.	张涛，蔚立元，苏海健，高亚楠，贺虎，魏江波. 基于多级力链网络分析的花岗岩压缩特性的矿物尺寸效应研究. 岩石力学与工程学报. 2023(08): 1988-2003 . 百度学术
11.	乔世范，刘钰，王刚，张细宝，张海凤，董常瑞，谭晶仁，檀俊坤. 考虑岩石细观结构的TBM滚刀破岩过程数值研究. 中国安全生产科学技术. 2023(07): 106-112 . 百度学术
12.	向衍斌. 煤系岩石单轴压缩损伤破坏演化规律与表征. 煤矿安全. 2023(09): 88-95 . 百度学术
13.	赵光明，高宇，吴旭坤. 岩石变刚度实验条件下力学与声发射特性. 安徽理工大学学报(自然科学版). 2023(06): 63-72 . 百度学术
14.	王桂林，王润秋，孙帆. 块体离散元颗粒模型细观参数标定方法及花岗岩细观演化模拟. 长江科学院院报. 2022(01): 86-93 . 百度学术
15.	张涛，蔚立元，鞠明和，李明，苏海健，季浩奇. 基于PFC3D-GBM的晶体–单元体尺寸比对花岗岩动态拉伸特性影响分析. 岩石力学与工程学报. 2022(03): 468-478 . 百度学术
16.	兰恒星，包含，孙巍锋，刘世杰. 岩体多尺度异质性及其力学行为. 工程地质学报. 2022(01): 37-52 . 百度学术
17.	李博，梁秦源，周宇，赵程，伍法权. 基于CT-GBM重构法的花岗岩裂纹扩展规律研究. 岩石力学与工程学报. 2022(06): 1114-1125 . 百度学术
18.	Tongzhao Zhang，Hongguang Ji，Xiaobo Su，Shuang You，Daolu Quan，Zhou Zhang，Jinzhe Li. Evaluation and classification of rock heterogeneity based on acoustic emission detection. International Journal of Minerals, Metallurgy and Materials. 2022(12): 2117-2125 . 必应学术
19.	郑强强，徐颖，胡浩，钱佳威，宗琦，谢平. 单轴荷载作用下砂岩的破裂与速度结构层析成像. 岩土工程学报. 2021(06): 1069-1077 . 本站查看
20.	李博，朱强，张丰收，赵程，伍法权. 基于矿物晶体模型的非均质性岩石双裂纹扩展规律研究. 岩石力学与工程学报. 2021(06): 1119-1131 . 百度学术