PhaseFieldFEniCS FEniCS Python script with a staggered implementation of the phase field fracture method, suitable for 2D and 3D case studies. Includes a document with detailed instructions. Paper: Hirshikesh, S. Natarajan, R. K. Annabattula, E. Martínez-Pañeda. Phase field modelling of crack propagation in functionally graded materials. Composites Part B: Engineering 169: 239-248 (2019)