Numerical Study on Shock Wave Profiles in 6061-T6 Aluminum Using Shock Viscous Stresses

In order to accurately predict particle velocity profiles of steady shock wave fronts propagating in solid materials, a new numerical calculation method was proposed. The present method was based on one-dimensional Lagrangian finite difference wave code, and shock viscous stresses evaluated theoretically were introduced into the present numerical code. The shock viscous stress that is one of the important parameters to shape the rising profile of the shock wave front was calculated by the inside temperature estimate method for the steady shock wave fronts. The present calculation method was applied to the uniaxial strain problem of 6061-T6 aluminum at some stress levels below 8.86 GPa, where the elastic and plastic steady shock waves will appear. The results of the present calculation could reproduce more accurately the experimental data measured by the velocity interferometer system (VISAR). In addition, the present calculation was applied to a shock stress level of 20 GPa to investigate the shock viscous effects in the overdriven shock wave rising profile.