Adaptive phase field simulation of quasi-static crack propagation in rocks

Abstract In this study, we present an adaptive phase field method (APFM) for modeling quasi-static crack propagation in rocks. Crack initiation due to positive strains is considered, and a numerical simulation is implemented using a commercial software, COMSOL Multiphysics. Two benchmark tests are first examined, namely, a single-edge-notched square plate subjected to respective tension and shear loadings. The crack propagation in Brazil splitting tests, 2D notched semi-circular bend (NSCB) tests, and 3D NSCB tests are subsequently simulated and analyzed. All the numerical examples indicate that the propagation of the cracks is autonomous and external fracture criteria are not required for phase field modeling. Furthermore, the adaptive remeshing scheme reduces unnecessary global mesh refinement and exhibits good adaptability for fracture modeling. The simulations are in good agreement with the experimental observations, and thereby indicate the feasibility and practicability of the APFM in rocks (even in 3D cases).

[1]  Timon Rabczuk,et al.  A software framework for probabilistic sensitivity analysis for computationally expensive models , 2016, Adv. Eng. Softw..

[2]  B. Bourdin,et al.  The Variational Approach to Fracture , 2008 .

[3]  Ted Belytschko,et al.  Cracking particles: a simplified meshfree method for arbitrary evolving cracks , 2004 .

[4]  G. R. Johnson,et al.  Eroding interface and improved tetrahedral element algorithms for high-velocity impact computations in three dimensions , 1987 .

[5]  Timon Rabczuk,et al.  Fracture properties prediction of clay/epoxy nanocomposites with interphase zones using a phase field model , 2017 .

[6]  Khader M. Hamdia,et al.  Hybrid nonlinear surrogate models for fracture behavior of polymeric nanocomposites , 2017 .

[7]  Antonio Bobet,et al.  Coalescence of multiple flaws in a rock-model material in uniaxial compression , 2002 .

[8]  Christian Hesch,et al.  Thermodynamically consistent algorithms for a finite‐deformation phase‐field approach to fracture , 2014 .

[9]  H. Nguyen-Xuan,et al.  A geometrically non-linear three-dimensional cohesive crack method for reinforced concrete structures , 2008 .

[10]  Christian Miehe,et al.  A phase field model for rate-independent crack propagation: Robust algorithmic implementation based on operator splits , 2010 .

[11]  T. Anderson,et al.  Fracture mechanics - Fundamentals and applications , 2017 .

[12]  Antonio Bobet,et al.  Crack initiation, propagation and coalescence from frictional flaws in uniaxial compression , 2010 .

[13]  Timon Rabczuk,et al.  Phase-field analysis of finite-strain plates and shells including element subdivision , 2016 .

[14]  Timon Rabczuk,et al.  Finite strain fracture of plates and shells with configurational forces and edge rotations , 2013 .

[15]  Timon Rabczuk,et al.  Element-wise fracture algorithm based on rotation of edges , 2013 .

[16]  Timon Rabczuk,et al.  An h-adaptive thermo-mechanical phase field model for fracture , 2018 .

[17]  T. Rabczuk,et al.  Discontinuous modelling of shear bands using adaptive meshfree methods , 2008 .

[18]  Timon Rabczuk,et al.  Dual‐horizon peridynamics , 2015, 1506.05146.

[19]  Timon Rabczuk,et al.  A peridynamics formulation for quasi-static fracture and contact in rock , 2017 .

[20]  Ted Belytschko,et al.  A three-dimensional impact-penetration algorithm with erosion , 1987 .

[21]  Christian Miehe,et al.  Thermodynamically consistent phase‐field models of fracture: Variational principles and multi‐field FE implementations , 2010 .

[22]  L. Wong,et al.  Crack Coalescence in Molded Gypsum and Carrara Marble: Part 1. Macroscopic Observations and Interpretation , 2009 .

[23]  Antonio Bobet,et al.  Crack coalescence in specimens with open and closed flaws: A comparison , 2009 .

[24]  Timon Rabczuk,et al.  Detection of multiple flaws in piezoelectric structures using XFEM and level sets , 2014 .

[25]  T. Rabczuk,et al.  A three-dimensional meshfree method for continuous multiple-crack initiation, propagation and junction in statics and dynamics , 2007 .

[26]  Timon Rabczuk,et al.  Damage and fracture algorithm using the screened Poisson equation and local remeshing , 2016 .

[27]  Timon Rabczuk,et al.  A new crack tip element for the phantom‐node method with arbitrary cohesive cracks , 2008 .

[28]  Hao Cheng,et al.  An Experimental Study of Crack Coalescence Behaviour in Rock-Like Materials Containing Multiple Flaws Under Uniaxial Compression , 2014, Rock Mechanics and Rock Engineering.

[29]  T. Rabczuk,et al.  XLME interpolants, a seamless bridge between XFEM and enriched meshless methods , 2014 .

[30]  Seokwon Jeon,et al.  An experimental and numerical study of fracture coalescence in pre-cracked specimens under uniaxial compression , 2011 .

[31]  David J. Williams,et al.  Investigating the Effect of Cyclic Loading on the Indirect Tensile Strength of Rocks , 2011, Rock Mechanics and Rock Engineering.

[32]  Yuri Bazilevs,et al.  A higher-order stress-based gradient-enhanced damage model based on isogeometric analysis , 2016 .

[33]  Cv Clemens Verhoosel,et al.  A phase-field description of dynamic brittle fracture , 2012 .

[34]  Hung Nguyen-Xuan,et al.  A Phantom-Node Method with Edge-Based Strain Smoothing for Linear Elastic Fracture Mechanics , 2013, J. Appl. Math..

[35]  T. Belytschko,et al.  Extended finite element method for cohesive crack growth , 2002 .

[36]  T. Belytschko,et al.  A simplified mesh‐free method for shear bands with cohesive surfaces , 2007 .

[37]  Yun-Teng Wang,et al.  Numerical simulation of crack propagation and coalescence in pre-cracked rock-like Brazilian disks using the non-ordinary state-based peridynamics , 2016 .

[38]  C. Atkinson,et al.  Combined mode fracture via the cracked Brazilian disk test , 1982, International Journal of Fracture.

[39]  Jean-François Dubé,et al.  Non‐local damage model with evolving internal length , 2004 .

[40]  Timon Rabczuk,et al.  Dual-horizon peridynamics: A stable solution to varying horizons , 2017, 1703.05910.

[41]  Martin Entacher,et al.  Rock Failure and Crack Propagation Beneath Disc Cutters , 2015, Rock Mechanics and Rock Engineering.

[42]  Timon Rabczuk,et al.  Finite strain fracture of 2D problems with injected anisotropic softening elements , 2014 .

[43]  G. Gao,et al.  Application of Digital Image Correlation (DIC) in Dynamic Notched Semi-Circular Bend (NSCB) Tests , 2015 .

[44]  T. Rabczuk,et al.  A Meshfree Method based on the Local Partition of Unity for Cohesive Cracks , 2007 .

[45]  M. Cai,et al.  Fracture Initiation and Propagation in a Brazilian Disc with a Plane Interface: a Numerical Study , 2013, Rock Mechanics and Rock Engineering.

[46]  T. Belytschko,et al.  A three dimensional large deformation meshfree method for arbitrary evolving cracks , 2007 .

[47]  Timon Rabczuk,et al.  Steiner-point free edge cutting of tetrahedral meshes with applications in fracture , 2017 .

[48]  H. Nguyen-Xuan,et al.  A simple and robust three-dimensional cracking-particle method without enrichment , 2010 .

[49]  T. Belytschko,et al.  A comparative study on finite element methods for dynamic fracture , 2008 .

[50]  K. T. Chau,et al.  Analysis of crack coalescence in rock-like materials containing three flaws—Part I: experimental approach , 2001 .

[51]  Herbert H. Einstein,et al.  Fracture coalescence in rock-type materials under uniaxial and biaxial compression , 1998 .

[52]  Khader M. Hamdia,et al.  Uncertainty quantification of the fracture properties of polymeric nanocomposites based on phase field modeling , 2015 .

[53]  Khader M. Hamdia,et al.  Stochastic analysis of the fracture toughness of polymeric nanoparticle composites using polynomial chaos expansions , 2017, International Journal of Fracture.

[54]  Phill-Seung Lee,et al.  Phantom-node method for shell models with arbitrary cracks , 2012 .