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

TL;DR

This paper introduces a fast FFT-based iterative method for simulating non-local ductile failure in heterogeneous materials, effectively handling complex 3D microstructures with multiple damage variables and characteristic lengths.

Contribution

It develops a novel FFT framework that efficiently solves generalized Helmholtz equations for gradient damage modeling in multi-phase media, improving computational feasibility.

Findings

Validated the numerical implementation against known solutions.

Demonstrated robustness and efficiency in simulating failure of complex 3D composites.

Achieved significant computational speed-up over traditional finite element methods.

Abstract

The simulation of fracture using continuum ductile damage models attains a pathological discretization dependence caused by strain localization, after loss of ellipticity of the problem, in regions whose size is connected to the spatial discretization. Implicit gradient techniques suppress this problem introducing some inelastic non-local fields and solving an enriched formulation where the classical balance of linear momentum 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 Finite Element method 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 computational homogenization problems. In particular, the Helmholtz-type equation of the implicit gradient approach is properly generalized to model the regularization 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. 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.