A robust algorithm for configurational-force-driven brittle crack propagation with R-adaptive mesh alignment

• Published in: International journal for numerical methods in engineering Vol. 72; no. 2; pp. 127 - 155
• Main Authors: Miehe, C., Gürses, E.
• Format: Journal Article
• Language: English
• Published: Chichester, UK John Wiley & Sons, Ltd 08.10.2007; Wiley
• Subjects: Computational techniques; configurational forces; crack simulations; energy minimization; Exact sciences and technology; finite elements; fracture; Fracture mechanics (crack, fatigue, damage...); Fundamental areas of phenomenology (including applications); Mathematical methods in physics; Physics; Solid mechanics; Structural and continuum mechanics; Constitutive equation; Discontinuity; Energy release; Brittle fracture; Brittle material; Minimization; Modeling; Adaptive method; Crack propagation; crack simulations; fracture; Staggered arrangement; Finite element method; Thermodynamics; Alignment; configurational forces; Post critical range; Positive definite matrix; Variational calculus; Griffith crack; Crack tip; Mesh generation; finite elements; energy minimization
• ISSN: 0029-5981; 1097-0207; 1097-0207
• DOI: 10.1002/nme.1999

The paper considers a variational formulation of brittle fracture in elastic solids and proposes a numerical implementation by a finite element method. On the theoretical side, we outline a consistent thermodynamic framework for crack propagation in an elastic solid. It is shown that both the elastic equilibrium response as well as the local crack evolution follow in a natural format by exploitation of a global Clausius–Planck inequality in the sense of Coleman's method. Here, the canonical direction of the crack propagation associated with the classical Griffith criterion is the direction of the material configurational force which maximizes the local dissipation at the crack tip and minimizes the incremental energy release. On the numerical side, we exploit this variational structure in terms of crack‐driving configurational forces. First, a standard finite element discretization in space yields a discrete formulation of the global dissipation in terms configurational nodal forces. As a consequence, the constitutive setting of crack propagation in the space‐discretized finite element context is naturally related to discrete nodes of a typical finite element mesh. Next, consistent with the node‐based setting, the discretization of the evolving crack discontinuity is performed by the doubling of critical nodes and interface segments of the mesh. Critical for the success of this procedure is its embedding into an r‐adaptive crack‐segment reorientation procedure with configurational‐force‐based directional indicator. Here, successive crack releases appear in discrete steps associated with the given space discretization. These are performed by a staggered loading–release algorithm of energy minimization at frozen crack state followed by the successive crack releases at frozen deformation. This constitutes a sequence of positive‐definite discrete subproblems with successively decreasing overall stiffness, providing an extremely robust algorithmic setting in the postcritical range. We demonstrate the performance of the formulation by means of representative numerical simulations. Copyright © 2007 John Wiley & Sons, Ltd.