Atomistic Modeling of Spall Response in a Single Crystal Aluminum

Materials used in soldier protective structures, such as armor, vehicles and civil infrastructures, are being improved for performance in extreme dynamic environments. Accordingly, atomistic molecular dynamics simulations were performed to study the spall response in a single crystal aluminum atom system. A planar 9.6 picoseconds (ps) shock pulse was generated through impacts with a shock piston at velocities ranging from 0.6 km/s to 1.5 km/s in three <1,0,0>, <1,1,0>, and <1,1,1> crystal orientations. In addition to characterizing the transient spall region width and duration, the spall response was characterized interms of the traditional axial stress vs. axial strain response for gaining an understanding of the material failure in spall. Using an atom section averaging process, the snapshots, or the time history plots of the stress and strain axial distributions in the shock direction, were obtained from the MD simulations’ outputs of the atom level stresses and displacements. These snapshots guided the analyses to an estimation of the spall widths and spall transients. The results were interpreted to highlight the effects of crystal orientation and impact velocity on the spall width, spall duration, spall stress, strain rate, critical strain values at the void nucleation, and the void volume fraction at the void coalescence. For all the combinations of the crystal orientations and the impact velocities, the void nucleation was observed when the stress reached a peak hydrostatic state and the stress triaxiality factor reached a minimum, i.e. by the simultaneous occurring of these three conditions for the stress state: 1. pressure reaching a negative minimum, 2. axial stress reaching the magnitude value of the peak pressure, and 3. the effective stress reaching a zero value. At these conditions, void nucleation was mainly caused by atom de-bonding. In fact, the void nucleation strains were shown to have been preceded by the strains of the stress triaxiality condition in this study, thus confirming the stress triaxiality condition for the void nucleation in spall. Based on the observation that the axial stress reached a maximum value of ∼6 GPa during the void nucleation phase in spall and stayed approximately at that value for different crystal orientations and impact velocities, the value was proposed as a material spall strength.