🔢 mathematics

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

Dit artikel introduceert SoliDualSPHysics, een nieuwe, open-source, door GPU versnelde uitbreiding van DualSPHysics die een totale Lagrangiaanse formulering en een faseveldbenadering gebruikt om verenigde, grootschalige simulaties van hyperelasticiteit, eindvervormingsplasticiteit en brosse breuk in vervormbare vaste stoffen mogelijk te maken zonder expliciete scheurtracking.

Oorspronkelijke auteurs: Mohammad Naqib Rahimi, George Moutsanidis

Gepubliceerd 2026-06-23

📖 5 min leestijd🧠 Diepgaand

CC BY 4.0

Oorspronkelijke auteurs: Mohammad Naqib Rahimi, George Moutsanidis

Oorspronkelijk artikel gelicentieerd onder CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Dit is een AI-gegenereerde uitleg van het onderstaande artikel. Het is niet geschreven of goedgekeurd door de auteurs. Raadpleeg het oorspronkelijke artikel voor technische nauwkeurigheid. Lees de volledige disclaimer

Stel je voor dat je een gigantische, onzichtbare zak knikkers hebt. In de wereld van computersimulaties wordt deze zak met knikkers vaak gebruikt om te modelleren hoe water stroomt of hoe olie uitstort. Dit komt omdat de knikkers vrij kunnen bewegen, tegen elkaar aan kunnen botsen en samengedrukt kunnen worden zonder in de knoop te raken zoals een rommelige bol wol (wat wel gebeurt bij traditionele, op rasters gebaseerde computermodellen).

Dit artikel introduceert een nieuwe, krachtige tool genaamd SoliDualSPHysics. Zie dit als een grote software-upgrade die die "zak met knikkers" leert hoe het zich moet gedragen als vaste objecten (zoals metaal, rubber of bot) in plaats van alleen maar als een vloeistof.

Hier is een overzicht van wat deze software doet, met behulp van eenvoudige analogieën:

1. De Kern: Van Water naar Vast

Voorheen was de software (DualSPHysics) erg goed in het simuleren van vloeistoffen. De auteurs wilden het mogelijk maken om vaste stoffen te simuleren die kunnen rekken, vervormen en breken.

De Analogie: Stel je voor dat je een simulatie hebt van een zwembad. Nu wil je een rubberen band simuleren die knapt of een metalen balk die buigt. Deze nieuwe software stelt de "knikkers" in staat om hun vorm vast te houden als een solide object, terwijl ze nog steeds vrij kunnen bewegen als een vloeistof.

2. Hoe het omgaat met Rekken en Buigen (Hyperelasticiteit & Plasticiteit)

Vaste stoffen gedragen zich op verschillende manieren wanneer je aan trekt of erop drukt.

Hyperelasticiteit (De Rubberen Band): Sommige materialen, zoals rubber, rekken uit en springen dan perfect terug naar hun oorspronkelijke vorm. De software kan dit "verende" gedrag simuleren.
Plasticiteit (De Play-Doh): Andere materialen, zoals zacht metaal of klei, rekken uit en blijven uitgerekt. Ze veren niet terug. De software kan nu deze "permanente vervorming" simuleren waarbij het materiaal voor altijd van vorm verandert.

3. Hoe het omgaat met Breken (Fractuur)

Dit is het meest opwindende deel. In het verleden, als je een barst in een stuk glas wilde simuleren, moest je de computer precies vertellen waar de barst zou komen, zoals het tekenen van een lijn op een landkaart voordat de race begint. Als de barst ergens anders liep, faalde de simulatie.

De Nieuwe Aanpak (De "Mist"-methode): Deze software gebruikt een zogenaamde Phase-Field aanpak. Stel je voor dat het materiaal niet overal solide is; stel je in plaats daarvan een dun laagje "mist" voor dat door het materiaal verspreidt:
- Waar de mist dik is, is het materiaal gezond.
- Waar de mist dun is, is het materiaal beschadigd.
- Waar de mist volledig verdwijnt, is het materiaal gebroken.
Waarom het cool is: Je hoeft de computer niet te vertellen waar de barst zal komen. De computer ontdekt het op natuurlijke wijze. De barst kan ontstaan, groeien, in tweeën splitsen of met een andere barst versmelten, precies zoals in het echte leven gebeurt, zonder dat de computer in de war raakt of speciale instructies nodig heeft.

4. De Motor: Snelheid en Kracht

Het simuleren van miljoenen knikkers is extreem zwaar voor een computer.

De Analogie: Stel je voor dat je de baan van een miljoen biljartballen probeert te berekenen. Dat één voor één doen zou eeuwig duren.
De Oplossing: Deze software is gebouwd om te draaien op GPU's (de krachtige grafische kaarten die in gamingcomputers zitten). Het maakt gebruik van duizenden kleine processoren die tegelijkertijd samenwerken (parallelle verwerking). Dit maakt de simulatie ongelooflijk snel, waardoor wetenschappers enorme structuren met miljoenen deeltjes in een redelijke tijd kunnen modelleren.

5. Flexibiliteit: De "Afstandsbediening"

De software wordt geleverd met een speciale "afstandsbediening" (geschreven in een taal die XML wordt genoemd).

Wat het doet: Je kunt de simulatie vertellen om een specifiek deel van het object te duwen, te trekken of te schudden, met behulp van wiskundige formules. Je kunt zelfs de opdracht geven om kracht alleen op de linkerkant van een balk uit te oefenen, of alleen tijdens de eerste seconde van de simulatie.
Het Voordeel: Dit maakt het voor ingenieurs heel eenvoudig om complexe tests op te zetten zonder telkens het hele programma te hoeven herschrijven.

Wat het (nog) niet kan (De Kleine Lettertjes)

De auteurs zijn zeer eerlijk over de huidige beperkingen van hun tool:

Het werkt het beste met materialen die in alle richtingen hetzelfde zijn (isotroop), zoals puur metaal of rubber. Het kan complexe materialen zoals koolstofvezel (die sterker zijn in één specifieke richting) nog niet aan.
Het simuleert bros breken (zoals glas of keramiek) erg goed. Het simuleert nog niet het ductiele breken (zoals een metalen draad die uitrekt en dunner wordt tot hij knapt), omdat dit een andere soort fysica vereist.
Het behandelt momenteel plasticiteit (permanente vervorming) en fractuur (breken) als aparte gebeurtenissen. Het simuleert in deze specifieijke versie niet een materiaal dat tegelijkertijd vervormt én breekt.

Samenvatting

SoliDualSPHysics is een gratis, open-source tool die een vloeistofsimulatieprogramma transformeert in een krachtige simulator voor vaste mechanica. Het gebruikt een "mist"-methode waardoor barsten natuurlijk kunnen verschijnen zonder dat er een kaart nodig is, het draait razendsnel op moderne grafische kaarten, en stelt ingenieurs in staat om te testen hoe materialen rekken, buigen en breken onder extreme omstandigheden. Het is ontworpen als een flexibel, hoogwaardig laboratorium voor virtuele experimenten.

Technische Samenvatting: SoliDualSPHysics

Probleemstelling

Smoothed Particle Hydrodynamics (SPH) heeft zich gevestigd als een krachtige meshvrije methode voor vloeistofdynamica en de mechanica van vaste stoffen met grote vervormingen vanwege het Lagrangiaanse karakter. De huidige landschappen van high-performance, open-source SPH-frameworks zijn echter sterk scheefgetrokken naar computationele vloeistofmechanica (CFD) en vloeistof-structuurinteractie (FSI) waarbij rigide lichamen of kleine elastische vervormingen betrokken zijn. Er is een aanzienlijke kloof in verenigde frameworks die in staat zijn om autonome vervormbare vaste stofmechanica met geavanceerde constitutieve gedragingen te verwerken, specifiek eindige-vervormingsplasticiteit en fysisch consistente brosse breuk.

Bestaande SPH-breukformuleringen vertrouwen vaak op ad hoc empirische schade-modellen, deeltjes-cracking benaderingen of cohésieve zone-modellen, die lijden onder mesh-afhankelijkheid, niet-convergerende resultaten of de noodzaak voor expliciete scheurtrakking en lokale verrijking. Bovendien blijft het integreren van deze geavanceerde vaste stofmechanica-mogelijkheden in schaalbare CPU/GPU-architecturen die geschikt zijn voor grootschalige simulaties een uitdaging. De geïdentificeerde specifieke behoefte is een verenigde SPH-formulering die hyperelasticiteit, eindige-vervormingsplasticiteit en brosse breuk ondersteunt zonder expliciete scheurtrakking, terwijl de computationele efficiëntie behouden blijft door parallellisatie.

Methodologie

De auteurs introduceren SoliDualSPHysics, een open-source, GPU-versnelde uitbreiding van het DualSPHysics-framework. De methodologie is gebaseerd op een Total Lagrangian SPH (TLSPH) formulering, die directe toepassing van externe belastingen en randvoorwaarden in de onvervormde configuratie mogelijk maakt.

Theoretisch Kader

Impulsbalans: De solver maakt gebruik van een strong-form TLSPH-discretisatie van de impulsbalansvergelijking. Het maakt gebruik van een eerste-orde kernel-gradiëntcorrectie om consistentie te herstellen en bevat een kunstmatige viscositeitsterm om numerieke instabiliteiten en spurieuze nul-energie-modi te mitigeren.
Constitutieve Modellering:
- Hyperelasticiteit: Het framework ondersteunt de St. Venant–Kirchhoff en de samendrukbare neo–Hooke modellen. Deze zijn gekoppeld aan een phase-field formulering om brosse breuk te simuleren.
- Eindige-vervormingsplasticiteit: Een uitbreiding naar snelheidsonafhankelijke $J_2$ -elastoplasticiteit met isotrope verharding is geïmplementeerd. Dit maakt gebruik van een multiplicatieve decompositie van de vervormingsgradiënt ( $F = F^e F^p$ ) en een exponential-map return algoritme om frame-objectiviteit en thermodynamische consistentie onder grote rotaties te waarborgen.
- Koppelingsnotitie: In de huidige implementatie zijn plasticiteit en phase-field breuk niet gekoppeld; het plasticiteitsmodel beschrijft de elastoplastische respons zonder breuk, terwijl het phase-field model brosse breuk behandelt in hyperelastische materialen.
Brosse Breuk (Phase-Field): Breuk wordt gemodelleerd met behulp van een hyperbolische phase-field formulering. In tegenstelling tot elliptische formuleringen die vereisen dat globale lineaire stelsels worden opgelost (incompatibel met expliciete tijdstap-integratie), gebruikt deze benadering een hyperbolische evolutievergelijking voor de schadevariabele $s$ . Dit maakt het mogelijk dat crack initiatie, propagatie, vertakking en coalescentie natuurlijk ontstaan zonder expliciete tracking of verrijking. Het model dwingt onomkeerbaarheid af via een historie-functional en beperkt scheurpropagatie tot trekgebieden door de elastische energie te deconstrueren.
Tijdintegratie: Het systeem wordt voortgestuven met behulp van expliciete tijdintegratieschema's (Verlet of Symplectic/Leapfrog), waarbij tijdstappen worden beperkt door CFL-type condities gebaseerd op golfvoortplanting en versnelling.

Software Architectuur en Implementatie

Platform: De code is geschreven in C++ en CUDA, waarbij gebruik wordt gemaakt van de native CPU (OpenMP) en GPU (CUDA) parallelle architectuur van DualSPHysics.
Input/Output: Configuratie wordt afgehandeld via XML-inputbestanden. De uitbreiding introduceert nieuwe blokken (<deformstrucs>, <mathexpressions>) om vervormbare lichamen, materiaaleigenschappen, randvoorwaarden en gebruikersgedefinieerde expressies te definiëren.
Gebruikersgedefinieerde Expressies: Een kernkenmerk is de inclusie van een wiskundige expressie-parser die gebruikers in staat stelt om tijd- en ruimteafhankelijke grootheden (bijv. snelheidvelden, tractiebelastingen) voor te schrijven met behulp van ingebouwde variabelen (positie, tijd, verplaatsing) en standaard wiskundige/logische operatoren.
Randvoorwaarden: Flexibele snelheid- en tractie-randvoorwaarden kunnen globaal of op specifieke oppervlaktegroepen worden toegepast die zijn gedefinieerd door hulpgeometrieën.

Belangrijkste Bijdragen

Verenigd Framework: SoliDualSPHysics biedt een enkele, open-source omgeving voor het simuleren van hyperelasticiteit, eindige-vervormingsplasticiteit en brosse breuk binnen het SPH-paradigma.
Phase-Field Breuk in SPH: Het implementeert een hyperbolische phase-field formulering voor brosse breuk die de computationele kosten van elliptische solvers vermijdt, wat efficiënte expliciete dynamica voor scheurtakking en coalescentie mogelijk maakt zonder ad hoc criteria.
Eindige-vervormingsplasticiteit: De integratie van een multiplicatief $J_2$ plasticiteitsmodel met een exponential-map return mapping breidt de toepasbaarheid van de code uit naar grote vervormingsplastische problemen.
Schaalbaarheid: Door gebruik te maken van de DualSPHysics CPU/GPU infrastructuur bereikt de software een substantiële versnelling voor grootschalige simulaties met miljoenen deeltjes.
Flexibiliteit: De inclusie van een gebruikersgedefinieerde expressie-parser en flexibele definities van randvoorwaarden verhoogt de aanpasbaarheid van de code voor complexe belastingsscenario's en onderzoeksapplicaties.
Open Source Beschikbaarheid: De volledige broncode, documentatie en benchmarkgevallen worden gratis beschikbaar gesteld via een publieke GitHub-repository om reproduceerbaarheid te garanderen.

Resultaten en Validatie

De software is geverifieerd en gevalideerd tegen een reeks benchmark numerieke problemen en experimentele gegevens:

Elasto-Dynamica: Vrije oscillatie van cantilever balken en platen (2D en 3D) toonde uitstekende overeenstemming met analytische oplossingen en referentie FEM-resultaten. De simulaties vertoonden een bijna-conserveuze energie-uitwisseling en minimale numerieke dissipatie.
Grote Vervorming: Een 3D cantilever balk onderworpen aan een grote gedistribueerde oppervlaktetractie bevestigde het vermogen om geometrische nonlineariteit en complexe randvoorwaarden te verwerken zonder spurieuze oscillaties.
Torsiedynamica: Een draaiende 3D kolom casus valideerde het vermogen van het framework om torsie-golfvoortplanting en elastische herstel te resolveren, waarbij de resultaten overeenkwamen met referentiestudies.
Breuk: Een dynamische crack branching benchmark (pre-notched plaat onder spanning) demonstreerde het vermogen van het model om complexe breukpatronen, inclusioneel van vertakking, te vangen, consistent met de gevestigde phase-field literatuur.
Prestaties: De GPU-implementatie leverde een significante versnelling op vergeleken met CPU-only runs, wat schaalbare simulaties ondersteunt.

Betekenis en Omvang

Het artikel positioneert SoliDualSPHysics als een belangrijke stap voorwaarts in het overbruggen van de kloof tussen high-performance SPH vloeistofsolvers en geavanceerde vaste stofmechanica-mogelijkheden. De primaire betekenis ligt in het bieden van een verenigde, open-source en GPU-versnelde tool die de noodzaak voor ad hoc crack tracking in brosse breuksimulaties elimineert en eindige-vervormingsplasticiteit ondersteunt binnen een meshfree framework.

De auteurs merken expliciet de beperkingen van de huidige release op:

Het is beperkt tot isotrope materiaalrespons (isotrope hyperelasticiteit, isotrope $J_2$ plasticiteit).
Plasticiteit en breuk zijn niet gekoppeld; ductiele breukmechanismen zijn niet inbegrepen.
Anisotrope constitutieve modellen, composietmaterialen en specifieke hourglass-control operators zijn nog niet geïmplementeerd.

Ondanks deze beperkingen legt dit werk een robuuste fundering voor toekomstige ontwikkelingen in de vaste stofmechanica en breukmodellering binnen de SPH-gemeenschap, en biedt het een reproduceerbaar en uitbreidbaar platform voor onderzoekers.