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

1. Il Problema e il Contesto

Il lavoro si concentra sulla riduzione dimensionale variazionale di modelli di iperelasticità non locale, specificamente nell'ambito della peridinamica (una teoria non locale che gestisce naturalmente discontinuità come fratture).
L'obiettivo principale è derivare un modello efficace per membrane sottili (strutture tridimensionali con spessore $\varepsilon \to 0$) partendo da un modello tridimensionale non locale.

Le sfide principali affrontate sono:

• Anisotropia indotta dal ridimensionamento: Quando si scala un dominio sottile (riducendo lo spessore), i range di interazione non locali, che sono naturalmente isotropi (sferici) nel dominio fisico, diventano anisotropi (ellissoidali) nel dominio ridimensionato.
• Limiti dei modelli esistenti: La letteratura esistente sui gradienti non locali si basa spesso su kernel con supporto sferico (orizzonte isotropo). Questi modelli non coprono il caso in cui i raggi di interazione possono annullarsi indipendentemente nelle diverse direzioni (es. interazione nulla nella direzione dello spessore).
• Assenza di risultati rigorosi: Mentre esistono derivazioni formali o basate su simulazioni numeriche per membrane peridinamiche, manca una teoria variazionale rigorosa (basata sulla $\Gamma$-convergenza) per modelli non locali a orizzonte finito in regime di film sottile.

2. Metodologia

Gli autori sviluppano un approccio matematico rigoroso basato sulla $\Gamma$-convergenza e su nuove strutture analitiche per i gradienti non locali.

A. Teoria dei Gradienti Non Locali Anisotropi

Per gestire la geometria delle membrane sottili, gli autori introducono una teoria unificata per i gradienti non locali con orizzonte ellissoidale:

• Definizione del Gradiente: Sostituiscono l'orizzonte isotropo $\delta$ con una coppia di orizzonti $\delta = (\bar{\delta}, \delta_3)$, dove $\bar{\delta}$ controlla il raggio di interazione nel piano e $\delta_3$ quello fuori dal piano.
• Kernel Anisotropi: Il gradiente non locale $D_{\rho, \delta}$ è definito tramite un'ellissoide di interazione con assi principali $\bar{\delta}, \bar{\delta}, \delta_3$. Questo permette di interpolare tra modelli completamente non locali, parzialmente non locali e puramente locali (quando i raggi tendono a zero).
• Spazi Funzionali: Vengono definiti spazi di Sobolev non locali $H^{\rho, p}_{\delta}$ adattati a queste geometrie anisotrope.

B. Il Meccanismo di Traduzione (Translation Mechanism)

Un risultato tecnico fondamentale è l'estensione del "meccanismo di traduzione" (introdotto in lavori precedenti per orizzonti isotropi) al caso anisotropo e degenerato:

• Viene costruita un'isomorfismo tra lo spazio non locale $H^{\rho, p}_{\delta}$ e lo spazio locale classico $W^{1,p}$.
• Questo meccanismo permette di trasformare problemi non locali in problemi locali (e viceversa) tramite operatori di convoluzione ($Q_{\rho, \delta}$ e il suo inverso $P_{\rho, \delta}$), consentendo di utilizzare strumenti noti dell'analisi locale per dimostrare proprietà come disuguaglianze di Poincaré e compattezza, mantenendo stime uniformi rispetto al parametro $\delta$.

C. Analisi Asintotica ($\Gamma$-convergenza)

Il modello di energia elastica immagazzinata viene riscritto su un dominio fisso $\Omega = \omega \times (0,1)$ tramite un ridimensionamento lineare nella direzione fuori dal piano.

• Si studia il limite asintotico della famiglia di funzionali energetici $I_\varepsilon$ quando $\varepsilon \to 0$ e gli orizzonti $\delta(\varepsilon)$ convergono a un limite $\delta_0$.
• Viene affrontata la problematica delle funzioni con gradiente non locale nullo (che non sono costanti nel caso non locale), introducendo un termine di stabilizzazione per garantire la compattezza forte.

3. Contributi Chiave

1. Sviluppo di una teoria anisotropa: Creazione di un quadro teorico completo per i gradienti non locali con orizzonti ellissoidali, che include casi degeneri (raggi nulli) e generalizza i risultati esistenti per orizzonti sferici.
2. Dimostrazione della $\Gamma$-convergenza: Prima derivazione rigorosa di modelli di membrana per materiali iperelastici non locali a orizzonte finito.
3. Unificazione dei regimi: Il modello limite ottenuto copre automaticamente sia il caso non locale (se l'orizzonte rimane finito) sia il caso locale classico (se l'orizzonte tende a zero), mostrando come la teoria locale sia un caso particolare di quella non locale.
4. Riduzione dimensionale esplicita: Dimostrazione che il limite è un modello di membrana bidimensionale, dove l'energia dipende solo dalle variabili nel piano, anche partendo da un modello 3D non locale.

4. Risultati Principali

Il Teorema di $\Gamma$-convergenza (Teorema 3.1)

La famiglia di energie $I_\varepsilon$ $\Gamma$-converge (rispetto alla topologia forte in $L^p$) a un funzionale limite $I$:
$$I(u) = \begin{cases} \int_{\omega} W^{qc}(\bar{D}_{\rho, \bar{\delta}_0} \bar{u}) \, d\bar{x} & \text{se } u \in H^{\rho, p}_{\delta_0}(\Omega), (D_{\rho, \delta_0}u)|_{\Omega}e_3 = 0 \\ \infty & \text{altrimenti} \end{cases}$$
Dove:

• $\bar{D}_{\rho, \bar{\delta}_0}$ è il gradiente non locale bidimensionale (solo le prime due colonne).
• $W^{qc}$ è l'involucro quasiconvesso della densità di energia.
• La condizione $(D_{\rho, \delta_0}u)e_3 = 0$ impone che il gradiente non locale nella direzione dello spessore sia nullo, il che implica (a meno di funzioni nel nucleo) che la deformazione non dipende dalla coordinata verticale.

Casi Particolari di Interesse

Gli autori analizzano due casi di scaling fisico:

1. Orizzonte anisotropo scalato: Se $\delta(\varepsilon) = (1, \varepsilon)$, il limite è un modello non locale puro con orizzonte finito nel piano.
2. Orizzonte isotropo che scala con $\varepsilon$: Se $\delta(\varepsilon) = (\varepsilon, \varepsilon)$, il limite è un modello locale classico (iperelasticità di membrana), poiché l'orizzonte effettivo nel piano tende a zero.

Stabilizzazione e Compattezza

Poiché nel caso non locale esistono deformazioni non costanti con gradiente nullo (che rendono l'energia invariante per traslazioni in un sottospazio infinito-dimensionale), gli autori introducono un termine di penalizzazione $S_\varepsilon$.

• Questo termine garantisce la compattezza forte e la convergenza dei minimizzanti senza dover sottrarre arbitrariamente funzioni di "nucleo".
• Il risultato finale (Teorema 3.4) mostra la $\Gamma$-convergenza del funzionale stabilizzato verso un modello di membrana con un termine di penalizzazione aggiuntivo.

5. Significato e Impatto

Questo lavoro è significativo per diversi motivi:

• Rigore Matematico: Colma il divario tra le derivazioni formali/numeriche esistenti e una teoria variazionale rigorosa per le strutture sottili non locali.
• Generalità: La teoria sviluppata non è limitata alla peridinamica classica ma si applica a una vasta classe di modelli di iperelasticità non locale con kernel a supporto compatto.
• Flessibilità di Modellazione: Fornisce un framework unificato che permette di studiare transizioni continue tra comportamenti locali e non locali in strutture sottili, a seconda di come l'orizzonte di interazione scala rispetto allo spessore del materiale.
• Applicabilità: I risultati sono cruciali per la simulazione numerica e la progettazione di materiali avanzati (come elastomeri cristalli liquidi o materiali con microstruttura) dove gli effetti non locali e la geometria sottile giocano un ruolo fondamentale.

In sintesi, il paper stabilisce le fondamenta analitiche per derivare modelli di membrana efficaci a partire da teorie non locali complesse, gestendo con successo le sfide geometriche e analitiche introdotte dall'anisotropia indotta dal ridimensionamento.