Modelling and numerical simulation of contact and lubrication

Abstract

This thesis considers the simulation of a steel forming process including the lubricant between tool and workpiece. The model equations, that describe deformation, contact and hydrodynamic flow, are derived from fundamental physical laws. For the contact-simulation the two models of Signorini- and z-contact are introduced. The 3-dimensional Stokes problem for a thin fluid film is condensed to the 2-dimensional SubStokes model (for velocity and pressure) and further reduced to the Reynolds model, which is an equation for the pressure only. All fluid models are extended to variational inequations to consider cavitation. Basic mathematical concepts for the numerical treatment with the finite element method and the corresponding analysis are presented. Error estimation is done separate for the modelling error and discretisation error. An algorithm for model adaptivity is given and model error estimates are derived out of the physical models for z- and Signorini-contact aswell as SubStokes- and Reynolds-fluids. Error estimates for the discretisation error of elliptic problems - deformation and Reynolds-flow - are presented. The SubStokes problem is stabilized to apply linear finite elements. To estimate the discretisation error of the SubStokes inequation a Lagrangian multiplier is introduced. An algorithm for a simultanous refinement of model and grid is given. Finally the presented methods and estimates are validated by the application to various prototype examples. The results allow to compare the different models (z- and Signorini-contact, Reynolds- and SubStokes-fluid).

Diese Arbeit behandelt die Simulation eines Stahlformungsprozesses einschließlich des Schmieröls zwischen Werkzeug und Werkstück. Die Modellgleichungen, die Deformation, Kontakt und Hydrodynamik zugrunde liegen, werden aus physikalischen Gesetzen abgeleitet. Für die Kontakt-Simulation werden die beiden Modelle des Signorini- und z-Kontaktes vorgestellt. Das 3-dimensionale Stokes-Problem für eine dünne Schicht wird zum 2-dimensionalen SubStokes-Modell verdichtet und weiter vereinfacht zum Reynolds-Modell, das nur noch eine Gleichung für den Druck ist. Alle Fluid-Modelle werden zu Variationsungleichungen erweitert, um Kavitation abzubilden. Mathematische Grundlagen zur numerischen Simulation mittels Finiter Elemente und die zugehörige Analyse werden bereitgestellt. Modell- und Diskretisierungsfehler werden separat abgeschätzt. Es wird ein Algorithmus zur Modelladaptivität aufgestellt und aus den physikalischen Modellen für z- und Signorini-Kontakt sowie für SubStokes- und Reynolds-Fluide Modellfehlerschätzer hergeleitet. Fehlerschätzer für den Diskretisierungsfehler elliptischer Probleme - Deformation und Reynolds-Strömung - werden präsentiert. Das SubStokes-Problem wird stabilisiert, um lineare Finite Elemente darauf anzuwenden. Um den Diskretisierungsfehler der SubStokes-Ungleichung abzuschätzen, wird ein Lagrange-Multiplikator eingefügt. Ein Algorithmus zur simultanen Modell- und Gitter-Verfeinerung wird aufgestellt. Schließlich werden die vorgestellten Methoden und Schätzer durch die Anwendung auf vielfältige Beispiele validiert. Die Ergebnisse gewähren einen Vergleich der unterschiedlichen Modelle von Kontakt und Fluid.