A deeper understanding of delamination growth under fatigue loading in composite wind turbine blades is required. A new cohesive model for delamination fatigue initiation, which uses initiation S-N curve data to calculate the number of cycles to the introduction of delaminations, is developed and tailored to a state-of-the-art model for fatigue propagation. Both models are implemented in a novel user-defined cohesive element for ANSYS Mechanical APDL. In fatigue propagation-dominated test cases, including with multiple delaminations, the crack growth rate evolution is accurately predicted. In fatigue initiation-dominated cases, a satisfactory prediction is achieved, also showcasing the ability to model fatigue propagation after initiation.