An investigation on prevalent strategies for XFEM-based numerical modeling of crack growth in porous media

Abstract

Crack growth modeling has always been one of the major challenges in fracture mechanics. Among all numerical methods, the extended finite element method (XFEM) has recently attracted much attention due to its ability to estimate the discontinuous deformation field. However, XFEM modeling does not directly lead to reliable results, and choosing a strategy of implementation is inevitable, especially in porous media. In this study, two prevalent XFEM strategies are evaluated: a) applying reduced Young’s modulus to pores and b) using different partitions to the model and enriching each part individually. We mention the advantages and limitations of each strategy via both analytical and experimental validations. Finally, the crack growth is modeled in a natural porous media (Fontainebleau sandstone). Our investigations proved that although both strategies can identically predict the stress distribution in the sample, the first strategy simulates only the initial crack propagation, while the second strategy could model multiple cracks growths. Both strategies are reliable and highly accurate in calculating the stress intensity factor, but the second strategy can compute a more reliable reaction force. Experimental tests showed that the second strategy is a more accurate strategy in predicting the preferred crack growth path and determining the maximum strength of the sample.

Graphical abstract

Keywords

numerical modeling / extended finite element method / porous media / crack growth / stress intensity factor

References

[1] ChangS H, LeeC I, JeonS. Measurement of rock fracture toughness under modes I and II and mixed-mode conditions by using disc-type specimens. Engineering Geology, 2002, 66( 1−2): 79– 97

[2] HoekE, MartinC D. Fracture initiation and propagation in intact rock—A review. Journal of Rock Mechanics and Geotechnical Engineering, 2014, 6( 4): 287– 300

[3] LisjakA, KaifoshP, HeL, TatoneB S A, MahabadiO K, GrasselliG. A 2D, fully-coupled, hydro-mechanical, FDEM formulation for modelling fracturing processes in discontinuous, porous rock masses. Computers and Geotechnics, 2017, 81 : 1– 18

[4] JingL, HudsonJ A. Numerical methods in rock mechanics. International Journal of Rock Mechanics and Mining Sciences, 2002, 39( 4): 409– 427

[5] Cundall P A. A computer model for simulating progressive large scale movements in blocky rock systems. In: Proceedings of the International Symposium on Rock Mechanics. Nancy: International Society for Rock Mechanics, 1971

[6] LisjakA, GrasselliG. A review of discrete modeling techniques for fracturing processes in discontinuous rock masses. Journal of Rock Mechanics and Geotechnical Engineering, 2014, 6( 4): 301– 314

[7] WangS Y, SloanS W, ShengD C, YangS Q, TangC A. Numerical study of failure behaviour of pre-cracked rock specimens under conventional triaxial compression. International Journal of Solids and Structures, 2014, 51( 5): 1132– 1148

[8] KatoT, NishiokaT. Analysis of micro–macro material properties and mechanical effects of damaged material containing periodically distributed elliptical microcracks. International Journal of Fracture, 2005, 131( 3): 247– 266

[9] RezanezhadM, LajevardiS A, KarimpouliS. Effects of pore(s)-crack locations and arrangements on crack growth modeling in porous media. Theoretical and Applied Fracture Mechanics, 2020, 107 : 102529–

[10] Rodriguez-Florez N. Mechanics of cortical bone: Exploring the micro- and nano-scale. Dissertation for the Doctoral Degree. London: Imperial College London, 2015

[11] DuarteA P C, SilvaB A, SilvestreN, deBrito J, JúlioE. Mechanical characterization of rubberized concrete using an Image-Processing/XFEM coupled procedure. Composites Part B: Engineering, 2015, 78 : 214– 226

[12] DuarteA P C, SilvestreN, deBrito J, JúlioE. Numerical study of the compressive mechanical behaviour of rubberized concrete using the extended finite element method (XFEM). Composite Structures, 2017, 179 : 132– 145

[13] Supar K, Ahmad H. XFEM Modelling of Multi-holes Plate with Single-row and Staggered Holes Configurations. In: International Symposium on Civil and Environmental Engineering 2016 (ISCEE 2016). Wuhan: MATEC Web of Conferences, 2017

[14] RezanezhadM, LajevardiS A, KarimpouliS. Application of equivalent circle and ellipse for pore shape modeling in crack growth problem: A numerical investigation in microscale. Engineering Fracture Mechanics, 2021, 253 : 107882–

[15] HedjaziL, GuessasmaS, Della ValleG, BenseddiqN. How cracks propagate in a vitreous dense biopolymer material. Engineering Fracture Mechanics, 2011, 78( 6): 1328– 1340

[16] HedjaziL, MartinC L, GuessasmaS, Della ValleG, DendievelR. Application of the Discrete Element Method to crack propagation and crack branching in a vitreous dense biopolymer material. International Journal of Solids and Structures, 2012, 49( 13): 1893– 1899

[17] ChenM, WangH, JinH, PanX, JinZ. Effect of pores on crack propagation behavior for porous Si3N4 ceramics. Ceramics International, 2016, 42( 5): 5642– 5649

[18] Rodriguez-FlorezN, CarrieroA, ShefelbineS J. The use of XFEM to assess the influence of intra-cortical porosity on crack propagation. Computer Methods in Biomechanics and Biomedical Engineering, 2017, 20( 4): 385– 392

[19] RezanezhadM, LajevardiS A, KarimpouliS. Effects of pore-crack relative location on crack propagation in porous media using XFEM method. Theoretical and Applied Fracture Mechanics, 2019, 103 : 102241–

[20] RabczukT, BelytschkoT. Cracking particles: A simplified meshfree method for arbitrary evolving cracks. International Journal for Numerical Methods in Engineering, 2004, 61( 13): 2316– 2343

[21] RabczukT, BelytschkoT. A three-dimensional large deformation meshfree method for arbitrary evolving cracks. Computer Methods in Applied Mechanics and Engineering, 2007, 196( 29−30): 2777– 2799

[22] RabczukT, ZiG, BordasS, Nguyen-XuanH. A simple and robust three-dimensional cracking-particle method without enrichment. Computer Methods in Applied Mechanics and Engineering, 2010, 199( 37−40): 2437– 2455

[23] RenH, ZhuangX, RabczukT. Dual-horizon peridynamics: A stable solution to varying horizons. Computer Methods in Applied Mechanics and Engineering, 2017, 318 : 762– 782

[24] RenH L, ZhuangX Y, AnitescuC, RabczukT. An explicit phase field method for brittle dynamic fracture. Computers & Structures, 2019, 217 : 45– 56

[25] RenH, ZhuangX, RabczukT. A higher order nonlocal operator method for solving partial differential equations. Computer Methods in Applied Mechanics and Engineering, 2020, 367 : 113132–

[26] AreiasP, MsekhM A, RabczukT. Damage and fracture algorithm using the screened Poisson equation and local remeshing. Engineering Fracture Mechanics, 2016, 158 : 116– 143

[27] LehoucqR B, SillingS A. Force flux and the peridynamic stress tensor. Journal of the Mechanics and Physics of Solids, 2008, 56( 4): 1566– 1577

[28] KarimpouliS, TahmasebiP. A hierarchical sampling for capturing permeability trend in rock physics. Transport in Porous Media, 2017, 116( 3): 1057– 1072

[29] KarimpouliS, TahmasebiP, SaengerE H. Estimating 3D elastic moduli of rock from 2D thin-section images using differential effective medium theory. Geophysics, 2018, 83( 4): MR211– MR219

[30] HillerborgA, ModéerM, PeterssonP E. Analysis of crack formation and crack growth in concrete by means of fracture mechanics and finite elements. Cement and Concrete Research, 1976, 6( 6): 773– 781

[31] AsadpoureA, MohammadiS, VafaiA. Crack analysis in orthotropic media using the extended finite element method. Thin-walled Structures, 2006, 44( 9): 1031– 1038

[32] Mohammadi S. Extended Finite Element Method: For Fracture Analysis of Structures, Oxford: Blackwell Publishing Ltd, 2008

[33] SharafisafaM, NazemM. Application of the distinct element method and the extended finite element method in modelling cracks and coalescence in brittle materials. Computational Materials Science, 2014, 91 : 102– 121

[34] MoësN, BelytschkoT. Extended finite element method for cohesive crack growth. Engineering fracture mechanics, 2002, 69( 7): 813– 833

[35] GinerE, SukumarN, TarancónJ E, FuenmayorF J. An Abaqus implementation of the extended finite element method. Engineering Fracture Mechanics, 2009, 76( 3): 347– 368

[36] BelytschkoT, BlackT. Elastic crack growth in finite elements with minimal remeshing. International Journal for Numerical Methods in Engineering, 1999, 45( 5): 601– 620

[37] DolbowJ, MoësN, BelytschkoT. An extended finite element method for modeling crack growth with frictional contact. Computer Methods in Applied Mechanics and Engineering, 2001, 190( 51−52): 6825– 6846

[38] LiL, WangM Y, WeiP. XFEM schemes for level set based structural optimization. Frontiers of Mechanical Engineering, 2012, 7( 4): 335– 356

[39] MoësN, DolbowJ, BelytschkoT. A finite element method for crack growth without remeshing. International Journal for Numerical Methods in Engineering, 1999, 46( 1): 131– 150

[40] MoësN, CloirecM, CartraudP, RemacleJ F. A computational approach to handle complex microstructure geometries. Computer Methods in Applied Mechanics and Engineering, 2003, 192( 28−30): 3163– 3177

[41] ZhuQ Z. On enrichment functions in the extended finite element method. International Journal for Numerical Methods in Engineering, 2012, 91( 2): 186– 217

[42] AgathosK, ChatziE, BordasS P A. Multiple crack detection in 3D using a stable XFEM and global optimization. Computational Mechanics, 2018, 62( 4): 835– 852

[43] Sih G C. Methods of Analysis and Solution of Crack Problems. Leyden: Noordhoff International Publishing, 1973

[44] Anderson T L. Fracture Mechanics: Fundamentals and Applications. 3rd ed. Boca Raton: Taylor and Francis, 2005

[45] ArshadnejadS. Analysis of the first cracks generating between two holes under incremental static loading with an innovation method by numerical modelling. Mathematics and Computer Science, 2017, 2( 6): 120– 129

[46] ZhangZ. An empirical relation between mode I fracture toughness and the tensile strength of rock. International Journal of Rock Mechanics and Mining Sciences, 2002, 39( 3): 401– 406

[47] BažantZ P, KazemiM T. Size effect in fracture of ceramics and its use to determine fracture energy and effective process zone length. Journal of the American Ceramic Society, 1990, 73( 7): 1841– 1853

[48] MarshallG P, WilliamsJ G, TurnerC E. Fracture toughness and absorbed energy measurements in impact tests on brittle materials G. Journal of Materials Science, 1973, 8( 7): 949– 956

[49] Nasaj MoghaddamH, KeyhaniA, AghayanI. Modelling of crack propagation in layered structures using extended finite element method. Civil Engineering Journal, 2016, 2( 5): 180– 188

[50] ZhangC, CaoP, CaoY, LiJ. Using finite element software to simulation fracture behavior of three-point bending beam with initial crack. Journal of Software, 2013, 8( 5): 1145– 1150

[51] AbdellahM Y. Delamination modeling of double cantilever beam of unidirectional composite laminates. Journal of Failure Analysis and Prevention, 2017, 17( 5): 1011– 1018

[52] GrigoriuM, SaifM T A, El BorgiS, IngraffeaA R. Mixed mode fracture initiation and trajectory prediction under random stresses. International Journal of Fracture, 1990, 45( 1): 19– 34

[53] S. Moaveni, Finite Element Analysis: Theory and Application with ANSYS. Hoboken: Prentice Hall, 1999

[54] TroyaniN, PérezA, BaízP. Effect of finite element mesh orientation on solution accuracy for torsional problems. Finite Elements in Analysis and Design, 2005, 41( 14): 1377– 1383

[55] Logan D L. A First Course in the Finite Element Method. 4th ed. Toronto: Nelson, 2007

[56] SongJ, BelytschkoT. Cracking node method for dynamic fracture with finite elements. International Journal for Numerical Methods in Engineering, 2009, 77( 3): 360– 385

[57] LinderC, ArmeroF. Finite elements with embedded branching. Finite Elements in Analysis and Design, 2009, 45( 4): 280– 293

[58] LiX, KonietzkyH. Simulation of time-dependent crack growth in brittle rocks under constant loading conditions. Engineering Fracture Mechanics, 2014, 119 : 53– 65

[59] AndräH, CombaretN, DvorkinJ, GlattE, HanJ, KabelM, KeehmY, KrzikallaF, LeeM, MadonnaC, MarshM, MukerjiT, SaengerE H, SainR, SaxenaN, RickerS, WiegmannA, ZhanX. Digital rock physics benchmarks—Part I: Imaging and segmentation. Computers & Geosciences, 2013, 50 : 25– 32

[60] MadonnaC, QuintalB, FrehnerM, AlmqvistB S G, TisatoN, PistoneM, MaroneF, SaengerE H. Synchrotron-based X-ray tomographic microscopy for rock physics investigations. Geophysics, 2013, 78( 1): D53– D64

[61] Huang J Q, Huang Q A, Qin M, Dong W J, Chen X W. Experimental study on the dielectrostriction of SiO₂ with a micro-fabricated cantilever. In: IEEE Sensors 2009 Conference. Christchurch: IEEE, 2009

[62] KarimpouliS, TahmasebiP, SaengerE H. Coal cleat/fracture segmentation using convolutional neural networks. Natural Resources Research, 2020, 29( 3): 1675– 1685

RIGHTS & PERMISSIONS

2021 Higher Education Press 2021.