Discrete element simulation of shock wave propagation in polycrystalline copper

Abstract

The implementation of the characteristic of compressive plasticity into the Discrete Element Code, DM2, while maintaining its quasi-molecular scheme, is described. The code is used to simulate the shock compression of polycrystalline copper at 3.35 and 11.0 GPa. The model polycrystal has a normal distribution of grain sizes, with mean diameter 14 μm, and three distinct grain orientations are permitted with respect to the shock direction; 〈1 0 0〉, 〈1 1 0〉, and 〈1 1 1〉. Particle velocity dispersion (PVD) is present in the shock-induced flow, attaining its maximum magnitude at the plastic wave rise. PVD normalised to the average particle velocity of Δ u y / u ¯ P = 0.05 and Δ u y / u ¯ P = 0.087 are yielded for the 3.35 and 11.0 GPa shocks, respectively, and are of the same order as those seen in the experiment. Non-planar elastic and plastic wave fronts are present, the distribution in shock front position increasing with propagation distance. The rate of increase of the spread in shock front positions is found to be significantly smaller than that seen in probabilistic calculations on nickel polycrystals, and this difference is attributed, in the main, to grain interaction. Reflections at free surfaces yield a region of tension near to the target free surface. Due to the dispersive nature of the shock particle velocity and the non-planarity of the shock front, the tensile pressure is distributed. This may have implications for the spall strength, which are discussed. Simulations reveal a transient shear stress distribution behind the shock front. Such a distribution agrees with that put forward by Lipkin and Asay to explain the quasi-elastic reloading phenomenon. Simulation of reloading shocks show that the shear stress distribution can give rise to quasi-elastic reloading on the grain scale.