基于SPH的花岗岩高精度建模方法及矿物尺度裂纹扩展规律研究

doi:10.3724/1000-6915.jrme.2025.0857

摘要
图/表
参考文献
相关文章 (15)

全文: PDF (4750 KB) HTML (1 KB)
输出: BibTeX | EndNote (RIS)

摘要准确预测岩石材料的宏观力学行为，其关键在于深刻理解其矿物尺度的微裂纹演化规律及机制。提出一种高精度的多尺度建模框架，该框架耦合高精度X射线计算机断层扫描（CT）技术、纳米压痕技术与光滑粒子流数值方法（SPH）。通过该建模框架，成功构建能够精确反映花岗岩内部矿物晶体（石英、长石、云母）及晶界空间分布的SPH-GBM数值模型。模拟结果表明：花岗岩在单轴加载下的破坏主要由拉伸损伤驱动。微裂纹首先在易产生应力集中的晶界处以拉伸形式（晶界裂纹）萌生；随着荷载增加，裂纹沿晶界扩展并汇合，随后穿过晶体内部（晶内裂纹）。裂纹路径具有明显的矿物选择性，微裂纹优先在长石内部生成，其次是石英，而刚度最小的云母颗粒则在接近峰值强度时才发生断裂。研究证实了SPH-GBM作为一种连接岩石高精度微观结构与宏观失效行为的强大工具，在岩石细观力学研究中具有广阔的应用前景。

	服务

	把本文推荐给朋友
	加入我的书架
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	周宇1
	2
	3
	钟露露1
	郁舒阳4*
	梁秦源5
	吕文俊6
	韩观胜1
	2
	唐琼琼1
	2
	李博7
	伍法权1
	2
	3

关键词 ：岩石力学, 矿物尺度, SPH, 高精度建模方法, 单轴压缩试验, 裂纹扩展

Abstract：The key to accurately predicting the macro-mechanical behavior of rock materials lies in a profound understanding of the evolution laws and mechanisms governing microcracks at the mineral scale. This paper proposes a high-precision multi-scale modeling framework that integrates advanced X-ray computed tomography (CT) technology, nanoindentation techniques, and the Smoothed Particle Hydrodynamics (SPH) numerical method. Utilizing this modeling framework, an SPH-GBM numerical model that accurately represents the spatial distribution of mineral crystals (quartz, feldspar, and mica) and grain boundaries within granite was successfully developed. The simulation results indicate that the failure of granite under uniaxial loading is primarily driven by tensile damage. Microcracks initially emerge as tensile intragranular cracks at grain boundaries, where stress concentrations are likely to occur. As the load increases, these cracks propagate and converge along the grain boundaries and subsequently penetrate the interior of the crystals, forming grain boundary cracks. The crack propagation path exhibits distinct mineral selectivity, with microcracks preferentially forming within feldspar, followed by quartz, while mica particles, which possess the lowest stiffness, tend to fracture near peak strength. This study confirms that SPH-GBM is a robust tool for linking high-precision microstructural characteristics with macroscopic failure behavior in rocks, offering significant potential for advancing research in rock mesomechanics.

Key words： rock mechanics mineral-scale SPH high-precision modeling methods uniaxial compression test crack propagation

引用本文:

周宇1，2，3，钟露露1，郁舒阳4*，梁秦源5，吕文俊6，韩观胜1，2，唐琼琼1，2，李博7，伍法权1，2，3. 基于SPH的花岗岩高精度建模方法及矿物尺度裂纹扩展规律研究[J]. 岩石力学与工程学报, 2026, 45(10): 3016-3030.
ZHOU Yu1, 2, 3, ZHONG Lulu1, YU Shuyang4*, LIANG Qinyuan5, LYU Wenjun6, HAN Guansheng1, 2, TANG Qiongqiong1, 2, LI Bo7, WU Faquan1, 2, 3. Crack propagation law at the mineral scale of granite based on the high-precision modeling approach using SPH. , 2026, 45(10): 3016-3030.

链接本文:

https://rockmech.whrsm.ac.cn/CN/10.3724/1000-6915.jrme.2025.0857 或 https://rockmech.whrsm.ac.cn/CN/Y2026/V45/I10/3016

参考文献(References)：

[1] TIAN J J，XU D J，LIU T H. An experimental investigation of the fracturing behaviour of rock-like materials containing two V-shaped parallelogram flaws[J]. International Journal of Mining Science and Technology，2020，30(6)：777–783.

[2] 张睿，赵程，幸金权，等. 基于平节理模型的花岗岩固有微裂纹对裂纹演化影响分析[J]. 岩石力学与工程学报，2021，40(9)：1 857–1 867.(ZHANG Rui，ZHAO Cheng，XING Jinquan，et al. Influence of granite inherent microcracks on crack evolution based on the flat-joint contact model[J]. Chinese Journal of Rock Mechanics and Engineering，2021，40(9)：1 857–1 867.( in Chinese))

[3] 张国凯，李海波，王明洋，等. 单裂隙花岗岩破坏强度及裂纹扩展特征研究[J]. 岩石力学与工程学报，2019，38(增1)：2 760–2 771. (ZHANG Guokai，LI Haibo，WANG Mingyang，et al. Study on characteristics of failure strength and crack propagation of granite rocks containing a single fissure[J]. Chinese Journal of Rock Mechanics and Engineering，2019，38(Supp.1)：2 760–2 771.( in Chinese))

[4] 李博，梁秦源，周宇，等. 基于CT-GBM重构法的花岗岩裂纹扩展规律研究[J]. 岩石力学与工程学报，2022，41(6)：1 114–1 125. (LI Bo，LIANG Qinyuan，ZHOU Yu，et al. Research on crack propagation law of granite based on CT-GBM reconstruction method[J]. Chinese Journal of Rock Mechanics and Engineering，2022，41(6)：1 114–1 125.( in Chinese))

[5] 贾蓬，杨其要，刘冬桥，等. 高温花岗岩水冷却后物理力学特性及微观破裂特征[J]. 岩土力学，2021，42(6)：1 568–1 578.(JIA Peng，YANG Qiyao，LIU Dongqiao，et al. Physical and mechanical properties and related microscopic characteristics of high-temperature granite after water-cooling [J]. Rock and Soil Mechanics，2021，42(6)：1 568–1 578. (in Chinese))

[6] 梁昌玉，吴树仁，李晓. 中低应变率范围内单轴压缩下花岗岩断口细–微观特征研究[J]. 岩石力学与工程学报，2015，34(增1)：2 977–2 986.(LIANG Changyu，WU Shuren，LI Xiao. Research on micro-meso characteristics of granite fracture under uniaxial compression at low and intermediate strain rates[J]. Chinese Journal of Rock Mechanics and Engineering，2015，34(Supp.1)：2 977–2 986.(in Chinese))

[7] 张志镇，高峰，刘治军. 温度影响下花岗岩冲击倾向及其微细观机制研究[J]. 岩石力学与工程学报，2010，29(8)：1 591–1 602. (ZHANG Zhizhen，GAO Feng，LIU Zhijun. Research on rockburst proneness and its microcosmic mechanism of granite considering temperature effect[J]. Chinese Journal of Rock Mechanics and Engineering，2010，29(8)：1 591–1 602.( in Chinese))

[8] XUE D J，ZHOU H W，ZHAO Y W，et al. Real-time SEM observation of mesoscale failures under thermal-mechanical coupling sequences in granite[J]. International Journal of Rock Mechanics and Mining Sciences，2018，112：35–46.

[9] JING L. A review of techniques，advances and outstanding issues in numerical modelling for rock mechanics and rock engineering[J]. International Journal of Rock Mechanics and Mining Sciences，2003，40(3)：283–353.

[10] TANG C A，ZHANG Y B，LIANG Z Z，et al. Fracture spacing in layered materials and pattern transition from parallel to polygonal fractures[J]. Physical Review E，2006，73(5)：056120.

[11] TANG C A，HUDSON J A. Rock failure mechanisms：illustrated and explained[M]. [S. l.]：CRC Press，2010：58–203.

[12] POTYONDY D O. A grain-based model for rock：approaching the true microstructure[C]// Proceedings of rock mechanics in the Nordic Countries. [S. l.]：[s. n.]，2010：9–12.

[13] KONG L，SHANG J，RANJITH P G，et al. Grain-based DEM modelling of mechanical and coupled hydro-mechanical behaviour of crystalline rocks[J]. Engineering Geology，2024，339：107649.

[14] ZHANG Y H，WONG L N Y. A review of numerical techniques approaching microstructures of crystalline rocks[J]. Computers & Geosciences，2018，115：167–187.

[15] HU X J，GUO P P，XIE N，et al. Creep behavior and associated acoustic characteristics of heterogeneous granite containing a single pre-existing flaw using a grain-based parallel-bonded stress corrosion model[J]. Rock Mechanics and Rock Engineering，2023，56(7)：4 799–4 832.

[16] HU X J，XIE N，ZHU Q Z，et al. Modeling damage evolution in heterogeneous granite using digital image-based grain-based model[J]. Rock Mechanics and Rock Engineering，2020，53(11)：4 925–4 945.

[17] ZHENG Z，LI S X，QIN Z Y，et al. Dual-scale 3DEC-GBM Discrete element simulation on mechanical behavior and anisotropic fracture evolution mechanism of rock induced by true three-dimensional stress[J]. Rock Mechanics and Rock Engineering，2024，57(5)：3 885–3 915.

[18] 王桂林，王润秋，孙帆. 块体离散元颗粒模型细观参数标定方法及花岗岩细观演化模拟[J]. 长江科学院院报，2022，39(1)：86–93.(WANG Guilin，WANG Runqiu，SUN Fan. A discrete element GBM simulation method for meso?parameter calibration and granite meso?evolution simulation[J]. Journal of Yangtze River Scientific Research Institute，2022，39(1)：86–93.(in Chinese))

[19] LAN H，MARTIN C D，ANDERSSON J C. Evolution of in situ rock mass damage induced by mechanical–thermal loading[J]. Rock Mechanics and Rock Engineering，2013，46(1)：153–168.

[20] 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.

[21] 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.

[22] ZHOU Y，LV W J，ZHANG C，et al. Novel hard rock breaking technique using ultra-high-frequency particle impact induced by ultrasonic vibration field[J]. Energy，2024，288(C)：129747.

[23] ZHOU Y，LV W J，LI B，et al. A realistic 3D grain-based modeling approach for reproducing the mechanical and failure behavior of brittle granites[J]. International Journal of Rock Mechanics and Mining Sciences，2025，185：105981.

[24] ZHOU Y，LV W J，LI B，et al. Impact of soft minerals on crack propagation in crystalline rocks under uniaxial compression：A grain‐based numerical investigation[J]. International Journal for Numerical and Analytical Methods in Geomechanics，2024，48(8)：2 020–2 042.

[25] DAS R，CLEARY P W. Effect of rock shapes on brittle fracture using Smoothed Particle Hydrodynamics[J]. Theoretical and Applied Fracture Mechanics，2010，53(1)：47–60.

[26] ZHOU X P，ZHAO Y，QIAN Q H. A novel meshless numerical method for modeling progressive failure processes of slopes[J]. Engineering Geology，2015，192：139–153.

[27] ZHAO Y，ZHOU X P，QIAN Q H. Progressive failure processes of reinforced slopes based on general particle dynamic method[J]. Journal of Central South University，2015，22(10)：4 049–4 055.

[28] 周小平，赵毅，钱七虎. 单轴压缩条件下岩石破坏的光滑粒子流体动力学数值模拟[J]. 岩石力学与工程学报，2015，34(增1)：2 647– 2 658.(ZHOU Xiaoping，ZHAO Yi，QIAN Qihu. Numerical simulation of rock failure process in uniaxial compression using smoothed particle hydrodynamics[J]. Chinese Journal of Rock Mechanics and Engineering，2015，34(Supp.1)：2 647–2 658.(in Chinese))

[29] BI J，ZHOU X P. Numerical simulation of zonal disintegration of the surrounding rock masses around a deep circular tunnel under dynamic unloading[J]. International Journal of Computational Methods，2015，12(3)：1550020.

[30] BI J，ZHOU X P，QIAN Q H. The 3D numerical simulation for the propagation process of multiple pre-existing flaws in rock-like materials subjected to biaxial compressive loads[J]. Rock Mechanics and Rock Engineering，2016，49(5)：1 611–1 627.

[31] HU M，TAN Q T，FENG D L，et al. Simulation of rock crack propagation and failure behavior based on a mixed failure model with SPH[J]. Rock Mechanics and Rock Engineering，2024，57(11)：9 575–9 596.

[32] YU S Y，REN X H，ZHANG J X，et al. An improved form of smoothed particle hydrodynamics method for crack propagation simulation applied in rock mechanics[J]. International Journal of Mining Science and Technology，2021，31(3)：421–428.

[33] MA G W，WANG X J，REN F. Numerical simulation of compressive failure of heterogeneous rock-like materials using SPH method[J]. International Journal of Rock Mechanics and Mining Sciences，2011，48(3)：353–363.

[34] GHAREHDASH S，SAINSBURY B A L，BARZEGAR M，et al. Field scale modelling of explosion-generated crack densities in granitic rocks using dual-support smoothed particle hydrodynamics(DS-SPH)[J]. Rock Mechanics and Rock Engineering，2021，54(9)：4 419–4 454.

[35] LIBERSKY L D，PETSCHEK A G. Smooth particle hydrodynamics with strength of materials[C]// Advances in the Free-Lagrange Method Including Contributions on Adaptive Gridding and the Smooth Particle Hydrodynamics Method. Berlin Heidelberg：Springer，1990：248–257.

[36] LIBERSKY L D，PETSCHEK A G，CARNEY T C，et al. High strain Lagrangian hydrodynamics：A three-dimensional SPH code for dynamic material response[J]. Journal of Computational Physics，1993，109(1)：67–75.

[37] ZHOU X P，YAO W W. Smoothed bond-based peridynamics[J]. Journal of Peridynamics and Nonlocal Modeling，2022，4(4)：452–474.

[38] MONAGHAN J J. Simulating free surface flows with SPH[J]. Journal of Computational Physics，1994，110(2)：399–406.

[39] BI J，ZHOU X P，XU X M. Numerical simulation of failure process of rock-like materials subjected to impact loads[J]. International Journal of Geomechanics，2017，17(3)：04016073.

[40] ZHOU X P，BI J. Numerical simulation of thermal cracking in rocks based on general particle dynamics[J]. Journal of Engineering Mechanics，2018，144(1)：04017156.

[41] ZHANG Z Y，JIN X G，LUO W. Numerical study on the collapse behaviors of shallow tunnel faces under open-face excavation condition using mesh-free method[J]. Journal of Engineering Mechanics，2019，145(11)：04019085.

[42] 杨振琦，王述红，孟嫣然，等. 基于CT扫描的花岗岩三维数值试件重构模型及应用[J]. 固体力学学报，2017，38(6)：591–600. (YANG Zhenqi，WANG Shuhong，MENG Yanran，et al. The 3D numerical specimen of granite based on CT technology and application[J]. Journal of Solid Mechanics，2017，38(6)：591–600.(in Chinese))

[43] 尧俊凯，陈晓斌，蔡德钩，等. 基于X-ray CT粗粒土填料细观结构表征分析[J]. 铁道建筑，2021，61(5)：70–74.(YAO Junkai，CHEN Xiaobin，CAI Degou，et al. Characterization analysis of meso?structure of coarse?grained soil filling based on X?Ray CT[J]. Railway Engineering，2021，61(5)：70–74.(in Chinese))

[44] KOMADJA G C，STANISLAS T T，MUNGANYINKA P，et al. New approach for assessing uniaxial compressive strength of rocks using measurement from nanoindentation experiments[J]. Bulletin of Engineering Geology and the Environment，2022，81(8)：299.

[45] 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.

[46] GHAZVINIAN E，DIEDERICHS M S，QUEY R. 3D random Voronoi grain-based models for simulation of brittle rock damage and fabric-guided micro-fracturing[J]. Journal of Rock Mechanics and Geotechnical Engineering，2014，6(6)：506–521.

[47] ZHANG X P，JI P Q，PENG J，et al. A grain-based model considering pre-existing cracks for modelling mechanical properties of crystalline rock[J]. Computers and Geotechnics，2020，127：103776.

[48] ZHOU W，JI X，MA G，et al. FDEM Simulation of Rocks with Microstructure Generated by Voronoi Grain-Based Model with Particle Growth[J]. Rock Mechanics and Rock Engineering，2020，53(4)：1 909–1 921.

[49] ABDELAZIZ A，ZHAO Q，GRASSELLI G. Grain based modelling of rocks using the combined finite-discrete element method[J]. Computers and Geotechnics，2018，103：73–81.

[50] LI X F，LI H B，LIU L W，et al. Investigating the crack initiation and propagation mechanism in brittle rocks using grain-based finite-discrete element method[J]. International Journal of Rock Mechanics and Mining Sciences，2020，127：104219.

[51] MÜLLER C，FRÜHWIRT T，HAASE D，et al. Modeling deformation and damage of rock salt using the discrete element method[J]. International Journal of Rock Mechanics and Mining Sciences，2018，103：230–241.

[52] GUI Y L，ZHAO Z Y，JI J，et al. The grain effect of intact rock modelling using discrete element method with Voronoi grains[J]. Géotechnique Letters，2016，6(2)：136–143.

[53] LI H，YANG J，HAN Y，et al. Weibull grain-based model (W-GBM) for simulating heterogeneous mechanical characteristics of salt rock[J]. Engineering Analysis with Boundary Elements，2019，108：227–243.

[54] QUEY R，DAWSON P R，BARBE F. Large-scale 3D random polycrystals for the finite element method：Generation，meshing and remeshing[J]. Computer Methods in Applied Mechanics and Engineering，2011，200(17–20)：1 729–1 745.

[55] 谢和平，鞠杨，黎立云，等. 岩体变形破坏过程的能量机制[J]. 岩石力学与工程学报，2008，27(9)：1 729–1 740.(XIE Heping，JU Yang，LI Liyun, et al. Energy mechanism of deformation and failure of rock masses[J]. Chinese Journal of Rock Mechanics and Engineering，2008，27(9)：1 729–1 740.(in Chinese))

[56] MIKHALYUK A V，ZAKHAROV V V. Dissipation of dynamic-loading energy in quasi-elastic deformation processes in rocks[J]. Journal of Applied Mechanics and Technical Physics，1997，38(2)：312–318.

[57] RECHES Z，LOCKNER D A. Nucleation and growth of faults in brittle rocks[J]. Journal of Geophysical Research：Solid Earth，1994，99(B9)：18 159–18 173.

[58] LEE S E，KIM M I，PARK J H，et al. Damage process of intact granite under uniaxial compression：microscopic observations and contact stress analysis of grains[J]. Geosciences Journal，2006，10(4)：457–463.

[59] CUI Z D，ZHANG J Y，RANJITH P G，et al. Crystallographic constraints on microcrack dynamics in granite：Insights from in-situ microscopic characterization of crack network propagation[J]. International Journal of Rock Mechanics and Mining Sciences，2025，195(1)：106278.

[60] XIN X K，MENG Q B，PU H，et al. Numerical simulation analysis of energy evolution law and distribution characteristics of anchored body unidirectional loading[J]. Tunnelling and Underground Space Technology，2025，158：106419.