Finite Element Method (FEM) based modelling and simulation technique which can suitably interpret and analyze initiation and propagation of delamination damages in laminated curved shell structures, has been developed in the present analysis. The FE model yields result in close agreement with literature. Three-dimensional interlaminar out-of-plane stresses, responsible for damage initiation have been studied in details at different ply-interfaces of the composite spherical shell. Tsai-wu failure criterion has been used to interpret the initiation of delamination damages at different critical ply-interfaces. However, Strain Energy Release Rate (SERR), which is a Fracture Mechanics, based parameter, calculated through Virtual Crack Closure Technique (VCCT) has been used to quantify the damage growth rate. Delamination damages have been predominantly affecting the fundamental frequency of vibration of the curved shells. Effect of different possible ply-orientations on delaminated frequency and damage growth has been studied and pertinent conclusions have been drawn for an improved dynamic and structural response of the laminated curved shell structures. Keywords: Fibre Reinforced Polymer (FRP) composite, Finite Element Method (FEM), Spherical Shells, Strain Energy Release Rate (SERR), Virtual Crack Closure Technique (VCCT).