|
|
|
| Simulation of rock failure by Voronoi-based discontinuous deformation analysis |
| ZHANG Kaiyu1，2，XIA Kaiwen1，2，LIU Feng1，2 |
| (1. State Key Laboratory of Hydraulic Engineering Simulation and Safety，Tianjin University，Tianjin 300350，China；2. Department of Architectural Engineering，Tianjin University，Tianjin 300350，China) |
|
|
|
|
Abstract Discontinuous deformation analysis(DDA) is an implicit discrete element method for simulating deformable discrete block systems. Since the meso-grain structure of rocks is similar to Voronoi polygons，this paper will study the rock failure simulation with DDA using Voronoi discretizations，and compare its performance with DDA using triangle discretizations. The influence of the sub-block size on the elastic parameters of the material is studied，and the numbers of sub-blocks and the contact pairs using these two discretizations are compared. By simulating the uniaxial compression test，Brazilian test，Brazilian test with an eccentric circular hole and Kalthoff- Winkler impact experiment，the influences of these two types of discretizations on rock macroscopic parameters，failure modes and crack propagation paths are analyzed. Meanwhile，in order to reduce mesh-dependence，the influence of the randomness of the discretizations on the simulation results is discussed. The results show that an overly regular discretization will lead to an unreasonable failure pattern，while an overly random discretisation may lead to more disperse results and also worse convergence. That is to say，a proper degree of randomness is necessary for such simulations. It is also shown that DDA with Voronoi discretizations can effectively simulate the failure process of rocks，and that, in comparison with the triangle model，the Voronoi model leads to greater computational efficiency，higher elastic modulus，higher compression strength，higher tensile strength and more realistic compression-tensile ratio，and the crack propagation path is similar to the inter-granular fracture of the rock.
|
|
|
|
|
|
[1] TANG C A. Numerical simulation of progressive rock failure and associated seismicity[J]. International Journal of Rock Mechanics Mining Sciences，1997，34(2)：249–261.
[2] ZHU W C，TANG C A. Numerical simulation of Brazilian disk rock failure under static and dynamic loading[J]. International Journal of Rock Mechanics Mining Sciences，2006，43(2)：236–252.
[3] 马天辉，唐春安，杨天鸿，等. 基于非均匀网格的岩石破坏过程分析系统[J]. 东北大学学报：自然科学版，2007，(8)：1 182–1 185. (MA Tianhui，TANG Chun¢an，YANG Tianhong，et al. Rock failure process analysis based on non-uniform grids[J]. Journal of Northeastern University：Natural Science，2007，(8)：1 182–1 185.(in Chinese))
[4] CUNDALL P A，HART R D. Numerical modelling of discontinua[J]. Engineering Computations，1992，9(2)：101–113.
[5] POTYONDY D O，CUNDALL P A. A bonded-particle model for rock[J]. International Journal of Rock Mechanics and Mining Sciences，2004，41(8)：1 329–1 364.
[6] CHO N，MARTIN C D，SEGO D C. A clumped particle model for rock[J]. International Journal of Rock Mechanics and Mining Sciences.，2007，44(7)：997–1 010.
[7] POTYONDY D O. A flat-jointed bonded-particle material for hard rock[C]// Proceedings of the 46 textsuperscriptth U.S. Rock Mechanics/ Geomechanics Symposium. [S. l.]：[s. n.]，2012：24–27.
[8] 蒋明镜，陈 贺，刘 芳. 岩石微观胶结模型及离散元数值仿真方法初探[J]. 岩石力学与工程学报，2013，32(1)：15–23.(JIANG Mingjing，CHEN He，LIU Fang. A microscopic bond model for rock and preliminary study of numerical simulation method by distinct element method[J]. Chinese Journal of Rock Mechanics and Engineering，2013，32(1)：15–23.(in Chinese))
[9] 金 磊，曾亚武，李 欢，等. 基于不规则颗粒离散元的土石混合体大三轴数值模拟[J]. 岩土工程学报，2015，37(5)：829–838.(JIN Lei，ZEN Yawu，LI Huan，et al. Numerical simulation of large-scale triaxial tests on soil-rock mixture based on DEM of irregularly shaped particles[J]. Chinese Journal of Geotechnical Engneering，2015，37(5)：829–838.(in Chinese))
[10] CUNDALL P A. UDEC—a generalised distinct element program for modelling jointed rock[R]. London：European Research Office，1980.
[11] KAZERANI T，ZHAO J. Micromechanical parameters in bonded particle method for modelling of brittle material failure[J]. International Journal for Numerical and Analytical Methods in Geomechanics，2010，34(18)：1 877–1 895.
[12] KAZERANI T. Effect of micromechanical parameters of microstructure on compressive and tensile failure process of rock[J]. International Journal of Rock Mechanics and Mining Sciences，2013，64：44–55.
[13] GAO F Q，STEAD D. The application of a modified Voronoi logic to brittle fracture modelling at the laboratory and field scale[J]. International Journal of Rock Mechanics and Mining Sciences，2014，68：1–14.
[14] ZHAO G F，TOHID K，KE M，et al. Numerical study of the semi-circular bend dynamic fracture toughness test using discrete element models[J]. Science China Technological Sciences，2015，58(9)：1 587–1 595.
[15] GAO F，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.
[16] CHRISTIANSON M，BOARD M，RIGBY D. UDEC simulation of triaxial testing of lithophysal tuffs[C]// Proceedings of the 41st U.S. Rock Mechanics Symposium. [S. l.]：[s. n.]，2006：28–43.
[17] MAYER J M，STEAD D. Exploration into the causes of uncertainty in UDEC Grain Boundary Models[J]. Computers and Geotechnics，2017，82：110–123.
[18] LAN H，MARTIN C D，HU B. Effect of heterogeneity of brittle rock on micromechanical extensile behavior during compression loading[J]. Journal of Geophysical Research，2010，115(B1)：1–14.
[19] 严成增，郑 宏，孙冠华，等. 粗粒料多边形表征及二维FEM/DEM分析[J]. 岩土力学，2015，36(增2)：95–103.(YAN Chengzeng，ZHENG Hong，SUN Guanhua，et al. Polygon characterization of coarse aggregate and two-dimensional combined finite discrete element method analysis[J]. Rock and Soil Mechanics，2015，36(Supp.2)：95–103.(in Chinese))
[20] 周 伟，常晓林，周创兵，等. 堆石体应力变形细观模拟的随机散粒体不连续变形模型及其应用[J]. 岩石力学与工程学报，2009，28(3)：491–499.(ZHOU Wei，CHANG Xiaolin，ZHOU Chuangbin，et al. Stochastic granule discontinuous deformation model of rockfill and its application[J]. Chinese Journal of Rock Mechanics and Engineering，2009，28(3)：491–499.(in Chinese))
[21] 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.
[22] SHI G H. Manifold method of material analysis[C]// Transactions of the 9th Army Conference on Applied Mathematics and Computing. [S. l.]：[s. n.]，1992：57–76.
[23] WU Z，FAN L，LIU Q，et al. Micro-mechanical modeling of the macro-mechanical response and fracture behavior of rock using the numerical manifold method[J]. Engineering Geology，2017，225：49–60.
[24] LIU Q，JIANG Y，WU Z，et al. A Voronoi element based-numerical manifold method(VE-NMM) for investigating micro/macro- mechanical properties of intact rocks[J]. Engineering Fracture Mechanics，2018，199：71–85.
[25] 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.
[26] KE T C. Application of DDA to simulate fracture propagation in solid[C]// YUZO Ohnishi. Proceedings of the Second International Conference on Analysis of Discontinuous Deformation(ICADD–2). [S. l.]：[s. n.]，1997：155–185.
[27] LIN C T，AMADEI B，JUNG J，etc. 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.
[28] SHYU K. Nodal-based discontinuous deformation analysis[Ph. D. Thesis][D]. Berkeley：University of California，1993.
[29] BAO H R. Nodal-based discontinuous deformation analysis[Ph. D. Thesis][D]. Singapore：Nanyang Technological of University，2010.
[30] JIAO Y Y，HUANG G，ZHAO Z，et al. An improved three-dimensional spherical DDA model for simulating rock failure[J]. Science China Technological Sciences，2015，58：1 533–1 541.
[31] 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.
[32] 焦玉勇，张秀丽，刘泉声，等. 用非连续变形分析方法模拟岩石裂纹扩展[J]. 岩石力学与工程学报，2007，26(4)：682–691.(JIAO Yuyong，ZHANG Xiuli，LIU Quanshen，et al. Simulation of rock crack propagation using discontinuous deformation analysis method[J]. Chinese Journal of Rock Mechanics and Engineering，2007，26(4)： 682–691.(in Chinese))
[33] NING K，YANG J，NING Y，et al. A modified sub-block DDA fracturing modelling method for rock[J]. Engineering Analysis with Boundary Elements，2020，111：154–166.
[34] 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–49.
[35] 徐栋栋，邬爱清，李 聪，等. 破裂全过程模拟的改进非连续变形分析方法[J]. 岩土力学，2019，40(3)：1 169–1 178.(XU Dongdong，WU Aiqin，LI Cong，et al. An improved discontinuous deformation analysis method for simulation of whole fracturing process[J]. Rock and Soil Mechanics，2019，40(3)：1 169–1 178.(in Chinese))
[36] PEARCE C J，THAVALINGAM A，LIAO Z，et al. Computational aspects of the discontinuous deformation analysis framework for modelling concrete fracture[J]. Engineering Fracture Mechanics，2000，65(2)：283–298.
[37] ENGWIRDA D. Generalised primal-dual grids for unstructured co-volume schemes[J]. Journal of Computational Physics，2018，375： 155–176.
[38] 周 旺，张晨麟，吴建鑫. 一种基于Hartigan-Wong和Lloyd的定性平衡聚类算法[J]. 山东大学学报，2016，46(5)：37–44.(ZHOU Wang，ZHANG Chenlin，WU Jianxin. Qualitative balanced clustering algorithm based on Hartigan-Wong and Lloyd[J]. Journal of Shandong University：Engineering Science，2016，46(5)：37–44.(in Chinese))
[39] AMADEI B，LIN C T，STURE S，etc. Modeling Fracturing of Rock Masses With the DDA Method[C]// Proceedings of the 1st North American Rock Mechanics Symposium. [S. l.]：[s. n.]，1994.
[40] HOEK E，BROWN E T. Practical estimates of rock mass strength[J]. International Journal of Rock Mechanics and Mining Sciences，1997，34(8)：1 165–1 186.
[41] BAI Q S，TU S H，ZHANG C. DEM investigation of the fracture mechanism of rock disc containing hole(s) and its influence on tensile strength[J]. Theoretical and Applied Fracture Mechanics，2016，86：197–216.
[42] KALTHO F，WINKLER S. Failure mode transition at high rates of loading[C]// International conference on impact loading and dynamic behavior of materials. [S. l.]：[s. n.]，1987：185–195.
[43] SAKSALA T，BRANCHERIE D，IBRAHIMBEGOVIC A. Numerical modeling of dynamic rock fracture with a combined 3D continuum viscodamage-embedded discontinuity model[J]. International Journal for Numerical and Analytical Methods in Geomechanics，2016，40(9)： 1 339–1 357. |
|
|
|
2021, Vol. 40 
