GEOMETRICALLY NONLINEAR ANALYSIS OF THIN ARBITRARY SHELLS USING DISCRETE-KIRCHHOFF CURVED TRIANGULAR ELEMENTS (FINITE).

Abstract

The research work presented here deals with the problems of geometrically nonlinear analysis of thin shell structures. The specific objective was to develop geometrically nonlinear formulations, using Discrete-Kirchhoff Curved Triangular (DKCT) thin shell elements. The DKCT elements, formulated in the natural curvilinear coordinates, based on arbitrary deep shell theory and representing explicit rigid body modes, were successfully applied to linear elastic analysis of composite shells in an earlier research work. A detailed discussion on the developments of classical linear and nonlinear shell theories and the Finite Element applications to linear and nonlinear analysis of shells has been presented. The difficulties of developing converging shell elements due to Kirchhoff's hypothesis have been discussed. The importance of formulating shell elements based on deep shell theory has also been pointed out. The development of shell elements based on Discrete-Kirchhoff's theory has been discussed. The development of a simple 3-noded curved triangular thin shell element with 27 degrees-of-freedom in the tangent and normal displacements and their first-order derivatives, formulated in the natural curvilinear coordinates and based on arbitrary deep shell theory, has been described. This DKCT element has been used to develop geometrically nonlinear formulation for the nonlinear analysis of thin shells. A detailed derivation of the geometrically nonlinear (GNL) formulation, using the DKCT element based on the Total Lagrangian approach and the principles of virtual work has been presented. The techniques of solving the nonlinear equilibrium equations, using the incremental methods has been described. This includes the derivation of the Tangent Stiffness matrix. Various Newton-Raphson solution algorithms and the associated convergence criteria have been discussed in detail. Difficulties of tracing the post buckling behavior using these algorithms and hence the necessity of using alternative techniques have been mentioned. A detailed numerical evaluation of the GNL formulation has been carried out by solving a number of standard problems in the linear buckling and GNL analysis. The results compare well with the standard solutions in linear buckling cases and are in general satisfactory for the GNL analysis in the region of large displacements and small rotations. It is concluded that this simple and economical element will be an ideal choice for the expensive nonlinear analysis of shells. However, it is suggested that the element formulation should include large rotations for the element to perform accurately in the region of large rotations.