Study on the effect of rock mesoscopic heterogeneity on its compressive mechanical properties based on GB-DDA method
GAO Shanhua1, 2, ZHANG Kaiyu1, 2, 3, YANG Mei1, 2, LIU Feng4, ZHU Kanyuan5
(1. Department of Hydraulic Science and Engineering, Taiyuan University of Technology, Taiyuan, Shanxi 030024, China; 2. Shanxi Key Laboratory of Collaborative Utilization of River Basin Water Resources, Taiyuan, Shanxi 030024, China; 3. China Institute for Radiation Protection, Taiyuan, Shanxi 030006, China; 4. Department of Architectural Engineering, Tianjin University, Tianjin 300350, China; 5. Chang Jiang Geotechnical Engineering Co., Ltd., Wuhan, Hubei 430010, China)
Abstract：Discontinuous deformation analysis (DDA) is an implicit discrete element method used to simulate the deformation of discrete block systems. To investigate the influence of mineral grains on the initiation, propagation, and coalescence of microcracks during rock loading, a Voronoi mineral grain model (GB-DDA) that adheres to statistical laws was established based on the DDA method. Compared to laboratory experiments, the GB-DDA, calibrated by meso parameters, can accurately simulate the mechanical response of Barre granite under tension and compression, enabling the quantitative analysis of intragranular and intergranular damage evolution processes. Using quasi-static uniaxial compression experiments, several GB-DDA numerical models with varying average grain sizes, grain roundness, and mineral content were developed to study the influence of mesoscopic heterogeneity on the compressive mechanical properties of Barre granite. Research findings indicate that during quasi-static uniaxial compression, the number of cracks in the numerical Barre granite samples, generated by different factors, follows the order: intergranular tensile cracks, transgranular tensile cracks, intergranular shear cracks, and transgranular shear cracks. The uniaxial compressive strength, elastic modulus, and Poisson's ratio of the rock are significantly affected by average grain size, grain roundness, and mineral content. An increase in average grain size and grain roundness promotes the initiation of shear cracks but has little effect on intergranular and transgranular cracks. Additionally, an increase in biotite content promotes the generation of intergranular cracks but has minimal impact on tensile and shear cracks.
高善铧1，2，张开雨1，2，3，杨 玫1，2，刘 丰4，朱看远5. 基于GB-DDA方法的岩石细观非均质性对其压缩力学特性影响研究[J]. 岩石力学与工程学报, 2025, 44(6): 1612-1623.
GAO Shanhua1, 2, ZHANG Kaiyu1, 2, 3, YANG Mei1, 2, LIU Feng4, ZHU Kanyuan5. Study on the effect of rock mesoscopic heterogeneity on its compressive mechanical properties based on GB-DDA method. , 2025, 44(6): 1612-1623.
YAO W，XU Y，XIA K. Damage evolution during rock pulverization induced by dynamic compressive loading[J]. Journal of Geophysical Research：Solid Earth，2020，125(5)：1-22.
[2]
LI S G，HUO R K，WANG B，et al. Experimental study on physicomechanical properties of sandstone under acidic environment[J]. Advances in Civil Engineering，2018，2018(3)：1-15.
[3]
ZHANG K Y，LIU F，XIA K W，et al. On the calibration and verification of Voronoi-based discontinuous deformation analysis for modeling rock fracture[J]. Journal of Rock Mechanics and Geotechnical Engineering，2023，15(8)：2 025-2 038.
[4]
CUNDALL P A. A computer model for simulating progressive large scale movements in blocky rock systems[C]// Proceedings of the International Symposium on Rock Mechanics. Nancy：[s. n.]，1971：129-136.
[5]
SHI G H，GOODMAN R E. Discontinuous deformation analysis-a new method for computing stress，strain and sliding of block systems[J]. Internal Medicine Journal，1988，44(3)：309-311.
[6]
TAN X，KONIETZKY H，CHEN W. Numerical Simulation of Heterogeneous Rock Using Discrete Element Model Based on Digital Image Processing[J]. Rock Mechanics and Rock Engineering，2016，49(12)：4 957-4 964.
[7]
GAO F Q，STEAD D，ELMO D. Numerical simulation of microstructure of brittle rock using a grain-breakable distinct element grain-based model[J]. Computers and Geotechnics，2016，78：203-217.
[8]
POTYONDY D O，CUNDALL P A. A bonded-particle model for rock[J]. International Journal of Rock Mechanics and Mining Sciences，2004，41(8 Special)：1 329-1 364.
[9]
ZHAO G F，FANG J N，ZHAO J. A 3D distinct lattice spring model for elasticity and dynamic failure[J]. International Journal for Numerical and Analytical Methods in Geomechanics，2011，35(8)：859-885.
[10]
WANG X，CAI M. Modeling of brittle rock failure considering inter- and intra-grain contact failures[J]. Computers and Geotechnics，2018，101：224-244.
[11]
殷坤龙，姜清辉，汪 洋. 新滩滑坡运动全过程的非连续变形分析与仿真模拟[J]. 岩石力学与工程学报，2002，21(7)：959-962.(YIN Kunlong，JIANG Qinghui，WANG Yang. Numerical simulation on the movement process of Xintan landslide by DDA method[J]. Chinese Journal of Rock Mechanics and Engineering，2002，21(7)：959-962.(in Chinese))
[12]
张国新，武晓峰. 裂隙渗流对岩石边坡稳定性的影响——渗流、变形耦合作用的DDA法[J]. 岩石力学与工程学报，2003，22(8)：1 269-1 275.(ZHANG Guoxin，WU Xiaofeng. Influence of seepage on the stability of rock slope—coupling of seepage and deformation by DDA method[J]. Chinese Journal of Rock Mechanics and Engineering，2003，22(8)：1 269-1 275.(in Chinese))
[13]
邬爱清，丁秀丽，李会中，等. 非连续变形分析方法模拟千将坪滑坡启动与滑坡全过程[J]. 岩石力学与工程学报，2006，25(7)：1 297-1 303.(WU Aiqing，DING Xiuli，LI Huizhong，et al. Numerical simulation of startup and whole failure process of Qianjiangping landslide using discontinuous deformation analysis[J]. Chinese Journal of Rock Mechanics and Engineering，2006，25(7)：1 297-1 303.(in Chinese))
[14]
JIAO Y Y，HUANG G H，ZHAO Z Y，et al. An improved three-dimensional spherical DDA model for simulating rock failure[J]. Chinese Academy of Sciences，2015，58(9)：1 533-1 541.
[15]
BEYABANAKI S A R，JAFARI A，BIABANAKI S O R，et al. Nodal-based three-dimensional discontinuous deformation analysis (3-D DDA)[J]. Computers and Geotechnics，2009，36(3)：359-372.
[16]
BAO H R. Nodal-based discontinuous deformation analysis[Ph. D. Thesis][D]. Singapore：Nanyang Technological University，2010.
[17]
LIN C T，AMADEI B，JUNG J，et al. Extensions of discontinuous deformation analysis for jointed rock masses[J]. International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts，1996，33(7)：671-694.
[18]
PEARCE C J，THAVALINGAM A，LIAO Z. Computational aspects of the discontinuous deformation analysis framework for modelling concrete fracture[J]. Engineering Fracture Mechanics，2000，65(2/3)：283-298.
[19]
NING Y J，YANG J，AN X M，et al. Modelling rock fracturing and blast-induced rock mass failure via advanced discretisation within the discontinuous deformation analysis framework[J]. Computers and Geotechnics，2011，38(1)：40.
[20]
NI K，YANG J，NING Y J，et al. A modified sub-block DDA fracturing modelling method for rock[J]. Engineering Analysis with Boundary Elements，2020，111：154.
[21]
LI C Y，NING Y J，LIU X L. Discontinuous deformation analysis (DDA) simulations of crack propagation，branching and coalescence in brittle materials under dynamic loading[J]. Theoretical and Applied Fracture Mechanics，2022，122：103648.
[22]
JIAO Y Y，ZHANG X L，ZHAO J. Two-dimensional DDA contact constitutive model for simulating rock fragmentation[J]. Journal of Engineering Mechanics，2012，138(2)：199-209.
[23]
张开雨，夏开文，刘 丰. 基于Voronoi多边形离散的DDA方法模拟岩石破坏[J]. 岩石力学与工程学报，2021，40(4)：725-738. (ZHANG Kaiyu，XIA Kaiwen，LIU Feng. Simulation of rock failure by Voronoi-based discontinuous deformation analysis[J]. Chinese Journal of Rock Mechanics and Engineering，2021，40(4)：725-738.(in Chinese))
[24]
PAN C，LI X，LI J C，et al. Numerical investigation of blast-induced fractures in granite：insights from a hybrid LS-DYNA and UDEC grain-based discrete element method[J]. Geomechanics and Geophysics for Geo-Energy and Geo-Resources，2021，7(2)：1-18.
[25]
MORGAN S P，JOHNSON C A，EINSTEIN H H. Cracking processes in Barre granite： fracture process zones and crack coalescence[J]. International Journal of Fracture，2013，180(2)：177-204.
[26]
LI B Q，CASANOVA M，EINSTEIN H H. Laboratory study of fracture initiation and propagation in Barre granite under fluid pressure at stress state close to failure[J]. Journal of Rock Mechanics and Geotechnical Engineering，2023，15(3)：538-550.
BANADAKI M M D. Stress-wave induced fracture in rock due to explosive action[Ph. D. Thesis][D]. Toronto：University of Toronto，2011.
[29]
LAN H X，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(B1)：B01202.
[30]
PENG J，WONG L N Y，TEH C L. Effects of grain size-to-particle size ratio on micro-cracking behavior using a bonded-particle grain-based model[J]. International Journal of Rock Mechanics and Mining Sciences，2017，100：207-217.
[31]
张世瑞，邱士利，李邵军，等. 北山花岗岩细观非均质性对单轴压缩力学特性影响的FDEM数值研究[J]. 岩石力学与工程学报，2022，41(增1)：2 658-2 672.(ZHANG Shirui，QIU Shili，LI Shaojun，et al. Study on mechanical properties of Beishan granite meso-heterogeneity under uniaxial compression with FDEM[J]. Chinese Journal of Rock Mechanics and Engineering，2022，41(Supp.1)：2 658-2 672.(in Chinese))
[32]
LI L，LEE P K K，TSUI Y，et al. Failure process of granite[J]. International Journal of Geomechanics，2003，3(1)：84-98.
[33]
梁昌玉，吴树仁，李 晓. 中低应变率范围内单轴压缩下花岗岩断口细-微观特征研究[J]. 岩石力学与工程学报，2015，34(增1)：2 977-2 985.(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 985.(in Chinese))
[34]
ZHOU J，ZHANG L Q，YANG D X，et al. Investigation of the quasi-brittle failure of Alashan granite viewed from laboratory experiments and grain-based discrete element modeling[J]. Materials，2017，10(7)：835.
[35]
PENG J，WONG L N Y，TEH C I. Influence of grain size on strength of polymineralic crystalline rock：New insights from DEM grain-based modeling[J]. Journal of Rock Mechanics and Geotechnical Engineering，2021，13(4)：755-766.
[36]
王穗丰，赵先宇，张利平，等. 基于晶粒织构模型的花岗岩矿物晶粒形状及朝向对其力学特性影响研究[J]. 岩石力学与工程学报，2023，42(10)：2 515-2 527.(WANG Suifeng，ZHAO Xianyu，ZHANG Liping，et al. Influence of grain shape and orientation of granite on mechanical property based on the grain texture model[J]. Chinese Journal of Rock Mechanics and Engineering，2023，42(10)：2 515-2 527.(in Chinese))
[37]
PAN C，ZHAO G M，MENG X R，et al. Numerical investigation of the influence of mineral mesostructure on quasi-static compressive behaviors of granite using a breakable grain-based model[J]. Frontiers in Ecology and Evolution，2023，11：1288870.
[38]
YESILOGLU-GULTEKIN N，SEZER E A，GOKCEOGLU C，et al. An application of adaptive neuro fuzzy inference system for estimating the uniaxial compressive strength of certain granitic rocks from their mineral contents[J]. Expert Systems with Applications，2013，40(3)：921-928.
[39]
SAJID M，ARIF M. Reliance of physico-mechanical properties on petrographic characteristics：consequences from the study of Utla granites，north-west Pakistan[J]. Bulletin of Engineering Geology and the Environment，2015，74(4)：1 321-1 330.
2025, Vol. 44 
