A blended approach for evolving phase fields using peridynamics: Cyclic loading in quasi-brittle fracture

Dit artikel introduceert een nieuw veldtheoretisch model dat peridynamica combineert met een tweepunts-historieafhankelijke faseveldbenadering om schade en breuk in kwasibrosse materialen onder cyclische belasting te voorspellen, waarbij zowel plastische vervorming als degradatie van de elasticiteitsmodulus wordt meegenomen en de theorie voldoet aan de thermodynamische wetten.

De kern

Stel je voor dat je een stuk beton bekijkt, zoals een brug of een muur. Als je er te hard op duwt, ontstaan er barstjes. Vroeger waren computersimulaties van dit proces heel lastig, vooral als het materiaal niet alleen breekt, maar ook een beetje "plastisch" vervormt (zoals een stukje metaal dat buigt en niet meer helemaal recht komt) en als je het materiaal heen en weer belast (cyclisch), zoals bij trillingen.

Deze paper introduceert een nieuwe, slimme manier om dit te simuleren. Ze noemen het een "geblende" (samengevoegde) aanpak. Hier is hoe het werkt, vertaald naar alledaagse taal:

1. Het idee: Een web van elastiekjes

Stel je het materiaal voor als een enorm groot web van elastiekjes die aan elkaar zijn geknoopt.

• De oude manier: Vaak gebruikten wetenschappers twee aparte regels: één voor hoe het materiaal beweegt (Newton's wetten) en één aparte, ingewikkelde regel om te zeggen "hier is nu een barst". Dit was als het hebben van twee verschillende bestuurders voor dezelfde auto.
• De nieuwe manier: In deze paper gebruiken ze één enkele regel (Newton's tweede wet) voor alles. Ze zeggen: "Als twee punten in het web te ver uit elkaar worden getrokken, wordt het elastiekje zwakker en breekt het." Er is geen aparte regel voor de barst; de barst ontstaat vanzelf uit het gedrag van de elastiekjes.

2. De "Geheugen-Elastiekjes" (Peridynamica)

Normale elastiekjes in een computermodel kijken alleen naar wat er nu gebeurt. Deze nieuwe elastiekjes hebben echter geheugen.

• Als je een elastiekje uitrekt en het laat los, onthoudt het: "Ik ben al een keer uitgerekt, ik ben nu een beetje versleten."
• Dit is cruciaal voor materialen zoals beton. Als je beton belast, ontstaan er microscopische barstjes. Als je de last wegneemt, komt het niet helemaal terug in zijn originele vorm. Het materiaal is "vermoeid".
• De simulatie houdt bij hoe ver elk elastiekje ooit is uitgerekt (de "geschiedenis"). Als het te vaak of te ver is uitgerekt, verliest het zijn kracht.

3. Twee soorten elastiekjes: Trek en Druk

Materiaal gedraagt zich anders als je erop trekt dan als je erop duwt.

• Trek (Tension): Als je aan een elastiekje trekt en het is beschadigd, breekt het. Het kan de kracht niet meer dragen.
• Druk (Compression): Als je op een beschadigd elastiekje duwt, werkt het nog steeds! Het kan nog steeds kracht overbrengen. Denk aan een gebroken baksteen: als je erop duwt, blijft hij staan, maar als je eraan trekt, valt hij uit elkaar.
• De nieuwe methode gebruikt twee aparte "geheugen-registers": één voor trek en één voor druk. Dit zorgt ervoor dat de simulatie heel realistisch is voor cyclische belasting (heen-en-weer bewegen).

4. De "Verdwijnende Kracht" (De Barst)

Hoe ziet een barst eruit in dit model?

• In de oude methodes moest je een lijn tekenen waar de barst zou komen.
• In deze methode gebeurt het vanzelf. Als de elastiekjes tussen twee punten te zwak worden, verdwijnt de kracht tussen hen.
• Als je kijkt naar een groot aantal van deze "gebroken" elastiekjes die in een lijn liggen, zie je plotseling een barst ontstaan. De computer hoeft niet te weten waar de barst komt; hij volgt gewoon de regels van de elastiekjes en de barst "ontdekt zichzelf".

5. Waarom is dit belangrijk? (De "Grootte-effect" en Energie)

• Energiebalans: Het model houdt perfect rekening met energie. De energie die je in het materiaal stopt, gaat naar beweging, elasticiteit of wordt "verslind" door het maken van nieuwe barstoppervlakken. Dit voldoet aan de natuurwetten.
• Grootte-effect: Als je een klein betonblok en een groot betonblok hebt, breekt het grote vaak relatief "makkelijker" dan je zou denken op basis van de sterkte van het materiaal alleen. Dit model pakt dit fenomeen automatisch op zonder dat je handmatig regels moet toevoegen.
• Mesh-vrij: Normaal gesproken moet je een model indelen in een rooster (een raster van vierkantjes). Als een barst door een vierkantje loopt, moet je het rooster opnieuw tekenen. Deze methode is "mesh-vrij". Het werkt met losse punten die met elkaar communiceren. Als er een barst komt, hoeft het rooster niet te worden aangepast; de punten communiceren gewoon niet meer met elkaar.

Samenvattend

De auteurs hebben een manier bedacht om brekende materialen te simuleren die werkt als een gigantisch, slim web van elastiekjes met geheugen.

• Ze gebruiken één simpele wet (beweging) in plaats van twee complexe regels.
• Ze houden geschiedenis bij (hoeveel is het materiaal al vervormd?).
• Ze onderscheiden tussen trekken en duwen.
• De barst ontstaat vanzelf uit het gedrag van de elastiekjes.

Dit maakt het mogelijk om heel nauwkeurig te voorspellen hoe beton en andere kwetsbare materialen falen onder zware of herhaalde belasting, zonder ingewikkelde handmatige aanpassingen. Het is alsof je een virtueel laboratorium hebt waar je kunt zien hoe een brug langzaam uit elkaar valt, punt voor punt, volledig gebaseerd op de natuurwetten.

Technische samenvatting

Titel: Een geblokte aanpak voor het evolueren van fasevelden met behulp van peridynamica: Cyclische belasting in quasi-sprongbrekende breuk

1. Het Probleem

Het identificeren van een geschikte veldtheorie voor het formuleren van het initiële randwaardeprobleem (IBVP) voor breuk in continue media is een langdurige uitdaging. Traditionele methoden voor het modelleren van schade en breuk in quasi-sprongbrekende materialen (zoals beton) hebben vaak te kampen met de volgende beperkingen:

• Koppeling van vergelijkingen: Traditionele faseveldmethoden koppelen een tijdsafhankelijke faseveldvergelijking (vaak gebaseerd op de Ginzburg-Landau vergelijking) aan een bewegingsvergelijking voor verplaatsing (Newton's tweede wet). Dit vereist vaak extra aannames of variatiele beperkingen.
• Plasticiteit en hysterese: Het modelleren van irreversibele vervorming (plasticiteit) en hysterese bij cyclische belasting binnen een dynamisch breukmodel is complex. Veel bestaande modellen kunnen de degradatie van de elasticiteitsmodulus in coördinatie met de accumulatie van irreversibele rek niet adequaat beschrijven.
• Mesh-afhankelijkheid: Veel numerieke methoden zijn mesh-afhankelijk, wat leidt tot artefacten in de voorspelling van kraspatronen.

De auteurs streven naar een theorie die uitsluitend gebaseerd is op Newton's tweede wet, maar die toch in staat is om schade, plasticiteit en breuk te lokaliseren zonder een aparte evolutievergelijking voor het faseveld.

2. Methodologie

De auteurs presenteren een geblendeerde (blended) aanpak die peridynamica combineert met een tweepunts geschiedenis-afhankelijk faseveld. De kern van de methodologie omvat:

• Peridynamische Basis: De verplaatsing van het materiaal evolueert volgens Newton's tweede wet, waarbij krachten worden berekend via een niet-lokale constitutieve wet. De interactie tussen punten $x$ en $y$ binnen een horizon $\epsilon$ wordt beschreven door een "bond" (verbinding).
• Tweepunts Fasevelden: In plaats van een apart veld voor schade, wordt het faseveld $\gamma$ geïntroduceerd als een functie van de tweepunts rek tussen $x$ en $y$. Er worden twee onafhankelijke fasevelden gebruikt:
  • $\gamma^+$ voor trekkrachten (tensile).
  • $\gamma^-$ voor drukkrachten (compressive).
    Deze velden variëren van 1 (onbeschadigd) tot 0 (volledig beschadigd).
• Constitutieve Wet met Geheugen: De spanning-rek relatie omvat:
  • Een sterkteomhulsel (strength envelope) dat het gedrag bij belasting beschrijft (lineair elastisch tot vloeigrens, dan niet-lineair met verharding).
  • Een familie van ontlastingswetten (unloading laws) gebaseerd op schade en plasticiteit.
  • Irreversibele rek: De totale rek wordt opgesplitst in elastische en irreversibele (plastic) rek. De plasticiteit wordt gemodelleerd via een ontlastingsverhouding ($\beta$), die de verhouding aangeeft tussen plastic rek met materiaalverzachting en plastic rek zonder verzachting.
• Niet-lokale en Tijdsafhankelijke Eigenschappen: De niet-lokale interactie is zowel ruimtelijk (via de horizon $\epsilon$) als tijdsafhankelijk (via de maximale historische rek $S^*$).
• Mesh-vrije Numerieke Implementatie: De IBVP wordt opgelost met een mesh-vrije methode gebaseerd op kwadratuur, waarbij geen elementen nodig zijn tussen de knooppunten.

3. Belangrijkste Bijdragen

• Unieke Formulering: De paper toont aan dat materialdegradatiemodellen die consistent zijn met plastische dissipatie, volledig kunnen worden beschreven door een tweepunts geschiedenis-afhankelijk faseveld, zonder een aparte evolutievergelijking voor het faseveld.
• Thermodynamische Consistentie: De theorie voldoet aan de energiebalans. De dissipatiesnelheid is positief en volgt direct uit de bewegingsvergelijkingen en de constitutieve wetten, in overeenstemming met de wetten van de thermodynamica.
• Herstel van Breukenergie: De breukenergie voor vlakke krasen wordt direct uit het model afgeleid en komt overeen met het product van de energieafgifte-snelheid en het krasoppervlak (Griffith-theorie). Dit resultaat is onafhankelijk van de schaal van de niet-lokale interactie.
• Unieke Oplossing: Er wordt bewezen dat het initiële randwaardeprobleem een unieke oplossing heeft voor verplaatsing, elastische spanning en zowel elastische als irreversibele rekvelden, mits gebruik wordt gemaakt van de vastgestelde ontlastingswetten.
• Calibratie: De modelparameters (elasticiteitsmodulus, sterkte, kritieke energieafgifte) worden gekalibreerd op basis van $\Gamma$-convergentie en geometrische maattheorie, waarbij de niet-lokale schaal $\epsilon$ niet als een materiaaleigenschap wordt behandeld, maar als een discretisatieparameter.

4. Resultaten

De auteurs hebben de methode getest op diverse benchmarks en vergelijkingen met experimentele data:

• Mesh Convergentie: Er is aangetoond dat het model convergeert wanneer de verhouding tussen horizon en roosterafstand ($m$) toeneemt.
• Cyclische Belasting (Beton): Simulaties van driedubbele punt-buigtests op beton tonen een uitstekende overeenkomst (zowel kwalitatief als kwantitatief) met experimentele data, inclusief de hysterese in de krasmondopening (CMOD) versus kracht. Het model vangt de schade-accumulatie en plastische vervorming correct op.
• Gekoppelde Modus Breuk: Voor gemengde modus breuk (trek + schuif) in genotste betonbalken en L-vormige proefstukken, voorspelt het model de kraspaden nauwkeurig, in overeenstemming met experimentele observaties.
• Grootte-effect (Size Effect): Het model reproduceert het deterministische grootte-effect in betonstralen (waarbij grotere balken een lagere nominale sterkte vertonen) zonder expliciete regels voor sterkte-aanpassing op basis van de proefstukgrootte. Dit komt voort uit de initiële randwaardeprobleem-oplossing zelf.
• Vergelijking met andere methoden: De resultaten komen overeen met die van andere geavanceerde peridynamische modellen en experimentele literatuur, maar met een vereenvoudigde constitutieve benadering.

5. Significantie

Deze studie is significant omdat het een wiskundig goed gesteld raamwerk biedt voor het modelleren van dynamische schade en breuk in quasi-sprongbrekende materialen met plastisch gedrag.

• Eenvoud en Efficiëntie: Door te werken met één bewegingsvergelijking en geschiedenis-afhankelijke constitutieve wetten, wordt de complexiteit van gekoppelde systemen vermeden.
• Fysische Betekenis: Het model koppelt macroscopische observaties (zoals hysterese en grootte-effecten) direct aan microscopische mechanismen (plastische rek en materiaalverzachting) via de parameter $\beta$.
• Toepasbaarheid: De methode is breed toepasbaar op geomaterialen, metselwerk en beton, en biedt een solide basis voor toekomstige uitbreidingen naar anisotrope media en composietmaterialen.
• Validatie: De sterke overeenkomst met experimentele data voor zowel monotone als cyclische belasting bevestigt de praktische bruikbaarheid van de geblokte peridynamische-faseveld aanpak voor engineering-toepassingen.

Kortom, de paper levert een robuust, thermodynamisch consistent en numeriek efficiënt model dat de complexiteit van quasi-sprongbrekende breuk onder cyclische belasting succesvol simuleert zonder de beperkingen van traditionele mesh-gebaseerde of gekoppelde faseveldmethoden.