Converging shocks in elastic-plastic solids

We present an approximate description of the behavior of an elastic-plastic material processed by a cylindrically or spherically symmetric converging shock, following Whitham's shock dynamics theory. Originally applied with success to various gas dynamics problems, this theory is presently derived for solid media, in both elastic and plastic regimes. The exact solutions of the shock dynamics equations obtained reproduce well the results obtained by high-resolution numerical simulations. The examined constitutive laws share a compressible neo-Hookean structure for the internal energy e = e(s) (I-1) + e(h)(rho, zeta), where e(s) accounts for shear through the first invariant of the Cauchy-Green tensor, and e(h) represents the hydrostatic contribution as a function of the density rho and entropy zeta. In the strong-shock limit, reached as the shock approaches the axis or origin r = 0, we show that compression effects are dominant over shear deformations. For an isothermal constitutive law, i.e., e(h) = e(h) (rho), with a power-law dependence e(h) proportional to rho(alpha) , shock dynamics predicts that for a converging shock located at r = R(t) at time t, the Mach number increases as M proportional to [log(1/R)](proportional to), independently of the space index s, where s = 2 in cylindrical geometry and 3 in spherical geometry. An alternative isothermal constitutive law with p(rho) of the arctanh type, which enforces a finite density in the strong-shock limit, leads to M proportional to R-(s-1) for strong shocks. A nonisothermal constitutive law, whose hydrostatic part e(h) is that of an ideal gas, is also tested, recovering the strong-shock limit M proportional to R-(s-1)/n(gamma) originally derived by Whitham for perfect gases, where gamma is inherently related to the maximum compression ratio that the material can reach, (gamma + 1)/(gamma - 1). From these strong-shock limits, we also estimate analytically the density, radial velocity, pressure, and sound speed immediately behind the shock. While the hydrostatic part of the energy essentially commands the strong-shock behavior, the shear modulus and yield stress modify the compression ratio and velocity of the shock far from the axis or origin. A characterization of the elastic-plastic transition in converging shocks, which involves an elastic precursor and a plastic compression region, is finally exposed.