This study introduces a static analysis of laminated sandwich spherical shell panels applying a novel two-dimensional facet-shell finite element. This element combines the characteristics of membrane elements and plate bending element behaviors and draws upon the zigzag theory. Further development includes the discrete Kirchhoff quadrilateral element refined by the co-author for composite plate analysis. The incorporation of a matrix of transformation is crucial for transforming local actions and displacements into global actions and displacements. This element, notable for its inclusion of two fictitious degrees of freedom, supports nine local and nine global DOF at each node, effectively bypassing issues associated with ill-conditioned stiffness matrices. The efficacy of this element is verified through rigorous comparisons with existing analytical, 3D elasticity solutions from the literature. It is tested under various boundary conditions, material types, and geometrical configurations, demonstrating robust and accurate performance in evaluating thick and moderately thick spherical shell panels.