Door discrete scheuren in beton numeriek te simuleren, kunnen verschillende aspecten van het constructief gedrag van betonelementen worden bestudeerd. In een afstudeeronderzoek aan de TU Eindhoven is aan de hand van modellen, de invloed onderzocht van discrete scheuren op de onthechting van uitwendig opgelijmde wapening. Stap voor stap zijn verschillende materiaalmodellen onderzocht en op geschiktheid vergeleken. Auteur:ir. Steven Schoenmakers (Van Rossum Raadgevende Ingenieurs b.v.) ENCI Studieprijs 2014Dit is het vierde artikel in een serie met bijdragen van prijswinnaars van de ENCI Studieprijs 2014. De studie die in dit artikel wordt beschreven, ontving een eervolle vermelding in de categorie Universiteiten. De jury over deze studie: "[?] De jury is onder de indruk van het vermogen van de onderzoeker een eigen materiaalmodel (UMAT) te ontwikkelen, waarmee de prestatie van de beschikbare modellen is verbeterd. Het werk wordt gekenmerkt, door een systematische aanpak van de problematiek, met een kritische analyse van de verkregen resultaten en een duidelijke verantwoording van de uiteindelijk gekozen modellering. Het wetenschappelijk niveau is hoog en de onderzoeker toont een diepgaand inzicht in complexe materiaalmodellen. De presentatie is verzorgd en het werk is ondanks de diepgang vlot leesbaar." Lees ook het volledige afstudeerrapport van Steven Schoenmakers.
Modelleren discrete scheurvorming5 20152ModellerendiscretescheurvormingInvloed van discrete scheuren op deonthechting van uitwendig opgelijmdewapening onderzochtModelleren discrete scheurvorming 5 2015 3effectieve relatieve verplaatsingeffectievetractiemax0tmaxKK (1 - d)Gc4321In Nederland staan momenteel veel betonnen gebouwen leeg.Deze gebouwen kunnen mogelijk worden herbestemd. Vooreen deel van deze gebouwen is versterking van de betoncon-structie benodigd voordat het gebouw kan worden hergebruikt.Een beproefde methode voor die versterking is het aanbrengenvan extern gelijmde wapening, zoals bijvoorbeeld koolstof-vezelwapening (fig. 1). Bezwijken van een constructie die isversterkt met koolstofvezelwapening, vindt meestal plaatsdoordat de extern gelijmde wapening onthecht van het beton.Deze onthechting treedt op ter plaatse van de scheurvormingin het beton en resulteert in een bros bezwijkmechanisme.Indien men de onthechting van extern gelijmde wapening wilvoorkomen, kan het nodig zijn het onthechtgedrag te simulerenmet een analyse met behulp van een eindige-elementenmodel.Het is hiervoor noodzakelijk om het discrete scheurpatroonvan het betonelement te kunnen simuleren.OnderzoekIn het afstudeeronderzoek is het hechtgedrag tussen twee mate-rialen bestudeerd. Daarnaast zijn verschillende mogelijkhedenvoor het modelleren van discrete scheurvorming beschouwd.Zowel het hechtgedrag als het modelleren van discrete scheur-vorming is beschouwd voor 2D-configuraties.Scheurvorming kan worden beschouwd als een schade dielokaal optreedt in een materiaal, dat kan worden beschrevenmet materiaalmodellen in een numeriek model. Hiertoe kaneen tractie-separatiediagram worden gedefinieerd als materi-aalmodel, zoals weergegeven in figuur 2. Hierin kunnen vierpunten worden onderscheiden:1.initi?le lineair elastische tak;2.schade initiatie criterium;3.schade verloop;4.schade is volledig.ir. Steven Schoenmakers 1)Van Rossum RaadgevendeIngenieurs b.v.1 Vloer versterkt met extern gelijmde wapening2 Schadeverloop: ti= (1 ? d) Kijj1) Steven Schoenmakers is met het onderzoek Hetmodel≤renvandiscretescheurvorming∈be→n′af≥studeerdaandeTUE∈doven,fact̲eitBouwkunde,StructuralDesign.Hijwerddaarbijbe≥≤idd∞rprof.dr.ir.A.S.J.Suiker,prof.dr.ir.T.A.M.Sa≤tendr.ir.C.Veρosel(al≤nTUe).Eenl∈knaarhetrapportstaatopwww.cemen→nl∈e.nl.ENCIStudieprijs2014Ditishetvierdeartikel∈eenseriemetbijdra≥nvanprijsw∈naarsvandeENCIStudieprijs2014.Destudiedie∈ditartikelwordtbeschreven,ontv∈gee≠ervol≤vermeld∈g∈decategorieUniversiteiten.Dejuryoverdezestudie:[...] De jury is onder de indruk van het vermogenvan de onderzoeker een eigen materiaalmodel (UMAT) teontwikkelen, waarmee de prestatie van de beschikbare modellenis verbeterd. Het werk wordt gekenmerkt door een systematischeaanpak van de problematiek, met een kritische analyse van deverkregen resultaten en een duidelijke verantwoording van deuiteindelijk gekozen modellering. Het wetenschappelijk niveau ishoog en de onderzoeker toont een diepgaand inzicht incomplexe materiaalmodellen. De presentatie is verzorgd en hetwerk is ondanks de diepgang vlot leesbaar.Meer∈formatieopwww.cemen→nl∈e.nlencistudieprijs.D∞rdiscretescheuren∈be→nνmeriektesiμ≤ren,ku∩enversχl≤ndeaspectenvanhetconstructief≥dragvanbe→≠≤mentenwordenbestudeerd.IneenafstudeeronderzoekaandeTUE∈dhovenisaandehandvanmodel≤n,de∈vloedonderzochtvandiscretescheurenopdeonthecht∈gvanuitwendigop≥lijmdewapen∈g.Stapv∞rstapzijnversχl≤ndemateriaalmodel≤nonderzochtenop≥sχktheidver≥≤ken.12Mo∂≤rendiscretescheurvorming5201543H∞fdricht∈≥nbel∗∈≥n4Knopenverreiktmeteenaanv≤̲ndeverplaatsingunctiebijdeXFEM-methode(∈groenweer≥≥ven)een≥r∈≥mesh-afhankelijkheidheeft.Ditkomtd∞rdatnietdevol≤digaanv≤̲ndeverplaatsingunctiewordt→e≥voegdaandee≤menten∈hetscheurdome∈,waard∞rhetscheurontal≤enopeene≤mentrandkanoptreden.CohesiveZo≠Mo∂(CZM)Eencohesivezo≠modeliseenmodelwaarbijtussenal≤cont∈u?me≤menten,∫erface-e≤mentenzijnaan≥bracht(fig.5).Deze∫erface-e≤mentenzijn∈staatschadetebeschrijven,zodatdiscretescheurvorming∈hetmodelkanoptredenterplaatsevande∫erface-e≤menten.Erkanworden≥kozende∫erface-e≤mentenmetofzonderfysiekediktetedef∈i?ren.De∫erface-e≤mentenwordend∞rmid∂vaneenPython-script≥def∈ieerdenvervol≥ns≥?mporteerd∈depreprocessor.BenchmarkstudieIneenbenchmarkstudiezijndeXFEM-methodeenCZM-methodever≥≤ken.Indezestudieishetontstaanvandiscretescheurvorming≥comb∈eerdmetdeonthecht∈gtussentweemateria≤n.Hierbijwordtv∞rhetmodel≤renvandehecht∈>ussentweemateria≤n≥bruik≥maaktvan∫erface-e≤menten,zoalsdeze∞k∈hetcohesivezo≠modelworden→e≥p∗v∞rhetbeschrijvenvandehecht∈g.Bijdestudieiseenvezel-ep⊗ysysteembel∗optrekbeschouwd,zoalsweer≥≥ven∈fig∪r6.Deep⊗ywordtl∞drechtopdericht∈gvandevezelbel∗optrek,d∞rmid∂vaneenop≥≤gdeverplaatsing.Als≥volghiervanontstaatonthecht∈gtussendevezelendeep⊗y,evenalsscheurvorming∈deep⊗ymatrix.D∞rAlfaroetal.[1]iseerdereenstudieuit≥voerdnaarditsysteem,waarbijzowelnumeriekonderzoekd∞rmid∂vaneen∫erfacedama≥model,alsexperimenteelonderzoekisuit≥voerd.Ditonder-zoekisalsbenchmark≥bruiktv∞rdetweeonderzochtemethodes.Infig∪r7zijndeoptredendescheurpatro≠nweer-≥≥ven.Fig∪r8laathet≥mid∂despa∩∈gsverl∞pzienaanderandvanhetmodelwaardeverplaatsingwordtop≥≤gd.Hettractie-separatiediagramdientv∞reen3D-modelteworden≥def∈ieerdv∞rdedrieh∞fdricht∈≥n,namelijkde∥a∥al-richt∈gendetweeafschuifricht∈≥n.V∞reen2D-modelzijndits≤chtstweericht∈≥n,modeIenmodeII(fig.3).Mo∂≤r∈ghecht∈gInhetonderzoekisaan≥→onddathetmodel≤renvandehechtei≥nschappend∞rmid∂vaneenoppervlakte-ei≥n-schaps≤chtseenbeperkt→epassings≥biedheeft.OndermodeIIbel∗∈gblijkteene≤ment,bijoptredendeverplaatsin≥ngroterdande∈iti?≤e≤ment≤n>e,niet∈staatspa∩∈≥novertedra≥n.Eenmeernauwkeuri≥beschrijv∈gvanhechtei≥nschappenkanwordenbereiktd∞rhet→epassenvan∫erface-e≤menten.Bijhet→epassenhiervanwordtzowelv∞rmodeI,modeII,alsmixed-modebel∗∈geen≥dragovereenkomendaande≥model≤erdewaarden≥vonden.Mo∂≤r∈gscheurvormingTweemethodenv∞rhetbeschrijvenvandiscretescheurvor-mingku∩enwordenonderscheiden:deextendedf∈itee≤mentmethodenhetcohesivezo≠model.eXtendedF∈iteE≤mentMethod(XFEM)Bijdeextendedf∈itee≤mentmethodwordteenaanv≤̲ndeverplaatsingunctie→e≥voegdaandee≤menten∈hetscheur-dome∈.Dezeverplaatsingunctiemaakthetmo≥lijkdisconti-ν?teiten∈desterkteb∈≠neene≤menttebeschrijven.Dezediscont∈u?teitenontstaan∈debereken∈gopdezwakste' plaatsin de doorsnede. Op het moment dat `schade' aan elementenontstaat wordt aan deze elementen een aanvullende verplaat-singsfunctie toegekend (fig. 4). Deze aanvullende verplaatsings-functie wordt vooraf gedefinieerd, maar in de berekening uitslui-tend toegepast op de elementen waarvoor dit van toepassing is.In theorie is de XFEM-methode mesh-onafhankelijk. Het blijktechter dat de implementatie van deze methode in het toege-paste eindige-elementensoftwarepakket (Abaqus 6.11.2) toch43mode I (N) mode II (S) mode III (T)NI 1 i i i( ) ( )u N x u H x a= = + Modelleren discrete scheurvorming 5 2015 55 Plaatsing interface-elementen (in blauw weergegeven) in een cohesive zone model6 Vezel-epoxysysteem7 Optredende scheurpatronen (a) onderzoek Alfaro et al. [1]; (b) XFEM-methode;(c) cohesive zone modelelementen worden toegepast met een geringe fysieke dikte,doordat deze normaal van de interface-elementen eenduidig isgedefinieerd in het model. Er wordt ook altijd aangeraden eencohesive zone model toe te passen waarbij de interface-elemen-ten een fysieke dikte hebben.MateriaalmodelVoor het beschrijven van mixed-modegedrag is een gebruikers-materiaalmodel (UMAT) ge?mplementeerd. Dit UMAT is zelfgeprogrammeerd voor dit onderzoek, met als doel de viscositeiten de wijze waarop deze in de berekening is verwerkt eenduidigte beschrijven. Het toepassen van deze viscositeit (in eenbeperkte mate) is aanbevelingswaardig zodat het convergentie-gedrag van het model kan worden gewaarborgd.Het gebruikersmateriaalmodel kan worden toegekend aanelementen in het model. Het materiaalmodel in dit onderzoekModellering scheurpatroonZowel met de XFEM-methode als met het cohesive zone modelkan een discreet scheurpatroon beschreven worden. , de hoekwaaronder scheurvorming in het epoxy materiaal optreedt(fig. 7), kan analytisch worden afgeleid en bedraagt 68?. Hettheoretische bepaalde scheurpatroon wordt echter beter bena-derd door het toepassen van een cohesive zone model (fig. 7).Ook het optredende spanningsverloop in het model wordtbeter benaderd met het cohesive zone model ten opzichte vande XFEM-methode (fig. 8). Het optredende scheurpatroon bijtoepassing van een cohesive zone model is, in tegenstelling totde XFEM-methode, echter wel mesh-afhankelijk, aangezienhier scheurvorming alleen kan optreden tussen de continu?m-elementen. Voor het bepalen van een realistisch scheurpatroonis het van belang een voldoende willekeurige en fijne mesh tegenereren, zodat er geen voorkeursrichting aanwezig is in demesh.Het toepassen van de interface-elementen in een cohesive zonemodel met of zonder fysieke dikte heeft geen invloed op hetoptredende scheurpatroon en spanningsverloop. De benodigderekentijd blijkt echter aanzienlijk korter wanneer de interface-567a 7b 7cepoxy = 40 m14m9 m5 m14 m5m9mvezel = 56? = 36?Modelleren discrete scheurvorming5 201568 Optredene spanningen9 Experimentele en numerieke (in rood weergegeven)resultaten scheurpatroon betonliggersis geprogrammeerd met een bilineair materiaalgedag: een initi?lelineair-elastische tak, waarna lineaire softening kan optredenzoals het schade verloop (fig. 2) is weergegeven. Het UMATlevert voor de verschillende belastingtypen resultaten dienauwkeurig overeenkomen met de gevonden resultaten in eenbeschouwd referentieonderzoek.8910,90,80,70,60,50,40,30,20,1021,81,61,41,20,80,60,40,20 1XFEMCZMAlfaro et al.tucohucohVergelijk met experimenteel onderzoekHet cohesive zone model is toegepast voor het modelleren vanhet gedrag van een betonbalk waarbij de mate van slankheid isgevarieerd. Deze betonbalk is eerder experimenteel bepaald enbeschreven door Leonhardt en Walther [2]. Voor de verschil-lende slankheden van de ligger komt het gesimuleerde scheur-patroon nauwkeurig overeen met de experimentele resultaten.In figuur 9 is in het rood het numeriek bepaalde scheurpatroonweergegeven voor de experimentele resultaten van Leonhardten Walter [2], die in zwart zijn weergegeven. Voor de semi-gedrongen en slanke liggers komt de bezwijkcapaciteit van degesimuleerde liggers goed overeen met de experimentele resul-taten. Voor de meer gedrongen liggers (a/d 2) wordt eenafwijking in bezwijkcapaciteit gevonden.Conclusie en aanbevelingenOp basis hiervan kan worden geconcludeerd dat het cohesivezone model kan worden gebruikt voor het modelleren vandiscrete scheurvorming bij semi-gedrongen en slank gewa-pende betonnen liggers.Aanbevolen wordt in vervolgonderzoek de cohesive-zone-modelbenadering uit te breiden voor 3D-configuraties.Zodoende kunnen meer complexe problemen gesimuleerdworden zoals betonliggers met dwarskrachtwapening. Ook kannaar verwachting, de bezwijkcapaciteit van gedrongen beton-liggers (a/d 2) nauwkeuriger benaderd worden met3D-modellen.Als laatste dient in toekomstig onderzoek de invloed van lokaleimperfecties in de toegepaste materialen te worden opgenomenin de numerieke simulaties. De eigenschappen van de lokaleimperfecties, zoals grindkorrels in beton, kunnen wordenverwerkt in de numerieke simulaties door een Monte-Carlo-verdeling toe te passen op de materiaaleigenschappen in hetmodel. LITERATUUR1 Alfaro, M. V., Suiker, A. S. J. en Borst, R.,Transverse Failure Behavior of Fibre-epoxySystems. Journal of Composite Materials,Vol. 44, No. 12/2010, 2010.2 Leonhardt, F. en Walther, R., Schubversuchean einfeldrigen Stahlbetonbalken mit undohne Schubbewehrung. DeutscherAusschuss f?r Stahlbeton, nr. 151, 1962.
Reacties