Kontext
Individuelle Projektarbeit im ersten Masterjahr (MESI, Modul Finite Elemente) an der Fakultät für Physik & Ingenieurwesen der Universität Straßburg, betreut von Dr. Salah Elbarnaty und Prof. J.P.M. Correia. Der Untersuchungsfall: ein 2 m langer Kragträger (Querschnitt 200 mm, Dicke 1 mm, Stahl) unter einer Kraft von 250 N, unter der Annahme des ebenen Spannungszustands. Ziel: ein linear-elastisches Modell validieren, indem die ABAQUS-Simulation mit der analytischen Lösung der technischen Mechanik verglichen wird, und verstehen, wie Netzdichte und Elementformulierung die Genauigkeit bestimmen.
Ansatz / Lösung
Zuerst habe ich die analytische Referenz aufgestellt (Euler-Bernoulli): Flächenträgheitsmoment, maximale Normalspannung (±75 MPa, deutlich unter der Streckgrenze von 400 MPa) und Durchbiegung am Trägerende (4,762 mm). Diese Lösung ist der Maßstab für jede numerische Konfiguration. In ABAQUS habe ich das Modell sorgfältig aufgebaut: Geometrie, Werkstoff Stahl (E = 210.000 MPa), konsistente Einheiten (mm / MPa / N), die Einspannung bereits im Initial-Step definiert und Sets verwendet, um Randbedingungen und Last einmal aufzubringen und das Modell anschließend über die acht Konfigurationen zu duplizieren. Ich habe ebene-Spannungs-Elemente CPS4 (lineare Vierecke) und CPS8 (quadratische) verwendet. Anschließend habe ich zwei Faktoren kombiniert: die Netzdichte (1×6, 2×12, 4×12, 8×24) und die Interpolationsordnung (CPS4 vs. CPS8). Für jeden Fall habe ich die S11-Spannungskarte, die U2-Verformung, die Knotenreaktionen RF2 und den relativen Fehler der Durchbiegung extrahiert und in einer Vergleichstabelle zusammengefasst.
Ergebnisse
Das herausragende Ergebnis: Mit nur 6 quadratischen CPS8-Elementen erreiche ich einen Fehler von 0,33 % bei der Durchbiegung, besser als 192 lineare CPS4-Elemente (0,94 %). Bei CPS4 ist die Konvergenz monoton, aber langsam (40,7 % Fehler bei 1×6, bis auf 0,9 % bei 8×24); bei CPS8 ist sie bereits beim gröbsten Netz praktisch erreicht. Außerdem habe ich das Schubversteifen (Shear Locking) offengelegt und verstanden: Bei einem groben CPS4-Netz sind die linearen Elemente zu steif, um die Biegekrümmung abzubilden, und unterschätzen die Durchbiegung um 40 %. Zur Validierung beträgt die Auflagerreaktion stets exakt 250 N (globales Gleichgewicht), und die S11-Verteilung bleibt symmetrisch zur neutralen Faser, zwei Kontrollen, die ich jedes Mal überprüft habe. Mein Fazit: Eine an die Physik angepasste Elementformulierung zu wählen, ist wirksamer als blindes Netzverfeinern, eine wertvolle Denkweise im industriellen Umfeld, wo Rechenzeit zählt. Diese Projektarbeit war zugleich mein erster praktischer Umgang mit ABAQUS und hat mir verständlich gemacht, warum ein solcher Solver unverzichtbar wird, sobald man den linearen Fall verlässt (Nichtlinearitäten, Kontakt, Crash), wo CAD-Werkzeuge wie Inventor Nastran an ihre Grenzen stoßen. Wege, die ich noch erkunden könnte: ins Nichtlineare vorstoßen (große Verschiebungen, Plastizität, Kontakt), mit Balken- oder Schalenelementen vergleichen, die reduzierte Integration und ihre Hourglassing-Moden untersuchen und ABAQUS am selben Fall gegen einen anderen Solver benchmarken.