💻 computer science

Accurate simulation of delamination with a resin-rich layer-dependent penalty stiffness based on structural cohesive elements

Dit artikel verbetert een hogere-orde structureel cohesief element door de introductie van een harsrijke laag-afhankelijke penalty-stijfheidsformulering die onderscheid maakt tussen normale en afschuivingscomponenten, waardoor de mesh-convergentie en de nauwkeurigheid van compressie- en delaminatievoorspellingen in shell-gebaseerde simulaties aanzienlijk worden verbeterd in vergelijking met conventionele methoden.

Oorspronkelijke auteurs: Xiaopeng Ai, Christos Kassapoglou, Boyang Chen

Gepubliceerd 2026-08-07

📖 7 min leestijd🧠 Diepgaand

CC BY 4.0

Oorspronkelijke auteurs: Xiaopeng Ai, Christos Kassapoglou, Boyang Chen

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 huis bouwt van duizenden piepkleine, kleurrijke LEGO-steentjes die aan elkaar zijn gelijmd om een supersterke muur te maken. In de echte wereld bouwen ingenieurs soortgelijke structuren met lagen koolstofvezel en kunststof hars om vleugels voor vliegtuigen of onderdelen voor raketten te maken. Deze materialen zijn ongelooflijk sterk, maar ze hebben een geheim zwak punt: de lagen kunnen soms van elkaar loskomen, zoals pagina's in een boek die eruit worden gescheurd. Dit afbladderen wordt "delaminatie" genoemd, en als dit onverwacht gebeurt, kan de hele structuur falen. Om dit te voorkomen, gebruiken wetenschappers computerprogramma's om te simuleren hoe deze lagen uit elkaar kunnen raken. Echter, deze programma's worstelen vaak met een lastig probleem: ze moeten extreem gedetailleerd zijn (met miljoenen piepkleine digitale steentjes) om de wiskunde kloppend te krijgen, wat ervoor zorgt dat de computer eeuwig nodig heeft om de simulatie te draaien. Als de digitale steentjes te groot zijn, raakt de computer in de war en voorspelt hij dat de lagen op onmogelijke manieren op elkaar worden samengedrukt, wat leidt tot foutieve antwoorden.

Dit artikel pakt deze verwarring aan door een slimmere manier uit te vinden om de computer te vertellen hoe de lagen zich moeten gedragen wanneer ze uit elkaar worden getrokken. De auteurs, Xiaopeng Ai, Christos Kassapoglou en Boyang Chen, stellen een nieuwe wiskundige "regelset" voor voor de lijm tussen de lagen. In plaats van te gokken hoe stijf de lijm moet zijn op basis van de dikte van de hele muur, kijken zij naar de piepkleine, rommelige laag extra hars (de lijm) die van nature tussen de vezellagen zit. Door een slimme truc te gebruiken die gebaseerd is op hoe balken buigen, hebben zij precies uitgewerkt hoe de spanning door die dunne lijmlaag zich verspreidt. Hun nieuwe methode stelt de computer in staat om veel grotere, grovere digitale steentjes te gebruiken zonder dat de wiskunde fout gaat. Dit betekent dat ingenieurs complexe afbladderingsproblemen veel sneller en met grotere nauwkeurigheid kunnen simuleren, waardoor de vliegtuigen en raketten die van deze materialen zijn gemaakt, veilig en betrouwbaar zijn.

Het Probleem: De "Snotterige" Lijm-fout

Denk aan een composietmateriaal als een enorme sandwich gemaakt van vele lagen brood (de vezel) en een klein beetje extra boter (de hars) tussen elke plak. Wanneer de sandwich begint te scheuren, is de boterlaag de eerste die uitrekt en breekt. In het verleden behandelden computersimulaties deze boterlaag alsof het een eenvoudige, uniforme veer was. Het probleem is dat deze "veer" vaak te stijf of te zwak was, afhankelijk van hoe de computer de sandwich in digitale stukjes sneed.

Wanneer de computer grote digitale stukken gebruikte, maakte hij een specifieke fout: hij voorspelde dat de lagen vlak voor de scheur met enorme kracht op elkaar werden samengedrukt, zoals bij een auto-ongeluk, terwijl ze in werkelijkheid gewoon voorzichtig uit elkaar zouden moeten bladderen. Deze "snotterige" fout betekende dat ingenieurs, om een correct antwoord te krijgen, minuscule, piepkleine digitale stukjes moesten gebruiken, wat de simulaties dagen liet duren. Het was alsof je een vloeiende curve probeert te tekenen met alleen maar kleine, hoekige stappen; je had miljoenen stappen nodig om het er goed uit te laten zien.

Het Nieuwe Idee: Luisteren naar de Spanning

De auteurs realiseerden zich dat de oude manier van het berekenen van de "stijfheid" van de lijm te simpel was. Ze wisten dat de spanning (de interne kracht) niet gelijkmatig door de dikte van de sandwich verspreid wordt; het verandert naarmate je van de bovenste laag naar de onderste laag beweegt. Om dit op te lossen, keken ze niet alleen naar de lijm; ze keken naar de fysica van hoe een balk buigt.

Stel je voor dat je een dik boek buigt. De pagina's aan de buitenkant rekken uit, de pagina's aan de binnenkant worden samengedrukt, en de pagina's in het midden blijven neutraal. De auteurs gebruikten dezezelfde logica om precies in kaart te brengen hoe de spanning door de dikte van de harslaag verandert. Ze leidden een nieuwe formule af die werkt als een "slimme vertaler". Het neemt het complexe, veranderende spanningspatroon en zet dit om in één enkel, accuraat stijfheidsgetal dat de computer kan gebruiken.

Cruciaal is dat deze nieuwe formule werkt, zelfs wanneer de computer veel lagen van hetzelfde type groepeert in één grote "superlaag" (een equivalent enkel-laag model). Eerdere methoden vereisten dat de computer elke laag individueel modelleerde, wat traag en inefficiënt was. De nieuwe methode stelt de computer in staat om lagen samen te voegen zonder nauwkeurigheid te verliezen, waardoor de simulatie veel sneller kan draaien.

Het Bewijs: Het Testen van de Nieuwe Regelset

Om te zien of hun nieuwe idee daadwerkelijk werkte, voerde het team een reeks tests uit met standaard "benchmark"-problemen die ingenieurs gebruiken om hun wiskunde te controleren. Deze omvatten:

De Double Cantilever Beam (DCB): Een test waarbij je twee lagen uit elkaar trekt, zoals het openen van een boek.
De End-Notched Flexure (ENF): Een test waarbij je de lagen langs elkaar laat glijden.
De Mixed-Mode Bending (MMB): Een test die trekken en glijden combineert.
De Reinforced DCB (R-DCB): Een complexere test met extra lagen aan de buitenkant, wat een echt versterkte structuur simuleert.

In elke test vergeleken ze hun nieuwe "slimme" stijfheid met de oude "domme" stijfheid. De resultaten waren duidelijk:

Mesh Convergentie: Wanneer ze grove (grote) digitale meshes gebruikten, faalde de oude methode. Het gaf verschillende antwoorden afhankelijk van hoe groot de stukken waren, en het voorspelde dat de lagen veel te hard op elkaar werden samengedrukt. De nieuwe methode gaf echter hetzelfde accurate antwoord, of de digitale stukken nu groot of klein waren. Het was alsof je een kaart had die er net zo goed uitzag of je nu inzoomde of uitzoomde.
Nauwkeurigheid: In de DCB-test verminderde de nieuwe methode de fout in de belasting-verplaatsingscurve (hoeveel kracht er nodig is om de lagen te scheiden) aanzienlijk. Voor een mesh van 7,5 mm had de oude methode een fout van meer dan 50%, terwijl de nieuwe methode deze zeer laag hield.
Vergelijking met de Werkelijkheid: In de complexe R-DCB-test, die werd vergeleken met echte fysieke experimenten, voorspelde de nieuwe methode het pad van de scheur en de kracht die nodig is om deze te breken veel beter dan de oude methode. Specifiek voorspelde het de tweede piek in de breukkracht 14% lager dan de oude methode, wat veel beter overeenkwam met de gegevens uit de echte wereld.

Wat Dit Betekent

Het artikel beweert niet dat het elk probleem in het universum van composietmaterialen heeft opgelost. Het richt zich specifelijk op het verbeteren van de "penalty stiffness" voor structurele cohesieve elementen, wat de digitale hulpmiddelen zijn die worden gebruikt om afbladdering te simuleren. De auteurs laten via simulaties en vergelijkingen met experimenten zien dat hun nieuwe formule superieur is aan de conventionele methode.

Door deze nieuwe methode te gebruiken, kunnen ingenieurs delaminatie in complexe 3D-structuren simuleren met grovere meshes. Dit betekent dat ze accurate resultaten kunnen krijgen in een fractie van de tijd die het vroeger kostte. Het is alsof je upgradet van een trage, gepixelde videogame naar een vloeiende high-definition versie zonder dat je een supercomputer nodig hebt. De auteurs suggereren dat deze aanpak in de toekomst kan worden uitgebreid om ook andere soorten schade te voorspellen, zoals scheuren binnen de lagen zelf, maar voor nu hebben ze succesvol het "snotterige lijm"-probleem opgelost, waardoor het digitale testen van composietmaterialen sneller en betrouwbaarder is geworden.

Technische Samenvatting: Nauwkeurige Simulatie van Delaminatie met een op de Resin-Rich Layer gebaseerde Penalty Stiffness op basis van Structurele Cohesieve Elementen

Probleemstelling
Glasvezelversterkte polymeercomposieten (FRP) zijn cruciaal voor luchtvaartstructuren, maar het voorspellen van delaminatie blijft een uitdaging. Hoewel cohesive zone modellen (CZM) met gebruik van cohesieve elementen (CE) de standaard zijn voor het simuleren van delaminatie, lijden ze aan een bekende beperking: de vereiste dat elementgroottes kleiner zijn dan de cohesieve zone lengte (CZL) om hoge spanningsgradiënten te vangen, wat leidt tot hoge computationele kosten.

Om dit aan te pakken, zijn eerder hogere-orde structurele cohesieve elementen (gebaseerd op de Kirchhoff-Love schelptheorie) ontwikkeld om solide elementen te vervangen, wat de efficiëntie aanzienlijk verbeterde. Deze elementen vertrouwen echter op een penalty stiffness-formulering om het constitutieve gedrag te definiëren. Conventionele formuleringen (bijv. Turon et al.) definiëren de penalty stiffness op basis van de totale laminaatdikte en de uit-het-vlak modulus. Deze aanpak faalt bij het nauwkeurig voorspellen van door-de-dikte spanningverdelingen, waarbij de compressie vóór de scheurtip wordt overschat en het niet lukt om mesh-convergentie te bereiken.

Eerdere pogingen om de nauwkeurigheid te verbeteren omvatten resin-rich layer-gebaseerde formuleringen. Hoewel accuraat, vereisten deze een layer-wise modelleringsaanpak waarbij elke ply-interface expliciet gemodelleerd moest worden. Dit maakte de winst in computationele efficiëntie van equivalent single-layer (ESL) modellen, waarbij meerdere plies met dezelfde oriëntatie worden gehomogeniseerd tot één enkel schelelement, teniet. Er was behoefte aan een penalty stiffness-formulering die de fysieke nauwkeurigheid van de resin-rich layer aanpak behoudt, maar ook toepasbaar is op ESL-modellen.

Methodologie
De auteurs stellen een nieuwe penalty stiffness-formulering voor voor structurele cohesieve elementen die de kloof overbrugt tussen de nauwkeurigheid van de resin-rich layer en de efficiëntie van ESL. De kern van de methode bestaat uit het afleiden van analytische expressies voor de door-de-dikte verdelingen van de normale spanning uit-het-vlak ( $\sigma_z$ ) en de transversale afschuifspanning ( $\tau_{xz}$ ) op basis van balktheorie en de Kirchhoff-Love hypothese.

Afleiding van de Spanningsverdeling:
- Het 3D-schelpprobleem wordt gereduceerd tot een 2D-balkprobleem.
- Onder de Kirchhoff-Love aannames varieert de axiale spanning lineair, terwijl de transversale afschuifspanning een kwadratische verdeling volgt.
- Evenwichtsvergelijkingen worden gebruikt om een cubische verdeling voor de normale spanning uit-het-vlak ( $\sigma_z$ ) af te leiden.
- Randvoorwaarden voor specifieke belastingsmodi (Mode I/DCB en Mode II/ENF) worden toegepast om de coëfficiënten van deze spanningsverdelingen te bepalen.
Penalty Stiffness Formulering:
- De totale opening (Mode I) of verschuiving (Mode II) van het cohesieve element wordt berekend door de rekken over de dikte van het gehomogeniseerde laminaat te integreren.
- Rekken worden afgeleid uit de spanningsverdelingen en de eigenschappen van de resin-rich layer (modulus $E_{rr}$ of $G_{rr}$ en dikte $h_{rr}$ ).
- De penalty stiffness ( $K_n$ voor normaal, $K_s$ voor afschuiving) wordt geformuleerd als een functie van de eigenschappen van de resin-rich layer en de som van de spanningsratio's bij elk interface-laag relatief aan de maximale spanning.
- Dit maakt de formulering in staat om rekening te houden met de cumulatieve compliantie van alle resin-rich lagen binnen een gehomogeniseerde ply-block, in plaats van de block te behandelen als een enkele homogene laag.

Belangrijkste Bijdragen

Nieuwe Formulering: Een op de resin-rich layer afhankelijke penalty stiffness-formulering die onderscheid maakt tussen normale en afschuivingsstijfheid.
ESL Compatibiliteit: De methode breidt de resin-rich layer aanpak uit van layer-wise modellering naar het equivalent single-layer framework, wat het gebruik van grovere meshes mogelijk maakt zonder verlies van fysieke nauwkeurigheid.
Analytische Afleiding: De afleiding van door-de-dikte spanningsverdelingen met behulp van balktheorie om de specifieke spanningsratio's vereist voor de stijfheidsberekening te bepalen.
Modus-afhankelijke Toepassing: De formulering wordt toegepast op Mode I, Mode II en mixed-mode scenario's met de juiste randvoorwaarden voor DCB, ENF en Mixed-Mode Bending (MMB) configuraties.

Resultaten en Validatie
De voorgestelde methode is gevalideerd tegen vier benchmark-problemen: Double Cantilever Beam (DCB), End-Notched Flexure (ENF), Mixed-Mode Bending (MMB) en Reinforced DCB (R-DCB).

Mesh Convergentie: In de DCB- en MMB-benchmarks vertoonde de voorgestelde formulering een aanzienlijk verbeterde mesh-convergentie vergeleken met de conventionele formulering. Nauwkeurige last-verplaatsingscurven en tractieverdelingen werden bereikt, zelfs met grove meshes (bijv. 7,5 mm), terwijl de conventionele methode niet convergeerde of veel fijnere meshes (2 mm) vereiste voor vergelijkbare nauwkeurigheid.
Spanningvoorspelling: De voorgestelde methode verminderde de overschatting van de compressieve spanning vóór de scheurtip, die wordt waargenomen bij conventionele schel gebaseerde cohesieve elementen, aanzienlijk. De voorspelde tractieprofielen vertoonden een betere overeenstemming met solide elementmodellen en analytische oplossingen.
Mode II Prestaties: In ENF-benchmarks presteerde de voorgestelde afschuivingsstijfheid beter dan andere schel-gebaseerde formuleringen (bijv. Bazilevs et al.), die de initiële stijfheid aanzienlijk onderschatten.
3D Complexiteit (R-DCB): In de Reinforced DCB benchmark, die complexe 3D delaminatie en versterking omvat, bood de voorgestelde methode een betere overeenstemming met experimentele gegevens met betrekking tot de tweede piekbelasting en de evolutie van het scheurfront. Het legde correct het verzachtingsgedrag in versterkte regio's vast, waar de conventionele methode dit als excessief stijf voorspelde.
Computationele Efficiëntie: Door het gebruik van grovere meshes te faciliteren terwijl de nauwkeurigheid behouden blijft, vermindert de methode de benodigde rekentijd voor delaminatiesimulaties.

Betekenis
Het artikel stelt dat dit werk een kritieke beperking in de toepassing van hogere-orde structurele cohesieve elementen oplost: het onvermogen om equivalent single-layer modellen te gebruiken zonder verlies van nauwkeurigheid in de door-de-dikte spanningsvoorspelling. Door de fysieke eigenschappen van de resin-rich layer in de penalty stiffness-afleiding op te nemen via analytische spanningsverdelingen, bereikt de methode een balans tussen computationele efficiëntie en voorspellende nauwkeurigheid.

De auteurs stellen dat deze aanpak een nauwkeurige voorspelling van uit-het-vlak delaminatie en spanningsverdelingen in 3D composiet laminaatmodellen op basis van de ESL-methode mogelijk maakt. De verbeterde mesh-convergentie en de verminderde discrepantie met solide elementmodellen maken de methode bijzonder waardevol voor het simuleren van complexe structuren met versterkingen. Het artikel concludeert door te merken dat toekomstig werk zich zal richten op het uitbreiden van deze formulering om in-plane schade mechanismen (matrixcracking en vezelbreuk) te voorspellen door in-plane schelspanningen te combineren met de uit-het-vlak spanningen geleverd door de cohesieve elementen.

Het Probleem: De "Snotterige" Lijm-fout

Het Nieuwe Idee: Luisteren naar de Spanning

Het Bewijs: Het Testen van de Nieuwe Regelset

Wat Dit Betekent

Meer zoals dit