Simple and effective approach to modeling crack propagation in the framework of extended finite element method

Abstract In this study, the extended finite element method (XFEM) coupled with the hierarchical mesh adaptation method and direct method is developed to model crack propagation. First, we observe that crack tip stresses can be approximated with high accuracy by using the stress intensity factor (SIF) terms and T-stress term within 0.1 % of the crack length away from the crack tip, without considering other high order terms of the crack-tip stresses. Next, based on the XFEM, the hierarchical mesh adaptation method is developed to obtain the stresses around the location which is very close to the crack tip. Then, the SIFs are derived directly from the crack-tip stresses based on the least square fitting without considering the special algorithms and procedures. Finally, the numerical results reveal that the proposed method can obtain SIFs accurately and it can predict a better crack trajectory for crack simulation. In addition, the proposed method has the advantage over the J-integral method and it is benefit for the simulation of large-scale engineering problems.

[1]  A. Needleman,et al.  A COMPARISON OF METHODS FOR CALCULATING ENERGY RELEASE RATES , 1985 .

[2]  M. Williams,et al.  On the Stress Distribution at the Base of a Stationary Crack , 1956 .

[3]  Surendra P. Shah,et al.  Mixed-mode fracture of concrete , 1988, International Journal of Fracture.

[4]  M. Ayatollahi,et al.  Mixed mode brittle fracture in PMMA—An experimental study using SCB specimens , 2006 .

[5]  Herbert Edelsbrunner,et al.  Geometry and Topology for Mesh Generation , 2001, Cambridge monographs on applied and computational mathematics.

[6]  Bostjan Zafosnik,et al.  Određivanje koeficijenta intenziteta naprezanja KI metodom ekstrapolacije , 2016 .

[7]  Haim Waisman,et al.  Direct evaluation of stress intensity factors for curved cracks using Irwin's integral and XFEM with high‐order enrichment functions , 2017 .

[8]  Grégory Legrain,et al.  Robust and direct evaluation of J2 in linear elastic fracture mechanics with the X‐FEM , 2008 .

[9]  Filippo Berto,et al.  A comparison among some recent energy- and stress-based criteria for the fracture assessment of sharp V-notched components under Mode I loading , 2014 .

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

[11]  J. Rice.  A path-independent integral and the approximate analysis of strain , 1968 .

[12]  Stefano Mariani,et al.  Extended finite element method for quasi‐brittle fracture , 2003 .

[13]  N. Moës,et al.  Improved implementation and robustness study of the X‐FEM for stress analysis around cracks , 2005 .

[14]  Harish Kumar,et al.  Progressive damage response and crack growth direction for multiple through cracks of laminated composite finite plate , 2018 .

[15]  James C. Newman,et al.  Stress-intensity factors for a wide range of semi-elliptical surface cracks in finite-thickness plates , 1979 .

[16]  M.R.M. Aliha,et al.  Size-dependent fracture behavior of Guiting limestone under mixed mode loading , 2014 .

[17]  Wataru Fujisaki,et al.  An Determination Method for Stress Intensity Factor by Extrapolation from an Arbitrary Angle under KI, KII Mixture Mode. , 1999 .

[18]  Jun-Hyok Ri,et al.  Limit analysis of cracked structure by combination of extended finite element method with linear matching method , 2018 .

[19]  Bijay K. Mishra,et al.  The numerical simulation of fatigue crack growth using extended finite element method , 2012 .

[20]  Timon Rabczuk,et al.  Modeling and simulation of kinked cracks by virtual node XFEM , 2015 .

[21]  Pierre Kerfriden,et al.  Minimum energy multiple crack propagation. Part III: XFEM computer implementation and applications , 2017 .

[22]  M.R.M. Aliha,et al.  Geometry and size effects on fracture trajectory in a limestone rock under mixed mode loading , 2010 .

[23]  Andrei Kotousov,et al.  Three-dimensional stress states at crack tip induced by shear and anti-plane loading , 2013 .

[24]  Shuodao Wang,et al.  A Mixed-Mode Crack Analysis of Isotropic Solids Using Conservation Laws of Elasticity , 1980 .

[25]  Per-Olof Persson,et al.  A Simple Mesh Generator in MATLAB , 2004, SIAM Rev..

[26]  Ted Belytschko,et al.  An extended finite element method for modeling crack growth with frictional contact , 2001 .

[27]  D. Kershaw.  The incomplete Cholesky—conjugate gradient method for the iterative solution of systems of linear equations , 1978 .

[28]  Toshiaki Murakami,et al.  Proportional extrapolation techniques for determining stress intensity factors , 1988 .

[29]  I. Harari,et al.  A direct analytical method to extract mixed‐mode components of strain energy release rates from Irwin's integral using extended finite element method , 2013 .

[30]  Surendra P. Shah,et al.  Fracture of Concrete Subjected to Impact Loading , 1986 .

[31]  Filippo Berto,et al.  The improvement of crack propagation modelling in triangular 2D structures using the extended finite element method , 2018, Fatigue & Fracture of Engineering Materials & Structures.

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

[33]  Bhushan Lal Karihaloo,et al.  XFEM for direct evaluation of mixed mode SIFs in homogeneous and bi‐materials , 2004 .

[34]  M. Ayatollahi,et al.  Analysis of a new specimen for mixed mode fracture tests on brittle materials , 2009 .

[35]  F. Berto,et al.  On Higher Order Terms in the Crack Tip Stress Field , 2010 .

[36]  H. M. Westergaard.  Stresses At A Crack, Size Of The Crack, And The Bending Of Reinforced Concrete , 1933 .

[37]  Leslie Banks-Sills,et al.  Comparison of methods for calculating stress intensity factors with quarter-point elements , 1986 .

[38]  Hyon-Sik Hong,et al.  XFEM fracture analysis by applying smoothed weighted functions with compact support , 2018 .

[39]  C. Buchalet,et al.  Stress Intensity Factor Solutions for Continuous Surface Flaws in Reactor Pressure Vessels , 1976 .

[40]  Hussain,et al.  Strain Energy Release Rate for a Crack Under Combined Mode I and Mode II , 1974 .

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

[42]  M. Aliha,et al.  Crack growth path prediction for the angled cracked plate using higher order terms of Williams series expansion , 2013 .

[43]  P. C. Paris,et al.  The Stress Analysis of Cracks Handbook, Third Edition , 2000 .

[44]  Naser A. Al-Shayea,et al.  Crack propagation trajectories for rocks under mixed mode I–II fracture , 2005 .

[45]  M. Pavier,et al.  The role of T-stress in brittle fracture for linear elastic materials under mixed-mode loading , 2001 .

[46]  S. Ahmad,et al.  Mixed mode crack growth in elasto-plastic-creeping solids using XFEM , 2018, Engineering Fracture Mechanics.

[47]  N. Subramonian,et al.  The inclined crack under biaxial load , 1978 .

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

[49]  S. Chan,et al.  On the Finite Element Method in Linear Fracture Mechanics , 1970 .

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

[51]  Per Christian Hansen,et al.  REGULARIZATION TOOLS: A Matlab package for analysis and solution of discrete ill-posed problems , 1994, Numerical Algorithms.

[52]  M.R.M. Aliha,et al.  Experimental and Theoretical Study of Fracture Paths in Brittle Cracked Materials Subjected to Pure Mode II Loading , 2011 .

[53]  Jaime Planas,et al.  KI evaluation by the displacement extrapolation technique , 2000 .

[54]  T. Rabczuk.  Computational Methods for Fracture in Brittle and Quasi-Brittle Solids: State-of-the-Art Review and Future Perspectives , 2013 .

[55]  M.R.T. Arruda,et al.  Energetic Convergence of a New Hybrid Mixed Finite Element , 2019, Engineering Solid Mechanics.

[56]  Hamid Hashemolhosseini,et al.  Determination of crack initiation and propagation in two disc shaped specimens using the improved maximum tangential stress criterion , 2017 .

[57]  S. K. Maiti,et al.  Comparison of the criteria for mixed mode brittle fracture based on the preinstability stress-strain field Part I: Slit and elliptical cracks under uniaxial tensile loading , 1983 .

[58]  S. Bordas,et al.  A robust preconditioning technique for the extended finite element method , 2011 .