Simulation of three-dimensional tension-induced cracks based on cracking potential function-incorporated extended finite element method

Abstract

In the finite element method, the numerical simulation of three-dimensional crack propagation is relatively rare, and it is often realized by commercial programs. In addition to the geometric complexity, the determination of the cracking direction constitutes a great challenge. In most cases, the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation. However, in the case of three-dimensional analysis, the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method. In this paper, based on the extended finite element method, the stress-related function field is introduced into the calculation domain, and then the boundary value problem of the function is solved. Subsequently, the envelope surface of all propagation directions can be obtained at one time. At last, the possible surface can be selected as the direction of crack development. Based on the aforementioned procedure, such method greatly reduces the programming complexity of tracking the crack propagation. As a suitable method for simulating tension-induced failure, it can simulate multiple cracks simultaneously.

Keywords

extended finite element method / crack / three-dimensional calculation / cracking potential function / tensile failure

References

[1] TangC-a, YangY-feng. Crack branching mechanism of rock-like quasi-brittle materials under dynamic stress [J]. Journal of Central South University, 2012, 19(11): 3273-3284

[2] FormicaG, MilicchioF. Crack growth propagation using standard FEM [J]. Engineering Crack Mechanics, 2016, 165: 1-18

[3] ZhangP, DuC-b, ZhangD-H. Crack propagation modelling in concrete based on scaled boundary finite element generalized shape function [J]. Journal of Hydraulic Engineering, 2019, 50(12): 1491-1501(in Chinese)

[4] ZhaoC-b, HobbsB E, ordAConvective and Advective heat transfer in geological systems [M], 2008, Berlin, Springer

[5] ZhaoC-binDynamic and transient infinite elements: Theory and geophysical, geotechnical and Geoenvironmental Aapplications [M], 2009, Berlin, Springer

[6] ZhaoC-binPhysical and chemical dissolution front instability in porous media: Theoretical analyses and computational simulations [M], 2014, Berlin, Springer

[7] ZhaoC-b, HobbsB E, ordA. Finite element modeling of convective pore-fluid flow in fluid-saturated porous rocks within upper crust: An overview [J]. Journal of Central South University, 2019, 26: 501-514

[8] DebD, DasK C. Extended finite element method for the analysis of discontinuities in rock masses [J]. Geotechnical & Geological Engineering, 2010, 28(5): 643-659

[9] RuZ-l, ZhaoH-b, YinS-de. Evaluation of mixed-mode stress intensity factors by extended finite element method [J]. Journal of Central South University, 2013, 20(5): 1420-1425

[10] WANG Xiang-nan, YU Peng, YU Jia-lin, YU Yu-zhen, LV He. Simulated crack and slip plane propagation in soil slopes with embedded discontinuities using XFEM [J]. International Journal of Geomechanics, 2018, 18(12). DOI: https://doi.org/10.1061/(ASCE)GM.1943-5622.0001290.

[11] VIGNERON L M, VERLY J G, WARFIELD S K. On extended finite element method (XFEM) for modelling of organ deformations associated with surgical cuts [C]// Medical Simulation: International Symposium. Cambridge, MA, USA, 2004. DOI: https://doi.org/10.1007/978-3-540-25968-8_15.

[12] KrsticB, RasuoB, TrifkovicD, RadisavlievicL, RajicZ, DinulovicM. An investigation of the repetitive failure in an aircraft engine cylinder head [J]. Engineering Failure Analysis, 2013, 34(8): 335-349

[13] AreiasP M A, BelytschkoT. Analysis of three-dimensional crack initiation and propagation using the extended finite element method [J]. International Journal for Numerical Methods in Engineering, 2005, 63(5): 760-788

[14] LoehnertS, Mueller-hoeppeD S, wriggersP. 3D corrected XFEM approach and extension to finite deformation theory [J]. International Journal for Numerical Methods in Engineering, 2011, 86(4): 431-452 5

[15] WangZ, YuT-t, TinhQ B, satoyukiT, ZhangC-z, SohichiH, JoseL, CurielS. 3-D Local mesh refinement XFEM with variable-node hexahedron elements for extraction of stress intensity factors of straight and curved planar cracks [J]. Computer Methods in Applied Mechanics and Engineering, 2017, 313: 375-405

[16] AgathosK, ChatziE, BordasS P A. Multiple crack detection in 3D using a stable XFEM and global optimization [J]. Computational Mechanics, 2018, 62: 835-852

[17] ZhaoC-b, HobbsB E, ordAFundamentals of computational geoscience: Numerical methods and algorithms [M], 2009, Berlin, Springer

[18] OLIVER X. On strategies for tracking strong discontinuities in computational failure mechanics [C]// Fifth World Conference on Computational Mechanics. Vienna, Austria, 2002. https://previa.uclm.es/profesorado/evieira/ftp/articulos/congreso/wccm2_2002.pdf.

[19] CHAVES E W V, OLIVER X. A three-dimensional setting for strong discontinuities modeling in failure mechanics [M]. Monograph CIMNE N-73, 2003. https://www.tdx.cat/handle/10803/6861

[20] PhilippeJ, PaulS, EllenK. Towards the treatment of boundary conditions for global crack path tracking in three-dimensional brittle fracture [J]. Computational Mechanics, 2009, 45(1): 91-107

[21] PeterD, GüntherM. Crack propagation criteria in the framework of X-FEM-based structural analyses [J]. International Journal for Numerical & Analytical Methods in Geomechanics, 2010, 31(2): 239-259

[22] ManzoliO L, claroG K S, RodriguesE A, lopesJ. A local-global scheme for tracking crack path in three-dimensional solids [J]. Computers & Concrete, 2013, 12(3): 261-283

[23] PengC, WuW, ZhangB-Y. Three-dimensional simulations of tensile cracks in geomaterials by coupling meshless and finite element method [J]. International Journal for Numerical & Analytical Methods in Geomechanics, 2015, 39(2): 135-154

[24] SongJ H, areiasP M A, BelytschkoT. A method for dynamic crack and shear band propagation with phantom nodes [J]. International Journal for Numerical Methods in Engineering, 2006, 67(6): 868-893

[25] HillerborgA, ModeerM, PetersonP E. Analysis of crack propagation and crack growth in concrete by means of crack mechanics and finite elements [J]. Cement Concrete Research, 1976, 6: 773-782

[26] PlanasJ, ElicesM, GunineaG V, gomezF J, cendonD A, arbillaI. Generalizations and specializations of cohesive crack models [J]. Engineering Crack Mechanics, 2003, 70(14): 1759-1776

[27] FangX-j, JinF, WangJ-ting. Cohesive crack model based on finite element method [J]. Journal of Tsinghua University (Natural Science Edition), 2007, 47(3): 344-347(in Chinese)

[28] WellsG N, sluysL J. A new method for modelling cohesive cracks using finite elements [J]. International Journal for Numerical Methods in Engineering, 2001, 50: 2667-2682

[29] ArreaM, IngraffeaA RMixed-mode crack propagation in mortar and concrete [R], 1982, Ithaca, Department of Structure Engineering, Cornell University

[30] BoccaP, CappinteriA, ValenteS. Mixed mode fracture of concrete [J]. International Journal of Solids and Structures, 1991, 27(9): 1139-1153

[31] ZhangB-y, ZhangM-c, SunX. Centrifugal modeling of transverse cracking in earth core dams [J]. Geomechanics, 2008, 29(5): 1254-1258(in Chinese)

[32] ZhangY, WangJ-g, ZhangB-y, LiQ. Meshless method for numerical simulation of crack propagation in earth dams [J]. Journal of Geotechnical Engineering, 2009, 31(5): 727-731(in Chinese)