Numerical Integration of Plane Elastic-Relaxing Plastic Shock Waves by a Two-Step Method

Hyperbolic equations with a differential constitutive relation of dynamic yielding are numerically integrated without the pseudoviscosity. The scheme is shown to be stable and the resulting shock profiles are smooth. The example of iron shows a good agreement with those of the viscosity method and suggests that a stable profile will emerge as mesh sizes decrease.