We present a series of plate-impact experiments on 17-4 stainless steel to study the effect of manufacturing method and heat treatment on the Hugoniot elastic limit (HEL), Hugoniot, phase transformation stress, and spallation strength. Two traditional manufacturing methods were considered, wrought processing and casting, as well as two additive manufacturing methods, laser powder-bed fusion (LPBF) and wire-fed electron beam (EBAM). For both LPBF and EBAM 17-4 stainless steel variants, two billets were printed, enabling the application of two unique heat treatments. The HEL stress depended heavily on the thermal history, with the HEL increasing after the formation of Cu-rich precipitates via heat treatment. The Hugoniot response both below and above the phase transition was unaffected by the manufacturing method or heat treatment. The phase transition stress depended heavily on the thermal history, with its variation being attributed to the presence of various microstructural features. This is supported by a marked increase in the phase transition stress after precipitation hardening. These results suggest that the notion of the phase transition stress being dictated by bulk composition is an oversimplification and the stress fields generated by the meso-scale structure are a dominant force. The spallation strength was lower in the cast material compared to all other 17-4 stainless steel variants due to the presence of brittle δ-ferrite inclusions. Additionally, a drop in the tensile strain-rate was observed in the spallation response above the phase transition stress, which was hypothesized to stem from the kinetics of the reversion to the low-pressure phase during spall.
Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of NixAly, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available NixAly and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.
Mesoscale modeling of shock waves in Ni+Al multilayers poses significant challenges that are due, in part, to shock-induced chemical reactions. Current modeling approaches utilize reactive molecular dynamics (MD), but they are limited to resolving domains of only a few hundred nanometers. In contrast, actual multilayer superlattices can be tens of micrometers thick, and they exhibit non-ideal (i.e., wavy) interfaces. The second part of our research builds upon previous work developing physically based, thermodynamically complete equations of state for various Ni and Al intermetallic compositions. Here, we introduce a novel workflow for high-fidelity mesoscale simulations of Ni+Al multilayers using a continuum hydrocode. By increasing the simulation domain size beyond MD limitations (e.g., 2×6μm2) and incorporating explicit interfacial roughness, we investigate the shock response of Ni+Al multilayers at previously unexplored scales. Our experimental design encompasses nine multilayer geometries with varying roughness amplitudes and tilt angles (θ=15°, 30°, and 45°), alongside 19 flyer impact velocities ranging from 0.3 to 3.0 km/s, resulting in a total of 171 high-fidelity simulations. The bulk shock state from inert 2D mesoscale simulations aligns with the law of mixtures, while temperature and pressure fluctuations strongly correlate with multilayer geometry types. A new metric dubbed the “hot spot probability integral” shows a greater dependence on a tilt angle than interfacial roughness.
We measured the austenitic FeCr18Ni12.5 stainless steel Hugoniot as a function of crystallographic direction to approximately 60 GPa. We shock-compressed FeCr18Ni12.5 samples oriented along ⟨100⟩, ⟨110⟩, and ⟨111⟩ to mean stresses ranging 30.5–58.1 GPa via Ta plate impact in a large-bore powder gun and measured the free-surface velocities with laser interferometry. We unambiguously observed the largest post-shock free-surface velocity along ⟨100⟩ in each experiment, which consequently produced the lowest shock velocity along that orientation. However, the propagation of experimental uncertainties through the impedance matching scheme used to compute the shock velocity produced sufficient uncertainty overlap to preclude definitive conclusion of Hugoniot anisotropy.
Both shock and shockless compression experiments were performed on laser powder bed fusion (LPBF) Ti–5Al–5V–5Mo–3Cr (Ti-5553) to peak compressive stresses near 15 GPa. Experiments were performed on the as-built material, containing a purely β (body centered cubic) microstructure, and two differing heat treatments resulting in a dual phase α (hexagonal close packed) and β microstructure. The Hugoniot, Hugoniot elastic limit (HEL), and spallation strength were measured and compared to wrought Ti-6Al-4V (Ti-64). The results indicate the LPBF Ti-5553 Hugoniot response is similar between heat treatments and to Ti-64. The HEL stress observed in the LPBF Ti-5553 was considerably higher than Ti-64, with the as-built, fully β alloy exhibiting the largest values. The spallation strength of the LPBF Ti-5553 was also similar to Ti-64. Clear evidence of initial porosity serving as initiation sites for spallation damage was observed when comparing computed tomography measurements before and after loading. Post-mortem scanning electron microscopy images of the recovered spallation samples showed no evidence of retained phase changes near the spall plane. The spall plane was found to have kinks aligned with the loading direction near areas with large concentrations of twin-like, crystallographic defects in the as-built condition. For the heat-treated samples, the concentrations of twin-like, crystallographic defects were absent, and no preference for failure at the interface between the α and β phases was observed.
Magnetically-driven, shockless-compression experiments were performed to peak stresses approaching 200 GPa on both a direct energy deposition, additively manufactured (AM) and conventionally, wrought-processed 304L stainless steel to compare their thermodynamic and constitutive responses. Velocimetry measurements were used to infer the response of the 304L stainless steel samples during shockless-compression and release from peak stress. A self-consistent, inverse Lagrangian analysis technique was used to determine the isentrope of each sample to peak compression, while a wave profile analysis method was used to estimate the flow strength, shear modulus, and bulk modulus from the unloading from peak stress. The thermodynamic response of both stainless steels were similar and consistent with current equation of state (EOS) formulations up to 200 GPa. The flow strength, shear modulus, and bulk modulus measurements were also similar between stainless steel variants. The flow strength measurements in the AM material showed evidence of more sample-to-sample variability, particularly at pressures above 100 GPa. This was hypothesized to result from the large grain size imparted during printing. The flow strength and shear modulus measurements of both stainless steels deviated from current calibrations of the Steinberg-Guinan-Cochran and Preston-Tonks-Wallace constitutive models. The bulk modulus values extracted for both stainless steels deviated from those predicted by a Vinet EOS fit to the extracted isentropes, particularly at pressures above 100 GPa. This was found to result from the assumption of linearity of the bulk modulus with pressure in the Vinet model. The observed response of the 304L stainless steel shows a quadratic dependence of the bulk modulus with pressure, which is hypothesized to result from the formation of martensite during loading.
We studied the orientation-dependent dynamic properties and plastic deformation mechanisms in single-crystal austenitic FeCr18Ni12.5 stainless steel shocked to peak stresses ranging 3.8-12.0 GPa. The 〈110〉 and 〈111〉 principal Hugoniots generally match that of polycrystalline 304L stainless steel, but the 〈100〉 Hugoniot exhibits lower shock velocities and comprises two distinct regions. The Hugoniot elastic limits (HELs) and spallation strengths range 0.5-1.2 GPa and 2.2-3.8 GPa, respectively, and exhibit trends previously observed in polycrystalline steels. We generally observed the largest HELs along 〈100〉 and the largest spallation strengths along 〈111〉. Microscopy revealed larger twin and martensite fractions in the 〈100〉 samples, correlations in twinning and martensite prevalence with dynamic response, and orientation-dependent differences in dislocation behavior near the spallation plane.
Pressure-shear plate impact experiments were performed to quantify flow strength of wrought, as-built additively manufactured (AM), and heat-treated and recrystallized AM 304 L stainless steel (SS304L) under combined loading. Impact velocities spanned between 0.03 and 0.24 mm/μs, resulting in corresponding pressures of 0.62–5.93 GPa. Flow strength measurements are comparable for the sample variants across the studied loading conditions; however, shear wave structures significantly differ between sample type. Microstructurally aware simulations indicate local strain differences attributed to anisotropic elastic constants of large grains ( ∼ 1 mm) in the as-built and heat-treated AM may impede the ability to uniformly transmit a shear wave.
The mineral CaF2 is the archetype of the α fluorite structure and its high-pressure phase transition to γ cotunnite is an ideal test bed for exploring the effects of kinetics. The inter-disciplinary topic of the kinetics of dynamically driven phase transitions is at the forefront of condensed matter physics, both for its theoretical importance and its relevance to technological applications at extreme conditions of pressure and temperature. Here we probe the α → γ → α structural transformations taking place over the nanosecond timescale of a dynamic event, beginning-to-end: from the principal shock Hugoniot state, followed by a quasi-steady off-Hugoniot release state, and finally the unsteady return to near-ambient conditions. We present quantitative, atomic-scale data of the unfolding of the dynamically driven phase transition and its subsequent reversal close to the α/γ phase boundary. Dynamic loading with a two-stage gas gun is coupled with in situ time-resolved synchrotron X-ray diffraction and with continuum scale velocimetry at the Dynamic Compression Sector (DCS), Advanced Photon Source, Argonne National Laboratory. Our results demonstrate the time dependence of phase transitions and highlight the need for modeling of transition kinetics in dynamically driven processes.
Uniaxial strain, reverse-ballistic impact experiments were performed on wrought 17-4 PH H1025 stainless steel, and the resulting Hugoniot was determined to a peak stress of 25 GPa through impedance matching to known standard materials. The measured Hugoniot showed evidence of a solid–solid phase transition, consistent with other martensitic Fe-alloys. The phase transition stress in the wrought 17-4 PH H1025 stainless steel was measured in a uniaxial strain, forward-ballistic impact experiment to be 11.4 GPa. Linear fits to the Hugoniot for both the low and high pressure phase are presented with corresponding uncertainty. The low pressure martensitic phase exhibits a shock velocity that is weakly dependent on the particle velocity, consistent with other martensitic Fe-alloys.