Hybrid finite-discrete element simulator based on GPGPU-parallelized computation for modelling crack initiation and coalescence in sandy mudstone with prefabricated cross-flaws under uniaxial compression

被引:31
作者
Liang, Dongxu [1 ,2 ]
Zhang, Nong [1 ,3 ]
Liu, Hongyuan [2 ]
Fukuda, Daisuke [4 ]
Rong, Haoyu [1 ]
机构
[1] China Univ Min & Technol, Sch Mines, Key Lab Deep Coal Resource Min, Minist Educ China, Xuzhou 221116, Jiangsu, Peoples R China
[2] Univ Tasmania, Coll Sci & Engn, Hobart, Tas 7001, Australia
[3] Jiangsu Normal Univ, Open Lab Large Scale Sci Instruments, Xuzhou 221116, Jiangsu, Peoples R China
[4] Hokkaido Univ, Fac Engn, Sapporo, Hokkaido 0608628, Japan
关键词
Crack initiation and coalescence; Hybrid FDEM; Initiation stress; Coalescence stress; Rocks; ROCK-LIKE MATERIALS; HYDRAULIC-FRACTURE PROPAGATION; NUMERICAL-SIMULATION; DYNAMIC FRACTURE; PRECRACKED ROCK; MULTIPLE FLAWS; BEHAVIOR; SPECIMENS; NETWORKS; CODE;
D O I
10.1016/j.engfracmech.2021.107658
中图分类号
O3 [力学];
学科分类号
08 ; 0801 ;
摘要
Discrete element method or extended finite element method is often adopted to simulate crack initiation, propagation, and coalescence. However, the transition of rock from continuous medium to the discontinuous medium cannot be simulated. Moreover, previous studies on propagation and coalescence of prefabricated flaws often focus on single or multiple parallel prefabricated flaws. In reality, the fractures in rock mass usually exist in the form of cross flaws. To solve the inconsistency, the uniaxial compression of rock samples with prefabricated crossflaws is modelled in this study through the GPGPU-parallelized Y-HFDEM IDE, and the modelling results are compared with those from laboratory experiments. The effects of four geometric parameters of the prefabricated cross-flaws on the crack initiation, propagation, and coalescence process are explored: (1) the joint continuity parameter k, (2) the angle ? between the axial load direction and the dip direction of the primary flaw, (3) the angle ? between the direction of the rock bridge and the dip direction of the primary flaw, and (4) the angle ? between primary and secondary flaws. It is found that numerical simulation results are in good agreement with those of laboratory experiments in terms of not only the qualitative crack initiation, propagation, and coalescence processes but also the quantitative crack initiation locations and stress, crack propagation mechanism, and coalescence locations. It is concluded that the presence of the prefabricated cross flaws promotes the rock samples fracture more seriously during uniaxial loading. Moreover, two new crack coalescence mechanisms that have not been reported in previous studies are observed through numerical simulations and laboratory experiments. Thus, compared with laboratory experiments, the GPGPU-parallelized Y-HFDEM IDE provides another cheaper but more flexible and robust tool for the study on the crack initiation, propagation, and coalescence process in rocks with prefabricated flaws.
引用
收藏
页数:24
相关论文
共 86 条
[1]  
Akara MEM, HYDROGEOL J
[2]   Hybrid finite-discrete element modelling of dynamic fracture and resultant fragment casting and muck-piling by rock blast [J].
An, H. M. ;
Liu, H. Y. ;
Han, Haoyu ;
Zheng, Xin ;
Wang, X. G. .
COMPUTERS AND GEOTECHNICS, 2017, 81 :322-345
[3]  
[Anonymous], 2015, 49 US ROCK MECH GEOM
[4]   Numerical investigation of the effect of joint geometrical parameters on the mechanical properties of a non-persistent jointed rock mass under uniaxial compression [J].
Bahaaddini, M. ;
Sharrock, G. ;
Hebblewhite, B. K. .
COMPUTERS AND GEOTECHNICS, 2013, 49 :206-225
[5]   Estimation of natural asbestos content in rocks by fracture network modeling and petrographic characterization [J].
Barale, Luca ;
Botta, Serena ;
Piana, Fabrizio ;
Tallone, Sergio ;
Fidelibus, Corrado ;
Avataneo, Chiara ;
Turci, Francesco ;
Compagnoni, Roberto ;
Cossio, Roberto ;
Alberto, Walter .
ENGINEERING GEOLOGY, 2020, 271
[6]  
Ben Y, 2013, P 11 INT C AN DISC I
[7]   The 3D Numerical Simulation for the Propagation Process of Multiple Pre-existing Flaws in Rock-Like Materials Subjected to Biaxial Compressive Loads [J].
Bi, J. ;
Zhou, X. P. ;
Qian, Q. H. .
ROCK MECHANICS AND ROCK ENGINEERING, 2016, 49 (05) :1611-1627
[8]   Fracture coalescence in rock-type materials under uniaxial and biaxial compression [J].
Bobet, A ;
Einstein, HH .
INTERNATIONAL JOURNAL OF ROCK MECHANICS AND MINING SCIENCES, 1998, 35 (07) :863-888
[9]   The initiation of secondary cracks in compression [J].
Bobet, A .
ENGINEERING FRACTURE MECHANICS, 2000, 66 (02) :187-219
[10]  
Brown E.T., 1981, Rock characterization, Testing and monitoring, ISRM suggested methods