Numerical simulation of shock wave propagation in spatially-resolved particle systems

The shock compression of spatially-resolved particle systems is studied at the mesoscale through a series of finite element simulations. The simulations involve propagating shock waves through aluminium–iron oxide thermite systems (Al+Fe2O3) composed of micron-size particles suspended in a polymer binder. Shock-induced chemical reactions are not considered in this work; the particle systems are modelled as inert mixtures. Eulerian formulations are used to accommodate the highly dynamic nature of particulate shock compression. The stress–strain responses of the constituent phases are modelled explicitly at high strain rates and elevated temperatures. Dynamic behaviour of the model system is computed for a set of mixtures (20% and 50% epoxy content by weight) subjected to a range of loading conditions (particle velocities that span 0.300–1.700 km s−1). Spatial profiles of pressure and temperature obtained from the numerical simulations provide insight into thermomechanical responses at the particle level; such resolution is not available in experiments. Finally, Hugoniot data are calculated for the particle mixtures. Stationary pressure calculations are in excellent agreement with experiments, while shock velocity calculations exhibit larger deviations due to the 2D approximation of the microstructure.

[1]  D. Benson,et al.  Numerical simulation of the shock compaction of copper powder , 1994 .

[2]  B. Zhou,et al.  Constitutive Behavior of Epon 828/T-403 at Various Strain Rates , 1998 .

[3]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[4]  L. S. Bennett,et al.  Constitutive model of shock‐induced chemical reactions in inorganic powder mixtures , 1994 .

[5]  T. Vreeland,et al.  Correlation of shock initiated and thermally initiated chemical reactions in a 1:1 atomic ratio nickel‐silicon mixture , 1991 .

[6]  M. Boyce,et al.  A constitutive model for the nonlinear viscoelastic viscoplastic behavior of glassy polymers , 1995 .

[7]  D. Benson.  Computational methods in Lagrangian and Eulerian hydrocodes , 1992 .

[8]  R. Bradt,et al.  Microhardness anisotropy and the indentation size effect on the basal plane of single crystal hematite , 2002 .

[9]  M. Ostoja-Starzewski.  Wavefront propagation in a class of random microstructures—1. bilinear elastic grains☆ , 1991 .

[10]  I ScottKirkpatrick.  Optimization by Simulated Annealing: Quantitative Studies , 1984 .

[11]  N. Thadhani.  Shock-induced chemical reactions and synthesis of materials , 1993 .

[12]  N. Thadhani.  SHOCK COMPRESSION PROCESSING OF POWDERS , 1988 .

[13]  M. Boslough.  A thermochemical model for shock‐induced reactions (heat detonations) in solids , 1990 .

[14]  David J. Benson,et al.  THE CALCULATION OF THE SHOCK VELOCITY: PARTICLE VELOCITY RELATIONSHIP FOR A COPPER POWDER BY DIRECT NUMERICAL SIMULATION , 1995 .

[15]  Farrokh Mistree,et al.  An Approach for Robust Design of Reactive Power Metal Mixtures Based on Non-deterministic Micro-scale Shock Simulation , 2005 .

[16]  David J. Benson,et al.  Eulerian finite-element simulations of experimentally acquired HMX microstructures , 1999 .

[17]  Henry Lee,et al.  Handbook of Epoxy Resins , 1967 .

[18]  D. Benson,et al.  Micromechanical modeling of shock-induced chemical reactions in heterogeneous multi-material powder mixtures , 2001 .

[19]  S. Marsh.  Lasl Shock Hugoniot Data , 1980 .

[20]  M. Meyers,et al.  Effect of shock pressure and plastic strain on chemical reactions in NbSi and MoSi systems , 1995 .

[21]  N. Thadhani,et al.  Shock compression of Al+Fe2O3 powder mixtures of different volumetric distributions , 1998 .

[22]  R. D. Dick,et al.  Equation of State of Aluminum — Iron Oxide (Fe2O3) — Epoxy Composite: Modeling and Experiment , 2006 .

[23]  K. Vecchio,et al.  Shock synthesis of silicides—I. experimentation and microstructural evolution , 1994 .

[24]  D. V. Widder,et al.  Note on Tschebycheff Approximation , 1927 .

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

[26]  L YoungsD,et al.  Time-dependent multi-material flow with large fluid distortion. , 1982 .

[27]  N. Thadhani,et al.  Time-resolved measurements of the shock-compression response of Mo+2Si elemental powder mixtures , 2003 .

[28]  Francis D. Murnaghan,et al.  Finite Deformation of an Elastic Solid , 1967 .

[29]  M. Meyers,et al.  Controlled high‐rate localized shear in porous reactive media , 1994 .

[30]  D. Benson.  A mixture theory for contact in multi-material Eulerian formulations , 1997 .

[31]  N. Thadhani,et al.  Shock-induced chemical reactions in titanium-silicon powder mixtures of different morphologies: Time-resolved pressure measurements and materials analysis , 1997 .

[32]  R. Williamson.  Parametric studies of dynamic powder consolidation using a particle-level numerical model , 1990 .

[33]  David J. Benson,et al.  Particle-level modeling of dynamic consolidation of Ti-SiC powders , 1995 .

[34]  A. M. Wahl.  Finite deformations of an elastic solid: by Francis D. Murnaghan. 140 pages, 15 × 23 cm. New York, John Wiley & Sons, Inc., 1951. Price, $4.00 , 1952 .

[35]  G. Tan,et al.  Modeling of Constitutive Behavior for Epon 828/T-403 at High Strain Rates , 2001 .

[36]  D. Benson.  A multi-material Eulerian formulation for the efficient solution of impact and penetration problems , 1995 .

[37]  K. Vecchio,et al.  SHOCK SYNTHESIS OF SILICIDES--II. THERMODYNAMICS AND KINETICS , 1994 .

[38]  Salvatore Torquato,et al.  Random Heterogeneous Media: Microstructure and Improved Bounds on Effective Properties , 1991 .

[39]  David J. Benson,et al.  An analysis by direct numerical simulation of the effects of particle morphology on the shock compaction of copper powder , 1994 .