Modeling mixed-mode dynamic crack propagation nsing finite elements: Theory and applications

Previous work in modeling dynamic fracture has assumed the crack will propagate along predefined mesh lines (usually a straight line). In this paper we present a finite element model of mixed-mode dynamic crack propagation in which this constraint is removed. Applying linear elasto-dynamic fracture mechanics concepts, discrete cracks are allowed to propagate through the mesh in arbitrary directions. The fracture criteria used for propagation and the algorithms used for remeshing are described in detail. Important features of the implementation are the use of triangular elements with quadratic shape functions, explicit time integration, and interactive computer graphics. These combine to make the approach robust and applicable to a broad range of problems.Example analyses of straight and curving crack problems are presented. Verification problems include a stationary crack under dynamic loading and a propagating crack in an infinite body. Comparisons with experimental data are made for curving propagation in a cracked plate under biaxial loading.

[1]  Ares J. Rosakis,et al.  The determination of dynamic fracture toughness of AISI 4340 steel by the shadow spot method , 1984 .

[2]  H. P. Rossmanith,et al.  How ‘mixed’ is dynamic mixed-mode crack propagation? — A dynamic photoelastic study , 1983 .

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

[4]  Hyun Moo Koh.  A Mixed Eulerian-Lagrangian Model for the Analysis of Dynamic Fracture (Finite Element Method, Moving Mesh, Virtual Energy Release Rate) , 1986 .

[5]  Samuel W. Key,et al.  HONDO II: a finite element computer program for the large deformation dynamic response of axisymmetric solids. [HONDO] , 1978 .

[6]  J. R. M. Radok.  On the solution of problems of dynamic plane elasticity , 1956 .

[7]  L. B. Freund,et al.  Crack propagation in an elastic solid subjected to general loading—I. Constant rate of extension , 1972 .

[8]  R. Barsoum.  Triangular quarter‐point elements as elastic and perfectly‐plastic crack tip elements , 1977 .

[9]  R. D. Henshell,et al.  CRACK TIP FINITE ELEMENTS ARE UNNECESSARY , 1975 .

[10]  E. Yoffe,et al.  The moving Griffith crack , 1951 .

[11]  L. Freund.  Crack propagation in an elastic solid subjected to general loading—III. Stress wave loading , 1973 .

[12]  L. B. Freund,et al.  On the uniqueness of plane elastodynamic solutions for running cracks , 1974 .

[13]  L. B. Freund,et al.  Crack propagation in an elastic solid subjected to general loading—II. Non-uniform rate of extension , 1972 .

[14]  Satya N. Atluri,et al.  Numerical studies in dynamic fracture mechanics , 1985 .

[15]  Stephen A. Thau,et al.  Transient stress intensity factors for a finite crack in an elastic solid caused by a dilatational wave , 1971 .

[16]  S. Atluri,et al.  Path-independent integrals, energy release rates, and general solutions of near-tip fields in mixed-mode dynamic fracture mechanics , 1983 .

[17]  Jan Drewes Achenbach,et al.  ELASTODYNAMIC NEAR-TIP STRESS AND DISPLACEMENT FIELDS FOR RAPIDLY PROPAGATING CRACKS IN ORTHOTROPIC MATERIALS. , 1975 .

[18]  H. Saunders,et al.  Finite element procedures in engineering analysis , 1982 .

[19]  E. Yoffe,et al.  LXXV. The moving griffith crack , 1951 .

[20]  Albert S. Kobayashi,et al.  Dynamic-finite-element and dynamic-photoelastic analyses of two fracturing homalite-100 plates , 1976 .

[21]  Fred Nilsson,et al.  Dynamic stress-intensity factors for finite strip problems , 1972 .

[22]  T. Belytschko,et al.  Computational Methods for Transient Analysis , 1985 .

[23]  L. B. Freund.  Crack propagation in an elastic solid subjected to general loading— IV. Obliquely incident stress pulse , 1974 .

[24]  R. Shaw,et al.  Modification to the Suhara‐Fukuda method of network generation , 1978 .