Dynamic brittle fracture using Lip-field approach in an explicit dynamics context

Technische Samenvatting: Dynamische Brosser Breuk met de Lip-veld Benadering in een Expliciete Dynamica Context

Probleemstelling
Dit artikel behandelt de numerieke simulatie van dynamische brosse breuk in materialen die worden blootgesteld aan belasting met een hoge vervormingssnelheid, zoals bij impacts of explosies. Terwijl klassieke breukmechanica (bijv. Griffiths model) moeite heeft met scheurinitiatie en vertakking, zijn er moderne benaderingen zoals Cohesive Zone Models (CZM) en diffuse schade modellen ontwikkeld. Echter, lokale diffuse schade modellen lijden aan pathologische mesh-gevoeligheid, wat regularisatie vereist om een lengteschaal te introduceren en mesh-onafhankelijke resultaten te garanderen. Bestaande regularisatietechnieken, zoals Phase-Field en Thick Level Set (TLS) methoden, impliceren vaak niet-lokale energiefunctionalen die computationeel duur kunnen zijn. De specifieat geadresseerde uitdaging is het uitbreiden van de Lip-veld benadering—een variationele methode die lengteschalen introduceert via Lipschitz-restricties op het schadeveld—van de eerdere quasi-statische en eendimensionale dynamische toepassingen naar tweedimensionale dynamische breuk binnen een expliciet tijdintegratiekader.

Methodologie
De auteurs stellen een variationele formulering voor voor dynamische schademecanica gebaseerd op de minimalisatie van een incrementeel potentiaal. De kern van de methodologie omvat:

1. Variationele Formulering: Het mechanische probleem wordt gedefinieerd door het minimaliseren van een incrementeel potentiaal $F$ dat kinetische energie, vrije energie en dissipatiepotentiaal bevat. De toestandsvariabelen zijn verplaatsing ($u$) en een scalaire schadevariabele ($d$). Het constitutieve model maakt gebruik van generalized standard materials met specifieke energie-degradatie ($g(d)$) en verweking ($h(d)$) functies.
2. Lip-veld Regularisatie: In tegenstelling tot Phase-Field modellen die de energiefunctional modificeren met een gradiëntterm, houdt de Lip-veld benadering de energiefunctional lokaal (niet-geregulariseerd). In plaats daarvan worden Lipschitz-restricties opgelegd aan het schadeveld, waarbij vereist wordt dat de gradiënt van de schade begrensd is door een kritische waarde ($1/l$). Dit beperkt het schadeveld tot een specifieke Lipschitz-ruimte $\mathcal{L}$, wat effectief een lengteschaal $l$ introduceert zonder de lokale energiedichtheid te wijzigen.
3. Expliciet Gestaffeld Schema: Het artikel gebruikt een expliciet Newmark tijdintegratieschema (central difference methode). Een gestaffelde aanpak wordt gebruikt waarbij de verplaatsings- en schadevelden sequentieel worden opgelost. Voor een gegeven tijdstap wordt eerst de verplaatsing bijgewerkt, gevolgd door de schade-update.
4. Efficiënte Schade Solver: Een belangrijke algoritmische bijdrage is het gebruik van grenswaarde-schatting (afgeleid van eerder werk) om de schadeberekening te optimaliseren. Het proces omvat:
   • Het oplossen van een puur lokaal minimalisatieprobleem om een lokaal schadeveld ($d_{loc}$) te vinden.
   • Het projecteren van $d_{loc}$ naar boven- en ondergrenzen op basis van de Lipschitz-restrictie.
   • Het identificeren van regio's waar deze grenzen verschillen.
   • Het oplossen van het niet-lokale (Lipschitz-geconstraineerde) minimalisatieprobleem alleen in het kleine deel van het domein waar de grenzen niet gelijk zijn. In regio's waar de grenzen samenvallen, wordt de lokale oplossing geaccepteerd. Dit vermindert de computationele kosten aanzienlijk vergeleken met het oplossen van het niet-lokale probleem over het gehele domein.
5. Discretisatie: De ruimtelijke discretisatie maakt gebruik van lineaire driehoekige elementen voor het verplaatsingsveld. Een apart "lip-mesh" (geconstrueerd vanuit elementcentra) wordt gebruikt om het schadeveld op te slaan en te beperken, waardoor het aantal Lipschitz-restricties wordt geminimaliseerd.

Belangrijkste Bijdragen

• Extensie naar 2D Dynamica: Het artikel slaagt erin de Lip-veld benadering, die eerder gevalideerd is in 1D dynamische fragmentatie en quasi-statische gevallen, uit te breiden naar tweedimensionale dynamische breukproblemen.
• Integratie van Expliciete Dynamica: Het schetst een specifieke implementatie van de Lip-veld benadering met behulp van een expliciet tijdintegratieschema, wat computationeel efficiënt is voor problemen met hoge snelheid impacts.
• Computationele Efficiëntie via Grenzen: De auteurs demonstreren dat het gebruik van de grenswaarde-schatting de dure niet-lokale minimalisatie kan beperken tot kleine subdomeinen, waardoor de Lip-veld benadering computationeel concurrerend wordt voor dynamische simulaties.
• Variationeel Onderscheid: Het werk verheldert het onderscheid tussen de Lip-veld en Phase-field benaderingen: de Lip-veld benadering houdt een lokale energiefunctional met niet-lokale restricties, terwijl Phase-field een niet-lokale energiefunctional gebruikt.

Simulatieresultaten
Het voorgestelde model werd gevalideerd door middel van tweedimensionale numerieke voorbeelden:

1. Single Edge Notched Tension (SENT): Het model simuleerde succesvol dynamische scheurvoortplanting en vertakking. De resultaten toonden mesh-convergentie in potentiaal en gedissipeerde energie. Het scheurvertakkingsfenomeen werd gevangen zonder ad-hoc criteria, en de scheurtip-snelheden bleven over het algemeen onder de Rayleigh-golfsnelheid, behalve tijdens transiënte pieken tijdens vertakking. De resultaten waren vergelijkbaar met Phase-field en TLS simulaties uit de literatuur.
2. Invloed van Lastintensiteit: Door de toegepaste tractie te variëren, demonstreerde het model het vermogen om meerdere scheurvertakkingen te vangen, waarbij vertakking dichter bij de inkeping (notch) optrad onder hogere belastingen.
3. Kalthoff-Winkler Test: Het model simuleerde scheurinitiatie en voortplanting onder dynamische afschuivingsbelasting.
   • Bij lage impact-snelheden ($33$ m/s) voorspelde het een Mode-I scheurvoortplantingshoek van $\approx 62.4^\circ$, consistent met experimentele observaties van $\approx 70^\circ$.
   • Bij hoge impact-snelheden ($100$ m/s) ving het model opeenvolgende scheurvertakking en meerdere scheuren op.
   • Beperking: De auteurs merken op dat hoewel experimenten een overgang naar ductiele breuk laten zien (afschuivingsbreuk bij $\approx -10^\circ$) bij hoge snelheden, het huidige elastisch-schademodel dit ductiele gedrag niet vangt, omdat het geen plasticiteit bezit.

Betekenis en Claims
Het artikel claimt dat de Lip-veld benadering een robuust variationeel kader biedt voor dynamische breuk dat de wiskundige relevantie van het probleem behoudt terwijl noodzakelijke lengteschalen worden geïntroduceerd. De primaire betekenis ligt in de computationele efficiëntie gerealiseerd door de bounds-gebaseerde solver, die het oplossen van niet-lokale vergelijkingen globaal vermijdt. De auteurs stellen dat de benadering effectief complexe breukverschijnselen zoals scheurinitiatie, vertakking en coalescentie afhandelt zonder ad-hoc criteria.

Het artikel concludeert bescheiden dat hoewel de Lip-veld benadering effectief is voor brosse breuk en fragmentatie, het beperkt is tot elastisch-schade gedrag. Er wordt erkend dat het simuleren van ductiele breuk bij hoge vervormingssnelheden een uitbreiding naar elastisch-plastische schade modellen zou vereisen, wat buiten de scope van het huidige werk valt. De methode wordt gepresenteerd als een levensvatbaar alternatief voor Phase-field en TLS methoden, met name waar computationele efficiëntie in expliciete dynamica cruciaal is.