Subject:  A New Family of Mixed Finite Elements for Elasticity 
treated by:  Dipl.Ing. Astrid Sinwel 
Supervisor:  Prof. Dr. Joachim Schöberl 
External Referee:  Prof. Dr. Ragnar Winther (University of Oslo) 
In our thesis, we introduce a new finite element method to discretize the equations of elasticity. It is analyzed thoroughly both in the infinitedimensional and in the discrete setting, where finite element schemes of arbitrary order are presented. As the main result of this work, we prove that the new method is lockingfree with respect to volume and shear locking, i.e. that it is applicable for both nearly incompressible materials and the discretization of the thin structures using flat elements. To date, several wellknown methods for the discretization and subsequent solution of the equationsof elasticity have been introduced: the primal method using continuous finite element functions for the displacement, the mixed method due to Hellinger and Reissner, and mixed methods with weak symmetry. However, each of these methods has its drawbacks, which motivates the need for yet another formulation. The primal method suffers from volume and shear locking phenomena. Mixed methods can be shown to be stable with respect to volume locking; however, the construction of finite elements can only be done at high polynomial orders, and therefore high computational cost. Mixed methods with weak symmetry are easier to construct, but the exact symmetry of the stress tensor is lost. The new method lies in between the primal method and the mixed HellingerReissner method. The vectorvalued displacement function is chosen in the space H(curl), which ensures continuity of the tangential component across interfaces. This choice implies to search for the stresses in the newly introduced space H(div div) consisting of symmetric tensor valued functions, whose divergence again allows for a distributional divergence lying in H^{1}. It is shown that such tensor fields have their normalnormal component continuous. The resulting mixed formulation is referred to as TangentialDisplacementNormalNormalStress (TDNNS) formulation. Its stability in the infinitedimensional setting is analyzed. To discretize the TDNNS formulation, a stable pair of finite elements is introduced. For the displacement space, Nédélec elements as standard choice for an H(curl) conforming discretization are used. For the stress space, a new family of symmetric tensorvalued finite elements of arbitrary order is constructed. Stability and optimal order approximation properties of the resulting finite element method are provided. A major drawback of mixed methods is the indefiniteness of the resulting system matrix. In order to obtain a positive definite matrix, we apply hybridization. The normalnormal continuity of the stresses is torn, and reenforced by Lagrangian multipliers resembling the normal displacement on interfaces. After local elimination of the stresses, one obtains a positive definite matrix. In this setup, stability with respect to volume locking is provided. An additive Schwarz block preconditioner is shown to be an optimal choice. Finally, the TDNNS method is used on thin structures, where flat prismatic elements are used. Stability and error estimates which are independent of the aspect ratio of the elements are derived. All results are confirmed by computational examples. 
Download

back to PhD theses 