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.