Abstract: Free vibration response of laminated composite shells with delamination is presented using the finite element method based on first order shear deformation theory. The shell theory used is the extension of dynamic, shear deformable theory according to the Sanders'' first approximation for doubly curved shells, which can be reduced to Love''s and Donnell''s theories by means of tracers. An eight-noded C 0 continuity, isoparametric quadrilateral element with five degrees of freedom per node is used in the formulation. For modeling the delamination, multipoint constraint algorithm is incorporated in the finite element code. The natural frequencies of the delaminated cylindrical (CYL), spherical (SPH) and hyperbolic paraboloid (HYP) shells are determined by using the above mentioned shell theories, namely Sanders'', Love''s, and Donnell''s. The validity of the present approach is established by comparing the authors'' results with those available in the literature. Additional studies on free vibration response of CYL, SPH and HYP shells are conducted to assess the effects of delamination size and number of layers considering all three shell theories. It is shown that shell theories according to Sanders and Love always predict practically identical frequencies. Donnell''s theory gives reliable results only for shallow shells. Moreover, the natural frequency is found to be very sensitive to delamination size and number of layers in the shell. [Copyright &y& Elsevier]