Shock compression of monocrystalline copper: Atomistic Simulations

BY Cao and EM Bringa and MA Meyers, METALLURGICAL AND MATERIALS TRANSACTIONS A-PHYSICAL METALLURGY AND MATERIALS SCIENCE, 38A, 2681-2688 (2007).

DOI: 10.1007/s11661-007-9248-9

Molecular dynamics (MD) simulations were used to model the effects of shock compression on 001 and 221 monocrystals. We obtained the Hugoniot for both directions, and analyzed the formation of a two-wave structure for the 221 monocrystal. We also analyzed the dislocation structure induced by the shock compression along these two crystal orientations. The topology of this structure compares extremely well with that observed in recent transmission electron microscopy (TEM) studies of shock-induced plasticity in samples recovered from flyer plate and laser shock experiments. However, the density of stacking faults in our simulations is 10(2) to 10(4) times larger than in the experimental observations of recovered samples. The difference between experimentally observed TEM and calculated MD results is attributed to two effects: (1) the annihilation of dislocations during post-shock relaxation (including unloading) and recovery processes and (2) a much shorter stress rise time at the front in MD ( < 1 ps) in comparison with flyer-plate shock compression (similar to 1ns).

Return to Publications page