3.1. Acquisition of Experimental Stress-Strain Data
The specific chemical composition (wt.%) of SAE 5137H steel, as provided through the manufacturer, was listed in
Table 1. The initial microstructure with an average grain size of 64.8 µm was exhibited in
Figure 3. The samples were cylinders from rolled billet with a diameter of 10 mm and a height of 12 mm. The stress-strain data of this steel was acquired from a type of isothermal compression test at a certain temperature and a certain strain rate. A computer-controlled servo-hydraulic thermos-mechanical machine, i.e., Gleeble 3500, with a stable heating system, was adopted for this isothermal compression. The experiment procedures were roughly shown in
Figure 4. First of all, tantalum sheets were padded at both ends of the specimen before the experiment to minimize the friction between the edges of the specimen and the dies. Then, a specimen was heated to the proposal temperature at a rate of 10 K/s and held at a fixed temperature for 3 min. The holding treatment was to ensure that specimen achieved uniform temperature at the onset of a compression process, which could reduce the anisotropy of the material in flow deformation behavior. In addition, two thermocouple wires were welded to two points in the middle of each specimen before the test and the two wires will feedback the changes of temperature of the specimen at any time to the control system of Gleeble 3500, which is convenient for the system to adjust Itself and accurately control the temperature during text. Referring to the actual range of forging process parameters used in the factory, the temperatures of the compression test were selected as 1123 K, 1213 K, 1303 K, 1393 K, and 1483 K, and the strain rates were selected as 0.01 s
−1, 0.1 s
−1, 1 s
−1, and 10 s
−1. According to the material properties of SAE 5137H, the compression ratio of height should reach 60% during the test. Finally, the deformed specimens were immediately quenched with water to room temperature to preserve the elevated temperature microstructures. More details related to the isothermal compression procedures can be obtained in references [
26].
The true stress-strain curves of SAE 5137H steel during the compression processes under different deformation conditions are shown as
Figure 5. Considering adiabatic heating during the deformation at high strain rates may significantly influence the true stress–true strain curves, the reliability of the compression experimental results was judged using the expansion coefficient in Equation (1). When the expansion coefficient
B > 0.9 is considered reliable. On the contrary, it is not reliable, and the stress value obtained is corrected with Equation (2). After the test, expansion coefficient of all specimens after this isothermal compression test were greater than 0.9. Therefore, the results obtained from this experiment were reliable and can be used directly.
where
B is expansion coefficient,
L0 and
Lf are the original height of the specimen and the height after compression deformation, respectively, d
0 and d
f are the original diameter of the specimen and the diameter after compression deformation, respectively.
where
σi is corrected true stress,
Fi,
di,
Li, and
are the pressure on the specimen, the average diameter, the average height, and the friction coefficient, respectively.
The variation in flow stress with true strain can be summarized in three stages. In the first stage, the flow stress increases rapidly to a critical value with the true strain increases, and work hardening (WH) is the predominant deformation mechanism. At the same time, the grain boundary storage energy increases rapidly to the dynamic recrystallization activation energy (DRX). In the second stage, DRX and dynamic recovery (DRV) occur and increase, and the rate of increase of flow stress decreases until the maximum stress is reached. At this point, thermal softening begins to exceed WH. In the last stage, two types of stress change patterns are observed. Flow stress continues to decrease and the DRX softens obviously at strain rates of 0.01 and 0.1. Different from the first type, true stress-strain curves at strain rates of 1 and 10 do not show significant stress peaks. With the increase of true strain, the true stress basically remains stable due to the dynamic equilibrium of the softening behaviors of WH and DRV. At higher temperatures and lower strain rates, the flow stress is relatively lower due to sufficient time for recrystallization gain for energy accumulation and nucleation. In addition, high temperature accelerates dislocation movement and grain boundary migration. The flow behaviors are nonlinear with complex deformation mechanism. Therefore, it is important to establish a constitutive model to characterize the flow behaviors of SAE 5137H steel.
3.2. Semi-Physical Model with Improved Arrhenius-Type for Flow Behavior of SAE 5137H Steel
As for the Arrhenius-type constitutive equation, the effects of the temperatures and strain rates on the deformation behaviors are represented by Zener-Hollomon parameter,
Z, in an exponent-type expression.
where
is strain rate (s
−1).
Q is the activation energy of hot deformation (kJ·mol
−1),
A is material constant,
R is the universal gas constant (8.31 J·mol
−1·K
−1),
T is the absolute temperature (K),
is the function of stress, it takes three forms, and it is expressed as Equation (4),
where
,
, and
are the material constants,
.
is the flow stress (MPa) for a given strain.
For the low stress level
, substituting
into Equation (3), respectively, gives (5). Similarly for the high stress level
, substituting
into Equation (3), gives Equation (6).
Taking the logarithm on both sides of Equations (5) and (6), we have,
Consequently, the slope of
versus
and
versus
gives the value of
and
. This means
,
, respectively. Then, substituting the values of the peak stress and corresponding strain rates into the logarithm Equation (7) gives the relationships of
as shown in
Figure 6a. It is not difficult to find that the natural logarithms of stresses at every temperature are linear and the slopes are approximated the same with each other. The average value of all the lines’ slopes can be regarded as the inverse of
, thus
. Meanwhile, substituting the values of the peak stress and corresponding strain rates into the logarithm Equation (8) gives the relationships of
as shown in
Figure 6b. For the same, stresses at every temperature are linear and the average value of all the lines’ slopes can be regarded as the inverse of
, thus
. Then,
.
Substituting the hyperbolic law of
into Equation (2) gives,
Taking the natural logarithm of both sides of Equation (7) gives,
By linear fit, Equation (10) can be rewritten as,
As shown in
Figure 7a, the distribution of all points is linear. By a linear fitting with an average error of 0.10, the relationships between
and 1/
T is linear at different strain rate, and the slopes are approximated the same with each other. The value of
Q can be obtained from the slope of
versus 1/
T. The average value of all the slope rates is accepted
Q/(
RT), furthermore, the value of
Q is obtained as 359.6095 kJ·mol
−1.
Equation (10) can also be expressed as following:
By substituting the values of the peak stress at different temperatures and strain rates into Equation (12), the linear relationships between
and
for different temperatures can be obtained as shown in
Figure 7b. The average value of all the intercepts of
versus
plot is obtained as the value of
A, thus
A value is calculated as 6.9633 × 10
13 s
−1. Submitting the values of material constants
,
,
Q, and
A into Equation (9) gives,
Substituting Equation (3) into Equation (4), then the flow stress can be expressed as Equation (14).
Substituting Equation (3) and
A into Equation (14), the constitutive equation of flow stress for SAE 5137H steel can be calculated as Equation (15).
Substituting the polynomial functions of
,
,
and
into Equation (9), and gives Equation (17).
Finally, the Arrhenius type equation of SAE 5137H steel can be developed as following:
The above is the calculation process of the Arrhenius-type constitutive equation for SAE 5137H steel, while this equation ignores the effect of strain on the flow stress, and then this equation is lack of the ability to predict the stresses at different strains. In order to solve this issue, the strain compensation is introduced by constructing a series of polynomials as Equation (16) representing the nonlinear relationships between the variables (including activation energy of deformation
Q, material constants
n, and
α, and structure factor
A) in Arrhenius-type constitutive equation and strains. In order to find the variation pattern of the variables, the values of the variables were fitted nonlinearly at 0.1 true strain interval. Such nonlinear relationships were shown as
Figure 8, and the coefficients of the fitted polynomials were listed in
Table 2.
3.3. BP-ANN for Flow Behavior of SAE 5137H Steel
In this investigation, the BP-ANN was developed by MATLAB software. The input variables include deformation temperatures and strains, and the output variables were flow stresses. The 20 curves were divided into two datasets, i.e., the training dataset and the test dataset, as shown in
Table 3. A total of 308 input-output pairs were selected from the stress-strain curves to train and test the BP-ANN. The 36 stress points on the test stress-strain curves in the strain range 0.075~0.875 with a distance of 0.1 were not used for training, but for testing the BP-ANN generation ability. The BP-ANN was trained using 272 stress points in the strain range of 0.05~0.85 and distance of 0.5 in the training stress-strain curves. Due to the different units of measurement from experimental data such as temperatures, strains, strain rates, and stresses, there were large differences between different types of data, and such differences would reduce the speed and accuracy of convergence within the network. Therefore, the input and output datasets measured in different units need to be normalized to dimensionless units before the networks are trained to eliminate the arbitrary effect of similarity between different data. The input and output data were normalized in the range of 0~1 according to the relationship given in Equation (19).
where
is the normalized value of
,
is the experimental data,
and
are the maximum and minimum value of
respectively.
As mentioned above, for a typical BP-ANN structure, one or more hidden layers were required, and two hidden layers were used here to ensure a high training accuracy. In addition, the number of nodes in the input layer was 3, the number of nodes in the hidden layer was 12, and the number of nodes in the output layer was 1. Considering the range of values, the tansig was used as hidden layers function, and the purelin was used as the output layer function. Other parameters were shown in
Table 4.