Abstract:The simulation of fracture using continuum ductile damage models attains a pathological discretization dependence caused by <a href="/topics/engineering/strain-localization">strain localization</a>, after loss of <a href="/topics/engineering/ellipticity">ellipticity</a> of the problem, in regions whose size is connected to the <a href="/topics/engineering/spatial-discretization">spatial discretization</a>. Implicit gradient techniques suppress this problem introducing some inelastic non-local fields and solving an enriched formulation where the classical balance of <a href="/topics/engineering/linear-momentum">linear momentum</a> is fully coupled with a Helmholtz-type equation for each of the non-local variable. Such Helmholtz-type equations determine the distribution of the non-local fields in bands whose width is controlled by a characteristic length, independently on the spatial discretization. The numerical resolution of this coupled problem using the <a href="/topics/engineering/finite-element-method">Finite Element method</a> is computationally very expensive and its use to simulate the damage process in 3D multi-phase microstructures becomes prohibitive.In this work, we propose a novel FFT-based iterative algorithm for simulating gradient ductile damage in <a href="/topics/engineering/computational-homogenisation">computational homogenization</a> problems. In particular, the Helmholtz-type equation of the implicit <a href="/topics/engineering/gradient-approach">gradient approach</a> is properly generalized to model the <a href="/topics/engineering/regularization">regularization</a> of damage in multi-phase media, where multiple damage variables and different characteristic lengths may come into play. In the proposed iterative algorithm, two distinct problems are solved in a staggered fashion: (i) a conventional mechanical problem via a FFT-Galerkin solver with mixed macroscopic loading control and (ii) the generalized Helmholtz-type equation using a Krylov-based algorithm combined with an efficient pre-conditioner. The numerical implementation is firstly validated on simple two-dimensional microstructures, showing identical responses for different spatial discretizations and reproducing a ductility change dependent on the characteristic length. Finally, the robustness and efficiency of the algorithm is demonstrated in the simulation of failure of complex 3D particle reinforced composites characterized by millions of degrees of freedom.

A simple and efficient implementation of localizing gradient damage method in COMSOL for fracture simulation

A Computationally Efficient Vectorized Implementation of Localizing Gradient Damage Method in MATLAB

GPU-accelerated Approach for 2D Fracture Analysis of Structures Combining Finite Particle Method and Cohesive Zone Model

A simple implementation of localizing gradient damage model in Abaqus

Phase field modeling of quasi-static and dynamic crack propagation: COMSOL implementation and case studies

An eXtended Finite Element Method implementation in COMSOL Multiphysics: Solid Mechanics

Thermodynamically consistent linear-gradient damage model in Abaqus

Efficient quasi-brittle fracture simulations of concrete at mesoscale using micro CT images and a localizing gradient damage model

Phase-field simulation of crack propagation in quasi-brittle materials: COMSOL implementation and parameter sensitivity analysis

An FFT framework for simulating non-local ductile failure in heterogeneous materials

An efficient nonlocal treatment for damage modeling based on the background mesh concept

A Cohesive Fracture Model for Discrete Element Method Based on Polyhedral Blocks

A Nonlocal Finite Element Model for Progressive Failure Analysis of Composite Laminates

A local to global (L2G) finite element method for efficient and robust analysis of arbitrary cracking in 2D solids

A Nonlocal Treatment Technique Based on the Background Cell Concept for Micro-Mechanical Damage Modeling

Image-based adaptive domain decomposition for continuum damage models

Moving from theory to application: Evaluating the numerical implementation of void shape effects and damage delocalization in the modeling of ductile fracture in porous plastic metals

An eXtended finite element method implementation in COMSOL multiphysics: Thermo-hydro-mechanical modeling of fluid flow in discontinuous porous media

A fully-coupled computational framework for large-scale simulation of fluid-driven fracture propagation on parallel computers

A Unified FFT Framework with Coupled Gradient Damage Model and Phase-Field Method for the Failure Analysis of Composites

A Non-Local Fracture Model for Composite Laminates and Numerical Simulations by Using the FFT Method