An eXtended Finite Element Method Implementation in COMSOL Multiphysics: Solid Mechanics

This paper presents the first time implementation of the eXtended Finite Element Method (XFEM) in the general purpose commercial software COMSOL Multiphysics. An enrichment strategy is proposed, consistent with the structure of the software. To this end, for each set of enrichment functions, an additional Solid Mechanics module is incorporated into the numerical framework, coupled with compatible modifications to the internal variables. The Linear Elastic Fracture Mechanics (LEFM) is exclusively adopted for the crack analysis. The model pre-processing, level set update, stress intensity factor calculation and crack propagation analysis are conducted by employing COMSOL’s built-in features in conjunction with external MATLAB functions through COMSOL LiveLink. All implementational aspects and suggested remedies for the treatment of enriched elements, framework setup, evaluation of stress intensity factors, and numerical integration are described in detail. The accuracy and robustness of the proposed method are examined by several numerical examples for stationary and propagating crack problems in 2D and 3D settings. The results represent excellent agreement with available analytical, numerical and experimental observations in the literature.

[1]  Amir R. Khoei,et al.  Modeling of crack propagation via an automatic adaptive mesh refinement based on modified superconvergent patch recovery technique , 2008 .

[2]  A. Khoei.  Extended Finite Element Method: Theory and Applications , 2015 .

[3]  Kamy Sepehrnoori,et al.  XFEM-Based CZM for the Simulation of 3D Multiple-Cluster Hydraulic Fracturing in Quasi-Brittle Shale Formations , 2015, Rock Mechanics and Rock Engineering.

[4]  Joseph M. Wallace,et al.  Mechanics of linear microcracking in trabecular bone. , 2019, Journal of biomechanics.

[5]  Timon Rabczuk,et al.  Phase field modeling of quasi-static and dynamic crack propagation: COMSOL implementation and case studies , 2018, Adv. Eng. Softw..

[6]  Amir R. Khoei,et al.  An augmented Lagrangian contact formulation for frictional discontinuities with the extended finite element method , 2015 .

[7]  Guirong Liu,et al.  Smoothed Point Interpolation Methods: G Space Theory and Weakened Weak Forms , 2013 .

[8]  Amir R. Khoei,et al.  The extended finite element method for large deformation ductile fracture problems with a non-local damage-plasticity model , 2013 .

[9]  I. Babuska,et al.  The design and analysis of the Generalized Finite Element Method , 2000 .

[10]  T. Anderson,et al.  Fracture mechanics - Fundamentals and applications , 2017 .

[11]  Guirong Liu,et al.  Extended finite element method with edge-based strain smoothing (ESm-XFEM) for linear elastic crack growth , 2012 .

[12]  S. Mohammadi,et al.  An adapting cohesive approach for crack-healing analysis in SMA fiber-reinforced composites , 2019, Computer Methods in Applied Mechanics and Engineering.

[13]  Nasser Khalili,et al.  Fully Coupled XFEM Model for Flow and Deformation in Fractured Porous Media with Explicit Fracture Flow , 2016 .

[14]  T. Belytschko,et al.  An Eulerian–Lagrangian method for fluid–structure interaction based on level sets , 2006 .

[15]  Amir R. Khoei,et al.  X-FEM modeling of large plasticity deformation; a convergence study on various blending strategies for weak discontinuities , 2015 .

[16]  A. Khoei,et al.  An extended finite element method for hydraulic fracture propagation in deformable porous media with the cohesive crack model , 2013 .

[17]  R. Gracie,et al.  Simulation of fracture propagation induced acoustic emission in porous media , 2020 .

[18]  T. Belytschko,et al.  New crack‐tip elements for XFEM and applications to cohesive cracks , 2003 .

[19]  K. Terada,et al.  Simulations of cohesive fracture behavior of reinforced concrete by a fracture-mechanics-based damage model , 2019, Engineering Fracture Mechanics.

[20]  R. Gracie,et al.  An enriched finite element model for wave propagation in fractured media , 2017 .

[21]  Zhigang Suo,et al.  Partition of unity enrichment for bimaterial interface cracks , 2004 .

[22]  Ted Belytschko,et al.  Elastic crack growth in finite elements with minimal remeshing , 1999 .

[23]  I. Babuska,et al.  The Partition of Unity Method , 1997 .

[24]  Soheil Mohammadi,et al.  Extended Finite Element Method: for Fracture Analysis of Structures , 2008 .

[25]  J. Prévost,et al.  Modeling quasi-static crack growth with the extended finite element method Part I: Computer implementation , 2003 .

[26]  T. Belytschko,et al.  Arbitrary branched and intersecting cracks with the eXtended Finite Element Method , 2000 .

[27]  I. Babuska,et al.  The partition of unity finite element method: Basic theory and applications , 1996 .

[28]  Ted Belytschko,et al.  Concurrently coupled atomistic and XFEM models for dislocations and cracks , 2009 .

[29]  S. Sundararajan,et al.  Investigating the effect of retained austenite and residual stress on rolling contact fatigue of carburized steel with XFEM and experimental approaches , 2018, Materials Science and Engineering: A.

[30]  A. Khoei,et al.  An enriched–FEM technique for numerical simulation of interacting discontinuities in naturally fractured porous media , 2018 .

[31]  Francisco Cruz,et al.  An XFEM implementation in Abaqus to model intersections between fractures in porous rocks , 2019, Computers and Geotechnics.

[32]  P. Alam.  ‘E’ , 2021, Composites Engineering: An A–Z Guide.

[33]  R. Gracie,et al.  Enriched mixed finite element models for dynamic analysis of continuous and fractured porous media , 2019, Computer Methods in Applied Mechanics and Engineering.

[34]  R. Gracie,et al.  XFEM simulation of a mixed-mode fracture experiment in PMMA , 2020 .

[35]  A. Jafari,et al.  Numerical analysis of multiple hydro-fracture growth in layered media based on a non-differentiable energy minimization approach , 2020 .

[36]  Amir R. Khoei,et al.  A mesh-independent finite element formulation for modeling crack growth in saturated porous media based on an enriched-FEM technique , 2014, International Journal of Fracture.

[37]  P. Bouchard,et al.  Numerical modelling of crack propagation: automatic remeshing and comparison of different criteria , 2003 .

[38]  Ted Belytschko,et al.  A finite element method for crack growth without remeshing , 1999 .

[39]  Kamran Goshtasbi,et al.  3D Numerical Modeling of the Propagation of Hydraulic Fracture at Its Intersection with Natural (Pre-existing) Fracture , 2017, Rock Mechanics and Rock Engineering.

[40]  Eugenio Giner,et al.  An Abaqus implementation of the extended finite element method , 2009 .

[41]  Matthew G. Pike,et al.  XFEM modeling of short microfiber reinforced composites with cohesive interfaces , 2015 .

[42]  M. Kabir,et al.  Experimental and numerical fatigue crack growth of an aluminium pipe repaired by composite patch , 2017 .

[43]  R. Borja.  Assumed enhanced strain and the extended finite element methods: A unification of concepts , 2008 .

[44]  N. Khalili,et al.  Robust simulation of dynamic fluid-driven fracture in naturally fractured impermeable media , 2019 .

[45]  T. Belytschko,et al.  Extended finite element method for three-dimensional crack modelling , 2000 .

[46]  Huang Yuan,et al.  On damage accumulations in the cyclic cohesive zone model for XFEM analysis of mixed-mode fatigue crack growth , 2009 .

[47]  Ted Belytschko,et al.  Discontinuous enrichment in finite elements with a partition of unity method , 2000 .

[48]  Jean-François Remacle,et al.  Substructuring FE-XFE approaches applied to three-dimensional crack propagation , 2008 .

[49]  Amir R. Khoei,et al.  Thermo-hydro-mechanical modeling of impermeable discontinuity in saturated porous media with X-FEM technique , 2012 .

[50]  T. Belytschko,et al.  A review of extended/generalized finite element methods for material modeling , 2009 .

[51]  Julien Réthoré,et al.  A Numerical Approach for Arbitrary Cracks in a Fluid-Saturated Medium , 2006 .

[52]  A. Khoei,et al.  An X-FEM technique in modeling hydro-fracture interaction with naturally-cemented faults , 2019, Engineering Fracture Mechanics.

[53]  A. Jafari,et al.  Fully coupled XFEM formulation for hydraulic fracturing simulation based on a generalized fluid leak-off model , 2021 .

[54]  Bernhard A. Schrefler,et al.  XFEM, strong discontinuities and second-order work in shear band modeling of saturated porous media , 2018, Acta Geotechnica.

[55]  René de Borst,et al.  Mesh-independent discrete numerical representations of cohesive-zone models , 2006 .

[56]  Amir R. Khoei,et al.  An efficient enrichment strategy for modeling stress singularities in isotropic composite materials with X-FEM technique , 2017 .

[57]  Ronaldo I. Borja,et al.  A contact algorithm for frictional crack propagation with the extended finite element method , 2008 .

[58]  Timon Rabczuk,et al.  Phase-field modeling of fluid-driven dynamic cracking in porous media , 2019, Computer Methods in Applied Mechanics and Engineering.

[59]  Vinh Phu Nguyen,et al.  Architecture tradeoffs of integrating a mesh generator to partition of unity enriched object-oriented finite element software , 2007 .

[60]  R. Gracie,et al.  On the undrained and drained hydraulic fracture splits , 2019, International Journal for Numerical Methods in Engineering.

[61]  Amir R. Khoei,et al.  General framework for dynamic large deformation contact problems based on phantom-node X-FEM , 2018 .