Homogenization of rod-like metamaterials as a special Cosserat rod

Dieser Beitrag stellt ein variationsbasiertes Homogenisierungsschema vor, das auf der geometrisch exakten speziellen Cosserat-Stabtheorie basiert, um die nichtlineare konstitutive Antwort und Steifigkeit periodisch angeordneter stabförmiger Metamaterialien abzuleiten, wobei dies durch numerische Beispiele von einfachen Gittern bis hin zu komplexen auxetischen und künstlichen Muskelstrukturen validiert wird.

Das Wesentliche

Stellen Sie sich vor, Sie haben einen langen, flexiblen Schlauch, der nicht aus einem einzigen festen Material besteht, sondern aus einem komplexen, sich wiederholenden Muster winziger Stäbchen, die miteinander verbunden sind, wie eine mikroskopische Leiter oder ein Maschendrahtzaun, der zu einem Zylinder gewickelt ist. Dies ist das, was die Autoren als „stabförmiges Metamaterial" bezeichnen.

Das Problem, das sie angehen, lautet: Wenn Sie wissen möchten, wie sich dieser gesamte lange Schlauch biegt, dehnt oder verdreht, können Sie nicht einfach nur einen einzelnen winzigen Stab betrachten. Sie müssen betrachten, wie das gesamte Netzwerk aus Tausenden von Stäbchen interagiert. Die Simulation jedes einzelnen Stabs für einen langen Schlauch ist wie der Versuch, jedes Sandkorn an einem Strand zu zählen, um zu verstehen, wie sich der Strand im Wind bewegt – es erfordert zu viel Rechenleistung.

Die Autoren schlagen einen cleveren Abkürzungsweg vor, ein „Rezept", um vorherzusagen, wie sich der gesamte Schlauch verhält, indem nur ein winziges, repräsentatives Stück davon untersucht wird. So gehen sie vor, erklärt mit einfachen Analogien:

1. Der „magische Zoom" (Homogenisierung)

Stellen Sie sich das Metamaterial als ein riesiges, sich wiederholendes Tapetenmuster vor. Anstatt die gesamte Wand zu analysieren, betrachten Sie nur ein einziges Quadrat der Tapete (das sogenannte RVE oder repräsentatives Volumenelement).

Der Trick der Autoren besteht darin anzunehmen, dass, wenn Sie den gesamten langen Schlauch dehnen oder verdrehen, auch dieses winzige quadratische Stück sich dehnt oder verdreht, jedoch auf eine sehr spezifische, spiralförmige Weise. Sie nennen dies eine „helikale" Verformung. Stellen Sie sich vor, Sie ziehen an einer Feder; die Windungen werden nicht nur länger, sie drehen sich auch leicht. Die Autoren erkannten, dass sie, indem sie dieses winzige quadratische Stück zwingen, genau diese spiralförmige Bewegung nachzuahmen, herausfinden können, wie der gesamte lange Schlauch reagieren würde, ohne das Ganze simulieren zu müssen.

2. Die „perfekt flexiblen" Stäbe

Die meisten Computermodelle behandeln Stäbe als steif und unveränderlich, wie ein Lineal aus Stahl. Aber im echten Leben, besonders bei diesen winzigen Metamaterialien, können die Stäbe gleichzeitig biegen, dehnen und scheren (seitwärts gleiten), selbst wenn die Verformung enorm ist.

Die Autoren verwenden ein spezielles mathematisches Modell, das als „Spezieller Cosserat-Stab" bezeichnet wird.

• Analogie: Stellen Sie sich ein Stück gekochte Spaghetti vor. Es kann biegen, es kann sich ein wenig dehnen und es kann sich verdrehen. Stellen Sie sich nun vor, diese Spaghetti besteht aus einem Material, das all diese Dinge perfekt und genau tun kann, selbst wenn Sie sie zu einem Kreis biegen oder auf das Doppelte ihrer Länge dehnen. Das ist es, was ihr Modell leistet. Es approximiert nicht nur; es erfasst die exakte Geometrie der Biegung und Verdrehung.

3. Die „Tanzboden"-Regeln (Randbedingungen)

Um das winzige quadratische Stück so zu verhalten zu lassen, als wäre es Teil eines riesigen, sich wiederholenden Schlauchs, mussten die Autoren eine Reihe von Regeln erfinden, wie die Ränder dieses Quadrats miteinander kommunizieren.

• Das Problem: Wenn Sie ein Stück einer Wendeltreppe abschneiden, passt die obere Kante nicht perfekt zur unteren Kante.
• Die Lösung: Sie schufen eine „helikale Randbedingung". Stellen Sie sich vor, die linke Seite Ihres winzigen quadratischen Stücks hält die rechte Seite an der Hand, aber die rechte Seite ist leicht gedreht und verschoben, genau wie Stufen auf einer Wendeltreppe.
• Die Innovation: Frühere Methoden konnten nur kleine, sanfte Bewegungen handhaben. Die neue Regel der Autoren funktioniert auch dann, wenn der Schlauch zu einem Brezel verdreht oder so weit gedehnt wird, dass er dünn wie ein Faden ist. Sie ist „geometrisch exakt", was bedeutet, dass sie niemals an Genauigkeit verliert, egal wie wild die Form wird.

4. Die „Gelenke" und der „Kleber"

Innerhalb dieses winzigen quadratischen Stücks sind die Stäbe an Gelenken verbunden.

• Steife Gelenke: Einige Gelenke sind wie superstarker Kleber; die Stäbe können sich am Verbindungspunkt nicht relativ zueinander bewegen.
• Die Mathematik: Die Autoren haben ein System aufgebaut, bei dem der Computer die beste Position jedes Stabs in diesem winzigen Quadrat berechnet, wobei sichergestellt wird, dass die Gelenke verbunden bleiben und die „Wendeltreppe"-Regeln eingehalten werden, während gleichzeitig der geringstmögliche Energieaufwand verwendet wird.

5. Was sie fanden (Die Ergebnisse)

Sobald sie die Mathematik für das winzige Stück gelöst hatten, konnten sie vorhersagen, wie sich der gesamte Schlauch verhalten würde. Sie testeten dies mit verschiedenen Formen:

• Die Kreuz- und Quadratformen: Sie begannen mit einfachen Formen (wie einem Pluszeichen oder einem Quadrat aus Stäben), um zu beweisen, dass ihre Mathematik funktionierte. Sie stellten fest, dass, wenn die winzigen Stäbe dick und kurz sind, es sehr wichtig ist, ob sie sich dehnen oder scheren können. Wenn sie sehr dünn und lang sind, funktioniert die alte, einfachere Mathematik gut.
• Die helikalen (federförmigen) Stäbe: Sie betrachteten ein Quadrat aus Stäben, die bereits wie Federn (Helices) gekrümmt sind.
  • Die „J-förmige" Dehnung: Als sie dieses Material zogen, war es zunächst weich (wie eine sich entrollende Feder), wurde aber sehr steif, sobald es sich gerade zog. Dies erzeugt eine „J-förmige" Kurve. Dies ist genau das Verhalten biologischer Gewebe (wie Muskeln), weshalb die Autoren erwähnen, dass diese für künstliche Muskeln verwendet werden könnten.
  • Die Erweichung beim Biegen: Als sie es bogen, wurde das Material weicher, je mehr sie es bogen. Dies geschah, weil der verbindende Federstab aus der Ebene heraus zu verdrehen begann und wie ein Scharnier wirkte.
• Das auxetische Rohr: Sie modellierten ein hohles Rohr, das sich verbreitert, wenn man es zieht (wie eine Wabe).
  • Sie zeigten, dass man durch Ändern des Winkels der Stäbe das Rohr so abstimmen kann, dass es sehr flexibel von Seite zu Seite ist (gut zum Biegen), aber sehr steif gegen das Zusammendrücken ist (gut zum Offenhalten von Blutgefäßen).
  • Sie stellten fest, dass diese Strukturen so abgestimmt werden können, dass sie „Verkürzung" vermeiden (kurzer werden, wenn sie erweitert werden), was ein häufiges Problem bei kardiovaskulären Stents (Gitterrohren, die verwendet werden, um Arterien offen zu halten) ist.

Zusammenfassung

Die Autoren bauten einen „universellen Übersetzer" für Metamaterialien. Sie entwickelten eine Methode, die ein komplexes, dreidimensionales Netzwerk winziger Stäbe in eine einfache, glatte mathematische Beschreibung eines einzelnen Stabs übersetzt. Dies ermöglicht es Ingenieuren, komplexe, flexible Materialien für Dinge wie robotische Arme, künstliche Muskeln und medizinische Stents zu entwerfen, indem sie die winzigen inneren Muster anpassen, wobei sie genau wissen, wie das Endprodukt sich biegen und dehnen wird, ohne für jede einzelne Designänderung eine Supercomputer-Simulation durchführen zu müssen.

Technische Zusammenfassung

Technische Zusammenfassung: Homogenisierung von stäbchenförmigen Metamaterialien als spezieller Cosserat-Stab

Problemstellung
Stäbchenförmige Metamaterialien sind Strukturen, die durch periodische Anordnung mikroskopischer Einheiten (typischerweise Netzwerke aus Stäben) in einer einzigen Richtung entstehen, wodurch eine makroskopische Geometrie resultiert, die in einer Dimension deutlich länger ist als in den anderen beiden. Diese Strukturen finden Anwendungen im Bereich von kardiovaskulären Stents und künstlichen Muskeln bis hin zu robotischen Aktuatoren. Während eine direkte numerische Simulation solcher heterogener Strukturen im Vollmaß aufgrund des hohen Rechenaufwands untragbar ist, stoßen Standard-Homogenisierungstechniken häufig an Grenzen. Insbesondere sind Homogenisierungsschemen erster Ordnung unzureichend, um große Verformungen mit signifikanten Dehnungsgradienten (wie Biegung und Verdrehung), die für stäbchenförmige Strukturen typisch sind, zu erfassen. Darüber hinaus versagen Standard-Ansätze der dreidimensionalen Homogenisierung, Größeneffekte zu berücksichtigen, wenn die Trennung von Längenskalen in den transversalen Richtungen nicht gegeben ist. Es besteht ein Bedarf an einem Rahmenwerk, das die nichtlineare konstitutive Antwort dieser Metamaterialien unter beliebigen großen Verformungen präzise modellieren kann, wobei endliche Scherung, Dehnung, Biegung und Verdrehung sowohl auf Mikro- als auch auf Makroebene berücksichtigt werden.

Methodik
Die Autoren schlagen ein computergestütztes Homogenisierungsschema vor, das stäbchenförmige Metamaterialien auf Makroebene als „speziellen Cosserat-Stab" modelliert, während die mikroskopische Einheit (repräsentatives Volumenelement oder RVE) als Netzwerk geometrisch exakter Cosserat-Stäbe modelliert wird.

1. Geometrisch exakte Kinematik: Sowohl der makroskopische effektive Stab als auch die mikroskopischen Komponentenstäbe werden unter Verwendung der geometrisch exakten Theorie spezieller Cosserat-Stäbe modelliert. Diese Formulierung berücksichtigt endliche Biegung, Verdrehung, Scherung und axiale Dehnung und vermeidet die Annahmen kleiner Dehnungen klassischer Balkentheorien.
2. Erweiterung der Helical Cauchy-Born (HCB)-Regel: Die Methode erweitert die Helical Cauchy-Born-Regel. Sie geht davon aus, dass das Metamaterial auf Makroebene einem homogenen Dehnungsfeld entlang seiner Bogenlänge unterliegt. Diese Annahme ermöglicht die Reduktion des gesamten Strukturproblems auf ein einzelnes RVE.
3. Helisch periodische Randbedingungen: Eine wesentliche Neuerung ist die Herleitung einer geometrisch exakten, nichtlinearen Abbildung zur Anwendung helisch periodischer Randbedingungen auf das RVE. Diese Abbildung verknüpft die translatorischen und rotatorischen Freiheitsgrade der rechten Grenze des RVE mit denen der linken Grenze basierend auf den makroskopischen Dehnungsparametern ($v_0, k_0$). Im Gegensatz zu früheren linearen Übergangsabbildungen bleibt diese Formulierung auch bei beliebig großen makroskopischen Dehnungen gültig.
4. Eingeschränkte Energie-Minimierung: Das RVE-Problem wird als eingeschränkte Minimierung der elastischen Energie formuliert. Die Einschränkungen umfassen:
   • Interne Gelenkbedingungen: Erzwingen einer starren Verbindung zwischen den Stäben innerhalb des RVE mittels eines Master-Slave-Ansatzes.
   • Helische Randbedingungen: Erzwingen der aus der HCB-Regel abgeleiteten Periodizität.
   • Globale Bedingungen: Einschränkung der Starrkörpertranslation und -rotation des RVE, um eine eindeutige Lösung sicherzustellen und während der Minimierung feste makroskopische Dehnungsparameter beizubehalten.
5. Diskretisierung und Lösung: Die RVE-Energie wird unter Verwendung „helischer Elemente" diskretisiert, wobei angenommen wird, dass jedes Segment einem homogenen Dehnungsfeld unterliegt, wodurch das Segment selbst helisch wird. Die resultierende nichtlineare schwache Form wird mit einer hausinternen Software auf der deal.II-Plattform gelöst.
6. Herleitung makroskopischer Eigenschaften: Geschlossene Ausdrücke für die homogenisierte innere Kontaktkraft, das Moment und den Steifigkeitstensor des Stabs werden durch Differentiation der minimierten RVE-Energiedichte nach den makroskopischen Dehnungsparametern abgeleitet.

Hauptbeiträge
Die Arbeit identifiziert drei primäre neuartige Beiträge:

1. Geometrische Exaktheit auf zwei Skalen: Sowohl Mikro- als auch Makrostäbe werden als geometrisch exakte spezielle Cosserat-Stäbe modelliert, wodurch neben Biegung und Verdrehung auch endliche Scherung und Dehnung auf beiden Skalen erfasst werden.
2. Geometrisch exakte Makro-zu-Mikro-Abbildung: Es wird eine vollständig nichtlineare Abbildung zur Anwendung helisch periodischer Randbedingungen auf translatorische und rotatorische Freiheitsgrade vorgeschlagen. Diese Abbildung ist auch bei großen makroskopischen Dehnungen gültig und adressiert eine Einschränkung früherer linearer Übergangsabbildungen.
3. Beliebiges mikroskopisches konstitutives Verhalten: Das Rahmenwerk ermöglicht beliebiges konstitutives Verhalten der Komponentenstäbe innerhalb des RVE und geht über lineare elastische Annahmen hinaus.

Ergebnisse und numerische Beispiele
Die Autoren validieren das Schema anhand mehrerer numerischer Beispiele mit variierender RVE-Komplexität:

• Kreuz- und quadratische RVEs: Einfache planare RVEs (Kreuz- und Quadratformen) wurden verwendet, um die Ergebnisse gegen bestehende analytische Formeln und die Literatur zu validieren. Die Studie zeigte, dass die Einbeziehung von Scherung und Dehnung in Mikrostäben die makroskopische Steifigkeit erheblich beeinflusst, insbesondere bei Stäben mit größeren Schlankheitsverhältnissen. Die Ergebnisse stimmten mit analytischen Grenzfällen für Kirchhoff-Stäbe (undehnbar/unverformbar) überein.
• Helische Stab-RVEs: Eine komplexere quadratische RVE, bestehend aus helischen Komponentenstäben, wurde analysiert. Die Ergebnisse zeigten eine J-förmige makroskopische Spannungs-Dehnungs-Antwort, die von einem biegungsdominierten zu einem dehungsdominierten Verhalten übergeht. Die Studie hob hervor, wie Parameter wie die Anzahl der helischen Windungen und die Steigung diese Antwort einstellen können. Eine Biegeanalyse offenbarte ein weichwerdendes Verhalten, das durch die außerebene Biegung der kreuzenden helischen Stäbe ausgelöst wird.
• Auxetische tubuläre Strukturen: Das Schema wurde auf auxetische tubuläre Metamaterialien (relevant für Stents) angewendet. Das Modell erfasste erfolgreich die Abhängigkeit der makroskopischen Biegesteifigkeit, der Poisson-Zahl (Verkürzungsverhalten) und der radialen Steifigkeit vom Designwinkel $\theta_a$. Es zeigte zudem die Erfassung des Bazier-Effekts (Ovalisierung des Querschnitts) unter großer Biegung.
• FE2-Validierung: Ein Vergleich zwischen dem vorgeschlagenen Homogenisierungsschema (FE2-Ansatz) und Vollmaßsimulationen einer Kreuz-RVE-Metamaterialstruktur zeigte eine gute Übereinstimmung und validierte die Genauigkeit der Methode selbst unter großen Verdrehmomenten, bei denen lineare konstitutive Modelle versagten.

Bedeutung und Behauptungen
Die Arbeit behauptet, ein robustes Rahmenwerk für die Homogenisierung stäbchenförmiger Metamaterialien zu liefern, das die Einschränkungen von Homogenisierungsschemen erster Ordnung und linearer Homogenisierung überwindet. Durch die Nutzung einer geometrisch exakten Theorie spezieller Cosserat-Stäbe auf beiden Skalen und die Herleitung einer nichtlinearen Makro-zu-Mikro-Übergangsabbildung sagt die Methode die nichtlineare konstitutive Antwort dieser Strukturen unter großen Verformungen präzise voraus.

Die Autoren betonen, dass ihr Ansatz allgemein ist und auf andere RVE-Typen (z. B. poröse Festkörper, Origami-Strukturen) erweitert werden kann und komplexe mikroskopische Phänomene wie Elastoplastizität, Viskoelastizität und Multiphysik-Kopplung (thermo-, chemo- oder elektroelastisch) einbeziehen kann. Die Arbeit dient als Grundlage für die Entwicklung maßgeschneiderter Metamaterialien mit spezifischen makroskopischen Eigenschaften, wie künstlichen Muskeln oder einsetzbaren Stents, durch Abstimmung mikroskopischer Parameter. Die Autoren vermerken bescheiden, dass die aktuelle Arbeit zwar davon ausgeht, das RVE sei ein Stabnetzwerk, die abgeleiteten Randbedingungen jedoch für andere mikroskopische Darstellungen angepasst werden könnten.