Use of higher-order Legendre polynomials in node-dependent kinematic shell elements