An efficient and robust staggered algorithm applied to the quasi-static description of brittle fracture by a phase-field approach

Abstract The phase field method has been widely adopted in brittle fracture analysis for its ability to handle complex crack topology. This paper presents a novel efficient and robust phase field algorithm for quasi-static brittle fracture analysis. This algorithm overcomes two major issues that affect significantly the numerical cost of the method: the treatment of discontinuous crack propagation and the inequality constraint associated with the irreversibility of the damage evolution. To handle discontinuous crack propagation, a semi-implicit scheme, which combines the usual explicit and implicit schemes, is proposed. Different from explicit schemes that require small time steps and purely implicit schemes that lose immediately efficiency when encountering discontinuous propagation, the proposed method can alleviate the steps constraint while keeping a good robustness with discontinuous cracking. Concerning the irreversibility constraint, this work proposes a practical and easy-to-implement method. It is shown that this method is extremely efficient and robust without any supplementary numerical coefficient. The efficiency of the method is demonstrated by means of representative numerical examples.

[1]  A. A. Griffith.  The Phenomena of Rupture and Flow in Solids , 1921 .

[2]  Gilles A. Francfort,et al.  Revisiting brittle fracture as an energy minimization problem , 1998 .

[3]  Jean-Jacques Marigo,et al.  Initiation of a periodic array of cracks in the thermal shock problem: a gradient damage modeling , 2014 .

[4]  Interacting Cracks Obey a Multiscale Attractive to Repulsive Transition. , 2018, Physical review letters.

[5]  Julien Yvonnet,et al.  A phase-field method for computational modeling of interfacial damage interacting with crack propagation in realistic microstructures obtained by microtomography , 2016 .

[6]  D. Mumford,et al.  Optimal approximations by piecewise smooth functions and associated variational problems , 1989 .

[7]  Loïc Daridon,et al.  Artificial compliance inherent to the intrinsic cohesive zone models: criteria and application to planar meshes , 2012, International Journal of Fracture.

[8]  B. Bourdin,et al.  Numerical experiments in revisited brittle fracture , 2000 .

[9]  M. Wheeler,et al.  An augmented-Lagrangian method for the phase-field approach for pressurized fractures , 2014 .

[10]  I. Babuska,et al.  The Partition of Unity Method , 1997 .

[11]  Anthony Gravouil,et al.  2D and 3D Abaqus implementation of a robust staggered phase-field solution for modeling brittle fracture , 2017 .

[12]  K. Saanouni,et al.  A 3D h-adaptive local remeshing technique for simulating the initiation and propagation of cracks in ductile materials , 2018 .

[13]  A. Combescure,et al.  A selective mass scaling method for shear wave propagation analyses in nearly incompressible materials , 2017 .

[14]  S. Melin.  Why do cracks avoid each other? , 1983 .

[15]  Bruno Michel,et al.  Introducing the open-source mfront code generator: Application to mechanical behaviours and material knowledge management within the PLEIADES fuel element modelling platform , 2015, Comput. Math. Appl..

[16]  Christopher J. Larsen,et al.  A time-discrete model for dynamic fracture based on crack regularization , 2011 .

[17]  Patrick E. Farrell,et al.  Linear and nonlinear solvers for variational phase‐field models of brittle fracture , 2015, 1511.08463.

[18]  Yongxing Shen,et al.  Massive parallelization of the phase field formulation for crack propagation with time adaptivity , 2016 .

[19]  Giovanni Lancioni,et al.  The Variational Approach to Fracture Mechanics. A Practical Application to the French Panthéon in Paris , 2009 .

[20]  Yulong Shao,et al.  Adaptive consistent element-free Galerkin method for phase-field model of brittle fracture , 2019, Computational Mechanics.

[21]  L. Ambrosio,et al.  Approximation of functional depending on jumps by elliptic functional via t-convergence , 1990 .

[22]  Christopher J. Larsen,et al.  EXISTENCE OF SOLUTIONS TO A REGULARIZED MODEL OF DYNAMIC FRACTURE , 2010 .

[23]  T. Belytschko,et al.  Non‐planar 3D crack growth by the extended finite element and level sets—Part I: Mechanical model , 2002 .

[24]  L. Lorenzis,et al.  Phase-field modeling of ductile fracture , 2015, Computational Mechanics.

[25]  Endre Süli,et al.  An Adaptive Finite Element Approximation of a Variational Model of Brittle Fracture , 2010, SIAM J. Numer. Anal..

[26]  L. Xia,et al.  Topology optimization for maximizing the fracture resistance of quasi-brittle composites , 2018 .

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

[28]  Thomas J. R. Hughes,et al.  A higher-order phase-field model for brittle fracture: Formulation and analysis within the isogeometric analysis framework , 2014 .

[29]  Gaël Gibert,et al.  A 3D automatic mesh refinement X-FEM approach for fatigue crack propagation , 2019, Finite Elements in Analysis and Design.

[30]  Thomas Wick,et al.  Modified Newton methods for solving fully monolithic phase-field quasi-static brittle fracture propagation , 2017 .

[31]  Isabelle Ramière,et al.  Iterative residual-based vector methods to accelerate fixed point iterations , 2015, Comput. Math. Appl..

[32]  Mary F. Wheeler,et al.  A primal-dual active set method and predictor-corrector mesh adaptivity for computing fracture propagation using a phase-field approach , 2015 .

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

[34]  Chengbao Jiang,et al.  A study of the mechanism of formation and numerical simulations of crack patterns in ceramics subjected to thermal shock , 2012 .

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

[36]  T. Belytschko,et al.  Non‐planar 3D crack growth by the extended finite element and level sets—Part II: Level set update , 2002 .

[37]  Ted Belytschko,et al.  A finite element method for crack growth without remeshing , 1999 .

[38]  L. De Lorenzis,et al.  On penalization in variational phase-field models of brittle fracture , 2018, Computer Methods in Applied Mechanics and Engineering.

[39]  Indra Vir Singh,et al.  An adaptive multiscale phase field method for brittle fracture , 2018 .

[40]  Fei Zhang,et al.  Moving mesh finite element simulation for phase-field modeling of brittle fracture and convergence of Newton's iteration , 2017, J. Comput. Phys..

[41]  Yongxing Shen,et al.  Variational h-adaption method for the phase field approach to fracture , 2019, International Journal of Fracture.

[42]  Jean-Jacques Marigo,et al.  Regularized formulation of the variational brittle fracture with unilateral contact: Numerical experiments , 2009 .

[43]  M. Elices,et al.  The cohesive zone model: advantages, limitations and challenges , 2002 .

[44]  I. Ramière,et al.  A new numerical methodology for simulation of unstable crack growth in time independent brittle materials , 2017 .