This paper presents a nonlinear model for piezoelastic laminated plates with damage effect of the intra-layers and inter-laminar interfaces. Discontinuity of displacement and electric potential on the interfaces are depicted by three shape functions. By using the Hamilton variation principle, the three-dimensional nonlinear dynamic equations of piezoelastic laminated plates with damage effect are derived. Then, by using the Galerkin method, a mathematical solution is presented. In the numerical studies, effects of various factors on the natural frequencies and nonlinear amplitude-frequency response of the simply-supported peizoelastic laminated plates with interfacial imperfections are discussed. These factors include different damage models, thickness of the piezoelectric layer, side-to-thickness
ratio, and length-to-width ratio.