Finite element implementation of virtual internal bond model for simulating crack behavior

The virtual internal bond (VIB) model has been recently proposed to describe material deformation and failure under both static and dynamic loading. The model is based on the incorporation of a cohesive type law in a hyperelastic framework, and is capable of fracture simulation as a part of the constitutive formulation. However, with an implicit integration scheme, difficulties are often encountered in the finite element implementation of the VIB model due to possible negative eigenvalues of the stiffness matrix. This paper describes the implementation of an explicit integration scheme of the VIB model. Issues pertaining to the implementation, such as mesh size and shape dependence, loading rate dependence, crack initiation and growth characteristics, and solution time are examined. Both quasi-static and dynamic loading cases have been studied. The experimental validation of the VIB model has been done by calibrating the model parameters using the experimental data of Andrews and Kim [Mech. Mater. 29 (1988) 161]. The simulations using the VIB model are shown to agree well with the experimental observations.

[1]  Huajian Gao,et al.  Crack nucleation and growth as strain localization in a virtual-bond continuum , 1998 .

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

[3]  Paulo B. Lourenço,et al.  A plane stress softening plasticity model for orthotropic materials , 1997 .

[4]  J. Hutchinson,et al.  The influence of plasticity on mixed mode interface toughness , 1993 .

[5]  J. Hutchinson,et al.  The relation between crack growth resistance and fracture process parameters in elastic-plastic solids , 1992 .

[6]  C. Shih,et al.  Ductile crack growth-I. A numerical study using computational cells with microstructurally-based length scales , 1995 .

[7]  E. W Andrews,et al.  Threshold conditions for dynamic fragmentation of glass particles , 1998 .

[8]  M. Ortiz,et al.  Quasicontinuum analysis of defects in solids , 1996 .

[9]  Huajian Gao,et al.  A theory of local limiting speed in dynamic fracture , 1996 .

[10]  F. Milstein.  Theoretical elastic behaviour of crystals at large strains , 1980 .

[11]  G. I. Barenblatt.  The formation of equilibrium cracks during brittle fracture. General ideas and hypotheses. Axially-symmetric cracks , 1959 .

[12]  Alan Needleman,et al.  Oscillatory crack growth in glass , 1999 .

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

[14]  A. Needleman.  A Continuum Model for Void Nucleation by Inclusion Debonding , 1987 .

[15]  R. Hill.  Acceleration waves in solids , 1962 .

[16]  K. Hsia,et al.  Fracture Simulation Using an Elasto-Viscoplastic Virtual Internal Bond Model With Finite Elements , 2004 .

[17]  R. H. Dodds,et al.  Three-dimensional modeling of ductile crack growth in thin sheet metals: computational aspects and validation , 1999 .

[18]  David H. Allen,et al.  Formulation of a three-dimensional cohesive zone model for application to a finite element algorithm , 2000 .

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

[20]  J. Willis,et al.  A comparison of the fracture criteria of griffith and barenblatt , 1967 .

[21]  Huajian Gao,et al.  Numerical simulation of crack growth in an isotropic solid with randomized internal cohesive bonds , 1998 .

[22]  R. Larsson,et al.  A generalized fictitious crack model based on plastic localization and discontinuous approximation , 1995 .

[23]  Huajian Gao.  Elastic waves in a hyperelastic solid near its plane-strain equibiaxial cohesive limit , 1997 .