基于Voronoi多边形离散的DDA方法模拟岩石破坏

doi:10.13722/j.cnki.jrme.2020.0745

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

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

摘要非连续变形分析(DDA)方法是为模拟可变形离散块体系统而提出的一种隐式离散元方法。考虑到真实岩石的细观晶粒结构更加接近于Voronoi多边形，将建立基于Voronoi多边形离散形式的DDA方法，并将其与基于三角形离散的DDA方法进行对比分析。首先研究不同子块体尺寸对材料弹性参数的影响，并对比2种离散方法在子块体数和接触对数量上的差异。通过模拟单轴压缩、巴西劈裂、带孔圆盘劈裂以及Kalthoff-Winkler冲击等试验，对比分析2种离散形式对岩石宏观参数、破坏形式以及裂纹扩展路径的影响。同时，为了减小网格的依赖性，考虑了3种程度的随机离散网格对结果的影响。结果表明，过于规则的网格离散形式会导致不合理的破坏形式，过于随机的网格离散形式可能导致结果离散性偏大，收敛性变差，而一定程度的随机性则是必要的；基于Voronoi多边形离散形式的DDA方法能够有效地模拟岩石的破坏过程；相比三角形离散，基于Voronoi多边形离散的DDA方法具有更高的计算效率，能得到更高的弹性模量、抗压、抗拉强度和更加真实的压拉比，且裂纹扩展路径也更接近于岩石的沿晶断裂。

	服务

	把本文推荐给朋友
	加入我的书架
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	张开雨1
	2
	夏开文1
	2
	刘丰1
	2

关键词 ：岩石力学, 岩石破坏, 非连续变形分析, Voronoi多边形, 三角形离散, 随机性, 网格依赖性

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.

Key words： rock mechanics')" href="#">

rock mechanics rock failure discontinuous deformation analysis Voronoi polygon triangle discretization randomness mesh-dependence

引用本文:

张开雨1，2，夏开文1，2，刘丰1，2. 基于Voronoi多边形离散的DDA方法模拟岩石破坏[J]. 岩石力学与工程学报, 2021, 40(4): 725-738.
ZHANG Kaiyu1，2，XIA Kaiwen1，2，LIU Feng1，2. Simulation of rock failure by Voronoi-based discontinuous deformation analysis. , 2021, 40(4): 725-738.

链接本文:

http://rockmech.whrsm.ac.cn/CN/10.13722/j.cnki.jrme.2020.0745 或 http://rockmech.whrsm.ac.cn/CN/Y2021/V40/I4/725

[10]	CUNDALL P A. UDEC—a generalised distinct element program for modelling jointed rock[R]. London：European Research Office，1980.
[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.
[29]	BAO H R. Nodal-based discontinuous deformation analysis[Ph. D. Thesis][D]. Singapore：Nanyang Technological of University，2010.
[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.
[4]	CUNDALL P A，HART R D. Numerical modelling of discontinua[J]. Engineering Computations，1992，9(2)：101-113.
[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.
[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.
[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.
[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.
[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.
[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.
[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.
[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.
[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.
[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.
[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.
[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.
[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.
[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.
[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))
[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.
[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.
[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.
[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))
[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.
[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.
[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))
[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.
[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))
[37]	ENGWIRDA D. Generalised primal-dual grids for unstructured co-volume schemes[J]. Journal of Computational Physics，2018，375： 155-176.
[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))
[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.
[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))
[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.
[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.
[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.
[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))
[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.
[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.
[28]	SHYU K. Nodal-based discontinuous deformation analysis[Ph. D. Thesis][D]. Berkeley：University of California，1993.
[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))
[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.