Adaptive finite elements for a contact problem in elastoplasticity with Lagrange techniques

Abstract

Das Thema dieser Dissertation ist die Herleitung und numerische Analyse von finiten Elementen für ein Problem in der Elastoplastizität mit Kontaktbedingungen. Die hergeleiteten finite Elemente Verfahren basieren auf einer Formulierung als Sattelpunktproblem und der Nutzung von Polynomen höherer Ordnung. Die Analyse der vorgestellten Verfahren beginnt mit dem Zeigen der Wohldefiniertheit und der Konvergenz. Im nächsten Schritt werden a priori Abschätzungen der Konvergenzraten gezeigt. Weiterhin führt die Einführung von Lagrange Multiplikatoren zu einem einheitlichen Ansatz zur a posteriori Abschätzung des Diskretisierungfehlers unter der Verwendung von Elementen höherer Ordnung. Zusätzlich ermöglicht es der Zugang über Lagrange Multiplikatoren die Äquivalenz der Diskretisierungsfehler in den Spannungen und in den Energien für finite Elemente niederer Ordnung zu zeigen, was insbesondere neu für Viereckselemente ist. Diese Äquivalenz wiederum erlaubt nun den Beweis der Konvergenz von adaptiven finiten Elementen niederer Ordnung. Für Dreieckselemente wird sogar die optimale Konvergenz bewiesen. Die theoretischen Erkenntnisse werden durch numerische Experimente bestätigt.

The topic of this thesis is the derivation and analysis of some finite element schemes for a contact problem in elastoplasticity. These schemes are based on the formulation of the models as saddle point problems and use finite element spaces of arbitrary polynomial degrees. In this thesis, these new approaches with higher-order finite elements are shown to be well defined and convergent. Moreover, some a~priori estimates on the rates of convergences are proven. The use of Lagrange multipliers in the saddle point formulation yields a coherent approach to reliable a~posteriori error estimates for the proposed higher-order schemes. Additionally, the Lagrange multipliers are used to show the equivalence of the errors of the stresses and the energies, for low order finite elements using triangular or quadrilateral cells. For the first time, this allows for a proof of convergence for quadrilateral-based adaptive finite elements. Furthermore, the approach based on triangular cells is shown to be of optimal convergence. The theoretical findings are confirmed by numerical experiments.