Nonlinear Incompressible Shear Wave Models in Hyperelasticity and Viscoelasticity Frameworks, with Applications to Love Waves

Questo lavoro presenta equazioni generali per onde di taglio non lineari in materiali iper-viscoelastici incompressibili, applicandole all'analisi numerica delle onde di Love interfacciali e superficiali, dimostrando che, nel caso pienamente non lineare, la velocità d'onda tende asintoticamente alla maggiore delle velocità dei materiali costituenti pur rimanendo inizialmente compresa tra i limiti lineari.

L'essenziale

Immagina di essere un geologo o un ingegnere che studia come le onde si muovono attraverso la Terra, o forse attraverso un tessuto biologico. Fino a poco tempo fa, la maggior parte dei modelli usati per descrivere queste onde (come le onde sismiche chiamate "onde di Love") si basava su una semplificazione: assumevano che i materiali si comportassero come molle perfette. Se le allunghi un po', tornano indietro; se le allunghi troppo, si rompono. È un modello lineare, semplice e comodo.

Ma la realtà è più complessa. Materiali come la gomma, le schiume, i polimeri e i tessuti del nostro corpo non sono molle perfette. Quando vengono sollecitati fortemente, si comportano in modo non lineare: si deformano in modo strano, cambiano rigidità e dissipano energia (come se avessero una sorta di "attrito interno").

Questo articolo scientifico, scritto da tre ricercatori dell'Università del Saskatchewan, si propone di creare una nuova mappa matematica per descrivere come si muovono le onde di taglio (onde che fanno scivolare gli strati di materiale l'uno sull'altro) in questi materiali "stravaganti" e complessi.

Ecco i concetti chiave spiegati con analogie semplici:

1. Il Problema: La Molla vs. Il Gomma da Masticare

Immagina due strati di terra: uno sopra (la crosta terrestre) e uno sotto (il mantello).

• Il vecchio modello (Lineare): Immagina che entrambi gli strati siano fatti di molle d'acciaio identiche. Se lanci un'onda, questa viaggia a una velocità fissa e prevedibile. Se le velocità sono diverse, l'onda rimane intrappolata all'interfaccia, come un'auto che corre su una corsia veloce ma non può uscire.
• Il nuovo modello (Non lineare): Immagina ora che lo strato superiore sia fatto di gomma da masticare e quello inferiore di gelatina. Quando l'onda passa, la gomma si allunga, si riscalda e cambia forma. Non è più una semplice molla. La velocità dell'onda non è più fissa: dipende da quanto è forte la spinta iniziale. Se spingi forte, l'onda diventa più veloce o più lenta in modo imprevedibile.

2. La Soluzione: Un'Equazione "Intelligente"

Gli autori hanno derivato una nuova equazione matematica (un'equazione differenziale) che tiene conto di questa complessità.

• I termini cubici e quintici: Nella loro equazione, compaiono termini matematici strani (come $x^3$ o $x^5$). In parole povere, questi termini dicono al modello: "Attenzione! Più l'onda è forte, più il materiale cambia comportamento". È come dire a un'auto: "Più premi l'acceleratore, più la macchina diventa pesante e lenta a reagire".
• La viscosità (L'attrito): Hanno anche aggiunto una parte che simula l'attrito interno (viscosità). Immagina di muovere la mano nell'acqua: più vai veloce, più senti resistenza. Questo aiuta a spiegare perché le onde sismiche perdono energia e si smorzano col tempo, invece di rimbalzare per sempre.

3. L'Esperimento Virtuale: La "Bomba" di Gomma

Per testare la loro teoria, i ricercatori hanno creato simulazioni al computer (come un videogioco fisico molto avanzato).

• Hanno creato un mondo virtuale con due strati di materiale diverso.
• Hanno simulato un'esplosione improvvisa (come un piccolo terremoto) in un punto preciso.
• Cosa hanno visto?
  • Le onde si propagano, rimbalzano e interagiscono con l'interfaccia tra i due strati.
  • Anche nel caso "non lineare" (gomma e gelatina), le onde tendono a comportarsi in modo simile alle onde classiche: rimangono intrappolate tra i due strati se le velocità sono giuste.
  • Tuttavia, c'è una differenza: nel lungo periodo, l'onda tende a viaggiare alla velocità del materiale più veloce, "catturando" l'energia e disperdendo il resto.

4. Perché è importante?

Questa ricerca è fondamentale per diversi motivi:

• Sismologia: Aiuta a capire meglio come i terremoti si propagano attraverso la crosta terrestre, che non è fatta di roccia perfetta ma di materiali complessi.
• Medicina: Potrebbe aiutare a studiare come le onde sonore o meccaniche si muovono attraverso i tessuti del corpo umano (che sono molto elastici e viscosi), utile per diagnosi mediche o terapie.
• Materiali: Aiuta a progettare materiali migliori, come gomme per pneumatici o schiume per assorbire gli urti, capendo come reagiscono a sollecitazioni forti.

In Sintesi

Gli autori hanno preso le vecchie regole della fisica delle onde (che funzionano bene per le molle leggere) e le hanno aggiornate per descrivere il mondo reale, dove i materiali sono "testardi", si deformano in modo strano e perdono energia. Hanno dimostrato che, anche quando le cose diventano molto complicate e non lineari, le onde di Love mantengono alcune delle loro regole fondamentali, ma con un comportamento più ricco e realistico.

È come passare da una mappa disegnata su un foglio di carta liscio a una mappa interattiva 3D che mostra anche le colline, le valli e il meteo in tempo reale: molto più difficile da calcolare, ma molto più vicina alla realtà.

Sommario tecnico

1. Problema e Contesto

Il lavoro affronta la modellazione matematica e la simulazione numerica delle onde di taglio (onde SH) non lineari che si propagano attraverso materiali iperelastici e iper-viscoelastici incompressibili.

• Contesto: Le onde sismiche, in particolare le onde di Love (onde di superficie polarizzate orizzontalmente), sono tradizionalmente descritte dalla teoria dell'elasticità lineare (Legge di Hooke). Tuttavia, materiali come tessuti biologici, polimeri, gomme e formazioni geologiche mostrano comportamenti non lineari significativi sotto grandi deformazioni.
• Limiti degli approcci attuali: I modelli lineari trascurano gli effetti non lineari nello stress-sforzo, che possono bilanciare gli effetti dispersivi. Inoltre, le grandi deformazioni vicino al fuoco di un terremoto o in materiali complessi richiedono modelli che vadano oltre l'approssimazione di piccole deformazioni.
• Obiettivo: Derivare equazioni generali per le onde di taglio non lineari in materiali incompressibili con una funzione di densità di energia di deformazione ($W^H$) arbitraria, specializzandole per il modello cubico di Yeoh, e includere gli effetti viscoelastici (dissipazione energetica).

2. Metodologia

Gli autori utilizzano un approccio basato sulla meccanica del continuo e sull'analisi matematica, combinato con simulazioni numeriche avanzate.

A. Formulazione Teorica

1. Quadro Iperelastico:
   • Si parte dalle equazioni del moto nel riferimento Lagrangiano per un corpo continuo incompressibile.
   • Si assume un campo di spostamento di taglio in direzione $Y$ che si propaga nel piano $(X, Z)$ (moto anti-piano).
   • Viene derivata un'equazione alle derivate parziali (PDE) generale per lo spostamento $v(X, Z, t)$ valida per qualsiasi funzione di energia $W^H$.
   • Restrizione: Per evitare condizioni di compatibilità complesse sulla pressione idrostatica quando $W^H$ dipende dal secondo invariante ($I_2$), il modello si restringe ai materiali neo-Hookeani generalizzati ($W^H = W^H(I_1)$), dove $I_1$ è il primo invariante del tensore di deformazione di Cauchy-Green.
2. Modello Specifico (Yeoh):
   • Si specializza il modello per la funzione di energia di Yeoh cubica: $W^H = \mu(I_1-3) + \frac{1}{2}M(I_1-3)^2 + \frac{1}{3}A(I_1-3)^3$.
   • Questo porta a un'equazione d'onda non lineare contenente termini polinomiali cubici e quintici nelle derivate dello spostamento.
3. Estensione Viscoelastica:
   • Viene introdotto un potenziale di dissipazione $W^V$ per modellare la viscosità.
   • Si deriva un'equazione di moto "iper-viscoelastica" che include termini di dispersione con derivate miste (es. $u_{xxt}$), permettendo di studiare la dissipazione dell'energia nel tempo.

B. Analisi Numerica

• Metodo: Vengono utilizzate simulazioni numeriche in 2D (spazio + tempo) basate sul Metodo delle Linee (MOL).
• Implementazione:
  • Il dominio spaziale è discretizzato con una griglia non uniforme (più fitta vicino alla sorgente e all'interfaccia per catturare meglio le non linearità).
  • Le derivate spaziali sono approssimate con differenze finite del secondo ordine.
  • Il sistema di ODE risultante è integrato nel tempo utilizzando ode23 di MATLAB (metodi Runge-Kutta), scelto per la sua capacità di gestire sistemi non rigidi e per l'efficienza computazionale rispetto ai solver rigidi che richiederebbero il calcolo di Jacobiani enormi.
• Configurazione: Si simula un'esplosione gaussiana improvvisa in un sistema a due strati (uno strato superiore su un semispazio inferiore) per generare onde di Love.

C. Analisi di Simmetria

• Per il caso monodimensionale (1D), viene applicata l'analisi di simmetria di Lie per derivare soluzioni esatte invarianti.
• Si identificano generatori di simmetria che portano a soluzioni di tipo onda viaggiante, sebbene queste risultino spesso non limitate (unbounded) su tutto il dominio reale.

3. Risultati Chiave

1. Equazioni Derivate

• È stata ottenuta una nuova PDE non lineare per le onde di taglio in materiali iperelastici incompressibili (Eq. 4.9 e 4.10).
• L'equazione include termini non lineari di ordine 3 e 5, accoppiati con la pressione idrostatica.
• L'estensione viscoelastica introduce termini di dissipazione che agiscono come termini di dispersione mista.

2. Comportamento delle Onde di Love Non Lineari

• Condizione di Esistenza: Le simulazioni confermano che, anche nel caso pienamente non lineare, la velocità di fase $v$ delle onde di interfaccia e di superficie soddisfa la condizione classica delle onde di Love lineari: $c_1 < |v| < c_2$, dove $c_1$ e $c_2$ sono le velocità delle onde di taglio nei due materiali.
• Evoluzione Temporale:
  • Inizialmente, le velocità delle onde sulla superficie e all'interfaccia possono differire.
  • Nel tempo, le velocità tendono ad assestarsi sul valore della velocità del materiale più veloce ($\max(c_1, c_2)$).
  • Questo comportamento è osservato sia nel caso puramente iperelastico che in quello iper-viscoelastico.

3. Effetti della Non Linearità e della Viscosità

• Iperelasticità: L'aumento del coefficiente non lineare ($M$) genera oscillazioni ad alta frequenza di piccola ampiezza, ma il comportamento qualitativo rimane simile al caso lineare.
• Iper-viscoelasticità: L'introduzione della viscosità ($\eta$) causa una rapida dispersione dell'energia. Le onde si attenuano più velocemente e il profilo d'onda diventa più prevedibile e stabile al crescere di $\eta$.
• Rottura Numerica (Appendice): Gli autori riportano la scoperta di soluzioni che mostrano "rottura" (breaking) numerica in configurazioni radialmente simmetriche, dove la derivata spaziale sviluppa discontinuità o oscillazioni spurie, un fenomeno non osservato nel caso lineare puro.

4. Soluzioni Esatte (1D)

• L'analisi di simmetria ha prodotto soluzioni esatte per il caso 1D, che mostrano un profilo d'onda stabile. Le simulazioni numeriche confermano che le soluzioni non lineari evolvono verso forme simili a queste soluzioni invarianti, con un'ampiezza che decade a causa della viscosità.

4. Significato e Contributi

• Generalizzazione Teorica: Il lavoro estende la teoria classica delle onde di Love (lineare) al regime di grandi deformazioni e materiali complessi, fornendo un quadro matematico rigoroso per materiali iper-viscoelastici.
• Validazione Numerica: Dimostra che la condizione di esistenza delle onde di Love ($c_1 < v < c_2$) è robusta e sopravvive anche in presenza di forti non linearità e dissipazione.
• Applicabilità: Il modello è rilevante per:
  • Sismologia: Migliore comprensione della propagazione delle onde superficiali vicino al fuoco di un terremoto (dove le deformazioni sono grandi).
  • Scienza dei Materiali e Biomedicina: Modellazione delle onde di taglio in tessuti biologici, gel e materiali polimerici, dove la viscosità e la non linearità sono dominanti.
• Metodologia Computazionale: L'uso efficace del metodo delle linee con griglie non uniformi e solver adattivi offre un approccio efficiente per simulare problemi di onde non lineari in domini stratificati complessi.

In conclusione, il paper fornisce un ponte fondamentale tra la teoria elastica lineare classica e la realtà fisica dei materiali non lineari, confermando che le caratteristiche fondamentali delle onde di Love persistono anche in regimi estremi, pur con dinamiche di ampiezza e velocità modificate dalla non linearità e dalla dissipazione.