This paper deals with the study of interaction of manufacturing thermal residual stresses and mechanical loading in penny-shaped delaminations embedded between dissimilar, anisotropic fiber composite layers by conducting two sets of three-dimensional thermoelastic finite element analyses with and without residual stress effects. Modified crack closure integral (MCCI) techniques based on the concepts of linear elastic fracture mechanics (LEFM) have been used to calculate the distribution of individual modes of strain energy release rates (SERR) to investigate the interlaminar delamination initiation and propagation characteristics. Asymmetric variations of strain energy release rates obtained along the delamination front are caused by the overlapping stress fields due to the coupling effect of thermal and mechanical loadings. It is found that parameters such as ply sequence and orientation, thermoelastic anisotropy and material heterogeneity. and ply properties of the delaminated interface dictate the interlaminar fracture behavior of multi-ply laminated FRP composites.