Numerische Schwingungsanalyse

Die von uns entwickelten Methoden zur numerischen Simulation gehen robust und effizient mit der Unstetigkeit der Kontaktvorgänge um und reduzieren den Rechenaufwand im Vergleich zu gebräuchlichen Finite-Elemente-Tools um mehrere Größenordnungen.

Um die Dynamik von Bauteilen und technischen Systemen genau zu berechnen, müssen die Vorgänge in den Kontakten räumlich und zeitlich hinreichend fein aufgelöst werden. Die gebräuchlichen Simulationsmethoden zeigen Artefakte wie künstliche Oszillationen, Instabilität und/oder numerische Dämpfung. Die von uns entwickelten numerischen Methoden ermöglichen die Simulation hoch-präziser Modelle. So haben wir voll-gekoppelte Methoden zur Strömung-Struktur-Wechselwirkung entwickelt. Neben den nichtlinearen Vorgängen in den Kontakten wird das inhärent nichtlineare Verhalten der Strömung abgebildet, um die aerodynamische Dämpfung und Anregung von Schaufeln in Triebwerken zu berechnen.

00:01

Harmonische Balance

Mit der Harmonischen Balance lassen sich periodische Lösungen nichtlinearer gewöhnlicher Differential- und Differential-Algebraischer Gleichungen berechnen. Als Ansatz wird eine abgebrochene Fourier-Reihe verwendet und es wird ausgenutzt, dass in vielen Fällen nur wenige Fourier-Terme nötig sind, um die Schwingungen sehr genaue zu beschreiben. Auch muss das bei der typischerweise schwachen Dämpfung langwierige Einschwingen nicht simuliert werden. Wir haben ein einfaches Tool veröffentlicht, das auch Möglichkeiten zur numerischen Pfadverfolgung bietet, siehe NLvib-Unterseite.

Die folgende Animation zeigt wie die Harmonische Balance mit zunehmender Fourier-Ordnung (H) gegen die Zeitbereichslösung (schwarze Linie) konvergiert (gerechnet mit unseren Tools an einem Knoten).

00:20

In der Animation der räumlichen Relativverschiebung der Kontaktfläche kann man die drei lokal wechselnden Kontaktzustände (Abheben, Haften und Gleiten) beobachten (gerechnet mit unseren Tools). Lila markiert dabei keine Relativverschiebung, rot bedeutet maximale Relativverschiebung.

00:04

Simulation im Zeitbereich

Neben periodischen Schwingungen gewinnen nicht-periodische und transiente Vorgänge zunehmend an Bedeutung. Selbst im periodischen Fall können stoßartige Vorgänge zu einer langsamen Konvergenz der Harmonischen Balance führen. Daher entwickeln wir Methoden zur numerischen Zeit-Schritt-Integration, welche speziell für reibungsfreien sowie reibungsbehafteten Kontakt geeignet sind. Übliche Finite-Elemente-Methoden verknüpfen eine endliche Masse mit jedem Knoten, auch mit denen an der Kontaktfläche. Wie veranschaulicht, führt dies zu künstlichen Oszillationen und großen Fehlern. Unsere Ansätze führen auf einen masselosen Rand, um solche Fehler zu vermeiden. Gleichzeitig reduziert sich damit der Rechenaufwand um mehrere Größenordnungen ohne die hervorragende Übereinstimmung mit den Messungen zu beeinträchtigen.

00:05

Parallelisierte Modell-Verfeinerung

Die numerische Pfadverfolgung wird häufig angewendet, um zu bestimmen, wie sich Schwingungen mit einem freien Parameter (z.B. der Anregungsfrequenz) entwickeln. Bei hoch-präzisen Modellen komplexer nichtlinearer Systeme kann die der Pfadverfolgung inhärente sequenzielle Natur zu einem unüberwindbaren Hindernis werden. Das von uns entwickelte, breit anwendbare Konzept ist sehr einfach und birgt dennoch das Potenzial für einen Durchbruch: Die Idee besteht darin, Pfadverfolgung und Modell-Verfeinerung auf neue und kreative Weise zu verknüpfen. Zunächst wird die numerische Pfadverfolgung auf ein Modell niedriger Genauigkeit angewendet (geringer Rechenaufwand). Anschließend wird eine Teilmenge interessanter Punkte ausgewählt, und ausgehend von diesen werden Punkte auf der Lösungskurve des Modells hoher Genauigkeit berechnet. Die iterative Korrektur kann für jeden Punkt unabhängig durchgeführt werden und ist somit prädestiniert für die parallele Berechnung. Dieses Konzept hat sich als numerisch sehr robust erwiesen und ermöglicht eine Beschleunigung von Berechnungen um mehrere Größenordnungen.

00:12
Dieses Bild zeigtJohann Groß

Johann Groß

Dr.-Ing.

Gruppenleiter „Numerische Methoden“ im Bereich Strukturmechanik

Zum Seitenanfang