A Physically Constrained Energy-Based Excess Pore Pressure Model for Saturated Sands with an Ordinary Differential Equation Interpretation

Sintesi Tecnica: Un Modello di Eccesso di Pressione Interstiziale Basato sull'Energia e Vincolato Fisicamente per Sabbie Sature

Problematica
La liquefazione del suolo rimane un rischio geotecnico critico, tradizionalmente valutata attraverso metodi basati sullo stress, sulla deformazione, sulla pressione interstiziale e sull'energia. Tra questi, i Metodi Basati sull'Energia (EBM) offrono un vantaggio distintivo, poiché collegano la generazione di eccesso di pressione interstiziale all'energia dissipata cumulativa, adattandosi a carichi sismici irregolari senza l'incertezza delle conversioni in cicli equivalenti. La formulazione EBM più influente è la legge di potenza di Berrill–Davis (1985) ($r_u = \alpha W^\beta$). Tuttavia, questo studio identifica due limitazioni fondamentali nel modello di Berrill–Davis:

1. Risposta Illimitata: La legge di potenza è matematicamente illimitata, permettendo al rapporto di eccesso di pressione interstiziale ($r_u$) di superare l'unità, il che è fisicamente impossibile per le sabbie sature dove $r_u$ dovrebbe tendere asintoticamente a 1,0 al verificarsi della liquefazione.
2. Interdipendenza dei Parametri: I parametri di calibrazione $\alpha$ e $\beta$ sono spesso trattati come indipendenti, sebbene studi empirici dimostrino che siano mutualmente correlati. Questa dipendenza complica la calibrazione, riduce l'interpretabilità fisica e ostacola l'estrapolazione.

Metodologia
Gli autori propongono un modello riformulato per sabbie sature, uniformi e pulite, che affronta queste limitazioni attraverso un processo di derivazione in più fasi:

1. Analisi e Decomposizione dei Dati: Utilizzando dati di triassiali ciclici di sabbia Ottawa F65 (LEAP-2017) e sabbia Akpınar, viene analizzata la relazione tra il rapporto di pressione interstiziale normalizzato ($r_u$) e l'energia dissipata normalizzata ($W_N$) su scala log–log. Gli autori osservano che i dati non seguono una singola linea retta, ma consistono in almeno tre segmenti lineari.
2. Approssimazione Polinomiale: I dati sono approssimati da un polinomio di sesto ordine senza termine costante. Questa funzione viene decomposta in un termine di crescita primario ($W_N^b$) e una funzione di mediazione $g(W_N)$.
3. Espansione di Taylor e Semplificazione: La funzione di mediazione viene ulteriormente analizzata. Un'espansione di Taylor della componente esponenziale rivela che un polinomio di terzo ordine rappresenta adeguatamente l'andamento dei dati. Ciò porta a una forma di modello candidata: $r_u = W_N^b (a_0 + a_1 W_N + a_2 W_N^2 + a_3 W_N^3)$.
4. Interpretazione dell'Equazione Differenziale Ordinaria (ODE): Riconoscendo la struttura matematica, il modello viene reinterpretato come la soluzione di un'equazione differenziale ordinaria di Cauchy–Euler non omogenea. Il lato sinistro rappresenta la risposta interna del sistema, mentre il lato destro rappresenta le condizioni di carico esterno.
5. Vincoli Fisici: Per eliminare fluttuazioni non fisiche e l'interdipendenza dei parametri, il modello è sottoposto a specifici vincoli al contorno nel punto di liquefazione ($W_N = 1, r_u = 1$):
   • $r_u = 1$ quando $W_N = 1$.
   • Una tangente orizzontale ($r_u' = 0$) a $W_N = 1$, garantendo che la curva si appiattisca al raggiungimento del limite.
   • Crescita monotona ($r_u' > 0$) durante tutto il processo di carico.
   • Una specifica relazione tra i coefficienti ($a_2 = -3a_3$) derivata dalla scomparsa del termine non omogeneo alla liquefazione.

Questi vincoli determinano analiticamente i coefficienti del polinomio di mediazione in funzione del singolo parametro $b$, riducendo il modello a una forma governata esclusivamente da $b$.

Contributi Chiave

• Formulazione Vincolata Fisicamente: Il modello proposto è costruito per soddisfare intrinsecamente $0 \le r_u \le 1$, prevenendo la sovrastima non fisica della pressione interstiziale comune nel modello di Berrill–Davis.
• Singolo Parametro di Calibrazione: Collegando analiticamente i coefficienti di mediazione all'esponente di crescita $b$, il modello elimina la necessità di calibrare indipendentemente molteplici parametri, risolvendo il problema dell'interdipendenza riscontrato nelle formulazioni precedenti.
• Framework ODE: Lo studio fornisce una nuova interpretazione della relazione empirica della pressione interstiziale come soluzione di un'equazione differenziale ordinaria di Cauchy–Euler non omogenea. Questo framework separa la risposta in una parte omogenea (che governa le caratteristiche interne del materiale, dominante nelle fasi iniziali) e una parte particolare (che governa le condizioni di carico esterno).
• Validazione: Il modello è validato rispetto ad alta qualità di dati triassiali ciclici per sabbie Ottawa F65 e Akpınar attraverso una gamma di densità relative (da sciolte a molto dense), raggiungendo un $R^2 \ge 0,99$.

Risultati

• Il modello semplificato (Eq. 20) riproduce con successo le storie di pressione interstiziale misurate per le sabbie Ottawa F65 e Akpınar.
• È stato riscontrato che il parametro $b$ varia con la densità relativa (ad es., $b \approx 0,36$ per una densità relativa del 51%, $b \approx 0,43$ per l'80–91% di densità relativa).
• Il modello cattura l'evoluzione non lineare della pressione interstiziale, inclusa la rapida salita iniziale e il successivo plateau man mano che ci si avvicina alla liquefazione.
• Il confronto con l'originale modello di Berrill–Davis dimostra che la formulazione proposta offre una descrizione più coerente dell'evoluzione della pressione interstiziale, in particolare vicino alla condizione terminale di liquefazione completa, dove il modello originale presenta una dispersione significativa e una crescita illimitata.
• L'applicazione del modello semplificato al restante dataset LEAP-2017 conferma la sua capacità di simulare la generazione di pressione interstiziale attraverso varie densità relative e rapporti di sforzo ciclico, a condizione che $b \ge 0,27$ per evitare massimi locali.

Significatività
L'articolo sostiene che questo lavoro offre un'alternativa valida, compatta e fisicamente coerente alla classica legge di potenza di Berrill–Davis per studi di liquefazione che coinvolgono sabbie uniformi e pulite. Limitando la risposta e rimuovendo l'interdipendenza dei parametri attraverso una derivazione basata su ODE, il modello semplifica la calibrazione mantenendo un'alta accuratezza. Gli autori sottolineano che, sebbene la formulazione sia empirica, la sua derivazione da una regressione sotto rigorosi vincoli fisici le permette di catturare il meccanismo sottostante della liquefazione: una fase iniziale governata dalle caratteristiche granulari interne (soluzione omogenea) seguita da una fase in cui il materiale perde la sua memoria meccanica e richiede un'energia sostanziale per raggiungere la piena liquefazione (soluzione particolare). L'ambito è esplicitamente limitato a sabbie sature, uniformi e pulite.