  Phase field fracture models have seen widespread application in the last
decade. Among these applications, its use to model the evolution of fatigue
cracks has attracted particular interest, as fatigue damage behaviour can be
predicted for arbitrary loading histories, dimensions and complexity of the
cracking phenomena at play. However, while cycle-by-cycle calculations are
remarkably flexible, they are also computationally expensive, hindering the
applicability of phase field fatigue models for technologically-relevant
problems. In this work, a computational framework for accelerating phase field
fatigue calculations is presented. Two novel acceleration strategies are
proposed, which can be used in tandem and together with other existing
acceleration schemes from the literature. The computational performance of the
proposed methods is documented through a series of 2D and 3D boundary value
problems, highlighting the robustness and efficiency of the framework even in
complex fatigue problems. The observed reduction in computation time using both
of the proposed methods in tandem is shown to reach a speed-up factor of 32,
with a scaling trend enabling even greater reductions in problems with more
load cycles.
