Hot Deformation Behavior and Constitutive Modeling of H13-Mod Steel

The H13-mod steel optimized by composition and heat treatment has reached the performance index of the shield machine hob. The hot deformation behavior of the H13-mod steel was investigated by compression tests in the temperature range from 900 to 1150 °C and the strain rate range from 0.01 to 10 s−1. The true stress-strain curve showed that the rising stress at the beginning of deformation was mainly caused by work hardening. After the peak stress was attained, the curve drop and the flow softening phenomenon became more obvious at low strain rates. The flow behavior of the H13-mod steel was predicted by a strain-compensated Arrhenius-type constitutive equation. The relationship between the material constant in the Arrhenius-type constitutive equation and the true strain was established by a sixth-order polynomial. The correlation coefficient between the experimental value and the predicted value reached 0.987, which indicated that the constitutive equation can accurately estimate the flow stress during the deformation process. A good linear correlation was achieved between the peak stress (strain), critical stress (strain) and the Zener‒Hollomon parameters. The processing maps of the H13-mod steel under different strains were established. The instability region was mainly concentrated in the high-strain-rate region; however, the microstructure did not show any evidence of instability at high temperatures and high strain rates. Combined with the microstructure and electron backscattered diffraction (EBSD) test results under different deformations, the optimum hot working parameters were concluded to be 998–1026 °C and 0.01–0.02 s−1 and 1140–1150 °C and 0.01–0.057 s−1.


Introduction
Shield machines are widely used in major projects such as urban subways, railway and highway transportation, energy transportation, underground passages, etc. They have the advantages of high excavation speed, high construction quality, low labor intensity, high cost effectiveness, and environmental protection [1]. The harsh working environment of shield machines necessitates that the tool have high hardness, high wear resistance, good impact toughness, certain anti-tempering performance, and good processing performance. At present, the international hob materials are mainly based on American AISI 4340 steel and German X50CrMoV51 steel [2]. The former exhibits higher toughness [3], whereas the latter exhibits higher hardness [1]. Because the domestic R&D of hobs started relatively late, most hob materials remain in the development stage; these materials include 9Cr2Mo [4], 5Cr5MoSiV1 [5], and 5Cr5MoV [6], which generally have shortcomings such as a poor combination of strength and toughness and short service life. Many international hob materials have been developed from AISI H13, including the aforementioned X50CrMoV51 steel. However, H13 steel itself has a high content of alloying elements, which tend to segregate during solidification, adversely affecting product performance. Numerous scholars have improved the  Figure 1 shows the true stress-true strain curve of the H13-mod steel at different strain rates. In the early stage of deformation, due to the work-hardening effect, the dislocation density increases continuously and the flow stress increases sharply with increasing strain. As the strain increases, the dislocation density and the dislocation rearrangement offset part of the increase in dislocation density due to work hardening and the flow stress slowly reaches a peak value. The peak stress gradually decreases as the strain rate decreases and the deformation temperature increases. The typical single peak of the curve illustrates that DRX is the primary softening mechanism [30]. During the subsequent deformation process, the flow stress gradually decreases and a substantial softening phenomenon occurs until the flow stress reaches a steady state. The steady-state stress reflects the strain hardening and strain softening caused by the formation of new grains and related grain-boundary migration [24]. Notably, the flow softening phenomenon is quite different in different deformation conditions. The occurrence of this phenomenon is related to deformation heating and microstructure changes such as DRV or DRX. To evaluate the flow softening behavior more intuitively, we use the extent of flow softening (∆σ/σ P , ∆σ = σ p − σ 0.6 ) according to the previous practice, where σ p is the peak stress and σ 0.6 is the flow stress at a strain of 0.6 [23]. The extent of flow softening under different deformation conditions is shown in Figure 2, which shows that at low strain rates, the extent of softening of the flow is greater because the low strain rate gives a sufficient incubation period for recrystallization; recrystallization promotes the softening process, and this phenomenon is more pronounced at high temperatures. As the strain rate increases, the extent of flow softening decreases. This phenomenon is attributed to higher work hardening speeds at high strain rates slowing the softening caused by DRX, thereby causing the curve to be stable [20]. This effect is also evident in the stress-strain curve. The curve in Figure 1a shows a dramatic drop after reaching the peak. In Figure 1d, especially at 1100 and 1150 • C, the stress shows only a slight decrease after the strain of 0.4 and the difference between the peak stress and the steady-state stress is small after the true strain of 0.8. At the same time, we observe that the flow stress tends to reach a steady state at high temperatures and low strain rates, which may be a consequence of higher temperatures and lower strain rates promoting the volume fraction of DRX [31]. These results are consistent with the results of flow softening. reaching the peak. In Figure 1d, especially at 1100 and 1150 °C, the stress shows only a slight decrease after the strain of 0.4 and the difference between the peak stress and the steady-state stress is small after the true strain of 0.8. At the same time, we observe that the flow stress tends to reach a steady state at high temperatures and low strain rates, which may be a consequence of higher temperatures and lower strain rates promoting the volume fraction of DRX [31]. These results are consistent with the results of flow softening.   Although the true stress-strain curve reflects the intrinsic relationship between flow stress and thermodynamics, we cannot elucidate the deformation mechanism of the material solely by the true stress-strain curve. For example, flow softening is often associated with microstructure mechanisms such as DRX, flow instability and cracking [16]. Therefore, we investigated the kinetics, processing maps and microstructures in detail, as discussed in the following subsections. reaching the peak. In Figure 1d, especially at 1100 and 1150 °C, the stress shows only a slight decrease after the strain of 0.4 and the difference between the peak stress and the steady-state stress is small after the true strain of 0.8. At the same time, we observe that the flow stress tends to reach a steady state at high temperatures and low strain rates, which may be a consequence of higher temperatures and lower strain rates promoting the volume fraction of DRX [31]. These results are consistent with the results of flow softening.   Although the true stress-strain curve reflects the intrinsic relationship between flow stress and thermodynamics, we cannot elucidate the deformation mechanism of the material solely by the true stress-strain curve. For example, flow softening is often associated with microstructure mechanisms such as DRX, flow instability and cracking [16]. Therefore, we investigated the kinetics, processing maps and microstructures in detail, as discussed in the following subsections. Although the true stress-strain curve reflects the intrinsic relationship between flow stress and thermodynamics, we cannot elucidate the deformation mechanism of the material solely by the true stress-strain curve. For example, flow softening is often associated with microstructure mechanisms such as DRX, flow instability and cracking [16]. Therefore, we investigated the kinetics, processing maps and microstructures in detail, as discussed in the following subsections.

Kinetic Analysis
DRX begins when the strain reaches a critical value; thus, the critical strain strongly affects DRX. In recent years, numerous methods for determining the critical strain have been reported [21,32]. The most widely used method is to fit the true stress-strain curve and then derive the stress to obtain the work-hardening rate curve. The critical stress is obtained from the inflection point of the curve. Figure 3 is a graph showing the work-hardening rate of the material at 900 • C and at a strain rate of 0.01 s −1 . From Figure 3, we can obtain the critical stress σ c , saturation stress σ s , stable stress σ ss , peak stress σ p under this deformation condition. The values of the various parameters under different deformation conditions are shown in Table 2. DRX begins when the strain reaches a critical value; thus, the critical strain strongly affects DRX. In recent years, numerous methods for determining the critical strain have been reported [21,32]. The most widely used method is to fit the true stress-strain curve and then derive the stress to obtain the work-hardening rate curve. The critical stress is obtained from the inflection point of the curve. Figure  3 is a graph showing the work-hardening rate of the material at 900 °C and at a strain rate of 0.01 s −1 . From Figure 3, we can obtain the critical stress σ , saturation stress , stable stress , peak stress under this deformation condition. The values of the various parameters under different deformation conditions are shown in Table 2.  The relationship between the mechanical and thermodynamic parameters of the metal involved in different deformation processes can be expressed using the flow stress equation. The most important data in the high-temperature deformation process are the flow stress; the maximum influence on the flow stress is the strain rate ε and the deformation temperature T. The relationship between and T can be expressed by the Garofalo formula. Based on the Garofalo formula, Sellars  The relationship between the mechanical and thermodynamic parameters of the metal involved in different deformation processes can be expressed using the flow stress equation. The most important data in the high-temperature deformation process are the flow stress; the maximum influence on the flow stress is the strain rate be expressed by the Garofalo formula. Based on the Garofalo formula, Sellars and Tgart [33] proposed a new constitutive equation that includes the thermal activation energy Q and the deformation temperature T, which is the Arrhenius equation: .
According to the stress level, this equation can be divided into the following three forms.
(1) When the stress is low (ασ < 0.8), it can be simplified into an exponential form: (2) When the stress is high (ασ < 1.2), it can be simplified into a power exponential form: (3) It can be simplified into a hyperbolic sinusoidal form over the entire stress range: In these formulas, .
ε is the strain rate, A is the structural factor of the material, α is the stress-level parameter (α = β/n), n is the stress index, σ is the flow stress, and Q is the thermal deformation activation energy. Whether the metal material is easily deformed can be estimated from the magnitude of the Q value. R is the molar gas constant, i.e., R = 8.31441 J·mol −1 ·K −1 , and T is the absolute temperature.
When Zener and Hollomon [34] studied the stress-strain relationship of steel, they found that the effect of deformation rate and temperature on flow stress can be expressed by the Z parameter: We linearly fitted the logarithm of Equations (2) and (3), where σ takes the peak stress as the calculation value. The relationship between ln σ p − ln . ε and σ p − ln . ε is shown in Figure 4a,b, respectively. Taking the average of the reciprocal of the slope as the values of n 1 and β, we found that n 1 = 9.894, β = 0.0594 and α = 0.0059 by calculation.
Furthermore, by taking the logarithm of the natural number for Equation (4), we find that ln [sin h(ασ p )] has a linear relationship with ln . ε and T −1 , as shown in Figure 4c,d. By linear regression analysis, we obtain n = 7.06 and Q = 563.62 kJ·mol −1 . At the same time, the relationship between lnZ and ln [sin h(ασ p )] is fitted and the intercept is the value of lnA, where A = 5.85 × 10 21 . The values of A, α, n and Q are calculated from the experimental data; thus, the flow stress equation of the material is The relationship between peak stress and Z parameter is as follows: Strain is generally assumed to have little effect on the thermal deformation behavior of metallic materials. However, a large number of studies have shown that the thermal deformation activation energy and the material constant are strongly affected by strain during the whole deformation process [23,[35][36][37]; we observe this effect from the true stress-strain curve. Therefore, strain should be introduced into the Arrhenius equation to improve the accuracy of the flow stress prediction. For this reason, based on the experimental data, we calculated the values of Q, n, α and lnA of H13-mod steel at a strain of 0.05 to 0.9 with an interval of 0.05. Table 3 shows the values of Q, n, α and lnA for H13-mod steel at different strains. These results show that the effect of strain on Q, n, α and lnA is substantial. These deformation constants and strains were subjected to the polynomial fitting shown in Equation (8), and the order was 1 to 9 [38,39]. Based on the correlation analysis, a six-step fitting was performed. The fitting coefficients are listed in Table 4. From Figure 5, the experimental data and the polynomial fitting results show a good correlation.
the effect of deformation rate and temperature on flow stress can be expressed by the Z parameter: We linearly fitted the logarithm of Equations (2) and (3), where takes the peak stress as the calculation value. The relationship between lnσ − lnε and σ − lnε is shown in Figure 4a,b, respectively. Taking the average of the reciprocal of the slope as the values of n1 and β, we found that n1 = 9.894, β = 0.0594 and α = 0.0059 by calculation.
Furthermore, by taking the logarithm of the natural number for Equation (4), we find that ln [sinℎ ( )] has a linear relationship with lnε and T −1 , as shown in Figure The relationship between peak stress and Z parameter is as follows: = 168.39ln   Strain is generally assumed to have little effect on the thermal deformation behavior of metallic materials. However, a large number of studies have shown that the thermal deformation activation energy and the material constant are strongly affected by strain during the whole deformation process [23,[35][36][37]; we observe this effect from the true stress-strain curve. Therefore, strain should be introduced into the Arrhenius equation to improve the accuracy of the flow stress prediction. For this reason, based on the experimental data, we calculated the values of Q, n, α and lnA of H13-mod steel at a strain of 0.05 to 0.9 with an interval of 0.05. Table 3 shows the values of Q, n, α and lnA for H13-mod steel at different strains. These results show that the effect of strain on Q, n, α and lnA is     After calculating the relationship between the material constant and strain, we can more accurately predict the stress value under a specific strain. According to Equation (3), the constitutive equation of the H13-mod steel under a specific strain can be written as A comparison between the experimental flow stress and flow stress predicted by the Arrhenius-type constitutive equation is shown in Figure 6. Through Figure 6, we observe that the experimental value has a high degree of coincidence with the predicted value, indicating that the proposed constitutive equation has higher prediction accuracy. To further illustrate the accuracy of the constitutive equation, we used the standard statistical parameters to calculate the correlation coefficient (R) and the mean absolute relative error (∆) [40,41] according to the equations: where σ E is the experimental flow stress, σ p is the predicted flow stress, σ E and σ p are the mean values of σ E and σ p , respectively. In addition, N is the total number of data employed in this work. The correlation coefficient (R) is commonly used to reflect the capability of the linear relationship between the experimental and predicted data. Figure 7 shows that a good correlation (R = 0.99) is obtained between the experimental and the predicted data. This correlation indicates that the proposed constitutive equations exhibit an excellent fit with the measured values. In the present work, the ∆ value is calculated to be approximately 6.3% for H13-mod steel deformed in the temperature range from 900 to 1150 • C and strain rate range from 0.01 to10 s −1 , which reflects the good prediction capabilities of the proposed constitutive equation. where is the experimental flow stress, is the predicted flow stress, and are the mean values of and , respectively. In addition, N is the total number of data employed in this work. The correlation coefficient (R) is commonly used to reflect the capability of the linear relationship between the experimental and predicted data. Figure 7 shows that a good correlation (R = 0.99) is obtained between the experimental and the predicted data. This correlation indicates that the proposed constitutive equations exhibit an excellent fit with the measured values. In the present work, the ∆ value is calculated to be approximately 6.3% for H13-mod steel deformed in the temperature range from 900 to 1150 °C and strain rate range from 0.01 to10 s −1 , which reflects the good prediction capabilities of the proposed constitutive equation.    Jorge and Balancin [42] suggested that the mechanical parameters can be described by simple Zener parameters. Moreover, this method is being cited in an increasing number of literature works. The Z parameters and these mechanical performance parameters can be described by the following power exponential equation: Jorge and Balancin [42] suggested that the mechanical parameters can be described by simple Zener parameters. Moreover, this method is being cited in an increasing number of literature works. The Z parameters and these mechanical performance parameters can be described by the following power exponential equation: The Z parameters and the changes in these parameters are shown in Figure 8a,b. This figure shows that the mechanical performance parameters have a good linear relationship with the Z parameters. By linear regression, we obtain the slopes of the different curves: At the same time, some studies have shown that there is a linear relationship between the critical strain ε c and the peak strain ε p , so we also plotted the relationship between ε c and ε p and Z parameters, as shown in Figure 8c. From the figure we can see that the two have a good linear relationship with the Z parameters:

Processing Map
Thus far, the most widely used model for metals in a thermal compression process is the dynamic material model (DMM) proposed by Prasad et al. [16,43,44]. The DMM links the large plastic deformation and the dissipation of the structure, possibly clarifying how the energy of the external action is dissipated through the plastic deformation of the workpiece. According to the principles of thermodynamics, the deformation process is regarded as a closed thermodynamic system. The entire deformation process is related to the dissipation of energy. According to the model, the total power

Processing Map
Thus far, the most widely used model for metals in a thermal compression process is the dynamic material model (DMM) proposed by Prasad et al. [16,43,44]. The DMM links the large plastic deformation and the dissipation of the structure, possibly clarifying how the energy of the external action is dissipated through the plastic deformation of the workpiece. According to the principles of thermodynamics, the deformation process is regarded as a closed thermodynamic system. The entire deformation process is related to the dissipation of energy. According to the model, the total power of the input system is P, and the energy consumed by the thermal deformation and the energy consumed by the evolution of the internal organizational structure are the dissipation amount G and the dissipation coordination J, respectively. The relationship is The strain-rate sensitivity coefficient m can be used to distribute the relationship between the dissipation amount G and the dissipation coefficient J during thermal deformation: According to irreversible thermodynamics, the relationship between ln σ and ln . ε can be approximately illustrated as Combining Equations (14) and (15), we can rewrite m as When the temperature and strain are constant, the stress and the strain rate of the specimen are constitutively related: σ = K . ε m or . ε = K 1 σ 1/m (K and K 1 are constants). Therefore, the dissipation J can be defined as For steady-state rheology of viscoplastic solids, m ranges from 0 to 1. Equation (14) shows that the larger the value of m, the larger G of the metal material, and that the smaller the value of m, the larger P of the metal material. When m = 1, the metal material is in a particularly ideal dissipative condition. At this time, the dissipation coefficient J reaches the maximum value J max = 0.5σ Equation (18) shows that η and m, which is a function of the strain ε, the strain rate . ε and the deformation temperature T, are closely related. The variation of η with T and . ε constitutes a power dissipation diagram and different regions in the map correspond to different microstructures. The value of η can reflect the thermal processing properties of the metal material to a certain extent. However, the power dissipation efficiency value of the instability region is sometimes also particularly high; thus, when the η value is large, the hot workability of the material is not necessarily good. Therefore, it is necessary for us to find the processing instability region.
According to the principle of entropy proposed by Ziegler [45], Prasad et al. applied the extreme principle of irreversible thermodynamics to continuum mechanics of large plastic flow and obtained the inequality criterion of plastic rheological instability: where . ε is the strain rate and ξ( . ε) is a dimensionless instability parameter. The variation of ξ( . ε) with temperature and strain rate constitutes an instability map and the flow instabilities are predicted to occur when ξ( . ε) becomes negative. We obtain the processing map of the material by superimposing the power dissipation map on the instability map. From the processing map, we obtain the strain rate and temperature range suitable for processing the material or some flow instability regions such as adiabatic shear bands or flow localization [23].
Based on the above theory, we have established processing maps with true strains of 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9, as shown in Figure 9. In the processing maps, the efficiency is represented by a contour line, the instability region is represented by colored shading and the instability level of the instability region is represented by a different color. Processing maps have a guiding role in the design of thermal processing. For example, processing in unstable areas should be avoided during processing [46]. It can be seen from the processing map that η and the unstable region of the flow vary with the strain. When the true strain is 0.4, the processing map exhibits two instability regions that occur at temperature ranges from 910 to 1020 • C and from 1030 • C to 1150 • C and strain-rate ranges from 0.03 to 10 s −1 and from 0.41 to 10 s −1 , respectively. In the processing stability region, two peaks of efficiency occur in this region: one is in the temperature range 998-1026 • C with a relatively high strain rate from 0.01 to 0.02 s −1 and the other is in the temperature range 1140-1150 • C with a strain rate from 0.01 to 0.057 s −1 (marked A and B in Figure 9a).
its instability levels reveals that a very low level of instability (>−0.2) in most areas of the instability region. This region has been reported to be a metastable region, and instability may not exist here [50]. This low-instability region (>−0.2) is more pronounced under low strain conditions. This material may have a larger processing range at low strains, which will be discussed in greater detail later in connection with the microstructure. Notably, the power dissipation efficiency values change substantially as A and B increase with increasing strain. We therefore calculated the variation in the efficiency of power dissipation with true strain at deformation temperatures of 1000 and 1150 °C and a strain rate of 0.01 s −1 . As shown in Figure 10, as the true strain increases, the peak efficiency of power dissipation increases gradually because the DRX volume fraction increases and consumes a large amount of energy. With further increases of the true strain, the peak efficiency of power dissipation decreases at a true strain of 0.6 at 1000 °C and a true strain of 0.5 at 1150 °C. This behavior has been previously attributed to the coarsening of recrystallized grains [51]. To verify this result, we conducted thermal compression simulation experiments at 1000 °C and true strains of 0.3 and 0.6 and analyzed the microstructure, size and grain boundaries and the average grain size by EBSD. Figure 11 shows the microstructure of three different true strains. Figure 11a shows that at a true strain of 0.3, a significant shear band is caused by flow instability, which affects the mechanical properties of the product after forming. As the true strain increases, this phenomenon disappears and the microstructure tends to be uniform. Figure 12a-c shows the type of grain boundary under these three deformation conditions. We calculated the grain boundary distribution under different deformations, as shown in Figure 12d. From Figure 12d, when the true strain is 0.3, the proportion of low-angle grain boundaries (LAGBs) is larger and the proportion of high-angle grain boundaries (HAGBs) is smaller. This result indicates that the recrystallization volume fraction is small at a true strain of 0.3. When the true strain is 0.6, the proportion of LAGBs is substantially reduced and the proportion of HAGBs is significantly increased. Moreover, the number of LAGBs decreases and is approximately the same as the number of HAGBs. This phenomenon can be interpreted on the basis of the mechanism of strain-induced In general, a higher value of η indicates that the microstructure involved in the deformation process consumes more energy, as is the case for DRV, DRX or a phase transition, and high η contributes to thermal processing [47]. DRX has been reported occur when the value of η is between 0.3 and 0.5 [48,49]. However, the values of the A and B regions in the present study are just within this range, indicating that DRX occurred under the deformation conditions; we also observe this phenomenon from the stress-strain curve. With increasing strain, the instability region gradually moves to the high strain rate; however, the positions of the A and B regions remain basically unchanged. At the same time, the unstable area gradually merges from two to one. Observation of its instability levels reveals that a very low level of instability (>−0.2) in most areas of the instability region. This region has been reported to be a metastable region, and instability may not exist here [50]. This low-instability region (>−0.2) is more pronounced under low strain conditions. This material may have a larger processing range at low strains, which will be discussed in greater detail later in connection with the microstructure.
Notably, the power dissipation efficiency values change substantially as A and B increase with increasing strain. We therefore calculated the variation in the efficiency of power dissipation with true strain at deformation temperatures of 1000 and 1150 • C and a strain rate of 0.01 s −1 . As shown in Figure 10, as the true strain increases, the peak efficiency of power dissipation increases gradually because the DRX volume fraction increases and consumes a large amount of energy. With further increases of the true strain, the peak efficiency of power dissipation decreases at a true strain of 0.6 at 1000 • C and a true strain of 0.5 at 1150 • C. This behavior has been previously attributed to the coarsening of recrystallized grains [51]. To verify this result, we conducted thermal compression simulation experiments at 1000 • C and true strains of 0.3 and 0.6 and analyzed the microstructure, size and grain boundaries and the average grain size by EBSD. Figure 11 shows the microstructure of three different true strains. Figure 11a shows that at a true strain of 0.3, a significant shear band is caused by flow instability, which affects the mechanical properties of the product after forming. As the true strain increases, this phenomenon disappears and the microstructure tends to be uniform. Figure 12a-c shows the type of grain boundary under these three deformation conditions. We calculated the grain boundary distribution under different deformations, as shown in Figure 12d. From Figure 12d, when the true strain is 0.3, the proportion of low-angle grain boundaries (LAGBs) is larger and the proportion of high-angle grain boundaries (HAGBs) is smaller. This result indicates that the recrystallization volume fraction is small at a true strain of 0.3. When the true strain is 0.6, the proportion of LAGBs is substantially reduced and the proportion of HAGBs is significantly increased. Moreover, the number of LAGBs decreases and is approximately the same as the number of HAGBs. This phenomenon can be interpreted on the basis of the mechanism of strain-induced boundary migration, in which LAGBs migrate to random HAGBs and are counteracted or swallowed by HAGBs during the process of DRX [52]. Figure 12 shows that with increasing true strain, the recrystallization volume fraction increases obviously; however, the size and grain-boundary distribution of true strain at 0.9 are not substantially different from those at 0.6. At the same time, the average grain size under the strains of 0.3, 0.6, and 0.9 was 2.044 µm, 1.73 µm, and 1.84 µm, respectively; thus, the grain size appears to first decrease and then increase. This behavior shows that the DRX is almost complete when the true strain is 0.6 and that the recrystallized grains grow when the true strain is 0.9. Through the aforementioned analysis, we can further explain that the growth of the grain at the true strain of 0.9 leads to a decrease in the peak power dissipation efficiency. Because grain-coarsening reduces the mechanical properties of the material, the true strain should not exceed 0.6 at 1000 and 1150 • C and a strain rate of 0.01 s −1 .

Microstructure
The DRX of steel depends on the strain, temperature and strain rate [53,54]. To study the deformation mechanism under different deformation conditions, we observed the microstructure of the material in the flow stability region and the flow instability region. Figure 13 shows the microstructure of the material in the flow stability region. Figure 13a shows that complete DRX occurred at 1050 °C and a strain rate of 0.01 s −1 and that the η in this region is high. The results

Microstructure
The DRX of steel depends on the strain, temperature and strain rate [53,54]. To study the deformation mechanism under different deformation conditions, we observed the microstructure of the material in the flow stability region and the flow instability region. Figure 13 shows the microstructure of the material in the flow stability region. Figure 13a shows that complete DRX occurred at 1050 • C and a strain rate of 0.01 s −1 and that the η in this region is high. The results obtained from the processing map are consistent with the microstructure results. DRX is prone to occur at lower strain rates due to the lower critical strain (i.e., critical dislocation density) required. The nucleation conditions of DRX expressed by dislocation density (ρ m ) and strain rate ( . ε) can be expressed by the formula [55]: where γ b is the boundary energy, M is the boundary mobility, L is the average slip distance of the dislocations, G is the shear modulus, b is the Burgers vector and K is a constant. According to Equation (20), at a given temperature, DRX nucleation requires a lower ρ m (i.e., critical strain) at lower strain-rate conditions. The critical strain of DRX has been reported to increase with increasing strain rate for various materials [56,57]. We observed a similar phenomenon in this work by comparing the critical strain at different strain rates in Table 2. The lower critical strain promotes nucleation of recrystallization, whereas the lower strain rate provides sufficient time for grain-boundary migration, which causes DRX grains to grow further, resulting in accelerated DRX kinetics. As the temperature increases, the recrystallized grains grow, as shown in Figure 13b,c. For the strain rate of 0.1 s −1 , the relatively fast deformation speed reduces the time of grain-boundary migration. At the same time, the energy storage during deformation is substantially lower than that at a high strain rate. The recrystallization process becomes slower, and the structure shown in Figure 13d appears. However, a uniform equiaxed grain appears at a strain rate of 1 s −1 , as shown in Figure 13e,f. We further analyzed this situation using a strain rate of 10 s −1 , as discussed later in this work. Figure 14 shows the microstructure of the flow instability region under different deformation conditions. As evident from Figure 14a, the lower temperature allows the deformed grains to coexist with the recrystallized grains and the value of η is low (17%) under this condition. Other authors have suggested that this microstrip structure is an unstable phenomenon at high strain rates [58].
This phenomenon is improved at 1000 • C. Figure 14c-e shows that the equiaxed crystal is uniform at a strain rate of 10 s −1 and that the recrystallized grain size increases with increasing deformation temperature. However, this region is marked as an instability region in the processing map. Under this deformation condition, the η value is 28%. This phenomenon has been previously attributed to a high density of dislocations promoting the nucleation of DRX grains [59]. At the same time, the increase of dislocation density and stored energy in a deforming material with increasing strain rate has been reported to be due to a strong dislocation-dislocation interaction and weak recovery [60]. Therefore, more energy is stored in the deformed sample at high strain rates. At the same time, the adiabatic heating effect produced at high strain rates is more pronounced. Moreover, at high strain rates, the time for grain-boundary migration is very short. These reasons together promote the nucleation of DRX. This situation, which occurs at high strain rates, has also appeared in several other reports [15,24,52,61]. At the same time, some scholars have questioned the instability parameters. They have speculated that if a large number of DRX grains are produced at high strain rates, the instability parameters will be difficult to predict under the deformation conditions [61]. Chiba et al. speculated that the ξ( . ε) value ranges from −1 to 0, that the processing region may be a metastable processing region, and that some instability phenomena may not occur in this region [50]. Meanwhile, we also marked the range of different ξ( . ε) values from the processing map. At high strain rates, most of the values correspond to metastable processing regions. We consider it a "safe haven" for processing because it helps to maximize production efficiency in industrial production. However, a uniform equiaxed grain appears at a strain rate of 1 s −1 , as shown in Figure 13e,f. We further analyzed this situation using a strain rate of 10 s −1 , as discussed later in this work.  Figure 14 shows the microstructure of the flow instability region under different deformation conditions. As evident from Figure 14a, the lower temperature allows the deformed grains to coexist with the recrystallized grains and the value of η is low (17%) under this condition. Other authors have suggested that this microstrip structure is an unstable phenomenon at high strain rates [58].
This phenomenon is improved at 1000 °C. Figure 14c-e shows that the equiaxed crystal is uniform at a strain rate of 10 s −1 and that the recrystallized grain size increases with increasing rates, the instability parameters will be difficult to predict under the deformation conditions [61]. Chiba et al. speculated that the ( ) value ranges from −1 to 0, that the processing region may be a metastable processing region, and that some instability phenomena may not occur in this region [50]. Meanwhile, we also marked the range of different ( ) values from the processing map. At high strain rates, most of the values correspond to metastable processing regions. We consider it a "safe haven" for processing because it helps to maximize production efficiency in industrial production.

Conclusions
The hot deformation behavior of H13-mod steel was analyzed by the flow curve, kinetics, processing map and microstructure under the deformation temperature range from 900 °C to 1150 °C and the strain rate range from 0.01 s −1 -10 s −1 .
The main conclusions are as follows: (1) In the initial stage of deformation, as the strain increases, the flow stress rapidly increases to a peak. After reaching the peak stress, the flow curve exhibits significant flow softening at low strain rates, mainly related to deformation heating and dynamic changes such as DRV and DRX. At high strain rates, flow softening is reduced by the higher work-hardening rate.

Conclusions
The hot deformation behavior of H13-mod steel was analyzed by the flow curve, kinetics, processing map and microstructure under the deformation temperature range from 900 • C to 1150 • C and the strain rate range from 0.01 s −1 -10 s −1 .
The main conclusions are as follows: (1) In the initial stage of deformation, as the strain increases, the flow stress rapidly increases to a peak. After reaching the peak stress, the flow curve exhibits significant flow softening at low strain rates, mainly related to deformation heating and dynamic changes such as DRV and DRX. At high strain rates, flow softening is reduced by the higher work-hardening rate. (2) The relationship between peak stress, strain rate, and deformation temperature can be expressed by . ε = 5.85 × 10 21 [sin h(0.0059lσ p )] exp − 563918.98 RT .
(3) A good linear correlation exists between peak stress (strain), critical stress (strain) stable stress, peak stress and Z parameters. (4) The relationship between material constants (α, Q, n, and lnA) and true strain in the Arrhenius-type constitutive equation is established by a sixth-order polynomial. The predicted flow stress values of the constitutive equations are in good agreement with the experimental values, indicating that the constitutive equation can accurately estimate the flow stress of H13-mod steel. (5) The processing maps under different strains were established. The instability region mostly corresponded to high-strain-rate conditions; however, the microstructure did not show any signs of instability at high temperatures and high strain rates. By studying the thermal compression test under different deformation amounts, we determined that the processing temperature is 998-1026 • C and 1140-1150 • C and the strain rate is 0.01-0.02 s −1 and 0.01-0.057 s −1 . Moreover, the true strain does not exceed 0.6.