🔢 mathematics

SoliDualSPHysics: An extension of DualSPHysics for solid mechanics with hyperelasticity, plasticity, and fracture

Questo articolo introduce SoliDualSPHysics, una nuova estensione open-source di DualSPHysics accelerata su GPU che utilizza una formulazione lagrangiana totale e un approccio a campo di fase per consentire simulazioni unificate e su larga scala di iperelasticità, plasticità a deformazione finita e frattura fragile in solidi deformabili senza il tracciamento esplicito delle crepe.

Autori originali: Mohammad Naqib Rahimi, George Moutsanidis

Pubblicato 2026-06-23

📖 5 min di lettura🧠 Approfondimento

CC BY 4.0

Autori originali: Mohammad Naqib Rahimi, George Moutsanidis

Articolo originale sotto licenza CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Questa è una spiegazione generata dall'IA dell'articolo qui sotto. Non è stata scritta né approvata dagli autori. Per precisione tecnica, consulta l'articolo originale. Leggi il disclaimer completo

Immagina di avere un sacco invisente di biglie gigante. Nel mondo delle simulazioni informatiche, questo sacco di biglia viene spesso usato per modellare come l'acqua scorre o come una fuoriuscita di petrolio si diffonde. Questo perché le biglie possono muoversi liberamente, scontrarsi tra loro e comprimersi senza incastrarsi come un groviglio di fili di lana (che è ciò che accade con i modelli informatici tradizionali basati su griglia).

Questo articolo presenta uno strumento potente e innovativo chiamato SoliDualSPHysics. Immagina che sia un grande aggiornamento software che prende questo "sacco di biglie" e insegna loro come comportarsi come oggetti solidi (come metallo, gomma o osso) invece che solo come liquidi.

Ecco una suddivisioni di ciò che fa questo software, utilizzando analogie semplici:

1. L'idea centrale: dall'acqua al solido

Precedentemente, il software (DualSPHysics) era ottimo nel simulare i fluidi. Gli autori volevano renderlo capace di simulare solidi che possono allungarsi, schiacciarsi e rompersi.

L'analogia: Immagina di avere una simulazione di una piscina. Ora, vuoi simulare un elastico che si spezza o una trave di metallo che si piega. Questo nuovo software permette alle "biglie" di mantenere la loro forma come un oggetto solido pur continuando a muoversi liberamente come un fluido.

2. Come gestisce l'allungamento e la flessione (Iperelasticità e Plasticità)

I solidi si comportano in modi diversi quando vengono tirati o spinti.

Iperelasticità (L'elastico): Alcuni materiali, come la gomma, si allungano e poi tornano perfettamente alla loro forma originale. Il software può simulare questo comportamento "elastico".
Plasticità (Il pongo): Altri materiali, come il metallo tenero o l'argilla, si allungano e rimangono allungati. Non tornano indietro. Il software può ora simulare questa "deformazione permanente" dove il materiale cambia forma per sempre.

3. Come gestisce la rottura (Frattura)

Questa è la parte più entusiasmante. In passato, se volevi simulare una crepa in un pezzo di vetro, dovevi dire al computer esattamente dove sarebbe andata la crepa, come disegnare una linea su una mappa prima che la gara inizi. Se la crepa fosse andata altrove, la simulazione sarebbe fallita.

Il nuovo approccio (Il metodo della "Nebbia"): Questo software utilizza un approccio chiamato Phase-Field. Immagina che il materiale non sia solido ovunque; immagina invece un sottile strato di "nebbia" che si diffonde attraverso il materiale.
- Dove la nebbia è densa, il materiale è sano.
- Dove la nebbia è sottile, il materiale è danneggiato.
- Dove la nebbia scompare completamente, il materiale si è rotto.
Perché è fantastico: Non devi dire al computer dove andrà la crepa. Il computer la individua naturalmente. La crepa può iniziare, crescere, dividersi in due rami o fondersi con un'altra crepa, proprio come accade nella realtà, senza che il computer si confonda o abbia bisogno di istruzioni speciali.

4. Il motore: Velocità e Potenza

Simulare milioni di biglie è incredibilmente difficile per un computer.

L'analogia: Immagina di dover calcolare la traiettoria di un milione di palle da biliardo. Farlo una per una richiederebbe un tempo infinito.
La soluzione: Questo software è costruito per girare su GPU (le potenti schede grafiche presenti nei computer da gioco). Utilizza migliaia di piccoli processori che lavorano insieme contemporaneamente (elaborazione parallela). Questo rende la simulazione incredibilmente veloce, permettendo agli scienziati di modellare enormi strutture con milioni di particelle in un tempo ragionevole.

5. Flessibilità: Il "Telecomando"

Il software viene fornito con un "telecomando" speciale (scritto in un linguaggio chiamato XML).

Cosa fa: Puoi dire alla simulazione di spingere una parte specifica dell'oggetto, tirarla o scuoterla, usando formule matematiche. Puoi persino dire di applicare una forza solo sul lato sinisto di una trave o solo durante il primo secondo della simulazione.
Il vantaggio: Questo rende molto facile per gli ingegneri impostare test complessi senza dover riscrivere l'intero programma ogni volta.

Cosa non può ancora fare (Le note tecniche)

Gli autori sono molto onesti riguardo ai limiti attuali del loro strumento:

Funziona meglio con materiali che sono uguali in tutte le direzioni (isotropi), come il metallo puro o la gomma. Non gestisce ancora materiali complessi come la fibra di carbonio (che è più resistente in una direzione rispetto a un'altra).
Simula molto bene la rottura fragile (come il vetro o la ceramica). Non simula ancora la rottura duttile (come un filo di metallo che si allunga e si assottiglia fino a spezzarsi), perché questo comporta un tipo diverso di fisica.
Attualmente tratta la plasticità (piegatura permanente) e la frattura (rottura) come eventi separati. Non simula un materiale che si piega e si rompe contemporaneamente in questa specifica versione.

Riassunto

SoliDualSPHysics è uno strumento gratuito e open-source che trasforma un programma di simulazione di fluidi in un potente simulatore di meccanica dei solidi. Utilizza un metodo a "nebbia" per permettere alle crepe di apparire naturalmente senza bisogno di una mappa, gira velocemente sulle moderne schede grafiche e permette agli ingegneri di testare come i materiali si allungano, si piegano e si rompono sotto condizioni estreme. È progettato per essere un laboratorio virtuale flessibile e ad alta velocità per esperimenti reali.

Sintesi Tecnica: SoliDualSPHysics

Definizione del Problema

La tecnica Smoothed Particle Hydrodynamics (SPH) si è affermata come un potente metodo meshfree per la dinamica dei fluidi e la meccanica dei solidi a grandi deformazioni grazie alla sua natura lagrangiana. Tuttavia, l'attuale panorama dei framework SPH open-source ad alte prestazioni è fortemente sbilanciato verso la fluidodinamica computazionale (CFD) e l'interazione fluido-struttura (FSI) che coinvolge corpi rigidi o elasticità a piccole deformazioni. Esiste una lacuna significativa in termini di framework unificati capaci di gestire la meccanica dei solidi deformabili autonomi con comportamenti costitutivi avanzati, specificamente la plasticità a deformazione finita e la frattura fragile fisicamente consistente.

Le formulazioni di frattura SPH esistenti spesso si affidano a modelli di danno empirici ad hoc, approcci di cracking delle particelle o modelli di zona coesiva, che soffrono di dipendenza dal mesh, risultati non convergenti o la necessità di tracking esplicito delle crepe o arricchimento locale. Inoltre, l'integrazione di queste capacità avanzate di meccanica dei solidi in architetture scalabili CPU/GPU adatte a simulazioni su larga scala rimane una sfida. La specifica necessità identificata è un framework SPH unificato che supporti l'iperelasticità, la plasticità a deformazione finita e la frattura fragile senza il tracking esplicito delle crepe, mantenendo al contempo l'efficienza computazionale attraverso la parallelizzazione.

Metodologia

Gli autori introducono SoliDualSPHysics, un'estensione open-source del framework DualSPHysics accelerata tramite GPU. La metodologia si basa su una formulazione Total Lagrangian SPH (TLSPH), che consente l'applicazione diretta di carichi esterni e condizioni al contorno nella configurazione indeformata.

Framework Teorico

Bilancio del Momento: Il solver impiega una discretizzazione TLSPH in forma forte dell'equazione del bilancio del momento. Utilizza una correzione del gradiente del kernel del primo ordine per ripristinare la consistenza e include un termine di viscosità artificiale per mitigare le instabilità numeriche e i modi di energia zero spurii.
Modellazione Costitutiva:
- Iperelasticità: Il framework supporta i modelli di St. Venant–Kirchhoff e neo-Hookeano comprimibile. Questi sono accoppiati a una formulazione phase-field per simulare la frattura fragile.
- Plasticità a Deformazione Finita: È implementata un'estensione alla plasticità $J_2$ indipendente dalla velocità con incrudimento isotropo. Questa utilizza una decomposizione moltiplicativa del gradiente di deformazione ( $F = F^e F^p$ ) e un algoritmo di ritorno tramite mappa esponenziale per garantire l'oggettività di frame e la consistenza termodinamica sotto grandi rotazioni.
- Nota sull'Accoppiamento: Nella corrente implementazione, plasticità e frattura phase-field non sono accoppiate; il modello di plasticità descrive la risposta elastoplastica senza frattura, mentre il modello phase-field gestisce la frattura fragile in materiali iperelastici.
Frattura Fragile (Phase-Field): La frattura è modellata utilizzando una formulazione phase-field iperbolica. A differenza delle formulazioni ellittiche che richiedono la risoluzione di sistemi lineari globali (incompatibili con l'integrazione temporale esplicita), questo approccio utilizza un'equazione di evoluzione iperbolica per la variabile di danno $s$ . Ciò consente l'iniziazione, la propagazione, la ramificazione e la coalescenza delle crepe di emergere naturalmente senza il tracking esplicito o l'arricchimento. Il modello impone l'irreversibilità tramite un funzionale di storia e limita la propagazione delle crepe alle regioni di trazione scomponendo l'energia di deformazione elastica.
Integrazione Temporale: Il sistema è avanzato utilizzando schemi di integrazione temporale espliciti (Verlet o Symplectic/Leapfrog), con passi temporali vincolati da condizioni di tipo CFL basate sulla propagazione delle onde e l'accelerazione.

Architettura Software e Implementazione

Piattaforma: Il codice è scritto in C++ e CUDA, sfruttando l'architettura parallela nativa CPU (OpenMP) e GPU (CUDA) di DualSPHysics.
Input/Output: La configurazione è gestita tramite file di input XML. L'estensione introduce nuovi blocchi (<deformstrucs>, <mathexpressions>) per definire corpi deformabili, proprietà dei materiali, condizioni al contorno ed espressioni definite dall'utente.
Espressioni Definite dall'Utente: Una caratteristica chiave è l'inclusione di un parser di espressioni matematiche che permette agli utenti di prescrivere quantità dipendenti dal tempo e dallo spazio (ad esempio, campi di velocità, carichi di trazione) utilizzando variabili integrate (posizione, tempo, spostamento) e operatori matematici/logici standard.
Condizioni al Contorno: Condizioni al contorno di velocità e trazione flessibili possono essere applicate globalmente o a specifici gruppi superficiali definiti da geometrie ausiliarie.

Contributi Chiave

Framework Unificato: SoliDualSPHysics fornisce un ambiente unico, open-source e accelerato da GPU per simulare l'iperelasticità, la plasticità a deformazione finita e la frattura fragile all'interno del paradigma SPH.
Frattura Phase-Field in SPH: Implementa una formulazione phase-field iperbolica per la frattura fragile che evita il costo computazionale dei solver ellittici, consentendo una dinamica esplicita efficiente per la ramificazione e la coalescenza delle crepe senza criteri ad hoc.
Plasticità a Deformazione Finita: L'integrazione di un modello di plasticità $J_2$ moltiplicativo con ritorno tramite mappa esponenziale estende l'applicabilità del codice a problemi plastici a grandi deformazioni.
Scalabilità: Utilizzando l'infrastruttura CPU/GPU di DualSPHysics, il software ottiene una sostanziale accelerazione per simulazioni su larga scala che coinvolgono milioni di particelle.
Flessibilità: L'inclusione di un parser di espressioni definite dall'utente e la definizione flessibile delle condizioni al contorno ne potenziano l'adattabilità per scenari di carico complessi e applicazioni di ricerca.
Disponibilità Open Source: Il codice sorgente completo, la documentazione e i casi di benchmark sono resi liberamente disponibili tramite un repository pubblico su GitHub per garantire la riproducibilità.

Risultati e Validazione

Il software è stato verificato e validato rispetto a una serie di problemi numerici di benchmark e dati sperimentali:

Elasto-Dinamica: La libera oscillazione di travi e piastre cantilever (2D e 3D) ha dimostrato un eccellente accordo con le soluzioni analitiche e i risultati di riferimento del Metodo degli Elementi Finiti (FEM). Le simulazioni hanno mostrato uno scambio di energia quasi conservativo e una dissipazione numerica minima.
Grandi Deformazioni: Una trave cantilever 3D soggetta a una grande trazione superficiale distribuita ha confermato la capacità di gestire la non linearità geometrica e le condizioni al contorno complesse senza oscillazioni spurie.
Dinamica Torsionale: Un caso di colonna torsionale 3D ha validato la capacità del framework di risolvere la propagazione delle onde torsionali e il recupero elastico, con risultati corrispondenti agli studi di riferimento.
Frattura: Un benchmark di ramificazione dinamica delle crepe (piastra pre-intagliata sotto trazione) ha dimostrato la capacità del modello di catturare pattern di frattura complessi, inclusa la ramificazione, in modo coerente con la letteratura phase-field stabilita.
Performance: L'implementazione GPU ha dimostrato di fornire un'accelerazione significativa rispetto alle esecuzioni solo CPU, supportando simulazioni scalabili.

Significato e Ambito

L'articolo posiziona SoliDualSPHysics come un passo significativo nel colmare il divario tra i solver fluidi SPH ad alte prestazioni e le avanzate capacità della meccanica dei solidi. Il suo significato primario risiede nel fornire uno strumento unificato, open-source e accelerato da GPU che elimina la necessità di tracking ad hoc delle crepe nelle simulazioni di frattura fragile e supporta la plasticità a deformazione finita all'interno di un framework meshfree.

Gli autori notano esplicitamente le limitazioni dell'attuale rilascio:

È limitato alla risposta dei materiali isotropi (iperelasticità isotropa, plasticità $J_2$ isotropa).
Plasticità e frattura non sono accoppiate; i meccanismi di frattura duttile non sono inclusi.
Modelli costitutivi anisotropi, materiali compositi e operatori dedicati al controllo dell'hourglass non sono ancora implementati.

Nonostante queste limitazioni, il lavoro stabilisce una base robusta per futuri sviluppi nella meccanica dei solidi e nella modellazione della frattura all'interno della comunità SPH, offrendo una piattaforma riproducibile ed estendibile per i ricercatori.