Comparison of Modiﬁed Johnson–Cook Model and Strain–Compensated Arrhenius Constitutive Model for 5CrNiMoV Steel during Compression around Austenitic Temperature

: Isothermal compression behaviors of 5CrNiMoV steel were studied at temperatures of 870, 800, 750, and 700 ◦ C, with strain rates of 0.001, 0.005, 0.01, 0.05, and 0.1 s − 1 , the compression temperatures 870 and 800 ◦ C are above Ac3, as well as 750 and 700 ◦ C below Ac3 temperature. The Modiﬁed Johnson–Cook (MJC) model and the Strain–Compensated Arrhenius (SCA) model were employed to demonstrate the relationship between the ﬂow stress and the compression parameters. The correlation coefﬁcient ( R ) and average absolute relative error ( AARE ) between the calculational and experimental ﬂow stress were used to evaluate the accuracy of the two models. The results show that the effect of dynamic softening on ﬂow stress is much more signiﬁcant at higher temperatures and lower strain rates, while this effect is not obvious when the strain rate exceeds 0.005 s − 1 with the temperature below Ac3. The MJC model has a good accuracy close to the reference conditions (0.001 s − 1 and 700 ◦ C), and it is suitable to predict the plastic behavior when the ﬂow stress is lower than 200 Mpa. The unbiased AARE values were 6.82 and 5.71 for MJC model and SCA model, respectively, which implied the SCA model has a higher accuracy than the MJC model. The SCA model was believed to be capable of being used to illustrate the thermomechanical behavior of 5CrNiMoV tool steel in a wide range of plastic deformation conditions.


Introduction
Hot working tool steels are mainly used in the production of various dies that are employed at elevated temperature, such as hammer forging die, hot extrusion die, machine forging die, and casting die. Because its working conditions are extremely harsh and complex, high strength, toughness, and good wear resistance are required [1][2][3]. To satisfy these requirements, searching for new alloys with additional alloying elements and using a reasonable hot working process have always attracted significant attention [4][5][6].
The stress-strain curve of metallic materials plays an important role, not only for designing deformation parameters but also for performing nonlinear finite element numerical simulation to solve an elastoplastic deformation problem [7,8]. Generally, there are two methods to obtain a stress-strain curve at a predefined temperature and strain rate. Conventional methods including tensile or compression tests, and then converting from engineering stress-strain data to true stress-strain data [7,9]. The other uses finite element simulations or calculational/artificial neural network methods to predict the stress-strain curves [10][11][12][13]. For the sake of obtaining more reliable results, experimental measurements, in conjunction with finite element analysis, are usually performed and verify each other [10,14]. The methods are mainly dependent on the material characteristics

Materials and Experiments
The chemical composition of 5CrNiMoV steel used in this study is given in Table 1. The microstructure of the as-received forged steel mainly consists of ferrite and cementite. Thermal expansion tests were carried out using a DIL 805A dilatometer machine (TA Instruments, Hüllhorst, Germany), and cylindrical samples with a diameter of 4 mm and length of 10 mm were used. The single-pass isothermal compression tests were conducted using a thermomechanical test machine DIL 805D (TA Instruments, Hüllhorst, Germany), and the specimens have a diameter of 5 mm and length of 10 mm. A type S (Pt-Pt/Rh10%) thermocouple was spot welded in the center of the specimen surface and used for a closedloop control of the specimen temperature. Cylindrical samples in thermal expansion and isothermal compression tests were heated by induction and cooled using inert gas. The cooling gas was ejected from the holes of the induction coil and directly meet the surface of the sample, meanwhile the temperature of the sample was controlled by adjusting the heating power and gas intensity automatically. The isothermal compression procedure and experimental setup are shown in Figure 1. All specimens were heated to austenite homogenization temperature (870 • C) with a heating rate of 2 K/s, hold for 10 min to obtain a uniform microstructure, and then quenched to the deformation temperature with a cooling rate of 40 K/s. The specimens were individually compressed uniaxially up to a total true strain of 1.1 at 870, 800, 750, and 700 • C, with a strain rate of 0.001, 0.005, 0.01, 0.05, and 0.1 s −1 . The deformed specimens were quenched to room temperature after isothermal compression, with a cooling rate of 40 K/s and cooling media of helium. isothermal compression tests were heated by induction and cooled using inert gas. The cooling gas was ejected from the holes of the induction coil and directly meet the surface of the sample, meanwhile the temperature of the sample was controlled by adjusting the heating power and gas intensity automatically. The isothermal compression procedure and experimental setup are shown in Figure  1. All specimens were heated to austenite homogenization temperature (870 °C) with a heating rate of 2 K/s, hold for 10 min to obtain a uniform microstructure, and then quenched to the deformation temperature with a cooling rate of 40 K/s. The specimens were individually compressed uniaxially up to a total true strain of 1.1 at 870, 800, 750, and 700 °C, with a strain rate of 0.001, 0.005, 0.01, 0.05, and 0.1 s −1 . The deformed specimens were quenched to room temperature after isothermal compression, with a cooling rate of 40 K/s and cooling media of helium. The isothermal compression experimental data including time, temperature, strain, true strain, force, true stress, change of length, induction heating power, and gas cooling intensity, etc., were recorded simultaneously alongside the test. It should be noted that the deformation stamps that DIL 805D used are made of Si3N4, and the maximum force is 20 KN for safety reasons [38]. The load provided by DIL 805D is relatively limited compared to the commonly used thermomechanical test machines Gleeble 3500/3800, whose maximum compressive loads are 100 KN and 200 KN, respectively [39]. That is why when the compression test performed at lower deformation temperature and higher strain rate, such as 700 °C and 0.1 s −1 in the test matrix, only a true strain of 0.8 was reached.
The compressed samples were sectioned, grinded, polished, and etched using the standard metallographic method. TESCAN VEGA3 scanning electron microscopic (SEM) (Tescan, Brno-Kohoutovice, Czech Republic) with secondary electron (SE) and back-scattering electron (BSE) mode was used for microstructure observation.

Results and Discussion
The dilatometer test results show that the austenitic transformation start temperature (Ac1) is 729 °C and the end temperature (Ac3) is 758 °C, when the heating rate is 0.05 K/s. The detailed temperature-dilatometer curves corresponding to different heating rates ranging from 0.05 to 50 K/s are not shown here. For specimens without external stress, temperature is the main driving force of phase transformation, the heating rate is a major factor affecting the phase transformation kinetics. Different heating rates and chemical The isothermal compression experimental data including time, temperature, strain, true strain, force, true stress, change of length, induction heating power, and gas cooling intensity, etc., were recorded simultaneously alongside the test. It should be noted that the deformation stamps that DIL 805D used are made of Si3N4, and the maximum force is 20 KN for safety reasons [38]. The load provided by DIL 805D is relatively limited compared to the commonly used thermomechanical test machines Gleeble 3500/3800, whose maximum compressive loads are 100 KN and 200 KN, respectively [39]. That is why when the compression test performed at lower deformation temperature and higher strain rate, such as 700 • C and 0.1 s −1 in the test matrix, only a true strain of 0.8 was reached.
The compressed samples were sectioned, grinded, polished, and etched using the standard metallographic method. TESCAN VEGA3 scanning electron microscopic (SEM) (Tescan, Brno-Kohoutovice, Czech Republic) with secondary electron (SE) and back-scattering electron (BSE) mode was used for microstructure observation.

Results and Discussion
The dilatometer test results show that the austenitic transformation start temperature (Ac1) is 729 • C and the end temperature (Ac3) is 758 • C, when the heating rate is 0.05 K/s. The detailed temperature-dilatometer curves corresponding to different heating rates ranging from 0.05 to 50 K/s are not shown here. For specimens without external stress, temperature is the main driving force of phase transformation, the heating rate is a major factor affecting the phase transformation kinetics. Different heating rates and chemical compositions result in different Ac1 and Ac3. As is shown by the findings reported by Huang [1], where the transformation temperatures of Ac1 and Ac3 for 5CrNiMoV steel were 730 • C and 765 • C, respectively. Thus, the isothermal compression temperatures 870 • C and 800 • C are above the Ac3, as well as 750 • C and 700 • C below the Ac3. The microstructure of steel alloys is mainly austenite and phase transformation does not happen when the holding temperature above Ac3. Alloy element diffusion, grain growths, and dislocation density decreases may happen, but austenite is stable and the lattice type of the metal matrix does not change in this stage from the point of view of thermodynamics [15]. When the holding temperature is below Ac3, the microstructure is supercooled austenite which is not stable anymore and is prone to decomposition and transforming into pearlite, bainite, or martensite, resulting in the thermomechanical analysis being more complicated [5,15]. The isothermal phase transformation curves obtained using DIL 805A dilatometer confirmed that the undercooled-austenite was stable and the incubation time was enough for a compression test at different strain rates, which means the supercooled austenite did not decompose when the compressed temperature was 750 • C or 700 • C.

True Stress-Strain Curves of the 5CrNiMoV Steel
The true stress-strain curves of 5CrNiMoV steel under different compression conditions are shown in Figure 2. It can be seen that the deformation temperature and strain rate have a significant effect on the flow stress, similar results have also been reported by other scholars [12,28,33], despite being above or below Ac3 temperature. The flow stress increases with the increase in the strain rate or the decrease in the deformation temperature. The strain rate and true strain are the main influencing factors for work hardening, as well as temperature for dynamic softening which includes dynamic recovery and dynamic recrystallization. The true stress-strain curve can be roughly divided into three parts. The flow stress increases steeply until the true strain reached 0.1, and then gradually increased to the maximum point is the first part. Within this part, the work hardening effect dominates the flow stress in the early stage and then is counterbalanced by the dynamic softening. As the true strain increases continuously, the flow stress decreased from the maximum point gradually before the true strain reached nearly 0.8. The effect of strain-induced work hardening on the flow stress faded while the effect of dynamic softening enhanced at the same time. The strain larger than about 0.8 belongs to the third part and the flow stress increases slowly again, until the total strain reached the predefined value of 1.1.
The stress-strain curve of the isothermal compression test depends on two competing factors, work hardening and dynamic softening. The competing mechanisms can be found in many metallic materials, such as steel [19,22,28], aluminum [31,32] and titanium alloys [21,40], etc. That is why the true stress-strain curves of 5CrNiMoV steel are similar to the sine curve. It should be noted that the effect of dynamic softening on flow stress is much more significant at a higher temperature and lower strain rate, especially when the compression temperature is above Ac3, but the dynamic softening effect is not obvious when the strain rate is higher than 0.005 s −1 with the compression temperature below Ac3. Even at a constant compression temperature, the true strain that corresponds to the maximum flow stress decreases with the decreasing strain rate. In contrast to the work hardening effect, which is much more obvious at lower temperatures and higher strain rates. The softening mechanism of metallic alloys usually contains dynamic recovery and continuous/discontinuous recrystallization, which depends on the alloy grade and the microstructure within it [41]. Thermal vibration and thermal diffusivity of the alloying element atom are enhanced at an elevated temperature, in addition to the stored energy accumulated with the strain, the driving force is adequate for dislocation cross-slip and climb, which lead to dynamic softening accompanied by lower flow stress. The dynamic softening is influenced by the time required to complete the corresponding true strain. If the strain rate is large enough and higher than 100 s −1 , that usually can be seen in Split Hopkinson Pressure Bar tests, dynamic softening during deformation is not obvious any more [23,34,37]. With the increase in strain, dislocations generation, slip, multiplication, and pile-up can be found in the deformed alloy, which leads to work hardening and an increase in flow stress [18].  [23,34,37]. With the increase in strain, dislocations generation, slip, multiplication, and pile-up can be found in the deformed alloy, which leads to work hardening and an increase in flow stress [18].

Microstructures of the 5CrNiMoV Steel after Compression
The microstructure of the as-compressed samples that deformed with a constant strain rate 0.05 s −1 and different temperatures are shown in Figure 3. At temperatures below 750 °C, it is observed that the grains were severely deformed and presented as a continuous strip, few recrystallizations were found. With the increase in deformation temperature, fine recrystallized grains can be observed at the grain boundaries ( Figure 3c) and their size increases gradually ( Figure 3d). It should be mentioned that the cooling velocity after compression was 40 K/s, which is much higher than the critical cooling velocity for martensite transformation. If the effect of stress/strain-induced austenite decomposition can be ignored, the microstructure of the 5CrNiMoV steel should be martensite after cooling, no precipitation occurred in this process. It is implied that the dynamic softening phenomena were not significant when the compression temperature is lower than 750 °C, which is consistent with the results derived from the red stress-strain curves in Figure 2.
The SEM morphologies of the 5CrNiMoV steel compressed at 870 °C with different strain rates can be seen in Figures 4 and 3d. The corresponding true stress-strain curves are shown in Figure 2a. The grain size decreases with the increase in the strain rate because the grains have enough time to grow up when the specimen was compressed at a lower strain rate, which results in dynamic softening and lower flow stress. When the strain rate was 0.1 s −1 , fine recrystallized grains can be found (Figure 4d), recrystallization has been completed, and the flow stress raised again. The flow stress is not only influenced by the temperature, but also the grain size. That is why the stress-strain curve is dependent on the work hardening and dynamic softening.

Microstructures of the 5CrNiMoV Steel after Compression
The microstructure of the as-compressed samples that deformed with a constant strain rate 0.05 s −1 and different temperatures are shown in Figure 3. At temperatures below 750 • C, it is observed that the grains were severely deformed and presented as a continuous strip, few recrystallizations were found. With the increase in deformation temperature, fine recrystallized grains can be observed at the grain boundaries ( Figure 3c) and their size increases gradually ( Figure 3d). It should be mentioned that the cooling velocity after compression was 40 K/s, which is much higher than the critical cooling velocity for martensite transformation. If the effect of stress/strain-induced austenite decomposition can be ignored, the microstructure of the 5CrNiMoV steel should be martensite after cooling, no precipitation occurred in this process. It is implied that the dynamic softening phenomena were not significant when the compression temperature is lower than 750 • C, which is consistent with the results derived from the red stress-strain curves in Figure 2.
The SEM morphologies of the 5CrNiMoV steel compressed at 870 • C with different strain rates can be seen in Figures 3d and 4. The corresponding true stress-strain curves are shown in Figure 2a. The grain size decreases with the increase in the strain rate because the grains have enough time to grow up when the specimen was compressed at a lower strain rate, which results in dynamic softening and lower flow stress. When the strain rate was 0.1 s −1 , fine recrystallized grains can be found (Figure 4d), recrystallization has been completed, and the flow stress raised again. The flow stress is not only influenced by the temperature, but also the grain size. That is why the stress-strain curve is dependent on the work hardening and dynamic softening. Metals 2022, 12, x FOR PEER REVIEW 6 of 18

Modified Johnson-Cook Model
The JC constitutive model was first proposed by Johnson and Cook in 1983 to illustrate the relationship among stress, strain and strain rate, especially for metallic materials under deformation at a relatively high temperature and high strain rate [16,42]. It has been widely used for describing the plastic deformation behavior and the model can be expressed as follows:

Modified Johnson-Cook Model
The JC constitutive model was first proposed by Johnson and Cook in 1983 to illustrate the relationship among stress, strain and strain rate, especially for metallic materials under deformation at a relatively high temperature and high strain rate [16,42]. It has been widely used for describing the plastic deformation behavior and the model can be expressed as follows:

Modified Johnson-Cook Model
The JC constitutive model was first proposed by Johnson and Cook in 1983 to illustrate the relationship among stress, strain and strain rate, especially for metallic materials under deformation at a relatively high temperature and high strain rate [16,42]. It has been widely used for describing the plastic deformation behavior and the model can be expressed as follows: .
where σ is the equivalent flow stress, ε is the strain, . ε * is the equivalent plastic strain, . ε is the strain rate, ε r is the reference strain rate, T * is the homologous temperature, T r is the reference temperature, and T m is the melting temperature. A, B, C, n, and m are material dependent constants, where A is the yield stress at reference temperature and reference strain rate, B is the coefficient of strain hardening, C is the coefficient of strain rate hardening, n is the strain hardening exponent and m is the thermal softening exponent, respectively. The JC constitutive model considered the main factors that affect the flow stress are isotropic hardening, strain rate hardening, and thermal softening, and their total effect can be simply calculated by multiplying these three factors [17], which can be seen in the right part of Equation (1). According to the previous research, the JC model does not consider the coupled effects of strain rate and temperature, which means that the calculational flow stress could not agree well with the experimental results [16,19]. Several MJC models [28,34,36] have been proposed in the past decades to improve the accuracy of the JC model, and the equation used in this study can be expressed as follows: where A 1 , B 1 , B 2 , B 3 , C 1 , λ 1 , and λ 2 are material dependent constants, the other parameters have the same meaning as the JC model. It can be seen in Equation (4) that the MJC model has a similar expression compared to JC model and the right part of the equal sign also can be divided into three terms. However, the second term is the same as that in the JC model which is related to strain rate hardening effect, the first term substitutes the monomial to a polynomial equation for better describing the isotropic hardening effect, the third term is much more complicated, and the coupled effects of strain rate and temperature was considered. The details for obtaining the material constants in the MJC model can be found elsewhere [25,29,36]. A quasi-static condition usually at lower temperature and strain rate was regarded as the reference condition [33], which means the sample compressed at 700 • C with a strain rate of 0.001 s −1 (Figure 2d) in the experimental range was selected. At the reference temperature and strain rate, the second and third terms are equal to 1 and then Equation (4) can be simplified as follows: The true stress-strain point data extracted from the 20 stress-strain curves were used to calculate the material constants, in which the true strain ranges from 0.05 to 1.05 with an interval of 0.05. That means 21 data were obtained that correspond to the reference temperature and strain rate. The parametric fitting expression of flow stress and strain were established by using non-linear numerical fitting analysis and the least square method, with the results of A 1 = 166.62, B 1 = 638.21, B 2 = −1150.07, B 3 = 581.86, and the square of correlation coefficient was 0.957, which can be seen in Figure 5.
When the compression temperature is equal to the reference temperature, 700 • C, Equation (4) can be simplified as Equation (6) and C 1 can be determined at each strain rate, ranging from 0.001 to 0.1 s −1 .
A few methods have been proposed for solving C 1 , Zhao et al. [17] suggested that the main factors that affecting C 1 were strain rate and temperature which are independent of each other, a sine fitting function that includes the above two parameters were obtained and found that the function agreed well with the experimental values. Other researchers [16,24,34] draw the fitting, in which the slope of the fitting line is the value of C 1 . The latter method was employed in this study and the slope in Figure 6 with the value of C 1 = 0.078. When the compression temperature is equal to the reference temperature, 700 °C, Equation (4) can be simplified as Equation (6) and can be determined at each strain rate, ranging from 0.001 to 0.1 s −1 .
( + + + ) = 1 + ln * A few methods have been proposed for solving , Zhao et al. [17] suggested that the main factors that affecting were strain rate and temperature which are independent of each other, a sine fitting function that includes the above two parameters were obtained and found that the function agreed well with the experimental values. Other researchers [16,24,34] draw the ( ) ~ ln * curve and carried out linear fitting, in which the slope of the fitting line is the value of . The latter method was employed in this study and the slope in Figure 6 with the value of = 0.078. Rearrange Equation (4) and take natural logarithm on each side and the equation was changed to the following form.
ln ( + + + )(1 + ln * ) = ( + ln * )( − ) (7) There were four deformation temperatures from 700-870 °C with five strain rates from 0.001-0.1s −1 in the test matrix, 20 true stress-strain curves were obtained and 21 data were extracted from each curve corresponding to the strain rate from 0.05 to 1.05 with an  When the compression temperature is equal to the reference temperature, 700 °C, Equation (4) can be simplified as Equation (6) and can be determined at each strain rate, ranging from 0.001 to 0.1 s −1 .
( + + + ) = 1 + ln * A few methods have been proposed for solving , Zhao et al. [17] suggested that the main factors that affecting were strain rate and temperature which are independent of each other, a sine fitting function that includes the above two parameters were obtained and found that the function agreed well with the experimental values. Other researchers [16,24,34] draw the ( ) ~ ln * curve and carried out linear fitting, in which the slope of the fitting line is the value of . The latter method was employed in this study and the slope in Figure 6 with the value of = 0.078. Rearrange Equation (4) and take natural logarithm on each side and the equation was changed to the following form. There were four deformation temperatures from 700-870 °C with five strain rates from 0.001-0.1s −1 in the test matrix, 20 true stress-strain curves were obtained and 21 data were extracted from each curve corresponding to the strain rate from 0.05 to 1.05 with an Figure 6. The relationship between σ/ A 1 + B 1 ε + B 2 ε 2 + B 3 ε 3 ~ln . ε/ε r , used for obtaining C 1 in the MJC model. (4) and take natural logarithm on each side and the equation was changed to the following form.

Rearrange Equation
There were four deformation temperatures from 700-870 • C with five strain rates from 0.001-0.1 s −1 in the test matrix, 20 true stress-strain curves were obtained and 21 data were extracted from each curve corresponding to the strain rate from 0.05 to 1.05 with an interval of 0.05. All these data were substituted in Equation (7)  interval of 0.05. All these data were substituted in Equation (7) and the value of + ln * at each strain rate can be calculated by linear fitting. All the results are shown in

Strain-Compensated Arrhenius Model
The flow stress of metallic material during deformation is not only related to the chemical composition and microstructure, but also the deformation temperature, strain rate, and true strain. The hyperbolic-sine Arrhenius equation and Zener-Hollomon parameter ( ) were widely used to predict the flow stress and illustrate the deformation behavior in the past few decades [26,27,36], it can be written as follows: The MJC model used to describe the relationship among flow stress σ, strain ε, strain rate . ε, and compression temperature T was established and expressed in the following form, where the reference strain rate ε r and reference temperature T r were 0.001 s −1 and 700 • C, respectively. σ = 166.62 + 638.21ε − 1150.07ε 2 + 581.86ε 3

Strain-Compensated Arrhenius Model
The flow stress of metallic material during deformation is not only related to the chemical composition and microstructure, but also the deformation temperature, strain rate, and true strain. The hyperbolic-sine Arrhenius equation and Zener-Hollomon parameter (Z) were widely used to predict the flow stress and illustrate the deformation behavior in the past few decades [26,27,36], it can be written as follows: where A, α, n, and Q are material dependent constants; . ε and σ have the same meaning with the MJC model; Q, R, and T denote the activation energy, the universe gas constant, and the absolute temperature, respectively.
The detailed procedures for solving the material constants also can be found elsewhere [26,27,31]. The calculation order of the constants or parameters related to the Arrhenius Model is as follows. n , which is the inverse of the slope of the fitting line ln . ε ∼ ln σ at each temperature. β, which is the inverse of the slope of the fitting line ln . ε ∼ σ at each temperature. α, α = β/n . n, which is the inverse of the slope of the fitting line ln[sinh(ασ)] ∼ ln . ε at each temperature. Q, which equals to the slope of the fitting line ln[sinh(ασ)] − 1/T at each strain rate multiply n and R. The calculational process of peak stress at different temperature and strain rates is shown in Figure 8.
where , , , and are material dependent constants; and have the same meaning with the MJC model; , , and denote the activation energy, the universe gas constant, and the absolute temperature, respectively.
The detailed procedures for solving the material constants also can be found elsewhere [26,27,31]. The calculation order of the constants or parameters related to the Arrhenius model is as follows.
, which is the inverse of the slope of the fitting line ~ at each temperature. , which is the inverse of the slope of the fitting line ~ at each temperature. , = ⁄ . , which is the inverse of the slope of the fitting line [ ℎ( )] ~ at each temperature. , which equals to the slope of the fitting line [ ℎ( )] − 1/ at each strain rate multiply and . The calculational process of peak stress at different temperature and strain rates is shown in Figure 8.  Based on the procedures mentioned above, the material constants of , , and in the SCA equation can be obtained at strain from 0.05 to 1.05 with an interval of 0.05, and the results are shown in Table 2 and red dots in Figure 9. It can be seen the strain has a significant effect on the material constants. As to improve the accuracy of the Arrhenius model, 6th order polynomial equations as shown in Equation 11 were used to demonstrate the relations between strain and the strain-dependent constants. The non-linear polynomial fitting curves and the square of the regression coefficients that correspond to , , , and also can be found in Figure 9. Each square of the regression coefficient is close to 1, which indicates the proposed 6th order polynomial equations have a higher prediction accuracy and the parameters in Equation (11) are shown in Table 3.  Table 2 and red dots in Figure 9. It can be seen the strain has a significant effect on the material constants. As to improve the accuracy of the Arrhenius Model, 6th order polynomial equations as shown in Equation (11) were used to demonstrate the relations between strain and the strain-dependent constants. The non-linear polynomial fitting curves and the square of the regression coefficients that correspond to α, n, Q, and lnA also can be found in Figure 9. Each square of the regression coefficient is close to 1, which indicates the proposed 6th order polynomial equations have a higher prediction accuracy and the parameters in Equation (11) are shown in Table 3.     Table 3. Polynomial fitting parameters of α, n, Q, and lnA in SCA constitutive model. The SCA constitutive equation related to Zener-Hollomon parameter can be rewritten in the following form, and the flow stress during compression can be predicted with the material constants α, n, Q, and A, which are a function of 6th order polynomial equations with the parameters in Table 3.

Evaluation of the Constitutive Models
Comparisons between the experimental and calculational flow stress using MJC and SCA models are shown in Figures 10 and 11, respectively. The MJC model presents a good accuracy close to the reference conditions (0.001 s −1 and 700 • C), the discrepancy between the predicted and measured stress values are negligible when the strain rate equal to 0.001 s −1 . As the compression temperature and strain rate increased from the reference conditions to a higher level, the discrepancy seems to expand. Moreover, when the true strain was larger than 0.6, the discrepancy was also gradually expanding with the increase in strain. The SCA model seems to have a relatively higher accuracy when the compression temperature is below Ac3 with a strain rate lower than 0.01 s −1 , and when it is higher than Ac3 with a strain rate larger than 0.01 s −1 . However, from the point of the shape of the predicted curves, the SCA model is much better than the MJC model. strain. The SCA model seems to have a relatively higher accuracy when the compression temperature is below Ac3 with a strain rate lower than 0.01 s −1 , and when it is higher than Ac3 with a strain rate larger than 0.01 s −1 . However, from the point of the shape of the predicted curves, the SCA model is much better than the MJC model. The deviation of the calculational flow stress from the experimental curves was investigated quantitatively by using two parameters, the correlation coefficient ( ) and the absolute average error ( ), and their equations can be expressed as follows [13,16,17]: where is the experimental flow stress and is the calculational flow stress predicted using the above two models, and are the average of experimental and calculational flow stresses, respectively.
is the number of data (equal to 21) used in this study for a single true stress-strain curve. The and values under different compression temperatures and strain rates are shown in Table 4. Generally, a higher value of and lower value of are usually considered as better parameters for prediction accuracy for the constitutive models. It is interesting that the values in the MJC model, which are varied in a wide range from 0.19 to 0.98 at different conditions, but have a relatively higher value of 0.97 when all data were considered. The correlation coefficient reflects the degree of linear correlation between the calculational and experimental flow stress, the greater the absolute value of indicates the stronger the correlation. When the changing trend of the flow stress, with the increase in strain, is similar for the calculational and experimental values, the absolute value of will increase correspondingly. The maximum value of could reach is 1, which means there exists a linear relationship between the calculational and experimental values. Each value in Table 4 that corresponds to a given compression condition can not represent the overall trend of flow stress data. The value of depends on the stress-strain data that is selected, that is why the value was much higher when all data were considered in the MJC model. Other researchers also mentioned that a higher value does not always mean better accuracy of the model [13,17], because compared to the experimental curves, the calculational flow stress may deviate upward or downward. The unbiased values can also be found in Table 4, and the values are 6.82% and 5.71% when all the stress-strain data were considered, which also indicates that the SCA model is more accurate. It should be noted that it is much more difficult to obtain the material dependent constants of the SCA constitutive model in terms of computational cost and complexity. The deviation of the calculational flow stress from the experimental curves was investigated quantitatively by using two parameters, the correlation coefficient (R) and the absolute average error (AARE), and their equations can be expressed as follows [13,16,17]: where E i is the experimental flow stress and P i is the calculational flow stress predicted using the above two models, E and P are the average of experimental and calculational flow stresses, respectively. N is the number of data (equal to 21) used in this study for a single true stress-strain curve. The R and AARE values under different compression temperatures and strain rates are shown in Table 4. Generally, a higher value of R and lower value of AARE are usually considered as better parameters for prediction accuracy for the constitutive models. It is interesting that the R values in the MJC model, which are varied in a wide range from 0.19 to 0.98 at different conditions, but have a relatively higher value of 0.97 when all data were considered. The correlation coefficient reflects the degree of linear correlation between the calculational and experimental flow stress, the greater the absolute value of R indicates the stronger the correlation. When the changing trend of the flow stress, with the increase in strain, is similar for the calculational and experimental values, the absolute value of R will increase correspondingly. The maximum value of R could reach is 1, which means there exists a linear relationship between the calculational and experimental values. Each R value in Table 4 that corresponds to a given compression condition can not represent the overall trend of flow stress data. The value of R depends on the stress-strain data that is selected, that is why the R value was much higher when all data were considered in the MJC model. Other researchers also mentioned that a higher R value does not always mean better accuracy of the model [13,17], because compared to the experimental curves, the calculational flow stress may deviate upward or downward. The unbiased AARE values can also be found in Table 4, and the values are 6.82% and 5.71% when all the stress-strain data were considered, which also indicates that the SCA model is more accurate. It should be noted that it is much more difficult to obtain the material dependent constants of the SCA constitutive model in terms of computational cost and complexity. The correlation between the experimental and calculational flow stress which is predicted using the MJC and SCA models is shown in Figure 12. The black line with a slope of 1 indicates the exactly suitable position between the experimental and calculational values. It can be seen that the MJC model has many data points that deviated from the black line, while the SCA model has fewer. Moreover, the degree of deviation is relatively small when the flow stress is lower than about 200 Mpa for the MJC model and 325 Mpa for the SCA model, which indicates that the two proposed models are suitable to predict the flow stress at higher temperatures and lower strain rates, but for higher flow stress usually obtained at lower temperatures or higher strain rates, the SCA models seem much better.  The correlation between the experimental and calculational flow stress which is predicted using the MJC and SCA models is shown in Figure 12. The black line with a slope of 1 indicates the exactly suitable position between the experimental and calculational values. It can be seen that the MJC model has many data points that deviated from the black line, while the SCA model has fewer. Moreover, the degree of deviation is relatively small when the flow stress is lower than about 200 Mpa for the MJC model and 325 Mpa for the SCA model, which indicates that the two proposed models are suitable to predict the flow stress at higher temperatures and lower strain rates, but for higher flow stress usually obtained at lower temperatures or higher strain rates, the SCA models seem much better.

Summary
Based on the above study, some conclusions could be drawn as follows: (1) The true stress-strain curve is much more like a sine curve when the compression temperature is above Ac3 or the strain rate is small enough, because the dynamic softening effect is more significant at a higher temperature and a lower strain rate. The softening effect is not obvious when the temperature is below Ac3 with a strain rate that exceeds 0.005 s −1 .
(2) The MJC model has a good accuracy close to the reference condition (0.001 s −1 and 700 • C). With the increase in compression temperature or increase in strain rate, the discrepancy between the MJC predicted flow stress and the measured flow stress expanded gradually, especially when the true strain larger than 0.6.
(3) The correlation coefficient R reflects the degree of linear correlation between the calculational and experimental flow stress under different conditions. The value of R depends on the stress-strain data that are selected, a higher R value does not imply a better prediction accuracy. The unbiased AARE values were 6.82 and 5.71 for MJC model and SCA model, respectively, which implied the SCA model has higher accuracy than the MJC model.
(4) The MJC and SCA models are suitable to predict the thermomechanical behavior of 5CrNiMoV steel when the flow stress lower than 200 Mpa. For a higher flow stress state, the SCA model is more preferred.