Physical and Numerical Simulations of Closed Die Hot Forging and Heat Treatment of Forged Parts

The paper describes physical and numerical simulations of a manufacturing process composed of hot forging and controlled cooling, which replace the conventional heat treatment technology. The objective was to investigate possibilities and limitations of the heat treatment with the use of the heat of forging. Three steels used to manufacture automotive parts were investigated. Experiments were composed of two sets of tests. The first were isothermal (TTT) and constant cooling rate (CCT) dilatometric tests, which supplied data for the identification of the numerical phase transformation model. The second was a physical simulation of the sequence forging-cooling on Gleeble 3800, which supplied data for the validation of the models. In the numerical part, a finite element (FE) thermal-mechanical code was combined with metallurgical models describing recrystallization and grain growth during forging and phase transformations during cooling. The FE model predicted distributions of the temperature and the austenite grain size in the forging, which were input data for further simulations of phase transformations during cooling. A modified JMAK equation was used to calculate the kinetics of transformation and volume fraction of microstructural constituents after cooling. Since the dilatometric tests were performed for various austenitization temperatures before cooling, it was possible to include austenite grain size as a variable in the model. An inverse algorithm developed by the authors was applied in the identification procedure. The model with optimal material parameters was used for simulations of hot forging and controlled cooling in one of the forging shops in Poland. Distributions of microstructural constituents in the forging after cooling were calculated. As a consequence, various cooling sequences during heat treatment could be analyzed and compared.


Introduction
Problems of the heat treatment of hot forged parts have been intensively studied and numerous papers have been published; see, for example, [1][2][3] or a paper dedicated to automotive parts [4]. On the other hand, due to its constraints, the heat treatment of automotive parts with the use of the heat of forging is still a challenge and published material is less numerous. Advantages of this process are discussed in [5]. They can be summarized as follows: (i) hot forgings with outstanding mechanical properties can be produced, (ii) economic and ecological advantages can be achieved in contrast to conventional hot forging process chains, which consist of various heating and cooling cycles, and (iii) different forging properties and microstructures can be set up by an adjusted temperature control. The majority of publications dealing with this topic are focused on experiments. Some papers consider annealing before forging [6]. The authors of [7] provided general information on direct-forge quenching (DFQ) and direct heat treatment (DHT) processes that are used in automotive and other mechanical industries. Technological advances in these processes and their ability to produce high-quality components at a low production cost were discussed. A spray quenching unit has been designed in [7] for integrated heat treatment in the forging process chain for steel components. Several quenching strategies with varying process conditions have been simulated based on time-temperature-transformation (TTT) diagrams. The design of a quenching unit for the control of the microstructural evolution was based on the simulation results. Few papers have dealt with evaluation of the effects of hot forging and post-forging heat treatment on the mechanical properties of products [8]. Authors of [9] designed a spray quenching unit for integrated heat treatment in the forging process chain for steel components. The simulations helped the design of a quenching unit integrating a sensor device for the control of the microstructural evolution in the workpieces during the quenching. The possibility to perform homogeneous bainite transformation from the forging heat using a stepped quenching strategy (combining gas jets and sprays as well as isothermal heat treatment) was demonstrated for such types of parts as a stepped shaft and a common rail. Authors of [10] used a finite element (FE) model in a macro scale and a material model based on the additive strain decomposition in the micro scale.
In a majority of the discussed papers, constant cooling rate (CCT) and/or isothermal (TTT) diagrams were used to predict phase transformations. Exploring the possibility of more reliable and still fast modelling of phase transformations during complex cooling schedules after hot forging was the motivation of the present work. A review of phase transformation models [11] divides these models into two groups: mean field models and full field models. The latter are based on such methods as phase field, Cellular Automata or Monte Carlo. As long as the latter models account for the geometrical features of the material microstructure, the former describes metallurgical phenomena in a single material point using closed-form equations. Since the hot forging process involves significant heterogeneity of strains, stresses and temperatures, the finite element (FE) method is needed to describe macro scale phenomena. Accounting for the influence of these heterogeneities on gradients of properties in the final product requires a solution of the phase transformation model in each Gauss point of the FE mesh [12], which practically eliminates the full field models. Thus, the emphasis in the present work was put on the upgrade of the mean field models to extend their predictive capabilities while the computing costs remain very low.
The general objective of the present paper was to further investigate possibilities and limitations of heat treatment with the use of the heat of forging. This goal was reached by both numerical and physical simulations of the process. The particular focus was on development of a new, reliable and fast numerical multiscale model of the thermal, mechanical and metallurgical phenomena during hot forging and cooling. This was done by combining basic material tests with physical and numerical simulations of the manufacturing cycle composed of hot forging and controlled cooling.

FE Model
The finite element (FE) program Forge NxT 3.0 (version number) developed by Transvalor [13] was used in simulations of thermal and mechanical phenomena in the macro scale during forging and cooling. Meshing was calculated using the package's tool from the curvature of initial mesh. The parts were re-meshed between stages. However, since the geometry changes significantly during forging, periodic remeshing every 20 computational steps was applied to avoid numerical errors. To make computations faster, tool geometries were simulated as contact surfaces.
During forging, the time step was set by built-in routines, which adapt the step according to previous deformation or temperature results. However, since during the cooling simulation the recalescence heat is returned to the solver, a different time step adaptation scheme was used in the thermal simulation. A standard adaptation based on previous values was used until the recalescence heat became non-zero. Then, a reduced time step was introduced as a constant value set, used for the entire length of transformations.
The heat transfer coefficient for the air cooling was taken from program's database. For accelerated cooling on the cooling line, the coefficient was identified by comparing the simulated temperature profile with measurements. The material models supplied in the Forge ® package were used in mechanical (flow stress) and thermal (conductivity, specific heat, density, expansion coefficient) parts. The flow stress, which is the material parameter in the mechanical part, was calculated based on the Hansel-Spittel equation with coefficients, which were identified using inverse analysis for the hot compression tests results. The inverse algorithm described in [11] was used. Equations describing phase transformations during cooling were implemented in the user's procedures in the FE program.
The objectives of the FE simulations were twofold. The first was calculation of the strain, strain rate and temperature history in the selected locations in the forging. This history was reproduced in the physical simulation on the Gleeble 3800. Simulations of the distribution of the microstructural parameters in the whole volume of the forgings was the second objective.

Microstructure Evolution Model
Equations describing processes of recrystallization and grain growth during hot forming were based on the fundamental works of Sellars at the University of Sheffield [14], and these equations are in the Forge database. Coefficients in the equations for the investigated steels were determined in the present project on the basis of compression tests, stress relaxation tests and grain growth simulations (see chapter 4). All equations are given in Table 1, where: SRX-static recrystallization, DRX-dynamic recrystallization, t-time, ε-strain, . ε-strain rate, T-temperature in K, R-gas constant, D 0 -grain size prior to deformation, Z-Zener-Holomon parameter, Q SRX , Q DRX , Q DSRX , Q GROWTH -activation energy for the relevant process, n, a 0 −a 3 , b 0 −b 3 , p 1 −p 10 , q, K-coefficients determined by the inverse analysis of the experimental tests. Table 1. Equations describing microstructure evolution during hot forming.

Parameter Equation
Kinetics of SRX Grain size after SRX Critical strain for DRX DRX volume fraction Saturation strain Grain size after DRX The FE model with microstructure evolution equations implemented into the FE Forge code predicted distributions of the temperature and the grain size in the part after hot forging. These data were a starting point for further simulations of phase transformations during cooling.

Phase Transformation Model
The modified JMAK (Johnson-Mehl-Avrami-Kolmogorov) equation was used to calculate the kinetics of transformation and volume fraction of microstructural constituents after cooling. The main equation of this model is where: X-volume fraction of a new phase, t-time, k, n-coefficients.
The numerical implementation of this model is described in [15]. Since the dilatometric tests were performed for various austenitization temperatures before cooling, it was possible to account for the effect of the austenite grain size on the kinetics of transformations. The austenite grain size at the beginning of transformations (D γ ) was included as a variable in the model. The main equations in the model are given in Table 2. The following upgrades were introduced in the JMAK equation [16]: (1) Coefficient n is constant for each transformation and in the identification it is referred to as a 4 , a 16 and a 24 for ferrite, pearlite and bainite transformations, respectively. (2) According to [15], coefficient k is temperature dependent. The modified Gauss function proposed in [17] was used. (3) Calculations of carbon concentration in the austenite during both ferrite and bainite transformations were added. (4) The T 0 temperature concept proposed in [18] was added. It allowed prediction of the return of the pearlitic transformation after bainitic has started during holding at the constant temperature. Beyond this, prediction of the occurrence of the retained austenite became possible. Table 2. Equations in the phase transformation models.

Parameter Equation
Incubation time of perlite for the remaining transformation τ P = Coefficients k for ferrite Coefficients k for perlite Coefficients k for bainite Average carbon content in austenite Bainite start temperatures B s = a 20 (15) Martensite start temperatures Martensite volume fraction Equations used for the calculation of the current carbon content in the austenite (c γ ), temperatures of bainite start (B s ) and martensite start (M s ), as well as martensite volume fraction (F m ) are shown in Table 2, where F f , F p and F b represent volume fractions of ferrite, pearlite and bainite, respectively. Parameter p in this table is explained in [19] and the numerical solution of the present model is described in [20]. The phase transformation model allows the effect of the recalescence to be taken into account. The heat generated during transformation per unit volume is calculated as where: ρ-density, H-enthalpy of the phase transformation. The phase transformation model was implemented into the FE Forge software using a user subroutine. The model was compiled to library form and connected with the FE package using its user routines. This technique allows the single implementation to be used in various FE programs only by adapting interface characteristics to the specific application. Source data is obtained from computation, intermediate variables are stored in mesh and results are passed as typical software-specific fields to be post-processed or used in other routines. In Forge software, two additional routines were added to obtain model coupling: (1) A routine which prepares storage for the model's intermediate variables, initializes or obtains simulation parameters passed to the model and calls the model main routine from a library file. The results are similarly stored in mesh fields and can be used in further simulation steps. (2) A routine which converts model results, especially transformation energy, to values accepted by the solver's FE thermal coupling (in the analysed case, thermal power). This routine allows a full thermal coupling of transformation on an energy level to be obtained without using additional models to compensate the transformation with modifying material parameters and can be turned off to make computations faster, but ignores thermal coupling.
Coupling of the FE code with the phase transformation model allowed for feedback from the micro scale to macro scale. As a consequence, the kinetics of transformations was calculated for the current, local temperatures, and the heat generated in the microstructure due to recalescence, calculated from Equation (18), was used in the macro FE calculations of the temperature.

Methodology
Uniaxial compression tests supplied data for identification of the flow stress model, and stress relaxation tests supplied data for identification of the microstructure evolution model. All the tests were carried on Gleeble 3800. The cylindrical samples having dimensions Ø10 mm × 12 mm were compressed at different temperatures, strain and strain rates. The methodology of these tests is well known and is not described here.
Dilatometric experiments were performed with a DIL 805 A/D/T dilatometer (New Castle, DE, USA). The experiments were intended to supply isothermal (TTT) and constant cooling rate (CCT) phase transformations diagrams, which were the input data for identification of phase transformation models. Prior to the cooling stage of the experiment, the samples (Ø5 mm × 10 mm) were subject to the deformations reproducing the forging process.
The physical simulation of the forging-cooling sequence was carried out on the Gleeble 3800 (Poestenkill, NY, USA) using blocky samples having dimensions 15 mm × 20 mm x 35 mm, which were subject to three-stage compression. Strain, strain rate and temperature variations calculated using finite element (FE) code Forge were reproduced in the experiments. This history for the two representative points in each forging (one located in a massive part and the second in the thin part of the forging) was reproduced on the Gleeble 3800. The physical simulations supplied data for the validation of the models and were used to design the industrial cooling sequence.
The microstructural observations were performed with light optical-digital microscopy (LOM) using an Olympus DSX500i (Tokyo, Japan) and with scanning electron microscopy (SEM) using INSPECT F (Tokyo, Japan). Quantitative analysis of grain size and volume fractions was carried out using MetiIo software. Quantitative evaluation of grain size was done on 10 micrographs recorded with 1000 times magnification, using an LOM microscope. In turn, the quantitative analysis of volume fractions of structural components was carried out on 20 micrographs recorded with 2000 times magnification, using the SEM microscope.

Materials
Manufacturing of three products, subject to different heat treatment involving controlled cooling after hot forging, was considered: (1) adapter made of 42CrMo4 steel (steel A in Table 3) subject to isothermal annealing, (2) flange made of the C22.8 steel (steel B in Table 3) subject to normalization, (3) fork made of the C45 steel (steel C in Table 3) subject to quenching and tempering. These steel grades were chosen in the present investigation due to their diversified behaviour during forging and subsequent cooling.

Deformation and Thermal Cycles
Dilatometric tests were performed for various austenitization parameters; as a consequence, the effect of the austenite grain size after forging on the kinetics of transformation could be included in the model. The cooling rates varied between 0.1 • C/s and 100 • C/s, which covered different investigated controlled cooling processes.

Physical Simulations
Physical simulations covered the whole process including three step forging and controlled cooling. As has been mentioned, two characteristic points in each forging were selected. One point (P) was located in the massive part of the forging and the second (Q) in the thinnest part. Schematic illustration of the investigated forgings with the locations of sensors is shown in Figure 1. Strain and temperature history, calculated with the Forge program, in these points during forging and cooling were reproduced on the Gleeble 3800 thermomechanical simulator. The samples were water quenched at various stages of the cooling process. The microstructure of the samples was analysed after each test and volume fractions of phases were estimated. Time-deformation-temperature data for the physical simulation are given in Table 4. The temperatures in this table refer to the beginning of the deformation. The strain rate varied during deformations and the maximum value was about 25 s −1 . Variations of the temperature versus time during cooling of the samples subject to the three-stage deformation on Gleeble 3800 are shown schematically in Figure 2. The following three variants of cooling before holding were considered: (1) Variant 0-cooling to the temperature above A c3 , heating to the holding temperature (all steels). (2) Variant 1-cooling until the austenite decomposition is completed, heating to the holding temperature (steels B and C).  The objective was to investigate the effect of the microstructure at the beginning of holding on the phase transformations.
Holding time for steel A varied between 1800 s and 3600 s. The objective was to investigate progress of pearlitic and bainitic transformations during cooling. The holding time for steels B and C was 1 h.

Results for Steel A
Grain size at the beginning of transformations was at the level of 50 ± 4 µm for the massive part of the forging and 40 ± 3 µm for the thinner part. Due to the small difference, all the results presented below are for the massive part only. In the tests with the holding temperature above 600 • C, pearlitic microstructures, with some ferrite at the level of around 5%, were obtained. Only one example of the microstructure of the pearlite for the holding temperature 620 • C is presented in Figure 3, and these results are not reproduced further in the paper. A more complex situation was observed at lower temperatures, in which bainitic transformation was dominating but it was controlled by the T 0 line. Microstructures of the samples after holding at temperatures of 560 • C and 520 • C are shown in Figures 4 and 5, respectively. Analysis of the pictures for 520 • C allows suggestion that pearlitic and bainitic transformations could occur alternately, however, the bainitic transformation occurs first and proceeds the pearlitic one until the moment when carbide content in untransformed austenite reaches the T 0 line. The possibility of modelling of this phenomenon is described in [19], and this approach was applied in the present work. Microstructures show also that even for very long holding times the austenite decomposition is not completed and the remaining austenite is transferred into bainite or martensite during subsequent cooling. Volume fractions of phases were determined for each sample and these results are presented in Section 5.1, where they are compared with the model predictions.

Results for Steel B
Physical simulations involved three step deformation followed by cooling according to schedules shown in Figure 2b. The microstructures in point P before the first deformation and after the last deformation are shown in Figure 6. In this case, the samples were cooled with water prior to and after deformations. The microstructure in that case contains thin ferrite layers decorating austenite grain boundaries and martensite + bainite mixture inside grains. The measured austenite grain size was 90 ± 7µm and 50 ± 4 µm, respectively. The grain size was similar in both points P and Q.
Three cooling sequences were simulated. The difference was in the lowest temperature after the first cooling following forging. This temperature was 344 • C, 697 • C and 726 • C for point P and 344 • C, 662 • C and 697 • C for point Q, for variants 1, 2 and 0, respectively. The temperatures were selected to obtain complete austenite decomposition (variant 1), only ferritic transformation without pearlite (variant 2) and less than 20% of the ferritic transformation (variant 0) in point P. The temperatures in point Q were slightly lower. The holding temperature was 920 • C for all variants. Cooling in air was applied after the holding, and the cooling rate varied between 3 • C/s at the beginning of cooling and below 1 • C/s at the end of cooling. Selected microstructures after cooling are shown in Figure 7. Quantitative analysis of the microstructural constituents is presented in Section 4.

Results for Steel C
Physical simulations involved four step deformation followed by cooling according to the schedules shown in Figure 2c. The microstructure in point P before the first deformation is shown in Figure 8. As in the case of A steel, fast cooling was applied to cause the precipitation of ferrite at the austenite grain boundaries. The grain size was similar in both points. It was non-uniform, with an average value of 100 ± 12 µm.     The results for cooling are presented for the holding temperature of 860 • C. The first set of results concerns point P in Figure 1c. Grain size during cooling according to the variant 0 was 70 ± 3 µm at the beginning of holding and 74 ± 4 µm at the end of holding. Microstructures during cooling according to variants 1 and 2 are shown in Figures 9 and 10. Microstructure after the first cooling to 540 • C (variant 1) is shown in Figure 9a. It is composed of about 12% of ferrite and pearlite. Ferrite and pearlite transferred into austenite during heating and fully austenitic microstructures at the beginning and end of holding were obtained. During further cooling the austenite grain boundaries were decorated by ferrite grains, which is seen in Figure 9b,c for the beginning and end of holding, respectively. The measured austenite grain size was 39 ± 2 µm in Figure 9b and 45 ± 2 µm in Figure 9c. The microstructure in point P after the fast cooling to 670 • C is shown in Figure 10a. It is composed of about 5% of pearlite and ferrite. Fully austenitic microstructures at the beginning and end of holding are shown in Figure 10b,c, respectively. The measured austenite grain size was 63 ± 3 µm and 85 ± 5 µm.
Similar results were obtained for point Q in Figure 1c. Grain size during cooling according to the variant 0 was 83 ± 5 µm at the beginning of holding and 87 ± 4 µm at the end of holding. The differences between austenite grain size after reheating and holding are not significant, which means that the grain growth slowed down after reaching the holding temperature. Microstructures during cooling according to variants 1 and 2 are shown in Figures 11 and 12. The microstructure after the first cooling to 500 • C (variant 1) is shown in Figure 11a. It is composed of ferrite (about 14%). Fully austenitic microstructures at the beginning and end of holding are shown in Figure 11b,c, respectively. The measured austenite grain size was 54 ± 3 µm and 61 ± 5 µm.   The microstructure in point Q after the first cooling to 609 • C followed by fast cooling in water is shown in Figure 12a. It is composed of about 40% of pearlite, 10% of ferrite and martensite/bainite, which was formed during water cooling. Microstructures at the beginning and end of holding are shown in Figure 12b,c, respectively. It is seen that pearlite transformation continued during heating prior to holding. On the contrary, the rate of cooling after holding was too small to prevent against the pearlite formation. Despite this, austenite grain size could be estimated by the measurement conducted on that sample. The measured austenite grain size was, respectively, 55 ± 3 µm and 71 ± 2 µm.
It is noteworthy that the fraction of pearlite in the sample reheated to the holding temperature contains more pearlite, which means that the austenite kept decomposing to this phase during heating of the sample (Figure 13).

Identification and Verification of the Austenite Microstructure Evolution Model
An inverse algorithm developed by the authors [11] was applied in the identification procedure. Identification of the microstructure evolutions model in Table 1 was divided into two parts. A static recrystallization model was identified on the basis of stress relaxation tests [21]. A dynamic recrystallization model was identified on the basis of compression tests at various temperatures and various strain rates. The coefficients in the microstructure evolution model for the three steels are given in Table 5. Identification of the phase transformation model was performed using inverse analysis for the dilatometric tests. In the analysis, the phase transformations start and finish temperatures, as well as volume fractions of structural constituents, determined in the dilatometric experiments, were used. Basic principles of application of this method to phase transformation models are described in [22]. In a majority of earlier applications of this method, results of the constant cooling rate (CCT) tests only were satisfactory to obtain reliable results. In the present project, however, heat treatment processes involve holding at the constant temperature in the range of pearlitic or bainitic transformation. Therefore, to improve the identification, the results of both CCT and TTT tests were combined. The following objective function was used: where: w T , w F , w t -weights, T m , T c -measured and calculated start and end temperatures of phase transformations in the CCT tests, F m , F c -measured and calculated volume fractions of phases, t m , t c -measured and calculated times to the start and end of phase transformations in the TTT tests, k, l, m-number of temperatures, volume fractions and times, respectively. The identification of the phase transformation model was performed using high performance computing (HPC) infrastructure. It was dictated by a large number of model coefficients and, as a consequence, requires huge computational costs to obtain reliable model. Over one hundred optimizations were performed using 50 nodes, each particular node had 24 computational cores. In this particular problem, global optimization algorithms were applied (PSO and Genetic), using from 24 to 96 specimens. Finally, the best model was selected. This approach allowed problems with local minima to be resolved and finally a reliable model to be obtained.
Coefficients in all models were determined by searching for the minimum of the objective Equation (19) with respect to these coefficients. The results are given in Table 6 for steel A, in Table 7 for steel B and in Table 8 for steel C.  The quality of the identification of the phase transformation model was evaluated via comparison of measured and calculated start and end temperatures for the phase transformations in the CCT tests and times to the start and end of transformations in the TTT tests. Selected results for steel B and steel C are shown in Figure 15 here. Reasonably good agreement was obtained. It should be noted that one model is used for different austenite grain size prior to transformations.   The model predicts very well volume fractions of structural components for ferritic/pearlitic microstructures. Discrepancies are larger when bainite and martensite are dominating. It is caused by much larger magnifications of the microstructure, which were necessary to distinguish these components. As a consequence, the area of analysis was constrained and was less representative for the whole material.

Numerical Simulation of Forging-Cooling Sequences
The validated model with optimal material parameters was used for simulations of hot forging and controlled cooling in one of the forges in Poland. Distributions of microstructural constituents in the forging after cooling were calculated. Selected results are described in the following sections.

Steel A
Calculated distributions of the temperature and austenite grain size after forging of the adapter are shown in Figure 18. Temperatures are at the level of 1100-1230 • C, and the austenite grain size varies in the range 30 ± 5 µm in the flash and 70 ± 4 µm in the massive part of the forging. These data were used as a starting point for simulation cooling according to the scheme shown in Figure 2a. Further simulations were performed for various holding temperatures and holding times, and distributions of phase volume fractions were calculated. Selected results for holding at the temperatures 560 • C and 520 • C are shown in Figure 19. It is seen that pearlite dominates for holding temperature 560 • C and bainite dominates for holding temperature 520 • C. At 560 • C, more pearlite was predicted for the thick end of the forging, which was due to larger austenite grains in that area. At 520 • C, distribution of the bainite volume fraction was more uniform. The martensite was a remaining phase in both cases. Figure 19. Cooling of the adapter from steel A-calculated distributions of the pearlite volume fraction after holding at 560 • C for 1 h (a) and of the bainite volume fraction after holding at 520 • C for 1 h (b).

Steel B
Calculated distributions of the temperature and austenite grain size after forging of the flange are shown in Figure 20. Temperatures are at the level of 950-1110 • C, and the austenite grain size varies in the range of 15 ± 2 µm in the thin part and 45 ± 3 µm in the massive part of the forging. These data were used as a starting point for simulation cooling according to the scheme shown in Figure 2b. Further simulations were performed for the cooling sequence shown in Figure 2b. Calculated distribution of the ferrite volume fraction after cooling to the room temperature is shown in Figure 21. The level of the ferrite fraction is in the range 72-75%. Pearlite is the remaining phase.  Figure 2c.
The first stage of cooling according to the schedule in Figure 2c was simulated next. Figure 23 shows the calculated distribution of the austenite volume fraction after cooling to the temperature at which ferritic transformation in the arms is completed and pearlitic transformation has not begun yet. It is seen that at this stage in the massive part of the forging the transformation is at the beginning stage and the volume fraction of the remaining austenite is in the range 80-95%.

Conclusions
The aim of the project was to develop the numerical tools capable of generating alternative cooling routes after forging as compared to conventional heat treatment methods. For this, an upgrade of the JMAK model was developed and applied to simulate various cooling strategies for hot forged parts. The upgrade of the model was composed of: (1) Introduction of the modified Gauss function to describe the relation of the Avrami coefficient k to the temperature. (2) Introduction of the austenite grain size prior to transformation as an independent variable in the model. Dilatometric tests were performed and supplied data for the identification of the phase transformation model for the three steels. Following this, the physical simulations of the selected forging processes were performed and the results were compared with the numerical simulations. The following conclusions were drawn: (1) The deformation process affects substantially the developed austenite microstructure and affects the kinetics of the decomposition of this phase to ferrite, pearlite and to a lesser extent to bainite and martensite. Therefore, it is important that this effect should be accounted for in the phase transformation model.
(2) The model predicts the start and end temperatures of transformations with good accuracy.
(3) The model predicts very well volume fractions of structural components for ferritic/pearlitic microstructures. (4) Despite the limitations, the methodology of the model validation and tuning adapted in this investigation using the results of the physical simulation of the forging process with Gleeble 3800 proved to be a very effective tool for improving the predictive capabilities of the mathematical models.
Larger discrepancies between measured and calculated volume fractions occurred when bainite and martensite were dominating components. This can be partly assigned to difficulties with reliable measurements of the bainite and martensite volume fraction. The other phenomenon not accounted for in the model relates to the chemical composition inhomogeneities related to the steel casting process.
Moreover, the global approach adapted in this investigation does not account for local phenomena, specifically those relating to non-equilibrium conditions prevailing at the physical interfaces created during phase transformations. However, this drawback of the developed methodology is counterbalanced by its simplicity and short time of the numerical calculations.