Implementation of a phase field damage model with a nonlinear evolution equation in an FFT-based solver

Abstract

This paper focuses on the numerical implementation of phase-field models of fracture using the Fast Fourier Transform (FFT) based numerical method. Recent studies on phase-field models focus on the discussions of the choice of the value of regularization length proposed to smear the discontinuity of the sharp crack. Some studies argue that it should be considered as a material property because it has a significant impact on the mechanical behavior of a material in some phase-field models, for instance, in the model proposed by Miehe. However, our results in this study for heterogeneous materials have shown that the choice of regularization length not only affects the macroscopic mechanical behavior but also the local crack propagation patterns. As a result, it can be challenging to select an appropriate value that produces both accurate macroscopic responses and local crack patterns for certain phase-field models, such as Miehe’s model. Thus, the phase-field model proposed by Wu, which has been proven to reduce the length sensitivity for homogeneous material, has been successfully implemented in an FFT-based solver with the application of Newton–Krylov algorithm.The length sensitivity of Wu’s phase-field model for heterogeneous materials is also investigated in this paper. Our tests show that Wu’s phase-field model has partial length sensitivity for heterogeneous materials. However, the main reason for this sensitivity is that one phase enters the damage zone of other phase, which is different from Miehe’s model. As a result, a set of criteria for accurately choosing the regularization value has been established for Wu’s model in this work. Meanwhile, it has also been found that Wu’s model may be more suitable than Miehe’s model for brittle failure due to the introduction of an elastic stage.