3D Finite Element-Based Multiphysics Simulation of a Shape Memory Alloy Hybrid Composite Module

Titel: 3D-Finite-Elemente-basierte Multiphysik-Simulation eines Hybridverbundmoduls aus Formgedächtnislegierungen (SMAHC)

Autoren: L. Handl, M. Kaiser, M. Duhovic, M. Gurka (Leibniz-Institut für Verbundwerkstoffe GmbH, Deutschland)

---

1. Problemstellung und Motivation

Formadaptive Hybridverbundwerkstoffe aus Formgedächtnislegierungen (SMAHCs) sind Verbundmaterialien, die SMA-Drahte integrieren, um reversible und kontinuierliche Formveränderungen zu ermöglichen. Diese werden zunehmend für Anwendungen wie aktive Spoiler, adaptive Trailing Edges oder Robotergreifer eingesetzt.

Obwohl zahlreiche analytische und Finite-Elemente-(FE)-Modelle zur Vorhersage des transienten Verhaltens von SMAHCs existieren, weisen viele bestehende Ansätze erhebliche Einschränkungen auf:

• Mangelnde Berücksichtigung der elektro-thermomechanischen Kopplung (Wechselwirkung zwischen elektrischem Strom, Erwärmung und mechanischer Verformung).
• Unzureichende Validierung durch umfassende experimentelle Daten.
• Oft vereinfachte Modelle, die keine detaillierten 3D-Geometrien oder Randbedingungen (wie Umgebungstemperatur und mechanische Lasten) abbilden.

Das Ziel dieser Studie war es, eine vollständig gekoppelte Multiphysik-Simulation in 3D zu entwickeln, die diese Lücken schließt und das dynamische Verhalten eines SMAHC-Aktors präzise vorhersagt.

2. Methodik und Modellansatz

Simulationsumgebung:
Die Simulation wurde mit dem kommerziellen FE-Code ANSYS LS-DYNA durchgeführt. Dieser nutzt integrierte gekoppelte Solver für Mechanik, Thermik und Elektromagnetismus (EM).

Materialmodell:
Es wurde ein mikromechanisches konstitutives Modell nach Kelly et al. [25] implementiert. Dieses Modell beschreibt das Materialverhalten auf der Mikroebene durch die Minimierung der Helmholtz-Freien-Energie ($W$).

• Es berücksichtigt die Kinetik der Phasenumwandlung zwischen Martensit und Austenit (Initiation, Sättigung, Wachstum).
• Der Zustand wird durch innere Variablen wie den Martensit-Volumenanteil ($\lambda$) und die nominelle Martensit-Dehnung ($\varepsilon_m$) beschrieben.
• Das Modell erfasst Pseudoelastizität, Superelastizität und den Einweg-Formgedächtniseffekt.

Geometrie und Modellierung:

• Objekt: Ein kommerziell erhältlicher SMAHC-Aktor (Typ A3950 von CompActive GmbH).
• Geometrie: Der FE-Modellansatz nutzt Symmetrien und modelliert einen Querschnitt mit einem einzelnen SMA-Draht. Es wurden 32.360 hexaedrische Volumenelemente verwendet.
• Elektromagnetik: Ein vereinfachter EM-Solver berechnet den ohmschen Widerstand und die Joule'sche Erwärmung (ohne Induktionseffekte), die als Wärmequelle in den thermischen Solver eingespeist wird.

Besonderer Initialisierungs-Schritt (Pre-Strain):
Ein kritischer Aspekt des Modells ist die physikalisch fundierte Initialisierung des Materialzustands. Da der Standard-Materialcard in LS-DYNA den SMA-Draht initial als reinen Austenit definiert, wurde ein vorgeschalteter Simulations-Schritt eingeführt:

1. Der Draht wird bei Raumtemperatur (23 °C) mechanisch gedehnt.
2. Dies führt zu einer Umwandlung von Austenit in verzwillingten Martensit und anschließend zu einem partiell entzwillingten Martensitzustand.
3. Dieser Schritt ist essenziell, da das Kontraktionspotenzial des Drahtes auf der Umwandlung von entzwillingtem Martensit zu Austenit beruht. Ohne diesen Schritt wäre das reale Verhalten nicht abbildbar.

Validierung:
Die Ergebnisse wurden mit zwei Referenzen verglichen:

1. Experimentelle Daten (verschiedene Lasten, Umgebungstemperaturen und Ströme).
2. Ein vollständig gekoppeltes, transienter Stufen-Schema-Modell (SSM) [22], das als numerischer Referenzstandard diente.

3. Wichtige Ergebnisse

Thermisches Verhalten:

• Die Joule'sche Erwärmung wurde qualitativ gut vorhergesagt. Die Temperaturverteilung entlang des Drahtes ist inhomogen (Mitte heizt schneller auf), was experimentellen Beobachtungen entspricht.
• Abweichung beim Abkühlen: Die Simulation zeigt eine signifikant schnellere Abkühlung als im Experiment. Dies führt dazu, dass die Rückkehr in die Ausgangsposition in der Simulation früher einsetzt. Dies wird auf eine Unterschätzung der thermischen Trägheit (z. B. durch Konvektionseffekte) zurückgeführt.

Mechanisches Verhalten (Durchbiegung):

• Qualitative Übereinstimmung: Die Simulation reproduziert die charakteristische Hysterese der Durchbiegung in Abhängigkeit von der Temperatur erfolgreich.
• Quantitative Ergebnisse: Die vorhergesagten Durchbiegungen liegen in der richtigen Größenordnung. In den meisten Fällen liegen die Ergebnisse innerhalb des 95 %-Konfidenzintervalls der Experimente, weichen jedoch oft marginal davon ab.
• Einfluss der Last: Bei höheren mechanischen Lasten (bis 39,2 N) zeigt die Simulation eine zu hohe Steifigkeit des Aktuators. Die simulierten Hysteresekurven sind verschoben, und die pseudoplastische Verformung des Drahtes wird weniger stark ausgeprägt als im Experiment.
• Einfluss der Geometrie: Bei Aktoren mit dünneren Elastomer-Schichten (Typ A und B) traten größere Abweichungen auf. Die Simulation konnte den experimentell beobachteten Anstieg der maximalen Durchbiegung bei dünneren Schichten nicht vollständig abbilden, was vermutlich auf eine Überschätzung der Steifigkeit der weichen Zwischenschicht im linearelastischen Modell zurückzuführen ist.

Vergleich mit SSM:
Das Stufen-Schema-Modell (SSM) zeigte in den meisten Fällen eine bessere quantitative Übereinstimmung mit den experimentellen Daten als die 3D-FE-Simulation, insbesondere bei der Hysterese und den Durchbiegungswerten.

4. Hauptbeiträge und Signifikanz

1. Erste vollständige 3D-Multiphysik-Simulation: Dies ist einer der ersten Ansätze, der mechanische, thermische und elektromagnetische Kopplungen in einem detaillierten 3D-FE-Modell für SMAHCs integriert.
2. Physikalisch fundierte Initialisierung: Die Einführung des Pre-Strain-Schritts zur Erzeugung eines realistischen, entzwillingten Martensitzustands ist ein wesentlicher methodischer Fortschritt für die Genauigkeit der Vorhersage.
3. Validierung unter realen Bedingungen: Das Modell wurde nicht nur unter idealisierten Bedingungen, sondern unter variierenden Randbedingungen (Umgebungstemperatur, externe Lasten, verschiedene Ströme) validiert.
4. Potenzial für komplexe Systeme: Trotz der identifizierten Abweichungen (insbesondere beim Abkühlen und bei hohen Lasten) demonstriert die Studie, dass komplexe 3D-SMAHC-Systeme simuliert werden können. Dies eröffnet Möglichkeiten für den CAE-gestützten Entwurf (Computer-Aided Engineering) von intelligenten Strukturen.

5. Fazit und Ausblick

Die Studie belegt, dass der neue Materialansatz in ANSYS LS-DYNA geeignet ist, das dynamische thermo-mechanische Verhalten von SMAHC-Aktoren qualitativ gut und quantitativ in der richtigen Größenordnung vorherzusagen. Die größten Diskrepanzen liegen in der Modellierung der Abkühlung und der Steifigkeit der Polymerkomponenten bei hohen Lasten.

Zukünftige Arbeiten sollen sich auf folgende Punkte konzentrieren:

• Nutzung der räumlichen Auflösung des Martensit-Volumenanteils ($\lambda$) für detailliertere Analysen.
• Entwicklung verbesserter Materialmodelle für die Elastomer-Zwischenschicht (nicht-linear).
• Sensitivitätsanalysen von Materialparametern und geometrischen Abmessungen zur Optimierung der Modellgenauigkeit.