Homogenization of rod-like metamaterials as a special Cosserat rod

Dit artikel presenteert een variatieel homogenisatieschema gebaseerd op de geometrisch exacte speciale Cosserat-staaftheorie om de niet-lineaire constitutieve respons en stijfheid van periodiek samengestelde staafachtige metamaterialen af te leiden, gevalideerd door middel van numerieke voorbeelden variërend van eenvoudige roosters tot complexe auxetische en kunstspierstructuren.

De kern

Stel je een lange, flexibele buis voor die niet is gemaakt van één massief materiaal, maar van een complex, zich herhalend patroon van tiny stokjes die met elkaar verbonden zijn, zoals een microscopische ladder of een schuimdraadomheining die tot een cilinder is gewikkeld. Dit noemen de auteurs een "stokachtig metamateriaal".

Het probleem dat ze aanpakken is als volgt: als je wilt weten hoe deze hele lange buis buigt, rekt of draait, kun je niet gewoon naar één klein stokje kijken. Je moet kijken naar hoe het hele netwerk van duizenden stokjes met elkaar interacteert. Het simuleren van elk enkel stokje voor een lange buis is als proberen elke korrel zand op een strand te tellen om te begrijpen hoe het strand beweegt in de wind; het kost te veel rekenkracht.

De auteurs stellen een slimme afkorting voor, een "recept" om te voorspellen hoe de hele buis zich gedraagt door slechts een klein, representatief stukje ervan te bestuderen. Hier is hoe ze dat doen, uitgelegd met eenvoudige analogieën:

1. De "Magische Zoom" (Homogenisatie)

Stel je het metamateriaal voor als een gigantisch, zich herhalend behangpatroon. In plaats van de hele muur te analyseren, kijk je gewoon naar één enkel vierkant van het behang (de RVE, of Representatief Volume-element).

De truc van de auteurs is om aan te nemen dat als je de hele lange buis uitrekt of draait, dat kleine vierkante stukje ook uitrekt of draait, maar op een zeer specifieke, spiraalvormige manier. Ze noemen dit een "helixvormige" vervorming. Stel je een veer voor die je uitrekt; de spoelen worden niet alleen langer, ze draaien ook een beetje. De auteurs realiseerden zich dat door dit kleine vierkante stukje te dwingen om die exacte spiraalbeweging na te bootsen, ze kunnen uitrekenen hoe de gehele lange buis zou reageren zonder het hele ding te simuleren.

2. De "Perfect Flexibele" Stokjes

De meeste computermodellen behandelen stokjes als stijf en onveranderlijk, zoals een liniaal van staal. Maar in het echt, vooral bij deze tiny metamaterialen, kunnen de stokjes tegelijkertijd buigen, rekken en schuiven (zijwaarts glijden), zelfs wanneer de vervorming enorm is.

De auteurs gebruiken een speciaal wiskundig model dat een "Speciale Cosserat-stok" wordt genoemd.

• Analogie: Stel je een stuk gekookte spaghetti voor. Het kan buigen, het kan een beetje rekken en het kan draaien. Stel je nu voor dat die spaghetti is gemaakt van een materiaal dat al deze dingen perfect en nauwkeurig kan doen, zelfs als je het in een cirkel buigt of tot het dubbele van zijn lengte uitrekt. Dat is wat hun model doet. Het benadert niet alleen; het vangt de exacte geometrie van de buiging en de draaiing.

3. De "Dansvloer"-Regels (Randvoorwaarden)

Om het kleine vierkante stukje te laten gedragen alsof het deel uitmaakt van een gigantische, zich herhalende buis, moesten de auteurs een reeks regels bedenken voor hoe de randen van dat vierkant met elkaar communiceren.

• Het Probleem: Als je een stuk van een wenteltrap afsnijdt, sluit de bovenrand niet perfect aan op de onderrand.
• De Oplossing: Ze creëerden een "helixvormige randvoorwaarde". Stel je voor dat de linkerkant van je kleine vierkante stukje hand in hand gaat met de rechterkant, maar de rechterkant is iets gedraaid en verschoven, net als treden op een wenteltrap.
• De Innovatie: Eerdere methoden konden alleen kleine, zachte bewegingen aan. De nieuwe regel van de auteurs werkt zelfs als de buis tot een pretzel is gedraaid of uitgerekt tot hij dun is als een draad. Het is "geometrisch exact", wat betekent dat het nooit nauwkeurigheid verliest, hoe wild de vorm ook wordt.

4. De "Gewrichten" en "Lijm"

Binnen dat kleine vierkante stukje zijn de stokjes verbonden bij gewrichten.

• Stijve Gewrichten: Sommige gewrichten zijn als supersterke lijm; de stokjes kunnen zich niet ten opzichte van elkaar bewegen op het verbindingspunt.
• De Wiskunde: De auteurs hebben een systeem opgezet waarbij de computer de beste positie voor elke stok in dat kleine vierkant berekent, ervoor zorgend dat de gewrichten verbonden blijven en de "wenteltrap"-regels worden gevolgd, terwijl er zo min mogelijk energie wordt gebruikt.

5. Wat Ze Vonden (De Resultaten)

Zodra ze de wiskunde voor het kleine stukje hadden opgelost, konden ze voorspellen hoe de hele buis zou reageren. Ze testten dit met verschillende vormen:

• Het Kruis- en Vierkantvormige Profiel: Ze begonnen met eenvoudige vormen (zoals een plus-teken of een vierkant gemaakt van stokjes) om te bewijzen dat hun wiskunde werkte. Ze ontdekten dat als de tiny stokjes dik en kort zijn, het veel uitmaakt of ze kunnen rekken of schuiven. Als ze erg dun en lang zijn, werkt de oude, eenvoudigere wiskunde prima.
• De Helixvormige (Veer) Stokjes: Ze keken naar een vierkant gemaakt van stokjes die al gebogen zijn als veren (helixen).
  • De "J-vormige" Rek: Toen ze dit materiaal trokken, was het eerst zacht (zoals een veer die ontspoelt) maar werd het erg stijf naarmate het zich rechtte. Dit creëert een "J-vormige" curve. Dit is precies hoe biologische weefsels (zoals spieren) zich gedragen, daarom noemen de auteurs dat deze kunnen worden gebruikt voor kunstspieren.
  • De Verzachtende Buiging: Toen ze het bogen, werd het materiaal zachter naarmate ze meer bogen. Dit gebeurde omdat het verbindende veer-stokje begon te draaien uit het vlak, werkend als een scharnier.
• De Auxetische Buis: Ze modelleerden een holle buis die breder wordt als je er aan trekt (zoals een honingraat).
  • Ze toonden aan dat je door de hoek van de stokjes te veranderen, de buis kunt afstemmen om zeer flexibel zijwaarts te zijn (goed voor buigen) maar zeer stijf tegen het worden samengeperst (goed voor het openhouden van bloedvaten).
  • Ze merkten op dat deze structuren kunnen worden afgestemd om "verkorting" te voorkomen (korter worden bij uitbreiding), wat een veelvoorkomend probleem is bij cardiovasculaire stents (netbuisjes die worden gebruikt om slagaders open te houden).

Samenvatting

De auteurs bouwden een "universele vertaler" voor metamaterialen. Ze creëerden een methode die een complex, 3D-netwerk van tiny stokjes vertaalt naar een eenvoudige, gladde wiskundige beschrijving van een enkele stok. Dit stelt ingenieurs in staat om complexe, flexibele materialen te ontwerpen voor dingen zoals robotarmen, kunstspieren en medische stents door de tiny interne patronen aan te passen, wetende precies hoe het eindproduct zal buigen en rekken, zonder dat ze voor elke ontwerpverandering een supercomputersimulatie hoeven te draaien.

Technische samenvatting

Technische Samenvatting: Homogenisatie van staafachtige metamaterialen als een speciale Cosserat-staaf

Probleemstelling
Staafachtige metamaterialen zijn structuren die worden gevormd door periodiek microstructurele eenheden (typisch netwerken van staven) in één richting te assembleren, wat resulteert in een macroscopische geometrie die aanzienlijk langer is in één dimensie dan in de andere twee. Deze structuren worden toegepast in toepassingen variërend van cardiovasculaire stents en kunstspieren tot robotische actuatoren. Hoewel directe numerieke simulatie in volledige schaal van dergelijke heterogene structuren computationeel onhaalbaar is, ondervinden standaard homogenisatietechnieken vaak beperkingen. Specifiek zijn homogenisatieschema's van de eerste orde ontoereikend voor het vastleggen van grote vervormingen met aanzienlijke rekgradiënten (zoals buigen en wringen) die kenmerkend zijn voor staafachtige structuren. Bovendien slagen standaard driedimensionale homogenisatiebenaderingen er niet in om grootte-effecten in aanmerking te nemen wanneer de scheiding van lengteschalen niet geldt in de transversale richtingen. Er is behoefte aan een raamwerk dat de niet-lineaire constitutieve respons van deze metamaterialen onder willekeurige grote vervormingen nauwkeurig kan modelleren, waarbij rekening wordt gehouden met eindige schuiving, rek, buiging en wringing op zowel micro- als macroschaal.

Methodologie
De auteurs stellen een computationeel homogenisatieschema voor dat staafachtige metamaterialen op macroschaal modelleert als een "speciale Cosserat-staaf", terwijl de microstructurele eenheid (Representative Volume Element of RVE) wordt gemodelleerd als een netwerk van geometrisch exacte Cosserat-staven.

1. Geometrisch exacte kinematica: Zowel de macroscopische effectieve staaf als de microscopische constituentstaven worden gemodelleerd met behulp van de geometrisch exacte theorie van de speciale Cosserat-staaf. Deze formulering houdt rekening met eindige buiging, wringing, schuiving en axiale rek, en vermijdt de aannames van kleine rekken van klassieke balktheorieën.
2. Extensie van de Helicale Cauchy-Born (HCB) regel: De methode breidt de Helicale Cauchy-Born regel uit. Er wordt aangenomen dat het metamateriaal op macroschaal wordt blootgesteld aan een uniform rekveld langs zijn booglengte. Deze aanname maakt het mogelijk het volledige structurele probleem te reduceren tot dat van een enkele RVE.
3. Helaal periodieke randvoorwaarden: Een belangrijke innovatie is de afleiding van een geometrisch exacte, niet-lineaire afbeelding om helaal periodieke randvoorwaarden toe te passen op de RVE. Deze afbeelding relateert de translatie- en rotatievrijheidsgraden van de rechtergrens van de RVE aan die van de linkergrens op basis van de macroscopische rekparameters ($v_0, k_0$). In tegenstelling tot eerdere lineaire overgangsafbeeldingen, blijft deze formulering geldig voor willekeurig grote macroscopische rekken.
4. Gedwongen energie-minimalisatie: Het RVE-probleem wordt geformuleerd als een gedwongen minimalisatie van elastische energie. De beperkingen omvatten:
   • Interne verbindingbeperkingen: Het afdwingen van stijve connectiviteit tussen staven binnen de RVE met behulp van een master-slave-benadering.
   • Helaale randbeperkingen: Het afdwingen van de periodiciteit die voortvloeit uit de HCB-regel.
   • Globale beperkingen: Het beperken van de translatie en rotatie van het RVE als een star lichaam om een unieke oplossing te garanderen en de macroscopische rekparameters tijdens de minimalisatie vast te houden.
5. Discretisatie en oplossing: De RVE-energie wordt gediscretiseerd met behulp van "helaale elementen", waarbij wordt aangenomen dat elk segment een uniform rekveld ondergaat, waardoor het segment zelf helicaal wordt. De resulterende niet-lineaire zwakke vorm wordt opgelost met behulp van een eigen code op het deal.II-platform.
6. Afleiding van macroscopische eigenschappen: Gesloten uitdrukkingen voor de geïntegreerde interne contactkracht, het moment en de stijfheidstensor van de homogeniseerde staaf worden afgeleid door de geminimaliseerde RVE-energiedichtheid te differentiëren naar de macroscopische rekparameters.

Belangrijkste bijdragen
Het artikel identificeert drie primaire nieuwe bijdragen:

1. Geometrische exactheid op twee schalen: Zowel microscopische als macroscopische staven worden gemodelleerd als geometrisch exacte speciale Cosserat-staven, waarbij eindige schuiving en rek worden vastgelegd naast buiging en wringing op beide schalen.
2. Geometrisch exacte macro-naar-micro afbeelding: Een volledig niet-lineaire afbeelding wordt voorgesteld voor het toepassen van helaal periodieke randvoorwaarden op zowel translatie- als rotatievrijheidsgraden. Deze afbeelding is zelfs geldig bij grote macroscopische rekken, waarmee een beperking van eerdere lineaire overgangsafbeeldingen wordt aangepakt.
3. Willekeurig microscopisch constitutief gedrag: Het raamwerk staat willekeurig constitutief gedrag toe van de constituentstaven binnen de RVE, en gaat verder dan lineair elastische aannames.

Resultaten en numerieke voorbeelden
De auteurs valideren het schema met behulp van verschillende numerieke voorbeelden met variërende RVE-complexiteit:

• Kruis- en vierkante RVE's: Eenvoudige planaire RVE's (kruis- en vierkante vormen) werden gebruikt om de resultaten te valideren tegen bestaande analytische formules en literatuur. De studie toonde aan dat het opnemen van schuiving en extensie in microscopische staven de macroscopische stijfheid significant beïnvloedt, met name voor staven met grotere slankheidsverhoudingen. De resultaten kwamen overeen met analytische limieten voor Kirchhoff-staven (onrekbaar/onbuigbaar).
• Helaale staaf-RVE's: Een complexere vierkante RVE bestaande uit helicaal samengestelde staven werd geanalyseerd. De resultaten toonden een J-vormige macroscopische spannings-rekrespons, die overgaat van buigingsgedomineerd naar rek-gedomineerd gedrag. De studie benadrukte hoe parameters zoals het aantal helicale windingen en de steek deze respons kunnen afstemmen. Buigingsanalyse onthulde een verzachtend gedrag dat wordt geactiveerd door het uit het vlak buigen van de kruisende helicale staven.
• Auxetische buisvormige structuren: Het schema werd toegepast op auxetische buisvormige metamaterialen (relevant voor stents). Het model slaagde erin de afhankelijkheid van de macroscopische buigstijfheid, de Poisson-ratio (krommingsgedrag) en de radiale stijfheid van de ontwerphoek $\theta_a$ vast te leggen. Het toonde ook aan dat het Bazier-effect (ovalisatie van de dwarsdoorsnede) onder grote buiging kan worden vastgelegd.
• FE2-validatie: Een vergelijking tussen het voorgestelde homogenisatieschema (FE2-benadering) en simulaties in volledige schaal van een kruis-RVE-metamateriaal toonde goede overeenkomst, waarmee de nauwkeurigheid van de methode werd gevalideerd, zelfs onder grote wringmomenten waarbij lineaire constitutieve modellen faalden.

Betekenis en claims
Het artikel claimt een robuust raamwerk te bieden voor de homogenisatie van staafachtige metamaterialen dat de beperkingen van homogenisatieschema's van de eerste orde en lineaire homogenisatie overwint. Door gebruik te maken van een geometrisch exacte theorie van de speciale Cosserat-staaf op beide schalen en het afleiden van een niet-lineaire macro-naar-micro overgangsafbeelding, voorspelt de methode nauwkeurig de niet-lineaire constitutieve respons van deze structuren onder grote vervormingen.

De auteurs benadrukken dat hun benadering algemeen is en kan worden uitgebreid naar andere RVE-types (bijvoorbeeld poreuze vaste stoffen, origami-structuren) en complexe microscopische fenomenen kan incorporeren, zoals elastoplasticiteit, viscoelasticiteit en multifysische koppeling (thermo-, chemo- of elektro-elasticiteit). Het werk dient als fundament voor het ontwerpen van op maat gemaakte metamaterialen met specifieke macroscopische eigenschappen, zoals kunstspieren of uitschuifbare stents, door microstructurele parameters af te stemmen. De auteurs merken bescheiden op dat hoewel het huidige werk ervan uitgaat dat de RVE een stavennetwerk is, de afgeleide randvoorwaarden kunnen worden aangepast voor andere microstructurele representaties.