On an isoparametric finite-element for composite laminates with finite rotations

For composite laminates consisting of an arbitrary number of orthotropic laminae first a finite-rotation theory is presented as basis of isoparametric finite-element formulations. The derivation is achieved by a Reissner-Mindlin type kinematic assumption which allows a constant shear deformation distribution across the thickness. The constitutive equations are presented in a general form such that orthotropic material behaviour with material axes varying arbitrarily across the thickness may easily be considered in numerical implementation, also when using curvilinear coordinates. Special attention is taken to predict the force distribution in the deformed shell structure. This theory is then transformed into a four-node isoparametric assumed-strain finite element. Unlike in the degeneration approach, interpolation polynomials are introduced directly for rotation variables determining the deformed position of the unit normal vector. The capability of the finite element developed to deal with strongly nonlinear situations is demonstrated by many examples. Also numerical results are presented permitting a systematical comparison of classical and isoparametric approaches concerning the numerical efficiency.