A scaled boundary shell element formulation using Neumann expansion

This paper proposes a new shell element formulation using the scaled boundary finite element (SBFE) method. A shell element is treated as a three-dimensional continuum. Its bottom surface is approximated with a quadrilateral spectral element and the shell geometry is represented through normal scaling of the bottom surface. Neumann expansion is applied to approximate the inversions of the matrix polynomials of the thickness coordinate ξ, including the Jacobian matrix and the coefficient of the second-order term in the SBFE equation. This permits the solution along the thickness to be expressed as a matrix exponential function whose exponent is a high-order matrix polynomial of ξ. After introducing the boundary conditions on the top and bottom surfaces and evaluating the resulting matrix exponential via Padé expansion, we derive the element stiffness and mass matrices. Poisson thickness locking is avoided fundamentally. Numerical examples demonstrate the applicability and efficiency of the formulation.