A new finite element scheme for instability analysis of thin shells

This paper describes a new finite element scheme for the analysis of instability phenomena of arbitrary thin shells. A computationally efficient procedure is proposed for calculating the non-linear stiffness and tangential stiffness matrices for a doubly-curved quadrilateral element defined by co-ordinate lines. The essential feature is the explicit addition of the non-linear terms into the rigid-body motion of the element. Thus the non-linear and tangential element stiffness matrices can easily be generated by transforming the generalized element stiffness matrix for linear analysis, and the non-linear terms of these matrices are separated into a number of component terms multiplied by the rigid-body rotations. These component terms can be stored permanently and used to calculate efficiently the non-linear and tangential stiffness matrices at each iteration. Illustrative examples are presented which confirm the validity of the present approach in the analysis of instability phenomena of thin plates and shells.