Software: FEM - Tutorial - Diskretisierung - Ansatzfunktionen

↑

← →

Ansatzfunktionen

Grundlagen

• Für die jeweilige Elementform werden nur die Potential-Werte ihrer Knoten-Punkte als so genannte Primär-Ergebnisse berechnet.
• Zusätzlich muss nun eine geeignete Funktion definiert werden, welche es ermöglicht, aus diesen Knoten-Potentialen für alle Punkte innerhalb des Elements einen eindeutigen (und sinnvollen) Wert zu berechnen (Sekundär-Ergebnisse).
• Man kann sich die Knoten als Stützstellen dieser Funktion vorstellen.
• Diese Funktion definiert man für Finite Elemente auf der Basis von Ansatzfunktionen.
• Von einer "Ansatzfunktion" spricht man, wenn sie folgende Bedingungen erfüllt (Nach [1]):

1. Die Funktion ist auf dem ganzen Element definiert.
2. Jedem Knoten des Elements ist eine Funktion zugeordnet.
3. Am Knoten hat diese Funktion den Wert=1, an allen anderen Knoten des Elements hat sie den Wert=0.
4. Für jeden Punkt des Elements hat die Summe aller Funktionen den Wert=1.
5. An gemeinsamen Kanten (oder Flächen) zu benachbarten Elementen haben die Funktionen gemeinsamer Knoten an jedem Punkt den gleichen Wert.

• Üblich für den Begriff "Ansatzfunktion" sind auch:
  • Formfunktion (eng. form-function / shape-function)
  • Interpolationsfunktion (engl. interpolation-function)
  • Näherungsfunktion

Dies soll am Beispiel eines 1-dimensionalen Elements der Länge=1 für lineare Ansatzfunktionen verdeutlicht werden:

• Gegeben sind die Funktionswerte y₁ und y₂ an den Knoten des Stabes.
• Es werden zwei Ansatzfunktionen definiert

h₁(x)=1-x

h₂(x)=x

• Gewichtet mit den Knotenwerten y₁ und y₂ lautet die Elementfunktion

y(x)=y₁·h₁(x)+y₂·h₂(x)

y(x)=y₁·(1-x)+y₂·x

Analog ist die Vorgehensweise bei höherwertigen Ansatzfunktionen:

• Quadratische (=parabolische) Ansatzfunktionen erfordern in der Mitte des Stabes einen zusätzlichen Knoten.
• Es werden für den Stab drei Ansatzfunktionen entsprechend der Knotenzahl gebildet.
• Diese sind so formuliert, dass Sie obigen Bedingungen genügen.

• Kubische Ansatzfunktionen erfordern zwei Zwischenknoten auf dem Stab:

Je höherwertig die gewählte Ansatzfunktion ist, desto besser gelingt die Annäherung der interpolierten Werte an den realen Verlauf der physikalischen Größen. Dies betrifft sowohl die Spannungen (Potentiale) als auch die mechanischen Verformungen der Elemente unter Belastung. Allerdings steigt mit der erforderlichen Knotenzahl auch die Zeit für die Modellberechnung.

In der Praxis beschränkt man sich auf lineare und insbesondere quadratische Ansatzfunktionen. Zur Erreichung einer höheren Genauigkeit wählt man statt einer höheren Ansatzfunktion eine engere Vernetzung!

Für Spezialanwendungen haben sich neben den obigen Polynomansätzen weitere Ansatzfunktionen etabliert:

• Realisierung weicher Übergänge an den Elementgrenzen (ohne Knick), d.h. stetige Differenzierbarkeit auch über die Grenzen hinweg.
• trigonometrische Funktionen
• Im Element=konstant (Treppe über die Elemente) für schnell zu berechnende, riesige Netze (Strömungs- und Crash-Simulationen).
• usw. ...

Höherwertige Ansatzfunktionen erfordern auf den Elementen zusätzlich zu den Eck-Knoten noch Zwischenknoten auf den Kanten bzw. innerhalb der Elemente. Dies soll am Beispiel der quadratischen Ansatzfunktion für das Viereck-Element verdeutlicht werden:

• Der normale Polynom-Ansatz erfordert 5 zusätzliche Knoten.
• Der heute übliche Lagrang'sche Ansatz benötigt den Innenknoten nicht.
• Ein unter Belastung verformtes Viereck kann mit seinen gekrümmten Kanten bei Verwendung von Zwischenknoten besser approximiert werden:
  

Wirkung bei der Modellberechnung

• Die verwendeten Ansatzfunktionen fließen vollständig ein in das Gleichungssystem zur Berechnung der Knoten-Potentiale.
• D.h., je höherwertig die Ansatzfunktion gewählt wird, desto besser ist zwischen den Knoten die Approximation an die realen Potentialverläufe (Verschiebungen in der Mechanik).
• Man kann demzufolge gleiche Genauigkeiten mit einer gröberen Vernetzung erzielen. Diesem Vorteil wirkt jedoch die höhere Anzahl von Knoten pro Element wieder entgegen.
• Die benötigte Zeit für die Lösung des Gleichungssystems steigt quadratisch mit der Zahl der Unbekannten, welche proportional mit der Knotenanzahl steigt:
  • Für Potentialprobleme wird nur eine Unbekannte pro Knoten berechnet, in der Mechaniksimulation sind es 1 bis 6 Unbekannte (abhängig von den berücksichtigten Freiheitsgraden).
  • Demzufolge hat man bei Potentialproblemen mehr Spielraum für eine feinere Vernetzung.
  • Im Sinne der benötigten Berechnungszeit haben sich quadratische Ansatzfunktionen mit einer entsprechend feinen Vernetzung als optimal erwiesen.

FEMAP+MEANS

• Infolge der FEMAP-Begrenzung auf 300 Knoten sind sinnvolle Netze mit quadratischen Ansatzfunktionen im Rahmen der Lehrveranstaltung nicht realisierbar.
• MEANS bietet die Möglichkeit, nachträglich aus linearen Elementen durch Einfügen von Zwischenknoten die Berechnung mit quadratischen Ansatzfunktionen durchzuführen. Dazu muss in FEMAP eine Variable midnodes mit dem Wert=1 definiert werden (Tools - Variables). Allerdings funktioniert dies nur, wenn man die Belastungen mittels (Load - Elemental) definiert, weil vom FEM-Solver nur dann die Last anteilig an alle Element-Knoten aufgeteilt werden kann!
• Der Menüpunkt (Means - Parameter allgemein) bietet eine komfortable Möglichkeit, den Wert der Variablen für "Mittenknoten erzeugen" zu setzen.

Auswirkung auf das Postprocessing

• Vom FEM-Solver (hier MEANS) werden dem Postprocessor (hier FEMAP) nur die Werte für die Knoten zur Verfügung gestellt, wobei prinzipiell unterschiedliche Varianten für diese Wertübergabe existieren:
  • Die Primär-Ergebnisse der Knoten (die berechneten Unbekannten = Potentialgrößen) werden immer übergeben.
  • Von MEANS werden für die Knoten auch Sekundär-Ergebnisse berechnet (z.B. mechanische Spannungen bzw. die Flussgrößen). Diese Berechnung könnte aber auch erst im Postprocessing erfolgen.
  • Von MEANS werden zugehörig zu jedem Element die Knotenwerte übergeben. D.h., für gemeinsame Knoten benachbarter Elemente erfolgt die Wertübergabe mehrmals. Das hat Bedeutung, wenn benachbarte Elemente unterschiedliche Material-Eigenschaften besitzen. Dann sind zwar die Primär-Ergebnisse an gleichen Knoten identisch, aber die Sekundär-Ergebnisse sind abhängig vom betrachteten Element (Zur Erinnerung: Probleme an Grenzschichten!).
• Die Berechnung der Ergebnis-Visualisierung (z.B. Contour bzw. Deformation) sollte sinnvoller Weise unter Nutzung der Ansatz-Funktionen erfolgen. Das ist aber in FEMAP zumindest bis zur Version 9.1 nicht der Fall. Dort erfolgt scheinbar nur eine lineare Interpolation zwischen den Knotenwerten, was in Abhängigkeit von der Elementform zu einer geringen Verfälschung des dargestellten Feldverlaufes führen kann.
• An Grenzschichten profitiert die FEMAP-Visualisierung von der Datenübergabe des MEANS-Solvers, da für die aneinander grenzenden Elemente der Sprung in den Flussgrößen (z.B. mech. Spannungen) auch dargestellt werden kann. Dies ist nicht bei allen FEM-Solvern üblich!

===>>> Hier geht es bald weiter !!!

← →