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 was supported and validated by tension and fatigue (non-interrupted and interrupted) tests. Fatigue delamination growth was simulated using a mixed-mode fatigue crack growth model, which was based on the cohesive zone modeling method. Quasi-static tension analyses on pristine and fatigued coupons were performed using a progressive damage model. These analyses were implemented using a set of Hashin-type strain-based failure criteria and a damage mechanics-based material property degradation module. Utilizing the fatigue model, we accurately foretold the expansion of delamination concerning the cycle count across 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 layup. Moreover, the combined models capture the decrease in the residual tensile strength of the coupons. During the quasi-static tension analysis of the fatigued coupons, we observed that the primary driving failure mechanisms were the rapid spread of existing delamination and the consequential severe matrix cracking.