This work presents a new plate theory for out-of-plane loaded thick plates where the static unknowns are those of the Love-Kirchhoff theory, to which six components are added representing the gradient of the bending moment. The Bending-Gradient theory is an extension to arbitrary multilayered plates of the Reissner-Mindlin theory which appears as a special case when the plate is homogeneous. The new theory is applied to multilayered plates and its predictions are compared to full 3D Pagano's exact solutions and other approaches. It gives good predictions of both deflection and shear stress distributions in any material configuration. Moreover, under some symmetry conditions, the Bending-Gradient model coincides with the second-order approximation of the exact solution as the slenderness ratio L/h goes to infinity
