Model-Informed Joint Material-Structural Optimization of Hard-Magnetic Soft Materials

Technische Samenvatting: Model-geïnformeerde Gezamenlijke Materiaal-Structurele Optimalisatie van Hard-Magnetische Zachte Materialen

Probleemstelling
Hard-magnetische zachte materialen (hMSM's) zijn composieten bestaande uit hard-magnetische deeltjes die zijn ingebed in een compliant elastomeer-matrix. Ze maken contactloze, veldgestuurde vervorming mogelijk die geschikt is voor zachte robotica en adaptieve structuren. De voorspellende ontwerping van hMSM's wordt echter geconfronteerd met twee onderling verbonden uitdagingen:

1. Constitutief Modellering: Nauwkeurige voorspelling vereist effectieve structuur-eigenschap-relaties die beschrijven hoe het volumefractie van magnetische deeltjes de effectieve elastische moduli modificeert. Bestaande literatuur biedt diverse modellen (rigid-inclusion, self-consistent, constrained-kinematics) die significant verschillende stijfheidswaarden kunnen voorspellen voor dezelfde deeltjesconcentratie. Bovendien kan de keuze van de rek-energie-dichtheidsfunctie de numerieke robuustheid en de voorspelde vervorming beïnvloeden.
2. Gezamenlijke Optimalisatie: Optimale ontwerpen vereisen de simultane controle van de structurele dichtheid, de verdeling van magnetische deeltjes en de richting van de remanente magnetisatie. Eerdere topologie-optimalisatiebenaderingen voor hMSM's behandelen de aanwezigheid van magnetisch materiaal vaak als een discrete keuze (gemagnetiseerd versus niet-gemagnetiseerd) in plaats van een continu volumefractieveld dat zowel de stijfheid als de actuatiekracht moduleert. Bovendien is het fysiek lastig om deze velden te koppelen omdat magnetische variabelen alleen betekenis hebben waar structureel materiaal aanwezig is.

Methodologie
De auteurs ontwikkelen een verenigd framework dat zowel de keuze van het constitutieve model als de gezamenlijke optimalisatie aanpakt.

• Verenigd Constitutief Framework: De auteurs herformuleren klassieke rigid-inclusion relaties (Guth, Mooney, Kerner), Hill's self-consistent theorie (gespecialiseerd voor rigide inclusies en bijna-incompressibele matrices), en constrained-kinematics modellen (Lopez-Pamies, Alshammari) in een gemeenschappelijke effectieve afschuifmodulusvorm, $G(\phi)$, waarbij $\phi$ de volumefractie van de deeltjes is. Dit levert zeven kandidaat-afschuifmodulus-relaties op. Deze worden gecombineerd met drie rek-energie-dichtheidsfuncties (twee compressibele neo-Hooke vormen, NH1 en NH2, en het St. Venant–Kirchhoff model, SVK) om 21 kandidaat-constitutieve modellen te creëren.
• Model Sensitiviteitsanalyse: Numerieke benchmarks (uitkragende balken, wielen, grijpers en 3D-balken) worden gebruikt om de gevoeligheid van de voorspelde vervorming te beoordelen ten opzichte van de keuze van de rek-energie-dichtheidsfunctie versus de keuze van $G(\phi)$.
• Constitutieve Model Selectie: Experimentele uniaxiale spanning–rekgegevens van Garcia-Gonzalez et al. (2021) worden gebruikt om de meest representatieve $G(\phi)$ relatie te selecteren. Twee metrieken worden gehanteerd: (1) relatieve fout ten opzichte van geëxtraheerde afschuifmodulus-waarden, en (2) relatieve wortel-gemiddelde kwadratische fout (rRMSE) ten opzichte van volledige ware spanning–rek curves.
• Gezamenlijke Materiaal–Structurele Optimalisatie Formulering: Een topologie-optimalisatie framework wordt ontwikkeld waarbij drie ontwerpvelden worden geparametriseerd: structurele dichtheid ($\rho$), de volumefractie van magnetische deeltjes ($\phi$) en de richting van de remanente magnetisatie ($\theta_{Br}$).
  • Koppelingsmechanisme: Om fysieke consistentie te waarborgen, worden de magnetische energie en de effectieve moduli geïnterpoleerd zodat magnetische bijdragen verdwijnen in holle gebieden ($\rho \approx 0$). De volumefractie van magnetische deeltjes $\phi$ fungeert als een continu veld dat zowel de magnetische interactie-energie als de effectieve afschuifmodulus schaalt via de geselecteerde $G(\phi)$ relatie.
  • Implementatie: Het framework maakt gebruik van de FEniTop bibliotheek en FEniCSx voor de finite-element implementatie, waarbij een SIMP-type penalisatie wordt toegepast voor dichtheid, Helmholtz filtering voor alle velden, en een adjoint-gebaseerde gradiëntberekening voor de Method of Moving Asymptotes (MMA) optimizer.

Belangrijkste Resultaten

• Constitutieve Sensitiviteit: De rek-energie-dichtheidsfunctie (NH1, NH2, SVK) heeft een relatief kleine invloed op de voorspelde vervorming voor de overwogen actuatieproblemen. In contrast hiermee heeft de keuze van de effectieve afschuifmodulus-relatie $G(\phi)$ een significante invloed op de vervorming, vooral wanneer magnetisch materiaal overlapt met sterk vervormende regio's. Wanneer magnetisch materiaal gelokaliseerd is buiten de primaire buigregio's (bijv. tip-gemagnetiseerde balken), is de modelafhankelijkheid verwaarloosbaar.
• Model Selectie: Op basis van experimentele gegevens vertoonde de Mooney-relatie voor $G(\phi)$ de beste algemene overeenkomst met zowel de geëxtraheerde afschuifmodulus-waarden als de volledige spanning–rek responsen. De NH2 rek-energie-dichtheidsfunctie werd geselecteerd vanwege de numerieke robuustheid.
• Optimalisatie Prestaties: Het gezamenlijke materiaal–structurele optimalisatie framework produceerde succesvol niet-intuïtieve ontwerpen voor drie verschillende gevallen:
  1. Rotationele Actuator: Het optimaliseren van $\phi$ en $\theta_{Br}$ in een vaste wielgeometrie verhoogde de tangentiële verplaatsing met 59,7% vergeleken met een intuïtief uniform ontwerp.
  2. Translationele Actuator: Het optimaliseren van $\phi$ en $\theta_{Br}$ in een vaste schaar-achtige geometrie bereikte een doelgerichte horizontale beweging terwijl verticale parasitaire verplaatsing werd onderdrukt.
  3. Herstellende Balk: Gelijktijdige optimalisatie van $\rho$, $\phi$ en $\theta_{Br}$ creëerde een balk die in staat is zijn vorm te herstellen onder tegenovergestelde mechanische belastingen met behulp van magnetische actuatie. Het ontwerp beperkte het magnetische materiaal van nature tot de structureel bezette regio's.

Significantie en Claims
Het artikel claimt een systematische, model-geïnformeerde route te bieden voor het ontwerpen van hMSM's. De primaire bijdragen zijn:

1. Een verenigd effectief afschuifmodulus-framework dat directe vergelijking van diverse structuur–eigenschap-relaties mogelijk maakt.
2. Bewijs dat effectieve afschuifmodulus-relaties cruciaal zijn voor voorspellende nauwkeurigheid wanneer magnetische en mechanische velden overlappen, terwijl rek-energie-vormen minder gevoelig zijn voor de geteste gevallen.
3. Een gezamenlijke materiaal–structurele optimalisatie formulering die de volumefractie van magnetische deeltjes behandelt als een continue ontwerpvariabele gekoppeld aan structurele dichtheid en magnetiseringsrichting, waarbij het fysieke koppelingsprobleem in gebieden met lage dichtheid wordt opgelost door middel van energie-interpolatie.
4. Demonstratie dat het gelijktijdig optimaliseren van de materiaalsamenstelling en de structurele topologie superieure prestaties oplevert ten opzichte van vast-geometrisch magnetisch ontwerp.

De auteurs merken op dat hoewel het framework gedemonstreerd wordt op hMSM's, de computationele structuur toepasbaar is op andere veld-responsieve elastomeren (bijv. diëlektrische, magnetorheologische of vormgeheugen-composieten) door de relevante energietermen te substitueren. De implementatie is beschikbaar via de open-source CEADpx/top_optim repository om reproduceerbaarheid en uitbreiding te ondersteunen.