Kontaktanalyse dünnwandiger Strukturen bei großen Deformationen

Abstract

Die vorliegende Arbeit beschäftigt sich mit der numerischen Simulation von Kontaktproblemen dünnwandiger Strukturen unter Verwendung der Methode der Finiten Elemente. Dazu wird eine mortar-basierte Kontaktformulierung vorgestellt und mit geeigneten räumlichen und zeitlichen Diskretisierungsstrategien verknüpft.

Aufbauend auf die 7-Parameter-Schalenformulierung von Büchter und Ramm (1992) wird im Hinblick auf eine sinnvolle Kopplung mit der elementunabhängigen Kontaktbeschreibung ein trilineares oberflächenorientiertes hybrides Volumen-Schalen-Element hergeleitet. Ergänzend wird auf der Basis des Prinzips von Hu-Washizu ein trilineares geometrisch nichtlineares Volumenelement ausgearbeitet. Numerische Untersuchungen dokumentieren die Leistungsfähigkeit beider FE-Formulierungen.

Für die zeitliche Diskretisierung werden zwei implizite Zeitintegrationsalgorithmen eingesetzt. Neben dem bestehenden "Generalized-alpha"-Verfahren findet vor allem die "Generalized-Energy-Momentum-Method" Anwendung. Diese erweist sich in allen durchgeführten numerischen Analysen als unbedingt stabil.

Hauptbestandteil dieser Arbeit ist die Erweiterung der in Hüber und Wohlmuth (2005) vorgestellten Mortar-Kontaktformulierung auf den geometrisch nichtlinearen Fall. Durch die Einführung von kontinuierlich approximierten Lagrange-Multiplikatoren, die physikalisch den Kontaktdruck repräsentieren, wird die einzuhaltende Inpenetrabilitätsbedingung in einem schwachen integralen Sinne formuliert. Die Wahl von dualen Ansatzfunktionen (Wohlmuth (2000)) zur Interpolation der Lagrange-Multiplikatoren ermöglicht die knotenweise Entkopplung der zu erfüllenden geometrischen Randbedingungen. In Kombination mit einer Aktiven-Mengen-Strategie entsteht ein Algorithmus, der die Elimination der diskreten Knotenwerte der Lagrange-Multiplikatoren erlaubt. Diese lassen sich in einer Nachlaufrechnung variationell konsistent in Abhöngigkeit der Verschiebungen berechnen. Der resultierende Kontaktalgorithmus vereint, im Gegensatz zu vielen anderen Formulierungen, zwei wesentliche Vorteile: Lediglich die diskreten Knotenverschiebungen treten als primäre Unbekannte auf, wodurch die Größe des abschließend zu lösenden Gleichungssystems konstant bleibt; es sind keinerlei benutzerdefinierte Parameter, wie beispielsweise ein Penalty-Parameter, notwendig.

Detaillierte numerische Untersuchungen dynamischer Kontaktprobleme verdeutlichen die Notwendigkeit zusätzlicher algorithmischer Energieerhaltungsstrategien. Die von Laursen und Love (2002) vorgestellte "Velocity-Update"-Methode zeichnet sich dadurch aus, dass sie die exakte Energieerhaltung bei gleichzeitiger Erfüllung der geometrischen Nichtdurchdringungsbedingung garantiert. Sie wird entsprechend der vorgestellten Kontaktformulierung überarbeitet und für die Kombination mit der "Generalized-Energy-Momentum-Method" verallgemeinert.

Anhand von numerischen Beispielen wird die Leistungsfähigkeit der vorgestellten Lösungsstragegie analysiert und bewertet.

The present thesis is concerned with the numerical simulation of contact problems of thin-walled structures using the finite element method. A mortar-based contact formulation is presented and combined with suitable strategies for the discretization in space and time.

In view of a useful coupling with the element independent contact description, a trilinear surface oriented hybrid shell element is derived on the basis of the 7-parameter shell model by Büchter and Ramm (1992). Additionally, a trilinear geometric nonlinear brick element based on the principle of Hu-Washizu is devised. Numerical tests demonstrate the performance of both element formulations.

For the discretization in time, two implicit time integration algorithms are used. In addition to the existing "Generalized-alpha"-Method especially the "Generalized-Energy-Momentum-Method" is applied. The latter is proven to be unconditionally stable in all performed numerical analyses.

The essential part of this thesis is the extension of the mortar contact formulation presented by Hüeber and Wohlmuth (2005) to the geometrically nonlinear regime. Introducing continuously approximated Lagrange Multipliers, physically representing the contact pressure, the geometric impenetrability condition is formulated in a weak, integral sense. Using dual shape functions (Wohlmuth (2000)) for the interpolation of the Lagrange Multipliers allows for a nodal decoupling of the geometric constraints. The combination with an active set strategy results in an algorithm which allows for elimination of the discrete nodal values of the Lagrange Multipliers. They can be easily recovered from the displacements in a variational consistent way. In contrast to many other formulations, the resulting contact algorithm combines two main advantages: Only the discrete nodal displacements appear as primal unknowns, thereby the size of the system of equations to be solved remains constant; there is no need for any user defined paramters like a penalty parameter.

Detailed numerical analyses of dynamic contact problems illustrate the necessity of additional, algorithmic energy-conserving strategies. The "Velocity-Update"-method by Laursen and Love (2002) is characterized by the fact that it guarantees the exact conservation of energy while simultaneously satisfying the geometric impenetrability condition. This method is revised according to the presented contact formulation and generalized for combination with the "Generalized-Energy-Momentum-Method".

Numerical examples are investigated to analyze and judge the effectiveness of the proposed solution strategy.