A crack-length control technique for phase field fracture in FFT homogenization

TL;DR

This paper introduces a novel crack-length control technique for phase field fracture in FFT homogenization, enabling stable, monotonic crack propagation simulations and improving numerical performance in micromechanical fracture modeling.

Contribution

It proposes a new crack control method compatible with FFT and finite element methods, resolving unstable crack propagation and snap-backs in phase field fracture simulations.

Findings

The method achieves stable crack propagation with superposed stress-strain curves.

J-integral calculations show constant energy release rate during propagation.

Application to composites and porous materials estimates effective fracture toughness.

Abstract

Modeling the propagation of cracks at the microscopic level is fundamental to understand the effect of the microstructure on the fracture process. Nevertheless, microscopic propagation is often unstable and when using phase field fracture poor convergence is found or, in the case of using staggered algorithms, leads to the presence of jumps in the evolution of the cracks. In this work, a novel method is proposed to perform micromechanical simulations with phase field fracture imposing monotonic increases of crack length and allowing the use of monolithic implementations, being able to resolve all the snap-backs during the unstable propagation phases. The method is derived for FFT based solvers in order to exploit its very high numerical performance n micromechanical problems, but an equivalent method is also developed for Finite Elements (FE) showing the equivalence of both implementations. It is shown that the stress-strain curves and the crack paths obtained using the crack control method are superposed in stable propagation regimes to those obtained using strain control with a staggered scheme. J-integral calculations confirm that during the propagation process in the crack control method, the energy release rate remains constant and equal to an effective fracture energy that has been determined as function of the discretization for FFT simulations. Finally, to show the potential of the method, the technique is applied to simulate crack propagation through the microstructure of composites and porous materials providing an estimation of the effective fracture toughness.