Wie funktioniert eine Kontaktanalyse in der FEM-Simulation?

Michael ·
Zwei präzisionsgefertigte Stahlbauteile mit sichtbarer Spannungskonzentration an der Kontaktzone auf einer Werkbank.

Eine Kontaktanalyse in der FEM-Simulation beschreibt das mechanische Verhalten an den Berührungsflächen zweier oder mehrerer Bauteile, wenn diese unter Last miteinander in Wechselwirkung treten. Die Finite-Elemente-Methode löst dabei die Gleichgewichtsbedingungen an der Kontaktzone iterativ, weil Kontakt eine geometrische und physikalische Nichtlinearität einführt. Die folgenden Abschnitte beantworten die häufigsten Fragen, die Projektingenieure bei der Vorbereitung und Auswertung einer Kontaktanalyse stellen.

Welche Kontaktarten werden in der FEM-Simulation unterschieden?

In der FEM-Simulation werden typischerweise vier grundlegende Kontaktarten unterschieden: gebundener Kontakt (bonded), reibungsfreier Kontakt (frictionless), Kontakt mit Reibung (frictional) sowie Kontakt ohne Trennung (no separation). Jede dieser Varianten beschreibt ein anderes physikalisches Verhalten an der Fügestelle und hat direkten Einfluss auf Rechenaufwand und Konvergenzverhalten.

Beim gebundenen Kontakt werden die Kontaktflächen starr miteinander verbunden, als wären sie verschweißt oder verklebt. Relativbewegungen sind nicht möglich, was die Analyse linear hält und den Rechenaufwand gering hält. Dieser Kontakttyp eignet sich für Pressverbindungen oder Klebstellen, bei denen keine Relativverschiebung zu erwarten ist.

Reibungsfreier Kontakt erlaubt Gleiten entlang der Kontaktfläche ohne Widerstand, verhindert aber das Durchdringen der Körper ineinander. Kontakt mit Reibung geht einen Schritt weiter und berücksichtigt den Reibungskoeffizienten nach Coulomb. Sobald die Tangentialspannung den Grenzwert überschreitet, wechselt der Kontaktzustand von Haften zu Gleiten, was die Nichtlinearität erheblich erhöht.

Wie wird Kontakt in der FEM mathematisch formuliert?

Kontakt in der FEM wird mathematisch über Ungleichungsbedingungen formuliert, die sicherstellen, dass sich zwei Körper nicht gegenseitig durchdringen. Die zwei verbreitetsten Methoden sind die Penalty-Methode und die Lagrange-Multiplikator-Methode, wobei viele Solver hybride Ansätze wie die Augmented-Lagrange-Methode verwenden.

Bei der Penalty-Methode wird eine virtuelle Federsteifigkeit zwischen den Kontaktpartnern eingeführt. Dringt ein Knoten in die Gegenfläche ein, erzeugt diese Federsteifigkeit eine rücktreibende Kontaktkraft proportional zur Eindringtiefe. Der Vorteil liegt in der numerischen Robustheit, der Nachteil in einer kleinen, kontrollierten Restpenetration, die bei zu weicher Penalty-Steifigkeit das Ergebnis verfälscht.

Die Lagrange-Multiplikator-Methode erzwingt die Kontaktbedingung exakt ohne Penetration, indem sie zusätzliche Freiheitsgrade in das Gleichungssystem einführt. Das Gleichungssystem wird dadurch größer und potenziell schlechter konditioniert, liefert aber physikalisch präzisere Kontaktkräfte. Die Augmented-Lagrange-Methode kombiniert beide Ansätze und ist in modernen Solvern wie ANSYS Mechanical der Standard für nichtlineare Kontaktanalysen.

Was ist der Unterschied zwischen symmetrischem und asymmetrischem Kontakt?

Bei asymmetrischem Kontakt wird eine Fläche als Kontaktfläche (contact) und die andere als Zielfläche (target) definiert. Nur die Knoten der Kontaktfläche werden auf Penetration in die Zielfläche geprüft. Bei symmetrischem Kontakt übernehmen beide Flächen gleichzeitig die Rolle von Kontakt- und Zielfläche, was eine beidseitige Penetrationsprüfung bewirkt.

Asymmetrischer Kontakt ist recheneffizienter und in den meisten Fällen ausreichend, solange die Netzfeinheit der Kontaktfläche geringer ist als die der Zielfläche und die Geometrien klar definiert sind. Die Zuweisung eines weichen Körpers als Kontaktfläche und eines steifen Körpers als Zielfläche ist dabei eine bewährte Faustregel.

Symmetrischer Kontakt wird empfohlen, wenn beide Körper ähnliche Steifigkeiten und vergleichbare Netzauflösungen aufweisen oder wenn unklar ist, welche Seite stärker deformiert. Der erhöhte Rechenaufwand ist in solchen Fällen gerechtfertigt, weil asymmetrischer Kontakt sonst Penetrationen auf der nicht geprüften Seite übersehen kann.

Warum konvergiert eine FEM-Kontaktanalyse nicht?

Eine FEM-Kontaktanalyse konvergiert nicht, wenn das nichtlineare Gleichungssystem innerhalb der erlaubten Iterationsschritte keinen Gleichgewichtszustand findet. Die häufigsten Ursachen sind sprunghafte Kontaktzustandswechsel, eine zu steife oder zu weiche Penalty-Steifigkeit, unzureichende Lastinkrementgrößen oder geometrisch schlecht definierte Kontaktflächen.

Sprunghafte Statuswechsel entstehen, wenn Knoten innerhalb eines Iterationsschritts mehrfach zwischen Haften, Gleiten und Abheben wechseln. Der Solver findet keinen stabilen Zwischenzustand und bricht ab. Abhilfe schaffen kleinere Lastinkremente, eine Stabilisierungsoption oder die Überprüfung, ob der Kontakt physikalisch sinnvoll modelliert ist.

Eine falsch gewählte Kontaktsteifigkeit ist eine weitere häufige Ursache. Ist sie zu hoch, erzeugt der Solver große Ungleichgewichtskräfte bei minimaler Penetration, was Konvergenzprobleme verursacht. Ist sie zu niedrig, penetrieren die Körper zu stark, und das Ergebnis verliert physikalische Aussagekraft. Die automatische Steifigkeitsanpassung moderner Solver löst dieses Problem in vielen Fällen, erfordert aber eine Plausibilitätsprüfung der Kontaktergebnisse.

Wann sollte eine nichtlineare Kontaktanalyse eingesetzt werden?

Eine nichtlineare Kontaktanalyse ist notwendig, wenn sich der Kontaktzustand während der Belastung ändert, also wenn Flächen abheben, gleiten oder neu in Kontakt treten. Sobald Reibung, große Verformungen oder wechselnde Kontaktflächen eine Rolle spielen, reicht ein linearer gebundener Kontakt nicht mehr aus.

Typische Anwendungsfälle sind Schraubenverbindungen unter Betriebslast, Lager- und Gelenkverbindungen, Dichtflächen unter Druckbelastung sowie Pressverbindungen mit Übermaß. In all diesen Situationen hängt die Spannungsverteilung direkt davon ab, welcher Bereich der Kontaktfläche tatsächlich trägt und welcher abhebt.

Für Bauteile, bei denen Kontaktflächen unter Last vollständig geschlossen bleiben und keine Relativbewegung auftritt, ist ein gebundener linearer Kontakt oft ausreichend und deutlich effizienter. Die Entscheidung zwischen linearem und nichtlinearem Kontakt sollte auf einer physikalischen Einschätzung des erwarteten Verhaltens basieren, nicht allein auf numerischer Vorsicht. Mehr zu den Grundlagen der FEM-Berechnung und Strukturanalyse zeigt ein Überblick über das Leistungsspektrum moderner Ingenieurdienstleister.

Wie beeinflusst die Netzqualität das Ergebnis einer Kontaktanalyse?

Die Netzqualität hat direkten Einfluss auf die Genauigkeit einer Kontaktanalyse, weil Kontaktspannungen und Kontaktkräfte an den Knoten der Kontaktfläche berechnet werden. Ein zu grobes Netz in der Kontaktzone führt zu einer Mittelung der Spannungen über zu große Bereiche, was lokale Spannungsspitzen unterschätzt und die Kontaktdruckverteilung verfälscht.

Besonders kritisch ist die Netzkompatibilität zwischen den Kontaktpartnern. Wenn die Elementgröße der Kontaktfläche deutlich größer ist als die der Zielfläche, können Knoten der Zielfläche zwischen den Kontaktknoten liegen und werden bei asymmetrischem Kontakt nicht auf Penetration geprüft. Eine vergleichbare Netzfeinheit auf beiden Seiten der Kontaktzone ist deshalb eine Grundvoraussetzung für verlässliche Ergebnisse.

Elementform und Elementtyp spielen ebenfalls eine Rolle. Quadratische Elemente zweiter Ordnung bilden den Spannungsverlauf in der Kontaktzone besser ab als lineare Elemente erster Ordnung, weil sie innerhalb eines Elements einen nichtlinearen Verschiebungsverlauf zulassen. Bei stark gekrümmten Kontaktflächen empfiehlt sich zudem eine Netzverfeinerung, die der Geometrie folgt, um Diskretisierungsfehler an der Krümmung zu minimieren.

FEM-Kontaktanalysen professionell durchführen lassen

Kontaktanalysen gehören zu den anspruchsvollsten Aufgaben in der FEM-Simulation, weil sie mehrere Nichtlinearitäten gleichzeitig erfordern und das Ergebnis stark von Modellierungsentscheidungen abhängt. CE-SYS Engineering unterstützt technische Projektleiter bei der Durchführung solcher Analysen, von der Auswahl der richtigen Kontaktformulierung bis zur Interpretation der Ergebnisse. Das Leistungsangebot umfasst:

  • Nichtlineare FEM-Berechnungen mit reibungsbehaftetem und statuswechselndem Kontakt
  • Spannungsnachweise und Verformungsanalysen unter statischer und dynamischer Belastung
  • Netzerstellung und Netzoptimierung speziell für Kontaktzonen
  • Konvergenzanalyse und Fehlerbehebung bei nicht konvergierenden Modellen
  • Nachweisführung nach europäischen und internationalen Regelwerken

Für Rückfragen zu laufenden Projekten oder zur Einschätzung, ob eine nichtlineare Kontaktanalyse für Ihre Anwendung notwendig ist, steht das Ingenieurteam direkt zur Verfügung. Jetzt Kontakt aufnehmen und eine unverbindliche Ersteinschätzung erhalten.

Ähnliche Artikel