Materially and geometrically nonlinear analysis of laminated anisotropic plates by p-version of FEM

A p-version finite element model based on degenerate shell element is proposed for the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with moderately large deflections and small rotations being accounted for in the sense of von Karman hypothesis. The material model is based on the Huber–Mises yield criterion and Prandtl–Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized for anisotropic materials by introducing the parameters of anisotropy. The model is also based on the equivalent-single layer laminate theory. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. Gauss–Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed p-version finite element model is demonstrated through several comparative points of view in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic zone.