In this article, a simple and accurate quadrilateral thin flat layered shell element formulation for the nonlinear analysis of reinforced concrete (RC) wall systems under static and cycling loads is presented. The 4 node shell element, with 6 degree of freedom (DOF) per node (3 displacements and 3 rotations) is created by superposing the quadrilateral layered membrane element with drilling degrees of freedom (12 DOF, 2 displacement and 1 rotation per node) developed by Rojas et al. (Eng Struct 124:521-538, 2016), and the Discrete Kirchhoff Quadrilateral Element (12 DOF, 1 displacement and 2 rotations) formulated by Batoz and Tahar (Int J Numer Methods Eng 18(11): 1982), to model the inplane and the out of plane bending behavior of the shell element, respectively. In addition, to model the complex behavior and coupling of the axial, flexural and shear behavior, observed in complex RC wall structures, the transversal section of the shell element consists of a layered system of fully bonded, smeared steel reinforcement and smeared orthotropic concrete material with the rotating angle formulation. The formulation used a tangent stiffness matrix approach, which include the coupling of membrane and bending effects. For verification, the shell element formulation is used to model a set of experimental results for T-shaped RC walls that are available in the literature. The proposed element is robust, simple to implement, and it can predict the global results (load vs. displacement and maximum capacity) and also the local behavior (vertical strain at the base level along the web and the flange) observed in RC wall structures.