Cohesive element approach for dynamic crack propagation: Artificial compliance and mesh dependency

Abstract Zero thickness cohesive element approach for arbitrary crack propagation has a deficiency of introducing artificial compliance to the model, especially when cohesive elements are inserted into every element interfaces. For dynamic problems, the artificial compliance decreases the stress wave speed and makes the result less accurate. In this paper, the reason of the artificial compliance is examined, and the bilinear and exponential cohesive law are compared. The work shows that by choosing the right cohesive stiffness, element size and using bilinear cohesive law rather than exponential cohesive law, the artificial compliance issue can be limited to a negligible level without greatly increasing the computational time. Three numerical simulations are used to support the argument. The bilinear cohesive law also shows more robust behavior in terms of cyclic loading than exponential law, and a new cyclic loading formulation is proposed to help fix the discontinuity observed in exponential law. According to our finding, in order to limit the artificial compliance and computational time, large element size is recommended, however, in fracture problems, small element size around crack tip is essential to capture the cohesive zone behavior. To escape the dilemma, we modify the cohesive zone enlargement approach presented in the literature (Dugdale, 1960) and adopt it for arbitrary crack propagation. The modified methodology enlarges cohesive zone size by reducing the cohesive strength only around crack tips to allow more cohesive elements inside the cohesive zone. Two benchmark numerical simulations are carried out to verify the modified methodology.

[1]  Robert O. Ritchie,et al.  Evaluation of toughness in AISI 4340 alloy steel austenitized at low and high temperatures , 1976, Metallurgical and Materials Transactions. A.

[2]  K. Volokh Comparison between cohesive zone models , 2004 .

[3]  C. Suryanarayana,et al.  Effect of clustering on the mechanical properties of SiC particulate-reinforced aluminum alloy 2024 metal matrix composites , 2003 .

[4]  T. Belytschko,et al.  Dynamic crack propagation based on loss of hyperbolicity and a new discontinuous enrichment , 2003 .

[5]  van den Mj Marco Bosch,et al.  An improved description of the exponential Xu and Needleman cohesive zone law for mixed-mode decohesion , 2006 .

[6]  J. G. Williams,et al.  The use of a cohesive zone model to study the fracture of fibre composites and adhesively-bonded joints , 2003 .

[7]  M. Ortiz,et al.  Computational modelling of impact damage in brittle materials , 1996 .

[8]  D. S. Dugdale Yielding of steel sheets containing slits , 1960 .

[9]  H. Nguyen-Xuan,et al.  A simple and robust three-dimensional cracking-particle method without enrichment , 2010 .

[10]  Loïc Daridon,et al.  Artificial compliance inherent to the intrinsic cohesive zone models: criteria and application to planar meshes , 2012, International Journal of Fracture.

[11]  Glaucio H. Paulino,et al.  Simulation of Crack Propagation in Asphalt Concrete Using an Intrinsic Cohesive Zone Model , 2006 .

[12]  Toshio Nakamura,et al.  Simulations of Crack Propagation in Porous Materials , 2001 .

[13]  Pedro P. Camanho,et al.  A damage model for the simulation of delamination in advanced composites under variable-mode loading , 2006 .

[14]  Loïc Daridon,et al.  Micromechanical-based criteria for the calibration of cohesive zone parameters , 2013, J. Comput. Appl. Math..

[15]  P. Camanho,et al.  Mixed-Mode Decohesion Finite Elements for the Simulation of Delamination in Composite Materials , 2002 .

[16]  John R. Rice,et al.  A Critical Evaluation of Cohesive Zone Models of Dynamic Fracture , 2001 .

[17]  Vinh Phu Nguyen An open source program to generate zero-thickness cohesive interface elements , 2014, Adv. Eng. Softw..

[18]  Alojz Ivankovic,et al.  The prediction of dynamic fracture evolution in PMMA using a cohesive zone model , 2005 .

[19]  Laxmikant V. Kalé,et al.  Parallel Simulations of Dynamic Fracture Using Extrinsic Cohesive Elements , 2009, J. Sci. Comput..

[20]  Xiaopeng Xu,et al.  Numerical simulations of fast crack growth in brittle solids , 1994 .

[21]  T. Belytschko,et al.  A comparative study on finite element methods for dynamic fracture , 2008 .

[22]  A. Hillerborg,et al.  Analysis of crack formation and crack growth in concrete by means of fracture mechanics and finite elements , 1976 .

[23]  G. Paulino,et al.  Cohesive Zone Models: A Critical Review of Traction-Separation Relationships Across Fracture Surfaces , 2011 .

[24]  M. Benzeggagh,et al.  Measurement of mixed-mode delamination fracture toughness of unidirectional glass/epoxy composites with mixed-mode bending apparatus , 1996 .

[25]  Carlos Armando Duarte,et al.  The role of cohesive properties on intergranular crack propagation in brittle polycrystals , 2011 .

[26]  H. Espinosa,et al.  A grain level model for the study of failure initiation and evolution in polycrystalline brittle materials. Part I: Theory and numerical implementation , 2003 .

[27]  Julien Réthoré,et al.  Efficient explicit time stepping for the eXtended Finite Element Method (X‐FEM) , 2006 .

[28]  G. Alfano On the influence of the shape of the interface law on the application of cohesive-zone models , 2006 .

[29]  M. Kolluri,et al.  Irreversible mixed mode interface delamination using a combined damage-plasticity cohesive zone enabling unloading , 2013, International Journal of Fracture.

[30]  V. Tomar,et al.  Bounds for element size in a variable stiffness cohesive finite element model , 2004 .

[31]  Ludovic Noels,et al.  A scalable 3D fracture and fragmentation algorithm based on a hybrid, discontinuous Galerkin, Cohesive Element Method , 2011 .

[32]  Stephen R Hallett,et al.  Cohesive zone length in numerical simulations of composite delamination , 2008 .

[33]  Huajian Gao,et al.  Physics-based modeling of brittle fracture: Cohesive formulations and the application of meshfree methods , 2000 .

[34]  Bert Sluys,et al.  Numerical simulation of quasi-brittle fracture using damaging cohesive surfaces , 2000 .

[35]  Lorenzo Iannucci,et al.  Formulation and implementation of decohesion elements in an explicit finite element code , 2006 .

[36]  P. Camanho,et al.  Numerical Simulation of Mixed-Mode Progressive Delamination in Composite Materials , 2003 .

[37]  J. Rice,et al.  The growth of slip surfaces in the progressive failure of over-consolidated clay , 1973, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[38]  F. Bobaru,et al.  Studies of dynamic crack propagation and crack branching with peridynamics , 2010 .

[39]  G. I. Barenblatt THE MATHEMATICAL THEORY OF EQUILIBRIUM CRACKS IN BRITTLE FRACTURE , 1962 .

[40]  Pedro P. Camanho,et al.  An engineering solution for mesh size effects in the simulation of delamination using cohesive zone models , 2007 .