On modeling fracture of soft polymers

De auteurs presenteren een damage-model gebaseerd op een kritieke waarde van stresswerk dat de breukweerstand van verschillende soorten zachte polymeren, variërend van elastomeer tot viskeus materiaal, onder diverse belastingcondities en geometrieën succesvol voorspelt.

De kern

Stel je voor dat je een stukje kauwgom, een zachte gel of een elastiekje hebt. Soms trek je er heel langzaam aan en rekt het zich uit als een lange, plakkerige draad. Soms trek je er heel snel aan en pats, het breekt in tweeën.

Deze materialen heten zachte polymeren. Ze zitten overal: in onze huid, in medische implants, in zachte robots en in lijm. Maar er is een groot probleem: wetenschappers hebben tot nu toe geen goede manier gevonden om te voorspellen wanneer en hoe deze materialen breken, vooral niet als je ze snel of langzaam trekt.

Dit paper van Aditya Konale en Vikas Srivastava van Brown University lost dit probleem op. Hier is de uitleg in simpele taal, met een paar creatieve vergelijkingen.

1. Het probleem: De "breukkracht" is niet eenduidig

Voor harde materialen (zoals staal) weten we precies hoeveel kracht nodig is om ze te breken. Voor zachte, plakkerige materialen werkt dat niet.

• De analogie: Stel je voor dat je een stuk taart hebt. Als je er heel langzaam een mes doorheen duwt, zakt het mes erin en breekt de taart niet. Duw je heel snel, dan breekt de taart misschien wel.
• De huidige meetmethodes zeggen: "Het hangt ervan af hoe snel je trekt en hoe het stuk eruit ziet." Dat is vervelend voor ingenieurs die iets willen bouwen dat veilig is, ongeacht hoe snel je erop duwt.

2. De oplossing: Een nieuwe "breuk-meter"

De auteurs hebben een nieuwe manier bedacht om te meten hoeveel energie er nodig is om een punt in het materiaal te laten breken. Ze noemen dit $W_{cr}$ (kritieke spanningsarbeid).

• De analogie: Denk aan een emmer die je vult met water (energie).
  • Bij een hard materiaal is de emmer klein en breekt hij snel.
  • Bij een zacht materiaal is de emmer groot en kan hij veel water opnemen voordat hij lekt.
  • De oude methodes keken naar de vorm van de emmer. De auteurs zeggen: "Nee, kijk gewoon naar het totaal aantal liters water dat erin past voordat hij breekt."
  • Ze ontdekten dat voor een bepaald type zacht materiaal, dit "totaal aantal liters" (de energie) bijna altijd hetzelfde is, of je nu langzaam of snel trekt. Dit is hun nieuwe, betrouwbare maatstaf.

3. Het model: Een auto met twee versnellingen

Om te voorspellen hoe het breken verloopt, hebben ze een computermodel gemaakt. Ze zien het materiaal niet als één blok, maar als een systeem met twee mechanismen die samenwerken:

1. Het elastische netwerk (De veer): Dit is het skelet van het materiaal. Het werkt als een veer die terugveert.
2. De dynamische verbindingen (De lijm): In veel zachte polymeren zijn de verbindingen tussen de moleculen niet permanent. Ze breken en plakken weer aan elkaar, net als magneetjes die loslaten en weer vastklikken.

• De analogie: Stel je een touw voor dat gemaakt is van duizenden kleine magneetjes.
  • Trek je langzaam: De magneetjes hebben tijd om los te laten en weer vast te klikken op een nieuwe plek. Het touw wordt lang en dun, maar breekt niet snel. Het gedraagt zich als een vloeistof.
  • Trek je snel: De magneetjes hebben geen tijd om te bewegen. Ze worden allemaal tegelijk uitgerekt en breken. Het touw knapt plotseling. Het gedraagt zich als een elastiek.

Het model van de auteurs houdt rekening met beide: de veerkracht én de tijd die de magneetjes nodig hebben om te bewegen.

4. Wat kan dit model?

Dit nieuwe model is als een "Crystal Ball" voor zachte materialen. Het kan voorspellen:

• Het "Trompet"-effect: Bij zeer snelle trekking van een zacht, plakkerig materiaal (zoals Polyborosiloxane of PBS), ziet de breuklijn eruit als een trompet. Eerst wordt het gat rond, dan wordt er een tweede, kleiner gat voor gevormd voordat het helemaal breekt. Het model ziet dit precies zo.
• De "Parabool"-effect: Bij elastische materialen (zoals rubber) is de breuklijn gewoon een mooie, enkele boog (een parabool).
• Het stoppen van een scheur: Soms begint een scheur, maar stopt hij halverwege. Het model kan dit ook voorspellen, afhankelijk van hoe snel je trekt.

5. Waarom is dit belangrijk?

Tot nu toe konden wetenschappers alleen kijken naar materialen die al een scheurtje hadden. Dit model kan voorspellen waar een scheur ontstaat in een compleet nieuw stuk materiaal, zelfs als het een complexe vorm heeft.

Bovendien kunnen ze met dit grote model terugrekenen naar de microscopische wereld. Ze hebben bijvoorbeeld berekend hoeveel energie nodig is om één enkele "magneet" (een moleculaire binding) te breken. Toen ze dit in het lab testten met een ander experiment, klopte het precies!

Samenvatting in één zin

De auteurs hebben een nieuwe, universele "breuk-meter" bedacht en een slim computermodel gemaakt dat precies kan voorspellen hoe zachte, plakkerige materialen breken – of je ze nu langzaam uitrekt als een kauwgom of snel als een rubberen band – en dit werkt voor zowel natuurlijke weefsels als kunstmatige lijm.

Technische samenvatting

Titel: Op modellering van breuk in zachte polymeren

Auteurs: Aditya Konale en Vikas Srivastava (Brown University)

---

1. Het Probleem

Zachte polymeren (zoals elastomeer, viskeuze polymeren, polymergele en biologische weefsels) zijn veelvoorkomende materialen met eigenschappen die variëren van snelheidsonafhankelijk (elastisch) tot snelheidsafhankelijk (viskeus). Een groot probleem in de huidige mechanica is het ontbreken van een universeel maatstaf voor de breuktaaiheid van deze materialen:

• Niet-unieke maatstaven: Bestaande maatstaven voor breuktaaiheid (zoals de energie-afgifte-rate of spanningsintensiteit aan de scheurtip) zijn niet uniek voor snelheidsafhankelijke materialen. Ze variëren sterk afhankelijk van de belastinggeschiedenis en de geometrie van het proefstuk.
• Beperkte modellen: Bestaande modellen voor het modelleren van breuk zijn vaak beperkt tot specifieke voorgekraakte geometrieën of alleen van toepassing op elastomeren (permanente netwerken). Er ontbreekt een algemeen model dat zowel elastomeren als viskeuze, dynamisch verweven polymeren (zoals vitrimeren en polyborosiloxanen) kan beschrijven, inclusief het voorspellen van het ontstaan (nucleatie), de groei en de volledige breuk in willekeurige geometrieën.
• Complexiteit van het fenomeen: Zachte polymeren vertonen complexe fenomenen zoals notch-verdoving (blunting), overgangen van ductiel naar bros gedrag, en specifieke breukprofielen (bijv. "trompet"-vormige profielen bij viskeuze materialen versus paraboolvormige profielen bij elastomeren).

2. Methodologie

De auteurs stellen een nieuw, unificerend raamwerk voor dat gebaseerd is op de methode van virtuele macht en een gegeneraliseerd gradiënt-schade-model.

• Kernconcept: Kritieke Spanningsarbeid ($W_{cr}$):
  De auteurs introduceren $W_{cr}$ als een uniek fysiek parameter om de breukweerstand te kwantificeren. Dit is gedefinieerd als de totale energie per eenheid referentievolume die nodig is voor de volledige falen van een willekeurig materiaalpunt onder voldoende snelle belasting. In tegenstelling tot traditionele maatstaven (J/m²), heeft $W_{cr}$ de eenheid J/m³ en wordt geacht onafhankelijk te zijn van de proefstukgeometrie en de belastingsnelheid (binnen een bepaalde klasse van materialen).

• Schade-initiatie en -groei:

  • Initiatie: Schade begint wanneer de energetische bijdrage aan de spanningsarbeid ($\psi^+_{cr}$) een kritieke drempelwaarde bereikt. Dit criterium houdt rekening met de materiaals sterkte.
  • Groei: De groei van schade wordt gemodelleerd via een gradiënt-schade raamwerk. Dit zorgt voor een diffuse schadezone en voorkomt numerieke instabiliteiten. Een kinetische modulus ($\zeta$) bepaalt de tijdschaal van de schadegroei.

• Constitutief Model (Multi-mechanisme):
  Het model gebruikt een multiplicatieve decompositie van de vervormingsgradiënt en beschouwt twee parallelle mechanismen:

  1. Intermoleculair weerstand (Mechanisme A): Beschrijft elastisch-viskoplastisch gedrag en intermoleculaire interacties.
  2. Netwerkweerstand (Mechanisme B): Beschrijft het rekken van subketens binnen dynamische, reversibele crosslinks (zoals B:O-coördinatiebindingen). Dit mechanisme houdt rekening met de kinetiek van crosslink-ontbinding en -herformatie.
     De vrije-energie-dichtheid wordt zo geformuleerd dat deze thermodynamisch consistent is en geen expliciete snelheidsafhankelijkheid bevat in de vrije energie zelf (in plaats daarvan wordt dit via de kinetiek van de crosslinks gemodelleerd).

• Numerieke Implementatie:
  Het model is geïmplementeerd in de finite-element software ABAQUS via een expliciete formulering. Dit vereiste de ontwikkeling van een VUEL (User-Defined Element) voor de schadevariabele en een VUMAT voor het materiaalgedrag, waarbij de gegevens tussen beide worden uitgewisseld.

3. Belangrijkste Bijdragen

1. Unieke Breukweerstand: Het aantonen dat de kritieke spanningsarbeid ($W_{cr}$) een bijna unieke maatstaf is voor de breukweerstand van een bepaalde klasse zachte polymeren, ongeacht de belastingsnelheid of geometrie.
2. Unificerend Model: De ontwikkeling van een enkel model dat zowel elastomeren (rate-onafhankelijk) als viskeuze polymeren (rate-afhankelijk) kan simuleren, inclusief het voorspellen van schade-initiatie, groei en volledige breuk.
3. Microstructuur-koppeling: Het vermogen om microstructurele grootheden (zoals de dissociatie-energie van subketens) af te leiden uit macroscopische modelparameters.
4. Voorspelling van Fysische Fenomenen: Het model kan complexe experimentele observaties voorspellen, zoals de overgang van een "trompet"-vormig breukprofiel (bij viskeuze materialen) naar een enkele paraboolvorm (bij elastomeren), en het fenomeen van notch-verdoving.

4. Resultaten en Validatie

Het model is gevalideerd tegen experimentele data voor een breed scala aan materialen en belastingcondities:

• Polyborosiloxaan (PBS):

  • Het model voorspelde nauwkeurig de kracht-verplaatsingscurves voor verschillende geometrieën (gekraakt, met inkepingen, met gaten) en belastingsnelheden (2 mm/s tot 60 mm/s).
  • Het voorspelde correct het verschil in gedrag: bij hoge snelheden (60 mm/s) treedt breuk op met lage rekbaarheid en een trompet-profiel; bij lage snelheden (6 mm/s) vertoont het materiaal hoge rekbaarheid zonder breuk (door netwerk-reconfiguratie).
  • Microstructurele validatie: Uit de parameter $\psi^+_{cr}$ werd de dissociatie-energie van subketens geschat op ~15 kJ/mol. Dit werd bevestigd door onafhankelijke oscillatoire reologie-experimenten, waarbij de activeringsenergie voor de dissociatie werd bepaald op ~24 kJ/mol (zelfde orde van grootte).

• EPDM (Elastomeer):

  • Het model voorspelde de breukverplaatsing voor V-vormige inkepingen met verschillende dieptes en stralen, in overeenstemming met experimentele data.

• EPS25 Vitrimer (Dynamisch verweven):

  • Het model voorspelde de overgang van geen scheurgroei (bij lage snelheid) naar scheurgroei en vervolgens scheurgroei-stops (arrestatie) bij intermediaire snelheden, wat overeenkomt met experimentele waarnemingen.

• SBR en NR (Styreen-butadieen en Natuurlijk Rubber):

  • Het model voorspelde de faaloppervlakken (failure surfaces) onder multiaxiale belasting nauwkeurig, in tegenstelling tot variatieve faseveldmodellen die hier vaak onfysische resultaten opleveren.

• Viskeuze vs. Elastische Effecten:

  • Simulaties waarbij viskeuze effecten werden verwijderd, toonden een duidelijke overgang in het breukprofiel van een trompet-vorm (viskeus) naar een parabool-vorm (elastisch), wat de rol van snelheidsafhankelijkheid in het breukmechanisme bevestigt.

5. Betekenis en Conclusie

Dit werk biedt een aanzienlijke stap voorwaarts in de mechanica van zachte materialen door een universeel, continuüm-gebaseerd model te bieden dat de complexiteit van rate-afhankelijke breuk adresseert.

• Praktische Toepassing: Het model kan worden gebruikt voor het ontwerpen van veiligheidskritieke componenten van zachte polymeren (bijv. in biomedische toepassingen, impactbescherming en flexibele elektronica) zonder te hoeven vertrouwen op specifieke proefstukgeometrieën.
• Fundamenteel Inzicht: Het koppelt macroscopische mechanische respons direct aan microstructurele processen (zoals crosslink-dissociatie), wat inzicht geeft in de fundamentele oorzaken van falen.
• Toekomstperspectief: De auteurs suggereren dat $W_{cr}$ kan dienen als een robuuste, geometrie-onafhankelijke eigenschap voor het karakteriseren van nieuwe zachte materialen, mits de belastingssnelheden voldoende hoog zijn om schade-initiatie te veroorzaken.

Samenvattend biedt deze studie een krachtige, thermodynamisch consistente en computationally haalbare oplossing voor het voorspellen van breuk in een breed spectrum van zachte polymeren, van volledig elastisch tot hoogviskeus.