Search algorithm, and simulation of elastodynamic crack propagation by modified smoothed particle hydrodynamics (MSPH) method

We first present a nonuniform box search algorithm with length of each side of the square box dependent on the local smoothing length, and show that it can save up to 70% CPU time as compared to the uniform box search algorithm. This is especially relevant for transient problems in which, if we enlarge the sides of boxes, we can apply the search algorithm fewer times during the solution process, and improve the computational efficiency of a numerical scheme. We illustrate the application of the search algorithm and the modified smoothed particle hydrodynamics (MSPH) method by studying the propagation of cracks in elastostatic and elastodynamic problems. The dynamic stress intensity factor computed with the MSPH method either from the stress field near the crack tip or from the J-integral agrees very well with that computed by using the finite element method. Three problems are analyzed. One of these involves a plate with a centrally located crack, and the other with three cracks on plates’s horizontal centroidal axis. In each case the plate edges parallel to the crack are loaded in a direction perpendicular to the crack surface. It is found that, at low strain rates, the presence of other cracks will delay the propagation of the central crack. However, at high strain rates, the speed of propagation of the central crack is unaffected by the presence of the other two cracks. In the third problem dealing with the simulation of crack propagation in a functionally graded plate, the crack speed is found to be close to the experimental one.

[1]  Y. M. Chen,et al.  Numerical computation of dynamic stress intensity factors by a Lagrangian finite-difference method (the HEMP code) , 1975 .

[2]  L. Lucy.  A numerical approach to the testing of the fission hypothesis. , 1977 .

[3]  Shigeru Aoki,et al.  Elastodynamic analysis of crack by finite element method using singular element , 1978 .

[4]  S. Atluri,et al.  Numerical Modeling of Dynamic Crack Propagation in Finite Bodies, by Moving Singular Elements—Part 1: Formulation , 1980 .

[5]  J. W. Eischen,et al.  Fracture of nonhomogeneous materials , 1987, International Journal of Fracture.

[6]  L. Hernquist,et al.  TREESPH: A Unification of SPH with the Hierarchical Tree Method , 1989 .

[7]  J. Bonet,et al.  An alternating digital tree (ADT) algorithm for 3D geometric searching and intersection problems , 1991 .

[8]  S. L. Burchell,et al.  Variation of spherical element packing angle and its influence on computer simulations of blasting induced rock motion , 1992 .

[9]  J. Ballmann,et al.  Re-consideration of Chen's problem by finite difference method , 1993 .

[10]  L. Libersky,et al.  High strain Lagrangian hydrodynamics: a three-dimensional SPH code for dynamic material response , 1993 .

[11]  T. Belytschko,et al.  Element‐free Galerkin methods , 1994 .

[12]  D. Rooke,et al.  The dual boundary element method:Ĵ-integral for dynamic stress intensity factors , 1994 .

[13]  Larsgunnar Nilsson,et al.  The position code algorithm for contact searching , 1994 .

[14]  Wing Kam Liu,et al.  Reproducing kernel particle methods , 1995 .

[15]  Brian Mirtich,et al.  Impulse-based dynamic simulation of rigid body systems , 1996 .

[16]  Mark A Fleming,et al.  Continuous meshless approximations for nonconvex bodies by diffraction and transparency , 1996 .

[17]  Zhongmin Jin,et al.  Some basic fracture mechanics concepts in functionally graded materials , 1996 .

[18]  D. Physics,et al.  A Comparison between the Fast Multipole Algorithm and the Tree-Code to Evaluate Gravitational Forces in 3-D , 1997, astro-ph/9703122.

[19]  A. Munjiza,et al.  NBS contact detection algorithm for bodies of similar size , 1998 .

[20]  Carsten Franke,et al.  Convergence order estimates of meshless collocation methods using radial basis functions , 1998, Adv. Comput. Math..

[21]  S. Atluri,et al.  A new Meshless Local Petrov-Galerkin (MLPG) approach in computational mechanics , 1998 .

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

[23]  Romesh C. Batra,et al.  DETERMINATION OF CRACK TIP FIELDS IN LINEAR ELASTOSTATICS BY THE MESHLESS LOCAL PETROV-GALERKIN (MLPG) METHOD , 2001 .

[24]  S. Atluri,et al.  The Meshless Local Petrov-Galerkin (MLPG) Method: A Simple \& Less-costly Alternative to the Finite Element and Boundary Element Methods , 2002 .

[25]  Romesh C. Batra,et al.  ANALYSIS OF ELASTODYNAMIC DEFORMATIONS NEAR A CRACK-NOTCH TIP BY THE MESHLESS LOCAL PETROV–GALERKIN (MLPG) METHOD , 2002 .

[26]  S. Atluri,et al.  The meshless local Petrov-Galerkin (MLPG) method , 2002 .

[27]  John F. Wallin,et al.  A Performance Comparison of Tree Data Structures for N-Body Simulation , 2002 .

[28]  John A. Nairn,et al.  Material Point Method Calculations with Explicit Cracks , 2003 .

[29]  Paulo Sollero,et al.  Dual reciprocity boundary element method in Laplace domain applied to anisotropic dynamic crack problems , 2003 .

[30]  Romesh C. Batra,et al.  Modified smoothed particle hydrodynamics method and its application to transient problems , 2004 .

[31]  Toshio Nakamura,et al.  Simulations of crack propagation in elastic-plastic graded materials , 2004 .

[32]  R. Ogden,et al.  On formulas for the Rayleigh wave speed , 2004 .

[33]  A. Shukla,et al.  Dynamic damage growth in particle reinforced graded materials , 2004 .

[34]  Young-Cheol Yoon,et al.  Numerical prediction of crack propagation by an enhanced element-free Galerkin method , 2004 .

[35]  Romesh C. Batra,et al.  Analysis of adiabatic shear bands in elasto-thermo-viscoplastic materials by modified smoothed-particle hydrodynamics (MSPH) method , 2004 .

[36]  R. Batra,et al.  Crack propagation due to brittle and ductile failures in microporous thermoelastoviscoplastic functionally graded materials , 2005 .

[37]  A. Shukla,et al.  Mixed Mode Dynamic Fracture in Particulate Reinforced Functionally Graded Materials , 2006 .

[38]  R. Batra,et al.  Crack Tip Fields in Functionally Graded Materials with Temperature-Dependent Properties , 2006 .

[39]  Romesh C. Batra,et al.  Wave propagation in functionally graded materials by modified smoothed particle hydrodynamics (MSPH) method , 2007, J. Comput. Phys..