Thermoplastic composites are continuously replacing thermosetting composites in lightweight structures. However, the accomplished work on the fatigue behavior of thermoplastics is quite limited. In the present work, we propose a numerical modeling approach for simulating fatigue delamination growth and predicting the residual tensile strength of quasi-isotropic TC 1225 LM PAEK thermoplastic coupons. The approach is supported and validated by tension and fatigue (non-interrupted and interrupted) tests. Fatigue delamination growth is simulated using a mixed-mode fatigue crack growth model, which is based on the cohesive zone modeling method. All interfaces between the thermoplastic plies were modelled using cohesive elements. Quasi-static tension analyses on pristine and fatigued coupons were performed using a progressive damage model implementing a set of Hashin-type strain-based failure criteria and a damage mechanics-based material property degradation module. Using the fatigue model, delamination growth was predicted as a function of number of cycles at all interfaces. The results agree well with C-scan images taken on fatigued coupons during interruptions of fatigue tests. An unequal and unsymmetric delamination growth was predicted due to the quasi-isotropic lay-up. Moreover, the combined models capture the decrease in the residual tensile strength of the coupons. In the quasi-static tension analysis of the fatigued coupons, the driving failure mechanisms are the fast propagation of the pre-existing delamination and the severe matrix cracking.