Abstract: In this paper, a new computational framework for the analysis of incompressible Electro Active Polymer (EAP) shells subjected to large strains and large electric fields is presented. Two novelties are incorporated in this work. First, the variational and constitutive frameworks developed by the authors in recent publications (Gil and Ortigosa, 2016; Ortigosa and Gil, 2016; Ortigosa et al., 2016) in the context of three-dimensional electromechanics are particularised/degenerated to the case of geometrically exact shell theory. This formulation is computationally very convenient as EAPs are typically used as thin shell-like components in a vast range of applications. The proposed formulation follows a rotationless description of the kinematics of the shell, enhanced with extra degrees of freedom corresponding to the thickness stretch and the hydrostatic pressure, critical for the consideration of incompressibility. Different approaches are investigated for the interpolation of these extra fields and that of the electric potential across the thickness of the shell. Crucially, this allows for the simulation of multilayer and composite materials, which can display a discontinuous strain distribution across their thickness. As a second novelty, a continuum degenerate approach allows for the consideration of complex three-dimensional electromechanical constitutive models, as opposed to those defined in terms of the main strain measures of the shell. More specifically, convex multi-variable (three-dimensional) constitutive models, complying with the ellipticity condition and hence, satisfying material stability for the entire range of deformations and electric fields, are used for the first time in the context of shell theory.