Next Article in Journal
A Simplified Ductile Fracture Model for Predicting Ultra-Low Cycle Fatigue of Structural Steels
Previous Article in Journal
A New Wood Adhesive Based on Recycling Camellia oleifera Cake-Protein: Preparation and Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Physical and Numerical Simulations for Predicting Distribution of Microstructural Features during Thermomechanical Processing of Steels

1
Łukasiewicz Research Network-Institute for Ferrous Metallurgy, ul. K. Miarki 12, 44-100 Gliwice, Poland
2
Department of Applied Computer Science and Modeling, Faculty of Metals Engineering and Industrial Computer Science, AGH University of Science and Technology, al. Mickiewicza 30, 30-059 Kraków, Poland
*
Author to whom correspondence should be addressed.
Materials 2022, 15(5), 1660; https://doi.org/10.3390/ma15051660
Submission received: 26 January 2022 / Revised: 14 February 2022 / Accepted: 17 February 2022 / Published: 23 February 2022
(This article belongs to the Topic Metallurgical and Materials Engineering)

Abstract

:
The design of modern construction materials with heterogeneous microstructures requires a numerical model that can predict the distribution of microstructural features instead of average values. The accuracy and reliability of such models depend on the proper identification of the coefficients for a particular material. This work was motivated by the need for advanced experimental data to identify stochastic material models. Extensive experiments were performed to supply data to identify a model of austenite microstructure evolution in steels during hot deformation and during the interpass times between deformations. Two sets of tests were performed. The first set involved hot compressions with a nominal strain of 1. The second set involved hot compressions with lower nominal strains, followed by holding at the deformation temperature for different times. Histograms of austenite grain size after each test were measured and used in the identification procedure. The stochastic model, which was developed elsewhere, was identified. Inverse analysis with the objective function based on the distance between the measured and calculated histograms was applied. Validation of the model was performed for the experiments, which were not used in the identification. The distance between the measured and calculated histograms was determined for each test using the Bhattacharyya metric and very low values were obtained. As a case study, the model with the optimal coefficients was applied to the simulation of the selected industrial hot-forming process.

1. Introduction

Enhancing the strength–ductility synergy of materials has been an objective of research on structural materials for years. It was shown in many publications that a significant improvement of this synergy can be obtained by tailoring heterogeneous microstructures. Multiphase steels (advanced high strength steels (AHSS)) are a leading example of the enhanced balance between strength and ductility [1]. The special mechanical properties of multiphase steels are due to the heterogeneity in strength between their structural components [2,3]. The heterogeneity of mechanical properties may come from differences in various microstructural features [4,5]. The authors of [6] investigated compositional, microstructural and local hardness variations in a complex-phase (CP) steel. They showed that the hierarchically heterogeneous microstructure with smooth gradients of properties was promising for maintaining the balance between strength and workability. This problem has been recently investigated by many researchers. Various gradient structures with varying ferrite grain size were investigated in [7], and a significant influence of grain size distribution on mechanical properties was shown. Simultaneous improvement of strength and plasticity was obtained in [8] using a novel strategy that emphasised the importance of work hardening and thickness of gradient layer (avoiding sharp gradients characteristic for thin layers). A similar hypothesis, which combines improved local fracture resistance with smooth gradients of properties, was put forward in [9]. Although the successful control or application of chemical or microstructural heterogeneity to achieve the desired properties was achieved by some researchers (e.g., for pipeline [10], multiphase [11] and ferritic/martensitic steels [12], as well as for other fields in engineering [13]), it seems that numerical modelling can still be useful support for the design of these materials. Advanced numerical models with the capability to predict distributions of various microstructural features are needed to reach this goal.
The stochastic model, which describes the microstructure evolution during multi-step hot forming of metallic materials, was developed in [14], and it is described in [15] in the deformation part (dynamic processes) and in [16] for interpass times (static processes). The correct evaluation of coefficients in the model for a particular material has a significant influence on the accuracy and reliability of the model. This goal is usually reached by performing experimental tests and the application of inverse analysis to find optimal coefficients in the model. The inverse solution for hot forming is well described in the literature [17,18,19], which also includes the extensive discussion of the existence and uniqueness of the solution in [20]. However, all published examples of inverse analysis concern deterministic models. The situation is different for the stochastic models when measurements of histograms of the output parameters have to be considered. The mathematical background of the inverse approach to the stochastic variables is described in [15]. Furthermore, the optimal numerical parameters of the model (number of Monte Carlo points, number of bins), which provide the best convergence and accuracy of optimisation, are proposed in that paper. Beyond this, the existence and convergence of the model are discussed and various measures of the distance between the predicted and experimental histograms are analysed. The practical application of inverse analysis to the stochastic model is described in [21]. To the best of the authors’ knowledge, research on the identification of stochastic microstructure evolution models is scarce. The published papers focused on the statistical inverse problem for the identification of a non-Gaussian tensor-valued random field, e.g., in [22], or the uncertainty of inverse analysis [23]. The analysis of the published works shows that there is a need for the development of the inverse approach to the identification of the stochastic microstructure models of steels. This approach requires advanced experiments, which supply information about the heterogeneity of microstructural parameters instead of the average values. Our objective was to perform experiments that would supply reliable and accurate data for the identification of the stochastic model described in [15]. In the investigation, a thermomechanical Gleeble 3800 simulator was used to perform experiments of controlled deformation to gain the data required for the model development. The experiments allowed for identifying the kinetics of structural changes occurring during and after deformation. Using this information, further experiments were planned with water quenching to “freeze” the austenite microstructure at different stages of processes occurring in the samples. Following this, the samples were quenched to reveal austenite grain boundaries, and the measurement of grain size distribution was performed on microstructure images.

2. Model

The model was developed in [14] and the subsequent steps of the development are described in the literature [15,16,21]. The general idea of the model, which is based on the stochastic internal variable (dislocation density) representing microstructure evolution, is presented in [15]. In Ref. [21], the grain size was included as the second stochastic variable, and the identification of the model based on the experimental data available in the literature for medium-carbon steel was performed. In Ref. [16], the model was extended by accounting for metadynamic and static recrystallisation during interpass times, which allowed for a simulation of multi-step forming processes. The main equations of the model are repeated below to complete the present paper. The evolution of dislocation density as a stochastic internal variable is governed by the following equation [15]:
ρ ( t i ) = ρ ( t 0 ) [ 1 ξ ( t i ) ] + { ρ ( t i 1 ) + [ A 1 ε ˙ A 2 ρ ( t i 1 ) ε ˙ 1 a 7 ] Δ t } ξ ( t i )
where t—time; ε ˙ —strain rate; A1 and A2—parameters of the model responsible for athermal storage (hardening) of dislocations and recovery, respectively; ξ(ti)—parameter responsible for a random character of recrystallisation; and a7—coefficient.
The coefficients in Equation (1), which are in general based on the Kocks–Estrin–Mecking (KEM) model [24,25], are defined in Table 1, where b—Burgers vector module; Z—Zener–Hollomon parameter; l—average free path for dislocations; T—temperature in K; R—universal gas constant; D—grain size; τ—energy per unit of dislocation length; G—shear modulus; and a1, a2, a3, a4, a5, a9 and a10—coefficients.
The parameter ξ(ti) accounts for the random character of recrystallisation, where its distribution is described by the following conditions [15]:
P [ ξ ( t i ) = 0 ] = { p ( t i ) i f p ( t i ) < 1 1 o t h e r w i s e P [ ξ ( t i ) = 1 ] = 1 P [ ξ ( t i ) = 0 ]
In Equation (2), p(ti) is a function that combines the probability that the material point recrystallises in a current time step and the present state of material:
p ( t i ) = a 4 × 10 10 ρ ( t i 1 ) a 6 3 γ ( t i ) τ D ( t i 1 ) exp ( a 5 R T ) Δ t
where γ—mobile fraction of the recrystallised grain boundary area, which depends on the (already known) distribution of ξ in the previous step (see [26]), and a4, a5, a6 and a17—coefficients.
γ ( t i ) = 1 exp { P [ ξ ( t i 1 ) = 0 ] q } a 8 { 1 P [ ξ ( t i 1 ) = 0 ] }
where a5, a6 and a8—coefficients, and q—a small number representing a nucleus of recrystallised grain, which is added to avoid a zero value of γ(ti) in the case of P[ξ(ti−1) = 0] = 0.
Details of the numerical solution are given in [15]. In [21], the grain size was introduced as a stochastic variable. The initial grain size is D(t0) ≡ D0, which is a random variable. During heating, a non-uniform distribution of grain size is observed (see Section 3), and there is a tendency for large grains to begin to consume small grains [27,28]; therefore, the Weibull distribution was assumed for D0. This distribution is described by the following probability density function:
f ( D 0 ) = k D ¯ 0 ( D 0 D ¯ 0 ) k 1 exp [ ( D 0 D ¯ 0 ) k ]
where k and D ¯ 0 —the shape parameter and the scale parameter of the distribution, respectively.
A shape parameter of k = 10 was assumed in [21], which was determined based on measured grain size distributions. The scale parameter D ¯ 0 of the distribution was established as the average grain size measured after preheating before deformation. Grain growth was calculated based on the fundamental research of Sellars [29], who proposed the following equation:
D ( t i ) = [ D ( t i 1 ) a 11 + a 12 exp ( a 13 R T ) Δ t ] 1 a 11
where Δt = titi−1—time step; D(ti−1) and D(t)—grain size at the beginning and the end of the time step, respectively; and a11, a12 and a13—coefficients.
During the calculation random parameter ξ(ti) = 0, the considered point recrystallises, and its new grain size D(ti) is drawn from the Gauss distribution:
f [ D ( t i ) ] = 1 2 π σ 2 exp { [ D ( t i ) D ¯ ( t i ) ] 2 2 σ 2 }
where D ¯ ( t i ) —the expected grain size value, calculated as either dynamically or statically recrystallised grain size; see details in [16].
The whole model contains 22 coefficients grouped in vector a. Some of these coefficients have physical meaning: a3—activation energy for self-diffusion, a5—activation energy for recrystallisation, a10—activation energy in the Zener–Hollomon parameter and a13—activation energy for grain growth. Other coefficients were introduced as the results of approximation in inverse analysis in previous research works. The application of the model to real materials and processing methods requires the identification of the model coefficients. The identification was performed in [16] using the inverse approach for the experimental data available in the literature. This identification, however, requires a special set of experimental data, including histograms of microstructure parameters at various stages of the process. These data were obtained in the experiments performed in the present work and described in Section 3.

3. Experiment

3.1. Material, Methodology and Parameters of the Tests

The material for the research was steel S355J2, which is an unalloyed, low-carbon, welded structural steel used for the hot forging and rolling of long products. As mentioned in the Introduction, two sets of tests were performed. All tests were carried out on the thermomechanical simulator Gleeble 3800. Conventional inverse analysis [19] was applied to eliminate the effects of friction and deformation heating. The first set consisted of hot uniaxial compression of φ10 × 12 mm cylindrical samples with a total strain of 1. Different strain rates and temperatures were applied. The objective was to supply data for the identification of the coefficients in the equations in Table 1 and Equations (2)–(4). The procedure of identification was based on the average dislocation density, as described in Section 3.2, and on the histograms of grain size after deformation. The parameters of the experiments in the first set of tests are given in Table 2. The second set of experiments consisted of hot uniaxial compression with lower nominal strains, followed by holding at the deformation temperature for different times th after the completion of the recrystallisation. The time for the completion of the static recrystallisation was determined using the stress relaxation technique [30]. The objective was to supply data for the identification of static recovery and static recrystallisation parts in the model, see [19] for details The parameters of the experiments in the second set of tests are given in Table 3. As with previous experiments, two preheating temperatures were used: 1200 °C and 1100 °C. The preheating time was 120 s in all tests.
The images used as the input for the analysis and measurement of austenite grain sizes were taken using an Olympus DSX500i light microscope at 693× magnification via the bright field technique. The equivalent diameter was used as a measure of grain size, and the measurements were performed in the central part of the sample. To reveal the grain boundaries of austenite, the samples were etched in saturated picric acid with the addition of detergent in the amount of 5 mL per 100 mL of acid at 56 °C for the time necessary to reveal the boundaries of prior austenite. This time was dependent on the deformation history. For the clear etching of grain boundaries, the material for testing was additionally subject to tempering in a Carbolite Chamber Furnace C-1600. This treatment created the desired effect, but some of the grains remained undisclosed (Figure 1a), which made further analysis much more difficult. The determination of grain size distributions and their measurement was carried out using the MetIlo program. On average, about 200 grain-size measurements were taken for each sample. In the course of the measurement, manual austenite grain correction was applied to avoid errors connected with the wrong automatic identification of boundaries. An example image of austenite grain detection after the manual correction is shown in Figure 1b.

3.2. Dislocation Density

Measurement of the distribution of dislocation density during hot forming is not possible for steels subject to allotropic transformations. Therefore, as it was suggested in [21], the average dislocation density calculated from the measured compression forces was used in the identification procedure. Uniaxial compression tests supplied data in the form of force vs. die displacement curves, which are shown in Figure 2. An analysis of the results showed that differences between loads for various preheating temperatures were negligible. An increase in the slope of the plots during deformation was caused by the effect of friction, which increased with an increase in the tool–workpiece contact surface. This effect was accounted for by the application of conventional inverse analysis [19].
Calculations of the flow stress vs. strain curves required accounting for the effect of heterogeneities due to the effect of friction, heat generation due to plastic work and heat transfer due to the die. The inverse approach proposed in [19] was applied. The flow stress vs. strain curves were calculated by searching for the minimum of the following objective function:
Φ = 1 N t i = 1 N t [ 1 N s j = 1 N s ( F i j m F i j c F i j m ) 2 ]
where Nt—number of tests; Ns—number of load measurement sampling points in one test; and F i j m and F i j c —measured and calculated forces, respectively.
The flow stress determined using inverse analysis was used to calculate the average dislocation density as a function of time. The changes in the average dislocation density were calculated using the following equation based on the KEM [24,25] model:
ρ m ( t ) = [ σ m ( t ) α M b G ] 2
where σm—flow stress calculated from the measured compression forces, b—length of the Burgers vector, G—shear modulus, M—the Taylor factor (≈3 for an FCC structure) and α—a constant called the Taylor coefficient hereafter. Following the discussion in [31], α = 0.5 was assumed in the present work.
The accuracy of the calculations of the average dislocation density (ρm) depends on the correctness of the evaluation of the shear modulus G, which is a function of temperature. The dependence of elastic modulus E and shear modulus G on the temperature was investigated by several studies but, unfortunately, there are large discrepancies between the published data (cf. [32,33,34]). The authors of [34] recapitulated the results for various steels and proposed the following approximation function for the relation E = E(T).
E = E 0 { b 1 exp [ ( T b 2 b 3 ) 2 ] + c 1 exp [ ( T c 2 c 3 ) 2 ] }
where E0—elastic modulus at room temperature.
The coefficients b1, b2, b3, c1, c2 and c3 that appear in Equation (10) were determined in [34] and they are given in Table 4. In the present work, to make the results closer to the conditions of hot forming of steels, coefficient c2 was slightly modified [21].
Plots of the average dislocation density vs. logarithmic strain, which were calculated from the measurements of forces during compression, are shown in Figure 3. The logarithmic strain is defined as ε = ln(h1/h), where h1 and h—height of the sample before the compression and during the compression, respectively. The effect of the dynamic recrystallisation of the evolution of the dislocation density is clearly visible in this figure. The dislocation density data were used in the identification of the model.

3.3. Histograms of Grain Size after Preheating (Prior to Deformation)

After each test, the samples were etched and austenite grain boundaries were revealed. Micrographs of the samples after preheating are shown in Figure 4. Histograms of the grain size prior to deformation for the two preheating temperatures (Tp) are shown in Figure 5. The average grain size for the histograms was 36.5 μm and 22.5 μm for Tp = 1200 °C and Tp = 1100 °C, respectively.

3.4. Grain Size after Hot Deformation (Dynamic Recrystallisation)

All the samples were subjected to metallographic investigations; selected microstructures of samples preheated at 1200 °C followed by deformation are shown in Figure 6, Figure 7, Figure 8 and Figure 9.
Generally, a substantial spread of the measured grain diameters was observed at a higher deformation temperature (1200 °C vs. 1100 °C). Moreover, the spread was larger for lower strain rates. For the applied deformation temperatures, dynamic recrystallisation was initiated in the steels, which is coherent with the fact that grain size in the steady-state regime of the deformation depends on the Zner–Hollomon parameter.

3.5. Grain Size during Interpass Times (Metadynamic and Static Recrystallisation)

Selected examples of grain size histograms measured at various times after the end of deformation are shown in Figure 10. The parameters of the tests are given in Table 3. It can be seen that a substantial spread of the distribution of equivalent diameter was obtained. This can be attributed to the spread of the initial size distribution, as well as to the partially stochastic nature of the static recrystallisation phenomenon. As opposed to the samples subject to dynamic recrystallisation, the histograms of equivalent diameter were strain-dependent.

4. Identification and Validation of the Model

The identification procedure is described in [16]. Inverse analysis was applied to determined the coefficients a in the model. The aim of this analysis was to find the optimal values of these coefficients by searching for the minimum of the following objective function:
Φ ( a ) = Φ ρ ( a ) + Φ D D R X ( a ) + Φ D S R X ( a )
where
Φ ρ ( a ) = w ρ i = 1 N t c d ( ρ c i ( a ) , ρ m i )
Φ D D R X ( a ) = w D i = 1 N t c d ( H c i ( a ) , H m i )
Φ D S R X ( a ) = w D i = 1 N t i d ( H c i ( a ) , H m i )
Here, ρc(a)—expected average value of dislocation density calculated for the model coefficients a, ρm—average dislocation density determined from the compression tests (see Section 3.2), Hc(a)—distribution of grain size calculated for the model using the coefficients a, Hm—distribution of grain size measured in the experiments, Ntc—number of compression tests, Nti—number of measurements of histograms during interpass times and wρ and wD—weighted coefficients.
In Equation (12), d(ρci(a),ρmi) was defined as the mean square root error (MSRE) between the measured and the calculated average dislocation density in the ith experiment:
d ( ρ c i ( a ) , ρ m i ) = 1 N s j = 1 N s ( ρ c i j ( a ) ρ m i j ρ m i j ) 2
where Ns—number of sampling points for measurements of average dislocation density in the ith test.
In parts of the objective Equations (13) and (14), d(Hci(a),Hmi) was defined as the distance between the measured and the calculated histograms of grain size in the ith experiment. Following the analysis in [15], the Bhattacharyya [35] metric was used to calculate this distance. This metric is calculated as follows:
d ( H c , H m ) = log i = 1 n P m ( i ) P c ( i )
where n—number of histogram bins; Hm and Hc—measured and calculated histograms, respectively; and
P k ( i ) = H k ( i ) i = 1 n H k ( i )
The experiments were divided into two sets. Hot deformation tests (T1, T3, T5, T7, T13, T15, T19 and T20) and static recrystallisation tests (T25, T28, T32 and T35) were used in the identification procedure, and the remaining tests were used for the validation of the model. The coefficients in the model determined for the objective Equation (11) are given in Table 5.

5. Model Validation

Predictions of the model with the optimal coefficients were compared with the experiments, and selected results are presented in Figure 11. Some selected results for static recrystallisation tests are shown in Figure 12. The Bhattacharaya distance is given in the bottom-left corner of the plot. The visual comparison of the measured and calculated histograms confirmed the qualitatively good predictive capability of the model.
The accuracy of the model can be better evaluated using a quantitative comparison of the distance between the measured and calculated histograms defined by Equation (16). The values of the Bhattacharyya distance for all hot deformation tests are shown in Figure 13 and for all static recrystallisation tests in Figure 14.
The comparison between the measured and calculated average grain sizes was an additional evaluation of the accuracy of the model. These results are shown in Figure 15 for all hot deformation tests and in Figure 16 for all static recrystallisation tests. The comparisons of the measured and calculated grain sizes generated conclusions that qualitatively agreed with the observations based on the Bhattacharyya distance shown in Figure 13 and Figure 14.

6. Discussion of Results

Advanced stochastic models, which predict the distribution of microstructural features instead of average values, can describe the heterogeneity of microstructures of modern multiphase steels. The accuracy and reliability of such models depend on the proper identification of coefficients for a particular material. The need for advanced experimental data for the identification of the stochastic material model was a motivation for this work.
Extensive experiments were performed to supply data for the identification of the model of austenite microstructure evolution in steels during multi-step hot deformation and during interpass times between deformations. The applied experimental methodology allowed for the precise identification and quantitative description of the processes involved in austenite microstructure restoration in the process of deformation using controlled deformation parameters. The crucial condition to generate an accurate model is to precisely perform austenite boundary etching. Boundaries that are not clearly etched need to be manually reconstructed, which may lead to a significant error in the process of generation of austenite grain size distribution histograms.
The analysis of the optimisation procedure in the inverse solution showed that the Bhattacharyya metric was an efficient and reliable measure of the distance between the measured and calculated histograms. A good convergence of the optimisation was observed, and reasonably low values of the Bhattacharyya metric were obtained for all the tests. The results confirmed the good accuracy and reliability of the inverse stochastic approach. It can be assumed that a Bhattacharyya metric below 0.3 is acceptable for the comparison of microstructures. The predictions of the model with optimal coefficients were in good agreement with the measurements. The Bhattacharyya metric was below 0.2 for the experiments used in the identification, and below 0.33 for the experiments used in the validation. The value of 0.3 was slightly exceeded in a few tests, which may have been due to experimental errors. In the tests with the largest Bhattacharyya metrics, the standard deviations in the experiments were large (e.g., T1, T26, T31). This problem will be the subject of further research.
In the hot deformation experiments, better accuracy was obtained for tests at higher temperatures and lower strain rates (e.g., T1–T4) in which dynamic recrystallisation dominated. In the static recrystallisation experiments, better accuracy was obtained for the tests with larger strains and at higher temperatures (e.g., T28–T30) in which recrystallisation was faster.

7. Case Study

Simulations of hot forging industrial processes were performed to demonstrate the predictive capabilities of the stochastic model. Forge 3D finite element (FE) software was used to simulate the thermal and mechanical phenomena during deformation and during interpass times. Forge is a computational package dedicated to the simulation of metal forming processes. The software consists of a loosely coupled set of programs, bound by a common launcher user interface. The package also consists of a database with material models. Forge 3D offers automatic adaptive remeshing, which is the main advantage of this software. The post-processing application takes a project input file and uses its data for estimating results and showing them to the user through the graphical user interface (GUI). Three-step forging of the part shown in Figure 17a was considered. The details of the industrial forging process for this part are given in [36], where more information on the application of Forge 3D to simulations of hot forging can be found. The results of FE simulations in the form of temperature distribution after subsequent stages of forging are presented in Figure 17b–d.
The solution of the stochastic evolution Equation (1) can be performed for each Gauss point of the FE mesh, accounting for the current local temperatures and strain rates. However, such detailed information is usually not needed, and it is satisfactory if the histograms are calculated in the selected critical points in the forging. IT was done by using the “sensors” in the Forge software. Forge allows for selecting a few points (sensors) in the volume of the forging in which the time–temperature–deformation history is recorded. This history was used to solve Equation (1) and to calculate histograms of dislocation density and grain size in the selected points. The results for points 1 and 2 in Figure 17a are presented in Figure 18. These points represented the massive part and thin part of the forging and were subjected to different time–temperature–deformation histories.
A simulation of the manufacturing cycle composed of three-stage forging followed by accelerated cooling with a rate of 25 °C/s to the temperature 820 °C was performed. The calculated distributions of dislocation density and grain size at the end of forging and after cooling to 820 °C are shown in Figure 19. The analysis of the results confirmed the model’s capability to supply quantitative information about the distribution of microstructural parameters in the complex industrial hot-forming processes. The predictions were in qualitative agreement with the authors’ knowledge about the forging process. For point 2, which was subjected to larger deformations, the dynamic recrystallisation was faster and the grain size was smaller. This tendency was maintained during cooling. Finer grains and a larger static recrystallisation volume fraction were obtained for point 2.
The representation of the dislocation density histograms was not convincing. This model was based on the assumption that the bins in the histograms were of equal length. Consequently, the majority of the material either had a very large dislocation density in the last bin or, after recrystallisation, had very low dislocation density and was located in the first bin. The volume fraction of the material with intermediate dislocation densities was negligible. In the future, the authors will change the model and introduce a different length of the bins.

8. Conclusions

The identification of the stochastic model of hot forming for metallic materials requires experimental data, which supplies information about distributions of microstructural parameters for different process conditions. Tests composed of hot deformation and holding after deformation for different times were performed. The following conclusions were drawn:
  • The objective function in the inverse analysis was formulated as a distance between the measured and the calculated histograms. It was shown that using the Bhattacharyya metric as a measure of this distance was a very efficient approach. The results confirmed the good accuracy and reliability of the inverse stochastic approach. The predictions of the model with optimal coefficients agreed with the measurements. The Bhattacharyya metric was below 0.2 for the experiments used in the identification and below 0.33 for the experiments used in the validation, which is a reasonable accuracy when steel microstructures are compared.
  • In the hot deformation experiments, better accuracy was obtained for the tests at higher temperatures and lower strain rates (e.g., T1–T4) in which dynamic recrystallisation dominated. In the static recrystallisation experiments, better accuracy was obtained for the tests with larger strain (T28–T30) in which recrystallisation was faster.
  • It was observed that the accuracy of the model identification was directly dependent on the quality of the etching of austenite grain boundaries. Since, typically, not all boundaries were clearly etched, a manual reconstruction of the missing boundaries segments needed to be done, which could have resulted in significant errors.
  • The lower accuracy of the developed model for low strains could be connected with a significant overlapping of the recovery and recrystallisation processes, which needs the improved methodology of the separation of their effect in plastometric studies. This will be the subject for future activities.
  • The model with optimal coefficients was applied to the simulation of the three-stage hot-forging process. The predictions of the model were in qualitative agreement with the authors’ knowledge regarding hot forging, and the extensive predictive capabilities of the model were confirmed. Dislocation density (when the recrystallisation of austenite is not completed) and grain size at the beginning of phase transformation have a strong influence on the kinetics of these transformations. Thus, calculated distributions of the microstructural features can be used as a starting point for stochastic modelling of the phase transformation and prediction of histograms of phase composition in the final product. This problem will be the scope of the authors’ future research.

Author Contributions

Conceptualisation, all Authors; methodology, Ł.P., R.K., D.S. and M.P.; software, D.S., J.K. and M.P.; validation, All Authors; formal analysis, D.S., M.P. and J.K.; investigation, Ł.P., R.K. and M.P.; resources, Ł.P., R.K. and V.P.; data curation, Ł.P., R.K. and V.P.; writing—original draft preparation, Ł.P., R.K. and M.P.; writing—review and editing, All Authors; visualisation, D.S. and M.P.; supervision, Ł.P.; project administration, Ł.P., R.K. and M.P.; funding acquisition, Ł.P. and R.K. All authors have read and agreed to the published version of the manuscript.

Funding

The financial support of NCBiR, project No. Techmatstrateg 1/348491/10/NCBR/2017, is acknowledged.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kuziak, R.; Kawalla, R.; Waengler, S. Advanced high strength steels for automotive industry. Arch. Civ. Mech. Eng. 2008, 8, 103–117. [Google Scholar] [CrossRef]
  2. Wu, X.; Zhu, Y. Heterogeneous materials: A new class of materials with unprecedented mechanical properties. Mater. Res. Lett. 2017, 5, 527–532. [Google Scholar] [CrossRef]
  3. Chang, Y.; Haase, C.; Szeliga, D.; Madej, L.; Hangen, U.; Pietrzyk, M.; Bleck, W. Compositional heterogeneity in multiphase steels: Characterization and influence on local properties. Mater. Sci. Eng. A 2021, 827, 142078. [Google Scholar] [CrossRef]
  4. Motaman, S.A.H.; Roters, F.; Haase, C. Anisotropic polycrystal plasticity due to microstructural heterogeneity: A multi-scale experimental and numerical study on additively manufactured metallic materials. Acta Mater. 2020, 185, 340–369. [Google Scholar] [CrossRef]
  5. Motaman, S.A.H.; Haase, C. The microstructural effects on the mechanical response of polycrystals: A comparative experimental-numerical study on conventionally and additively manufactured metallic materials. Int. J. Plast. 2021, 140, 102941. [Google Scholar] [CrossRef]
  6. Chang, Y.; Lin, M.; Hangen, U.; Richter, S.; Haase, C.; Bleck, W. Revealing the relation between microstructural heterogeneities and local mechanical properties of complex-phase steel by correlative electron microscopy and nanoindentation characterization. Mater. Des. 2021, 203, 109620. [Google Scholar] [CrossRef]
  7. Wang, L.; Li, B.; Shi, Y.; Huang, G.; Song, W.; Li, S. Optimizing mechanical properties of gradient-structured low-carbon steel by manipulating grain size distribution. Mater. Sci. Eng. A 2019, 743, 309–313. [Google Scholar] [CrossRef]
  8. Shao, C.W.; Zhang, P.; Zhu, Y.K.; Zhang, Z.J.; Tian, Y.Z.; Zhang, Z.F. Simultaneous improvement of strength and plasticity: Additional work-hardening from gradient microstructure. Acta Mater. 2018, 145, 413–428. [Google Scholar] [CrossRef]
  9. Szeliga, D.; Chang, Y.; Bleck, W.; Pietrzyk, M. Evaluation of using distribution functions for mean field modelling of multiphase steels. Proc. Manuf. 2019, 27, 72–77. [Google Scholar] [CrossRef]
  10. Hassan, S.F.; Al-Wadei, H. Heterogeneous microstructure of low-carbon microalloyed steel and mechanical properties. J. Mater. Eng. Perform. 2020, 29, 7045–7051. [Google Scholar] [CrossRef]
  11. Heibel, S.; Dettinger, T.; Nester, W.; Clausmeyer, T.; Tekkaya, A.E. Damage mechanisms and mechanical properties of high-strength multi-phase steels. Materials 2018, 11, 761. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Li, S.; Vajragupta, N.; Biswas, A.; Tang, W.; Wang, H.; Kostka, A.; Yang, X.; Hartmaier, A. Effect of microstructure heterogeneity on the mechanical properties of friction stir welded reduced activation ferritic/martensitic steel. Scr. Mater. 2022, 207, 114306. [Google Scholar] [CrossRef]
  13. Ding, R.; Yao, Y.; Sun, B.; Liu, G.; He, J.; Li, T.; Wan, X.; Dai, Z.; Ponge, D.; Raabe, D.; et al. Chemical boundary engineering: A new route toward lean, ultrastrong yet ductile steels. Sci. Adv. 2020, 6, eaay1430. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Pietrzyk, M.; Czyżewska, N.; Klimczak, K.; Kusiak, J.; Morkisz, P.; Oprocha, P.; Przybyłowicz, P.; Rauch, Ł.; Szeliga, D. Stochastic approach to modelling microstructure evolution and properties of metals subjected to thermomechanical processing. NCN project no. 2017/25/B/ST8/01823, 2018-2021.
  15. Klimczak, K.; Oprocha, P.; Kusiak, J.; Szeliga, D.; Morkisz, P.; Przybyłowicz, P.; Czyżewska, N.; Pietrzyk, M. Inverse problem in stochastic approach to modelling of microstructural parameters in metallic materials during processing. Math. Probl. Eng. 2022. in review. [Google Scholar]
  16. Szeliga, D.; Czyżewska, N.; Klimczak, K.; Kusiak, J.; Kuziak, R.; Morkisz, P.; Oprocha, P.; Poloczek, Ł.; Pietrzyk, M.; Przybyłowicz, P. Stochastic through process model describing evolution of microstructural parameters during multi-step hot forming processes. Arch. Civ. Mech. Eng. 2022. in review. [Google Scholar]
  17. Gavrus, A.; Massoni, E.; Chenot, J.L. An inverse analysis using a finite element model for identification of rheological parameters. J. Mater. Proc. Technol. 1996, 60, 447–454. [Google Scholar] [CrossRef]
  18. Forestier, R.; Massoni, E.; Chastel, Y. Estimation of constitutive parameters using an inverse method coupled to a 3D finite element software. J. Mater. Proc. Technol. 2002, 125, 594–601. [Google Scholar] [CrossRef]
  19. Szeliga, D.; Gawąd, J.; Pietrzyk, M. Inverse analysis for identification of rheological and friction models in metal forming. Comput. Methods Appl. Mech. Eng. 2006, 195, 6778–6798. [Google Scholar] [CrossRef]
  20. Pietrzyk, M.; Madej, Ł.; Rauch, Ł.; Szeliga, D. Computational Materials Engineering: Achieving High Accuracy and Efficiency in Metals Processing Simulations; Butterworth-Heinemann: Oxford, UK; Elsevier: Amsterdam, The Netherlands, 2015. [Google Scholar]
  21. Szeliga, D.; Czyżewska, N.; Klimczak, K.; Kusiak, J.; Kuziak, R.; Morkisz, P.; Oprocha, P.; Pidvysotsk’yy, V.; Pietrzyk, M.; Przybyłowicz, P. Identification and validation of the stochastic model describing evolution of microstructural parameters during hot forming of metallic materials. Int. J. Mater. Form. 2022. in review. [Google Scholar]
  22. Soize, C.; Desceliers, C.; Guilleminot, J.; Le, T.-T.; Nguyen, M.-T.; Perrin, G.; Allain, H.; Gharbi, J.-M.; Duhamel, D.; Funfschilling, C. Stochastic representations and statistical inverse identification for uncertainty quantification in computational mechanics. In Proceedings of the 1st ECCOMAS Thematic International Conference on Uncertainty Quantification in Computational Sciences and Engineering, Athens, Greece, 21–23 June 2021; pp. 1–26. [Google Scholar]
  23. Ekwaro-Osire, S.; Dhorje, M.; Khandaker, M. Accounting for uncertainty in the inverse problem. In Proceedings of the XIth International Congress and Exposition, Orlando, FL, USA, 2–5 June 2008. [Google Scholar]
  24. Mecking, H.; Kocks, U.F. Kinetics of flow and strain-hardening. Acta Metall. 1981, 29, 1865–1875. [Google Scholar] [CrossRef]
  25. Estrin, Y.; Mecking, H. A unified phenomenological description of work hardening and creep based on one-parameter models. Acta Metall. 1984, 32, 57–70. [Google Scholar] [CrossRef]
  26. Sandstrom, R.; Lagneborg, R. A model for hot working occurring by recrystallization. Acta Metall. 1975, 23, 387–398. [Google Scholar] [CrossRef]
  27. Kaonda, M.K.M. Prediction of the Recrystallised Grain Size Distribution after Deformation for the Nb Free and Model Steel. Ph.D. Thesis, University of Birmingham, Birmingham, UK, 2017. [Google Scholar]
  28. Wang, J.; Kang, Y.; Yang, C.; Wang, Y. Effect of heating temperature on the grain size and titanium solid-solution of titanium microalloyed steels. Mater. Sci. Appl. 2019, 10, 558–567. [Google Scholar] [CrossRef] [Green Version]
  29. Sellars, C.M. Physical metallurgy of hot working. In Hot Working and Forming Processes; Sellars, C.M., Davies, G.J., Eds.; The Metals Society: London, UK, 1979; pp. 3–15. [Google Scholar]
  30. Somani, M.C.; Porter, D.A.; Karjalainen, L.P.; Kantanen, P.K.; Kömi, J.I.; Misra, D.K. Static recrystallization characteristics and kinetics of high-silicon steels for direct quenching and partitioning. Int. J. Mater. Res. 2019, 110, 183–193. [Google Scholar] [CrossRef] [Green Version]
  31. Kang, J.-H.; Kim, S.-J. Critical Assessment 33: Dislocation density-based constitutive modelling for steels with austenite. Mater. Sci. Technol. 2019, 35, 1128–1132. [Google Scholar] [CrossRef]
  32. Fukuhar, M.; Sanpei, A. Elastic steels moduli and internal friction as a function of temperature. ISIJ Int. 1993, 33, 508–512. [Google Scholar] [CrossRef]
  33. Outinen, J.; Mäkeläinen, P. Mechanical properties of structural steel at elevated temperatures and after cooling down. In Proceedings of the 2nd International Workshop Structures in Fire, Christchurch, New Zealand, 18–19 March 2002; pp. 273–298. [Google Scholar]
  34. Maraveas, C.; Fasoulakis, Z.C.; Tsavdaridis, K.D. Mechanical properties of high and very high steel at elevated temperatures and after cooling down. Fire Sci. Rev. 2017, 6, 1–13. [Google Scholar]
  35. Bhattacharyya, A. On a measure of divergence between two multinomial populations. Indian J. Stat. 1946, 7, 401–406. [Google Scholar] [CrossRef] [Green Version]
  36. Poloczek, Ł.; Rauch, Ł.; Wilkus, M.; Bachniak, D.; Zalecki, W.; Rozmus, R.; Kuziak, R.; Pietrzyk, M. Physical and numerical simulations of closed die hot forging and heat treatment of forged parts. Materials 2021, 14, 15. [Google Scholar] [CrossRef]
Figure 1. Example of the detection of austenite grain boundaries’ revealed microstructure (a) and the detection obtained using the MetIlo software (b) for sample T35.
Figure 1. Example of the detection of austenite grain boundaries’ revealed microstructure (a) and the detection obtained using the MetIlo software (b) for sample T35.
Materials 15 01660 g001
Figure 2. Results of force vs. die displacement in uniaxial compression tests for nominal strain rates 0.1 s−1 (a), 1 s−1 (b) and 10 s−1 (c).
Figure 2. Results of force vs. die displacement in uniaxial compression tests for nominal strain rates 0.1 s−1 (a), 1 s−1 (b) and 10 s−1 (c).
Materials 15 01660 g002aMaterials 15 01660 g002b
Figure 3. Average dislocation density vs. logarithmic strain calculated from the measurements of forces during compression for nominal strain rates 0.1 s−1 (a), 1 s−1 (b) and 10 s−1 (c).
Figure 3. Average dislocation density vs. logarithmic strain calculated from the measurements of forces during compression for nominal strain rates 0.1 s−1 (a), 1 s−1 (b) and 10 s−1 (c).
Materials 15 01660 g003
Figure 4. Micrographs of samples after preheating for 120 s at 1200 °C (a) and 1100 °C (b).
Figure 4. Micrographs of samples after preheating for 120 s at 1200 °C (a) and 1100 °C (b).
Materials 15 01660 g004
Figure 5. Histograms of grain size prior to deformation for the preheating temperatures 1200 °C and 1100 °C.
Figure 5. Histograms of grain size prior to deformation for the preheating temperatures 1200 °C and 1100 °C.
Materials 15 01660 g005
Figure 6. Micrographs after deformation at 1200 °C with strain rates 0.1 s−1 (a), 1 s−1 (b) and 10 s−1 (c), preheating temperature 1200 °C.
Figure 6. Micrographs after deformation at 1200 °C with strain rates 0.1 s−1 (a), 1 s−1 (b) and 10 s−1 (c), preheating temperature 1200 °C.
Materials 15 01660 g006
Figure 7. Micrographs after deformation at 1100 °C with strain rates 0.1 s−1 (a), 1 s−1 (b) and 10 s−1 (c), preheating temperature 1200 °C.
Figure 7. Micrographs after deformation at 1100 °C with strain rates 0.1 s−1 (a), 1 s−1 (b) and 10 s−1 (c), preheating temperature 1200 °C.
Materials 15 01660 g007
Figure 8. Micrographs after deformation at 1000 °C with strain rates 0.1 s−1 (a), 1 s−1 (b) and 10 s−1 (c), preheating temperature 1200 °C.
Figure 8. Micrographs after deformation at 1000 °C with strain rates 0.1 s−1 (a), 1 s−1 (b) and 10 s−1 (c), preheating temperature 1200 °C.
Materials 15 01660 g008
Figure 9. Micrographs after deformation at 900 °C with strain rates 0.1 s−1 (a), 1 s−1 (b) and 10 s−1 (c), preheating temperature 1200 °C.
Figure 9. Micrographs after deformation at 900 °C with strain rates 0.1 s−1 (a), 1 s−1 (b) and 10 s−1 (c), preheating temperature 1200 °C.
Materials 15 01660 g009
Figure 10. Examples of histograms measured at the end of selected static recrystallisation tests.
Figure 10. Examples of histograms measured at the end of selected static recrystallisation tests.
Materials 15 01660 g010
Figure 11. Selected examples of the comparison between measured and calculated histograms of the grain size after hot deformation tests: T1 (a), T2 (b), T13 (c) and T14 (d).
Figure 11. Selected examples of the comparison between measured and calculated histograms of the grain size after hot deformation tests: T1 (a), T2 (b), T13 (c) and T14 (d).
Materials 15 01660 g011
Figure 12. Selected examples of the comparison between measured and calculated histograms of the grain size measured at the end of the static recrystallisation tests T25 (a) and T29 (b).
Figure 12. Selected examples of the comparison between measured and calculated histograms of the grain size measured at the end of the static recrystallisation tests T25 (a) and T29 (b).
Materials 15 01660 g012
Figure 13. Bhattacharyya distance calculated for all hot deformation tests for preheating temperatures 1200 °C (a) and 1100 °C (b).
Figure 13. Bhattacharyya distance calculated for all hot deformation tests for preheating temperatures 1200 °C (a) and 1100 °C (b).
Materials 15 01660 g013
Figure 14. Bhattacharyya distance calculated for all static recrystallisation tests.
Figure 14. Bhattacharyya distance calculated for all static recrystallisation tests.
Materials 15 01660 g014
Figure 15. Comparison of the measured and calculated average grain sizes for all hot deformation tests for preheating temperatures 1200 °C (a) and 1100 °C (b).
Figure 15. Comparison of the measured and calculated average grain sizes for all hot deformation tests for preheating temperatures 1200 °C (a) and 1100 °C (b).
Materials 15 01660 g015
Figure 16. Comparison of the measured and calculated average grain sizes for all static recrystallisation tests.
Figure 16. Comparison of the measured and calculated average grain sizes for all static recrystallisation tests.
Materials 15 01660 g016
Figure 17. View of the hot-forged part investigated in the present work (a), distributions of temperature after stage 1 (b), stage 2 (c) and stage 3 (d) of the forging process.
Figure 17. View of the hot-forged part investigated in the present work (a), distributions of temperature after stage 1 (b), stage 2 (c) and stage 3 (d) of the forging process.
Materials 15 01660 g017
Figure 18. Time–temperature–deformation histories for points 1 and 2 in Figure 17a.
Figure 18. Time–temperature–deformation histories for points 1 and 2 in Figure 17a.
Materials 15 01660 g018
Figure 19. Calculated distributions of dislocation density (a,c) and grain size (b,d) at the end of the three-stage forging (a,b) and after cooling to 820 °C with a rate of 25 °C/s (c,d).
Figure 19. Calculated distributions of dislocation density (a,c) and grain size (b,d) at the end of the three-stage forging (a,b) and after cooling to 820 °C with a rate of 25 °C/s (c,d).
Materials 15 01660 g019
Table 1. Relationships describing the coefficients in Equation (1).
Table 1. Relationships describing the coefficients in Equation (1).
Hardening A 1 = 1 b l where l = a 1 Z a 9 ,  Z = ε ˙ exp ( a 10 R T )
Recovery A 2 = a 2 exp ( a 3 R T )
Table 2. Parameters of hot deformation tests for the preheating temperatures of 1200 °C and 1100 °C. The samples were water quenched directly after deformation.
Table 2. Parameters of hot deformation tests for the preheating temperatures of 1200 °C and 1100 °C. The samples were water quenched directly after deformation.
Tp = 1200 °C, ε = 1Tp = 1100 °C, ε = 1
TestTd, °C ε ˙ , S−1TestTd, °C ε ˙ , S−1
T112000.1T1611000.1
T212001T1711001
T3120010T18110010
T411000.1T1910000.1
T511001T2010001
T6110010T21100010
T710000.1T229000.1
T810001T239001
T9100010T2490010
T109000.1
T119001
T1290010
Table 3. Parameters of static recrystallisation tests for the preheating temperatures of 1200 °C and 1100 °C.
Table 3. Parameters of static recrystallisation tests for the preheating temperatures of 1200 °C and 1100 °C.
Tp = 1200 °CTp = 1100 °C
TestTd, °C ε ˙ , S−1εth, sTestTd, °C ε ˙ , S−1εth, s
T25100010.112T32100010.28
T261000100.27T331000100.27
T2790010.242T3490010.238
T28100010.41.7T35900100.217
T29900100.48.5
T3090010.412
T31900100.244
Table 4. Coefficients in Equation (10) for the investigated steel.
Table 4. Coefficients in Equation (10) for the investigated steel.
b1b2, °Cb3, °Cc1c2, °Cc3, °C
0.1871−12.961320.9199181.6670
Table 5. Optimal coefficients in the stochastic model.
Table 5. Optimal coefficients in the stochastic model.
a1a2a3a4a5a6a7
0.000846302256334,9361.145290,4761.79250.295
a8a9a10a11a12a13a14
0.7860.283,740,9454.35 × 1021410,00067,484
a15a16a17a18a19a20a21
0.27420.910.5450.290.196430
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Poloczek, Ł.; Kuziak, R.; Pidvysots’kyy, V.; Szeliga, D.; Kusiak, J.; Pietrzyk, M. Physical and Numerical Simulations for Predicting Distribution of Microstructural Features during Thermomechanical Processing of Steels. Materials 2022, 15, 1660. https://doi.org/10.3390/ma15051660

AMA Style

Poloczek Ł, Kuziak R, Pidvysots’kyy V, Szeliga D, Kusiak J, Pietrzyk M. Physical and Numerical Simulations for Predicting Distribution of Microstructural Features during Thermomechanical Processing of Steels. Materials. 2022; 15(5):1660. https://doi.org/10.3390/ma15051660

Chicago/Turabian Style

Poloczek, Łukasz, Roman Kuziak, Valeriy Pidvysots’kyy, Danuta Szeliga, Jan Kusiak, and Maciej Pietrzyk. 2022. "Physical and Numerical Simulations for Predicting Distribution of Microstructural Features during Thermomechanical Processing of Steels" Materials 15, no. 5: 1660. https://doi.org/10.3390/ma15051660

APA Style

Poloczek, Ł., Kuziak, R., Pidvysots’kyy, V., Szeliga, D., Kusiak, J., & Pietrzyk, M. (2022). Physical and Numerical Simulations for Predicting Distribution of Microstructural Features during Thermomechanical Processing of Steels. Materials, 15(5), 1660. https://doi.org/10.3390/ma15051660

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop