Summary
We present a finite element (FE) formulation for the free vibration analysis of doubly curved laminated composite and sandwich shells having multiple delaminations, employing a facet shell element based on the efficient third‐order zigzag theory and the region approach of modeling delaminations. The methodology, hitherto not attempted, is general for delaminations occurring at multiple interfacial and spatial locations. A recently developed hybrid method is used for satisfying the continuity of the nonlinear layerwise displacement field at the delamination fronts. The formulation is shown to yield very accurate results with reference to full‐field three‐dimensional FE solutions, for the natural frequencies and mode shapes of delaminated shallow and deep, composite and highly inhomogeneous soft‐core sandwich shells of different geometries and boundary conditions, with a significant computational advantage. The accuracy is sensitive to the continuity method used at the delamination fronts, the usual point continuity method yielding rather poor accuracy, and the proposed hybrid method giving the best accuracy. Such efficient modeling of laminated shells with delamination damage will be of immense use for model‐based techniques for structural health monitoring of laminated shell‐type structures.