Acessibilidade / Reportar erro

Numerical and Experimental Study of an Industrial Case for Grain Size Evolution in Bainitic Steel in Controlled Hot Forging and its Influence on Mechanical Performance

Abstract

Controlling the recrystallization is an important way to reach grain size refinement and outstanding strength and toughness on alloy metals. This study sets out the application and investigation of mathematical microstructure modeling of a newly designed bainitic steel for hot forging industrial applications. The macro-scale model was used to observe and predict the austenitic grain size behavior during the controlled forging of a gear. Arrhenius grain growth kinetic and recrystallization model for a new class of bainitic steel was established for the given strain rate ranges and temperatures. This model was calibrated through microscopic analysis and used to simulate the unpublished constants of low alloyed bainitic forging steel DIN 18MnCrSiMo6-4 microstructure module using DEFORM® commercial finite element code. The increased temperature due to the adiabatic effect was investigated by numerical analysis, demonstrating its influence on grain coarsening. Local tensile test and Charpy-V notch were compared at different industrial hot forging temperatures and local plastic strain. Changes in yield strength and ductility have demonstrated the grain size influence on the processing parameters. The employed numerical model was an efficient tool to predict and present an alternative path to develop robust industrial forging using semi-empirical models.

Keywords:
Hot Forging; Grain Size; JMAK; Mechanical Properties; Bainitic Steel

1. Introduction

Engineers and scientists have transformed how low alloyed carbon steels are designed and processed in the last decades. The relationship between tailored metals processing, microstructure, and mechanical properties remains challenging, especially when more energy-efficient and environmentally friendly products are on the global manufacturing agenda. Advanced bainitic steels have stood out in the last decades11 Shimanov M, Korpala G, Terzic A, Kawalla R. Bainitic steels, their characteristics and applications. KEM. 2016;684:104-10. https://doi.org/10.4028/www.scientific.net/KEM.684.104.
https://doi.org/10.4028/www.scientific.n...

2 Caballero FG, Capdevila C, Chao J. The microstructure of continuously cooled tough bainitic steels. In: 2nd International Conference of Super-High Strength Steels; 2010 Oct 27-30; Verona, Italy. Proceedings. Verona, Italy: TIB; 2010. p. 10-8.

3 Roelofs H, Hasler S, Urlau U, Lembke MI, Olschewski G. Continuously cooled bainitic steel HSX ® Z12: one decade of experience. In: 4th International Conference on Steels in Cars and Trucks; 2014, Jun 15-19; Braunschweig, Germany. Proceedings. Braunschweig, Germany.
-44 Caballero FG, Garcia-Mateo C, Cornide J, Allain S, Puerta J, Crouvizer M, et al. New advanced ultra high strength bainitic steels: ductility and formability (DUCTAFORM) [Internet]. Brussels: European Commission; 2013 [cited 2021 Nov 24]. Available from: https://op.europa.eu/en/publication-detail/-/publication/9f8dc6ba-d882-48da-9b7b-aa8d7f00dd24
https://op.europa.eu/en/publication-deta...
due to the excellent combinations of mechanical properties such as yield strength and toughness caused by their microconstituents. One way to obtain such mechanical properties is prior austenitic grain size modification. It is known that coalescence of the bainite laths is avoided by grain refinement, increasing the volume fraction and stability of retained austenite in isothermal heat treatment55 Lan HF, Du LX, Li Q, Qiu CL, Li JP, Misra RDK. Improvement of strength-toughness combination in austempered low carbon bainitic steel: the key role of refining prior austenite grain size. J Alloys Compd. 2017;710:702-10.. Although it is an excellent method for this diffusional transformation, its mechanical performance is achieved at the cost of high energy consumption, making it less attractive in hot forged parts manufacturing chains.

As an alternative, the thermomechanical process in continuous cooling bainitic steels can lead to energy savings and tailored mechanical performance by controlling forging and cooling conditions, replacing long isothermal treatments such as quenching and tempering (Q&T)22 Caballero FG, Capdevila C, Chao J. The microstructure of continuously cooled tough bainitic steels. In: 2nd International Conference of Super-High Strength Steels; 2010 Oct 27-30; Verona, Italy. Proceedings. Verona, Italy: TIB; 2010. p. 10-8.,44 Caballero FG, Garcia-Mateo C, Cornide J, Allain S, Puerta J, Crouvizer M, et al. New advanced ultra high strength bainitic steels: ductility and formability (DUCTAFORM) [Internet]. Brussels: European Commission; 2013 [cited 2021 Nov 24]. Available from: https://op.europa.eu/en/publication-detail/-/publication/9f8dc6ba-d882-48da-9b7b-aa8d7f00dd24
https://op.europa.eu/en/publication-deta...
. Therefore, understanding how these steels perform in hot forging applications is undoubtedly necessary to achieve excellent mechanical properties. Efforts in understanding the bainitic transformation after the lower plastic strain was reported by66 Lemos Bevilaqua W, Epp J, Meyer H, Da Silva Rocha A, Roelofs H. In situ investigation of the bainitic transformation from deformed austenite during continuous cooling in a low carbon Mn-Si-Cr-Mo steel. Metall Mater Trans, A Phys Metall Mater Sci. 2020;51(7):3627-37. using in situ XRD experiments coupled with a dilatometer. Microstructure transformations in continuous cooling were accelerated after compression, and simultaneously, the granular bainite and retained austenite were refined. Another work77 Hatwig RA, Dong J, Epp J, Rocha A S. Effect of compressive deformations on the final microstructure of a low carbon high silicon bainitic steel thermomechanically processed. Mater Res. 2021;24(1):1-7. showed the effect of the plastic strain of metastable austenite on resulting bainite morphology by physical simulation. It has presented the doorway for the forging application using the DIN 18MnCrSiMo6-4 steel in different thermomechanical routes.

Silveira and co-workers88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
found that, besides the forging having a significant impact on the bainite morphology, at 1200 °C, a predominantly bainitic microstructure was formed from the coarse austenitic grains. They also summarized that plastic deformation at high temperatures inhibits the formation of pro-eutectoid ferrite, thus favoring the austenite transformation into bainite77 Hatwig RA, Dong J, Epp J, Rocha A S. Effect of compressive deformations on the final microstructure of a low carbon high silicon bainitic steel thermomechanically processed. Mater Res. 2021;24(1):1-7.,88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
. When the temperature was reduced to 1100 °C or 1000 °C, in addition to displaying more refined grains, there was an increase in the amount of polygonal ferrite and a reduction in retained austenite fraction. In summary, the austenitic grain significantly impacts the final microstructure of DIN 18MnCrSiMo6-4 steel.

In this way, and as acquired in an early work99 Ivaniski TM, de Souza TF, Épp J, Rocha AS. Constitutive modelling of high temperature flow behaviour for a low carbon high silicon bainitic steel. Mater Res. 2020;23(5):e20200264., the modeling of the flow curves of this material, as well as the (JMAK) "Johnson-Mehl-Avrami-Kolmogorov" parameters1010 Jonas JJ, Quelennec X, Jiang L, Martin É. The Avrami kinetics of dynamic recrystallization. Acta Mater. 2009;57(9):2748-56., are necessary for assessing its microstructural evolution such as in recrystallization kinetics. Thus, enabling the application of "semi-empirical" numerical models, based on methods previously proposed by1111 Irani M, Lim S, Joun M. Experimental and numerical study on the temperature sensitivity of the dynamic recrystallization activation energy and strain rate exponent in the JMAK model. J Mater Res Technol. 2019;8(2):1616-27. https://doi.org/10.1016/j.jmrt.2018.11.007.
https://doi.org/10.1016/j.jmrt.2018.11.0...

12 Razali MK, Joun M. A new approach of predicting dynamic recrystallization using directly a flow stress model and its application to medium Mn steel. J Mater Res Technol. 2021;11(18):81-94. http://dx.doi.org/10.1016/j.jmrt.2021.02.026.
http://dx.doi.org/10.1016/j.jmrt.2021.02...

13 Chamanfar A, Chentouf SM, Jahazi M, Lapierre-Boire L-P. Austenite grain growth and hot deformation behavior in a medium carbon low alloy steel. J Mater Res Technol. 2020;9(6):12102-12114. https://doi.org/10.1016/j.jmrt.2020.08.114.
https://doi.org/10.1016/j.jmrt.2020.08.1...

14 Moraes ALI, Balancin O. Numerical simulation of hot closed die forging of a low carbon steel coupled with microstructure evolution. Mater Res. 2015;18(1):92-7.

15 Ivaniski TM, Epp J, Zoch H-W, Da Silva Rocha A. Austenitic grain size prediction in hot forging of a 20mncr5 steel by numerical simulation using the JMAK model for industrial applications. Mater Res. 2019;22(5)

16 Sellars CM, Whiteman JA. Recrystallization and grain growth in hot rolling. Met Sci. 1978;13(3-4):187-94. https://doi.org/10.1179/msc.1979.13.3-4.187.
https://doi.org/10.1179/msc.1979.13.3-4....
-1717 Henke T, Hirt G, Bambach M. Optimization of a closed die forging process to manufacture a gear wheel by the use of a response surface model. Adv Mat Res. 2014;922:254-9., which can simulate the complexity of the hot forging process and the grain size evolution.

Although the new methodology1111 Irani M, Lim S, Joun M. Experimental and numerical study on the temperature sensitivity of the dynamic recrystallization activation energy and strain rate exponent in the JMAK model. J Mater Res Technol. 2019;8(2):1616-27. https://doi.org/10.1016/j.jmrt.2018.11.007.
https://doi.org/10.1016/j.jmrt.2018.11.0...
,1313 Chamanfar A, Chentouf SM, Jahazi M, Lapierre-Boire L-P. Austenite grain growth and hot deformation behavior in a medium carbon low alloy steel. J Mater Res Technol. 2020;9(6):12102-12114. https://doi.org/10.1016/j.jmrt.2020.08.114.
https://doi.org/10.1016/j.jmrt.2020.08.1...
is excellent for predicting grain size evolution and flow curves through the JMAK equation, it has not yet been applied in industrial hot forging in high strain rate presses (above 50 s-1) in this bainitic steel. Also, the kinetics of austenitic grain growth industrial closed die hot forging after load release has not been discussed using a finite element FEM method approach. A higher strain rates increase dislocation density and stored energy during the plastic deformation, affecting recrystallization kinetics and possibly the grain growth1010 Jonas JJ, Quelennec X, Jiang L, Martin É. The Avrami kinetics of dynamic recrystallization. Acta Mater. 2009;57(9):2748-56.. In closed die forging, such behavior is expected. Therefore, in the present research, the suitability of these proposed methods was addressed in industrial hot forging conditions in a high strain rate press.

This work aimed at validating a macroscale model of the grain size evolution using numerical simulation for a new continuously cooled bainitic steel DIN 18MnCrSiMo6-4 hot forged. Some aspects regarding recrystallization during forging were addressed and what effects the process caused on grain growth. The synergetic effect between closed die forging parameters in response to austenitic grain size at the end of the process has been verified. Finally, it has compared two temperature levels on the mechanical properties.

2. Material and Methods

2.1. Grain growth kinetic of annealing test

All experiments were performed with the steel DIN 18MnCrSiMo6-4, classified as a low-carbon carbide-free bainitic steel grade. The chemical composition is given in Table 1.

Table 1
Chemical composition for steel 18MnCrSiMo6-4, mass percentage.

The evolution of the austenitic grain size of this steel was studied by88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
, considering the same sample mass of 515.47 g in a cylindrical geometry under isothermal conditions. The authors heated the samples at temperatures of 1000 °C, 1100 °C, and 1200 °C. It was found that the abnormal grain growth occurred at 1200 °C, which was attributed to the dissolution of precipitates formed by microalloying elements such as Ti, N, and Al in this steel. The reported results significantly explain the grain growth equation values, where the methodology is shown below to be discussed in the results section. Figure 1 shows the as received prior austenitic grain size PAGS and the microstructure of the steel revealed as the same methodology written in section 2.3.

Figure 1
(a) 1 PAGS as received 23,5 μm, (b) 2 Initial Microstructure BF – Bainitic Ferrite, GB – Granular Bainite, PF – Polygonal Ferrite, A – Retained Austenite, M – Martensite.

The constitutive Arrhenius equation described grain growth kinetics1313 Chamanfar A, Chentouf SM, Jahazi M, Lapierre-Boire L-P. Austenite grain growth and hot deformation behavior in a medium carbon low alloy steel. J Mater Res Technol. 2020;9(6):12102-12114. https://doi.org/10.1016/j.jmrt.2020.08.114.
https://doi.org/10.1016/j.jmrt.2020.08.1...
. In addition, the activation energy for this bainitic steel was determined according to Equation 1 implemented in DEFORM® V. 11.0 subroutine.

d = d 0 m + a 9 t e x p Q g g R T 1 / m (1)

The initial and final grain sizes are given by d0 and d, while Qgg describes the activation energy necessary for grain growth. The inverse of the time t exponent for grain growth is given by m, and a9 is a constant that depends on the chemical composition of the steel1818 Janssens KGF, Raabe D, Kozeschnik E, Miodownik MA, Nestler B. Computacional materials engineering: an introduction to microstructure evolution. Cambridge: Academic Press; 2007. 359 p.. The R-value is the universal gas constant (8.314 J/mol.K), and T the absolute temperature in Kelvin (K). By applying natural logarithm to both sides of Equation 1, we can determine m (Equation 2) and Q9 (Equation 3) values, and it can be expressed as1313 Chamanfar A, Chentouf SM, Jahazi M, Lapierre-Boire L-P. Austenite grain growth and hot deformation behavior in a medium carbon low alloy steel. J Mater Res Technol. 2020;9(6):12102-12114. https://doi.org/10.1016/j.jmrt.2020.08.114.
https://doi.org/10.1016/j.jmrt.2020.08.1...
:

1 m = l n d 1 t T (2)
Q 9 = m R l n d l n T t (3)

In this case, it is essential to emphasize that the equation describing growth does not consider prior plastic deformation. Moreover, the grain size distribution was assumed to be homogeneous since the grain size measurements were the same in distinct parts of the cylindrical rod obtained.

2.2. Flow curves and microstructure analysis

Hot isothermal compression tests for flow curves determination were conducted in a DIL 805A/T dilatometer. Cylindrical samples with 10 mm x 5 mm were heated at a rate of 10 °C/s and held for five (5) minutes for homogenization before compression. Strain rates of 0.1 s-1, 0.5 s-1, 1 s-1, and 5 s-1 and temperatures of 850 °C, 950 °C, 1050 °C, and 1150 °C with a fixed true plastic strain of 0.69 were set up for testing.

2.3. Metallography analysis

After isothermal compression and quenching, the austenitic grain size samples were prepared and analyzed by optical microscopy according to DIN EN ISO 643 and ISO 14250 using the line-intercept method. A saturated picric acid solution revealed the grain boundary (AGB) with 42 mL of wetting agent, 58 mL of distilled water, and 2.3 g of picric acid. Wetted cotton was swabbed into the steel surface for 5 min to reveal the (AGB). The overall microstructure was etched by immersion in a Nital 2% solution for 10 s.

Electron backscatter diffraction (EBSD) was performed using a Philips XL 30 to determine the distribution of the crystallographic units. The instrument operated at 25 keV, with a working distance of 8.4 mm, a sample tilt angle of 75 degrees, and a step size of 0.1 μm. The data were processed in Orientation Imaging Microscopy for two quenched samples in temperatures of 950 °C and 1150 °C at ε˙= 0.1 s-1 of plastic strain rate.

2.4. Recrystallization model

The modeled flow curves obtained by compression tests and the JMAK recrystallization coefficients were implemented in DEFORM® finite element code V11.0. Figure 2 shows the flowchart used for the acquisition of the coefficients.

Figure 2
Flowchart of modeling interaction to calibrate the experiment.

The Zener-Hollomon parameter Z was calculated using Equation 4, and the overall equation that depends on it is shown in Equation 5.

Where Z is:

Z = ε ˙ . e x p Q R T (4)

In Equation 5, n, α, and A are the material-dependent constants, determined by experimental flow curve results, and Q is the activation energy in (J/mol.K) for the onset of dynamic recrystallization.

The parameters determined from flow stress curves for the DIN 18MnCrSiMo6-4 steel were α = 0.009833 MPa-1, n = 5.3, A(T. ε¯) = 7.6 x1011and QTotal = 276684 J/mol.K99 Ivaniski TM, de Souza TF, Épp J, Rocha AS. Constitutive modelling of high temperature flow behaviour for a low carbon high silicon bainitic steel. Mater Res. 2020;23(5):e20200264.. The peak stress σp was defined by deriving the stress and strain from the flow curves (θ=σε).

ε ˙ = A sinh α σ p n e x p Q R T (5)

The modeled flow stress in the subroutine of the finite element code is shown in Equation 6, where the work hardening σWH represents the stress caused by the increase in dislocation motion of the material up to σp.

σ W H = σ 2 s a t σ 2 s a t σ 0 2 exp r ε (6)

Where σWH is derived from the saturation stress σsat, which is the first stage in which flow stress increases. Where r is the dynamic recovery coefficient that decreases as the strain rate increases in DIN 18MnCrSiMo6-4 bainitic steel99 Ivaniski TM, de Souza TF, Épp J, Rocha AS. Constitutive modelling of high temperature flow behaviour for a low carbon high silicon bainitic steel. Mater Res. 2020;23(5):e20200264., this coefficient has an essential function in this work, as it improves the model accuracy at high strain rates in the hot forging process.

Since these coefficients are obtained by linear fitting from the experimental results already determined by99 Ivaniski TM, de Souza TF, Épp J, Rocha AS. Constitutive modelling of high temperature flow behaviour for a low carbon high silicon bainitic steel. Mater Res. 2020;23(5):e20200264., the stress model can be calculated via Equation 7. The σss represents the stress at a steady-state when complete dynamic recrystallization XDRX occurs.

σ m o d = σ W H σ s a t σ s s X D R X (7)

A more detailed description of these hyperbolic constitutive equations and how these results were obtained for this steel can be found in99 Ivaniski TM, de Souza TF, Épp J, Rocha AS. Constitutive modelling of high temperature flow behaviour for a low carbon high silicon bainitic steel. Mater Res. 2020;23(5):e20200264..

The mathematical model of the critical strain before reaching a steady-state is presented in the following equation (Equation 8) and has been proposed in several previous works1616 Sellars CM, Whiteman JA. Recrystallization and grain growth in hot rolling. Met Sci. 1978;13(3-4):187-94. https://doi.org/10.1179/msc.1979.13.3-4.187.
https://doi.org/10.1179/msc.1979.13.3-4....
,1818 Janssens KGF, Raabe D, Kozeschnik E, Miodownik MA, Nestler B. Computacional materials engineering: an introduction to microstructure evolution. Cambridge: Academic Press; 2007. 359 p.,1919 Lin J, Balint D, Pietrzyk M. Microstructure evolution in metal forming processes: modelling and applications. Sawston, Reino Unido: Woodhead Publishing; 2012. p. 67-74. for C-Mn steels. When the critical strain (εc) is exceeded, as described in Equation 7, new grains nucleate and grow during the strain step only for ε¯>εc.

ε c = a 2 ε p (8)

In Equation 8 and Equation 9, a1,a2, n1, m1, and εp are the dependent material parameters and the peak strain, respectively, which is equal to the peak stress σp derived from the flow curves.

The activation energy Q1 (J.mol-1), gas constant R (J.mol-1. K-1), and the initial grain size d0 depend on the effective strain rate ε¯, and the temperature T. Thereby, Equation 8, Equation 9, and Equation 10 describe the overall kinetics of dynamic recrystallization XDRX.

ε p = a 1 d 0 n 1 ε ¯ m 1 e x p Q 1 R T (9)
X D R X = 1 e x p β d ε ε p a 10 ε ¯ 0.5 k d (10)

In Equation 10, the exponent (βd) represents the Avrami coefficient which is sensitive to variations in temperature. The Avrami exponent (kd) is determined by polymorphic changes, discontinuous precipitation, and interfacial growth control, among other factors1818 Janssens KGF, Raabe D, Kozeschnik E, Miodownik MA, Nestler B. Computacional materials engineering: an introduction to microstructure evolution. Cambridge: Academic Press; 2007. 359 p.,1919 Lin J, Balint D, Pietrzyk M. Microstructure evolution in metal forming processes: modelling and applications. Sawston, Reino Unido: Woodhead Publishing; 2012. p. 67-74.. The ε¯0.5 described in Equation 11 is the strain when 50% of the recrystallized fraction is reached, where a5, h5, m5, and n5 are material parameters.

ε ¯ 0.5 = a 5 d 0 h 5 ε ¯ m 5 ε n 5 e x p Q 5 R T + C 5 (11)

Finally, Equation 12 and Equation 13 provide the recrystallized grain size DDRX after full XDRX and the average grain size DAVG, respectively.

D D R X = a 8 d 0 h 8 ε ¯ m 8 ε n 8 exp Q 8 R T (12)
D A V G = X D R X D D R X + ( 1 X D R X ) d 0 (13)

If, after the interruption of deformation, the cumulative strain during the industrial forging process is higher than the critical strain εc then metadynamic recrystallization occurs1818 Janssens KGF, Raabe D, Kozeschnik E, Miodownik MA, Nestler B. Computacional materials engineering: an introduction to microstructure evolution. Cambridge: Academic Press; 2007. 359 p.. In MDRX, new nuclei appear during forming, but grain growth occurs after the onset of plastic strain. Equation 14 and Equation 15 describe the kinetics:

X M D R X = 1 exp ( β m t sec t 0.5 k m ) (14)
t 0.5 = a 4 ε ¯ n 4 d o h 4 ε ¯ m 4 e x p Q 4 R T (15)

Where t0.5 is the time required for 50% metadynamic recrystallization, which is dependent on the effective strain ε¯n4, initial grain size doh4 and strain rate ε¯˙m4 exponents. Equation 16 describes the grain size diameter for a complete metadynamic recrystallization after plastic strain:

D M D R X = a 7 ε ¯ h 7 d 0 h 7 ε ¯ m 7 e x p Q 7 R T (16)

Static recrystallization, described in Equation 17, occurs when the cumulative strain is lower than the critical strain εc. During this phenomenon, nuclei of recrystallized grains appear and grow after the forging step. The equation is given by:

X S R X = 1 e x p β s t t 0.5 k s (17)

Where t0.5, Equation 18 is:

t 0.5 = a 3 ε ¯ n 3 d o h 3 ε ¯ m 3 e x p Q 3 R T (18)

Grain growth for a full SRX is given in Equation 19:

D s = a 6 ε ¯ n 6 d 0 h 6 ε ¯ m 6 e x p Q 6 R T + C (19)

After full dynamic, static and metadynamic recrystallization, equiaxed grains may coarsen due to grain growth, as described in Equation 20, which is considered temperature-dependent.

d f m T = d 0 d r x m T + a 9 e x p Q 9 R T . t (20)

Here, d0 is the full-recrystallized grain size, t is the time after complete recrystallization, Q9 is the apparent activation energy, and m and a9 are the material and process-dependent constants.

2.5. Forging trial and numerical modeling

The industrial closed die forging tests were performed using a friction press with an energy of 240 kJ. Table 2 shows the forging conditions and targets used in this work. A cylindrical billet was forged at two different temperature conditions. The thermomechanical route and a forged specimen are shown schematically in Figure 3.

Table 2
Processing conditions in the industrial forging experiments.
Figure 3
(a) Thermomechanical routes on industrial environment steps, (b) Experiment forged samples indicating the regions of metallography analysis.

After forging, two dwell times, which represent the interval before quenching, were estimated with the aid of recorded videos. The first dwell time represents the interval where the lower and upper dies touch, 2 seconds. The second dwell time represents the interval in which the upper die was lifted, allowing the workpiece to be manipulated/moved. The workpiece was then quenched with moderate agitation at a water temperature of 23 °C.

After the first methodology was used for forging and quenching, other samples were forged and cooled to room temperature (25 °C) in calm air. Three representative sample regions were analyzed to determine the average grain size after hot forging. Prior austenitic grain size (PAGS) after quenching with water was revealed using a 3% picric acid solution and quantified by optical microscopy at 500x magnification. The circular intercept method was used for grain size determination of the hot-forged samples according to the ASTM E112 standard2020 ASTM: American Society for Testing and Materials. ASTM E112. West Conshohocken: ASTM; 2014..

The evolution of the grain size of steels by induction coil heating at 1030 °C and 1130 °C after 300 seconds was studied in annealing tests, which gave values of 26 μm for 1030 °C and 54 μm for 1130 °C. After removing the samples from the induction furnace, the surface temperature was measured using a Fluke®, Ti - 400 thermal camera. The forging temperatures in this work were chosen to be 1030 °C and 1130 °C based on the results of88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
to avoid abnormal austenitic grain growth.

The coupled simulation of heat transfer and microstructure evolution was performed using a rigid-viscoelastic code in DEFORM® V.11.0. The multi-operation code was used to represent the industrial conditions in the numerical simulation. Figure 4 shows the proposed forging model corresponding to the industrial conditions and the tracked points representing the average grain size of twelve points in the measured regions (R1, R2, and R3).

Figure 4
The axisymmetric model implemented in DEFORM for gear forging and a cross-section of the simulated forging profile of the preform showing the analyzed areas. 1 - Billet, 2 - Upper die, 3- Lower die, 4 - Knock-off.

Hexahedral elements were used to discretize the components using a coupled thermomechanical process. Elasto-plastic bodies were considered in the upper, lower, and knock-off dies, with 250 °C as the initial temperature in the forming dies. Since high plastic strains lead to severe mesh distortion, the automatic remeshing technique was used. This method solves problems related to the inability to model geometry and state variables1919 Lin J, Balint D, Pietrzyk M. Microstructure evolution in metal forming processes: modelling and applications. Sawston, Reino Unido: Woodhead Publishing; 2012. p. 67-74..

The boundary conditions of the process are listed in Table 3, according to the industrial environment and the steel properties. The thermal conductivity and heat capacity were taken from Castro et al.2121 Castro PJ, Silveira AC F, Ivaniski TM, Turra CJ, Epp J, Rocha A S. Two-step continuous cooling heat treatment applied in a low carbon bainitic steel. Mater Res. 2021;24(2) inverse heat transfer analysis for this steel. The friction factor µ was set at 0.3, corresponding to the use of graphite and water in the closed die hot forging condition1414 Moraes ALI, Balancin O. Numerical simulation of hot closed die forging of a low carbon steel coupled with microstructure evolution. Mater Res. 2015;18(1):92-7..

Table 3
Thermal properties in the numerical model.

2.6. Strength and toughness analysis

This study performed tensile testing of the miniature specimen to determine strength in specific forged areas, as shown in Figure 5. A hydraulic testing machine MTS 810 250 kN with a laser strain gage, was used to perform the tests. For each analyzed region, seven (7) samples were prepared for testing to achieve good repeatability of the experiments. The samples were produced by wire EDM, as suggested by2222 Kashaev N, Horstmann M, Ventzke V, Riekehr S, Huber N. Comparative study of mechanical properties using standard and micro-specimens of base materials Inconel 625, Inconel 718 and Ti-6Al-4 v. J Mater Res Technol. 2013;2(1):43-7., using the standards DIN EN 2002:2006, DIN EN ISO 6892-:2009-12, and DIN 50100. Four specimens were taken from two forged test pieces for each temperature proposed in this study to determine the bar impact energy by the Charpy notched bar impact test. The specimens were tested at a room temperature of 23.5 °C in a standard atmosphere. The specimens were processed following ASTM E22982323 ASTM: American Society for Testing and Materials. ASTM E2298. West Conshohocken: ASTM; 2018.

Figure 5
Representative and machined zone of collected samples for the miniature tensile test specimen and Charpy V-notch.

3. Results and Discussion

3.1. Grain growth model of annealing

Figure 6a and 6b show a parallel trend in grain growth with increasing time and temperature. Linear regression was used to estimate the mean values of the coefficients m, a9, and Q9 given in Equation 21.

Figure 6
Modeled Annealing curves for 18MnCrSiMo6-4 steel; (a) Natural logarithm of grain size and time, (b) Grain size and inverse of temperature, (c) Comparison between experiment88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
and simulation grain growth evolution.

From these data fitting, Figure 6c shows the results of the isothermal numerical simulation of austenitic grain growth compared to the experiments performed by88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
. In the current study, it was necessary to consider a homogeneous grain size distribution to find a reasonable agreement between the experiment and predicted results.

d = d 0 5,25 + 4.56 10 16 t e x p 299500 R T 0.1904 (21)

Figure 6c also shows that grain growth exhibits exponential behavior, with final size strongly affected by temperature after a short time. At a specific annealing time, say 1200 seconds, a variation in grain size deviates from the exponential trend. It is possible that at this point, a change in activation energy occurs, possibly changing the exponent (m) in Equation 1 of grain growth. In particular, this is caused by dissociating precipitates that reduce grain growth's boundary motion or coalescence due to Zenner pinning88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
.

It is well known that carbon diffuses rapidly because of its interstitial characteristics, facilitating dissolution in the austenitic matrix. Therefore, it requires less self-diffusion activation energy. The opposite is true in the presence of substitutional elements such as Mn and Ti, which require higher self-diffusion activation energy2424 Vasilyev AA, Sokolov SF, Kolbasnikov NG, Sokolov DF. Effect of alloying on the self-diffusion activation energy in γ-iron. Phys Solid State. 2011;53(11):2194-200..

The abrupt increase in austenite grain size at 1200 °C is the result of the dissolution of carbides formed by microalloying elements in this steel. These elements are responsible for forming precipitates such as manganese sulfides MnS and TiC, which causes pinning effects to be a barrier to grain growth88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
. Molybdenum carbides (MoC) may have formed since a small amount (~0.27% Mo) is present, affecting grain growth stabilization.

The theoretical implications may suggest that the total grain boundary area, in this case, is higher than the 1200 °C austenitization condition. Therefore, the system's total energy may be higher2424 Vasilyev AA, Sokolov SF, Kolbasnikov NG, Sokolov DF. Effect of alloying on the self-diffusion activation energy in γ-iron. Phys Solid State. 2011;53(11):2194-200.. In other words, it can be assumed that the activation energy increased.

The value of the activation energy of 18MnCrSiMo6-4 is (299 kJ/mol.K), which is slightly higher than the results calculated by1313 Chamanfar A, Chentouf SM, Jahazi M, Lapierre-Boire L-P. Austenite grain growth and hot deformation behavior in a medium carbon low alloy steel. J Mater Res Technol. 2020;9(6):12102-12114. https://doi.org/10.1016/j.jmrt.2020.08.114.
https://doi.org/10.1016/j.jmrt.2020.08.1...
(230.992 kJ/mol.K), which was determined for a steel grade with a higher carbon content of 0.3% C and lower amount of manganese 0.81% Mn. However, after the annealing test, the obtained grain size in this work is smaller than the steel studied by1313 Chamanfar A, Chentouf SM, Jahazi M, Lapierre-Boire L-P. Austenite grain growth and hot deformation behavior in a medium carbon low alloy steel. J Mater Res Technol. 2020;9(6):12102-12114. https://doi.org/10.1016/j.jmrt.2020.08.114.
https://doi.org/10.1016/j.jmrt.2020.08.1...
(about 7.8 times lower).

3.2. Model validation

This work's industrial application model has been based on validating the numerical model of grain growth after the recrystallization kinetics from the compression test. The interactive mode presented in Figure 2 determines the coupling of the flow curve and grain size model throughout the experimental results. Table 4 displays the overall constants calibrated for the microstructure modeling.

Table 4
Constants used in DEFORM® after calibration of the experiment for microstructural modeling of DIN 18MnCrSiMo6-4 steel.

The accuracy of the results of the FEM model can be seen in the flow stress shown in (Figure 7a), which shows what would be the flow curves at high strain and strain rate levels (100 s-1). These curves were converted into time (s) format using the ratio between the strain and strain rate (s-1)99 Ivaniski TM, de Souza TF, Épp J, Rocha AS. Constitutive modelling of high temperature flow behaviour for a low carbon high silicon bainitic steel. Mater Res. 2020;23(5):e20200264.,1010 Jonas JJ, Quelennec X, Jiang L, Martin É. The Avrami kinetics of dynamic recrystallization. Acta Mater. 2009;57(9):2748-56., as shown in Equation 22. In (Figure 7b), the predicted values of recrystallization fraction in different strain rates are shown.

t = ε ε c ε ˙ (22)

With the increase of strain rate at constant T (1150 °C), the XDRX time decreased1818 Janssens KGF, Raabe D, Kozeschnik E, Miodownik MA, Nestler B. Computacional materials engineering: an introduction to microstructure evolution. Cambridge: Academic Press; 2007. 359 p.. As reported by1010 Jonas JJ, Quelennec X, Jiang L, Martin É. The Avrami kinetics of dynamic recrystallization. Acta Mater. 2009;57(9):2748-56.,2525 Guo-Zheng Q, Yang W, Ying-Ying L, Jie Z. Effect of temperatures and strain rates on the average size of grains refined by dynamic recrystallization for as-extruded 42CrMo steel. Mater Res. 2013;16(5):1092-105., a high strain rate increases the migration rate of dislocations and grain boundaries, therefore, accelerating the XDRX process after the εc is reached2525 Guo-Zheng Q, Yang W, Ying-Ying L, Jie Z. Effect of temperatures and strain rates on the average size of grains refined by dynamic recrystallization for as-extruded 42CrMo steel. Mater Res. 2013;16(5):1092-105..

Figure 7
(a) Modeled flow curves and predicted for high strain rates on temperature dependence99 Ivaniski TM, de Souza TF, Épp J, Rocha AS. Constitutive modelling of high temperature flow behaviour for a low carbon high silicon bainitic steel. Mater Res. 2020;23(5):e20200264., (b) XDRX kinetics on strain rate dependence, modeled and predicted, (c) Comparison between experiment and prediction of a load of a compression test, (d) FEM recrystallization results.

Figure 7c and 7d show the predicted results of the experimental punch load of the numerical model during the compression test and FEM recrystallization results, respectively.

A good approximation of the computational simulation results with the force and displacement experiments has been shown, using the model proposed by99 Ivaniski TM, de Souza TF, Épp J, Rocha AS. Constitutive modelling of high temperature flow behaviour for a low carbon high silicon bainitic steel. Mater Res. 2020;23(5):e20200264. for a continuous cooling bainitic steel. A slight variation of the experimental forces compared to the simulation can be seen at the temperature and strain rate of 950 °C and 5 s-1, respectively. It is assumed that this condition is more prone to adiabatic heating variation for the calculated model when compared to the phenomena occurring during the experiment. The chemical interactions and the activation energy vary with the degree of plastic strain and strain rate. Therefore, this can sensitively vary the force and displacement results. Moreover, friction was kept constant during the computer simulation, not depicting the natural behavior during the hot forging process. A dynamic friction coefficient could be preferred to solve this issue.

Figure 8 shows the cross-section view of the experimentally measured and numerically calculated fully recrystallized austenitic microstructure after being compressed at 950 °C (Figure 8a) and 1150 °C (Figure 8b), respectively. Three points on the specimens were measured and compared to the simulation. What stands out in this figure is that small nucleate recrystallized grains can be observed at different temperatures, with grains being more refined at 950 °C.

Figure 8
Simulation of austenitic grain size in isothermal compression test (a) 950 °C, 0.1 (1/s), (b) 1150 °C , 0.1 (1/s).

Figure 9 illustrates the experimental and simulated austenitic grain size variation at different temperatures and strain rates. The most interesting aspect of this graph is that the temperature influenced the austenitic grain refinement and the strain rate. The values show that increasing the strain rate by 10 times at 1150 °C caused a higher variation in grain size when compared to 950 °C. The variation was 40.74% at a constant temperature of 1150 °C. A percentage variation of 82% in the grain size refinement was obtained from 1150 °C to 950 °C. However, at 950 °C, the strain rate did not influence the grain refinement considerably. That behavior could be related to a decrease in driving force where the temperature dominates the recrystallization mechanism activation99 Ivaniski TM, de Souza TF, Épp J, Rocha AS. Constitutive modelling of high temperature flow behaviour for a low carbon high silicon bainitic steel. Mater Res. 2020;23(5):e20200264.,1515 Ivaniski TM, Epp J, Zoch H-W, Da Silva Rocha A. Austenitic grain size prediction in hot forging of a 20mncr5 steel by numerical simulation using the JMAK model for industrial applications. Mater Res. 2019;22(5).

Figure 9
The predicted austenitic grain size of hot compressed samples ε = 0.69.

The EBSD map in Figure 10 shows these results in the crystallographic unit sizes (UC) of quenched samples, and Figure 11 shows a histogram of the UC. Smaller crystallographic units were seen in the sample compressed at 950 °C (Figure 10a) than in the samples at 1150 °C (Figure 10b). At 950 °C, less UC size distribution was observed, as shown in (Figure 11a) in comparison to the higher forming temperatures (Figure 11b), which showed a more extensive UC distribution. Also, a morphological variation, where more granular packages were seen at lower temperatures whereas, at higher temperatures, large packages with laths were observed.

Figure 10
EBSD map showing the temperature dependence in different sizes of crystallographic units. (a) P3 - T = 950 °C / ε˙= 0.1 s-1, (b)P3 - T = 1150 °C / ε˙= 0.1 s-1.
Figure 11
Histogram of measured crystallographic units grain sizes, (a) T = 950 °C / ε˙= 0.1 s-1, (b) T = 1150 °C / ε˙= 0.1 s-1.

The smaller austenitic grain sizes were assumed to be responsible for the smaller crystallographic UC, which contribute to the increase in high-angle grain boundaries (HAGB), which, according to the literature2626 Li X, Ma X, Subramanian SV, Shang C. EBSD characterization of secondary microcracks in the heat affected zone of a X100 pipeline steel weld joint. Int J Fract. 2015;193(2):131-9.,2727 Białobrzeska B, Konat Ł, Jasiński R. The influence of austenite grain size on the mechanical properties of low-alloy steel with boron. Metals (Basel). 2017;7(1):1-20., improves the material's toughness and strength.

3.3. Industrial application

3.3.1. Hot forging simulation results

A point-tracked section of DIN 18MnCrSiMo6-4 helical gear was analyzed in the simulation results. The average temperature of the tracked points displayed in (Figure 4), representing each analyzed region (R1, R2, and R3) in the simulation, is shown in Figure 12.

Figure 12
Temperature evolution of forging and cooling simulation in the forged preform; (a) Austenitization at 1030 °C; (b) Austenitization at 1130 °C.

The temperature rises due to the energy dissipation in the form of heat (adiabatic heating) caused by strain rate and plastic strain, showed in (Figure 13a) and (Figure 13b). The simulation results show that the energy dissipation is higher during forging with the initial temperature of 1030 °C than 1130 °C due to mechanical energy available to be converted into heat being higher in lower temperatures. There was an increase of 50 °C from the temperature of 1030 °C and 30 °C at the temperature of 1130 °C.

Figure 13
Simulated evolution of (a) average strain rate (ε˙) and (b) effective strain (ε¯) on the industrial preform forged gear.

The conduction between the workpiece and dies provokes a decrease in temperature on both surfaces, which is more pronounced for the lower die, as shown after a dwelling of 2.023 seconds. Moreover, the heat flux from region R3 in direction to region R2 and R2 to R1 maintains R2's temperature constant for a longer time and, therefore, causes an increase in the austenitic grain size by grain coarsening. This behavior is shown in (Figure 12a) for forging at 1030 °C and (Figure 12b) for forging at 1130 °C.

Lower temperatures also increase the force required for the plastic flow and, therefore, the system's energy1919 Lin J, Balint D, Pietrzyk M. Microstructure evolution in metal forming processes: modelling and applications. Sawston, Reino Unido: Woodhead Publishing; 2012. p. 67-74.. Even with the increase in energy dissipation in higher strain rates, forging at lower temperatures (1030 °C) does not exceed the value of 1100 °C-in this way, avoiding the excessive grain coarsening.

Although these behaviors are expected, it is only through FEM that it can accurately quantify the local distribution of strain throughout the billet, allowing evaluation of mechanical properties after the tensile and impact tests. The numerical simulation showed a significant difference between the effective strain rate fields in the R2 region and R3 caused by the upper die piercing in the billet, changing the plastic flow rate.

Even though the increase in strain rate should reduce the recrystallization fraction due to the increase in critical strain for XDRX1616 Sellars CM, Whiteman JA. Recrystallization and grain growth in hot rolling. Met Sci. 1978;13(3-4):187-94. https://doi.org/10.1179/msc.1979.13.3-4.187.
https://doi.org/10.1179/msc.1979.13.3-4....
, we must take into account the adiabatic heating in high-speed closed die forging, which interacts in two ways: First, through the heterogeneous plastic strain, which was higher in R1 than in R2 and R3, and second, through the heterogeneous strain rate, which changes the parameter (Z) responsible for the magnitude of DDRX. Thus, the simultaneous increase in plastic strain and strain rate in R1 resulted in a higher XDRX fraction due to the temperature increase, which delayed the time of grain growth by increasing the dynamically recrystallized grains DDRX.

Different mechanisms such as XDRX, MDRX, and SRX can be qualitatively analyzed through FEM by the JMAK macro-scale model results. Figure 14 provides the results obtained from the numerical simulation of these phenomena, and, Figure 15 shows the grain size evolution. According to these results, when (ε¯>εc), there is a predominance of XDRX (Figure 14a), which is proportional to the plastic strain in each selected region.

Figure 14
Recrystallization kinetics of the modeled preform hot forged; (a) XDRX; (b) MDRX; (c) SRX.
Figure 15
X-Y Plane view of the grain size evolution computed on the tracked points: (a) Forged at 1030 °C, (b) Forged at 1130 °C.

Nonetheless, compared to the analyzed regions R1, R2 and R3, the strain rate may also be responsible for the lower XDRX values found in R2 and R3. It is well known that strain rate affects the dynamic recrystallization kinetics due to the increase in dislocation density1616 Sellars CM, Whiteman JA. Recrystallization and grain growth in hot rolling. Met Sci. 1978;13(3-4):187-94. https://doi.org/10.1179/msc.1979.13.3-4.187.
https://doi.org/10.1179/msc.1979.13.3-4....
. Higher strain rates may decrease the mobility of these dislocations in R2 and R3, which consequently reduces annihilation or rearrangement. In contrast, in the R1 region, which may increase mobility due to higher temperature, over 95% of complete recrystallization was achieved.

The regions that were submitted to minor plastic strain and lower than critical strain εc was recrystallized Metadynamically (MDRX) and statically (SRX), respectively, as shown in Figures 14b and 14c. An apparent phenomenon is that the MDRX is slower than the DRX's kinetics and is observed in R1. Therefore, the significant difference between each recrystallization phenomenon is that the final nucleated grain size may be larger than the XDRX at the same deformation time as holding. The main reason for this phenomenon is the lack of dislocation multiplication during the MDRX process2828 Cho SH, Kang KB, Jonas JJ. Dynamic, static and metadynamic recrystallization of a Nb-microalloyed steel. ISIJ Int. 2001;41(1):63-9..

In the case of closed-die forging, such a mechanism could be compensated by the temperature rise caused by the stored energy due to the plastic strain, strain rate, and the heat dissipation itself. The higher forging temperature (1130 °C) and strain rate may provide more thermal energy for grain boundary migration2929 Jiang H, Dong J, Zhang M, Yao Z. A study on the effect of strain rate on the dynamic recrystallization mechanism of alloy 617B. Metall Mater Trans, A Phys Metall Mater Sci. 2016;47(10):5071-87.. In R2, the temperature was higher than in R1 after 8.23 s, as shown in Figure 12b, and therefore more thermally activated for the grain growth mechanism (Figure 15b).

3.3.2. Grain growth simulation after hot forging

Based on the results in Figure 10 and the discussion in the previous section, we have considered the grain growth exponent m (Equation 20) as a variable, which depends on the grain size distribution1818 Janssens KGF, Raabe D, Kozeschnik E, Miodownik MA, Nestler B. Computacional materials engineering: an introduction to microstructure evolution. Cambridge: Academic Press; 2007. 359 p., plastic deformation degree, and temperatures.

Hodgson and Gibbs3030 Hodgson PD, Gibbs RK. A Mathematical model to predict hot rolled C-Mn and microalloyed the mechanical properties of steels. ISIJ Int. 1992;32(12):329-38. also stated that high values of m are strongly dependent on annealing temperatures and a weak function of the initial recrystallized grain size and holding time. However, in the case of short holding times in closed die forging, the exerted heterogeneous plastic strain distribution, strain rate, and adiabatic heating result in different grain size distribution and recrystallization mechanisms (Figure 11, Figure 12, Figure 13, and Figure 14). Therefore, a constant value of the growth exponent is insufficient to predict the austenite grain growth. Such consideration about this behavior can also be associated with the high-temperature peak condition caused by the adiabatic heating, in which all fine precipitates are dissolved3131 Banerjee K, Militzer M, Perez M, Wang X. Nonisothermal austenite grain growth kinetics in a microalloyed x80 linepipe steel. Metall Mater Trans, A Phys Metall Mater Sci. 2010;41(12):3161-72..

To distinguish between these parameters with Hillert's theory for ideal growth mechanism3232 Hillert M. On the theory of normal and abnormal grain growth. Acta Metall. 1965;13:227., we assumed that above 1100 °C, exponent m > 4 represents the increase of the growth rate, while m between 3-4 would represent a decrease in growth rate. This makes the magnitude of the collision factor of the atoms decrease in terms of the reactions based on Arrhenius's theory. The modified empirical grain growth exponent values are shown in Table 5.

Table 5
Values of m exponent implemented in DEFORM® computed as a temperature function in non-isothermal conditions after DRX, MDRX, and SRX.

As shown in Figure 15, in the first 2 seconds, most of the heat transfer is between the billet's surface and the lower and upper dies. The pressure of the mechanical system also raises, increasing the temperature in regions R1 and R2 and, therefore, may increase the growth rate for temperatures above 1100 °C.

Figure 16 compares the predicted grain size results with m as a constant value and m as a function of the temperature in the different regions. Modifying the growth exponent as a temperature function improved the model compared to the constant m exponent, closer to the experimental results. The temperature was mainly responsible for the austenitic grain size variation when comparing both cases. However, for the forging done at 1130 °C, the local strain also significantly influenced the grain size variation.

Figure 16
Local austenitic grain size after industrial hot forging comparing experimental and numerical simulation results.

Silveira et al.88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
showed that in a hot forging experiment above 1100 °C at a low strain rate of about 0.17 s-1, the austenitic grain size remained large due to the extended periods the continuous cooling required. In contrast, this work employed high strain rates reaching up to 50 s-1 showing the influence of this parameter when combined with high plastic strain, which resulted in grain refinement when forging above 1100 °C due to adiabatic heating. This can be seen by comparing the R1 and R3 regions. The increase in the strain rate delays the onset of dynamic recrystallization kinetics. i.e., the critical strain εc is the trigger for this reaction; therefore, this variation can be explained by the local strain rate gradient.

3.3.3. Mechanical properties response

Experimental results of austenitic grain size were used to validate the accuracy of the numerical simulation. It is essential to remind that, in this present study, the bainitic phase transformation was not modeled. Hence, only the grain size model during hot forging and its influence on mechanical properties is addressed.

Figure 17 shows the experimental austenitic grains in the R2 region and the final microstructure after quenching the forged part at temperatures of 1030 °C (Figure 17a) and 1130 °C (Figure 17b), respectively.

Figure 17
Austenitic grain and continuously cooled microstructure in the R2 region: (a,b) and (c) Forged at 1030 °C, quenched and cooled in calm air; (d,e) and (f) forged at 1130 °C quenched and cooled in calm air.

As already mentioned, the dislocation motion is influenced by temperature variation, which is a lower temperature likely decreases the growth rate of the dynamically recrystallized grains. As can be seen, more refined austenitic grain forms when forged at 1030 °C, which shows the importance of controlling the forging temperature of this steel. This grain size variation resulted in different bainitic morphologies in the microstructure after continuous air cooling, as seen at lower magnification in Figures 17b and 17e forged at 1030 °C and 1130 °C, respectively. Visually, it is possible to observe a considerable amount of polygonal ferrite PF at the 1030 °C (Figure 17c) compared to 1130 °C (Figure 17f), which showed less PF fraction in this region and more Lath-like bainite LB.

Austenitic grain refinement increases the amount of ferrite nucleation sites due to the increase in the grain boundary area. The increase in the dislocation density caused by plastic deformation, temperature, and strain rate, are responsible for retained austenite in the microstructure88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
,3333 Aranda MM, Kim B, Rementeria R, Capdevila C, De Andrés CG. Effect of prior austenite grain size on pearlite transformation in A hypoeuctectoid Fe-C-Mn steel. Metall Mater Trans, A Phys Metall Mater Sci. 2014;45(4):1778-86.. The silicon content in this steel can also be responsible for more carbon enrichment in the retained austenite, which is higher in temperature above 1100 °C88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
. The increase of defects caused by plastic strain increases the number of possible nucleation sites due to the growth of A3 austenite-ferrite transformation88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
.

Even at high strain rates, the metallurgical evolution of the ferritic bainite transformation was the same as reported by66 Lemos Bevilaqua W, Epp J, Meyer H, Da Silva Rocha A, Roelofs H. In situ investigation of the bainitic transformation from deformed austenite during continuous cooling in a low carbon Mn-Si-Cr-Mo steel. Metall Mater Trans, A Phys Metall Mater Sci. 2020;51(7):3627-37.

7 Hatwig RA, Dong J, Epp J, Rocha A S. Effect of compressive deformations on the final microstructure of a low carbon high silicon bainitic steel thermomechanically processed. Mater Res. 2021;24(1):1-7.
-88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
. The temperature significantly influenced the transformations. By forging at 1130 °C, a more significant fraction of the Lath-like bainite can be observed. At 1030 °C, granular bainite rises more.

Finally, the mechanical properties are analyzed. Figure 18a shows the boxplot of the yield YS, and Figure 18b shows the tensile strength US results from the miniaturized tensile tests for the three selected regions of the industrial hot forged component. Figure 19a shows the toughness response of the material to the Charpy impact tests, whereas the elongation (EL) is shown in Figure19b.

Figure 18
Boxplot results of the (a) Yield stress YS and (b) Ultimate stress US in MPa.
Figure 19
(a) Impact energy through Charpy- V notch tests; (b) Local Elongation in %.

The results suggest that the forging temperature was mainly responsible for the increase in the YS (Figure 18a) and US (Figure 18b), where the calculated average stress ratios (YSavg/TUavg) were 0.70 at 1030 °C and 0.68 at 1130 °C. However, by taking the R2 region into account, the plastic strain rate at low temperature also exerted a considerable effect on the YS and TS. This result can be explained through the numerical simulation, which showed that, by the time of forging unloading, the energy accumulated in the R1 region was greater than R2, which accelerated grain growth due to adiabatic heating, as shown in 2 and Figure 15, respectively.

Figure 20 shows the microstructure from the micro-tensile samples on a magnified scale. The different bainitic morphology revealed by scanning electron microscopy (SEM) in the three regions is shown.

Figure 20
SEM results of the micro-tensile specimen regions show the final continuous cooled microstructure of preform gear. (a) R1, (b) R2 and (c) R3 Forged at 1030 °C; (d) R1, (e) R2 and (f) R3 Forged at 1130 °C, respectively.

The physical simulations performed by66 Lemos Bevilaqua W, Epp J, Meyer H, Da Silva Rocha A, Roelofs H. In situ investigation of the bainitic transformation from deformed austenite during continuous cooling in a low carbon Mn-Si-Cr-Mo steel. Metall Mater Trans, A Phys Metall Mater Sci. 2020;51(7):3627-37.,77 Hatwig RA, Dong J, Epp J, Rocha A S. Effect of compressive deformations on the final microstructure of a low carbon high silicon bainitic steel thermomechanically processed. Mater Res. 2021;24(1):1-7. with controlled temperature in continuous cooling rate and low strain rate88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
have shown a similar bainitic structure66 Lemos Bevilaqua W, Epp J, Meyer H, Da Silva Rocha A, Roelofs H. In situ investigation of the bainitic transformation from deformed austenite during continuous cooling in a low carbon Mn-Si-Cr-Mo steel. Metall Mater Trans, A Phys Metall Mater Sci. 2020;51(7):3627-37.

7 Hatwig RA, Dong J, Epp J, Rocha A S. Effect of compressive deformations on the final microstructure of a low carbon high silicon bainitic steel thermomechanically processed. Mater Res. 2021;24(1):1-7.
-88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
. However, in this work, it has been shown that the influence of the temperature variation in high strain rate on the microstructure evolution in an industrial hot forged condition. Figure 20a (R1), Figure 20b (R2), and Figure 20c (R3), forged at 1030 °C, revealed different lengths of grain boundaries compared to the samples forged at 1130 °C (Figure 20d, Figure 20e, and Figure 20f). Finite element analysis showed the local temperature, plastic strain, and strain rate variation for each region of both hot forged conditions. These results encouraged us to interpret that grain boundary changes have been influenced by the synergetic effect in the closed die forging parameters. The temperature variation due to the energy dissipation through plastic strain and strain rate has possibly made the austenite grains in R1 and R2 smaller than R3. At 1030 °C, the grains are more homogeneous and refined than the 1130 °C forging temperature.

The toughness of the component forged at 1030 °C increased compared to the material received and forged at 1130 °C. This result is shown in Figure 19a and confirms the decrease in crystallographic unity (UC), previously shown by the EBSD in Figure 10a. The UC decrease provides a less heterogeneous microstructure package, which could significantly influence the toughness increase of the bainitic steel, as seen in Figure 11a.

The hot forging at 1030 °C promoted a greater fraction of polygonal ferrite, which, according to88 Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601.
https://doi.org/10.3390/met10050601...
, resulted in a chemical heterogeneity in the remaining austenite. This phenomenon causes the formation of a higher fraction of martensite in regions with lower concentrations of carbon and alloying elements3434 Lambert A, Lambert A, Drillet J, Gourgues AF, Sturel T, Pineau A. Microstructure of martensite–austenite constituents in heat affected zones of high strength low alloy steel welds in relation to toughness properties. Sci Technol Weld Join. 2000;5(3):168-73. http://dx.doi.org/10.1179/136217100101538164.
http://dx.doi.org/10.1179/13621710010153...
. This composition heterogeneity is more predominant in R2, which could have resulted in a higher fraction of martensite, which would lead to an increase in strength, as can be seen in Figure 18b and Figure 20b.

4. Conclusion

This work employed the semi-empirical JMAK models to simulate the recrystallization and grain growth after closed die forging for new continuously cooled bainitic steels. The influence of hot forging temperatures on local tensile strength and ductility (toughness) by Charpy V-notch tests was also verified. The following conclusions could be drawn:

  • Grain growth exponent, which depends on activation energy, geometric and frequency factors based on the Arrhenius relationship, can vary during the heating and forging steps. Due to the variation of the recrystallized grain size distribution (Figure 11), different exponent values (m) were determined according to temperature. Therefore, after recrystallization, the grain growth results were closer considering the temperature variation.

  • From the numerical simulation results, adiabatic heating, caused by high strain rates and plastic deformation, may be responsible for the rapid grain growth in forging. It was also observed that in hot forging below 1130 °C, the austenitic grain size was refined, and abnormal grain growth was absent. This led us to believe that there is no influence on the variation of the exponent m at temperatures below 1100 °C.

  • The hot forging results below 1130 °C led to higher austenitic grain refinement and no abnormal grain growth. This confirmed previous results in forging tests using a hydraulic press.

  • The austenitic grain size favored the formation of different morphologies of bainite at forging temperatures below 1100 °C.

  • The localized variations of the bainitic microstructure morphology resulted in different mechanical properties. The temperature of 1030°C in R2 achieved the highest yield strength and the highest ratio between tensile strength and maximum stress, while the toughness could be improved compared to forging at 1130 °C.

  • These results showed that, even with a temperature variation of around 100 °C, the control of the temperature forging has a significant impact on the mechanical properties. The numerical simulation was an efficient tool to determine and optimize the austenitic grain evolution.

5. Acknowledgments

The authors acknowledge the Brazilian agency CAPES (Project 1844/2017) and the Deutsch Forschungsgemeinschaft (German Research Foundation – DFG) via project EP 128/6-2 (Project number 327887503) for the financial support within the BRAGECRIM program (Brazilian German Cooperation Research in Manufacturing). The authors also thank the INPEL S.A industry for the research partnership. To Metal Forming Institute IBF- RWTH Aachen in Germany for the support during experimental data acquisition. Thiago M. Ivaniski and Alexandre da S. Rocha also acknowledge CNPq (process numbers 167948/2017-2 and 308773/2018-7) for the grants.

6. References

  • 1
    Shimanov M, Korpala G, Terzic A, Kawalla R. Bainitic steels, their characteristics and applications. KEM. 2016;684:104-10. https://doi.org/10.4028/www.scientific.net/KEM.684.104
    » https://doi.org/10.4028/www.scientific.net/KEM.684.104
  • 2
    Caballero FG, Capdevila C, Chao J. The microstructure of continuously cooled tough bainitic steels. In: 2nd International Conference of Super-High Strength Steels; 2010 Oct 27-30; Verona, Italy. Proceedings. Verona, Italy: TIB; 2010. p. 10-8.
  • 3
    Roelofs H, Hasler S, Urlau U, Lembke MI, Olschewski G. Continuously cooled bainitic steel HSX ® Z12: one decade of experience. In: 4th International Conference on Steels in Cars and Trucks; 2014, Jun 15-19; Braunschweig, Germany. Proceedings. Braunschweig, Germany.
  • 4
    Caballero FG, Garcia-Mateo C, Cornide J, Allain S, Puerta J, Crouvizer M, et al. New advanced ultra high strength bainitic steels: ductility and formability (DUCTAFORM) [Internet]. Brussels: European Commission; 2013 [cited 2021 Nov 24]. Available from: https://op.europa.eu/en/publication-detail/-/publication/9f8dc6ba-d882-48da-9b7b-aa8d7f00dd24
    » https://op.europa.eu/en/publication-detail/-/publication/9f8dc6ba-d882-48da-9b7b-aa8d7f00dd24
  • 5
    Lan HF, Du LX, Li Q, Qiu CL, Li JP, Misra RDK. Improvement of strength-toughness combination in austempered low carbon bainitic steel: the key role of refining prior austenite grain size. J Alloys Compd. 2017;710:702-10.
  • 6
    Lemos Bevilaqua W, Epp J, Meyer H, Da Silva Rocha A, Roelofs H. In situ investigation of the bainitic transformation from deformed austenite during continuous cooling in a low carbon Mn-Si-Cr-Mo steel. Metall Mater Trans, A Phys Metall Mater Sci. 2020;51(7):3627-37.
  • 7
    Hatwig RA, Dong J, Epp J, Rocha A S. Effect of compressive deformations on the final microstructure of a low carbon high silicon bainitic steel thermomechanically processed. Mater Res. 2021;24(1):1-7.
  • 8
    Silveira ACF, Bevilaqua WL, Dias VW, Castro PJ, Epp J, Rocha AS. Influence of hot forging parameters on a low carbon continuous cooling bainitic steel microstructure. Metals (Basel). 2020;10(5):601. https://doi.org/10.3390/met10050601
    » https://doi.org/10.3390/met10050601
  • 9
    Ivaniski TM, de Souza TF, Épp J, Rocha AS. Constitutive modelling of high temperature flow behaviour for a low carbon high silicon bainitic steel. Mater Res. 2020;23(5):e20200264.
  • 10
    Jonas JJ, Quelennec X, Jiang L, Martin É. The Avrami kinetics of dynamic recrystallization. Acta Mater. 2009;57(9):2748-56.
  • 11
    Irani M, Lim S, Joun M. Experimental and numerical study on the temperature sensitivity of the dynamic recrystallization activation energy and strain rate exponent in the JMAK model. J Mater Res Technol. 2019;8(2):1616-27. https://doi.org/10.1016/j.jmrt.2018.11.007
    » https://doi.org/10.1016/j.jmrt.2018.11.007
  • 12
    Razali MK, Joun M. A new approach of predicting dynamic recrystallization using directly a flow stress model and its application to medium Mn steel. J Mater Res Technol. 2021;11(18):81-94. http://dx.doi.org/10.1016/j.jmrt.2021.02.026
    » http://dx.doi.org/10.1016/j.jmrt.2021.02.026
  • 13
    Chamanfar A, Chentouf SM, Jahazi M, Lapierre-Boire L-P. Austenite grain growth and hot deformation behavior in a medium carbon low alloy steel. J Mater Res Technol. 2020;9(6):12102-12114. https://doi.org/10.1016/j.jmrt.2020.08.114
    » https://doi.org/10.1016/j.jmrt.2020.08.114
  • 14
    Moraes ALI, Balancin O. Numerical simulation of hot closed die forging of a low carbon steel coupled with microstructure evolution. Mater Res. 2015;18(1):92-7.
  • 15
    Ivaniski TM, Epp J, Zoch H-W, Da Silva Rocha A. Austenitic grain size prediction in hot forging of a 20mncr5 steel by numerical simulation using the JMAK model for industrial applications. Mater Res. 2019;22(5)
  • 16
    Sellars CM, Whiteman JA. Recrystallization and grain growth in hot rolling. Met Sci. 1978;13(3-4):187-94. https://doi.org/10.1179/msc.1979.13.3-4.187
    » https://doi.org/10.1179/msc.1979.13.3-4.187
  • 17
    Henke T, Hirt G, Bambach M. Optimization of a closed die forging process to manufacture a gear wheel by the use of a response surface model. Adv Mat Res. 2014;922:254-9.
  • 18
    Janssens KGF, Raabe D, Kozeschnik E, Miodownik MA, Nestler B. Computacional materials engineering: an introduction to microstructure evolution. Cambridge: Academic Press; 2007. 359 p.
  • 19
    Lin J, Balint D, Pietrzyk M. Microstructure evolution in metal forming processes: modelling and applications. Sawston, Reino Unido: Woodhead Publishing; 2012. p. 67-74.
  • 20
    ASTM: American Society for Testing and Materials. ASTM E112. West Conshohocken: ASTM; 2014.
  • 21
    Castro PJ, Silveira AC F, Ivaniski TM, Turra CJ, Epp J, Rocha A S. Two-step continuous cooling heat treatment applied in a low carbon bainitic steel. Mater Res. 2021;24(2)
  • 22
    Kashaev N, Horstmann M, Ventzke V, Riekehr S, Huber N. Comparative study of mechanical properties using standard and micro-specimens of base materials Inconel 625, Inconel 718 and Ti-6Al-4 v. J Mater Res Technol. 2013;2(1):43-7.
  • 23
    ASTM: American Society for Testing and Materials. ASTM E2298. West Conshohocken: ASTM; 2018
  • 24
    Vasilyev AA, Sokolov SF, Kolbasnikov NG, Sokolov DF. Effect of alloying on the self-diffusion activation energy in γ-iron. Phys Solid State. 2011;53(11):2194-200.
  • 25
    Guo-Zheng Q, Yang W, Ying-Ying L, Jie Z. Effect of temperatures and strain rates on the average size of grains refined by dynamic recrystallization for as-extruded 42CrMo steel. Mater Res. 2013;16(5):1092-105.
  • 26
    Li X, Ma X, Subramanian SV, Shang C. EBSD characterization of secondary microcracks in the heat affected zone of a X100 pipeline steel weld joint. Int J Fract. 2015;193(2):131-9.
  • 27
    Białobrzeska B, Konat Ł, Jasiński R. The influence of austenite grain size on the mechanical properties of low-alloy steel with boron. Metals (Basel). 2017;7(1):1-20.
  • 28
    Cho SH, Kang KB, Jonas JJ. Dynamic, static and metadynamic recrystallization of a Nb-microalloyed steel. ISIJ Int. 2001;41(1):63-9.
  • 29
    Jiang H, Dong J, Zhang M, Yao Z. A study on the effect of strain rate on the dynamic recrystallization mechanism of alloy 617B. Metall Mater Trans, A Phys Metall Mater Sci. 2016;47(10):5071-87.
  • 30
    Hodgson PD, Gibbs RK. A Mathematical model to predict hot rolled C-Mn and microalloyed the mechanical properties of steels. ISIJ Int. 1992;32(12):329-38.
  • 31
    Banerjee K, Militzer M, Perez M, Wang X. Nonisothermal austenite grain growth kinetics in a microalloyed x80 linepipe steel. Metall Mater Trans, A Phys Metall Mater Sci. 2010;41(12):3161-72.
  • 32
    Hillert M. On the theory of normal and abnormal grain growth. Acta Metall. 1965;13:227.
  • 33
    Aranda MM, Kim B, Rementeria R, Capdevila C, De Andrés CG. Effect of prior austenite grain size on pearlite transformation in A hypoeuctectoid Fe-C-Mn steel. Metall Mater Trans, A Phys Metall Mater Sci. 2014;45(4):1778-86.
  • 34
    Lambert A, Lambert A, Drillet J, Gourgues AF, Sturel T, Pineau A. Microstructure of martensite–austenite constituents in heat affected zones of high strength low alloy steel welds in relation to toughness properties. Sci Technol Weld Join. 2000;5(3):168-73. http://dx.doi.org/10.1179/136217100101538164
    » http://dx.doi.org/10.1179/136217100101538164

Publication Dates

  • Publication in this collection
    01 July 2022
  • Date of issue
    2022

History

  • Received
    24 Nov 2021
  • Reviewed
    28 Apr 2022
  • Accepted
    08 June 2022
ABM, ABC, ABPol UFSCar - Dep. de Engenharia de Materiais, Rod. Washington Luiz, km 235, 13565-905 - São Carlos - SP- Brasil. Tel (55 16) 3351-9487 - São Carlos - SP - Brazil
E-mail: pessan@ufscar.br