Theory and numerics of thin elastic shells with finite rotations

Summary A bending theory for thin shells undergoing finite rotations is presented, and its associated finite element model is described. The kinematic assumption is based on a Reissner-Mindlin theory. The starting point for the derivation of the strain measures is the polar decomposition of the material deformation gradient. The work-conjugate stress resultants and stress couples are integrals of the Biot stress tensor. This tensor is invariant with respect to rigid body motions and therefore appropriate for the formulation of constitutive equations. The rotations are described by using Eulerian angles. The finite element discretization of arbitrary shells is performed using isoparametric elements. The advantage of the proposed shell formulation and its numerical model is shown by application to different non-linear plate and shell problems. Finite rotations can be calculated within one load increment. Thus the step size of the load increment is only imited by the local convergence behaviour of Newton's method or the appearance of stability phenomena.

---

Theorie und numerische Analysis dünner elastischer Schalen mit endlichen Drehungen

Übersicht Es wird eine Biegetheorie elastischer Schalen mit endlichen Drehungen dargestellt und die numerische Behandlung mit der Methode der finiten Elemente beschrieben. Dazu müssen die Verzerrungen und die schwache Form des Gleichgewichts angegeben werden. Es wird eine Reissner-Mindlin-Kinematik zugrunde gelegt. Die Verwendung des Greenschen Verzerrungstensors führt im Schalenraum bei endlichen Drehungen auf sehr komplizierte Ausdrücke. Daher wird von der polaren Zerlegung des materiallen Deformationsgradienten ausgegangen. In diesem Fall sind die Kräfte und Momente Spannungsresultierende des Biot-Spannungstensors, der invariant gegenüber Starrkörperbewegungen und somit für ein Stoffgesetz geeignet ist. Als Drehkinematen werden Eulerwinkel verwandt. Um allgemeine Schalengeometrien beschreiben zu können, wird eine isoparametrische Elementformulierung gewählt. Endliche Drehungen können in den berechneten Platten- und Schalenproblemen in einem Lastschritt aufgebracht werden. Die Lastschrittweite wird nur durch das Lösungsverfahren und das Auftreten von Stabilitätspunkten beschränkt.

---

---

Dedicated to the memory of Dr. Dieter Withum who died on June 25, 1979. At the time of his death, Dr. Withum was Professor of Engineering Mechanics at the University of the German Federal Armed Forces in Munich. He would have been 60 years old on September 29, 1988