Shock compression modeling of metallic single crystals: comparison of finite difference, steady wave, and analytical solutions
- US Army, Aberdeen Proving Grounds, MD (United States). Impact Physics Branch; Georgia Inst. of Technology, Atlanta, GA (United States). Woodruff School of Mechanical Engineering
- US Army, Aberdeen Proving Grounds, MD (United States). Impact Physics Branch
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States). Materials Modeling and Simulation Group
- Georgia Inst. of Technology, Atlanta, GA (United States). Woodruff School of Mechanical Engineering and School of Materials Science
Background: The shock response of metallic single crystals can be captured using a micro-mechanical description of the thermoelastic-viscoplastic material response; however, using a such a description within the context of traditional numerical methods may introduce a physical artifacts. Advantages and disadvantages of complex material descriptions, in particular the viscoplastic response, must be framed within approximations introduced by numerical methods. Methods: Three methods of modeling the shock response of metallic single crystals are summarized: finite difference simulations, steady wave simulations, and algebraic solutions of the Rankine-Hugoniot jump conditions. For the former two numerical techniques, a dislocation density based framework describes the rate- and temperature-dependent shear strength on each slip system. For the latter analytical technique, a simple (two-parameter) rate- and temperature-independent linear hardening description is necessarily invoked to enable simultaneous solution of the governing equations. For all models, the same nonlinear thermoelastic energy potential incorporating elastic constants of up to order 3 is applied. Results: Solutions are compared for plate impact of highly symmetric orientations (all three methods) and low symmetry orientations (numerical methods only) of aluminum single crystals shocked to 5 GPa (weak shock regime) and 25 GPa (overdriven regime). Conclusions: For weak shocks, results of the two numerical methods are very similar, regardless of crystallographic orientation. For strong shocks, artificial viscosity affects the finite difference solution, and effects of transverse waves for the lower symmetry orientations not captured by the steady wave method become important. The analytical solution, which can only be applied to highly symmetric orientations, provides reasonable accuracy with regards to prediction of most variables in the final shocked state but, by construction, does not provide insight into the shock structure afforded by the numerical methods.
- Research Organization:
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- Sponsoring Organization:
- USDOE
- Grant/Contract Number:
- AC52-07NA27344
- OSTI ID:
- 1249149
- Report Number(s):
- LLNL-JRNL-663743
- Journal Information:
- Advanced Modeling and Simulation in Engineering Sciences, Vol. 2, Issue 1; ISSN 2213-7467
- Publisher:
- SpringerCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Dynamic Strength of AZ31B-4E and AMX602 Magnesium Alloys Under Shock Loading
|
journal | January 2020 |
Similar Records
Shock-induced chemical decomposition and overdriven detonation in hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) single crystals
Shock wave compression of hexagonal-close-packed metal single crystals: Time-dependent, anisotropic elastic-plastic response of beryllium