A 4-node co-rotational quadrilateral composite shell element is presented. The local coordinate system of the element is a co-rotational framework defined by the two bisectors of the diagonal vectors generated from the four corner nodes and their cross product. Thus, the element rigid-body rotations are excluded in calculating the local nodal variables from the global nodal variables. Compared with other existing co-rotational finite-element formulations, the present element has two features: (i) The two smallest components of the mid-surface normal vector at each node are defined as the rotational variables, leading to the desired additive property for all nodal variables in a nonlinear incremental solution procedure; (ii) both element tangent stiffness matrices in the local and global coordinate systems are symmetric owing to the commutativity of the nodal variables in calculating the second derivatives of strain energy with respect to the local nodal variables and, through chain differentiation with respect to the global nodal variables. In the modeling of composite structures, the first-order shear deformable laminated plate theory is adopted in the local element formulation, where both the thickness deformation and the normal stress in the direction of the shell thickness are ignored, and an assumed strain method is employed to alleviate the membrane and shear locking phenomena. Several examples involving composite plates and shells with large displacements and large rotations are presented to testify to the reliability and convergence of the present formulation