Derivation of variational membrane models in the context of anisotropic nonlocal hyperelasticity

1. Probleemstelling en Motivatie

Het artikel richt zich op de wiskundige afleiding van effectieve membraanmodellen voor dunne structuren (zoals films of membranen) binnen het kader van niet-lokale hyperelasticiteit.

• Niet-lokale continuümmechanica: In tegenstelling tot klassieke elasticiteit, die lokale afgeleiden (gradienten) gebruikt en faalt bij singulariteiten zoals breuken, maakt niet-lokale mechanica (zoals peridynamica) gebruik van interacties tussen deeltjes binnen een bepaald bereik (horizon). De energie wordt bepaald door paarsgewijze krachten binnen dit bereik.
• De uitdaging: Bestaande theorieën voor niet-lokale gradienten gaan vaak uit van isotrope interactiegebieden (bolvormig met een vaste straal). Wanneer men echter een driedimensionaal lichaam afplat tot een dunne film (dikte $\varepsilon \to 0$), verandert de geometrie van het interactiegebied fundamenteel. Een oorspronkelijk bolvormig interactiegebied wordt in de gerescaleerde coördinaten anisotroop (ellipsoïdaal), waarbij de straal in de dikte-richting anders schaalt dan in het vlak.
• Gaps in de literatuur: Er is een gebrek aan rigoureuze variatie-theoretische afleidingen voor dimensiereductie in niet-lokale modellen, vooral voor modellen met een eindige horizon (finite horizon). Bestaande werken focussen vaak op heuristische benaderingen of oneindige horizons (fractionele gradienten).

2. Methodologie

De auteurs ontwikkelen een nieuwe wiskundige framework om dit probleem aan te pakken, bestaande uit drie hoofdfasen:

A. Theorie van Anisotrope Niet-Lokale Gradienten (Sectie 2)

Om de effecten van het afvlakken van het domein te modelleren, definiëren de auteurs een nieuwe klasse van niet-lokale gradienten met anisotrope horizons.

• Definitie: In plaats van één isotrope straal $\delta$, gebruiken ze een vector $\delta = (\bar{\delta}, \delta_3)$, waarbij $\bar{\delta}$ de straal in het vlak ($x_1, x_2$) en $\delta_3$ de straal in de dikterichting ($x_3$) controleert. Deze kunnen onafhankelijk van elkaar naar nul gaan.
• Interactiegebied: Het interactiegebied is een ellipsoïde met hoofdradii $\bar{\delta}, \bar{\delta}, \delta_3$.
• Kernconcept: Ze introduceren een translatiemechanisme (Propositie 2.15). Dit is een isomorfisme dat de niet-lokale ruimte $H^{\rho,p}_\delta$ koppelt aan de klassieke Sobolev-ruimte $W^{1,p}$. Dit stelt hen in staat om bewijzen in het complexe niet-lokale kader terug te brengen naar het welbekende lokale geval.
• Technische uitdagingen: Ze bewijzen dichtheidsresultaten voor gladde functies en een Leibniz-regel die uniform geldig blijft, zelfs wanneer de horizoncomponenten naar nul convergeren (degenererende gevallen).

B. Dimensiereductie en $\Gamma$-convergentie (Sectie 3)

De auteurs analyseren de asymptotische limiet van de energie voor een dunne film met dikte $\varepsilon$.

• Rescaling: Het fysieke domein $\Omega_\varepsilon = \omega \times (0, \varepsilon)$ wordt gerescaleerd naar een referentiedomein $\Omega = \omega \times (0, 1)$.
• Gedrag van de horizon: De fysieke horizon $\delta(\varepsilon)$ wordt gerescaleerd naar $\delta_\varepsilon = (\bar{\delta}, \delta_3/\varepsilon)$. Afhankelijk van hoe $\delta_3$ schaalt ten opzichte van $\varepsilon$, ontstaan verschillende regimes:
  1. Anisotrope schaling: Als $\delta_3 \sim \varepsilon$, blijft de horizon in de gerescaleerde ruimte eindig en anisotroop.
  2. Isotrope schaling: Als $\delta_3 \ll \varepsilon$ (of $\delta \to 0$), convergeert het model naar het klassieke lokale model.
• $\Gamma$-convergentie: Ze bewijzen dat de geregelde energie $I_\varepsilon$ $\Gamma$-convergeert naar een limietfunctioanale $I$.

C. Stabilisatie

Omdat niet-lokale gradienten met eindige horizon niet-triviale functies hebben met een nul-gradient (die geen energie kosten), is de energie niet coercief. De auteurs introduceren een stabilisatieterm $S_\varepsilon$ die grote vervormingen met nul-gradient straft en de convergentie naar functies die constant zijn in de $x_3$-richting afdwingt.

3. Belangrijkste Resultaten

1. Unificatie van Modellen: De afgeleide theorie vormt een brug tussen volledig niet-lokale, gedeeltelijk niet-lokale en puur lokale modellen. Het is de eerste rigoureuze afleiding voor niet-lokale hyperelasticiteit met een eindige horizon in dunne structuren.
2. De Limietfunctioanale:
   • De $\Gamma$-limiet $I(u)$ heeft de structuur van een klassiek membraanmodel, maar met een niet-lokale in-vlak gradient $\bar{D}_{\rho, \bar{\delta}_0}$.
   • De energie dichtheid wordt de quasi-convexe omhulling ($W^{qc}$) van de oorspronkelijke energie, geminimaliseerd over de uitval-richting.
   • De limietfunctie $u$ is onafhankelijk van de dikte-coördinaat $x_3$ (binnen de ruimte van functies met nul-gradient).
3. Twee Specifieke Regimes:
   • Regime 1 (Anisotroop): Als de fysieke horizon in de verticale richting evenredig is met de dikte ($\delta_3 \sim \varepsilon$), resulteert dit in een niet-lokaal membraanmodel met een effectieve horizon $\bar{\delta}_0 > 0$.
   • Regime 2 (Lokaal): Als de horizon in alle richtingen schaalt met $\varepsilon$ ($\delta \sim \varepsilon$), dan is $\bar{\delta}_0 = 0$. Het model reduceert exact tot het klassieke lokale hyperelastische membraanmodel.
4. Compactheid en Convergentie van Minimizers:
   • Zonder stabilisatie: Convergentie geldt alleen "op tot een representant" in de ruimte van functies met nul-gradient.
   • Met stabilisatie: Er is sterke compactheid en convergentie van de minimizers zelf, zonder de noodzaak om representanten te kiezen.

4. Significatie en Impact

• Wiskundige Strenge: Het artikel levert een fundamentele bijdrage aan de variatie-analyse van niet-lokale problemen. Het oplossen van het probleem van anisotrope horizons en het bewijzen van uniforme schattingsresultaten (onafhankelijk van de horizon-grootte) is een belangrijke doorbraak.
• Fysische Relevantie: Het biedt een theoretische basis voor het modelleren van dunne materialen (zoals membranen of coatings) waarbij niet-lokale effecten (zoals microstructuur of langereafstandsinteracties) belangrijk zijn. Het verduidelijkt hoe de schaal van de interactie ten opzichte van de dikte van het materiaal bepaalt of het materiaal zich lokaal of niet-lokaal gedraagt.
• Toepassingsgebied: De resultaten zijn relevant voor de ontwikkeling van numerieke methoden (zoals peridynamica) voor dunne structuren en voor het begrijpen van overgangen tussen schaalniveaus in materialenwetenschap.

Kortom, dit werk generaliseert de klassieke theorie van dimensiereductie (zoals Kirchhoff-plate theorieën) naar het domein van niet-lokale hyperelasticiteit, waarbij het de complexiteit van anisotrope interactiegebieden oplost en een brug slaat tussen niet-lokale en lokale continuümmechanica.