MPM simulation of dynamic material failure with a decohesion constitutive model

An initiation criterion is developed for a discrete constitutive equation model of material failure. The decohesion constitutive model is combined with elasticity and implemented in the numerical material-point method. The complete numerical proceedure is used to study spall failure in brittle materials. A one-dimensional analytical solution for a bar subjected to a sinusoidal pulse is derived and used to verify the numerical method. The sinusoidal pulse leads to a finite region of the bar undergoing failure, rather than a single plane, as occurs with a square or triangular wave pulse. This observation may be related to experimental findings of multiple spall planes for impacted bars. The numerical method is not restricted to one-dimensional problems, thus we further study the additional features that arise for an impacted bar assuming that the bar deforms under plane stress conditions. Reflections of stress waves from the lateral surfaces of the bar result in the appearance of curved and secondary spall surfaces.

[1]  J. C. Simo,et al.  An analysis of strong discontinuities induced by strain-softening in rate-independent inelastic solids , 1993 .

[2]  D. Sulsky,et al.  A particle method for history-dependent materials , 1993 .

[3]  Francisco Armero,et al.  On the characterization of localized solutions in inelastic solids: an analysis of wave propagation in a softening bar , 2001 .

[4]  Howard L. Schreyer,et al.  Modeling delamination as a strong discontinuity with the material point method , 2002 .

[5]  U. S. Lindholm,et al.  Shock Wave and High-Strain-Rate Phenomena in Materials , 1992 .

[6]  Kenneth Runesson,et al.  Element-Embedded Localization Band Based on Regularized Displacement Discontinuity , 1996 .

[7]  Lynn Seaman,et al.  Dynamic failure of solids , 1987 .

[8]  D. R. Curran,et al.  A Continuum Model for Dynamic Tensile Microfracture and Fragmentation , 1985 .

[9]  J. Oliver.  On the discrete constitutive models induced by strong discontinuity kinematics and continuum constitutive equations , 2000 .

[10]  K. Kimsey,et al.  Material Modeling for Terminal Ballistic Simulation , 1992 .

[11]  Alberto Corigliano,et al.  Formulation, identification and use of interface models in the numerical analysis of composite delamination , 1993 .

[12]  D. Sulsky.  Erratum: Application of a particle-in-cell method to solid mechanics , 1995 .

[13]  J. Oliver.  MODELLING STRONG DISCONTINUITIES IN SOLID MECHANICS VIA STRAIN SOFTENING CONSTITUTIVE EQUATIONS. PART 2: NUMERICAL SIMULATION , 1996 .

[14]  Howard L. Schreyer,et al.  Axisymmetric form of the material point method with applications to upsetting and Taylor impact problems , 1996 .

[15]  C. Mariotti,et al.  Experimental and numerical study of concrete at high strain rates in tension , 2001 .

[16]  F. Armero,et al.  An analysis of strong discontinuities in multiplicative finite strain plasticity and their relation with the numerical simulation of strain localization in solids , 1996 .

[17]  Donald R. Curran,et al.  Fragmentation of rock under dynamic loads , 1974 .

[18]  Arunachalam M. Rajendran,et al.  Modeling the impact behavior of AD85 ceramic under multiaxial loading , 1994 .

[19]  U. Mayer,et al.  Crack Propagation and Arrest Study at Stress Pulse Loading , 1991 .

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

[21]  L. J. Sluys,et al.  Three-dimensional embedded discontinuity model for brittle fracture , 2001 .

[22]  R. de Borst,et al.  Numerical Study on Crack Dilatancy Part 2: Applications , 1991 .

[23]  Ares J. Rosakis,et al.  Three Dimensional Cohesive-Element Analysis and Experiments of Dynamic Fracture in C300 Steel , 2000 .

[24]  R. de Borst,et al.  Numerical Study on Crack Dilatancy Part 1: Models and Stability Analysis , 1991 .

[25]  M. Ortiz,et al.  FINITE-DEFORMATION IRREVERSIBLE COHESIVE ELEMENTS FOR THREE-DIMENSIONAL CRACK-PROPAGATION ANALYSIS , 1999 .

[26]  F. Armero,et al.  Large‐scale modeling of localized dissipative mechanisms in a local continuum: applications to the numerical simulation of strain localization in rate‐dependent inelastic solids , 1999 .

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

[28]  J. Oliver.  MODELLING STRONG DISCONTINUITIES IN SOLID MECHANICS VIA STRAIN SOFTENING CONSTITUTIVE EQUATIONS. PART 1: FUNDAMENTALS , 1996 .

[29]  L. J. Sluys,et al.  Application of embedded discontinuities for softening solids , 2000 .

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

[31]  F. J. Mello,et al.  Modeling the Initiation and Growth of Delaminations in Composite Structures , 1996 .

[32]  E. Dvorkin,et al.  Finite elements with displacement interpolated embedded localization lines insensitive to mesh size and distortions , 1990 .

[33]  M. E. Kipp,et al.  THE MICROMECHANICS OF IMPACT FRACTURE OF ROCK , 1979 .

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

[35]  Christian Miehe,et al.  Post-critical discontinuous localization analysis of small-strain softening elastoplastic solids , 1994, Archive of Applied Mechanics.