Simulation of high explosive explosion using adaptive material point method

Numerical simulation of high explosive explosion problems is a big challenge to traditional numerical methods because explosion usually involves ex- tremely large deformation and multi-material interaction of different phases. Re- centlydevelopedmeshfreemethodsshowmuchadvantagesovermesh-basedmethod for problems associated with very large deformation. Some of them have been successfully applied to impact and explosion problems, such as smoothed particle hydrodynamics (SPH). Similar to SPH, material point method (MPM) is an effi- cient meshfree particle method solving continuum problems. With combination of the advantages of Eulerian and Lagrangian methods, MPM is a promising numeri- cal tool for solving large deformation problems, such as high explosive detonation and consequent demolishment to the structures. A three dimensional MPM code, MPM3DPP, is developed by using C++ programming language. With adaptive particle splitting scheme proposed in this paper, MPM3DPP is capable of simulat- ing different explosion problems. Johnson-Cook material model is implemented in order to take strain rate effect and thermal softening effect into consideration. Mie- Gruneisen equation of state is used to treat volumetric response of metal under high pressure. Jones-Wilkins-Lee (JWL) equation of state is used for describing the ex- pansion process of detonation products. Artificial viscosity is added to pressure term to stabilize and capture the shock wave. The MPM3DPP code is validated by simulating TNT slab detonation and shock tube problem, and then is used to simulate different explosion problems including explosively driven flyer problem and shaped charge problem. The computational results are in good agreement with empirical formula and experimental results.

[1]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .

[2]  Xiong Zhang,et al.  Comparison study of MPM and SPH in modeling hypervelocity impact problems , 2009 .

[3]  Jianshi Zhu,et al.  Three-dimensional multi-mesh material point method for solving collision problems , 2007, 0708.3532.

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

[5]  Guirong Liu,et al.  Computer simulation of high explosive explosion using smoothed particle hydrodynamics methodology , 2003 .

[6]  Ranga Komanduri,et al.  Simulation of dynamic crack growth using the generalized interpolation material point (GIMP) method , 2007 .

[7]  William P. Walters,et al.  Explosive effects and applications , 1998 .

[8]  Hans Muhlhaus,et al.  A Lagrangian integration point finite element method for large deformation modeling of viscoelastic geomaterials , 2003 .

[9]  Wing Kam Liu,et al.  Nonlinear Finite Elements for Continua and Structures , 2000 .

[10]  James E. Guilkey,et al.  Improved Velocity Projection for the Material Point Method , 2007 .

[11]  Xiong Zhang,et al.  An explicit material point finite element method for hyper‐velocity impact , 2006 .

[12]  John A. Nairn,et al.  Three-Dimensional Dynamic Fracture Analysis Using the Material Point Method , 2006 .

[13]  S. Bardenhagen,et al.  The Generalized Interpolation Material Point Method , 2004 .

[14]  S. Atluri,et al.  Meshless Local Petrov-Galerkin (MLPG) Approaches for Solving Nonlinear Problems with Large Deformations and Rotations , 2005 .

[15]  James E. Guilkey,et al.  An Improved Contact Algorithm for the Material Point Method and Application to Stress Propagation in Granular Material , 2001 .

[16]  K. Y. Sze,et al.  Smoothed molecular dynamics for large step time integration , 2007 .

[17]  James J. Mason,et al.  On the strain and strain rate dependence of the fraction of plastic work converted to heat: an experimental study using high speed infrared detectors and the Kolsky bar☆ , 1992 .

[18]  Ranga Komanduri,et al.  Multiscale Simulation of Nanoindentation Using the Generalized Interpolation Material Point (GIMP) Method, Dislocation Dynamics (DD) and Molecular Dynamics (MD) , 2006 .

[19]  Deborah Sulsky,et al.  An unconditionally stable, energy–momentum consistent implementation of the material-point method , 2006 .

[20]  Ronald E. Brown,et al.  Building characteristics into a shaped charge to achieve unique performance requirements , 1995 .

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

[22]  J. Brackbill,et al.  Flip: A low-dissipation, particle-in-cell method for fluid flow , 1988 .

[23]  Wen Ho Lee.  Computer Simulation of Shaped Charge Problems , 2006 .

[24]  Mhamed Souli,et al.  High explosive simulation using multi-material formulations , 2006 .

[25]  Wenqing Hu,et al.  Model-based simulation of the synergistic effects of blast and fragmentation on a concrete wall using the MPM , 2006 .

[26]  H. Udaykumar,et al.  A particle-level set-based sharp interface cartesian grid method for impact, penetration, and void collapse , 2004 .

[27]  J. Brackbill,et al.  The material-point method for granular materials , 2000 .

[28]  Jean-François Molinari,et al.  Finite element simulation of shaped charges , 2002 .

[29]  Guirong Liu,et al.  Meshfree particle simulation of the detonation process for high explosives in shaped charge unlined cavity configurations , 2003 .

[30]  Mark A Fleming,et al.  Meshless methods: An overview and recent developments , 1996 .

[31]  Wing Kam Liu,et al.  Meshfree and particle methods and their applications , 2002 .

[32]  Ning Jianguo,et al.  Fuzzy interface treatment in Eulerian method , 2004 .

[33]  S. Atluri.  The meshless method (MLPG) for domain & BIE discretizations , 2004 .

[34]  Ning Jianguo,et al.  Multi-material Eulerian Formulations and Hydrocode for the Simulation of Explosions , 2008 .

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

[36]  Satya N. Atluri,et al.  The Applications of Meshless Local Petrov-Galerkin (MLPG) Approaches in High-Speed Impact, Penetration and Perforation Problems , 2006 .

[37]  G. R. Johnson,et al.  Evaluation of cylinder‐impact test data for constitutive model constants , 1988 .

[38]  L. Libersky,et al.  Smoothed Particle Hydrodynamics: Some recent improvements and applications , 1996 .

[39]  G. R. Johnson,et al.  SPH for high velocity impact computations , 1996 .

[40]  M. Meyers.  Dynamic Behavior of Materials , 1994 .

[41]  James E. Chisum,et al.  Modeling and Simulation of Underwater Shock Problems Using a Coupled Lagrangian-Eulerian Analysis Approach , 1997 .

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