Phase-Field Peridynamics

Titel: Phase-Field Peridynamics

Auteurs: Kai Partmann, Christian Wieners, Michael Ortiz, en Kerstin Weinberg.

1. Het Probleem

De modellering van breuk (fractuur) in structurele componenten is een fundamentele uitdaging in de computationele mechanica. Traditionele methoden zoals de Finite Element Methode (FEM) vereisen expliciete mesh-herindeling of speciale technieken om discontinuïteiten (scheuren) te behandelen, wat rekenkundig kostbaar en complex is.

Peridynamiek (PD) lost dit op door de balans van lineaire impuls te formuleren als een integro-differentiaalvergelijking, waardoor geen ruimtelijke afgeleiden van het verplaatsingsveld nodig zijn. Dit maakt PD van nature geschikt voor het modelleren van breuk. Echter, de meest nauwkeurige formulering, de bond-associated correspondence formulation, lijdt onder een specifiek probleem bij het modelleren van breuk:

• Stabiliteitsproblemen bij bond-deletie: In traditionele PD-modellen wordt een "bond" (interactie tussen twee punten) permanent verwijderd zodra een kritieke vervorming wordt overschreden.
• Verslechtering van de kinematische benadering: In correspondence-formuleringen wordt de vervormingsgradiënt benaderd via een vormtensor (shape tensor) die afhankelijk is van de configuratie van de buren. Het progressief verwijderen van buren rondom een voortplantende scheur degradeert deze vormtensor, wat leidt tot slecht geconditioneerde of singuliere matrices.
• Gevolg: Dit veroorzaakt numerieke instabiliteiten, vooral bij complexe belastingen zoals schuifkrachten (Mode II), waardoor de meest nauwkeurige PD-formulering onbetrouwbaar wordt voor breukproblemen.

2. Methodologie

De auteurs introduceren een nieuw Phase-Field Peridynamics (PFPD) kader dat de bond-associated correspondence-formulering combineert met een faseveld-benadering voor breuk, zonder bonds te verwijderen.

Kerncomponenten van de methode:

• Continuïteit in plaats van deletie: In plaats van bonds abrupt te verwijderen, wordt de energetische bijdrage van elke bond continu gedegradeerd via een bond-faseveldparameter $s \in [0, 1]$.
  • $s=0$: Onbeschadigde bond.
  • $s=1$: Volledig beschadigde bond.
• Twee degradatiefuncties:
  1. Energetische degradatie ($g(s)$): Vermindert de spanning (kracht) die de bond kan overbrengen. Dit volgt een kwadratische wet $g(s) = (1-s)^2$.
  2. Kinematische degradatie ($h(s)$): Een gescheiden functie die de bijdrage van de bond aan de vormtensor (shape tensor) en de kinematische benadering reguleert. Deze functie blijft gelijk aan 1 tot een drempelwaarde $s_c$ wordt bereikt, en degradeert daarna pas.
  • Doel: Zelfs als een bond energetisch "dood" is ($s \approx 1$), blijft hij kinematisch actief tot $s_c$ wordt overschreden. Dit behoudt de stabiliteit van de vormtensor en voorkomt numerieke instabiliteiten.
• Evolutie van schade: De schadeontwikkeling wordt afgeleid uit een variationeel principe dat de totale energie minimaliseert. Een geschiedenisvariabele zorgt voor irreversibiliteit (schade kan niet herstellen).
• Normalisatieconstante ($c_0$): Een cruciaal aspect is het analytisch afleiden van de normalisatieconstante voor de faseveldenergie. Deze moet ervoor zorgen dat de totale energie-dissipatie overeenkomt met de Griffith-breuktheorie (energievrijgave per oppervlakte-eenheid).
  • De auteurs leiden een exacte, gesloten vorm af voor algemene sferische kernfuncties als een verhouding van twee één-dimensionale integralen.

3. Belangrijkste Bijdragen

1. Nieuw PFPD Kader: Een robuust model dat de nauwkeurigheid van de bond-associated correspondence-formulering behoudt terwijl het de instabiliteiten van bond-deletie elimineert.
2. Kinematische Degradatie: De introductie van een gescheiden kinematische degradatiefunctie ($h(s)$) die de stabiliteit van de vormtensor garandeert, zelfs bij uitgebreide schade.
3. Analytische Normalisatie: Een wiskundig exacte afleiding van de normalisatieconstante $c_0$ voor willekeurige sferische kernfuncties, wat de thermodynamische consistentie met Griffith's theorie waarborgt.
4. Algoritmische Implementatie: Een gedetailleerde vergelijking van het algoritme met de standaard bond-based benadering, waarbij wordt aangetoond dat PFPD geen binaire deletie vereist maar continue degradatie.

4. Resultaten

De methode is gevalideerd via diverse numerieke voorbeelden:

• Mode I en Mode II Breuk:
  • Het model voorspelde correcte scheurpaden voor zowel trek- als schuifbelasting.
  • Cruciaal resultaat: Bij Mode II (schuif) werd de standaard BAQP-model met bond-deletie instabiel en faalde de simulatie. Het PFPD-model bleef stabiel en leverde accurate resultaten, dankzij de behouden kinematische benadering.
• Boundary Tension Test (BTT):
  • Gebruikt om de afgeleide normalisatieconstante te verifiëren.
  • Simulaties met verschillende horizon-ratio's ($m_\delta$) en kernfuncties (kubische B-spline en lineair) leverden identieke scheurvertakkingspatronen op. Dit bevestigt dat de normalisatie correct is en onafhankelijk van de discretisatieparameters.
  • De scheurtip-snelheid bleef onder de theoretische limiet van de helft van de Rayleigh-golfsnelheid ($0.5 c_R$), wat fysiek correct is.
• Kalthoff-Winkler Experiment:
  • Een simulatie van impactbreuk in een stalen plaat.
  • Het model voorspelde een scheurhoek van ongeveer 70°, wat overeenkomt met experimentele waarnemingen.
  • De snelheidsprofielen waren gladder dan bij traditionele modellen (zoals OSB), en de initiële scheurgroei was fysisch consistenter.

5. Betekenis en Conclusie

Deze studie biedt een doorbraak in de toepassing van peridynamiek op complexe breukproblemen.

• Stabiliteit: Door bonds niet te verwijderen maar continu te degraderen, wordt de vormtensor stabiel gehouden, wat de toepassing van de meest nauwkeurige PD-formulering (bond-associated correspondence) op breukproblemen mogelijk maakt.
• Thermodynamische Consistentie: De analytisch afgeleide normalisatieconstante zorgt ervoor dat het model strikt voldoet aan de energieprincipes van Griffith.
• Toepasbaarheid: Het model is geschikt voor dynamische breukproblemen, inclusief scheurvertakking en complexe belastingscenario's, zonder de numerieke instabiliteiten die traditionele PD-methoden plagen.

Kortom, de PFPD-methode combineert de nauwkeurigheid van klassieke continuümmechanica (via correspondence) met de flexibiliteit van peridynamiek voor breuk, en lost het langdurige stabiliteitsprobleem op dat de toepassing van deze geavanceerde formuleringen beperkte.