Software: FEM - Tutorial - 2D-Mechanik - manuell - strukturiertes Netz
In einem FEM-Programm werden verschiedene Methoden für die Erzeugung von Netzen bereitgestellt. Wir haben bisher nur die sogenannte freie Vernetzung (engl: free mesh) verwendet.
Free Mesh:
- Innerhalb eines FEM-Programms werden meist unterschiedliche Free-Mesher bereitgestellt, die man jeweils umfangreich konfigurieren kann (im Autodesk Simulation Multiphysics z.B. 3D-Mesher und 2D-Mesher).
- In Hinblick auf Flächennetze kann man mit einem Free-Mesher beliebige Oberflächen automatisch vernetzen. Netzverfeinerungen lassen sich mittels unterschiedlicher Parameter steuern.
- Die erzeugten Netze sind jedoch nicht immer optimal an die Fläche und die Spannungsverläufe angepasst. In unserem Beispiel äußert sich das z.B. in einer asymmetrischen Element-Anordnung trotz einer achsensymmetrischen Fläche.
Ziel dieses Abschnitts soll ein optimiertes 2D-Netz sein, welches mit möglichst geringem Berechnungsaufwand hinreichend genaue Ergebnisse liefert:
- Da unser Bauteil symmetrisch zur Mittelachse ist und nur in Richtung dieser Mittelachse belastet wird, genügt das Netz einer Symmetriehälfte.
- Wir verwenden die sogenannte strukturierte Vernetzung (engl: mapped mesh), um ein effizientes 2D-Netz zu erzeugen. Die strukturierte Vernetzung von Flächen wird meist mit Viereck-Elementen durchgeführt.
Mapped Mesh:
Angelehnt an die Kartografie versucht man bei dieser Vernetzungsmethode, Oberflächen mit einem Gitternetz zu überziehen, dessen Rechteck-Raster möglichst wenig verzerrt wird. Das gelingt bei komplexeren Oberflächen nur, wenn man die Gesamtfläche in einfachere 4-seitige Teilflächen zerlegt, welche auch gekrümmte Kanten besitzen können. Das Ziel besteht darin, über alle Teilflächen hinweg fluchtende Reihen- und Spalten von Elementen zu erhalten, z.B.:
Um bereits eingegebene Informationen möglichst weiter zu nutzen, erzeugen wir eine Kopie vom vorhandenen Free-Mesh-Szenarium und konfigurieren daraus ein Mapped-Mesh-Szenarium. Nach dem 2D-Netz-Löschen besteht die Hauptarbeit im Ändern der Geometrie mit den Mitteln der FEM-eigenen Zeichen-Funktionen:
Geometrische Vorbereitung
Hinweise:
- Auswahl-Modus "Konstruktionsobjekte" muss vor einer Wahl von Objekten aktiv sein, damit nicht aus Versehen z.B. Rasterpunkte gefangen werden.
- "Für Konstruktion verwenden" muss für neue Linien aktiviert sein, damit Konstruktionsobjekte entstehen. Nur solche sind für die strukturierte Vernetzung nutzbar.
- Geometrie auf Symmetriehälfte "zurechtstutzen"
- YZ-Ebene als Skizze aktivieren
- Mittels zweier Konstruktionspunkte eine Linie entlang der Symmetrie-Achse zeichnen
- Stutzen aller Skizzen-Elemente, welche durch die Symmetrie-Linie geschnitten werden:
Randelement=Symmetrie-Linie
zu änderndes Objekt=zu stutzendes Element
zu trimmendes Segment=zu löschender Elementabschnitt - Löschen der isolierten Elemente der unteren Symmetriehälfte einschließlich des Konstruktionspunktes (0,0,0). Den Kreismittelpunkt benötigen wir noch für die Konstruktion der Netz-Teilflächen!
- Zerlegung in 4-seitige Teilflächen
Leider existieren viele Funktionen, welche man aus der Skizzier-Umgebung eines CAD-Systems kennt, nicht in den Zeichen-Möglichkeiten eines FEM-Programms. Im Beispiel wird man z.B. Skizzierabhängigkeiten oder Fangen von Mittelpunkten vermissen. Man muss sich jeweils mit den vorhandenen Funktionen irgendwie behelfen:
- Teilen Halbkreis in 2 Abschnitte, Linie als Radius im Halbkreis, Dehnen der Radius-Linie bis zur oberen Kante (Fang der Rasterlinien vermeiden!), Anschließend Stutzen des Radius-Abschnitts
- Teilen linker Viertelkreis, Linie wie zuvor beschrieben
- Rechteck bis Endpunkt vorherige Linie
- Bereinigen der Teilflächen-Geometrie
- Stutzen von Linien, denn Linien dürfen sich nur zwischen zwei Eckpunkten eines Vierecks erstrecken. Überflüssige Segmente (z.B. unter dem Rechteck) kann man durch Stutzen beseitigen
- Schneiden von Linien in Teilstücke, wenn die getrennten Segmente noch benötigt werden
- Löschen übereinander liegender Linien nach Auswahl. Im Beispiel betrifft dies nur die linke Kante des Rechtecks. Die Funktion Doppelte Linien löschen ist hierfür nicht geeignet, denn sie entfernt diese übereinander liegenden Linien nicht.
Die geometrische Vorbereitung für die strukturierte Vernetzung bot die Gelegenheit, uns vertieft in die Funktionen des 2D-Geometrie-Editors eines FEM-Programms einzuarbeiten. Leider handelt sich bei solchen Editoren meist um individuelle Lösungen, auch wenn das FEM-Programm Bestandteil einer CAD-Suite ist. Das bedingt einen relativ hohen Einarbeitungsaufwand.
Strukturiertes Vernetzen
Die Zerlegung in geeignete 4-seitige Teilflächen basierte bereits auf qualitativen Vorstellungen, wie das Netz über die gesamte Fläche des Modells gespannt werden soll. Diese qualitativen Vorstellungen müssen nun quantitativ präzisiert werden:
- Innerhalb der 4-seitigen Teilflächen müssen gegenüberliegende Seiten (im Sinne des Netzverlaufs) mit jeweils der gleichen Anzahl von Elementen vernetzt werden.
- Die Planung des Netzes sollte an den kritischen Stellen des FEM-Modells beginnen. Dazu muss man jedoch bereits qualitative Vorstellungen besitzen, wo die größten Gradienten in den berechneten Feldverläufen auftreten. Im Beispiel ist dies am linken Rand des Loches der Fall:
- Für diesen Rand trifft man die grundlegende Entscheidung, mit wievielen Elementen er zu vernetzen ist. Zur Demonstration soll das am Beispiel einer ziemlich groben Vernetzung erläutert werden.
- Entscheidet man sich an der linken Seite des Bolzens für z.B. jeweils 6 Elemente auf den anliegenden Kanten, so zieht sich diese Feinheit der Vernetzung durch das gesamte Modell:
- Die Elementezahl auf den anliegenden Kanten der Teilfläche kann man unabhängig davon wählen. Dabei sollte man jedoch darauf achten, dass die Elemente des Netzes möglichst rhombenförmig werden. Deshalb fiel die Entscheidung hier auch auf eine Unterteilung in 6 Elemente.
- Im Beispiel wurden die noch nicht spezifizierten Kanten wie folgt in Elemente unterteilt:
- Damit ergäbe sich dann die bereits zu Beginn vorgestellte Vernetzung:
In der Registerkarte MFL > Netz wird deutlich, dass für die 2D-Geometrie keine Funktionen der 3D-Vernetzung zur Verfügung stehen. Die Funktion 2D-Netz erstellen ist nur verfügbar, wenn aktuell eine Bauteil-Ebene im Browser des FEM-Editors ausgewählt ist.
Verfügbar ist in der MFL die strukturierte Vernetzung einschließlich der Möglichkeiten einer Netzverfeinerung:
Bei unserem strukturierten Netz soll die Lochkante ähnlich fein vernetzt werden, wie zuvor bei der freien Vernetzung mit dem 2D-Netzgenerator:
- Die benutzte Teilung=250 für den gesamten Kreisumfang entspricht ungefähr einer Elementzahl von 30 auf dem linken Kreis-Segment des Loches.
- Daraus würde bereits eine ziemlich feine Vernetzung des Rechteck-Abschnittes der Lasche mit 30 Element-Zeilen resultieren.
- Wir planen deshalb nur mit einer Teilung in 20 Elemente auf dem linken Kreis-Segment.
Wir haben zwar systematisch die einzelnen 4-eckigen Teilflächen des Netzes in der Geometrie nachgebildet. Abhängig von der Vernetzungsfunktion werden wir jedoch nur ein Teil der geometrischen Objekte benötigen (Auswahlmodus immer "Konstruktionsobjekte"!).
benötigt nur die 4 Eckpunkte der Teilfläche:
- Voraussetzung für die Nutzung dieser Funktion sind gerade Begrenzungslinien der Teilfläche.
- Dies ist nur auf dem geraden Teil der Lasche gegeben. Dieses Rechteck vernetzen wir zuerst.
- Die Wahl der erforderlichen Punkte A, B, C und D muss entweder im oder gegen den Uhrzeigersinn erfolgen. Man sollte sich die mathematisch positive Richtung "gegen den Uhrzeigersinn" angewöhnen!
- Eine Teilung von 20:60 würde im Beispiel quadratische Elemente bewirken. Da ein Seiten-Verhältnis der Viereck-Elemente von 1:2 akzeptabel ist, genügt eine Teilung von 20:30.
- Wir verwenden die standardmäßig eingestellte konstante Teilung entlang der Netzkanten. Prinzipiell kann man dort variable Elementgrößen konfigurieren (Button "..." hinter Teilungswert):
spannt ein 4-eckiges Teilnetz zwischen zwei gegenüberliegenden Konstruktionsobjekten AA' und BB' (Linien, Bögen, Splines):
- Damit können wir Teilnetze zwischen Bögen des Lochrandes und der Außenkante des Bauteils erstellen:
- Leider gelingt die für die rechte Teilfläche geplante Vernetzung nur teilweise, weil keine Möglichkeit gefunden wurde, den äußeren Bogen mit den zwei angrenzenden Linienstücken zu einem Objekt zusammenzufassen. Wir vernetzen deshalb zuerst das Viereck zwischen den beiden Bögen mit einer gröberen Teilung von 20:10. Es bleiben vorläufig zwei 3-eckige Teilflächen unvernetzt:
- Hinweis: Ob eine Teilfläche mit mehreren Objekten zwischen 2 Ecken mit einem strukturierten Teilnetz gefüllt werden kann, ist abhängig von den Möglichkeiten des benutzten FEM-Programms. Das einleitende Demo-Vernetzung mit 6 Elementen an den Lochsegmenten konnte z.B. mit dem Programm FEMAP realisiert werden (Siehe: FEMAP-Übung).
spannt ein Netz zwischen die Punkte A, B und C:
- Die drei Punkte wählen wir wieder entgegen dem Uhrzeigersinn jeweils innerhalb einer der Netzlücken.
- Im Punkt A werden Dreieck-Elemente erstellt. Punkt A sollte deshalb am Lochrand liegen. Das übrige Teilnetz besteht dann aus viereckigen Elementen:
Alle Teilnetze werden im Browser des FEM-Editors als Bestandteil des Szenariums aufgelistet. Man kann sie nachträglich bearbeiten bzw. löschen.
Hinweise zur Aufteilung in Teilflächen:
- Insbesondere die Vernetzung in den 3-eckigen Teilflächen führt zu ungünstigen Element-Proportionen.
- Mit dem jetzigen Wissen, dass man Vierecke mit gekrümmten Seiten nur mit der Funktion
vernetzen kann, würde man die Auteilung in Teilflächen auf der rechten Lochseite etwas anders gestalten.
- Die senkrechte Trennlinie in der Lochmitte müsste nach rechts geneigt werden, um die "senkrechte" Dreieck-Teilfläche einzusparen.
- Die "waagerechte" Dreieck-Teilfläche kann man einsparen, wenn man das rechte Segment der Lochkante ebenfalls teilt und dadurch ein weiteres Viereck gewinnt.
Entsprechend dieser Hinweise ist die rechte Seite des strukturierten Netzes weitestgehend selbstständig umzugestalten:
- Löschen der zu ändernden Teilnetze
- Ändern der Geometrie (im Skizzen-Modus der YZ-Ebene, da ansonsten 3D-Objekte entstehen!)
- Vernetzung der verbesserten 4-eckigen Teilflächen (ähnlich wie im folgenden Bild)
Lasten und Abhängigkeiten
Da wir nur eine Symmetriehälfte der Lasche als FEM-Netz modellieren, muss man die Symmetrie im Beispiel wie folgt beachten:
- Es darf auf die Modellhälfte nur die Hälfte der realen Lastkraft wirken.
- Entlang der Symmetrielinie benötigt man zusätzliche Abhängigkeiten, welche eine Bewegung der Randknoten in Z-Richtung verhindern.
Kräfte und Einspannungen müssen in Betrag und Richtung den richtigen Knoten des Netzes zugewiesen werden. Man vermeidet nach Möglichkeit, diese Zuweisung manuell durch Auswahl der betroffenen Knoten vorzunehmen:
- Eine Gesamtkraft müsste auf alle betroffenen Knoten aufgeteilt werden. Bei ungleichmäßiger Vernetzung kann man die notwendigen Einzelkräfte praktisch nicht mehr berechnen.
- Eine Neuvernetzung (z.B. zur Netzverfeinerung) würde die Zuweisung von Kräften und Randbedingungen erneut erforderlich machen.
Wichtig:
Nur durch Zuweisen von Lasten und Abhängigkeiten an Kanten oder Flächen ist automatisch die richtige Belegung der betroffenen Knoten gewährleistet. Dieser Automatismus wirkt dann auch bei einer Neuvernetzung.
Problem:
- In der aktuellen Version 2013 von Autodesk Simulation Multiphysics kann man an strukturierten Netzen keine Kanten auswählen!
- Das kann man einfach überprüfen, indem man z.B. ein Rechteck einmal mit dem 2D-Free-Mesher vernetzt (Kanten auswählbar) und dann das gleiche Rechteck nach dem Löschen des freien 2D-Netzes mit einem strukturierten Netz überzieht (Kanten nicht auswählbar).
- Wir müssen uns also mit der Zuweisung direkt an die Knoten behelfen. Das hat im Rahmen der Übung den Vorteil, dass man das Prinzip der Kraftaufteilung versteht!
1. Zugkraft von 10000 N
Aufgrund der Symmetrie reduziert sich die erforderliche Gesamtkraft auf 5000 N.
- An der (nicht wählbaren) Kante liegen 20 gleichgroße Elemente mit ihren 21 Knoten.
- Die Kraft von 5000 N darf man nicht gleichmäßig auf diese 21 Knoten verteilen (entspräche á 238,1 N):
- Die meisten Rand-Knoten des Bauteils sind Eckpunkte zweier Elemente. Die Kraft solcher Knoten verteilt sich je zur Hälfte auf seine zwei Nachbar-Elemente.
- Die Kraft der beider Eck-Knoten des Bauteils wirkt jedoch vollständig auf das zugehörige Eck-Element.
- Deshalb dürfen Eck-Knoten im Beispiel nur die halbe Kraft erhalten, wie die restlichen Rand-Knoten der Kante. Ansonsten würde an den Eck-Knoten zu stark gezogen!
- Richtig ist bei gleicher Elementgröße folgende Lastverteilung bei n Elementen (nicht n Knoten!):
- Gesamtlast/n auf jeden "normalen" Randknoten
- Gesamtlast/(2n) auf jeden Eck-Knoten
- Hinweis:
Man kann einen Knoten nicht auswählen, wenn an gleicher Stelle ein Konstruktionspunkt existiert. Deshalb muss man vor der Kraft-Zuweisung z.B. den Konstruktionspunkt in der unteren linken Ecke löschen!
2. Fixierung des Lochrandes
- Auch hier das Problem, dass Knoten unter Konstruktionspunkten nicht wählbar sind. Man muss alle Konstruktionspunkte des Lochrandes löschen!
- Es genügt die Fixierung aller Randknoten durch Einschränkung der Freiheitsgrade Ty und Tz.
3. Symmetrie-Bedingungen
- Konstruktionspunkte beachten!
- Die Knoten auf der Symmetrie-Achse dürfen sich nicht in Z-Richtung bewegen (Tz einschränken)
- Die beiden Knoten des Lochrandes, welche auf der Symmetrie-Achse liegen, besitzen bereits die Einschränkung Tz. Wird beiden Knoten erneut eine Abhängigkeit zugewiesen, so erfolgt eine Kombination aller vorhandenen Einschränkungen. Dies stört hier nicht, man könnte die beiden Knoten aber auch beim Zuweisen der Symmetrie-Bedingung aussparen.
- Bei der Auswahl der Knoten entlang der Symmetrie-Achse werden diejenigen Knoten hervorgehoben, welche bereits eine Last bzw. eine Abhängigkeit besitzen. Dies dient nur zur Information und hat keine Auswirkung auf die Zuweisung der neuen Abhängkeit.
Ergebnisse bei fixiertem Lochrand
Die simulierten Ergebnisse des strukturierten 2D-Netzes entsprechen denen aus dem freien 2D-Netz. Die wenige Prozent Abweichung beim berechneten Maximalwert resultieren aus der unterschiedlich feinen Vernetzung des Lochrandes:
Verschweisst auf Bolzen
Wird der Lochrand der Lasche mit einem Bolzen verschweißt, so ist der Lochrand nicht ideal fixiert, wie wir es vereinfachend angenommen haben:
- Infolge der Belastung verformt sich in der Realität der Bolzen.
- Infolge der Verformung des Bolzens entfallen die Spannungsspitzen am Rand des Loches.
Mit dem bisher erworbenem Wissen und einigen Hinweisen sollte es kein Problem darstellen, einen Titan-Bolzen als zweites Bauteil zu definieren und als strukturiertes Netz zu modellieren:
- Die Belastung des Titan-Bolzens soll uns im Beispiel nicht interessieren. Auf ein "Voll-Material" können wir verzichten, es genügt für den Bolzen ein Materialring der Wandstärke 1 mm ebenfalls mit der Element-Dicke=3 mm. Die Bohrung des Titan-Ringes fixieren wir, um die Bolzenbefestigung nachzubilden.
- In Autodesk Simulation Multiphysics gilt:
1 Bauteil = 1 Elementkonfiguration = 1 Material. - Geometrie mit einem zusätzlichen Material (z.B. den Titan-Bolzen) muss man als weiteres, neues Bauteil definieren (Bolzen-Material: Other > Metal > Titanium.Ti)
- Beim Skizzieren der neuen Geometrie muss man darauf achten, dass diese dem Bolzen-Bauteil zugeordnet wird.
- Als Bolzengeometrie genügt für die strukturierte Vernetzung ein Halbkreis (als Nachbildung des Bolzen-Loches), den wir in vier Segmente teilen (für vier strukturierte Teilnetze zum Lochrand der Lasche). Die beiden Linien entlang der Symmetrielinie werden für den Bolzen nicht benötigt.
- Wichtig: Für die Bolzen-Teilnetze muss man die vorhandene Loch-Geometrie der Lasche benutzen. Zeichnet und verwendet man im neuen Bauteil zusätzliche Segmente für den Außenrand des Bolzens, so entsteht ein Bolzennetz, welches nicht mit dem Netz der Lasche verbunden ist!
Aktiviert man nach der Simulation für die Ergebnisse Bolzen-Bauteil > Durchsichtig zeichnen, damit man den Lochrand der Lasche besser erkennt, so ergibt sich der folgende Mises-Spannungsverlauf in der Lasche:
Die Ergebnisse an der Materialgrenze zwischen Stahl und Titan wollen wir etwas genauer betrachten:
- Die Sprünge im Spannungsverlauf resultieren überwiegend aus den unterschiedlichen E-Modulen der beiden Materialien (Faktor 1,74).
- Die Knoten der Grenzschicht gehören gleichzeitig zu Elementen mit unterschiedlichen Eigenschaften.
- Die Knoten einer Grenzschicht besitzen zwar eine eindeutige Position, aber es lässt sich ihnen keine eindeutige Spannung zuordnen (Spannungssprung im Knoten!).
- Im benutzten FEM-Programm werden alle Spannungswerte eines Knoten an den Post-Prozess (Ergebnisauswertung) übergeben. Im Beispiel sind es zwei Spannungen, weil diese Randknoten jeweils an zwei unterschiedliche Elemente grenzen.
- Diese Spannungswerte kann man an der Cursor-Position abfragen (Prüfen), wie im Bild dargestellt:
- Der berechnete Spannungsverlauf innerhalb der Randelemente ist abhängig von der gewählten Glättungsoption für die Farbkontur:
- "Über Bauteilränder glätten" darf nicht aktiviert sein (standardmäßig gewährleistet), damit der real vorhandene Spannungssprung abgebildet wird.
Frage 3:
- Wie groß darf die maximale Zugkraft bei verschweißtem Lochrand auf nachgiebigem Titan-Bolzen sein, damit ein Sicherheitsfaktor von 2 innerhalb der Lasche nicht unterschritten wird?
- Die Lösungsdatei ist konfiguriert mit dieser zulässigen Kraft einzusenden.
Frage 4:
- Welchen Wert sollte man als maximal zuzlässige Kraft für die Einhaltung von Sicherheitsfaktor 2 benutzen, wenn die Art der Befestigung der Lasche auf dem Bolzen beliebig sein kann.
- Diese Entscheidung ist nachvollziehbar zu begründen!