Use of higher-order Legendre polynomials for multilayered plate elements with node-dependent kinematics