Abstract
In this study, the impact property of 3D printed 15-5 PH stainless steel was investigated at low (77 K), room (298 K), and high temperatures (723 K) using integrated experimental and modelling studies. The finite element model was based on the Johnson-Cook phenomenological material model and fracture parameters. The experimentally measured impact energies are 0.01, 6.78 ± 4.07, and 50.84 ± 3.39 J cm−2, at the low, room, and high temperatures, respectively. The experimental and modelling predicted impact energies are in good agreement. The microstructures show that the steel exhibits a brittle behaviour at low and room temperatures as indicated by a transgranular fracture, but changes to a more ductile behaviour at high temperatures as illustrated by microvoid coalescence induced facture morphology.
Keywords
Introduction
Temperature-dependent impact properties are important because of the ductile-to-brittle response of metals from high to low temperature as exemplified by the failure of the Liberty ships in the Second World War [1]. The impact performance of the materials is also critical to determining applications in high strain rate deformation applications, e.g. ballistics or crashworthiness [2,3]. To date, there is limited published data on the impact properties of 3D printed 15-5 precipitation-hardening (PH) stainless steel; especially at cryogenic and high temperatures. The study by Pitrmuc [4] and Hendrickson [5] on 3D printed 316L stainless steel sample at room temperature shows the impact energy of 136.0 J cm−2 for horizontal built and 41.2 J cm−2 for a vertical build. Furthermore, at room temperature, Sistiaga [6] has found that the impact energy to be 119.5 J cm−2 using horizontal build using selective laser sintering, whereas the conventionally manufactured equivalent is 103 J cm−2. Therefore, previous research has been focused on the effects of 3D printing technique, process parameters and printing orientation rather than specimen temperature effect on the impact energy of 3D printed stainless steel.
For 15-5 PH stainless steel, the as-received wrought 15-5 PH stainless steel has the impact energy values of 57.2 and 46.7 J cm−2 at the room-temperature and 219 K, respectively. For the wrought steel heat treated at 900 F (755 K), the impact energy values are 50.0 and 41.3 J cm−2 for room-temperature and 219 K, respectively [7]. The DBTT diagrams for 15-5 stainless steel were measured in Ref. [8]. For the unaged steel, the DBTT is around 273 K, and the impact energy values are 25 and 150 J cm−2 at the brittle and ductile regions, respectively.
For the previous impact modelling studies, Amini et al. developed an LS-DYNA model to simulate the response of monolithic DH-36 steel plates and bilayer steel-polyurea plates under impulsive loads [9]. The transient response of the plates under impulsive pressure loads was studied and compared well with the experimentally observed results [9]. Barauskas et al. conducted the computational analysis of the impact of a bullet against the multilayer fabrics in LS-DYNA [10]. To avoid the direct modelling of filaments comprising the yarns, the yarns were modelled by using thin shell elements [10].
This paper presents a combined experimental and finite element modelling study to investigate the temperature-dependent Charpy impact property of 3D printed 15-5 PH stainless steel. The paper is arranged as follows: Section 2 provides the details of the specimen preparation and Charpy impact experiment. Section 3 presents the details of the finite element model for simulating the Charpy impact test. Section 4 gives the experimental and modelling results, and Section 5 summarises the conclusions.
Experimental procedure
Sample preparation
The geometry of the Charpy impact specimens follows the ASTM E23-16b standard [11]. EOS stainless steel PH1 powder (average particle size of 20 μm) was used in this study to produce the 3D printed Charpy test specimens. The specimens were fabricated by the laser powder bed fusion (PBF) process using an EOSINT M270 machine (EOS GmbH Electro Optical Systems, Germany) equipped with a 200 W single-mode Ytterbium fibre laser at 1070 nm wavelength. The processing setting parameters include the laser power of 170 W, scanning speed of 1250 mm s−1 in the continuous wavelength mode, hatch spacing of 100 μm, and layer thickness of 30 μm. The laser beam spot size is ∼50 μm in diameter. Three samples were printed for each of the 77, 298, and 723 K specimens that were prepared similar to the procedure detailed in Ref. [12].
Charpy impact test
Before the Charpy impact test, the samples were placed at room temperature, high-temperature furnace, and liquid nitrogen to achieve the desired testing temperatures. For the low-temperature test, the specimens were subjected to cryogenic temperature using liquid nitrogen for 10 min until the temperature was reduced to 77 K before testing. Similarly, the high-temperature specimens were placed in an oven for 15 min and heated to 1000 K to ensure that the target temperature of 723 K was achieved immediately before testing. The Charpy impact test was undertaken similar to [12].
Numerical modelling
Finite element model
The geometry of the Charpy impact specimens follows the ASTM E23-16b standard. The CAD file in iges format consisting of the striker, specimen, and anvil parts were imported into the structural simulation software LS-PREPOST/LS-DYNA [13] and meshed using solid elements. Further mesh refinement was executed on the specimen at the impact zone as shown in Figure 1.
Finite element model and mesh of striker, specimen, and anvil.
Rigid body properties were assigned on the striker and anvil in order to simulate experimental test conditions and to mitigate against kinetic energy loss through the impacting and support structures. Sliding contact conditions were assigned to the specimen/anvil and striker/specimen interfaces.
Given that the mass of the striker used in the Charpy impact experiment was 30.24 kg, an adjustment of the density based on the revised geometry was necessary. The volume of the striker is 6212.99 mm3; hence the adjusted density is 4.87e−3 kg mm−3.
By equating the potential energy of the dropping striker to the kinetic energy at impact, the initial velocity was calculated using
. The height, h was determined from the experimental Charpy impact test, i.e. 1.366 m, and g is the acceleration due to gravity thus yielding the initial velocity as 5.18 m s−1.
The mesh refinement was undertaken by creating a shell element out of a solid face (side), followed by further refining the mesh within the impact zone. Triangular meshes were created to transition the connection between the larger mesh with the refined zone. A solid mesh was achieved by shell dragging through a compartmentalised extrusion through the thickness of the specimen. The brick solid mesh was used in preference to tetrahedral mesh because of better accuracy per mesh size for high strain rate application. Some tetra mesh types exhibit volumetric locking which results in a more rigid response, further refinement results in no additional increased accuracy, but with the consequence of higher demand for computing resources.
The distance between the striker and specimen was assigned a relatively close distance, i.e. 1.0 mm in order to reduce analysis time. The Flanagan-Belytschko integration hourglass-control algorithm is used to damp out the overall zero-energy modes.
Johnson-Cook constitutive material model
In order to capture the high strain rate response of the specimen, the Johnson-Cook (J-C) phenomenological constitutive model [14] was used. The J-C is a flow stress model that consists of strain hardening, strain-rate, and thermal softening components as shown in Equation (1).
effect stress, A plastic yield stress, B plastic hardening parameter, C strain rate coefficient, n plastic hardening exponent, m temperature softening exponent,
effective plastic strain,
dimensionless plastic strain rate where
and
represents reference and effective the plastic strain rates respectively, and
the dimensionless homologous temperature, in which,
,
and
are the initial, reference, and melting temperatures respectively.
The temperature softening expression shown in Equation (1), represents the effect on flow stress due to the adiabatic condition caused by plastic work of deformation. The reference temperature,
represents the temperature of the specimen before testing, i.e. 77, 298, and 723 K, whereas the initial temperature,
represents the temperature after impact and is due to adiabatic heating. The maximum reference temperature,
was significantly lower than the melting temperature of the specimen,
because of the uncertainty of the accuracy of the J-C constitutive model in predicting the behaviour of the specimen at the near-melting point.
values of 77 K, 298 K and 723 K in the LS-DYNA model to match the behaviour of the experimental Charpy impact test; however, the effect of each constitutive model parameter on the impact energy was found to be negligible. More accurate values of each parameter may be explored using the torsion test, dynamic Hopkinson bar test, and tensile tests over a range of temperatures [14]. The J-C parameters used for defining the material properties of the specimen are presented in Table 1.
Johnson-Cook constitutive model parameters for 15-5 PH stainless steel.
Johnson-Cook fracture parameters
In order to model the fracture and failure characteristics beyond plastic deformation, the Johnson-Cook (J-C) failure criteria [17] were used. A fracture occurs when the following fracture parameter,
reaches unity as shown in Equation (2):
is the increment of equivalent plastic strain, and
is the equivalent strain at failure under given conditions defined by Equation (3):
initial failure strain,
exponential factor,
triaxiality factor,
strain rate factor,
temperature factor,
triaxiality ratio, where
and
represent the hydrostatic pressure and effective stress respectively. The input parameters
are material failure dimensionless parameters that are experimentally determined from the torsion test over a range of strain rates, Hopkinson bar test over a series of temperatures, and quasi-static tensile test over an assortment of notched specimens.
The J-C fracture model parameters,
for the 3D printed 15-5 PH stainless steel were extracted from Ref. [18]. The fracture parameters
and
were derived from [17] given the similar fracture mechanism as AISI 4340 steel.
By keeping the J-C constitutive material parameters and fracture parameters
constant, the numerical iterations indicated that the impact energy was mostly affected by
. The optimised values of
corresponding to
values of 77, 298, and 723 K were −0.600, −0.193, and −0.143 respectively. The damage parameters along with their references used for the Charpy impact test simulation are presented in Table 2.
A linear polynomial equation of state was selected for the Johnson-Cook material model. Also, in order to activate the material model, the bulk modulus, C1 (169.1 GPa) were placed in the input card while C2 ∼ C5 were set to zero.
Physical and mechanical properties
Results and discussion
Microstructures of the printed samples at various temperatures
Before the Charpy impact test, the sample was placed at room temperature, furnace, and liquid nitrogen to achieve the desired testing temperatures. Different temperatures affect the microstructure and the mechanical properties of the stainless steel including its ductility and fracture toughness. Microstructures of the specimen cross-sections were characterised using optical microscopy (OM) and scanning electron microscopy (SEM).
The optical microscopic images of three samples at 77, 298, and 723 K are shown in Figure 2. All of the three samples show dense microstructure with high relative densities (>99%), with very few pores.
Optical micrographs of sample cross-sections of 3D printed 15-5 PH stainless steel at (a) 77 K, (b) 298 K, and (c) 723 K.
As shown in Figure 2, the austenite phase is evident. The average grain size is determined to be ∼5 μm, following ASTM E112-13 standard. The grain size in this study is similar to the reported grain size 6.7∼5.6 μm in Ref. [19].
The surface textures of room temperature and low-temperature samples are very similar, except that the room temperature one has a few visible pores (Figure 2(b)). When heat-treated at high temperatures, large grains are formed (Figure 2(c)), an indication of grain growth at high temperatures. The high-temperature heat treatment acted as an annealing process so that the crystal lattice may relax and more ordered grains are formed.
Detailed microstructures are revealed by the SEM images as shown in Figure 3. At room temperature (Figure 3(b)), irregularly shaped lack-of-fusion pores are observed due to the insufficient heating by laser source during the powder bed fusion process. When the as-printed sample is subjected to furnace heat treatment, the high temperature acts in two ways: increase material diffusion rate and lead to thermal expansion. The surface of the high-temperature sample (Figure 3(c)) is smoother, which indicates that the rough surface is flattened by material diffusion at high temperatures. The black spots in Figure 3(c) are pores, as the result of the residual pores and unmelt particles in the powder bed fusion metal 3D printing process. On the other hand, due to the thermal expansion, the compressed materials near the pore move toward the pore centres, which makes the pore size smaller, and the pore becomes to be circular. At low temperatures, because the diffusion rate is much lower, the structural change is minimum, therefore the configuration is similar to the room temperature.
SEM images of the sample cross-sections at (a) 77 K, (b) 298 K, and (c) 723 K. (Scale bars represent 5 μm).
The microstructural evaluation suggests that the 3D printed samples are martensitic, which is consistent with the images reported in Ref. [20]. After high-temperature homogenisation treatment, the austenite microstructure undergoes a martensitic transformation during the subsequent quench giving rise to the formation of martensite laths inside the former austenite grains [20]. This martensitic structure is present after the precipitation treatment [20], similar to the one shown in Figure 2(c) in this study.
Elements in the 3D printed 15-5 PH stainless steel samples.
Fracture surfaces after Charpy impact test
Microstructure analysis of the fracture surface
After the Charpy impact test, the microstructures of the fracture surfaces were analysed to investigate the fracture mechanisms and the material behaviour under impact loading. The fracture surface micrographs of the low-temperature sample are shown in Figure 4. For the low-temperature specimen shown in Figure 4, a dominant transgranular or cleavage fracture morphology is evident which is usually associated with low-energy, brittle fracture caused by sources of weakness in the material [21]. The transgranular fracture or cleavage occurred along specific crystallographic planes, as shown in Figure 4. The transgranular fracture is commonly observed in certain BCC and HCP metals, but can also occur in FCC metals when they are subjected to severe environmental conditions such as extremely high strain rates or very low temperatures [21], which is the low-temperature case in this study. The transgranular fracture or cleavage is characterised by a relatively flat fracture surface with small converging ridges [21], as shown in Figure 4. Such a low toughness means that temperature 77 K is below the ductile-brittle transition temperature for the material. The movement of nearby atoms in the crystal lattice is suppressed by the low temperature, so that long cracks can form easily, resulting in a brittle fracture behaviour, which is consistent with the recorded impact energy at low temperature is close to 0 J cm−2.
Scanning electron microscopic images of 15-5 PH stainless steel samples test at low temperature (77 K). (2000× magnification. The scale bars represent 10 μm).
Figure 5 shows the fracture surface microstructure of the sample tested at room temperature (298 K). It illustrates a primary transgranular fracture with an increasingly small amount of intergranular fracture compared to the low-temperature specimen was observed. Such a microstructure suggests a brittle fracture, but with improved ductility compared to the low-temperature specimen. At room temperature, the tested impact energy is 6.78 ± 4.07 J cm−2. Although it is still brittle fracture at room temperature, the higher temperature increases the mobility of the lattice, and long crack propagation is partially suppressed. Therefore the impact energy is higher than that of low temperature.
Scanning electron microscopic images of 15-5 PH stainless steel samples test at room temperature (298 K). (2000× magnification. The scale bars represent 10 μm).
In contrast to the low temperature specimen in Figure 4, for the high-temperature specimen shown in Figure 6, a more ductile fracture morphology than the low-temperature specimens was observed. The plastic deformation tends to fail by a process called microvoid coalescence (MVC) [21]. The deformation and fracture took place by the nucleation of microvoids, followed by their growth and eventual coalescence to form cracks. Figure 6 shows the microvoid coalescence under shear loading, which leads to an elongated dimple morphology [21]. The tested impact energy, in this case, is 50.84 ± 3.39 J cm−2, which is the largest for three temperatures. This is because, at higher temperatures, the resistance to dislocation motion is reduced. So more impact energy can be absorbed.
Scanning electron microscopic images of 15-5 PH stainless steel samples test at high temperature (723 K). (2000× magnification. The scale bars represent 10 μm).
Comparison of fracture surface between with modelling results
The fracture surfaces of the broken specimens are compared with the modelling results. A qualitative assessment of both the experimental and numerical modelling indicates that a predominantly brittle failure at 77 K (shown in Figure 7(a,b)). The fracture surfaces are almost flat and perpendicular to the tensile direction, indicating a brittle failure as the striker moves through the specimen. Also, the maximum effective plastic strain of 0.451 m m−1 is the lowest among the three testing temperatures. When tested at room temperature (298 K), the fracture is still brittle (shown in Figure 7(c,d)), but the fracture surface is not as flat as the 77 K specimen, indicating that there has been an increase in ductility, and a small portion of the cracks propagates along the 45° maximum shear stress direction. The maximum effective plastic strain of 0.721 m m−1 is higher than the value at 77 K but lower than that at 723 K. A ductile failure is observed from the specimen tested at 723 K (shown in Figure 7(e,f)) as the fracture surfaces show a significant level of convexity. Most of the cracks propagate along the maximum shear direction. The effective plastic strain has a maximum value of 1.136 m m−1. This represents the highest effective plastic strain, which confirms the highest ductility in three cases.
Fractured specimens from experiment at (a) 77 K, (c) 298 K, and (e) 723 K, and the corresponding simulated effective plastic strain results at (b) 77 K, (d) 298 K, and (f) 723 K.
Comparison of Charpy impact energy between experiment and modelling results
The experimental tested impact energy at 77, 298, and 723 K are 0.01, 6.78 ± 4.07 and 50.84 ± 3.39 J cm−2 respectively. The numerical modelling results shown in Figure 8 indicate the impact energy at 77, 298, and 723 K are 1.05, 10.46 and 47.07 J cm−2 respectively.
Predicted impact energies as a function of striker displacement at 77, 298, and 723 K.
Both the experimental and numerical modelling results are in good agreement, as presented in Figure 9. It indicates a general trend of increasing ductility with the temperature increase. At 77 K, the measured impact energy in the experimental and numerical model were 0.01 and 1.05 J cm−2 respectively, which are relatively close to zero. At room temperature and high temperature, the modelling results are in the same range as that of the experiment data, suggesting the high-fidelity of the finite element model.
Comparison of impact energy between experimental and numerical modelling.
It is noted that due to the limited precision of the used impact tester, it is challenging to get a precise measurement at lower temperatures, which have low impact energy. Also, tester settings such as the friction were not fully included in the FE model, which assumed a frictionless contact surface; hence may predict a larger discrepancy between the FE model and experimentally measured values at lower temperatures.
Our study show that the 3D printed 15-5 stainless steel is generally more brittle than the wrought counterpart, and low impact energy values were acquired both in experiment and modelling (Figure 9). The brittleness of the printed materials may be because a powder is used as the feedstock in the powder bed fusion 3D printing process. The external defects (pores and unmelted particles) act as sites for crack initiation [22], as shown in Figure 3 in this study. Additionally, the powder may be exposed to humidity and oxygen, which cause impurity and defects in the printed materials and contribute to their brittleness [19]. A similar low toughness in the Charpy impact test for Ti-6A1-4V, AISI 316L, and maraging steel 300 samples prepared by the selective laser melting (SLM) process was observed. It was attributed to the presence of defects like pores, pick-up of impurities like oxygen and nitrogen especially for titanium alloys, and the presence of more brittle nonequilibrium phases [23]. In the case of aging reported in Ref. [24], Transmission Electron Microscopy (TEM) analysis was conducted to show the formation of fine spherical Cu precipitates in stainless steel, which solution strengthens but makes the specimen brittle [24]. As a result, relatively lower impact toughness is obtained as compared to the overaged condition where a combined effect of coarsening of Cu precipitates and increased retained austenite makes the specimen ductile [24].
Conclusions
Charpy impact test was used to investigate the impact properties of 3D printed 15-5 PH stainless steel at three different temperatures. A finite element model based on the phenomenological John-Cook theory was developed to simulate the Charpy test process. The results are summarised below.
The experimentally measured impact energies are 0.01, 6.78 ± 4.07, and 50.84 ± 3.39 J cm−2, at the low, room, and high temperatures, respectively. Experimental and numerically modelled impact energy show a good agreement for all three temperatures. The results show a trend of increasing impact energy with increasing temperature, indicating a transition from brittle to ductile failure as the temperature is increased. Fracture surface micrographs show that temperature plays an important role in the impact property of the 3D printed 15-5 PH stainless steel. When tested at low temperature and room temperature, the material showed a brittle behaviour with low fracture toughness, as indicated by a transgranular fracture manner. When the specimen was tested at high temperature, the impact energy was significantly increased, and a primarily brittle but with a slightly ductile fracture behaviour was observed, as illustrated by a microvoid coalescence induced fracture morphology.
Footnotes
Acknowledgements
The authors also thank Professor David R. Johnson in the School of Materials Engineering at Purdue University for assistance with the Charpy impact tests. S. Sagar is supported by the IUPUI University Diversity Fellowship.
Disclosure statement
No potential conflict of interest was reported by the author(s).
