**The Influence of Vanadium Additions on Isothermally Formed Bainite Microstructures in Medium Carbon Steels Containing Retained Austenite**

#### **Irina Pushkareva 1,\*, Babak Shalchi-Amirkhiz 1, Sébastien Yves Pierre Allain 2, Guillaume Geandier 2, Fateh Fazeli 1, Matthew Sztanko <sup>1</sup> and Colin Scott <sup>1</sup>**


Received: 25 February 2020; Accepted: 16 March 2020; Published: 19 March 2020

**Abstract:** The influence of V additions on isothermally formed bainite in medium carbon steels containing retained austenite has been investigated using in-situ high energy X-ray diffraction (HEXRD) and ex-situ electron energy loss spectroscopy (EELS) and energy dispersive X-ray analysis (EDX) techniques in the transmission electron microscope (TEM). No significant impact of V in solid solution on the bainite transformation rate, final phase fractions or on the width of bainite laths was seen for transformations in the range 375–430 ◦C. No strong influence on the dislocation density could be detected, although quantitative analysis was impeded by ferrite tetragonality. A reduction in the carbon content of retained austenite *C*γ that is not believed to be due to competition with VC or cementite precipitation was observed. No influence of V on the carbon supersaturation in bainitic ferrite Cb could be directly measured, although carbon mass balance calculations suggest Cb slightly increases. A beneficial refinement of blocky MA and a corresponding size effect induced enhancement in austenite stability were found at the lowest transformation temperature. Overall, V additions result in a slight increase in strength levels.

**Keywords:** bainite; vanadium microalloying; austenite stability; HEXRD; EELS

#### **1. Introduction**

Vanadium additions are known to be beneficial for improving the strength and/or the in-use performance of several classes of advanced high strength steel (AHSS) products. These include low carbon hot rolled bainitic steels [1], medium carbon air-cooled bainitic forging steels [2–4], intercritically annealed low and medium carbon transformation induced plasticity (TRIP) and dual phase (DP) sheets [5], medium carbon ferritic hot strips [6,7] and high carbon austenitic twinning induced plasticity (TWIP) steels [8,9]. Classically, the vanadium strengthening effect in carbon-manganese steels is attributed to a combination of microstructure refinement and precipitation hardening [10]. However, vanadium can also provide positive effects while still in solid solution. For example, in bainitic steels that are fast cooled from temperatures above Ac3 and transformed below Bs almost all of the added vanadium remains in solid solution [3,4]. Nevertheless, for low carbon alloys it is clear that there can be a beneficial strengthening effect from vanadium due to a notable decrease in the bainite transformation temperature [1]. Further, bainite softening during coiling is reduced in a manner analogous to the well known effect of vanadium on the temper softening of martensite [11]. While the situation for low carbon (i.e., containing < 0.1 wt. % C) bainite now seems reasonably clear, few detailed studies on

technologically important medium carbon bainitic alloys (0.2–0.4 wt. % C) are available. For example, in 0.25 wt. % C alloys with V additions of up to 0.4 wt. % Sourmail and co-authors [2,3] reported no significant precipitation or clustering of V and no effect on the bainite reaction for isothermal transformations in the range 375–450 ◦C. In fact, the dominant effect of vanadium in these alloys was a strong reduction in the critical cooling rate for pearlite formation, leading to a significant improvement in hardenability for large sections. Wang et al. [4] found that the addition of 0.13 wt. % V in bainitic forging steels was advantageous for both strength and toughness under as-forged slow air-cooled conditions, but that the benefits were reduced as the cooling rate increased.

In low carbon steels the bainitic transformation is very rapid and, depending on the composition and cooling rates applied on the hot strip mill run out table, can begin during fast cooling before the coiling temperature is reached. The situation is more complicated in the medium carbon steels considered here as the bainitic transformation is much slower and carbon partitioning will result in the presence of a second phase of carbon-stabilized retained austenite (RA). Often, silicon additions are employed to delay cementite formation in bainitic ferrite and increase the amount of carbon available for retained austenite stabilization [12]. Further, the high dislocation densities formed in fresh bainitic ferrite may accelerate vanadium carbide and nitride nucleation. In this work, a powerful combination of in-situ high energy synchrotron X-ray diffraction (HEXRD) and ex-situ electron back scattered diffraction (EBSD) and electron energy loss spectroscopy (EELS) and energy dispersive X-ray mapping (EDX) techniques to determine the effect of vanadium additions on the evolution of the microstructure before, during and after the bainitic transformation in two model alloys is applied. An important aspect of this study is that all of the experimental work including the heat treatment, in-situ diffraction, ex-situ microscopy and mechanical testing was carried out on the same samples. This was done in order to eliminate, as far as possible, the effects of heterogeneous composition and/or processing.

#### **2. Materials and Methods**

In this work the effect of vanadium on the bainitic transformation in two laboratory alloys, Ref and Ref+V was studied, using a range of in-situ diffraction and ex-situ microscopy techniques. The in-situ work was done using synchrotron X-ray diffraction to determine the evolution of the phase fractions, the average carbon content in austenite and the dislocation density in bainitic ferrite. Ex-situ optical and scanning electron microscopy (including EBSD) were applied to compare the microstructural parameters such as bainite lath size and MA island sizes. Transmission electron microscopy was also employed to characterize precipitate size and chemistry using EDX mapping and to measure the local carbon content in austenite, martensite and bainitic ferrite.

The chemical compositions of the two alloys, Ref and Ref+V, are given in Table 1. The nitrogen content was deliberately kept low in order to maximize the amount of vanadium in solution during the bainite transformation. The vanadium content was based on previous work on TRIP alloys [13] and is typical of bainitic forging applications. The steels were cast into 50 kg ingots from a single 200 kg heat in a vacuum induction melting (VIM) furnace at CanmetMATERIALS (CMAT). The chemical compositions of the solidified ingots were verified by optical emission spectroscopy (OES) and LECO combustion analysis for C and N. Hot rolling was carried out on the CMAT pilot hot strip mill. The ingots were reheated to 1220 ◦C, held at temperature for 3 h and then hot rolled from 75 to 15 mm in 6 passes. The final rolling pass was carried out above 900 ◦C and the finished plates were allowed to air cool.


**Table 1.** Chemical compositions of the studied steels, wt. %.

Hollow tube specimens with an outer diameter of 4 mm, an inner diameter of 2 mm and a length of 10 mm were directly machined from the hot rolled plates. They were cut with the cylindrical axis parallel to the rolling direction and centered at <sup>1</sup> <sup>4</sup> plate thickness. The tube geometry has two important advantages for this study; firstly the amount of X-ray absorption is reduced by having thin wall sections and secondly induction heating and gas cooling is much more efficient and homogeneous with hollow samples. A Bahr DIL 805 deformation dilatometer was used to determine the alloy transformation temperatures Ac1, and Ac3 using a heating rate of 5 ◦C/s. Optical microscopy (OM), EBSD and Vickers microhardness testing were all done on wall sections obtained by cutting the tubes into two equal half cylinders and polishing the cut faces. Specimens for OM were first polished with 1 μm diamond paste, then finished with Struers OPS solution before etching with a Klemm reagent. EBSD maps were acquired from the same OM samples after repolishing. The EBSD data was acquired using an EDAX TSL system on a Nova NanoSEM 650 field emission gun scanning electron microscope (FEG-SEM). The spot spacings used were 200 nm and 80 nm. Vickers microhardness was measured using a Clemex MMT-X7B tester under 300 gf force and 10 s dwell time according to the ASTM E384 standard.

The in-situ HEXRD experiments were carried out at the Petra (Hamburg) P07 beam line operated in transmission mode under powder diffraction conditions. A high-energy monochromatic beam with 100 keV energy and 0.12 Å wavelength permitted data acquisition in the transmission mode at a frequency of 10 Hz. A 2-D Perkin-Elmer charge coupled device (CCD) detector was positioned 1 m behind the sample, giving access to full Debye-Scherrer rings with a maximum 2θ angle of 12◦. The X-ray beam cross section was 400 μm × 400 μm so that the analyzed volume contained at least 10,000 austenitic grains at the start of each experiment. The applied thermal cycles were controlled using another Bahr dilatometer available at the P07 beam line. Specimens were heated by induction and cooled by argon gas. The temperature was regulated and recorded by a thermocouple spot-welded at the center of the sample. The setup was the same as that used in [14].

The 2-D diffraction patterns produced during the synchrotron experiments were circularly integrated using the Fit2D software [15] after calibration using Si powder and the resulting 1-D diffractograms were analyzed with a full Rietveld refinement procedure. The diffraction peaks were fitted by pseudo-Voigt functions using the FullProf code [16] with 16 degrees of freedom (background, lattice parameters, peak shape and temperature effects) to determine the phase fractions and the lattice parameters. In all the following experiments only two phases were detected in the diffraction patterns: a face-centered cubic phase, corresponding to austenite; and a body-centered phase. The latter could correspond to either body-centered tetragonal (bct-martensite or supersaturated bct-bainitic ferrite) or body-centered cubic (bcc) ferrite. It was not possible to isolate the contribution of each based on tetragonal distortion. Therefore, during the Rietveld refinement procedures the lattices of the two phases were considered to be cubic (Fm3m for austenite and Im3m for martensite/bainite). It is also important to note that no peaks from carbides of any type (including cementite) were detected on the diffractograms.

The austenite carbon content *C*γ starts to change from the nominal value during the bainitic transformation when carbon diffuses from fresh bainite into the parent austenite phase. This causes the austenite lattice parameter *a*γ to increase linearly with the amount of carbon rejected from bainitic ferrite. In these trials the transformation is isothermal, so the austenite carbon content was determined by calculating the relative lattice dilatation Δ*a*<sup>γ</sup> between each measured lattice parameter value, *am* and a reference lattice parameter determined at the transformation temperature before any bainite has formed, *atr*. This difference is proportional to the austenite carbon enrichment Δ*C*<sup>γ</sup> as follows [14]:

$$
\Delta a\_{\mathcal{V}} = a\_{\rm m} - a\_{\rm tr} = k\_{\mathbb{C}} \Delta \mathcal{C}\_{\mathcal{V}} \tag{1}
$$

The coefficient of proportionality *kC* is taken to be 0.00105 nm/(at. %)C [17,18]. The absolute carbon content in austenite is then obtained by adding Δ*C*<sup>γ</sup> to the nominal carbon content.

Planar TEM thin foils were prepared from the same specimens that were analyzed in the synchrotron. To do this, the tube samples were cut in half in the direction parallel to the rotation axis and one of the resulting half cylinders was gently flattened and then mechanically thinned to 100 μm. The other half cylinder was used for microscopy and microhardness testing. Standard 3 mm disks were punched out of the flattened tubes and perforated by twin-jet polishing using a Struers Tenupol 5 and Struers A8 electrolyte. The specimens were analyzed in a Technai Osiris 200 keV X-FEG (S)TEM equipped with a 4 element windowless Super-X EDX detection system and a Gatan Enfina EELS analysis system. The carbon content in martensite, bainite and retained austenite phases was measured by EELS using a standards based second difference acquisition technique developed by the authors [18]. All equilibrium and kinetic thermodynamic calculations were carried out using the MatCalc 6.02 (rel 1.003) software [19] using the mc\_fe\_v2.011 diffusion data base, the mc\_fe\_1.003 physical database and the mc\_fe\_2.059 thermodynamic database.

#### **3. Results**

Section 3.1 is concerned with the determination of the optimum experimental conditions. Section 3.2 contains the OM, SEM and microhardness data. Section 3.3 provides an exploitation of the in-situ synchrotron results. Complementary data obtained from TEM analyses are presented in Section 3.4.

#### *3.1. Definition of the Thermal Cycles*

The solution treatment employed here was designed as a compromise to best satisfy several mutually incompatible requirements:


The equilibrium transformation temperatures Ae1 and Ae3 for Ref and Ref+V were calculated using the MatCalc 6 code and compared with experimental dilatometry data [20] in Table 2. The bainitic start temperature Bs and the martensite start temperature Ms calculated from the empirical equations of Kircaldy and Venugopalan [21] and Andrews [22] respectively are also shown.


**Table 2.** Experimental and calculated transformation temperatures.

As expected, the addition of vanadium tends to increase the equilibrium transformation temperatures and stabilize the ferrite phase. The most important effect is on the kinetics—during continuous heating vanadium greatly reduces the austenite nucleation and growth rate in the intercritical region, resulting in a >40 ◦C increase in Ac1 and Ac3 at a constant heating rate of 5 ◦C/s. This behavior has been previously observed by the authors in microalloyed dual phase (DP) steels [5]. The calculated solution temperature for V(C,N) in Ref+V is 1067 ◦C. In the interests of optimizing the austenite grain size for the synchrotron experiment it was decided to use the lowest possible austenitizing temperature and soaking time and tolerate the inevitable presence of a small fraction of precipitated V(C,N). This will be considered further in the discussion section. Therefore a combination of 900 ◦C/120 s for the austenitization parameters was chosen. It was verified by dilatometry that this was sufficient for the complete austenitization of both alloys. Not much reliable information on the effect of vanadium on Bs and Ms exists for this class of steels. Therefore, three isothermal bainitic transformation temperatures spanning the range of the theoretical values; 375 ◦C, 400 ◦C and 430 ◦C were studied. The optimum cooling rate after austenitization was determined experimentally from the initial synchrotron trials. The cooling rate should be fast enough to avoid the formation of bainite before the isothermal temperature is attained, but slow enough so that there is no significant temperature undershoot at the transition point. It was found that a cooling rate of 40–50◦ C\*s−<sup>1</sup> satisfied these requirements. The experimental thermal cycles are shown in Figure 4a.

#### *3.2. Microstructure and Mechanical Properties*

The microstructure of the Ref and Ref+*V* alloys transformed at 375 ◦C can be seen after Klemm etching in Figure 1. The matrices (colored dark) are fully bainitic with second phase retained austenite and/or martensite (MA) colored light. The second phase consists of two distinct morphologies—alignments of thin films or small islands along interlath boundaries and prior austenite grain (PAG) boundaries and coarser blocky islands with a more random distribution. The area fraction and the mean size of second phase islands were determined using the ImageJ analysis software for all three transformation temperatures. The results are listed together with the corresponding synchrotron HEXRD, EBSD and TEM data in Tables 3 and 4.

**Figure 1.** Optical Micrographs (Klemm etch) of Ref and Ref+V alloys isothermally transformed at 375 ◦C in the synchrotron (**a**) Ref and (**b**) Ref+V. The light colored phase is martensite/retained austenite.

The EBSD data from specimens transformed at 375 ◦C are presented in Figure 2a–f. Similar results were obtained for the 400 ◦C transformation (not shown). The 430 ◦C transformation samples were not studied. The Euler orientation maps are shown in Figure 2a,b, the size distributions determined from low angle boundaries (θ < 5◦) are compared in Figure 2c, the grain misorientation frequencies in Figure 2d and the image quality (IQ) maps with superimposed austenite islands (red color) in Figure 2e,f. The frequency distributions of grain boundary misorientation angles were typical of lower bainitic microstructures with dominant fractions of low misorientation (θ < 15◦) and high misorientation (θ > 50◦) boundaries. Vanadium additions clearly had little or no effect on the distributions at 375 ◦C. The effective grain size for low angle boundaries (θ < 5◦) thought to define the effective lath size [23] is reported in Table 4. An estimate of the MA fractions and mean MA island area was determined from the poorly indexed regions in the IQ maps. This information, together with the RA fraction is reported in Tables 3 and 4.

**Figure 2.** Electron back scattered diffraction (EBSD) data from Ref and Ref+V alloys transformed at 375 ◦C: (**a,b**) Euler orientation maps with high angle (θ > 15◦) boundaries shown; (**c**) low angle (θ < 5◦) boundary size distributions; (**d**) misorientation frequencies and (**e,f**) image quality maps with superimposed austenite islands in red.


**Table 3.** Phase fractions from synchrotron high energy synchrotron X-ray diffraction (HEXRD), EBSD and optical microscopy data.

**Table 4.** Microstructure parameters from EBSD, TEM and OM data.


Vickers microhardness measurements were obtained from the as-rolled plates and from the tube samples after thermal treatment in the synchrotron (Table 4, Figure 3). The hardness of the air-cooled hot strips was found to be 343 ± 15 HV for the Ref alloy and 361 ± 12 HV for Ref+V after air cooling. This is equivalent to an increase of ~+60 MPa in the YS [24]. As expected, fast cooled and isothermally transformed bainitic samples were all stronger than the as-rolled condition, with the highest hardness (422 HV) occurring for the Ref+V alloy transformed at 375 ◦C. Vanadium additions improved the alloy strength at all transformation temperatures, but the difference at 400 ◦C was very small. This could indicate that the presence of vanadium introduces different strengthening mechanisms operating at low and high transformation temperatures. The maximum vanadium-related hardness increase after isothermal transformation occurred at 430 ◦C (+15 HV or ~+50 MPa). The gain in strength was therefore proportionately much lower than that achieved under slow air cooling conditions.

**Figure 3.** Vickers microhardness data from air-cooled as-rolled and isothermally transformed samples.

#### *3.3. Synchrotron XRD—Phase Fractions, Carbon Partitioning and Dislocation Densities.*

The evolution of the bcc phase fraction with time and temperature can be readily determined from the in-situ XRD data. Figure 4a shows the results obtained for the Ref and Ref+V alloys isothermally transformed at 375 ◦C, 400 ◦C and 430 ◦C. Time t = 0 s on the horizontal axis corresponds to the start of the fast cooling ramp after soaking at 900 ◦C. The isothermal section ends at t = 612 s. At this point, any unstable austenite transforms to martensite during the final cooling step. The amount of martensite formed during secondary cooling is very limited except for the Ref+V alloy transformed at 430 ◦C. Although the HEXRD detection limit for ferrite is excellent (0.1% volume fraction) there is a large mathematical uncertainty in the calculated volume fractions for ferrite fractions < 5% during the refinement procedure. This is partly due to the highly textured nature of the very first ferrite nuclei. In consequence only data with >5% bcc fraction is reported. The phase fractions are summarized in Table 3 (uncertainty ± 1%).

**Figure 4.** HEXRD data from Ref and Ref+V alloys as a function of the transformation temperature. (**a**) Body-centered cubic (bcc) phase fraction and (**b**) carbon content in austenite.

Following current theory, the total amount of bainite formed increased as the isothermal transformation temperature decreased (Figure 4a) and the RA fraction followed the opposite trend (Table 3). From the HEXRD data, the addition of vanadium had virtually no effect on the final bainite fraction after transformation at 375 ◦C and 400 ◦C and introduced a slight decrease at 430 ◦C. It had no influence on the final fraction of retained austenite at any transformation temperature, although it did increase the unstable austenite content at 430 ◦C. This slight excess austenite transformed to give 5.9% martensite fraction during final cooling, compared to only 2.5% for the Ref alloy under identical conditions. The bainite start transformation temperature, *Bs*, defined here as the temperature at which 5% of the final bainite fraction is detected, occurred either at or very close to the isothermal temperature for all of the trials, except for Ref+V at 400 ◦C. The time required to complete 95% of the bainite transformation decreased as the transformation temperature is lowered—this is expected as the thermodynamic driving force for austenite decomposition increases. The addition of vanadium appears to very slightly accelerate the transformation, but the effect is not significant.

The magnitude of the RA fraction detected by EBSD was lower than the HEXRD data by a factor of >2 and the MA fraction was much higher. Other authors have reported analogous discrepancies with similar bainitic steels [4]. Either more than 50% of all RA is contained in films and islands that are below the EBSD resolution limit (80 nm step size) or it is unstable and transforms to martensite during the specimen preparation. Nevertheless, all techniques show that the same trends are followed. For both steels, the mean MA island size decreased as the transformation temperature was lowered. However, a significant additional refinement from the vanadium addition was apparent at 375 ◦C where the mean EBSD MA island size in Ref+V was 0.38 μm2 compared to 0.53 μm<sup>2</sup> for the Ref alloy.

The evolution of *C*γ as the bainite transformation progresses can be seen in Figure 4b. Comparing Figure 4b with Figure 4a, it is interesting to note that partitioning of carbon to austenite follows a sigmoidal type kinetic that is much slower than the bainite formation rate. Indeed, at the lowest transformation temperature partitioning is still active at the end of the isothermal plateau (600 s), even although the bainite transformation stopped more than 300 s before. The final RA carbon content in both alloys (Table 5) increased monotonically with decreasing transformation temperature in agreement with the T0 model [25].

From Figure 4b and Table 5, the addition of vanadium changed the austenite carbon enrichment, i.e., Cγ was systematically lower in Ref+V at the end of transformation, for all transformation temperatures. The difference was large at 430 ◦C and at 375 ◦C but very small at 400 ◦C. As all of the Ref+V samples were machined from the same plate and processed identically, there is unlikely to be any variation in the fraction of carbon trapped in vanadium carbide precipitates formed prior to the fast cooling step. Now, for a fixed transformation temperature Table 3 shows that the RA fraction is the same for both alloys, so these data suggest that vanadium is in some way restricting the amount of carbon available for partitioning to austenite. Further, the severity of the effect may have a complex temperature dependence. This could be related to carbon trapping by undissolved prior VC precipitates or the formation of new VC/VxCy clusters in bainitic ferrite. Alternatively, vanadium may promote faster cementite formation in bainite. It is also possible that vanadium could directly consume some of the carbon partitioned to austenite by precipitating as VC inside highly enriched austenite islands. Another theory is that solute vanadium retards recovery in bainitic ferrite leading to a higher defect (dislocation) density that in turn increases the effective carbon solubility. These points are discussed further in the next sections.

#### Measuring Dislocation Densities

The dislocation density in bainitic ferrite at the end of the isothermal transformation is a critical parameter in evaluating the structure/properties relationship of these steels. TEM imaging can be used to estimate dislocation densities on a microscopic scale (as long as the foil thickness is known), however the highly disorganized nature of these microstructures (see the next section) would require extensive sampling to avoid serious errors in determining the average dislocation density. HEXRD can provide an accurate measurement of global dislocation densities [26]. Normally, Williamson–Hall (W–H) or Warren–Averbach (W–A) type analyses are used to isolate the contribution to the full width half maximum (FWHM) peak broadening from microstrain and hence determine the dislocation density. An attempt was made to calculate the dislocation densities in bainitic ferrite from the {110}α, {200} α, {211} α, {220} α and {310} α peak line broadening (after subtraction of the contribution from instrumental broadening) using the Williamson–Hall (W–H) method [27,28]. Unfortunately, the W–H values for dislocation densities were unphysically high for both alloys, almost certainly because tetragonality-induced peak splitting was present. However, the W–H plots did show a strong linear dependence on sin θ, indicating that strain broadening and not crystallite size was the major constituent. The degree of influence of vanadium on the dislocation density in bainitic ferrite was estimated by comparing the angular width at half maximum of the {011}α peaks from both alloys at the end of the isothermal hold (Figure 5a,b). The {222}α reflection, whose angular width is less sensitive to tetragonality than {011}α, would be a better choice for this. However, as there is almost perfect superposition between {222}α and {311}γ reflections, these peaks cannot be separated. Based on the trends seen in the W–H plots, it is reasonable to interpret Figure 5b as a qualitative plot of the variation of dislocation density with transformation temperature. For both steels, the dislocation density was highest during transformation at 375 ◦C and linearly decreased as the temperature increases. From Figure 5b the absolute difference in the {011}α FWHM of Ref and Ref+V was never more than 1.7% at any transformation temperature (measured after 600 s). As this is comparable with the uncertainty in the FWHM measurement, it is not possible to demonstrate that there is any significant difference in the dislocation density of bainitic ferrite due to vanadium additions.

**Figure 5.** (**a**) HEXRD {011}α reflection for Ref and Ref+V transformed for 600 s at 375 ◦C. (**b**) Comparison of {011}α peak broadening for Ref and Ref+V transformed for 600 s at 375 ◦C, 400 ◦C and 430 ◦C.

#### *3.4. TEM—Microstructure*

Thin foil samples from all of the specimens listed in Table 3 were examined in the TEM. Consistent with Figure 1, the matrix phase was lath bainite (Figure 6a) with thin (<100 nm) films of retained austenite situated between the bainite laths (Figure 6b). Some martensite was detected at higher transformation temperatures, usually associated with the RA films (Figure 6c). The linear defects visible in the austenite lath of Figure 6c are most likely epsilon martensite. Other areas showed plate or granular bainite (Figure 7). Overall, very few cementite precipitates were observed, in agreement with the HEXRD findings.

**Figure 6.** (**a**) Bright field STEM image of the Ref alloy transformed at 430 ◦C; (**b**) dark field STEM image showing austenite films (γ(110) reflection) from Ref+V transformed at 430 ◦C and (**c**) bright field STEM image of Ref transformed at 430 ◦C showing RA partially transformed to martensite.

The bainitic lath size was estimated by selecting blocks with lath boundaries that were normal to the foil surface and measuring the average distance between the lath boundaries (Figure 6a). The results can be seen in Table 4 and in Figure 7. The heterogeneous nature of these microstructures makes it very difficult to determine a reliable value for the lath width and it is impossible to say from the data that there is any significant difference between the Ref and Ref+V specimens. The bainite lath width actually appears to increase slightly with decreasing transformation temperature for the Ref alloy, which is contrary to current theory [25]. However, the scatter in the measurement is large and the trend could be misleading.

**Figure 7.** STEM measurement of bainite lath widths.

Figure 8a shows a bright field image of plate bainite containing RA films in Ref+V transformed at 400 ◦C. Figure 8b is the same area superimposed with a vanadium EDX chemical map and Figure 8c shows a carbon EDX map of the same region. Even although EDX carbon maps are not quantitative, carbon enrichment in retained austenite films is qualitatively obvious. Some coarse and fine precipitates are visible in Figure 8a. The chemical map in Figure 8b confirms that the majority of these particles contain vanadium. Interestingly, there is little evidence from carbon or nitrogen (not shown) maps that these are classical carbide or nitride precipitates. The largest of these particles are aligned, probably along austenite grain boundaries formed during hot rolling as their position is clearly uncorrelated with either the bainitic structure or the retained austenite films. In fact, the V precipitate density in retained austenite appears to be quite low, compared to the surrounding ferrite.

**Figure 8.** (**a**) Bright field STEM image of Ref+V transformed at 400 ◦C; (**b**) same as Figure 8a with a superimposed energy dispersive X-ray analysis (EDX) vanadium map (red) and (**c**) EDX carbon map (qualitative) from the same area.

Most of the particles observed in Ref+V specimens transformed at 400 ◦C and 430 ◦C contained vanadium with small amounts of carbon and nitrogen being occasionally detected, especially for the larger precipitates (Figure 9e). However, at 375 ◦C some of the particles contained V+Mo (Figure 9a–d). An equilibrium calculation (not shown) predicted that Mo enrichment of V(C,N) is only significant at temperatures below 800 ◦C, suggesting that these particles formed in ferrite.

**Figure 9.** Ref+V alloy transformed at 375 ◦C, (**a**) STEM bright field image, (**b**) EDX carbon map, (**c**) EDX vanadium map, (**d**) EDX molybdenum map from the same region, (**e**) STEM bright field image and EDX maps of a V(C,N) particle and (**f**) STEM dark field imaging and convergent beam electron diffraction (CBED) pattern confirming the crystal structure of the particle in e).

#### Measuring the Local Carbon Contents

The local carbon content, *C*γ, in a number of retained austenite and martensite islands was measured using a standards-based EELS technique [18]. At a carbon concentration level of 0.8 wt. % the EELS measurement precision was ±4% and the detection limit for solute carbon was 0.03 wt. %. The ex-situ EELS results for *C*γ were averaged and compared in Figure 10a with the in-situ HEXRD synchrotron data from Figure 4b and Table 5. Note that the error in the HEXRD carbon measurement (±0.011 wt. %) is too small to show in Figure 10a. The error bars in the EELS data are not measurement

errors but rather represent the standard deviation in the carbon contents obtained from different austenite islands. Clearly, there is a wide island-to-island scatter in *C*γ, as previously noted for TRIP steels [18].

From Figure 10a, the global HEXRD measurement and the local EELS values both predict the same trend of decreasing *C*<sup>γ</sup> with increasing transformation temperature. For Ref+V (yellow and brown points), the mean EELS data show very close agreement with HEXRD, especially bearing in mind the scatter in RA island carbon contents. The agreement for Ref (light blue and dark blue points) was slightly less good, but the overall trend was the same. In Figure 10b the mean carbon content in martensite Cα is compared with C<sup>γ</sup> from the EELS data. As expected, martensite islands contained less carbon than stable retained austenite. However, it is interesting to note that on average Ref+V martensite contained less carbon than Ref, i.e., Ref+V austenite appeared to be more stable, especially at lower transformation temperatures, and both alloys show *C*α (375 ◦C) < *C*α (400 ◦C). These observations clearly show that RA stability is not solely determined by carbon content [29].

If all the available carbon diffused into retained austenite, then the lever rule would predict significantly higher concentrations, *C*γ than those found in Table 5. In the absence of any significant precipitation this implies that bainitic ferrite is supersaturated in carbon. EELS analyses were therefore carried out to try to measure *Cb*, the bainitic ferrite carbon content. These measurements were made in regions of bainitic ferrite where there were no obvious signs of carbide precipitation. It was observed that *Cb* was heterogeneous—in areas close to austenite laths usually no carbon was detected (i.e., *Cb* < 0.03 wt. %). The experimental values for *Cb* given in Table 5 and plotted in Figure 10b were determined by averaging EELS measurements from bainitic ferrite far from austenite laths. The experimental *Cb* results can be compared with values calculated from a simple carbon mass balance using the phase fraction data in Table 3 and the experimental values for *C*γ and *C*α in Table 5. The two show very good agreement (Table 5). Evidently, bainitic ferrite remained highly supersaturated in carbon for both alloys, hence the tetragonality problem encountered with the W–H analysis. *Cb* increased as the transformation temperature decreased, in accordance with theory [30]. There appeared to be a slight influence of vanadium on *Cb* at 375 ◦C and 400 ◦C but the difference was within the experimental error. The carbon mass balance calculations in Table 5 suggest that vanadium additions increased *Cb*. More experimental work is required to confirm this.

**Figure 10.** (**a**) Comparison of HEXRD and electron energy loss spectroscopy (EELS) carbon contents in austenite and (**b**) EELS carbon contents of austenite, martensite and bainitic ferrite. The error in the HEXRD carbon measurements is ±0.011 wt. % The error bars in the EELS data show the standard deviation in the range of the measured values.


**Table 5.** Final carbon contents in austenite, martensite and bainitic ferrite.

#### **4. Discussion**

The yield stress σ*<sup>e</sup>* of a single phase fully bainitic structure can be expressed as the sum of contributions:

$$
\sigma\_{\ell} = \sigma\_0 + \sigma\_{\aleph} + \sigma\_d + \sigma\_p + \sigma\_L \tag{2}
$$

where σ<sup>0</sup> is the Peierls stress for pure iron, σ*ss* is the contribution from solid solution hardening, σ*<sup>d</sup>* is the strengthening due to forest dislocations, σ*<sup>p</sup>* is due to precipitation strengthening and σ*<sup>L</sup>* is the hardening due to grain or lath size dependence. As the interest here lies in explaining the difference in strength between Ref and Ref+V shown in Figure 3, the first two terms on the right side of Equation (1) can be neglected by assuming there is no difference in *Cb* and that small amounts of vanadium in solid solution have a negligible influence so that:

$$
\Delta \sigma\_{\mathfrak{e}} = \Delta \sigma\_{\mathfrak{d}} + \Delta \sigma\_{\mathfrak{p}} + \Delta \sigma\_{\mathfrak{L}} \tag{3}
$$

Figure 5a,b show that the dislocation term Δσ*<sup>d</sup>* is very small. This leaves two terms to evaluate; the contribution from precipitation strengthening, Δσ*<sup>p</sup>* and a size dependent term, Δσ*L*. The precipitation term consists of two distinct components; strengthening from (a) precipitates inherited from prior processing steps and (b) precipitates formed at the isothermal transformation temperature. The two cases are considered separately below.

#### *4.1. Strengthening from V Precipitation in Austenite and Ferrite*

The presence of 1.8 wt. % silicon strongly inhibits cementite formation in these steels [12]. This was confirmed by the lack of evidence for Fe3C in the TEM images and the absence of any cementite reflections in the HEXRD data. In view of the scarcity of these precipitates, there is no evidence to suggest that vanadium additions have a significant influence on their size, volume fraction or distribution or that they make a measurable contribution to the mechanical properties of these alloys. Therefore, cementite precipitation can be omitted from further discussion.

In order to restrict austenite grain growth the choice was made to austenitize at 900 ◦C, which is below the calculated solution temperature of 1067 ◦C for V(C,N) in Ref+V. Consequently, some undissolved vanadium-containing particles are present before the bainitic transformation occurs (Figures 8 and 9). In those images a majority of rather coarse (20–40 nm) vanadium particles can be seen, along with some finer (5–10 nm) precipitates. It must be emphasized that many areas showed little or no evidence of V precipitation, so that the average particle density observed in the TEM was actually quite low. Further, a close examination of the carbon and nitrogen maps suggests that many of these vanadium-rich regions are not carbides or nitrides and in fact could be due to local segregation.

An anisothermal kinetic precipitation calculation was carried out with the MatCalc 6 code in order to determine the amount of V(C,N) remaining undissolved before the bainitic transformation. Figure 11a shows the thermal path of the simulation. Starting at 1220 ◦C, hot rolling is represented by a uniform cooling rate of 10 ◦C/s in the austenite domain. This is decreased to 2 ◦C/s to represent air cooling in ferrite from 700 ◦C to room temperature. The material is then reheated to 900 ◦C at a

constant rate of 3 ◦C/s and held for 120 s before quenching. Two types of vanadium precipitates were simulated; the first class nucleates in austenite on grain boundaries and triple points (assuming a 100 μm austenite grain size), the second nucleates in ferrite on dislocations and defects. The dislocation density in ferrite is assumed to be low, 1 <sup>×</sup> 1012 m−<sup>2</sup> due to the slow air cooling after finish rolling. No attempt was made to include strain induced precipitation in austenite as previous experimental studies indicate that very little vanadium is actually precipitated before this point [3,4]. Figure 11b shows the evolution of the phase fractions, Figure 11c the number density and Figure 11d the mean radii of the two types of precipitate.

At the end of re-austenitization, the mean diameter of the GB V(C,N) precipitates was predicted to be 120 nm with a very low number density of 1 <sup>×</sup> 1015/m3. The particles nucleated in ferrite were much smaller with a mean diameter of 9 nm and higher number density of 1.6 <sup>×</sup> <sup>10</sup>19/m3. The total amount of V and C precipitated was actually very small: 0.0006 wt. % V and 0.0002 wt. % C. Compared to the TEM images, the calculated GB particle size appeared to be overestimated. However, the size of the dislocation nucleated precipitates correlated well with the thin foil images. Given their large size and scarcity, the strengthening effect of the GB particles can be neglected.

The contribution to strengthening made by the smaller precipitates can be estimated using the modified Ashby–Orowan equation proposed by Gladman [31]:

$$
\sigma\_p = 0.538Gb \left(\frac{f^{0.5}}{X}\right) \ln\left(\frac{X}{2b}\right) \tag{4}
$$

where *G* is the shear modulus (80,000 MPa), *b* is the Burgers vector (2.5 <sup>×</sup> 10−<sup>4</sup> <sup>μ</sup>m), *f* is the volume fraction of precipitates and *X* is the mean particle diameter in μm. Substituting the appropriate values from Figure 11 gives σ*<sup>p</sup>* = 8.5 MPa. Thus, the strength contribution from undissolved vanadium precipitates formed before the end of the re-austenitization step is expected to be very small.

**Figure 11.** (**a**) Simulation of the thermal path for V(C,N) precipitation in austenite and ferrite. (**b**) Calculated phase fractions of V(C,N) nucleated on GB and dislocation sites. (**c**) Precipitate number densities. (**d**) Mean radii.

#### *4.2. Strengthening from V Precipitation in Bainite*

Here the question of whether or not new vanadium precipitates can nucleate and grow during the fast cooling and isothermal bainitic transformation stage is addressed. It can be assumed that any further growth of the pre-existing vanadium particles is limited, as this would require low temperature long-range diffusion of vanadium across the (presumably) depleted zones surrounding those precipitates. Figure 12a is a calculated precipitation-time-temperature (PTT) diagram for V(C,N) formation in Ref+V calculated using the MatCalc 6 code [32] assuming a constant dislocation density of 1 <sup>×</sup> 1014/m<sup>2</sup> in ferrite. The blue dots show the cooling path of the Ref+V alloy. The high defect density decreased the nucleation time for V(C,N) by several orders of magnitude, but it is still not possible to nucleate vanadium particles in any ferrite that might form during cooling.

A PTT section from a lower temperature region corresponding to the isothermal bainitic transformations is shown in Figure 12b. Here the dislocation density was increased to 5 <sup>×</sup> 1014/m2, which is a reasonable value for lower bainite [1]. The simulation predicts that no V(C,N) could nucleate during the transformation times tested here. In fact, the temperature needs to be raised to almost 500 ◦C before any detectable precipitation can occur. The evidence points to the conclusion that no new vanadium precipitation occurred during the bainitic transformations studied here and that the precipitates observed in Figures 8 and 9 were all formed prior to bainite. The Δσ<sup>p</sup> strengthening term is thus a constant for all three transformation temperatures.

**Figure 12.** (**a**) V(C,N) PTT in Ref+V ferrite containing 1 <sup>×</sup> <sup>10</sup>14/m<sup>2</sup> dislocations and **(b**) low temperature PTT for V(C,N) in Ref+V bainitic ferrite containing 5 <sup>×</sup> 1014/m2 dislocations.

Even if classical precipitates do not form, the possibility that vanadium could influence the mechanical properties by segregating or clustering at a finer scale during the bainitic transformation is not precluded. Sourmail et al. [3] could not find any vanadium segregation at all in a medium carbon forging steel transformed at 450 ◦C using 3D tomographic atom probe reconstructions. However, many vanadium-rich areas that did not appear to contain either C or N were observed and these could well be an indication of clustering or segregation to boundaries and defects. Further investigations, beyond the scope of this work, are required.

#### *4.3. Size-Dependence E*ff*ects*

There is no doubt that the bainite matrix of both alloys becomes significantly stronger as the transformation temperature decreases from 430 to 375 ◦C (Figure 3). Following the current literature, a combination of packet and/or block size and/or lath refinement should be observed as the transformation temperature decreases. There are many arguments as to what is the best effective length L to use and whether yielding has a Hall–Petch or simple reciprocal type dependence on L [23]. The EBSD analysis of low angle (θ < 5◦) misorientations representing lath boundaries and the TEM analysis of lath sizes (Table 4) did not show the expected monotonic temperature-dependence (although it should be stated that the size variations were subtle and the trends might be being distorted by experimental error). It can be concluded that, based on the experimental data, the behavior shown in Figure 3 could not be explained by variations in the Δσ<sup>L</sup> term.

#### *4.4. Strengthening from Minority Phases*

The hardness measurements in Figure 3 involved plastic deformation and hence depended to some extent on the work hardening behavior of the alloys. The latter is greatly influenced by the minority retained austenite and martensite phases. In particular, the stability of retained austenite islands will play an important role. Tables 3–5 show that vanadium additions introduced four microstructural changes in Ref+V that may have a bearing on this:


Points 1–3 indicate that retained austenite in Ref+V transformed at low temperatures must be being stabilized by a strong size effect, i.e., the RA stability is not governed uniquely by the local chemistry [29]. This is very striking at 375 ◦C where the mean Ref+V thermal martensite carbon content was much lower than the Ref alloy (Figure 10b). From this together with point 4 two different transformation temperature dependent strengthening mechanisms to qualitatively explain the hardness response in Figure 3 can be proposed:


#### **5. Conclusions**

The influence of vanadium additions on isothermally formed bainite in medium carbon steels containing retained austenite was investigated using in-situ HEXRD and ex-situ electron microscopy techniques. In the range of transformation temperatures investigated (375–430◦C) the main conclusions are:


It is the author's opinion that the main benefit of adding vanadium to these steels is to be found in the refinement of martensite and austenite islands at low transformation temperatures. This provides a strong size effect stabilization of RA that could result in a significant increase in the ductility. Further, refining of blocky martensite should improve the fracture toughness behavior. If more strengthening is required then a combination of vanadium with higher nitrogen levels should be studied.

**Author Contributions:** Conceptualization, C.S.; methodology, G.G., S.Y.P.A., C.S., F.F; formal analysis, I.P.; investigation, I.P., B.S.-A., G.G., S.Y.P.A., M.S.; writing—original draft preparation, I.P.; writing—review and editing, C.S., S.Y.P.A.; project administration, F.F. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Canadian Office of Energy Research and Development (OERD) and the LABEX DAMAS (ANR-11-LABX-008-01) from Lorraine, France.

**Acknowledgments:** The authors would like to express their gratitude to the Canadian Office of Energy Research and Development (OERD) and the LABEX DAMAS (ANR-11-LABX-008-01) from Lorraine, France for their support. Finally, we wish to acknowledge the help and professionalism of all the staff at CanmetMATERIALS who participated in the study and in particular Jian Li for aid with the EBSD analysis.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

MDPI St. Alban-Anlage 66 4052 Basel Switzerland Tel. +41 61 683 77 34 Fax +41 61 302 89 18 www.mdpi.com

*Metals* Editorial Office E-mail: metals@mdpi.com www.mdpi.com/journal/metals

MDPI St. Alban-Anlage 66 4052 Basel Switzerland

Tel: +41 61 683 77 34 Fax: +41 61 302 89 18