Fast implicit solvers for phase-field fracture problems on heterogeneous microstructures

Abstract We study fast and memory-efficient FFT-based implicit solution methods for small-strain phase-field crack problems for microstructured brittle materials. A fully implicit first order formulation of the problem coupling elasticity and damage permits using comparatively few, but large, time steps compared to semi-explicit schemes. We investigate memory-efficient FFT-based solution techniques, and identify the heavy ball scheme as particularly powerful. We discuss the memory-efficient implementation and present demonstrative numerical examples.

[1]  Klaus Hackl,et al.  A Method for Gradient Enhancement of Continuum Damage Models , 2008 .

[2]  C. M. Reeves,et al.  Function minimization by conjugate gradients , 1964, Comput. J..

[3]  Radhi Abdelmoula,et al.  A damage model for crack prediction in brittle and quasi-brittle materials solved by the FFT method , 2012, International Journal of Fracture.

[4]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[5]  Songhe Meng,et al.  A non-local fracture model for composite laminates and numerical simulations by using the FFT method , 2012 .

[6]  Felix Fritzen,et al.  Fourier-Accelerated Nodal Solvers (FANS) for homogenization problems , 2018 .

[7]  Gianni Royer-Carfagni,et al.  Regularized variational theories of fracture: A unified approach , 2010 .

[8]  J. Michel,et al.  Analysis of Inhomogeneous Materials at Large Strains using Fast Fourier Transforms , 2003 .

[9]  R. Dembo,et al.  INEXACT NEWTON METHODS , 1982 .

[10]  A. Hartmaier,et al.  Formulation of nonlocal damage models based on spectral methods for application to complex microstructures , 2015 .

[11]  Homer F. Walker,et al.  Choosing the Forcing Terms in an Inexact Newton Method , 1996, SIAM J. Sci. Comput..

[12]  Lázaro M. Sixto-Camacho,et al.  Asymptotic homogenization of periodic thermo-magneto-electro-elastic heterogeneous media , 2013, Comput. Math. Appl..

[13]  Laura De Lorenzis,et al.  A phase-field model for ductile fracture at finite strains and its experimental verification , 2016 .

[14]  Michael I. Jordan,et al.  Accelerated Gradient Descent Escapes Saddle Points Faster than Gradient Descent , 2017, COLT.

[15]  T. Rabczuk,et al.  Phase-field modeling of fracture in linear thin shells , 2014 .

[16]  Homer F. Walker,et al.  Anderson Acceleration for Fixed-Point Iterations , 2011, SIAM J. Numer. Anal..

[17]  Matti Schneider,et al.  On polarization-based schemes for the FFT-based computational homogenization of inelastic materials , 2019, Computational Mechanics.

[18]  John D. Clayton,et al.  A geometrically nonlinear phase field theory of brittle fracture , 2014, International Journal of Fracture.

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

[20]  D. Jeulin,et al.  A Fourier-based numerical homogenization tool for an explosive material , 2015 .

[21]  R. Müller,et al.  On phase field modeling of ductile fracture , 2016 .

[22]  M. Schneider,et al.  FFT‐based homogenization for microstructures discretized by linear hexahedral elements , 2017 .

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

[24]  Antonin Chambolle,et al.  Approximation of a Brittle Fracture Energy with a Constraint of Non-interpenetration , 2017, 1709.04208.

[25]  Irene Arias,et al.  Phase-field modeling of crack propagation in piezoelectric and ferroelectric materials with different electromechanical crack conditions , 2012 .

[26]  L. Dormieux,et al.  Combining Galerkin approximation techniques with the principle of Hashin and Shtrikman to derive a new FFT-based numerical method for the homogenization of composites , 2012 .

[27]  Boris Polyak.  Some methods of speeding up the convergence of iteration methods , 1964 .

[28]  S Torquato,et al.  Robust algorithm to generate a diverse class of dense disordered and ordered sphere packings via linear programming. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[29]  Milan Jirásek,et al.  Nonlocal models for damage and fracture: Comparison of approaches , 1998 .

[30]  Euhanna Ghadimi,et al.  Global convergence of the Heavy-ball method for convex optimization , 2014, 2015 European Control Conference (ECC).

[31]  Felix Ernesti,et al.  An FFT‐based solver for brittle fracture on heterogeneous microstructures , 2019 .

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

[33]  Tymofiy Gerasimov,et al.  A line search assisted monolithic approach for phase-field computing of brittle fracture , 2016 .

[34]  Hervé Moulinec,et al.  A computational scheme for linear and non‐linear composites with arbitrary phase contrast , 2001 .

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

[36]  Jean-Jacques Marigo,et al.  The issues of the uniqueness and the stability of the homogeneous response in uniaxial tests with gradient damage models , 2011 .

[37]  B. Nestler,et al.  Phase-field modeling of crack propagation in multiphase systems , 2016 .

[38]  Yang Chen,et al.  A FFT solver for variational phase-field modeling of brittle fracture , 2019, Computer Methods in Applied Mechanics and Engineering.

[39]  Laura De Lorenzis,et al.  A review on phase-field models of brittle fracture and a new fast hybrid formulation , 2015 .

[40]  Ulrich Langer,et al.  Matrix-free multigrid solvers for phase-field fracture problems , 2019, Computer Methods in Applied Mechanics and Engineering.

[41]  M. Schneider,et al.  Efficient fixed point and Newton–Krylov solvers for FFT-based homogenization of elasticity at large deformations , 2014 .

[42]  Shenyang Y. Hu,et al.  Simulation of damage evolution in composites : A phase-field model , 2009 .

[43]  F. Willot,et al.  Fourier-based schemes for computing the mechanical response of composites with accurate local fields , 2014, 1412.8398.

[44]  M. Schneider.  The sequential addition and migration method to generate representative volume elements for the homogenization of short fiber reinforced plastics , 2017 .

[45]  D. Gross,et al.  Fracture simulation of ferroelectrics based on the phase field continuum and a damage variable , 2010 .

[46]  Matti Schneider,et al.  An FFT-based fast gradient method for elastic and inelastic unit cell homogenization problems , 2017 .

[47]  Boris Polyak.  The conjugate gradient method in extremal problems , 1969 .

[48]  B. Bourdin.  Numerical implementation of the variational formulation for quasi-static brittle fracture , 2007 .

[49]  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 .

[50]  Cv Clemens Verhoosel,et al.  A phase‐field model for cohesive fracture , 2013 .

[51]  Charlotte Kuhn,et al.  Numerical and Analytical Investigation of a Phase Field Model for Fracture , 2013 .

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

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

[54]  Ralf Müller,et al.  On degradation functions in phase field fracture models , 2015 .

[55]  J. Marigo,et al.  Gradient Damage Models and Their Use to Approximate Brittle Fracture , 2011 .

[56]  M. Ortiz,et al.  Modeling fracture by material-point erosion , 2013, International Journal of Fracture.

[57]  Ya-Xiang Yuan,et al.  A Nonlinear Conjugate Gradient Method with a Strong Global Convergence Property , 1999, SIAM J. Optim..

[58]  P. Kowalczyk,et al.  Calculation of cancellous bone elastic properties with the polarization‐based FFT iterative scheme , 2017, International journal for numerical methods in biomedical engineering.

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

[60]  L. Lorenzis,et al.  Phase-field modeling of brittle and ductile fracture in shells with isogeometric NURBS-based solid-shell elements , 2016 .

[61]  A. Giacomini,et al.  A Discontinuous Finite Element Approximation of Quasi-static Growth of Brittle Fractures , 2003, math/0304224.

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

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

[64]  M. Diehl,et al.  A spectral method solution to crystal elasto-viscoplasticity at finite strains , 2013 .

[65]  M. Saunders,et al.  Solution of Sparse Indefinite Systems of Linear Equations , 1975 .

[66]  M. Schneider.  On the Barzilai‐Borwein basic scheme in FFT‐based computational homogenization , 2019, International Journal for Numerical Methods in Engineering.

[67]  Adrian Willenbücher,et al.  Phase field approximation of dynamic brittle fracture , 2014 .

[68]  Daniel Kienle,et al.  Phase field modeling of fracture in anisotropic brittle solids , 2017 .

[69]  Bin Li,et al.  Phase‐field modeling and simulation of fracture in brittle materials with strongly anisotropic surface energy , 2015 .

[70]  Hervé Moulinec,et al.  A numerical method for computing the overall response of nonlinear composites with complex microstructure , 1998, ArXiv.

[71]  G. Bonnet,et al.  A polarization‐based FFT iterative scheme for computing the effective properties of elastic composites with arbitrary contrast , 2012 .

[72]  Renald Brenner,et al.  Numerical computation of the response of piezoelectric composites using Fourier transform , 2009 .

[73]  Luc Dormieux,et al.  FFT-based methods for the mechanics of composites: A general variational framework , 2010 .

[74]  Fernando Fraternali,et al.  Eigenfracture: An Eigendeformation Approach to Variational Fracture , 2009, Multiscale Model. Simul..

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

[76]  G. Piero,et al.  A variational model for fracture mechanics - Numerical experiments , 2007 .

[77]  Jan Novák,et al.  Accelerating a FFT-based solver for numerical homogenization of periodic media by conjugate gradients , 2010, J. Comput. Phys..

[78]  M. Schneider,et al.  Computational homogenization of elasticity on a staggered grid , 2016 .

[79]  Graeme W. Milton,et al.  A fast numerical scheme for computing the response of composites using grid refinement , 1999 .

[80]  Jean-Jacques Marigo,et al.  From the onset of damage to rupture: construction of responses with damage localization for a general class of gradient damage models , 2013 .

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

[82]  Xue Zhang,et al.  Numerical evaluation of the phase-field model for brittle fracture with emphasis on the length scale , 2017 .

[83]  Eric Lorentz,et al.  Gradient damage models: Toward full-scale computations , 2011 .

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

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

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

[87]  L. Gélébart,et al.  Non-linear extension of FFT-based methods accelerated by conjugate gradients to evaluate the mechanical behavior of composite materials , 2013 .

[88]  Jan Zeman,et al.  Guaranteed upper-lower bounds on homogenized properties by FFT-based Galerkin method , 2014, 1404.3614.

[89]  M. Focardi,et al.  Phase field approximation of cohesive fracture models , 2014, 1405.6883.

[90]  H. Moulinec,et al.  A fast numerical method for computing the linear and nonlinear mechanical properties of composites , 1994 .

[91]  R. Müller,et al.  A multiscale approach for modeling progressive damage of composite materials using fast Fourier transforms , 2014 .

[92]  Vinh Phu Nguyen,et al.  A length scale insensitive phase-field damage model for brittle fracture , 2018, Journal of the Mechanics and Physics of Solids.

[93]  C. Sauder,et al.  Analysis of the damage initiation in a SiC/SiC composite tube from a direct comparison between large-scale numerical simulation and synchrotron X-ray micro-computed tomography , 2019, International Journal of Solids and Structures.

[94]  F. Willot,et al.  Fourier‐based schemes with modified Green operator for computing the electrical response of heterogeneous media with accurate local fields , 2013, 1307.1015.

[95]  E. Süli,et al.  An adaptive finite element approximation of a generalized Ambrosio-Tortorelli functional , 2013 .