岩石细观颗粒结构具有明显离散性, 颗粒尺寸会影响岩石损伤演化过程和变形路径, 决定了岩石宏观力学性质和破裂模式[1].抗拉强度和断裂韧度作为重要力学参数对于岩体工程稳定性意义重大[2-4], 考虑测试技术复杂性, 断裂韧度与抗拉强度一般均采用间接法测试, 例如巴西劈裂与三点弯曲试验, 已有部分经验公式用于推导抗拉强度、断裂韧度[5].
已有研究从岩石细观尺度角度, 探究内部矿物颗粒尺寸、含量及分布对宏观力学性质影响.周喻等[6]采用颗粒体模型与光滑节理模型, 构建具有岩石矿物细观结构特征的等效晶质模型, 揭示岩石在加载过程中破裂机制与强度特性.Wong等[7]研究了内部矿物颗粒大小、排布方式等对于岩石起裂压力、裂纹扩展方向的影响.Wang等[8]采用虚裂纹模型对花岗岩试样进行三点弯曲试验, 间接测试出试样的抗拉强度和断裂韧度, 方法新颖易行.
本文引入虚裂纹模型量化试样缺口尖端的断裂过程区域, 进行三点弯曲试验计算岩石试样断裂韧度和抗拉强度, 并结合岩石细观数值计算软件RFPA-DIP版, 建立数值模型, 进行数值模拟, 并通过分形维数对三种颗粒尺寸岩石细观结构进行表征.
1 考虑虚裂纹模型裂纹尖端应力分布 1.1 三点弯曲试样缺口尖端应力分布岩石裂缝扩展时, 裂尖出现的微裂纹区被称为断裂过程区(fracture process zone).引入虚拟裂纹模型表征该区域的应力分布, 图 1所示为岩石三点弯曲试验过程中预制缺口尖端应力场的分布模型, 用式(1) 描述尖端应力分布(单位mm):图中S代表试样下部两支撑点的跨度, W代表试件宽度, Pmax为峰值载荷, ao为试样预制缺口的长度, Δafic取值试样平均颗粒尺寸, σn与σc分别为缺口区域的名义强度与压应力, x和y分别代表试验达到峰值荷载点σn和σc的分布范围.
(1) |
名义强度与抗拉强度的关系已在文献[8]中给出了解析解, 等效裂纹长度ae, 试样几何尺寸B(α)可由式(3), (4) 得到.
(2) |
(3) |
(4) |
(5) |
(6) |
(7) |
式中:α=a0/W;ae为等效裂纹长度;a∞*为特征裂纹长度;ft为抗拉强度;KIC为断裂韧度.特征裂纹长度a∞*是一个常数, 由抗拉强度和断裂韧度决定.利用式(7), 在准静态加载情况下, 通过大量试验结果拟合出抗拉强度、断裂韧度.
2 三点弯曲数值模型建立 2.1 花岗岩非均匀性表征利用数字图像处理技术表征颗粒非均匀性, 同时对岩石内部的细观介质分类.图 2为花岗岩数字图像, 像素尺寸为280×70, 模型尺寸为280 mm×70 mm.利用图像处理的多阈值分割技术区分不同矿物种类.图 2中, 红色区域代表长石, 黑色代表云母, 白色代表石英.图 3所示为数值图像某截面上矿物介质亮度值的变化, 经过反复试验, 将分割阈值设为50和200, 亮度阈值被分为三个区间段(0~50, 50~200, 200~255), 因此试样内部矿物颗粒也被分为三类.
图 4为阈值分割得到的数值模型, 由矩形排列的像素点组成, 每个像素点灰度不同, 赋予不同材料参数, 进而在数值模型中计入非均匀性影响.
选取颗粒尺寸为1 mm的花岗岩试样, 并预制不同长度缺口进行数值模拟, 通过与室内试验结果对比验证所建立模型有效性.位移初始值为0.001 mm, 单步增量为0.001 mm, 加载直至试样破坏.图 5为数值模拟与物理试验结果[8]比较, 可以发现模拟结果与试验结果较为接近.
选取如图 6所示三种不同颗粒尺寸花岗岩进行试验, 对其进行数字图像处理, 统计试样内部细观颗粒分布规律如图 7所示.得到平均颗粒尺寸分别为1.02, 2.12, 3 mm.
矿物颗粒力学参数见表 1, 单一颗粒为均匀介质, 不计各自内部非均匀性.
根据三种不同颗粒尺寸建立花岗岩数值模型,长、宽分别为280, 70 mm, 预制长度分别为0,1, 5, 10, …, 40 mm缺口(每5 mm预制), 如图 8所示.
试样裂纹扩展模式如图 9所示.第300步, 缺口尖端以拉应力导致积聚损伤(红色圆圈表示), 符合Hillerborg[10]提出的主虚裂纹模型特征, 此阶段对应图 10中弹性段.第348步, 如图 9b所示损伤积聚形成细微裂纹在缺口尖端开始萌生并逐渐扩展, 图 10中可见此区域曲线处于塑性阶段, 并进入虚裂纹产生区域.第380步达到峰值载荷, 随着加载的继续, 裂纹发生失稳扩展, 试样完全破坏(图 9d所示).如图 10所示微裂纹从萌生到突然发生断裂实际发生在第348~380步很短的加载区间内, 可以看出花岗岩主要表现出脆性破坏.从裂纹扩展路径观察到, 裂纹数量少、分支不明显, 裂缝面无剧烈折曲.
将峰值载荷Pmax代入式(1), 求得名义强度σn, 结果如表 2所示, 等效裂纹长度ae可通过式(3) 求得.σn, ae代入式(7), 采用线性回归方法计算抗拉强度和断裂韧度, 结果如图 11所示.
根据文献[11], 颗粒粒度分布具有分形特征:
(8) |
其中: r为包含所有特征颗粒尺寸;m(r)为尺寸小于r矿物颗粒累计质量;M为试样内部颗粒总质量;rm代表分布参数, 当r=rm时, r取值为最大颗粒尺寸;b代表颗粒分布参数.分形维数可用D=3-b计算得出, D代表不同颗粒尺寸的分维值.
图 12中可看出lnr与ln[m(r)/M]线性相关, 这表明岩石内部颗粒符合分形统计分布, 具有自相似性.得到三类颗粒试样b值分别为0.982, 0.855, 0.796, 因此其分维值分别为2.018, 2.145, 2.204.
图 13表示三类花岗岩颗粒分形维数与其抗拉强度和断裂韧度关系, 得出:抗拉强度和断裂韧度与分维值正相关.因此分形维数可以用来定量描述岩石试样内部细观颗粒分布的非均匀性.
1) 本文建立虚裂纹与岩石颗粒尺寸关系的模型, 利用细观岩石力学数值模拟软件RFPA-Dip版, 计算出三种颗粒尺寸花岗岩试样抗拉强度和断裂韧度.
2) 建立不同颗粒尺寸花岗岩试样数值模型, 验证模型的合理性, 加载过程中可观察到其裂纹数量少、分支不明显, 且裂缝面无剧烈折曲, 此现象体现出花岗岩试样的脆性特性.
3) 计算三类颗粒尺寸岩石试样分维值, 由结果可知:颗粒尺寸越大, 分维值越小.利用分形维数建立试样细观非均匀性与宏观强度特性的关系, 并得出结论:随着分形维数增大, 花岗岩试样的抗拉强度和断裂韧度均减小.
[1] |
于庆磊, 杨天鸿, 郑超.
岩石细观结构对其变形强度影响的数值分析[J]. 岩土力学, 2011, 32(11): 3468–3473.
( Yu Qing-lei, Yang Tian-hong, Zheng Chao. Numerical analysis of influence of rock mesostructure on its deformation and strength[J]. Rock and Solid Mechanics, 2011, 32(11): 3468–3473. DOI:10.3969/j.issn.1000-7598.2011.11.044 ) |
[2] | Dai F, Wei M D, Xu N W, et al. Numerical investigation of the progressive fracture mechanisms of four ISRM-suggested specimens for determining the mode I fracture toughness of rocks[J]. Computers and Geotechnics, 2015, 69: 424–441. DOI:10.1016/j.compgeo.2015.06.011 |
[3] | Zhang Z X. An empirical relation between mode I fracture toughness and the tensile strength of rock[J]. International Journal of Rock Mechanics and Mining Sciences, 2016, 39(3): 401–406. |
[4] | Wei M D, Dai F, Xu N W, et al. Experimental and numerical study on the cracked chevron notched semi-circular bend method for characterizing the mode I fracture toughness of rocks[J]. Rock Mechanics and Rock Engineering, 2016, 49(5): 1595–1609. DOI:10.1007/s00603-015-0855-2 |
[5] | Zhu W C, Niu L L, Li S H, et al. Dynamic Brazilian test of rock under intermediate strain rate:pendulum hammer-driven SHPB test and numerical simulation[J]. Rock Mechanics and Rock Engineering, 2015, 48(5): 1867–1881. DOI:10.1007/s00603-014-0677-7 |
[6] |
周喻, 高永涛, 吴顺川.
等效晶质模型及岩石力学特征细观研究[J]. 岩石力学与工程学报, 2015, 34(3): 511–519.
( Zhou Yu, Gao Yong-tao, Wu Shun-chuan. An equivalent crystal model for mesoscopic behavior of rock[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(3): 511–519. ) |
[7] | Wong R H C, Chau K T, Wang P. Microcracking and grain size effect in Yuen long marbles[J]. International Journal of Rock Mechanics and Mining Sciences Geomechanic Abstract, 1996, 33(5): 479–485. DOI:10.1016/0148-9062(96)00007-1 |
[8] | Wang Y S, Hu X Z, Liang L. Determination of tensile strength and fracture toughness of concrete using notched 3-p-b specimens[J]. Engineering Fracture Mechanics, 2016, 160: 67–77. DOI:10.1016/j.engfracmech.2016.03.036 |
[9] | Mavko G, Mukerji T, Dvorkin J. The rock physics handbook:tools for seismic analysis in porous media[M]. Cambridge: Cambridge University Press, 1998. |
[10] | Hillerborg A. Application of the fctitious crack model to different types of materials[J]. International Journal of Fracture, 1991, 51(2): 95–102. |
[11] |
高峰, 谢和平, 赵鹏.
岩石块度分布的分形性质及细观结构效应[J]. 岩石力学与工程学报, 1994, 13(3): 240–246.
( Gao Feng, Xie He-ping, Zhao Peng. Fractal properties of size-frequency distribution of rock fragments and the influence of meso-structure[J]. Chinese Journal of Rock Mechanics and Engineering, 1994, 13(3): 240–246. ) |