Modeling of shock compression of isotropic, polycrystalline elastic-plastic solids that deform by dislocation glide is undertaken. Internal energy density of the material depends on a thermoelastic Eulerian strain tensor, entropy, and an internal state variable representative of dislocation density. Dislocation glide is incompressible, but inelastic volume changes arising from residual local strain fields and core effects of dislocations are captured. A semi-analytical method is advanced for extracting inelastic constitutive response information from particle velocity histories of polycrystalline samples under planar shock loading. The only parameters entering the procedure are fundamental thermoelastic properties and assumed bounds on the fraction of plastic work corresponding to energy storage of generated dislocations in the lattice. Densities of statistically stored and geometrically necessary dislocations, in addition to shear stress, plastic strain, plastic strain rate, and temperature, are outcomes of the analysis. The model is implemented for polycrystalline aluminum and copper. Certain results are compared with others in the literature obtained under different kinematic and thermodynamic assumptions. Effects of choices of nonlinear thermoelastic formulations and higher-order elastic constants on extracted deviatoric stress profiles can be substantial. Stored energy of dislocations and commensurate residual volume changes together can notably affect extracted shear stress. Stored energy significantly affects predicted entropy production and its contribution to temperature rise. The calculated density of geometrically necessary dislocations required to maintain compatibility of total strain across the structured steady plastic wave front is small relative to the line density of total dislocations.
Read full abstract