Research on rock crack contact model considering linked substances based on particle flow method

Abstract

The models constructed by particle flow simulation method can effectively simulate the heterogeneous substance characteristics and failure behaviors of rocks. However, existing contact models overlook the rock cracks, and the various simulation methods that do consider cracks still exhibit certain limitations. In this paper, based on Flat-Joint model and Linear Parallel Bond model, a crack contact model considering linked substance in the crack is proposed by splitting the crack contact into two portions: linked portion and unlinked portion for calculation. The new contact model considers the influence of crack closure on the contact force-displacement law. And a better compressive tensile strength ratio (UCS/T) was obtained by limiting the failure of the contact bond to be solely controlled by the contact force and moment of the linked portion. Then, by employing the FISH Model tool within the Particle Flow Code, the contact model was constructed and verified through contact force-displacement experiments and loading-unloading tests with cracked model. Finally, the contact model was tested through simulations of rock mechanics experiments. The results indicate that the contact model can effectively simulate the axial and lateral strain laws of rocks simultaneously and has a relatively good reproduction of the bi-modularity of rocks.

Keywords

Crack closure / Contact model / Compaction stage / Initial effective modulus / Lateral strain

Acknowledgement

This study was supported by the Natural Science Foundation of Heilongjiang Province (No. ZD2021E006) and the National Natural Science Foundation of China (Nos. 52174075 and 52074110).

Supplementary material

Supplementary data to this article can be found online at https://doi.org/10.1016/j.ijmst.2025.02.007.

References

[1] Walsh JB. The effect of cracks on the uniaxial elastic compression of rocks. J Geophys Res 1965; 70(2):399-411.

[2] Walsh JB. The effect of cracks on the compressibility of rock. J Geophys Res 1965; 70(2):381-9.

[3] Bristow JR. Microcracks, and the static and dynamic elastic constants of annealed and heavily cold-worked metals. Br J Appl Phys 1960; 11(2):81-5.

[4] Jaeger JC, Cook NGW, Zimmerman RW. Fundamentals of Rock Mechanics. 4th ed. Oxford: Blackwell; 2007.

[5] David EC, Brantut N, Schubnel A, Zimmerman RW. Sliding crack model for nonlinearity and hysteresis in the uniaxial stress-strain curve of rock. Int J Rock Mech Min Sci 2012;52:9-17.

[6] David EC, Brantut N, Hirth G. Sliding crack model for nonlinearity and hysteresis in the triaxial stress-strain curve of rock, and application to antigorite deformation. J Geophys Res Solid Earth 2020; 125(10): e2019JB018970.

[7] Batzle ML, Simmons G, Siegfried RW. Microcrack closure in rocks under stress: Direct observation. J Geophys Res Solid Earth 1980; 85(B12):7072-90.

[8] Zhou HW, Zhong JC, Ren WG, Wang XY, Yi HY. Characterization of porefracture networks and their evolution at various measurement scales in coal samples using X-ray lCT and a fractal method. Int J Coal Geol 2018;189:35-49.

[9] Meng QS, Wu K, Zhou HR, Qin QL, Wang C. Mesoscopic damage evolution of coral reef limestone based on real-time CT scanning. Eng Geol 2022;307:106781.

[10] Sufian A, Russell AR. Microstructural pore changes and energy dissipation in Gosford sandstone during pre-failure loading using X-ray CT. Int J Rock Mech Min Sci 2013;57:119-31.

[11] Zhang T, Xie ZZ, Yuan L, Tang M, Zheng KG, Zhang C, Mao JL. Study on coal damage and pore-fracture evolution ahead of the mining face under the disturbance of coal mining. J China Coal Soc 2024:1-16. in Chinese.

[12] Meng XX, Qin YY, Zhuge FM, Liu WT. Microstructure characteristics and the damage evolution rule of sandstone under triaxial compression test using Nuclear Magnetic Resonance (NMR). Mater Lett 2022;309:131355.

[13] Bi J, Ning L, Zhao Y, Wu ZJ, Wang CL. Analysis of the microscopic evolution of rock damage based on real-time nuclear magnetic resonance. Rock Mech Rock Eng 2023; 56(5):3399-411.

[14] Cao WG, Tan X, Zhang C, He M. Constitutive model to simulate full deformation and failure process for rocks considering initial compression and residual strength behaviors. Can Geotech J 2019; 56(5):649-61.

[15] Ren CH, Yu J, Liu XY, Zhang ZQ, Cai YY. Cyclic constitutive equations of rock with coupled damage induced by compaction and cracking. Int J Min Sci Technol 2022; 32(5):1153-65.

[16] Wang CP, Liu JF, Chen L, Liu J, Wang L, Liao YL. Creep constitutive model considering nonlinear creep degradation of fractured rock. Int J Min Sci Technol 2024; 34(1):105-16.

[17] Li LC, Xia YJ, Huang B, Zhang LY, Li M, Li AS. The behaviour of fracture growth in sedimentary rocks: A numerical study based on hydraulic fracturing processes. Energies 2016;9:169.

[18] Li ZC, Li LC, Huang B, Zhang LY, Li M, Zuo JQ, Li AS, Yu Q. Numerical investigation on the propagation behavior of hydraulic fractures in shale reservoir based on the DIP technique. J Petrol Sci Eng 2017;154:302-14.

[19] Zhang HQ, Shi H, Jing HW, Wu Y, Pu H. Numerical study of remote fracturing around a circular opening in rock. Eur J Environ Civ Eng 2020; 24(7):1032-50.

[20] Wang J, Liu J, Ouyang XS, Zou HB. Analysis of the fracture characteristics and the stress intensity factor of the symmetric coalescence crack in rock. Geotech Geol Eng 2019; 37(4):2535-44.

[21] Shi H, Zhang HQ, Song L, Yang Z, Yuan GT, Xue XR, Wang Y. Failure characteristics of sandstone specimens with randomly distributed pre-cracks under uniaxial compression. Environ Earth Sci 2020; 79(9):193.

[22] Wu N, Liang ZZ, Zhou JR, Zhang LZ. Energy evolution characteristics of coal specimens with preformed holes under uniaxial compression. Geomech Eng 2021;20:55-66.

[23] Ji PQ, Zhang XP, Zhang Q. A new method to model the non-linear crack closure behavior of rocks under uniaxial compression. Int J Rock Mech Min Sci 2018;112:171-83.

[24] Patel S, Martin CD. Impact of the initial crack volume on the intact behavior of a bonded particle model. Comput Geotech 2020;127:103764.

[25] Zhang XP, Ji PQ, Peng J, Wu SC, Zhang Q. A grain-based model considering preexisting cracks for modelling mechanical properties of crystalline rock. Comput Geotech 2020;127:103776.

[26] Chen HG, Zhao C, Zhang R, Huang L, Pan HY, Qian Y, Yang HW. Micromechanics-based model of rock considering compaction deformation: Analysis of microcrack influence from discrete and continuous perspectives. Comput Geotech 2023;164:105766.

[27] Peng SY, Li XB, Li CJ, Liang LS, Huang LQ. Crack-closure behavior and stresssensitive wave velocity of hard rock based on flat-joint model in particle-flowcode (PFC) modeling. Comput Geotech 2024;170:106320.

[28] Zhang BY, Zhang CS, Wang CL, Hao BQ. Calibration method of mesoparameters of PFC2D flat-joint model. Chin J Comput Mech 2021; 38 (5):665-73. in Chinese.

[29] Lum JE, Viljoen F, Cairncross B, Frei D. Mineralogical and geochemical characteristics of BERYL (AQUAMARINE) from the erongo volcanic complex. Namibia J Afr Earth Sci 2016;124:104-25.

[30] Li Y, Chen JQ, Elsworth D, Pan ZJ, Ma XT. Nanoscale mechanical property variations concerning mineral composition and contact of marine shale. Geosci Front 2022; 13(4):101405.

[31] Goral J, Walton I, Andrew M, Deo M. Pore system characterization of organicrich shales using nanoscale-resolution 3D imaging. Fuel 2019;258:116049.

[32] Li Y, Yang JH, Pan ZJ, Tong WS. Nanoscale pore structure and mechanical property analysis of coal: An insight combining AFM and SEM images. Fuel 2020;260:116352.

[33] Zhao DF, Guo YH, Wang G, Mao XX. Characterizing nanoscale pores and its structure in coal: Experimental investigation. Energy Explor Exploit 2019; 37 (4):1320-47.

[34] Potyondy DO. A Flat-Jointed Bonded-Particle Material For Hard Rock. In:Proceedings of the 46th U.S. Rock Mechanics/Geomechanics Symposium. Chicago: American Rock Mechanics Association; 2012. p. 15-501.

[35] Itasca. Linear Parallel Bond Model, 2024.

[36] Itasca. Contact Resolution, 2024.

[37] Yoon J. Application of experimental design and optimization to PFC model calibration in uniaxial compression simulation. Int J Rock Mech Min Sci 2007; 44(6):871-89.

[38] Tang HD, Zhu ML. Timely spread characteristics of micro- and meso-scale fractures based on SEM experiment. Environ Earth Sci 2020; 79(20):484.

[39] Huang YY, Li TB. On the relationship between the mechanical strength and mineral microstructure of biotite granite under temperature influence. Front Earth Sci 2023;11:1209433.

[40] Huang XL, Qi SW, Guo SF, Zheng BW, Zhao Q, Sha P, Wang TZ, Yao XL, Liang N. Effect of the crystal habit on micromechanical extensile behaviors of the crystalline rock during compression. Eng Geol 2022;310:106874.

[41] Ding XB, Zhang LY. A new contact model to improve the simulated ratio of unconfined compressive strength to tensile strength in bonded particle models. Int J Rock Mech Min Sci 2014;69:111-9.

[42] Xu ZH, Wang WY, Lin P, Xiong Y, Liu ZY, He S. A parameter calibration method for PFC simulation: Development and a case study of limestone. Geomech Eng 2020;22:97-108.

[43] Bass JD. In:Elasticity of minerals, glasses, and melts. Mineral Physics & Crystallography. Washington: American Geophysical Union; 2013. p. 45-63.

[44] Li XF, Zhang QB, Li HB, Zhao J. Grain-based discrete element method (GB-DEM) modelling of multi-scale fracturing in rocks under dynamic loading. Rock Mech Rock Eng 2018; 51(12):3785-817.

[45] Li YJ, Wang BQ, Ren ZX, Zuo JP, Liu SX, Liu DJ. Effect of microstructure heterogeneity on shale strength: Insights from micro-CT imaging and numerical simulations. Rock Mech Rock Eng 2025; 58(1):225-44.

[46] Cheng Y, Wong LNY. A study on mechanical properties and fracturing behavior of Carrara marble with the flat-jointed model. Int J Numer Anal Meth Geomech 2020; 44(6):803-22.

[47] Chen CC, Xu J, Okubo S, Peng SJ. Damage evolution of tuff under cyclic tension- compression loading based on 3D digital image correlation. Eng Geol 2020;275:105736.

[48] Wu GZ, Wang WS, Peng SC. Analytical solution of the stress field and plastic zone at the tip of a closed crack. Front Earth Sci 2024;12:1370672.

[49] Patel S, Martin CD. Evaluation of tensile Young’s modulus and Poisson’s ratio of a bi-modular rock from the displacement measurements in a Brazilian test. Rock Mech Rock Eng 2018; 51(2):361-73.

[50] Cao MH, Yang SQ, Tian WL, Huang YH, Huang M. Three-dimensional discrete element simulation of the triaxial cyclic loading of sandstone based on a nonlinear parallel-bonded stress corrosion model. Comput Geotech 2024;168:106167.