Crack growth modeling via 3D automatic adaptive mesh refinement based on modified-SPR technique

Abstract In this paper, the three-dimensional automatic adaptive mesh refinement is presented in modeling the crack propagation based on the modified superconvergent patch recovery technique. The technique is developed for the mixed mode fracture analysis of different fracture specimens. The stress intensity factors are calculated at the crack tip region and the crack propagation is determined by applying a proper crack growth criterion. An automatic adaptive mesh refinement is employed on the basis of modified superconvergent patch recovery (MSPR) technique to simulate the crack growth by applying the asymptotic crack tip solution and using the collapsed quarter-point singular tetrahedral elements at the crack tip region. A-posteriori error estimator is used based on the Zienkiewicz–Zhu method to estimate the error of fracture parameters and predict the crack path pattern. Finally, the efficiency and accuracy of proposed computational algorithm is demonstrated by several numerical examples.

[1]  G. Irwin.  Crack-Extension Force for a Part-Through Crack in a Plate , 1962 .

[2]  J. Z. Zhu,et al.  The superconvergent patch recovery and a posteriori error estimates. Part 1: The recovery technique , 1992 .

[3]  M. Kanninen,et al.  A finite element calculation of stress intensity factors by a modified crack closure integral , 1977 .

[4]  Jörn Mosler,et al.  On the modeling of highly localized deformations induced by material failure: The strong discontinuity approach , 2004 .

[5]  Ted Belytschko,et al.  A vector level set method and new discontinuity approximations for crack growth by EFG , 2002 .

[6]  Amir R. Khoei,et al.  Three-dimensional data transfer operators in large plasticity deformations using modified-SPR technique , 2009 .

[7]  F. A. Brogan,et al.  Some computational tools for the analysis of through cracks in stiffened fuselage shells , 1993 .

[8]  S. Atluri,et al.  CALCULATION OF FRACTURE MECHANICS PARAMETERS FOR AN ARBITRARY THREE-DIMENSIONAL CRACK, BY THE ‘EQUIVALENT DOMAIN INTEGRAL’ METHOD , 1987 .

[9]  Julien Yvonnet,et al.  An XFEM/level set approach to modelling surface/interface effects and to computing the size-dependent effective properties of nanocomposites , 2008 .

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

[11]  R. de Borst,et al.  Analysis of concrete fracture in “direct” tension , 1989 .

[12]  Saeed Asil Gharebaghi,et al.  SUT-DAM: An integrated software environment for multi-disciplinary geotechnical engineering , 2006, Adv. Eng. Softw..

[13]  de R René Borst,et al.  Combination of fracture and damage mechanics for numerical failure analysis , 1999 .

[14]  Guido Dhondt,et al.  Computational fracture analysis of different specimens regarding 3D and mode coupling effects , 2001 .

[15]  Jean-Baptiste Leblond,et al.  Coplanar propagation paths of 3D cracks in infinite bodies loaded in shear , 2006 .

[16]  Robert W. Zimmerman,et al.  Numerical simulation of multiple 3D fracture propagation using arbitrary meshes , 2011 .

[17]  T. K. Hellen.  On the method of virtual crack extensions , 1975 .

[18]  A. Bolukbasi,et al.  Analytical and Experimental Studies of Crushing Behavior in Composite Laminates , 1995 .

[19]  Y. R. Rashid,et al.  Ultimate strength analysis of prestressed concrete pressure vessels , 1968 .

[20]  Ercan Gürses,et al.  A computational framework of three-dimensional configurational-force-driven brittle crack propagation , 2009 .

[21]  Amir R. Khoei,et al.  Three-dimensional superconvergent patch recovery method and its application to data transferring in small-strain plasticity , 2007 .

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

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

[24]  Victor E. Saouma,et al.  Numerical evaluation of 3-D SIF for arbitrary finite element meshes , 1997 .

[25]  M. Kawahara,et al.  A three-step finite element method for unsteady incompressible flows , 1993 .

[26]  Ted Belytschko,et al.  Dislocations by partition of unity , 2005 .

[27]  Paul J. Fanning,et al.  Smeared crack models of RC beams with externally bonded CFRP plates , 2000 .

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

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

[30]  Amir R. Khoei,et al.  Error estimation, adaptivity and data transfer in enriched plasticity continua to analysis of shear band localization , 2007 .

[31]  Paul A. Wawrzynek,et al.  Automated 3‐D crack growth simulation , 2000 .

[32]  Hans Albert Richard,et al.  Development of a new software for adaptive crack growth simulations in 3D structures , 2003 .

[33]  Amir R. Khoei,et al.  The superconvergence patch recovery technique and data transfer operators in 3D plasticity problems , 2007 .

[34]  Alberto Carpinteri,et al.  Numerical analysis of fracture mechanisms and failure modes in bi-layered structural components , 2007 .

[35]  Amir R. Khoei,et al.  3D adaptive finite element modeling of non-planar curved crack growth using the weighted superconvergent patch recovery method , 2009 .

[36]  M. K. Kassir,et al.  Three-Dimensional Stress Distribution Around an Elliptical Crack Under Arbitrary Loadings , 1966 .

[37]  Jin Wu,et al.  Development of the DYNA3D simulation code with automated fracture procedure for brick elements , 2003 .

[38]  Amir R. Khoei,et al.  H-adaptive mesh refinement for shear band localization in elasto-plasticity Cosserat continuum , 2005 .

[39]  Sankaran Mahadevan,et al.  Threshold stress intensity factor and crack growth rate prediction under mixed-mode loading , 2007 .