Abstract. The geometrically non-linear mechanical response of multilayer composite shells is addressed via an innovative discontinuous Galerkin formulation. In the framework of the Carrera Unified formulation, equivalent single layer kinematics with different through-the-thickness accuracy is adopted. The variational statement governing the shell nonlinear behavior is derived. The corresponding governing equations are solved via a discontinuous Galerkin approach, which employs the pure penalty method to weakly enforce the connection between the mesh elements. Numerical tests are presented to show the capabilities of the proposed approach.