Nonlinear large amplitude vibration and dynamic deformations of sandwich plates with composite face sheets and hyperelastic cores that encounter large thickness changes during the vibration are investigated here for the first time. Moreover, a new global–local zigzag hyperelastic plate theory that considers the incompressibility and large transverse deformability of the core, and employs transverse normal strains whose order is higher than that of the in-plane strains is proposed. The plate is resting on a viscoelastic Winkler-Pasternak substrate, exposed to various spatial and time distributions of loads, and experiencing a diversity of edge conditions. The governing equations of motion are obtained by using Hamilton’s principle, von Kármán’s assumptions, and the series expansion of the volume-preserving strain energy density function. A 2D differential quadrature (2D-DQ) technique is employed for the spatial discretization of a structure with both constitutive and kinematic nonlinearities, for the first time. The time-dependent combination of the nonlinear coupled motion equations and boundary conditions is treated by an iterative/updating Runge-Kutta time-marching technique and the effects of the thickness and aspect ratio of the plate, viscoelastic and shear parameters of the foundation, types of the edge conditions, constitutive parameters of the composite materials, fiber orientation angle, and magnitude and distribution pattern of the loads on the resulting dynamic lateral deflections are studied. Results reveal that while the vibration frequencies and amplitudes of the face sheets are quite different, the damping is more remarkable when stiffer cores are utilized. Moreover, the effects of the higher vibration modes are more apparent in looser edge conditions, smaller fiber angles, and special excitation frequencies.