We present a generalised phase field-based formulation for predicting fatigue crack growth in metals. The theoretical framework aims at covering a wide range of material behaviour. Different fatigue degradation functions are considered and their influence is benchmarked against experiments. The phase field constitutive theory accommodates the so-called AT1, AT2 and phase field-cohesive zone (PF-CZM) models. In regards to material deformation, both non-linear kinematic and isotropic hardening are considered, as well as the combination of the two. Moreover, a monolithic solution scheme based on quasi-Newton algorithms is presented and shown to significantly outperform staggered approaches. The potential of the computational framework is demonstrated by investigating several 2D and 3D boundary value problems of particular interest. Constitutive and numerical choices are compared and insight is gained into their differences and similarities. The framework enables predicting fatigue crack growth in arbitrary geometries and for materials exhibiting complex (cyclic) deformation and damage responses. The finite element code developed is made freely available at www.empaneda.com/codes. •We present a generalised, phase field-based fatigue model for elasto-plastic solids.•Cyclic deformation is modelled by a combined non-linear kinematic/isotropic hardening law.•Three classes of phase field models are considered: AT1, AT2 and PF-CZM.•A quasi-Newton algorithm is used to enable a robust and efficient monolithic scheme.•The potential of the model is showcased with paradigmatic 2D and 3D case studies.