Modelling of shock waves in fcc and bcc metals using a combined continuum and dislocation kinetic approach

Abstract

Recent experimental data has revealed that, over short time scales (on the nanosecond scale), during formation of a shock in metals, amplitude of the ‘elastic’ precursor greatly exceeds the Hugoniot elastic limit (HEL), before decaying to the level of the HEL. Standard continuum scale material models are unable to reproduce this behaviour. To capture this aspect of material behaviour in metals, physical effects related to high rate dislocation mechanics must be taken into consideration (Mayer et al., 2013) and included into the continuum scale material model. The constitutive model developed here is defined at the continuum level, where the evolution of plastic deformation is controlled with a system of equations for three microscale state variables, for each slip system of a single crystal. These three state variables are the density of mobile dislocations, the density of immobile dislocations and the mobile dislocation velocity. The density evolutions of mobile and immobile dislocations are controlled by dislocation kinetic equations, which account for the generation of new dislocations, immobilisation of mobile dislocations and annihilation of dislocations. Dislocation velocity is determined by integration of the equations of motion of the mobile dislocations. The dislocation micromechanics is incorporated into the continuum model using the generalised Orowan equation, which relates plastic strain rate to the density of mobile dislocations and the velocity of mobile dislocations. Evolution of the yield surface is controlled by density of immobile dislocations.