This work presents the effectiveness of differential quadrature shape functions (i.e., Lagrange interpolation polynomial, Cardinal sine function, Delta Lagrange kernel and Regularized Shannon kernel) in the solution of nonlinear vibration of multilayers piezoelectric plates with nonlinear elastic support. A piezoelectric composite laminated plate is rested on nonlinear Winkler and Visco-Pasternak elastic foundations problems. Based on 3D elasticity theory and piezoelectricity, the governing equations of motion are derived. Differential quadrature methods based on four shape functions are presented as numerical techniques for solving this problem. The perturbation method is implemented to solve the obtained nonlinear eigenvalue problem. A MATLAB code is written for each technique for solving this problem and extract the numerical results. To validate these methods, the computed results are we compare with the previous exact results. In addition, parametric analyses are offered to investigate the influence of length to thickness ratio, elastic foundation parameters, various boundary conditions, and piezoelectric layers thickness on the natural frequencies and mode shapes. Consequently, it is discovered that the obtained results via the proposed schemes can be applied in structural health monitoring.