Numerical modeling of dynamic rock fracture with a combined 3D continuum viscodamage‐embedded discontinuity model

This paper deals with numerical modeling of dynamic failure phenomena in rate‐sensitive quasi‐brittle materials, such as rocks, with initial microcrack populations. To this end, a continuum viscodamage‐embedded discontinuity model is developed and tested in full 3D setting. The model describes the pre‐peak nonlinear and rate‐sensitive hardening response of the material behavior, representing the fracture‐process zone creation, by a rate‐dependent continuum damage model. The post‐peak response, involving the macrocrack creation accompanied by exponential softening, is formulated by using an embedded displacement discontinuity model. The finite element implementation of this model relies upon the linear tetrahedral element, which seems appropriate for explicit dynamic analyses involving stress wave propagation. The problems of crack locking and spreading typical of embedded discontinuity models are addressed in this paper. A combination of two remedies, the inclusion of viscosity in the spirit of Wang's viscoplastic consistency approach and introduction of isotropic damaging into the embedded discontinuity model, is shown to be effective in the present explicit dynamics setting. The model performance is illustrated by several numerical simulations. In particular, the dynamic Brazilian disc test and the Kalthoff–Winkler experiment show that the present model provides realistic predictions with the correct failure modes and rate‐dependent tensile strengths of rock at different loading rates. The ability of initial embedded discontinuity populations to model the initial microcrack populations in rocks is also successfully tested. Copyright © 2016 John Wiley & Sons, Ltd.     

扩展有限元法在岩石断裂力学中的应用

岩石裂纹的扩展是一个经典的不连续问题,常规有限元方法难以实现裂纹扩展过程的仿真模拟。扩展有限元法(XFEM)实现了计算网格与不连续面相互独立,因此模拟移动的不连续面时无需对网格进行重新剖分。本文介绍了XFEM基本原理和岩石断裂力学常用判据,尝试对岩石类材料单缝Ⅰ型三点弯曲、单缝剪切和双缝平板实验进行模拟。分析结果表明:扩展有限元模拟岩石类材料断裂问题不受网格划分限制,裂纹以实际应力场分布随机扩展;直观地给出岩样的微裂纹产生、演化,直至完全破坏的全过程,并与实验结果吻合。该方法能够应用到岩石断裂力学方面的研究,模拟岩石类材料的宏细观破坏过程,为解决复杂问题提供了方便的途径。     

基于扩展有限元的深部岩体分区破裂化现象初步数值模拟

通过轮廓法追踪岩石裂纹扩展轨迹,将其嵌入ABAQUS扩展有限元（EFE）平台,对单轴压力下裂隙试件受压过程进行了模拟验证,裂纹起裂扩展效果良好。以此为平台,将推导得到的弧形裂纹应力强度因子嵌入其中,以最大周向拉应力准则为开裂准则,认为当其大于岩体断裂韧度时,硐室围岩体内初始裂纹将开始扩展。在此基础上,以发现分区破裂现象的圆形隧道模型试验为背景开展了分区破裂的数值模拟试验。模拟结果发现,深部巷道围岩出现了3~4层破裂分区,证实深部巷道围岩存在分区破裂现象。将数值模拟结果与模型试验完成后模型围岩破坏状态对比,发现二者破裂区分布特征基本一致。数值模拟结果表明,EFEM方法在处理复杂岩体裂纹问题方面的有效性。     

岩石裂纹损伤及围压对剪切破坏影响研究

岩石内部细观裂纹的存在,对压缩作用下岩石剪切断裂的宏观现象有着重要的影响。然而,能够通过解析解阐释细观裂纹几何特性、围压等影响因素对压缩作用下剪切断裂面角度变化趋势的研究很少。基于Ashby模型中提出的裂纹尖端应力强度因子,提出了一种改进的考虑裂纹角度影响的应力强度因子表达式。利用该改进的应力强度因子表达式,推出了一个可以预测岩石峰值强度的裂纹扩展、应变与应力之间的本构关系。结合本构关系的峰值强度与摩尔-库仑失效准则,得到了岩石损伤与内摩擦角、黏聚力、剪切强度及失效断裂面角度之间的理论关系;讨论了围压、裂纹尺寸、角度及摩擦系数对岩石宏观剪切断裂面角度的影响,通过试验结果验证了模型合理性。结果表明：随着损伤增大,内摩擦角、黏聚力及剪切强度不断减小;随着围压增大、摩擦系数增大和初始裂纹尺寸减小,剪切断裂面角度不断增大;随着裂纹角度增大,剪切断裂纹面角度先减小后增大。     

Failure pattern of brittle rock governed by initial microcrack characteristics

Fracturing processes of rock are simulated using a micro- and fracture mechanical-based numerical modeling approach. The numerical model considers material heterogeneity and initial microcrack distributions. The gradual formation of macroscopic fractures by coalesced microcracks is successfully reproduced. Distinct failure modes are observed in the model under different loadings. Agreement is shown between the numerical results and laboratory observations. The influence of microcrack orientations on the fracture patterns is quantified by the numerical models. Possible mechanisms describing the splitting failure of hard brittle rocks under uniaxial compression are proposed and discussed.

     

Numerical modeling of rock failure under dynamic loading with polygonal elements

Polygonal finite elements are gaining an increasing attention in the computational mechanics literature, but their application in rock mechanics is very rare. This paper deals with numerical modeling of rock failure under dynamic loading based on polygonal finite elements. For this end, a damage-viscoplastic constitutive model for rock based on the Mohr-Coulomb criterion with the Rankine criterion as a tensile cutoff is employed and implemented with the polygonal finite element method. Moreover, the mineral mesostructure or rock is described by randomly mapping groups of polygonal elements representing the constituent minerals into a global mesh and assigning these groups with the corresponding mineral material properties. The performance of the polygonal elements is compared with that of the linear and quadratic triangular and bilinear quadrilateral elements in numerical simulations of controlled shear band formation under uniaxial compression and lateral splitting failure in the dogbone tension test. Numerical simulations of uniaxial tension and compression tests as well as dynamic Brazilian disc test under increasing loading rates demonstrate that the present approach predicts the correct failure modes as well as the dynamic increase in strength of rock.     

Life-time prediction for rocks under static compressive and tensile loads: a new simulation approach

The life time or time to failure of rocks under load is governed by microstructural defects, like microcracks, voids etc. The life time can be predicted either by empirical exponential laws or physical laws based on damage and fracture mechanics. The proposed numerical model is based on subcritical crack growth using the linear elastic fracture mechanical approach and is implemented as a numerical cellular automate. The algorithm considers both tensile and shear fracturing. Each cell contains a microcrack of random length according to a given probability function. Fracture growth is controlled by the Charles equation. Macroscopic cracks are the results of the coalescence of growing microcracks. Within the numerical approach elasto-plastic stress redistributions take place. If the stress intensity factors have reached the critical values or the microcrack has reached the zone dimension, the zone is considered as fractured and residual strength values are assigned. The proposed approach was applied to rock samples under uniaxial compressive and tensile loads (creep tests). Successful results were obtained in respect to the predicted life time, damage evolution and the fracture pattern. Conclusions for further improvements and extensions of this methodology were drawn.     

模拟岩石中裂纹扩展连接的近场动力学方法

脆性岩石材料在压应力作用下常出现两类裂纹：翼型张拉裂纹和次生剪切裂纹。近场动力学是一种新型的无网格数值计算方法。在近场动力学理论中,采用积分形式的控制方程代替微分形式的控制方程使得该数值算法在断裂问题上具有独特的优势。将Mohr-Coulomb准则和最大主应力准则引入非普通“态”基近场动力学理论中,分别用于模拟材料常见的压剪和张拉破坏。这种扩展的非普通“态”基近场动力学可以有效地模拟脆性岩石材料在多种受力状态下的裂纹起裂、扩展和连接问题。通过5个不同的数值算例说明该数值算法在处理脆性岩石材料断裂问题的有效性和准确性。首先,通过模拟含圆孔的弹性板拉伸数值试验说明该数值算法的有效性和准确性。其次,数值模拟了简单三点弯曲试验以及不使用其他外部准则条件下动荷载作用下裂纹的分叉试验,所得结果与其他试验结果或数值结果相吻合,从而验证了该理论的有效性。然后,模拟了包含斜裂纹的巴西圆盘试验,裂纹扩展路径和计算所得的断裂韧度同样吻合于试验结果。最后,模拟了单轴压缩状态下,预制裂纹试样的裂纹扩展和连接问题。将该数值算法与试验结果对比表明,所提出的数值方法可以模拟和预测岩石类材料的张拉和压剪裂纹的起裂、扩展和连接行为。     

三点弯曲条件下薄层状岩体单层厚度对裂纹扩展路径的影响

为了揭示薄层状岩体单层厚度对裂纹扩展路径的影响,利用ABAQUS软件进行数值模拟,对含预制裂纹的不同单层厚度的层状岩体进行三点弯曲数值试验,岩体试样长度和宽度一定（512 mm×128 mm）,分5种不同的单层厚度建模,单层厚度分别为:64 mm、32 mm、16 mm、8 mm、4 mm。提取破裂过程中的声发射信息及主裂纹路径参数,对比分析不同层理厚度时,岩体的峰值载荷、声发射事件数、累计声发射事件数、主裂纹长度等对单层厚度的变化规律。研究结果表明:（1）薄层状岩体的强度与单层厚度有一定关系,单层厚度越大,岩体的抗弯强度越小,并且这种变化是非线性的。（2）在不同层厚条件下,层理面均表现出了阻止裂纹向上扩展的效果。同时,裂纹在层状岩体中的扩展是不连续的,裂尖沿层理面发生了水平迁移,整体呈现阶梯状。（3）在三点弯曲条件下,随着单层厚度的减小,岩体破裂时的裂纹总长度先增大后减小。（4）裂纹在岩体层面扩展时破裂类型多为滑移与拉张破裂的组合,而在非层理面的岩体内部扩展时,以拉张破裂为主。（5）薄层状岩体在三点弯曲条件下,声发射事件在岩体的压密阶段几乎没有出现,直至接近峰值载荷时才大量涌现。到了载荷快速下降时,单位时间的声发射事件数也达到最大。研究结果将对预测层状岩石的裂缝演化提供依据。     

Three-dimensional elastic analysis of cracks with the g2 constant displacement discontinuity method

A new triangular element was created that could be used for the improvement of the accuracy of the constant displacement discontinuity method (CDDM). This element is characterized by three degrees of freedom in the three-dimensional space as in the classical CDDM approach. The element is based on strain gradient elasticity theory that accounts for the difference of the average value of stress with the local stress at surfaces with large curvature (eg, crack borders, corners, and notches) in elastic bodies. The new element is characterized by a strain gradient term in addition to the two Lamè constants that gives a more representative value of the stresses at the centroid of crack edge elements compared with the classical elasticity solution and thus an accurate stress intensity factor. In this approach, special crack border elements with square-root radius dependent displacements and numerical integrations are avoided. The extra strain gradient term is calibrated once only on the analytical solution for the penny-shaped crack. In a verification stage, the accuracy of the computational algorithm for the elliptic and rectangular crack problems is demonstrated. Then, the algorithm, which also accounts for crack closure in compression, is applied for the modeling of crack propagation and crack interaction in uniaxial tension and compression loading. It is illustrated that the numerical predictions are in accordance with experimental evidence pertaining to uniaxial compression of transparent precracked specimens in the lab.     

基于分形与损伤理论的岩石声-应力相关性理论模型研究

建立单一微裂纹在单轴压缩荷载作用下的力学模型,采用宏观复合型断裂理论中的最大周向应力理论,求出单一微裂纹达到稳定时的形态尺寸.利用分形损伤理论,求出不同荷载阶段的分形维数;由分形维数的定义导出了不同荷载阶段的裂纹个数的变化情况;结合单一微裂纹的形态尺寸,求出不同荷载阶段中岩石孔隙率与荷载值之间的关系;利用目前被广泛采用的超声波纵波速度与孔隙率的关系式,最终导出单轴压缩荷载作用下岩石超声波纵波速度与应力的理论关系式.进行岩石声-应力的试验测试,根据试验成果作出声-应力曲线,并与理论声-应力曲线进行比对,结果表明,理论模型比较适合于反映岩石受载全过程中的后半部分的声-应力关系.     

Analyses of Inclined Cracks Neighboring Two Iso-Path Cracks in Rock-Like Specimens Under Compression

Crack propagation process in pre-cracked rock like specimens has been studied experimentally and numerically considering three cracks in the middle part of each specimen. The rock-like specimens are specially prepared from Portland pozzolana cement, fine sands and water. These pre-cracked cylindrical specimens (each containing a single inclined crack in the neighborhood of two iso-path cracks) are experimentally tested under compressive loading. The same problems are numerically simulated by a modified displacement discontinuity method using higher order displacement discontinuity elements and higher order special crack tip elements for crack tip treatment to increase the accuracy of the Mode I and Mode II stress intensity factors obtained based on linear elastic fracture mechanics theory. The crack propagation and coalescence paths of the inclined crack are estimated by implementing a suitable iteration algorithm of incremental crack length extension in a direction predicted by using the maximum tangential stress criterion. The numerical and analytical crack extension analyses are compared which are in good agreement and show the validity, applicability and accuracy of the present work.     

单裂隙砂岩单轴压缩的中等应变率效应颗粒流模拟

加载速率对裂隙岩体的力学性质及变形破坏均有重要影响。利用二维颗粒流程序PFC2D开展了不同倾角非贯通单裂隙砂岩试件的单轴压缩试验,研究了中等应变率对裂隙砂岩应力-应变曲线特征、裂隙尖端应力状态、特征应力状态、岩体损伤及裂隙扩展等力学响应的影响规律。裂隙岩体应力-应变曲线呈现明显的波动性,定义应力突变指标 对应力突变型波动剧烈程度进行了定量统计分析：随应变率的增加,曲线应力突变波动越剧烈,且峰后明显大于峰前;随裂隙倾角的增大,波动幅度峰前增大,而峰后减小。裂隙尖端破裂应力随应变率增大均有所提高,随裂隙倾角的增大,切向剪应力 总体上呈增加变化,而法向应力 明显减小。尖端破裂时岩样加载应力 、岩样临界扩容应力 及峰值应力 均随应变率增大而增大。裂隙尖端的破裂可立即引起岩体扩容,一般应变率越低,岩体裂隙尖端破裂点 和扩容点 越接近峰值强度 。随着应变率的提高,损伤裂纹及宏观裂隙类型越多,岩体试件损伤破裂程度越强,特别是试件端部效应愈显著。裂隙首先以I型翼裂纹在其尖端起裂,而I型翼裂纹的扩展长度与加载速率与裂隙倾角具有较强的相关性。     

On the role of pre-existing discontinuities on the micromechanical behavior of confined rock samples: a numerical study

The deformation process and failure mechanism of rock mass with increased density of initial joints subjected to confined stress state are investigated in this study using discrete element method (DEM). A numerical model of standard size granite samples is developed and validated using experimental data for both intact and jointed rocks. The micro-parameters of the rock material are first determined, and the effects of the rock discontinuity on strength, deformability, stress–strain relationship, and failure modes are then investigated at the macro-scale level. Analyses are also performed to examine the tensile and shear crack distributions, fragmentation characteristics, particle kinematics, and energy dissipation to advance the current understanding of the deformation processes and failure mechanisms of jointed rock masses. The microscopic evolutions in the fabric and force anisotropy during loading and distributions of contact forces provide insights into the influence of increasing initial jointing on the macroscopic deformational behavior of the rock. The results show how the deceleration in the growth of fabric and contact force anisotropies develops and confirms that the increase in initial jointing and the associated changes in microstructure can restrain the development of anisotropy, thereby reducing significantly the strength of the rock samples.

     

模拟岩石破裂过程的块体单元离散弹簧模型

在变形体离散元的基础上建立块体单元离散弹簧模型,并应用于岩石破裂过程的数值模拟研究。该模型以连续介质力学理论为基础,将块体单元离散为具有明确物理意义的弹簧系统,通过对弹簧系统的能量泛函求变分获得各弹簧的刚度系数,进而可以直接利用弹簧刚度求解单元的变形和应力,提高计算效率。以重力作用下的岩质边坡计算为例,通过与传统的有限元进行对比,验证该模型弹性计算结果的正确性。在该基础上,引入Mohr-Coulomb与最大拉应力的复合破坏准则,判断单元的破坏状态及破裂方向。当单元的内部破坏面确定后,则通过块体切割的方式实现单元破坏,并建立单元边界和单元内部的双重破裂机制,实现块体由连续到非连续的破裂过程,进而显示的模拟裂纹的形成和扩展。最后,以巴西圆盘劈裂、单轴压缩破裂以及三点弯曲梁等典型算例验证该方法,结果表明该方法可以较好地模拟拉伸、压剪等应力状态下裂纹的形成和扩展,从而可模拟岩石介质由连续到非连续的破裂过程。     

自由边界形状与近边界裂纹相互作用模型研究

大量的研究表明, 岩爆发生条件与岩石的力学性能有关。岩体中微裂纹的形成、扩展和汇合对岩体的力学性能产生显著影响,可以导致材料的逐渐劣化直至破坏。对于单轴压缩的岩石试件,其破坏型式主要是微裂纹的摩擦滑移、自相似平面扩展和弯折扩展。基于断裂力学机制,对单轴压缩荷载下3种不同边界形状,即直边、凸形、凹形情况下近边界原生斜裂纹的扩展进行了理论研究,并运用FLAC数值模拟技术,建立了基于Interface界面单元的裂纹扩展数值模型。近直边界及凸形边界情况的裂纹扩展使用梁渐近模型来解析薄片岩层的屈曲断裂问题,对于近凹形边界情况给予简化,提出了裂纹－圆－椭圆相互作用的修正模型并进行理论分析。研究结果表明,基于Interface界面单元的裂纹扩展的数值模型正确,解析解与数值解吻合较好,相应的梁渐近模型及裂纹-圆-椭圆模型可以合理地分析应力调整过程中裂纹的扩展和裂纹与不同洞室表面形状的相互作用问题,为岩爆工程及理论研究提出了新的理论依据。     

基于扩展有限元的压剪复合断裂数值分析

针对岩石类材料压剪破坏,国内外学者做了大量真实岩石和模型材料的试验,并进行相应的数值模拟工作。但由于试验结果往往具有很强的离散性,且对影响试验结果的各因素缺乏系统地理论分析,加之传统常规有限元在处理裂纹这类强不连续问题时存在很大的困难,如在裂尖附近的高应力区需要剖分令人难以接受的密度的网格,同时在模拟裂纹生长时还需要对网格进行重新剖分,操作繁琐复杂,效率极低,因此,针对压剪复合断裂翼裂纹扩展规律的研究始终没有取得令人满意的成果。对原生裂纹的几何形态,主裂纹面与其间填充物的相互作用对翼裂纹起裂角的影响做了系统的理论分析,并利用扩展有限元对翼型裂纹扩展全过程进行数值模拟,与理论分析结果对比,吻合良好。相关成果为以后针对压剪复合断裂的研究提供了理论分析框架,同时对今后相关试验的开展具有理论指导作用     

单轴拉伸岩样破坏过程及尺寸效应数值模拟

由于从实验及理论角度研究岩样单轴拉伸条件下的破坏全过程及尺寸效应难度都很大。因此采用拉格朗日元法来研究这些问题。在峰值强度之前后,岩石材料的本构模型分别取为线弹性及拉破坏线性应变软化模型。为了使拉伸塑性区不出现在试样的端部,在试样的两侧面中部预制了2个凹槽。数值模拟结果表明,全程拉应力-拉应变曲线分为峰前和峰后阶段。在接近峰值的峰前阶段,由于两凹槽附近具有明显的拉应力集中现象,拉伸塑性区最先出现在两凹槽附近。随着轴向拉应变的增加,发生拉伸破坏的单元的数目增加,新发生拉伸破坏的单元越来越接近试样的中心,直到两块拉伸塑性区在应变软化阶段贯通。两凹槽连线上各单元拉应力的分布呈现3个阶段,“澡盆型”（“U型”）阶段,“双峰型”（“M型”）阶段及“单峰型”（“Π型”）阶段。“澡盆型”阶段对应于全程拉应力-拉应变曲线的弹性阶段。“双峰型”阶段及“单峰型”阶段对应于全程拉应力-拉应变曲线的非弹性阶段（包括峰值强度之前的一小段,即应变硬化阶段及峰后的应变软化阶段）。增加试样高度及降低试样宽度,拉应力-拉应变曲线的软化段变得越来越陡峭,因而试样越容易发生失稳破坏。由于试样宽度较大时,试样内部的单元并非处于单向拉应力状态,因此,增加试样宽度,全程拉应力-拉应变曲线的峰值强度增加。当试样宽度较小时,从出现塑性区,到塑性区贯通所需要的时间步较小,或应变范围较窄。这说明试样的脆性较强,前兆不明显。前兆不明显的脆性破坏对应常见的是洞室岩爆、冲击地压及地震等灾害。     

一种改进的光滑粒子动力学法在地震条件下岩质边坡中的验证与应用

开发一种改进损伤框架的粒子流算法,被称为核断裂的光滑粒子流法（kernel-broken smoothed particle hydrodynamics,KBSPH）,用于模拟地震条件下岩质边坡的裂纹扩展和变形破坏过程.在KBSPH中,提出一种改进的损伤框架,通过引入断裂标志来改进损伤粒子的核函数,使损伤粒子的虚拟应力键直接断裂,裂纹在断裂的应力键间生成,从而模拟岩石的裂纹扩展过程.在地震边界上采用了双层边界,将动力输入边界与黏滞边界分离.首先通过薄板振动实验验证KBSPH的动力特性.其次以单裂隙岩体单轴压缩试验验证KBSPH的断裂力学特性.最后模拟地震条件下多节理岩质边坡中裂纹扩展过程和动力响应.薄板振动实验验证了KBSPH的动力特性的准确性.单裂隙岩体单轴压缩试验,证明了KBSPH可以正确模拟预制裂隙尖端的翼型裂纹.通过对比以往数值模拟方法和现场案例,表明KBSPH正确揭示了加速度放大效应以及地震条件下岩质边坡的裂纹扩展过程.KBSPH避免了传统算法的网格畸变,损伤粒子应力分量重新分配的问题,降低了编程难度,提高了运行速率,可为SPH在地震条件下岩石力学中的应用和理解岩石断裂机理提...     

拉-压应力状态下脆性岩石强度及破坏机制颗粒流模拟

高地应力区地下岩体工程开挖常形成围岩拉-压应力状态,发生岩体张性破坏灾害。本文针对传统PFC平行黏结模型不能模拟脆性岩石高单轴压缩与拉伸强度比的问题,建立双抗拉强度参数的平行黏结强度准则,开展岩石拉-压数值模拟试验,得到了与物理试验接近的拉-压强度,实现了岩石高压拉强度比的模拟,并深入分析了破坏机制。研究结果表明随着围压的增加,破裂面倾角逐渐增大,由拉伸破裂转化为拉-剪破裂,发现了拉-压应力状态下破裂面处的雁行裂纹。根据细观颗粒位移场揭示了破裂面力学性质,随着围压的增加(破裂面倾角逐渐增大),破裂面张性逐渐减弱而剪性增强。可将拉-压应力状态下岩石损伤演化过程大致分为弹性变形阶段、稳定破裂发展阶段、不稳定破裂发展阶段和整体破裂阶段(峰后应力跌落及残余阶段)。围压较大时弹性变形和稳定破裂发展阶段相对较短,不稳定破裂发展阶段相对较长较剧烈,峰后残余阶段破裂面摩擦更强、应力波动较大。