Element Stiffness Matrix and Mass Matrix for Planar Triangular Shell Elements

Planar triangular shell elements are fundamental unit types in finite element analysis, particularly well-suited for modeling thin-plate structures. These elements combine characteristics of plane stress elements and plate bending elements, enabling simultaneous consideration of in-plane and out-of-plane deformations. In finite element analysis, element stiffness and mass matrices form the foundation for constructing global structural equations. For planar triangular shell elements, the stiffness matrix typically comprises two components: in-plane stiffness and bending stiffness. The in-plane stiffness component resembles plane stress problems, while the bending stiffness component relates to thin plate bending behavior. Mass matrix formulation commonly adopts two approaches: consistent mass matrix and lumped mass matrix. The consistent mass matrix, obtained through shape function integration, more accurately represents mass distribution but requires greater computational effort. The lumped mass matrix concentrates element mass at nodes, offering simpler computation with slightly reduced accuracy. When implementing these matrices in MATLAB, developers typically follow these steps: First define geometric parameters and material properties using structure arrays or separate variables. Then compute shape functions and their derivatives using barycentric coordinates or area coordinates. Subsequently construct stiffness and mass matrices through numerical integration techniques like Gaussian quadrature, where integration points and weights are calculated using standard quadrature rules. For shell elements, special attention must be paid to coordinate transformations between local and global coordinate systems, typically implemented through transformation matrices using direction cosine calculations. These elements find extensive engineering applications in analyzing thin-walled structures such as aircraft wings and ship hulls. Through appropriate selection of element types and parameters, accurate structural response predictions can be achieved.