We investigate the potential of quasi-Newton methods in facilitating convergence of monolithic solution schemes for phase field fracture modelling. Several paradigmatic boundary value problems are addressed, spanning the fields of quasi-static fracture, fatigue damage and dynamic cracking. The finite element results obtained reveal the robustness of quasi-Newton monolithic schemes, with convergence readily attained under both stable and unstable cracking conditions. Moreover, since the solution method is unconditionally stable, very significant computational gains are observed relative to the widely used staggered solution schemes. In addition, a new adaptive time increment scheme is presented to further reduces the computational cost while allowing to accurately resolve sudden changes in material behavior, such as unstable crack growth. Computation times can be reduced by several orders of magnitude, with the number of load increments required by the corresponding staggered solution being up to 3000 times higher. Quasi-Newton monolithic solution schemes can be a key enabler for large scale phase field fracture simulations. Implications are particularly relevant for the emerging field of phase field fatigue, as results show that staggered cycle-by-cycle calculations are prohibitive in mid or high cycle fatigue. The finite element codes are available to download from www.empaneda.com/codes.

We present a gradient-based theoretical framework for predicting hydrogen assisted fracture in elastic-plastic solids. The novelty of the model lies in the combination of: (i) stress-assisted diffusion of solute species, (ii) strain gradient plasticity, and (iii) a hydrogen-sensitive phase field fracture formulation, inspired by first principles calculations. The theoretical model is numerically

The phase field paradigm, in combination with a suitable variational structure, has opened a path for using Griffith’s energy balance to predict the fracture of solids. These so-called phase field fracture methods have gained significant popularity over the past decade, and are now part of commercial finite element packages and engineering fitness- for-service assessments. Crack paths can be predicted, in arbitrary geometries and dimensions, based on a global energy minimization—without the need for ad hoc criteria. In this work, we review the fundamentals of phase field fracture methods and examine their capabilities in delivering predictions in agreement with the classical fracture mechanics theory pioneered by Griffith. The two most widely used phase field fracture models are implemented in the context of the finite element method, and several paradigmatic boundary value problems are addressed to gain insight into their predictive abilities across all cracking stages; both the initiation of growth and stable crack propagation are investigated. In addition, we examine the effectiveness of phase field models with an internal material length scale in capturing size effects and the transition flaw size concept. Our results show that phase field fracture methods satisfactorily approximate classical fracture mechanics predictions and can also reconcile stress and toughness criteria for fracture. The accuracy of the approximation is however dependent on modelling and constitutive choices; we provide a rationale for these differences and identify suitable approaches for delivering phase field fracture predictions that are in good agreement with well-established fracture mechanics paradigms. This article is part of a discussion meeting issue ‘A cracking approach to inventing new tough materials: fracture stranger than friction’.

The phase field fracture method has emerged as a promising computational tool for modelling a variety of problems including, since recently, hydrogen embrittlement and stress corrosion cracking. In this work, we demonstrate the potential of phase field-based multi-physics models in transforming the engineering assessment and design of structural components in hydrogencontaining environments. First, we present a theoretical and numerical framework coupling deformation, diffusion and fracture, which accounts for inertia effects.Several constitutive choices are considered for the crack density function, including choices with and without an elastic phase in the damage response. The material toughness is defined as a function of the hydrogen content using an atomistically-informed hydrogen degradation law. The model is numerically implemented in 2D and 3D using the finite element method. The resulting computational framework is used to address a number of case studies of particular engineering interest. These are intended

No abstract

No abstract

No abstract

No abstract