Deformation Behavior and Constitutive Model of 34CrNi3Mo during Thermo-Mechanical Deformation Process

With higher creep strength and heat resistance, 34CrNi3Mo has been widely used in the production of engine rotors, steam turbine impellers, and turbine blades. To investigate the hot deformation behaviors of 34CrNi3Mo steel, hot compressive tests were conducted on a Gleeble-3500 thermomechanical simulator, under the temperature range of 1073 K–1373 K and strain rate ranges of 0.1 s−1–20 s−1. The results show that the flow stress of 34CrNi3Mo steel under high temperatures is greatly influenced by the deformation temperature and strain rate, and it is the result of the interaction between strain hardening, dynamic recovery, and recrystallization. Under the same deformation rate, as the deformation temperature increases, the softening effect of dynamic recrystallization and dynamic recovery gradually increases, and the flow stress gradually decreases. Under the same deformation temperature, with the increase of strain rate, the influence of strain hardening on 34CrNi3Mo steel is gradually in power, and the flow stress gradually increases. To predict the flow stress of 34CrNi3Mo steel accurately, a modified Arrhenius-type constitutive model considering the effects of strain, temperature, and strain rate at the same time was made based on the experiment data. On this basis, the evolution law of deformation activation and instability characteristics of 34CrNi3Mo steel were investigated, and the processing map of 34CrNi3Mo steel was established. The formability of 34CrNi3Mo steel under high temperature deformation was revealed, which provided a theoretical foundation of the equation of reasonable hot working process.


Introduction
With high creep strength and high heat resistance, 34CrNi3Mo steel is a kind of high-strength alloy structural steel widely used in the production of engine rotors, steam turbine impellers, and turbine blades, and these important force transmission parts must be produced by forging to meet its application requirements. In the process of forging, the rheological properties and forming properties of materials at high temperature are the preconditions for working out a reasonable forging process. Therefore, it is important to establish the flow stress model of 34CrNi3Mo at high temperature to reveal its deformation characteristics and forming properties, which is of guiding significance for formulating a reasonable forging process.
A reasonable flow stress model can accurately describe the influence of forming process parameters on the high temperature deformation conditions, such as deformation temperature, deformation amount, and deformation rate, which is also a prerequisite for accurately predicting and analyzing the deformation law of materials in the plastic forming process [1,2]. The macroscopic phenomenological constitutive model based on mechanics of continuous media and thermodynamics describes the mechanical properties of materials by simulating and fitting the macroscopic phenomena of materials in the deformation process [3]. It cannot reveal the characteristics and evolution mechanism of material deformation from an imperceptible and microscopic perspective, however, the modeling and

Experiment Method
In this paper, 34CrNi3Mo hot-rolled bar with a diameter of 180 mm was taken as the research object. In order to eliminate the difference in material properties caused by the difference in internal microstructure, a compression sample with a diameter of 10 mm and a height of 15 mm were cut at the position on the outer surface of the bar, as shown in Figure 1, under the deformation conditions of four different temperatures (1073 K, 1223 K, 1323 K, 1373 K) and four different strain rates (0.1 s −1 , 1 s −1 , 10 s −1 , 20 s −1 ). In the experiment process, the amount of total strain of the compressed sample was set as 0.65, and the isothermal compression tests were carried out on 34CrNi3Mo cylindrical specimens with the help of Gleebel-3500 thermal simulation test machine, as shown in Figure 2.

Experiment Method
In this paper, 34CrNi3Mo hot-rolled bar with a diameter of 180 mm was taken as the research object. In order to eliminate the difference in material properties caused by the difference in internal microstructure, a compression sample with a diameter of 10 mm and a height of 15 mm were cut at the position on the outer surface of the bar, as shown in Figure 1, under the deformation conditions of four different temperatures (1073 K, 1223 K, 1323 K, 1373 K) and four different strain rates (0.1 s −1 , 1 s −1 , 10 s −1 , 20 s −1 ). In the experiment process, the amount of total strain of the compressed sample was set as 0.65, and the isothermal compression tests were carried out on 34CrNi3Mo cylindrical specimens with the help of Gleebel-3500 thermal simulation test machine, as shown in Figure 2.  In forging process, the 34CrNi3Mo workpiece is usually heated to the initial forging temperature (1473 K). However, due to the heat exchange among the workpiece, the mold, and the environment, the temperature of the workpiece gradually decreases. In order to simulate the effect of material temperature change on the properties during forging process, the samples were heated to the temperature T0 = 1473 K, then cooled to the deformation temperature and kept there before deformation for 180 s to make the sample temperature uniform.

Results and Analysis
The true stress-strain curves of 34CrNi3Mo steel at different deformation conditions were obtained through isothermal compression tests, as shown in Figure 3. It can be seen

Experiment Method
In this paper, 34CrNi3Mo hot-rolled bar with a diameter of 180 mm was taken as the research object. In order to eliminate the difference in material properties caused by the difference in internal microstructure, a compression sample with a diameter of 10 mm and a height of 15 mm were cut at the position on the outer surface of the bar, as shown in Figure 1, under the deformation conditions of four different temperatures (1073 K, 1223 K, 1323 K, 1373 K) and four different strain rates (0.1 s −1 , 1 s −1 , 10 s −1 , 20 s −1 ). In the experiment process, the amount of total strain of the compressed sample was set as 0.65, and the isothermal compression tests were carried out on 34CrNi3Mo cylindrical specimens with the help of Gleebel-3500 thermal simulation test machine, as shown in Figure 2.  In forging process, the 34CrNi3Mo workpiece is usually heated to the initial forging temperature (1473 K). However, due to the heat exchange among the workpiece, the mold, and the environment, the temperature of the workpiece gradually decreases. In order to simulate the effect of material temperature change on the properties during forging process, the samples were heated to the temperature T0 = 1473 K, then cooled to the deformation temperature and kept there before deformation for 180 s to make the sample temperature uniform.

Results and Analysis
The true stress-strain curves of 34CrNi3Mo steel at different deformation conditions were obtained through isothermal compression tests, as shown in Figure 3. It can be seen  In forging process, the 34CrNi3Mo workpiece is usually heated to the initial forging temperature (1473 K). However, due to the heat exchange among the workpiece, the mold, and the environment, the temperature of the workpiece gradually decreases. In order to simulate the effect of material temperature change on the properties during forging process, the samples were heated to the temperature T 0 = 1473 K, then cooled to the deformation temperature and kept there before deformation for 180 s to make the sample temperature uniform.

Results and Analysis
The true stress-strain curves of 34CrNi3Mo steel at different deformation conditions were obtained through isothermal compression tests, as shown in Figure 3. It can be seen from the test results that the flow stress of 34CrNi3Mo steel under high temperature deformation is the result of the interaction of strain hardening, dynamic recovery, and recrystallization. The flow stress of 34CrNi3Mo steel exhibits obvious dynamic recrystallization under high temperature deformation, especially under high temperature and low strain rate deformation. In the initial stage of deformation, the flow stress increases gradually with the increase of the strain; when the stress reaches a peak value, the flow stress decreases with the increase of the strain. from the test results that the flow stress of 34CrNi3Mo steel under high temperature deformation is the result of the interaction of strain hardening, dynamic recovery, and recrystallization. The flow stress of 34CrNi3Mo steel exhibits obvious dynamic recrystallization under high temperature deformation, especially under high temperature and low strain rate deformation. In the initial stage of deformation, the flow stress increases gradually with the increase of the strain; when the stress reaches a peak value, the flow stress decreases with the increase of the strain. Under different compression temperature conditions, the flow stress of 34CrNi3Mo steel is greatly affected by the deformation temperature at the same strain rate. With the increase of deformation temperature, the dislocation activity and the kinetic energy of atoms increases at the same time. These microscopic changes will make dislocation slip more easily. Therefore, with the increase of temperature, the softening effect of dynamic recrystallization and dynamic recovery of 34CrNi3Mo steel gradually increases, and the flow stress gradually decreases.
At the same deformation temperature, the dislocation proliferation rate caused by strain hardening increases gradually with the increase of strain rate, while the rate of dislocation disappearance caused by dynamic recovery gradually decreases. Although dynamic recrystallization is easier to nucleate under high dislocation density and high temperature, dynamic recrystallization is not enough to eliminate strain hardening caused by increased dislocation density because of too fast strain rate to short deformation time. Therefore, under the condition of high temperature deformation, with the increase of strain rate, the influence of strain hardening on 34CrNi3Mo steel gradually dominates, and the flow stress gradually increases.  Under different compression temperature conditions, the flow stress of 34CrNi3Mo steel is greatly affected by the deformation temperature at the same strain rate. With the increase of deformation temperature, the dislocation activity and the kinetic energy of atoms increases at the same time. These microscopic changes will make dislocation slip more easily. Therefore, with the increase of temperature, the softening effect of dynamic recrystallization and dynamic recovery of 34CrNi3Mo steel gradually increases, and the flow stress gradually decreases.
At the same deformation temperature, the dislocation proliferation rate caused by strain hardening increases gradually with the increase of strain rate, while the rate of dislocation disappearance caused by dynamic recovery gradually decreases. Although dynamic recrystallization is easier to nucleate under high dislocation density and high temperature, dynamic recrystallization is not enough to eliminate strain hardening caused by increased dislocation density because of too fast strain rate to short deformation time. Therefore, under the condition of high temperature deformation, with the increase of strain rate, the influence of strain hardening on 34CrNi3Mo steel gradually dominates, and the flow stress gradually increases.

Arrhenius-Type Constitutive Model
It can be seen from the thermal compression test results that the flow stress of 34CrNi3Mo steel is greatly affected by strain, strain rate, and temperature under thermal deformation conditions. In order to accurately describe the flow stress under thermal deformation conditions, the relationship between stress and strain, temperature, and strain rate must be established. The Arrhenius-type constitutive model can characterize the relationship between the flow stress σ, strain rate . ε, and deformation temperature T when materials are deformed at high temperature, where A is the material parameter, R = 8.31J/mol·K, and f (σ) is the stress function. Under the condition of low stress levels, namely, ασ < 0.8, f (σ) = σ n 1 . Where α is the stress level parameter (MPa −1 ). Then, substituted into Equation (1), the relationship between flow stress σ and strain rate . ε can be approximately shown as, where n 1 is the material parameter. Under the condition of high stress levels, namely, ασ > 1.2, f (σ) = exp(βσ), the flow stress σ and strain rate . ε can be expressed as the Equation as follows.
where β is the constant relevant to the materials, and β = α·n 1 . Under all stress levels conditions, f (σ) = [sin h(ασ)] n , n is the material parameter. The Arrhenius-type hyperbolic sine form of constitutive equation containing deformation activation energy Q and deformation temperature T can be obtained, By taking logarithms of both sides of Equations (2) and (3), we can see that, From Equations (5) and (6), it can be seen that when the temperature is fixed, the ln . ε satisfies the linear equation with σ and lnσ, respectively, namely both ln  (5) and (6) is the stress index n 1 and the material constant β. Therefore, under certain strain conditions, the test results under the same deformation temperature and different strain rates were linearly fitted according to Equation (5) and Equation (6), respectively, then the stress index n 1 and material constant β under different stress levels can be obtained.
The experimental results of plastic strain ε = 0.3 were selected as the research object, and the experimental results were linearly fitted according to Equation (5) and Equation (6), respectively. The linear fitting results of lnσ − ln  Table 1.   The deformation activation energy Q can be expressed as, It can be seen from Equation (7) that at a certain deformation temperature T, the slope of the curve ln − ln sinh( ) is the material parameter n. When the strain rate is fixed, the slope of the curve ln sinh( ) − ln(1/ ) is material parameter K, and K can be expressed as = ⁄ .
Combining the obtained value α and the experimental results obtained under different deformation conditions, the linear fitting results of − sinh( ) and ln sinh( ) − ln (1 ) ⁄ can be plotted, as shown in Figures 6 and 7. Thus, the values of n under different temperatures can be calculated according to Figures 6 and 7, which is shown in Table 2, and the deformation activation energy Q at different strain rates is shown in Table 3.    The deformation activation energy Q can be expressed as, It can be seen from Equation (7) that at a certain deformation temperature T, the slope of the curve ln − ln sinh( ) is the material parameter n. When the strain rate is fixed, the slope of the curve ln sinh( ) − ln(1/ ) is material parameter K, and K can be expressed as = ⁄ .
Combining the obtained value α and the experimental results obtained under different deformation conditions, the linear fitting results of − sinh( ) and ln sinh( ) − ln (1 ) ⁄ can be plotted, as shown in Figures 6 and 7. Thus, the values of n under different temperatures can be calculated according to Figures 6 and 7, which is shown in Table 2, and the deformation activation energy Q at different strain rates is shown in Table 3.  Since Equation (5) is obtained under the condition of a low stress level, Equation (6) is obtained under the condition of a high stress level. Therefore, we take the average value of the four lines as the value of n 1 in the research process of this paper, the value n 1 = 10.231 can be obtained through calculation. In a similar way, taking the average of the four lines as the value of β, thus β = 0.0662. According to the relationship between the material parameters n 1 and β, β = α·n 1 , the stress level parameter α=0.0064MPa −1 can be obtained.
The deformation activation energy Q can be expressed as, It can be seen from Equation (7) that at a certain deformation temperature T, the slope of the curve ln . ε − ln[sin h(ασ)] is the material parameter n. When the strain rate . ε is fixed, the slope of the curve ln[sin h(ασ)] − ln(1/T) is material parameter K, and K can be expressed as K = Q/nR.
Combining the obtained value α and the experimental results obtained under different deformation conditions, the linear fitting results of ln  Table 2, and the deformation activation energy Q at different strain rates is shown in Table 3.    When the materials produce plastic deformation under high temperature conditions, the relationship between strain rate, deformation activation energy, and deformation temperature can be expressed by the Zener-Hollomon parameter, namely, Z parameter representing the strain rate factor of temperature compensation.
Taking the logarithm of both sides of Equation (8), the Equation (8) will be changed as, According to Equation (9), lnA is the intercept of the linear curve of ln − ln sinh( ) . By substituting the value of Q into the Equation (8)    When the materials produce plastic deformation under high temperature conditions, the relationship between strain rate, deformation activation energy, and deformation temperature can be expressed by the Zener-Hollomon parameter, namely, Z parameter representing the strain rate factor of temperature compensation.
Taking the logarithm of both sides of Equation (8), the Equation (8) will be changed as, According to Equation (9), lnA is the intercept of the linear curve of ln − ln sinh( ) . By substituting the value of Q into the Equation (8)  When the materials produce plastic deformation under high temperature conditions, the relationship between strain rate, deformation activation energy, and deformation temperature can be expressed by the Zener-Hollomon parameter, namely, Z parameter representing the strain rate factor of temperature compensation. Taking the logarithm of both sides of Equation (8), the Equation (8) will be changed as, According to Equation (9), lnA is the intercept of the linear curve of ln Z − ln[sin h(ασ)]. By substituting the value of Q into the Equation (8), the Z parameters under different deformation conditions can be obtained. By linear fitting the relationship between ln Z and ln[sin h(ασ)] under different deformation conditions, as shown in Figure 8, the material parameter A can be obtained, A = 3.51 × 10 16 .  According to the experimental results of 34CrNi3Mo steel, it can be seen that the flow stress curve of 34CrNi3Mo steel is not only related to the deformation temperature and strain rate but is also affected by the plastic strain. Therefore, all the parameters in the Arrhenius-type model will be affected by plastic strain. With the same method as above, a strain value is selected every 0.05 within the strain range of 0.1-0.6, and the values of model parameters under different strain conditions can be calculated, as shown in Table 4. In this paper, the sixth-degree polynomial was used to fit model parameters under different deformation conditions, the fitting curve is shown in Figure 9, and the functional relationship between model parameters and strain obtained through fitting is shown in Equation (10). The coefficients of each material parameter model in Equation (10) are shown in Table 5.
As shown in Table 4 and Figure 9, in the initial stage of plastic deformation, the stress level parameter α in the Arrhenius-type constitutive model decreases rapidly with the increase of deformation. When the strain is greater than 0.3, the stress level parameter α tends to be stable with the increase of deformation. When the strain is less than 0.2, the material parameter n, the deformation activation energy Q, and the parameter A all increase with the increase of deformation. However, the strain is greater than 0.2 and less than 0.5, the material parameters n, A, and the deformation activation energy Q gradually decrease with the increase of deformation. This change rule will not exist when the strain is greater than 0.5, the three parameters increase gradually with the increase of strain.  According to the experimental results of 34CrNi3Mo steel, it can be seen that the flow stress curve of 34CrNi3Mo steel is not only related to the deformation temperature and strain rate but is also affected by the plastic strain. Therefore, all the parameters in the Arrhenius-type model will be affected by plastic strain. With the same method as above, a strain value is selected every 0.05 within the strain range of 0.1-0.6, and the values of model parameters under different strain conditions can be calculated, as shown in Table 4. In this paper, the sixth-degree polynomial was used to fit model parameters under different deformation conditions, the fitting curve is shown in Figure 9, and the functional relationship between model parameters and strain obtained through fitting is shown in Equation (10). The coefficients of each material parameter model in Equation (10) are shown in Table 5.   According to Equation (4), the flow stress can be expressed as, By substituting the stress level parameters α(ε), the material parameters n(ε), A(ε), and the deformation activation energy Q(ε) into Equation (11), the prediction results of flow stress under the different conditions can be obtained, as shown in Figure 10.  As shown in Table 4 and Figure 9, in the initial stage of plastic deformation, the stress level parameter α in the Arrhenius-type constitutive model decreases rapidly with the increase of deformation. When the strain is greater than 0.3, the stress level parameter α tends to be stable with the increase of deformation. When the strain is less than 0.2, the material parameter n, the deformation activation energy Q, and the parameter A all increase with the increase of deformation. However, the strain is greater than 0.2 and less than 0.5, the material parameters n, A, and the deformation activation energy Q gradually decrease with the increase of deformation. This change rule will not exist when the strain is greater than 0.5, the three parameters increase gradually with the increase of strain.
According to Equation (4), the flow stress can be expressed as, By substituting the stress level parameters α(ε), the material parameters n(ε), A(ε), and the deformation activation energy Q(ε) into Equation (11), the prediction results of flow stress under the different conditions can be obtained, as shown in Figure 10. By comparing the predicted flow stress curve shown in Figure 10 with the experimental results, the Arrhenius-type constitutive model considering the effect of strain can accurately predict the flow stress of 34CrNi3Mo steel when the deformation temperature is lower (less than 1223 K). However, with the increase of deformation temperature, the interaction between dynamic recrystallization and strain rate hardening is gradually obvious, under the condition of low strain rate, the softening phenomenon of the Arrheniustype constitutive model considering strain only is not obvious, and the strain rate hardening effect cannot be reflected under the higher strain rate. Therefore, there is a large error between model prediction results and experimental results at high temperature.

Modified Arrhenius-Type Constitutive Model
It can be seen from the above research that, in order to accurately predict the flow stress of 34CrNi3Mo steel under high temperature deformation with the help of the Arrhenius-type constitutive model, the interaction between strain, deformation temperature, and strain rate on model parameters must be considered simultaneously. Therefore, the modified Arrhenius-type model can be rewritten as [18], where the material parameter n and the stress level parameter α are the functions of strain and deformation temperature, namely, ( , ), ( , ), the deformation activation energy Q and the material parameter A are the functions of strain, strain rate, and deformation temperature, namely, ( , , ), (ε, , ). By comparing the predicted flow stress curve shown in Figure 10 with the experimental results, the Arrhenius-type constitutive model considering the effect of strain can accurately predict the flow stress of 34CrNi3Mo steel when the deformation temperature is lower (less than 1223 K). However, with the increase of deformation temperature, the interaction between dynamic recrystallization and strain rate hardening is gradually obvious, under the condition of low strain rate, the softening phenomenon of the Arrhenius-type constitutive model considering strain only is not obvious, and the strain rate hardening effect cannot be reflected under the higher strain rate. Therefore, there is a large error between model prediction results and experimental results at high temperature.

Modified Arrhenius-Type Constitutive Model
It can be seen from the above research that, in order to accurately predict the flow stress of 34CrNi3Mo steel under high temperature deformation with the help of the Arrheniustype constitutive model, the interaction between strain, deformation temperature, and strain rate on model parameters must be considered simultaneously. Therefore, the modified Arrhenius-type model can be rewritten as [18], where the material parameter n and the stress level parameter α are the functions of strain and deformation temperature, namely, n(ε, T), α(ε, T), the deformation activation energy Q and the material parameter A are the functions of strain, strain rate, and deformation temperature, namely, A ε, ε, T . The value of stress level parameter α under different deformation temperatures and strains are shown in Figure 11a. It can be seen from the figure that stress level parameters are greatly affected by deformation temperature. After undergoing large plastic deformation, strain has little influence on stress level parameters. Through nonlinear surface fitting, as shown in Figure 11b, the function of stress level parameter α on deformation temperature T and strain ε can be expressed as follows, α(ε, T) = 0.03419 − 0.01172ε − 0.000065T + 0.02498ε 2 + 3.7 × 10 −8 T 2 − 6.9 × 10 −6 εT (13) fitting, as shown in Figure 11b, the function of stress level parameter α on deformation temperature T and strain can be expressed as follows, ( , ) = 0.03419 − 0.01172 − 0.000065 + 0.02498 + 3.7 10 − 6.9 10 According to the Equation of the correlation coefficient, the correlation coefficient R 2 of the fitted surface is 0.951. This proves the accuracy of fitting formula.
According to Equation (14), n(ε,T) is the slope of curves ln -ln sin ℎ ( , ) under different strains and deformation temperatures. Therefore, according to the flow stress curve of 34CrNi3Mo steel, the curves of n(ε,T) under different strain and temperature conditions can be drawn, as shown in Figure 12a. It can be seen from Figure 12a that the value of n(ε,T) is greatly affected by strain and temperature. Under the same strain condition, the value of n(ε,T) decreases obviously with the increase of temperature. At the same deformation temperature, when the deformation temperature is greater than 1223 K, the value of n(ε,T) decreases with the increase of strain.
Similary, the function of material parameter n on deformation temperature T and strain ε was obtained by 3D surface fitting, as shown in Figure 12b, and the expression of fitting function is shown in Equation (15), with the correlation coefficient R 2 = 0.951. According to the Equation of the correlation coefficient, the correlation coefficient R 2 of the fitted surface is 0.951. This proves the accuracy of fitting formula.
Taking the logarithm of both sides of Equation (12), we can see that, According to Equation (14), n(ε,T) is the slope of curves ln . ε-under different strains and deformation temperatures. Therefore, according to the flow stress curve of 34CrNi3Mo steel, the curves of n(ε,T) under different strain and temperature conditions can be drawn, as shown in Figure 12a. It can be seen from Figure 12a that the value of n(ε,T) is greatly affected by strain and temperature. Under the same strain condition, the value of n(ε,T) decreases obviously with the increase of temperature. At the same deformation temperature, when the deformation temperature is greater than 1223 K, the value of n(ε,T) decreases with the increase of strain.
K, the value of n(ε,T) decreases with the increase of strain.
Similary, the function of material parameter n on deformation temperature T and strain ε was obtained by 3D surface fitting, as shown in Figure 12b, and the expression of fitting function is shown in Equation (15), with the correlation coefficient R 2 = 0.951. Similary, the function of material parameter n on deformation temperature T and strain ε was obtained by 3D surface fitting, as shown in Figure 12b, and the expression of fitting function is shown in Equation (15), with the correlation coefficient R 2 = 0.951. n(ε, T) = 7.66322 + 14.48418ε + 0.00582T + 1.68993ε 2 − 0.00000446633T 2 − 0.01402εT (15) Under the certain strain and deformation rate, Equation (14) can be rewritten as, According to Equation (16), it can be seen that the deformation activation energy Q is the slope of the curve of ln{sin h[α(ε)σ]} − 1/T.
According to Equation (17), the deformation activation energy Q ε, . ε, T consists of two functions, n(ε, T) and K ε, . ε . The function of n(ε, T) can be obtained according to Equation (15), and the K ε, . ε is a function of strain and strain rate. According to the experimental results, the value of K ε, . ε under different strains and strain rates can be obtained through calculation, the calculation results are shown in Figure 13a. It can be seen from the figure that the strain and strain rate have great influence on the value of K ε, . ε . Under the same strain condition, the value of K ε, . ε decreases gradually with the increase of strain rate. When the strain rate is small ( ε decreases with the increase of deformation (ε > 0.1). However, when the strain rate is large ( . ε ≥ 10 s −1 ), the value of K ε, . ε increases gradually with the increase of deformation (ε > 0.1). The functional expression of the function of K ε, . ε can be obtained through 3D surface fitting, as shown in Figure 13b. The fitting function is shown in Equation (18), and the correlation coefficient R 2 of the fitting results is 0.937.
By substituting Equation (18) and Equation (15) into Equation (17), the functional expression of the deformation activation energy Q ε, . ε, T considering the effects of strain, strain rate, and deformation temperature can be obtained. tion (18), and the correlation coefficient R 2 of the fitting results is 0.937.
(ε, ) = 5.99546 + 3.06765ε − 0.16104ln − 2.92598ε2 + 0.04702(ln )2 − 0.48355ε(ln ) (18) By substituting Equation (18) and Equation (15) into Equation (17), the functional expression of the deformation activation energy (ε, , ) considering the effects of strain, strain rate, and deformation temperature can be obtained. According to Equation (14), assuming the function value of ( , ) is equal to the function value of lnA( , , ) − (ε, , )⁄ , thus, the value of ( , ) is the intercept of the curve lnln sinh ( , ) under the different deformation temperatures and strains. Therefore, the value of ( , ) could be calculated according to the experimental results, as shown in Figure 14a. It can be seen from the figure that under the same strain condition, the value of ( , ) increases gradually with the increase of deformation temperature. However, at the same deformation temperature, the variation of ( , ) value is just opposite to the previous rule, the value of ( , ) decreases with the increase of strain. The function of ( , ) can be obtained through 3D surface fitting, as shown in  Figure 14a. It can be seen from the figure that under the same strain condition, the value of b(ε, T) increases gradually with the increase of deformation temperature. However, at the same deformation temperature, the variation of b(ε, T) value is just opposite to the previous rule, the value of b(ε, T) decreases with the increase of strain. The function of b(ε, T) can be obtained through 3D surface fitting, as shown in Figure 14b. The expression of fitting function b(ε, T) is shown in Equation (19), with the correlation coefficient By substituting Equations (19) and (17)  Using a similar method as in Equation (4), the rheological prediction model considering the interaction between the stress with plastic strain, temperature, and strain rate can be obtained by reordering and transforming the Equation (12), the prediction model is shown in Equation (21), where the parameters of the model could be calculated through Equations (13) and (17) According to the constitutive model of Equation (21), the predicted results of the flow stress of 34CrNi3Mo steel under the high temperature condition are shown in Figure 15. Comparing the predicted results according to Equation (21) with the experimental results, as shown in Figure 15 By substituting Equations (19) and (17) Using a similar method as in Equation (4), the rheological prediction model considering the interaction between the stress with plastic strain, temperature, and strain rate can be obtained by reordering and transforming the Equation (12), the prediction model is shown in Equation (21), where the parameters of the model could be calculated through Equations (13) and (17)- (20).
According to the constitutive model of Equation (21), the predicted results of the flow stress of 34CrNi3Mo steel under the high temperature condition are shown in Figure 15. Comparing the predicted results according to Equation (21) with the experimental results, as shown in Figure 15, the predicted results have a large error with the experimental results. For exploring the causes of such errors, the effects of material parameters n(ε, T), A ε, . ε, T , stress level parameter α(ε, T), and deformation activation energy Q ε, . ε, T on the constructive model (as shown in Equation (21)) were investigated. Through comparison research, it was found that stress level parameters α(ε, T) were the main reasons for the large error. When the stress level parameter is only related to strain, and other parameters function of Equation (21) remain the same, that is to say, in the calculation of the flow stress by using Equation (21), just replace the function of level parameter ( , ) with ( ), the expression of the ( ) satisfies Equation (10). Through the transformation mentioned above, the predicted results of flow stress are shown in Figure 16. It can be seen that the modified model's predicted results are highly consistent with the experimental results.
Therefore, the modified Arrhenius-type constitutive model will be rewritten as, where the model parameters can be expressed as follow,  When the stress level parameter is only related to strain, and other parameters function of Equation (21) remain the same, that is to say, in the calculation of the flow stress by using Equation (21), just replace the function of level parameter α(ε, T) with α(ε), the expression of the α(ε) satisfies Equation (10). Through the transformation mentioned above, the predicted results of flow stress are shown in Figure 16. It can be seen that the modified model's predicted results are highly consistent with the experimental results.
Therefore, the modified Arrhenius-type constitutive model will be rewritten as, where the model parameters can be expressed as follow,  In order to quantificationally analyze the accuracy of the Arrhenius-type model and the modified Arrhenius-type model established in this paper, the relative error curves of the prediction results of these two models and the test results were obtained, as shown in Figure 17. It can be seen from the relative error results that the predicted accuracy of the modified Arrhenius-type model is significantly higher than the Arrhenius-type model. Therefore, the modified Arrhenius-type model, which comprehensively considers the effects of strain, temperature, and strain rate on model parameters, can better predict the flow stress of 34CrNi3Mo steel under high temperature deformation.

Evolution of Activation Energy
The evolution of activation energy Q of 34CrNi3Mo steel at different temperatures and strain rates was calculated by the Arrhenius-type model, as shown in Figure 18. It can be seen from the figure that the evolution of activation energy of 34CrNi3Mo steel varies from 242 kJ·mol −1 to −607 kJ·mol −1 under different deformation conditions and decreases with the increase of deformation temperature and strain rate. This indicates that the dislocation motion of 34CrNi3Mo steel is promoted by temperature under the high temper- In order to quantificationally analyze the accuracy of the Arrhenius-type model and the modified Arrhenius-type model established in this paper, the relative error curves of the prediction results of these two models and the test results were obtained, as shown in Figure 17. It can be seen from the relative error results that the predicted accuracy of the modified Arrhenius-type model is significantly higher than the Arrhenius-type model. Therefore, the modified Arrhenius-type model, which comprehensively considers the effects of strain, temperature, and strain rate on model parameters, can better predict the flow stress of 34CrNi3Mo steel under high temperature deformation. In order to quantificationally analyze the accuracy of the Arrhenius-type model and the modified Arrhenius-type model established in this paper, the relative error curves of the prediction results of these two models and the test results were obtained, as shown in Figure 17. It can be seen from the relative error results that the predicted accuracy of the modified Arrhenius-type model is significantly higher than the Arrhenius-type model. Therefore, the modified Arrhenius-type model, which comprehensively considers the effects of strain, temperature, and strain rate on model parameters, can better predict the flow stress of 34CrNi3Mo steel under high temperature deformation.

Evolution of Activation Energy
The evolution of activation energy Q of 34CrNi3Mo steel at different temperatures and strain rates was calculated by the Arrhenius-type model, as shown in Figure 18. It can be seen from the figure that the evolution of activation energy of 34CrNi3Mo steel varies from 242 kJ·mol −1 to −607 kJ·mol −1 under different deformation conditions and decreases with the increase of deformation temperature and strain rate. This indicates that the dislocation motion of 34CrNi3Mo steel is promoted by temperature under the high temper-

Evolution of Activation Energy
The evolution of activation energy Q of 34CrNi3Mo steel at different temperatures and strain rates was calculated by the Arrhenius-type model, as shown in Figure 18. It can be seen from the figure that the evolution of activation energy of 34CrNi3Mo steel varies from 242 kJ·mol −1 to −607 kJ·mol −1 under different deformation conditions and decreases with the increase of deformation temperature and strain rate. This indicates that the dislocation motion of 34CrNi3Mo steel is promoted by temperature under the high temperature. At the same time, the softening effect of temperature is obviously enhanced with the increase of deformation temperature. This will cause the dislocation density and the activation energy to decrease. Under the condition of high strain rate, the high deformation rate increases the shear stress and is favorable to the dislocation movement, thus reducing the dislocation density.

Processing Map
According to the dynamic material model (DMM), the total energy absorbed by the workpiece during the thermal process can be expressed as, where P represents the external input energy, G represents the energy consumed by the plastic deformation, and J represents the energy dissipation caused by the transformation of microstructures, and these two energy distribution relationships are affected by the constitutive model. The proportion of the energy parameters G and J is determined by the strain rate sensitivity index m.
When the material meets the satisfied linear power dissipation system, the strain rate sensitivity index m = 1. On this condition, the energy dissipation J reaches the maximum value Jmax.
When m > 1, the material tends to form twins or shear bands during deformation. η is the energy dissipation factor and its value is J/Jmax. The energy dissipation factor is used to represent the ratio of microstructure conversion dissipation energy to ideal linear dissipation energy. According to the definition of the energy dissipation factor η and the Murty criterion [26], the power dissipation coefficient can be obtained as follows,

Processing Map
According to the dynamic material model (DMM), the total energy absorbed by the workpiece during the thermal process can be expressed as, εdσ (24) where P represents the external input energy, G represents the energy consumed by the plastic deformation, and J represents the energy dissipation caused by the transformation of microstructures, and these two energy distribution relationships are affected by the constitutive model. The proportion of the energy parameters G and J is determined by the strain rate sensitivity index m.
When the material meets the satisfied linear power dissipation system, the strain rate sensitivity index m = 1. On this condition, the energy dissipation J reaches the maximum value J max .
When m > 1, the material tends to form twins or shear bands during deformation.
η is the energy dissipation factor and its value is J/J max . The energy dissipation factor is used to represent the ratio of microstructure conversion dissipation energy to ideal linear dissipation energy. According to the definition of the energy dissipation factor η and the Murty criterion [26], the power dissipation coefficient can be obtained as follows, When the constitutive model is exponential, the expression of energy dissipation factor η can be simplified as, The energy consumed in the thermal deformation process consists of two parts, one part is the energy consumed by the dynamic recovery and recrystallization of the material, and the other part is the energy consumed by the microscopic defects, such as holes, adiabatic shear bands, and wedge cracks. Therefore, the larger the value of the energy dissipation factor, the greater the change in microstructure. However, it is not the higher the energy dissipation rate, the better the thermoforming formability.
The flow instability criterion proposed by Prasad [25] can be used as the instability criterion for 34CrNi3Mo steel during thermal deformation, thus the instability factors ζ . ε can be expressed as, From Equation (29), it can seen that the function of strain rate and deformation temperature. When ζ . ε < 0, the fluidity instability will occur. Based on the experimental results and Equations (28) and (29), the energy dissipation diagram and plastic instability diagram can be drawn, respectively. Then, the processing map of 34CrNi3Mo steel can be obtained by superposition the energy dissipation diagram and the plastic instability diagram, as shown in Figure 19. In the processing map, the contour line represents the percentage of power dissipation efficiency, and the shaded part represents the plastic instability area.
The energy consumed in the thermal deformation process consists of two parts, one part is the energy consumed by the dynamic recovery and recrystallization of the material, and the other part is the energy consumed by the microscopic defects, such as holes, adiabatic shear bands, and wedge cracks. Therefore, the larger the value of the energy dissipation factor, the greater the change in microstructure. However, it is not the higher the energy dissipation rate, the better the thermoforming formability.
The flow instability criterion proposed by Prasad [25] can be used as the instability criterion for 34CrNi3Mo steel during thermal deformation, thus the instability factors ( ) can be expressed as, From Equation (29), it can seen that the function of strain rate and deformation temperature. When ( ) < 0, the fluidity instability will occur.
Based on the experimental results and Equations (28) and (29), the energy dissipation diagram and plastic instability diagram can be drawn, respectively. Then, the processing map of 34CrNi3Mo steel can be obtained by superposition the energy dissipation diagram and the plastic instability diagram, as shown in Figure 19. In the processing map, the contour line represents the percentage of power dissipation efficiency, and the shaded part represents the plastic instability area.
Comparing the processing map under different strain conditions, the area of the plastic instability zone increases gradually with the increase of strain. When the deformation temperature is constant, the lower the strain rate, the less possibility of plastic instability. Under the same strain rate condition, the power dissipation efficiency increases gradually with the increase of deformation temperature, at the same times, the driving force of microstructure change increases gradually. The comprehensive analysis shows that the optimum deformation temperature and strain rate of 34CrNi3Mo steel should be higher than 1273 K and less than 10 s −1 . Comparing the processing map under different strain conditions, the area of the plastic instability zone increases gradually with the increase of strain. When the deformation temperature is constant, the lower the strain rate, the less possibility of plastic instability. Under the same strain rate condition, the power dissipation efficiency increases gradually with the increase of deformation temperature, at the same times, the driving force of microstructure change increases gradually. The comprehensive analysis shows that the optimum deformation temperature and strain rate of 34CrNi3Mo steel should be higher than 1273 K and less than 10 s −1 .

Conclusions
(1) The flow stress of 34CrNi3Mo steel in the thermal deformation process is the result of the interaction of work hardening, dynamic recovery, and recrystallization. Under the same deformation rate, with the increase of deformation temperature, the softening effect of dynamic recrystallization and dynamic recovery increases gradually, and the flow stress decreases gradually. However, under the same deformation temperature, with the increase of strain rate, the influence of work hardening on 34CrNi3Mo steel gradually dominates, and the flow stress increases gradually.
(2) The modified Arrhenius model can accurately reflect the influence of temperature, strain rate, and strain on flow stress under thermal deformation conditions, and it can also accurately predict the flow stress of 34CrNi3Mo steel under thermal deformation process.
(3) The evolution of the activation energy indicates that the deformation process of 34CrNi3Mo steel is a thermal activation process, which is affected by temperature, strain rate, and strain. With the increase of deformation temperature and strain rate, the value of deformation activation energy decreases gradually.
(4) According to the dynamic material model, the processing map of 34CrNi3Mo steel was established. At the same deformation temperature, the lower the strain rate, the less prone to plastic instability of 34CrNi3Mo steel. At the same strain rate, with the increase of deformation temperature, the efficiency of power dissipation gradually increases during plastic deformation, and the driving force for microstructure change of 34CrNi3Mo steel gradually increases. The optimal deformation temperature of 34CrNi3Mo steel should be higher than 1273 K and the strain rate should be less than 10 s −1 .

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author, [Xiang-Dong Jia], upon reasonable request.

Conflicts of Interest:
No conflict of interest exists in the submission of this manuscript, and the manuscript is approved by all authors for publication. I would like to declare, on behalf of my co-authors, that the work submitted was original research, not being published previously and not under consideration for publication elsewhere, in whole or in part. All the authors listed have approved the manuscript.