On the Prediction of the Flow Behavior of Metals and Alloys at a Wide Range of Temperatures and Strain Rates Using Johnson–Cook and Modified Johnson–Cook-Based Models: A Review

This paper reviews the flow behavior and mathematical modeling of various metals and alloys at a wide range of temperatures and strain rates. Furthermore, it discusses the effects of strain rate and temperature on flow behavior. Johnson–Cook is a strong phenomenological model that has been used extensively for predictions of the flow behaviors of metals and alloys. It has been implemented in finite element software packages to optimize strain, strain rate, and temperature as well as to simulate real behaviors in severe conditions. Thus, this work will discuss and critically review the well-proven Johnson–Cook and modified Johnson–Cook-based models. The latest model modifications, along with their strengths and limitations, are introduced and compared. The coupling effect between flow parameters is also presented and discussed. The various methods and techniques used for the determination of model constants are highlighted and discussed. Finally, future research directions for the mathematical modeling of flow behavior are provided.


Introduction
Hot deformation is one of the most well-known ways to enhance the mechanical properties of metals and alloys [1][2][3][4][5][6] via grain refinements [7][8][9][10], in which temperature is raised above recrystallization temperature during plastic deformation. The enhancement to the mechanical properties is controlled by work hardening, dynamic recovery (DRV), and dynamic recrystallization (DRX), which have a huge effect on the microstructure as well as the flow stress behavior of the alloys [11][12][13][14][15]. Accordingly, proper hot working parameters such as strain, strain rate, and temperature must be carefully chosen to achieve the desired mechanical properties [16][17][18][19][20]. The hot working parameters can be optimized by employing finite element simulations of the hot deformation process [21][22][23][24][25][26], in which the constitutive modeling of the flow stress behavior of alloys plays a significant role [27][28][29]. In addition to the importance of hot deformation, modeling the flow stress and plastic deformation of metals and alloys that are employed in applications under severe conditions such as very high strain rates (dynamic loadings) and different temperatures is also essential [30][31][32].
The constitutive models for the prediction of flow stress behavior at different strain rates and temperatures can be categorized as physical-based models [33][34][35][36], in which the physical aspects of material behavior are taken into consideration; phenomenological-based models [37][38][39][40][41], where empirical observations with mathematical functions establish the flow stress; and intelligence-based models [42][43][44][45][46], in which machine learning is used The flow behaviors of the B07 [99] and GH5188 [100] superalloys at elevated temperatures and a strain rate of 0.1 s −1 are shown in Figure 2a and Figure 2b, respectively. As can be seen, the flow stress is extremely sensitive to temperature and passes the three regions that are defined in Figure 1. The obtained stress-strain curves for the B07 and GH5188 superalloys under hot deformation show that both alloys are greatly affected by The flow behaviors of the B07 [99] and GH5188 [100] superalloys at elevated temperatures and a strain rate of 0.1 s −1 are shown in Figures 2a and 2b, respectively. As can be seen, the flow stress is extremely sensitive to temperature and passes the three regions that are defined in Figure 1. The obtained stress-strain curves for the B07 and GH5188 superalloys under hot deformation show that both alloys are greatly affected by hot working parameters, i.e., strain, strain rate, and temperature, in which stress decreases as the temperature increases and the strain rate decreases. Similarly, the same flow stress behavior under hot deformation has been reported for many different alloys, such as iron [101][102][103], titanium [104][105][106], aluminum [107][108][109][110], copper [111,112], nickel [113][114][115], and magnesium [116,117] alloys. The flow behaviors of the B07 [99] and GH5188 [100] superalloys at elevated temperatures and a strain rate of 0.1 s −1 are shown in Figure 2a and Figure 2b, respectively. As can be seen, the flow stress is extremely sensitive to temperature and passes the three regions that are defined in Figure 1. The obtained stress-strain curves for the B07 and GH5188 superalloys under hot deformation show that both alloys are greatly affected by hot working parameters, i.e., strain, strain rate, and temperature, in which stress decreases as the temperature increases and the strain rate decreases. Similarly, the same flow stress behavior under hot deformation has been reported for many different alloys, such as iron [101][102][103], titanium [104][105][106], aluminum [107][108][109][110], copper [111,112], nickel [113][114][115], and magnesium [116,117] alloys. As a result of DRX, dislocation is formed due to the release of stored elastic energy in the hardening region, which results in homogeneous sub-grains. As the strain increases, the misorientations between sub-grains increase, and the sub-grains are turned into fine grains [118][119][120][121][122]. Figure 3 shows the formation of new fine grains in AA1421 aluminum alloy due to DRX, in which aluminum alloy was processed using equal channel angular extrusion (ECAP) repeatedly to a strain of 12 at 673 K. As a result of DRX, dislocation is formed due to the release of stored elastic energy in the hardening region, which results in homogeneous sub-grains. As the strain increases, the misorientations between sub-grains increase, and the sub-grains are turned into fine grains [118][119][120][121][122]. Figure 3 shows the formation of new fine grains in AA1421 aluminum alloy due to DRX, in which aluminum alloy was processed using equal channel angular extrusion (ECAP) repeatedly to a strain of 12 at 673 K.  Figure 4 shows the peak stress vs. strain rate for Fe-26Mn-6.2Al-0.05C steel allo 1000 °C [124] (cf. Figure 4a) and Duplex cast steel alloy at 700 °C [125] (cf. Figure 4b). figure shows that, as the strain rate increases, the peak stress increases for both allo Similar findings can be found in [126][127][128][129][130]. As a matter of fact, this increase is attribu to the speedy rate of strain hardening, especially in the initial stage. Thus, the generat and multiplication of dislocations are accompanied by higher strain rates, in which h stresses are required due to dislocation interactions [131][132][133]. On the contrary, ther not enough time for DRV and DRX as restoration mechanisms [134,135].  Figure 4 shows the peak stress vs. strain rate for Fe-26Mn-6.2Al-0.05C steel alloy at 1000 • C [124] (cf. Figure 4a) and Duplex cast steel alloy at 700 • C [125] (cf. Figure 4b).

Effect of Strain Rate
The figure shows that, as the strain rate increases, the peak stress increases for both alloys. Similar findings can be found in [126][127][128][129][130]. As a matter of fact, this increase is attributed to the speedy rate of strain hardening, especially in the initial stage. Thus, the generation and multiplication of dislocations are accompanied by higher strain rates, in which high stresses are required due to dislocation interactions [131][132][133]. On the contrary, there is not enough time for DRV and DRX as restoration mechanisms [134,135].  Figure 4 shows the peak stress vs. strain rate for Fe-26Mn-6.2Al-0.05C steel alloy at 1000 °C [124] (cf. Figure 4a) and Duplex cast steel alloy at 700 °C [125] (cf. Figure 4b). The figure shows that, as the strain rate increases, the peak stress increases for both alloys. Similar findings can be found in [126][127][128][129][130]. As a matter of fact, this increase is attributed to the speedy rate of strain hardening, especially in the initial stage. Thus, the generation and multiplication of dislocations are accompanied by higher strain rates, in which high stresses are required due to dislocation interactions [131][132][133]. On the contrary, there is not enough time for DRV and DRX as restoration mechanisms [134,135]. Shi et al. [136] studied the effect of strain rate on the hot deformation properties of the GH690 superalloy. Figure 5 shows that, with a lower strain rate (0.001 s −1 ), fine grains due to full DRX are developed. It also shows the formation of a few column grains in addition to the equiaxed grains as the strain rate increases (0.1 s −1 and 5 s −1 ). The effect of strain rate on the hot deformation of other alloys has been reported with the same findings, such as magnesium [137], CoCrFeMnNi high-entropy [138], nickel [139,140], aluminum [141], iron [142], and titanium [143] alloys. Shi et al. [136] studied the effect of strain rate on the hot deformation properties of the GH690 superalloy. Figure 5 shows that, with a lower strain rate (0.001 s −1 ), fine grains due to full DRX are developed. It also shows the formation of a few column grains in addition to the equiaxed grains as the strain rate increases (0.1 s −1 and 5 s −1 ). The effect of strain rate on the hot deformation of other alloys has been reported with the same findings, such as magnesium [137], CoCrFeMnNi high-entropy [138], nickel [139,140], aluminum [141], iron [142], and titanium [143] alloys.  Figure 6 shows peak stress vs. temperature for Fe-26Mn-6.2Al-0.05C steel alloy [124] (Figure 6a) and Duplex cast steel alloy [125] (Figure 6b) at 1 s −1 and 0.1 s −1 , respectively. As can be seen, the peak stress decreases as the temperature increases. Similar findings can be found in [144][145][146][147][148]. At high temperatures and low strain rates, the flow stress decreases due to the slow rate of both DRV and DRX [149,150], in which the interaction and annihilation of the dislocations can be guided by adequate thermal energy in plenty of time [151,152].  Figure 6 shows peak stress vs. temperature for Fe-26Mn-6.2Al-0.05C steel alloy [124] ( Figure 6a) and Duplex cast steel alloy [125] (Figure 6b) at 1 s −1 and 0.1 s −1 , respectively. As can be seen, the peak stress decreases as the temperature increases. Similar findings can be found in [144][145][146][147][148]. At high temperatures and low strain rates, the flow stress decreases due to the slow rate of both DRV and DRX [149,150], in which the interaction and annihilation of the dislocations can be guided by adequate thermal energy in plenty of time [151,152]. Figure 6 shows peak stress vs. temperature for Fe-26Mn-6.2Al-0.05C steel alloy [124] ( Figure 6a) and Duplex cast steel alloy [125] (Figure 6b) at 1 s −1 and 0.1 s −1 , respectively. As can be seen, the peak stress decreases as the temperature increases. Similar findings can be found in [144][145][146][147][148]. At high temperatures and low strain rates, the flow stress decreases due to the slow rate of both DRV and DRX [149,150], in which the interaction and annihilation of the dislocations can be guided by adequate thermal energy in plenty of time [151,152]. Song et al. [153] studied the effect of hot deformation and DRX on a new Ni-Cr-Co nickel-based superalloy. Figure 7 shows that increasing the temperature leads to the formation of fine grains. Elongated grains can be seen in Figure 7a, passing with DRX at grain boundaries in Figure 7b and, finally, complete DRX and the growth of fine grains, as in Figure 7c. The effect of temperature on different alloys has been noticed with the same findings, such as nickel [154,155], Al-Cu-Mg [156], titanium [157] Ni-Co [158], and aluminum [159] alloys.  Figure 6 shows peak stress vs. temperature for Fe-26Mn-6.2Al-0.05C steel alloy [124] ( Figure 6a) and Duplex cast steel alloy [125] (Figure 6b) at 1 s −1 and 0.1 s −1 , respectively. As can be seen, the peak stress decreases as the temperature increases. Similar findings can be found in [144][145][146][147][148]. At high temperatures and low strain rates, the flow stress decreases due to the slow rate of both DRV and DRX [149,150], in which the interaction and annihilation of the dislocations can be guided by adequate thermal energy in plenty of time [151,152]. Song et al. [153] studied the effect of hot deformation and DRX on a new Ni-Cr-Co nickel-based superalloy. Figure 7 shows that increasing the temperature leads to the formation of fine grains. Elongated grains can be seen in Figure 7a, passing with DRX at grain boundaries in Figure 7b and, finally, complete DRX and the growth of fine grains, as in Figure 7c. The effect of temperature on different alloys has been noticed with the same findings, such as nickel [154,155], Al-Cu-Mg [156], titanium [157] Ni-Co [158], and aluminum [159] alloys.

Constitutive Models
In this section, the well-known JC and modified JC-based models are presented and carefully reviewed. In addition, the associated methods to determine the models' constants are also considered.

Original Johnson-Cook Model
Johnson and Cook [47] introduced their famous model to predict the flow behavior of materials at high temperatures and different strain rates. The model is very simple, and three independent parts were presented. The three parts are (i) strain hardening, (ii) strain rate, and (iii) thermal softening. Only five constants are introduced to constitute the JC model. The model can be expressed as σ = (A + Bε n )(1 + C ln ε · * )(1 − T * m ) (1) where σ is the flow stress, and ε is the plastic strain. The strain-hardening part is introduced by the term A + Bε n , where A is yield stress, and B and n are the strain-hardening strength and strain-hardening exponent. The strain rate term is represented by 1 + C ln ε · * , in which the strain rate constant, C, is interrelated with the strain rate. Finally, the softening term is represented by 1 − T * m , in that a softening constant, m, is associated with temperature. The ε · * represents a dimensionless value of the strain rate, in which the strain rate, ε · , is scaled to a reference strain rate, ε · • , so ε · * = ε · /ε · • . T * is a dimensionless temperature term introduced by (T − T r )/(T m − T r ), in which T, T r , and T m are the tested, reference, and melting temperatures.
To determine the JC constants, reference values for both the strain rate and temperature must be chosen at the beginning. In both reference values, Equation (1) reduces to Taking the logarithm after performing some rearrangements, Equation (2) can be linearly expressed as The value of constant B is obtained from the intercept, and n is the slope of the equation, which is obtained by plotting ln(σ − A) vs. ln(ε). At the reference temperature, After performing some rearrangements, Equation (4) can be linearly expressed as Constant C is the slope of the equation, in that it is obtained by plotting σ/(A + Bε n ) vs. ln(ε · * ). To accurately compute the strain rate sensitivity of coefficient C in the JC model, it is recommended that the ratio between every two following strain rates is not less than 0.1 [160], since the exponent, C, has to be calculated from the curves of the logarithmic graph of the true stress vs. true plastic strain for different strain rates.
By taking the logarithm after performing some rearrangements, at reference strain rate, Equation (1) can be written as Constant m is the slope of the equation, which is obtained by plotting ln[1 − σ/(A + Bε n )] vs. T * .
Li et al. [80] determined the JC constants for T24 steel alloy at elevated temperatures and different strain rates. Figure 8 shows the determination of JC constants for the T24 steel alloy, in which A is the yield stress, which has been found to have a value of 100 MPa, and ln B = 4.3567, which provides a B of 78 MPa, while n has the value of the slope, which is found to be 0.2742 (cf. Figure 8a). Constant C is determined as 0.08 (cf. Figure 8b), while constant m has a value of 0.5847 (cf. Figure 8c). The ability of the JC model to accurately predict the flow behavior of T24 steel alloy is assessed using the correlation coefficient (R), which provides a value of 0.962 (cf. Figure 8d).
One of the most powerful advantages of the JC model is that it is implemented in finite element software packages, in which it can be used to mimic real hot working processes and applications in severe conditions using simulation analysis. Mosleh et al. [161] implemented the JC model and two other models with finite element simulations (FES) to optimize the superplastic-forming of Ti-6%Al-4%V titanium alloy. Using FES, the optimized hot working parameters, strain rate, and temperate were determined as 0.002 s −1 and 875 • C. The FES was also helpful in checking the behavior of the superplastic-forming of the tested alloy (cf. Figure 9). Many other simulation analyses have been achieved based on the JC model [162][163][164][165][166][167][168]. and different strain rates. Figure 8 shows the determination of JC constants for the T24 steel alloy, in which is the yield stress, which has been found to have a value of 100 MPa, and ln = 4.3567, which provides a of 78 MPa, while has the value of the slope, which is found to be 0.2742 (cf. Figure 8a). Constant is determined as 0.08 (cf. Figure 8b), while constant has a value of 0.5847 (cf. Figure 8c). The ability of the JC model to accurately predict the flow behavior of T24 steel alloy is assessed using the correlation coefficient (R), which provides a value of 0.962 (cf. Figure 8d). One of the most powerful advantages of the JC model is that it is implemented in finite element software packages, in which it can be used to mimic real hot working processes and applications in severe conditions using simulation analysis. Mosleh et al. [161] implemented the JC model and two other models with finite element simulations (FES) to optimize the superplastic-forming of Ti-6%Al-4%V titanium alloy. Using FES, the optimized hot working parameters, strain rate, and temperate were determined as 0.002 s −1 and 875 °C . The FES was also helpful in checking the behavior of the superplasticforming of the tested alloy (cf. Figure 9). Many other simulation analyses have been achieved based on the JC model [162][163][164][165][166][167][168].

Modified Johnson-Cook-Based Models
Regardless of the simplicity of the JC model, which contains only five constants, the three parts that constitute the model are independently connected. Many criticisms have been reported for the JC model, since strain hardening, strain rate, and softening are

Modified Johnson-Cook-Based Models
Regardless of the simplicity of the JC model, which contains only five constants, the three parts that constitute the model are independently connected. Many criticisms have been reported for the JC model, since strain hardening, strain rate, and softening are interrelated. In reality, especially with the more complex non-linear behavior of flow stress, hot working parameters, i.e., strain, strain rate, and temperature, are connected [73,76,84,85,90,96,97,169]. Therefore, many modified JC-based models have been introduced in order to precisely predict flow behavior. In the following subsections, σ, ε, ε · , ε · 0 , ε · * , T * , T, T r , and T m are defined as they were used and explained in the original JC model (see Section 3.1).

Meyers JC-Based Modification
Meyers et al. [67] modified the original JC model by replacing the softening term with an exponential term that depends on the temperature rise. The modified model that was presented by Meyers et al. [67] can be expressed as where A, B, n, C, and λ are material constants. Constants A, B, n, and C are determined in the same way that was explained for the original JC model (see Section 3.1). Taking the logarithm after performing some rearrangements, at different values of the strain rate, Equation (7) can be written as By plotting ln σ/[(A + Bε n )(1 + C ln ε · * )] vs. T − T r , different values of constant λ can be obtained as the slope of the equation. The average of these values can be taken to determine the value of constant λ.

Rule JC-Based Modification
Rule and Jones [68] presented a modification for the original JC model to capture the sudden rise in strength at very high strain rates (over 10 4 s −1 ). This modified JC model was applied to OFHC copper, 7075-T6 aluminum, wrought iron, and high-strength steel. The modified JC model that was introduced by Rule and Jones [68] can be expressed as where C 1 , C 2 , N, C 3 , C 4 , C 5 , and M are material constants. Rule and Jones [68] modified the well-known finite element code known as elastic-plastic impact computation (EPIC) and determined the whole material's constants (N, M and C 1 − C 5 ) at once. The predicted stresses obtained using the modified JC model that was presented by Rule and Jones [68] are compared to yield strength at different strain rate values, obtained using a quasi-static experiment, and yield strength anticipated from a one-dimensional Taylor specimen model at very high strain rates for the tested high-strength steel (cf. Figure 10). As can be seen, a good agreement is obtained.
well-known finite element code known as elastic-plastic impact computation (EPIC) and determined the whole material's constants ( , and 1 − 5 ) at once. The predicted stresses obtained using the modified JC model that was presented by Rule and Jones [68] are compared to yield strength at different strain rate values, obtained using a quasi-static experiment, and yield strength anticipated from a one-dimensional Taylor specimen model at very high strain rates for the tested high-strength steel (cf. Figure 10). As can be seen, a good agreement is obtained. Figure 10. Predicted stresses obtained using the modified JC (RJC) model that was presented by Rule and Jones [68] compared to yield strength in the quasi-static experiment and anticipated from a one-dimensional Taylor specimen model at very high strain rates.

Kang JC-Based Modification
Kang et al. [69] modified the strain rate term in the original JC model to predict the dynamic behavior of a vehicle body at very high strain rates. The modified JC model that was presented by Kang et al. [69] can be expressed as where , , , 1 , 2 , and are material constants. Constants , , , and are determined in the same way that was explained for the original JC model (see Section 3.1). At the reference temperature, after performing some rearrangements, Equation (10) turns into Figure 10. Predicted stresses obtained using the modified JC (RJC) model that was presented by Rule and Jones [68] compared to yield strength in the quasi-static experiment and anticipated from a one-dimensional Taylor specimen model at very high strain rates.

Kang JC-Based Modification
Kang et al. [69] modified the strain rate term in the original JC model to predict the dynamic behavior of a vehicle body at very high strain rates. The modified JC model that was presented by Kang et al. [69] can be expressed as where A, B, n, C 1 , C 2 , and m are material constants. Constants A, B, n, and m are determined in the same way that was explained for the original JC model (see Section 3.1). At the reference temperature, after performing some rearrangements, Equation (10) turns into Figure 11 shows the initial yield stress vs. the strain rate for both the original JC and the modified JC that was presented by Kang et al. [69] for steel plate cold-rolled with grade C (SPCC) (cf. Figure 11a) and steel plate rolled with steel grade C (SPRC) (cf. Figure 11b). The figure shows that quadratic fitting can be used with the modified JC model. Similar findings can be obtained using the modified JC model introduced by Kang et al. [69] to investigate the dynamic behavior of autobody sheet metal [70]. Consequently, by plotting σ/(A + Bε n ) vs. ln ε · * , constants C 1 and C 2 can be determined with quadratic fitting.

Couque JC-Based Modification
Couque et al. [71] presented another modification for the strain rate term in the original JC model, in which a power constant is correlated with the strain rate. The modified JC model that was presented by Couque et al. [71] can be written as where A, B, n, C 1 , k, and m are material constants. The strain rate, ε · 1 , defines a transition between thermally activated and viscous regimes and is stated to have 0.001 s −1 [71]. Constants A, B, n, and m are determined in the same way that was explained for the original JC model (see Section 3.1). At the reference temperature, after performing some rearrangements, Equation (12) can be written as Constants C 1 , E, and k can be determined using algorithms that are based on non-linear least square methods. Good fitting between experimental stresses and predicted stresses are obtained using the modified JC method that was presented by Couque et al. [71] for the prediction of flow stress in nickel at different strain rates, and the home temperature can be obtained (cf. Figure 12). ( + ) = 1 + 1 ln • * + 2 (ln • * ) 2 (11) Figure 11 shows the initial yield stress vs. the strain rate for both the original JC and the modified JC that was presented by Kang et al. [69] for steel plate cold-rolled with grade C (SPCC) (cf. Figure 11a) and steel plate rolled with steel grade C (SPRC) (cf. Figure 11b). The figure shows that quadratic fitting can be used with the modified JC model. Similar findings can be obtained using the modified JC model introduced by Kang et al. [69] to investigate the dynamic behavior of autobody sheet metal [70]. Consequently, by plotting /( + ) vs. ln • * , constants 1 and 2 can be determined with quadratic fitting.
. Figure 11. Initial yield stress vs. strain rate using both the original JC model and the modified JC model that was presented by Kang et al. [69] for (a) SPCC and (b) SPRC.

Couque JC-Based Modification
Couque et al. [71] presented another modification for the strain rate term in the original JC model, in which a power constant is correlated with the strain rate. The modified JC model that was presented by Couque et al. [71] can be written as where , , , 1 , , and are material constants. The strain rate, • 1 , defines a transition between thermally activated and viscous regimes and is stated to have 0.001 s −1 [71]. Constants , , , and are determined in the same way that was explained for the original JC model (see Section 3.1). At the reference temperature, after performing some rearrangements, Equation (12) can be written as Constants 1 , , and can be determined using algorithms that are based on nonlinear least square methods. Good fitting between experimental stresses and predicted stresses are obtained using the modified JC method that was presented by Couque et al. [71] for the prediction of flow stress in nickel at different strain rates, and the home temperature can be obtained (cf. Figure 12).  A similar modification for the strain rate term is presented by Johnson et al. [72], in which the strain rate term is introduced by Figure 12. Experimental stresses (markers) compared to predicted stresses (solid lines) using the modified JC model that was presented by Couque et al. [71] for nickel.
A similar modification for the strain rate term is presented by Johnson et al. [72], in which the strain rate term is introduced by where constants C 1 , C 2 , and C 3 can be determined by using the exponential fitting of σ/(A + Bε n ) vs. ln ε · * .

Lin (1) JC-Based Modification
Lin et al. [73] introduced one of the famous modifications of the original JC model to study the flow behavior of typical high-strength alloy steel. The coupling effect between the hot working parameters, strain, strain rate, and temperature is taken into consideration in the proposed modification. The modified JC model that was introduced by Lin et al. [73] can be expressed as where A, B 1 , B 2 , C 1 , λ 1 , and λ 2 are material constants. As can be seen from Equation (14), the first term in the JC model is replaced with a quadratic function of strain, with constants A, B 1 , and B 2 . In this modification, the strain rate term in the original JC is employed as it is, with constant C 1 . The third term in the original JC model is modified considering the coupling effect between the strain rate and temperature considering two constants, λ 1 and λ 2 . At the reference strain rate and reference temperature, Equation (14) reduces to where constants A, B 1 , and B 2 are determined by fitting the experimental data of σ vs. ε at the reference values with the quadratic function in the strain (cf. Figure 13a). At the reference temperature, Equation (14) lowers to After performing some rearrangements, Equation (16) can be written as Strain rate constant 1 is the slope of Equation (17) and is obtained by plotting /( + 1 + 2 2 ) vs. ln • * (cf. Figure 13b). To obtain the value of the two constants, 1 and 2 , a new parameter, , that is equal to 1 + 2 ln • * is introduced, after performing some rearrangements, Equation (14) can be expressed as By taking the logarithm of both sides, Equation (18) can be written as After performing some rearrangements, Equation (16) can be written as Strain rate constant C 1 is the slope of Equation (17) and is obtained by plotting σ/ A + B 1 ε + B 2 ε 2 vs. ln ε · * (cf. Figure 13b).
To obtain the value of the two constants, λ 1 and λ 2 , a new parameter, λ, that is equal to λ 1 + λ 2 ln ε · * is introduced, after performing some rearrangements, Equation (14) can be expressed as By taking the logarithm of both sides, Equation (18) can be written as By plotting ln σ/ A + B 1 ε + B 2 ε 2 (1 + C 1 ln ε · * ) vs. (T − T r ) at different strain rate and temperature values, different values for λ (slope of the equation; cf. Figure 14a-c) are obtained, which can be plotted vs. ln ε · * , in which the value of λ 1 is the intercept and the value of λ 2 is the slope (Figure 14d). Experimental stresses are compared to predicted stresses that were obtained by the modified JC model that was introduced by Lin et al. [73] for typical high-strength alloy steel at different strain rates and different temperatures, as shown in Figure 15. The figure shows that the predicted stresses agree very well with the experimental stresses, with a maximum relative error (RE) of 5.15%. The modified JC model that was introduced by Lin et al. [73] has been used with many different alloys [170][171][172][173][174], with accurate predictions and inaccurate predictions for other alloys [61,96,[175][176][177]. Experimental stresses are compared to predicted stresses that were obtained by the modified JC model that was introduced by Lin et al. [73] for typical high-strength alloy steel at different strain rates and different temperatures, as shown in Figure 15. The figure shows that the predicted stresses agree very well with the experimental stresses, with a maximum relative error (RE) of 5.15%. The modified JC model that was introduced by Lin et al. [73] has been used with many different alloys [170][171][172][173][174], with accurate predictions and inaccurate predictions for other alloys [61,96,[175][176][177]. Experimental stresses are compared to predicted stresses that were obtained by the modified JC model that was introduced by Lin et al. [73] for typical high-strength alloy steel at different strain rates and different temperatures, as shown in Figure 15. The figure shows that the predicted stresses agree very well with the experimental stresses, with a maximum relative error (RE) of 5.15%. The modified JC model that was introduced by Lin et al. [73] has been used with many different alloys [170][171][172][173][174], with accurate predictions and inaccurate predictions for other alloys [61,96,[175][176][177].  [74] modified the softening term in the original JC model, considering that the tested temperature might be higher or lower than the reference temperature. The modified JC model that was presented by Hou Q. Y. et al. [74] can be expressed as where constants A, B, n, and C are determined as was explained for the original JC model (see Section 3.1). At different strain rate values, after performing some rearrangements, Equation (20) can be written as By plotting 1 − σ/(A + Bε n ) vs. e T/T m − e T r /T m / e − e T r /T m , at different values of strain, different values of constant λ (slope of the equation) can be obtained, in which the average can be implemented.
A good agreement between experimental stresses and predicted stresses obtained by the modified JC that was presented by Hou Q. Y. et al. [74] for Mg-10Gd-2Y-0.5Zr alloy is achieved (cf. Figure 16).
Perez et al. [178] used the modified JC that was introduced by Hou Q. Y. et al. [74] for the prediction of the flow behavior of Ti6Al4V alloy at high temperatures. The modified model succeeded in predicting the flow behavior of the tested titanium alloy with an R of 0.9765 compared with other models. In fact, considering Equation (21), λ varies with strain, and, hence, Guoliang et al. [179] modified Equation (21) by replacing λ with λ ε , in that λ has a constant value. Figure 17 shows the regression of both λ and λ by plotting e T/T m − e T r /T m / e − e T r /T m vs. 1 − σ/(A + Bε n ) and e T/T m − e T r /T m / e − e T r /T m vs. (σ/(A + Bε n ))ε, respectively, which was performed by Guoliang et al. [179]. A good agreement between experimental stresses and predicted stresses obtained by the modified JC that was presented by Hou Q. Y. et al. [74] for Mg-10Gd-2Y-0.5Zr alloy is achieved (cf. Figure 16). Perez et al. [178] used the modified JC that was introduced by Hou Q. Y. et al. [74] for the prediction of the flow behavior of Ti6Al4V alloy at high temperatures. The modified model succeeded in predicting the flow behavior of the tested titanium alloy with an R of 0.9765 compared with other models. In fact, considering Equation (21), varies with strain, and, hence, Guoliang et al. [179] modified Equation (21) by replacing with ′ ⁄ , in that ′ has a constant value. Figure

Shin JC-Based Modification
Shin and Kim [75] modified the original JC model for the accurate prediction of materials at wide regimes of temperature and strain rate. The well-known Voce hardening [180] is implemented in the modified JC model since it is describing the saturation of strain hardening. Thus, the strain rate term is also modified, so the high increase in the stress due to large strain rates can be taken into consideration. Furthermore, the temperature softening constant is applied to the whole bracket instead of * . The modified JC model that was presented by Shin and Kim [75] can be expressed as where constants and represent yield and saturated stress, constant represents the strain-hardening exponent, and constants and correlate to the strain rate, while constant correlates to softening. At the reference strain rate and reference temperature, the Voce hardening is the only left term in Equation (22), as = [ + {1 − exp(− )}], in which can be measured as the yield stress, and and can be computed using a regression analysis that is based on the non-linear least square method. Similarly, constants , , and m can be determined. Figure 18 shows that the predicted stresses of the modified JC model that was introduced by Shin and Kim [75] have a very good agreement with the experimental stresses for tungsten heavy alloy at high temperatures and different strain rates, with R = 0.995 and AARE = 1.98% (cf. Figure 18a) and with R = 0.991 and AARE = 1.82% (cf. Figure 18b).

Shin JC-Based Modification
Shin and Kim [75] modified the original JC model for the accurate prediction of materials at wide regimes of temperature and strain rate. The well-known Voce hardening [180] is implemented in the modified JC model since it is describing the saturation of strain hardening. Thus, the strain rate term is also modified, so the high increase in the stress due to large strain rates can be taken into consideration. Furthermore, the temperature softening constant m is applied to the whole bracket instead of T * . The modified JC model that was presented by Shin and Kim [75] can be expressed as where constants A and B represent yield and saturated stress, constant C represents the strain-hardening exponent, and constants D and E correlate to the strain rate, while constant m correlates to softening. At the reference strain rate and reference temperature, the Voce hardening is the only left term in Equation (22), as σ = [A + B{1 − exp(−Cε)}], in which A can be measured as the yield stress, and B and C can be computed using a regression analysis that is based on the non-linear least square method. Similarly, constants D, E, and m can be determined. Figure 18 shows that the predicted stresses of the modified JC model that was introduced by Shin and Kim [75] have a very good agreement with the experimental stresses for tungsten heavy alloy at high temperatures and different strain rates, with R = 0.995 and AARE = 1.98% (cf. Figure 18a) and with R = 0.991 and AARE = 1.82% (cf. Figure 18b). constant correlates to softening. At the reference strain rate and reference temperature, the Voce hardening is the only left term in Equation (22), as = [ + {1 − exp(− )}], in which can be measured as the yield stress, and and can be computed using a regression analysis that is based on the non-linear least square method. Similarly, constants , , and m can be determined. Figure 18 shows that the predicted stresses of the modified JC model that was introduced by Shin and Kim [75] have a very good agreement with the experimental stresses for tungsten heavy alloy at high temperatures and different strain rates, with R = 0.995 and AARE = 1.98% (cf. Figure 18a) and with R = 0.991 and AARE = 1.82% (cf. Figure 18b).  In another published article, Shin and Kim [181] studied the capability of the Shin and Kim [75] modified JC model for the prediction of the flow behavior of copper at a wide range of strain rates and temperatures. Compared with the original JC model, the mechanical threshold (MTS) model, and the Preston-Tonks-Wallace (PTW) model, the Shin and Kim (SK) model [75] provided the best predictions for the flow behavior of copper, with R = 0.931 and AARE = 14.01%, as shown in Figure 19. In another published article, Shin and Kim [181] studied the capability of the Shin and Kim [75] modified JC model for the prediction of the flow behavior of copper at a wide range of strain rates and temperatures. Compared with the original JC model, the mechanical threshold (MTS) model, and the Preston-Tonks-Wallace (PTW) model, the Shin and Kim (SK) model [75] provided the best predictions for the flow behavior of copper, with R = 0.931 and AARE = 14.01%, as shown in Figure 19.

Maheshwari JC-Based Modification
Maheshwari et al. [76] modified the original JC material flow model during hot deformation. The modified JC model can be expressed as where σ y is reference stress, and σ m is the true stress upon melting, which has a reference value of zero. Constants P, Q, n, r, α, and β are six material constants that constitute the modified model. At the reference strain rate and reference temperature, Equation (23) lowers to σ = P + Qε n , with P equal to yield stress. Constants Q and n can be determined by fitting the obtained linear logarithmic equation, as explained in the original JC model (see Section 3.1). Taking logarithms for both sides of Equation (23), after performing some rearrangements, and at a reference temperature, the equation can be expressed as By plotting ln[σ/(P + Qε n )] vs. ln ε · /ε · • , constant r has the value of the slope. Finally, at different strain rate and temperature values, after performing some rearrangement and taking logarithms for both sides, Equation (23) can be expressed as , the two constants α and β can be computed by fitting the data with an exponential function.
Compared with the original JC model, the modified JC that was introduced by Maheshwari et al. [76] gave accurate predictions for some of the strain rate and temperature combinations, while it gave inaccurate predictions with other combinations. A comparison between experimental stresses and predicted stresses obtained by the original JC and by the modified JC that was presented by Maheshwari et al. [76] for Al-2024 alloy under hot deformation is shown in Figure 20. As can be seen, the modified model cannot guarantee precise predictions of stresses.
The predictability of the flow behavior using the modified JC that is presented by Maheshwari et al. [76], along with other models, was investigated for 2024Al alloy under hot deformations by Trimble and O'Donnell [182]. The results showed that the model provides good predictions of the flow behavior of the tested alloy, with an R of 0.992 and an average absolute relative error (AARE) of 5.9491%. However, it failed to accurately predict the flow behavior in some combinations between the strain rate and temperature. The same findings were obtained by Maheshwari [183].

Wang (1) JC-Based Modification
Wang et al. [77] modified the JC model to predict the flow behavior of 30Cr2Ni4MoV rotor steel alloy through a large range of strain rates and temperatures. The coupling effect between the strain, strain rate, and temperature is taken into account in this modification. In this modification, the Voce hardening [180] equation is used to describe the strain-hardening term. The modified JC model that was presented by Wang et al. [77] can be written as where T L is defined as the lowest temperature in the tested range of temperatures, and A, B 0 , B 1 , C 1 , C 2 , λ 1 , and λ 2 are the material constants. At the reference strain rate and lowest temperature, Equation (26) Taking into account that A represents the yield stress, constants B 0 and B 1 can be determined using regression analysis by plotting σ vs. ε. At the lowest temperature, after performing some rearrangements with C = C 1 + C 2 ε, Equation (26) can be expressed as rearrangement and taking logarithms for both sides, Equation (23) can be expressed as By plotting ln[ /(( + )( • /°•) )] vs. ( − )/( − ), the two constants and can be computed by fitting the data with an exponential function. Compared with the original JC model, the modified JC that was introduced by Maheshwari et al. [76] gave accurate predictions for some of the strain rate and temperature combinations, while it gave inaccurate predictions with other combinations. A comparison between experimental stresses and predicted stresses obtained by the original JC and by the modified JC that was presented by Maheshwari et al. [76] for Al-2024 alloy under hot deformation is shown in Figure 20. As can be seen, the modified model cannot guarantee precise predictions of stresses. The predictability of the flow behavior using the modified JC that is presented by Maheshwari et al. [76], along with other models, was investigated for 2024Al alloy under hot deformations by Trimble and O'Donnell [182]. The results showed that the model  Figure 21a). Hence, C can be linearly plotted vs. ε, in which the value of C 1 is the intercept and the value of C 2 is the slope (cf. Figure 21b).
At different strain rate values, after performing some rearrangements, Equation (26) can be written as Taking the logarithm of both sides and introducing a new parameter, λ = λ 1 + λ 2 ln ε · * , the equation turned out to be linear with a slope of λ. The two constants, λ 1 and λ 2 , can be determined as explained in Section 3.2.5.
A comparison between experimental stresses and predicted stresses obtained using the modified JC model that was presented by Wang et al. [77] for 30Cr2Ni4MoV rotor steel alloy at elevated temperatures and different strain rates is shown in Figure 22. As can be seen, the predicted stresses obtained by the modified JC model have very good agreements with the experimental stresses. This might be due to the use of Voce hardening, as well as taking the coupling effect between the strain rate and temperature into consideration.

that
represents the yield stress, constants 0 and 1 can be determined using regression analysis by plotting vs. . At the lowest temperature, after performing some rearrangements with = 1 + 2 , Equation (26) can be expressed as By plotting /[ − 0 exp(− 1 )] vs. ln( • * ) at different strain and strain rate values, different values of can be obtained (cf. Figure 21a). Hence, can be linearly plotted vs. , in which the value of 1 is the intercept and the value of 2 is the slope (cf. Figure 21b). At different strain rate values, after performing some rearrangements, Equation (26) can be written as Taking the logarithm of both sides and introducing a new parameter, = 1 + 2 ln • * , the equation turned out to be linear with a slope of . The two constants, 1 and 2 , can be determined as explained in Section 3.2.5. A comparison between experimental stresses and predicted stresses obtained using the modified JC model that was presented by Wang et al. [77] for 30Cr2Ni4MoV rotor steel alloy at elevated temperatures and different strain rates is shown in Figure 22. As can be seen, the predicted stresses obtained by the modified JC model have very good agreements with the experimental stresses. This might be due to the use of Voce hardening, as well as taking the coupling effect between the strain rate and temperature into consideration. Figure 22. Experimental stresses (markers) compared to predicted stresses (solid lines) using the modified JC that was presented by Wang et al. [77] for 30Cr2Ni4MoV rotor steel at (a) 1323 K and (b) 1423 K.

Lin (2) JC-Based Modification
Lin et al. [78] modified the original JC model to predict the flow behavior of Al-Zn-Mg-Cu alloy under hot deformation. The modification considers the coupling effect of the strain, strain rate, and temperature. The modified JC model that was presented by Lin et al. [78] can be written as In this modification, and in order to overcome the difficulty of finding yield stress, the yield stress is replaced by beak stress and defined using an Arrhenius-type equation as per [78]: where = • exp( / ); is the activation energy; and is the universal gas constant, which has a value of 8.31 Jmol −1 ·K −1 . , , , and are constants to be determined; see [78] for more information. By taking the logarithm after performing some rearrangements, at the reference temperature, Equation (29) can be expressed as By plotting ln[ − 0 ] vs. ln ε, the parameters of ( • ) can be obtained from the intercept, and the parameters of ( • ) can be obtained from the slope; then, using Figure 22. Experimental stresses (markers) compared to predicted stresses (solid lines) using the modified JC that was presented by Wang et al. [77] for 30Cr2Ni4MoV rotor steel at (a) 1323 K and (b) 1423 K.

Lin (2) JC-Based Modification
Lin et al. [78] modified the original JC model to predict the flow behavior of Al-Zn-Mg-Cu alloy under hot deformation. The modification considers the coupling effect of the strain, strain rate, and temperature. The modified JC model that was presented by Lin et al. [78] can be written as In this modification, and in order to overcome the difficulty of finding yield stress, the yield stress is replaced by beak stress and defined using an Arrhenius-type equation as per [78]: where Z = ε · exp(Q/RT); Q is the activation energy; and R is the universal gas constant, which has a value of 8.31 Jmol −1 ·K −1 . Q, α, A, and m are constants to be determined; see [78] for more information. By taking the logarithm after performing some rearrangements, at the reference temperature, Equation (29) can be expressed as By plotting ln[σ − σ 0 ] vs. ln ε, the parameters of B(ε · ) can be obtained from the intercept, and the parameters of n(ε · ) can be obtained from the slope; then, using polynomial fitting for B(ε · ) with ε · and for n(ε · ) with ε · the constants that constitute their equations can be obtained. By taking the logarithm after performing some rearrangements, at different strain rates, Equation (29) can be written as At different values for both the temperature and strain rate, different values of p(ε · ) can be obtained; then, the polynomial fitting of p(ε · ) with ε · can be used to determine the parameters of p(ε · ) (see Section 3.2.5).
Lin et al. [78] fitted B(ε · ), n(ε · ), and p(ε · ) vs. ε · with linear relationships for Zn-Mg-Cu alloy, which can be expressed as B(ε · ) = 5.69337 A comparison between experimental stresses and predicted stresses obtained using the modified JC model that was introduced by Lin et al. [78] for the prediction of Zn-Mg-Cu alloy is shown in Figure 23. As can be seen, very accurate predictions were obtained, which might be due to the coupling effect of strain, strain rate, and temperature that was taken into consideration. A comparison between experimental stresses and predicted stresses obtained using the modified JC model that was introduced by Lin et al. [78] for the prediction of Zn-Mg-Cu alloy is shown in Figure 23. As can be seen, very accurate predictions were obtained, which might be due to the coupling effect of strain, strain rate, and temperature that was taken into consideration.

Lin (3) JC-Based Modification
Lin et al. [79] presented another modification of the original JC model to predict the flow behavior of 7075 Al alloy under hot deformation. The modification considers the coupling effect of strain, strain rate, and forming temperature. The modified JC model that was presented by Lin et al. [79] can be written as 0 is represented by using an Arrhenius-type equation, as in Section 3.2.10, the same as Equation (30). By following the same procedures in Section 3.2.10, the parameters of ( ), ( ) and ( ) can be obtained using fitting analysis. Lin et al. [79]  A comparison between the experimental stresses and predicted stresses obtained using the modified JC model that was presented by Lin et al. [79] for the prediction of the flow behavior of 7075 Al alloy is shown in Figure 24. As can be seen, very accurate

Lin (3) JC-Based Modification
Lin et al. [79] presented another modification of the original JC model to predict the flow behavior of 7075 Al alloy under hot deformation. The modification considers the coupling effect of strain, strain rate, and forming temperature. The modified JC model that was presented by Lin et al. [79] can be written as σ 0 is represented by using an Arrhenius-type equation, as in Section 3.2.10, the same as Equation (30). By following the same procedures in Section 3.2.10, the parameters of B(T), n(T) and C(T) can be obtained using fitting analysis.
A comparison between the experimental stresses and predicted stresses obtained using the modified JC model that was presented by Lin et al. [79] for the prediction of the flow behavior of 7075 Al alloy is shown in Figure 24. As can be seen, very accurate predictions were obtained, which might be due to the coupling effect of strain, strain rate, and temperature that was taken into consideration.

Li JC-Based Modification
Li et al. [80] introduced a modification of the original JC model for the prediction of the flow behavior of T24 steel at elevated temperatures and different strain rates. In this modification, both strain-hardening and -softening terms are modified, while the strain rate term is kept the same as in the original JC model. At reference values, the experimental stresses vs. strains are fitted to second-order polynomial function, while the softening term takes the coupling effect of strain and temperature, as well as strain rate and temperature into account. The modified JC model that was presented by Li et al. [80] can be expressed as where , 1 , 2 , , , and are material constants, taking into account that constant depends on , and constant depends on • and . Constants , 1 , 2 , and can be determined in the same way that was explained in Section 3.2.5. At the reference strain rate, Equation (34) turns into Taking logarithms for both sides after performing some rearrangements, Equation (35) can be written as The slope of the equation is ( ); since it depends on the strain, different values for can be obtained at different values of strain. Then, a polynomial fit is used to determine as a function of . Constant V( • , ) can be determined by introducing a new parameter, , which is equal to ( ) + V( • , ) ln • * . After performing some rearrangements, at the remaining values of the strain rates, Equation (34) can be expressed as By taking logarithms for both sides, Equation (37) can be introduced as

Li JC-Based Modification
Li et al. [80] introduced a modification of the original JC model for the prediction of the flow behavior of T24 steel at elevated temperatures and different strain rates. In this modification, both strain-hardening and -softening terms are modified, while the strain rate term is kept the same as in the original JC model. At reference values, the experimental stresses vs. strains are fitted to second-order polynomial function, while the softening term takes the coupling effect of strain and temperature, as well as strain rate and temperature into account. The modified JC model that was presented by Li et al. [80] can be expressed as where A, B 1 , B 2 , C, Q, and V are material constants, taking into account that constant Q depends on ε, and constant V depends on ε · and ε. Constants A, B 1 , B 2 , and C can be determined in the same way that was explained in Section 3.2.5. At the reference strain rate, Equation (34) turns into Taking logarithms for both sides after performing some rearrangements, Equation (35) can be written as The slope of the equation is Q(ε); since it depends on the strain, different values for Q can be obtained at different values of strain. Then, a polynomial fit is used to determine Q as a function of ε. Constant V(ε · , ε) can be determined by introducing a new parameter, S, which is equal to Q(ε) + V(ε · , ε) ln ε · * . After performing some rearrangements, at the remaining values of the strain rates, Equation (34) can be expressed as By taking logarithms for both sides, Equation (37) can be introduced as Different values for the slope, S, are obtained with different values for the temperature, strain, and strain rate; then, different values for V(ε · , ε), can be determined as a function of the strain and strain rate by inputting V(ε · , ε) = [S − Q(ε)]/ ln ε · * . For T24 steel, Li et al. [80] determined the polynomial function of Q as −0.00344 + 0.00177ε − 0.00135ε 2 , while the values of V at different strain and strain rate values are shown in Table 1. Good predictability for the prediction of the flow stress of T24 steel using the modified JC model that was presented by Li et al. [80] is obtained (cf. Figure 25b,d) when compared with those obtained using the original JC model (cf. Figure 25a,c). The modified JC model has an R-value of 0.991 compared with 0.962 for the original JC model. In addition, it has an AARE of 5.37% compared with 9.41% for the original JC model.

Song JC-Based Modification
Song et al. [81] modified the original JC model to predict the flow behavior of titanium matrix composite under hot deformation. The softening term is modified to take the strain rate and temperature effect into account. The modified model can be expressed as where A, B 1 , B 2 , C, λ 1 , and λ 2 are material constants. Constants A, B 1 , B 2 , and C can be determined in the same way that was explained in Section 3.2.5. At different strain rate values, Equation (39) can be expressed as Taking logarithms for both sides, after performing some rearrangements, Equation (40) can be introduced as By plotting the left side vs. T * , different values for λ 1 and λ 2 can be obtained at different strain rate values, for which the average of each one can be obtained (cf. Figure 26). Good predictability for the prediction of the flow stress of T24 steel using the modified JC model that was presented by Li et al. [80] is obtained (cf. Figure 25b,d) when compared with those obtained using the original JC model (cf. Figure 25a,c). The modified JC model has an R-value of 0.991 compared with 0.962 for the original JC model. In addition, it has an AARE of 5.37% compared with 9.41% for the original JC model.

Song JC-Based Modification
Song et al. [81] modified the original JC model to predict the flow behavior of titanium matrix composite under hot deformation. The softening term is modified to take the strain rate and temperature effect into account. The modified model can be expressed as where , 1 , 2 , , 1 , and 2 are material constants. Constants , 1 , 2 , and can be determined in the same way that was explained in Section 3.2.5. At different strain rate values, Equation (39) can be expressed as Taking logarithms for both sides, after performing some rearrangements, Equation (40) can be introduced as By plotting the left side vs. * , different values for 1 and 2 can be obtained at different strain rate values, for which the average of each one can be obtained (cf. Figure  26). The modified JC model that was introduced by Song et al. [81] provided a good agreement between predicted stresses and experimental stresses when compared with those obtained using the original JC model (cf. Figure 27). The modified JC model that was introduced by Song et al. [81] provided a good agreement between predicted stresses and experimental stresses when compared with those obtained using the original JC model (cf. Figure 27). The modified JC model that was introduced by Song et al. [81] provided a good agreement between predicted stresses and experimental stresses when compared with those obtained using the original JC model (cf. Figure 27).

Wang (2) JC-Based Modification
Wang et al. [48] modified the strain rate term in the original JC model to predict the flow behavior of Inconel 718 at high strain rates and high temperatures. In this modification, the strain rate has been found to be correlated with temperature in a sine wave function. The modified JC model that was presented by Wang et al. [48] can be written as where A, B, n, and m are material constants that are determined in the same way that was explained for the original JC model (see Section 3.1). Constant C for Inconel 718 is represented by the sine wave function (cf. Figure 28a Wang et al. [48] modified the strain rate term in the original JC model to predict the flow behavior of Inconel 718 at high strain rates and high temperatures. In this modification, the strain rate has been found to be correlated with temperature in a sine wave function. The modified JC model that was presented by Wang et al. [48] can be written as where , , , and are material constants that are determined in the same way that was explained for the original JC model (see Section 3.1). Constant for Inconel 718 is represented by the sine wave function (cf. Figure 28a Good agreement between experimental stresses and predicted stresses using the modified model for the flow behavior of Inconel 718 is thus obtained. However, this method does not provide accurate or precise predictions (Figure 28b). Good agreement between experimental stresses and predicted stresses using the modified model for the flow behavior of Inconel 718 is thus obtained. However, this method does not provide accurate or precise predictions (Figure 28b).
Another sine wave approximation of constant C(ε · , T) was introduced by Xu et al. [88] to predict the flow behavior of SnSbCu alloy at a wide range of strain rates and different temperatures, with good predictions, in which C(ε · , T) is determined as C(ε · , T) = 0.06 − 0.0232 0.05554 + 0.01777 sin ε · − 1000 5000 π sin T − 20 300 π

Tan JC-Based Modification
Tan et al. [53] modified the JC model to predict the flow behavior of 7050-T7451 aluminum alloy under dynamic loading. Both strain hardening and softening are kept constant, as in the original JC model, while the strain rate term is modified to take the coupling effect of the strain and strain rate into account. The modified JC model that was introduced by Tan et al. [53] can be written as where A, B, n, and m are material constants that are determined in the same way that was explained for the original JC model (see Section 3.1). The constant, C(ε, ε · ), is obtained by plotting C = [σ/((A + Bε n ) − 1)]/ ln ε · * vs. both ε and ε · (cf. Figure 29) in the original JC model at reference temperature and using polynomial fitting.
where constants 0 , 1 , 2 , 3 , 4 , and 5 are determined using regression analysis. Compared to the predictions of the original JC model (cf. Figure 30a), the modified JC model that was presented by Tan et al. [53] gave very good predictions for the flow behavior of the studied alloy at elevated temperatures (cf. Figure 30b). The good prediction of the flow stress using the modified model might be due to taking the interaction between the strain and strain rate into account. Constant C(ε, ε · ) is determined by Tan et al. [53] as: where constants C 0 , C 1 , C 2 , C 3 , C 4 , and C 5 are determined using regression analysis. Compared to the predictions of the original JC model (cf. Figure 30a), the modified JC model that was presented by Tan et al. [53] gave very good predictions for the flow behavior of the studied alloy at elevated temperatures (cf. Figure 30b). The good prediction of the flow stress using the modified model might be due to taking the interaction between the strain and strain rate into account.

Chen JC-Based Modification
Chen et al. [82] modified the original JC model to predict the flow behavior of 7050-T745 aluminum alloy at high strain rates and different temperatures. The coupling effect of strain, strain rate, and temperature is taken into account in this modification. The modified JC model that was introduced by Chen et al. [82] can be expressed as Material constants A, B, n, C, m, p1, p2, p3, and q are determined using a generic algorithm developed by Chen et al. [184] to minimize the mean square error between experimental stresses and predicted stresses. ( , • ) = 0 + 1 + 2 2 + 3 ln • * + 4 (ln • * ) 2 + 5 ln • * (44) where constants 0 , 1 , 2 , 3 , 4 , and 5 are determined using regression analysis. Compared to the predictions of the original JC model (cf. Figure 30a), the modified JC model that was presented by Tan et al. [53] gave very good predictions for the flow behavior of the studied alloy at elevated temperatures (cf. Figure 30b). The good prediction of the flow stress using the modified model might be due to taking the interaction between the strain and strain rate into account.

Chen JC-Based Modification
Chen et al. [82] modified the original JC model to predict the flow behavior of 7050-T745 aluminum alloy at high strain rates and different temperatures. The coupling effect of strain, strain rate, and temperature is taken into account in this modification. The modified JC model that was introduced by Chen et al. [82] can be expressed as The predicted stresses obtained using the modified JC model that was presented by Chen et al. [82] were found to be in good agreement with experimental stresses for 7050-T745 aluminum alloy at high strain rates and different temperatures for the tested alloy with RE less than 5%, which might be attributed to the interrelation between the hot working parameters (cf. Figure 31).
Material constants , , , , , 1, 2, 3 , and are determined using a generic algorithm developed by Chen et al. [184] to minimize the mean square error between experimental stresses and predicted stresses.
The predicted stresses obtained using the modified JC model that was presented by Chen et al. [82] were found to be in good agreement with experimental stresses for 7050-T745 aluminum alloy at high strain rates and different temperatures for the tested alloy with RE less than 5%, which might be attributed to the interrelation between the hot working parameters (cf. Figure 31).

Wang (3) JC-Based Modification
Wang et al. [83] modified the JC model to predict the flow behavior of a nickel-based superalloy at high strain rates and different temperatures. The coupling effect of anomalous temperature and strain rate dependences on flow stress is taken into account in this modification, using the strain rate term that was first presented by Rule and Jones [68] and modifying the softening term. The modified JC model that was introduced by Wang et al. [83] can be expressed as

Wang (3) JC-Based Modification
Wang et al. [83] modified the JC model to predict the flow behavior of a nickelbased superalloy at high strain rates and different temperatures. The coupling effect of anomalous temperature and strain rate dependences on flow stress is taken into account in this modification, using the strain rate term that was first presented by Rule and Jones [68] and modifying the softening term. The modified JC model that was introduced by Wang et al. [83] can be expressed as Material constants A, B, n, C 3 , C 4 , C 5 , m, D, d , and T p are determined by creating a generic algorithm that minimizes the mean square error between experimental stresses and predicted stresses, as reported by Rule and Jones [68].
The modified JC model that was introduced by Wang et al. [83] provided good predictions for the flow behavior of the tested alloy when predicted stresses were compared to experimental stresses (cf. Figure 32).

Shokry (1) JC-Based Modification
Shokry [84] introduced a modification of the JC model for the prediction of the flow behavior of alloy 800H at elevated temperatures and intermediate strain rates. The modification takes the coupling effect between strain and both strain rate and temperature into account. The modified JC that was presented by Shokry [84] can be expressed as where constants , 1 , 2 , and 3 are determined by fitting the experimental data of stress and strain at the reference strain rate and reference temperature with three-order polynomial functions (see Section 3.2.5). Shokry [84] determined both strain rate constants 1 and 2 and both softening constants 1 and 2 at once using the Kalman filter technique, a mathematical method for the determination of constants that minimizes the mean square error between experimental and predicted stresses [185].
A comparison between experimental stresses and predicted stresses obtained by the original JC and the modified JC that was presented by Shokry [84] is shown in Figure 33. The figure shows that the modified model can predict the flow behavior of the tested alloy better than the original JC model, with an R-value of 0.98 compared with 0.91 for the original JC model. The modified JC model that was presented by Shokry [84] provided good predictions of the flow behavior of a powder metallurgy nanoquasicrystalline Al93Fe3Cr2Ti2 alloy at high temperatures and different strain rates [61] when compared with other models, with R and AARE values of 0.98 and 7.8%, respectively, compared with 0.92 and 12.7% for the original JC model, respectively.

Shokry (1) JC-Based Modification
Shokry [84] introduced a modification of the JC model for the prediction of the flow behavior of alloy 800H at elevated temperatures and intermediate strain rates. The modification takes the coupling effect between strain and both strain rate and temperature into account. The modified JC that was presented by Shokry [84] can be expressed as where constants A, B 1 , B 2 , and B 3 are determined by fitting the experimental data of stress and strain at the reference strain rate and reference temperature with three-order polynomial functions (see Section 3.2.5). Shokry [84] determined both strain rate constants C 1 and C 2 and both softening constants m 1 and m 2 at once using the Kalman filter technique, a mathematical method for the determination of constants that minimizes the mean square error between experimental and predicted stresses [185]. A comparison between experimental stresses and predicted stresses obtained by the original JC and the modified JC that was presented by Shokry [84] is shown in Figure 33. The figure shows that the modified model can predict the flow behavior of the tested alloy better than the original JC model, with an R-value of 0.98 compared with 0.91 for the original JC model. The modified JC model that was presented by Shokry [84] provided good predictions of the flow behavior of a powder metallurgy nanoquasicrystalline Al9 3 Fe 3 Cr 2 Ti 2 alloy at high temperatures and different strain rates [61] when compared with other models, with R and AARE values of 0.98 and 7.8%, respectively, compared with 0.92 and 12.7% for the original JC model, respectively.

Zhao JC-Based Modification
Zhao et al. [85] used the modification of the original JC that was introduced by Lin et al. [73] and established a sine wave function for the strain rate term, as performed in [48], to predict the dynamic flow behavior of FeCr alloy manufactured using laser-additive manufacturing. The strain-hardening term is also modified by fitting stress with strain at the reference conditions, with three-order polynomial functions instead of two-order polynomial functions. The modified JC model that was introduced by Zhao et al. [85] can be written as where , 1 , 2 , 3 , 1 , and 2 are material constants that are determined in the same way that was explained for the modified JC model that was introduced by Lin et al. [ Good agreement between experimental stresses and predicted stresses is obtained using the modified JC model that was introduced by Zhao et al. [85], with R: 0.87-0.99 and AARE: 7.6%-22.7% when compared with those obtained using the original JC model with R: 0.83-0.88 and AARE: 21.8%-44.8%. However, very accurate predictions were not found (cf. Figure 34).

Zhao JC-Based Modification
Zhao et al. [85] used the modification of the original JC that was introduced by Lin et al. [73] and established a sine wave function for the strain rate term, as performed in [48], to predict the dynamic flow behavior of FeCr alloy manufactured using laseradditive manufacturing. The strain-hardening term is also modified by fitting stress with strain at the reference conditions, with three-order polynomial functions instead of twoorder polynomial functions. The modified JC model that was introduced by Zhao et al. [85] can be written as where A, B 1 , B 2 , B 3 , λ 1 , and λ 2 are material constants that are determined in the same way that was explained for the modified JC model that was introduced by Lin et al. [48] (see Section 3.2.5). Constant C for the FeCr alloy is represented by C(ε · , T) = 0.0031 + 0.0273 sin ε · − 900 800 π sin T − 100 800 π Good agreement between experimental stresses and predicted stresses is obtained using the modified JC model that was introduced by Zhao et al. [85], with R: 0.87-0.99 and AARE: 7.6-22.7% when compared with those obtained using the original JC model with R: 0.83-0.88 and AARE: 21.8-44.8%. However, very accurate predictions were not found (cf. Figure 34).

Iturbe JC-Based Modification
Iturbe et al. [49] introduced a modification of the JC model to predict the flow behavior of Inconel 718 nickel superalloy under hot deformations. The Lurdos model [186] is used for the strain-hardening term, and the coupling effect between the temperature and strain rate is taken into account. The modified JC model that was introduced by Iturbe et al. [49] can be written as where is saturation stress, and 0 is yield stress. , , , , , , , and are material constants. At the reference strain rate and reference temperature, constants , , and are determined by fitting the experimental stress and strain data. At the reference temperature, constants , , and are determined using exponential fitting between the softening term in the original JC model with (cf. Figure 35a), and exponential fitting between the strain rate term and is used to compute constants and at the reference strain rate (cf. Figure 35b).

Iturbe JC-Based Modification
Iturbe et al. [49] introduced a modification of the JC model to predict the flow behavior of Inconel 718 nickel superalloy under hot deformations. The Lurdos model [186] is used for the strain-hardening term, and the coupling effect between the temperature and strain rate is taken into account. The modified JC model that was introduced by Iturbe et al. [49] can be written as where σ s is saturation stress, and σ 0 is yield stress. A, n, r, m, B, C, D, and x are material constants. At the reference strain rate and reference temperature, constants A, n, and r are determined by fitting the experimental stress and strain data. At the reference temperature, constants C, D, and x are determined using exponential fitting between the softening term in the original JC model with T (cf. Figure 35a), and exponential fitting between the strain rate term and T is used to compute constants m and B at the reference strain rate (cf. Figure 35b).

Iturbe JC-Based Modification
Iturbe et al. [49] introduced a modification of the JC model to predict the flow behavior of Inconel 718 nickel superalloy under hot deformations. The Lurdos model [186] is used for the strain-hardening term, and the coupling effect between the temperature and strain rate is taken into account. The modified JC model that was introduced by Iturbe et al. [49] can be written as where is saturation stress, and 0 is yield stress. , , , , , , , and are material constants. At the reference strain rate and reference temperature, constants , , and are determined by fitting the experimental stress and strain data. At the reference temperature, constants , , and are determined using exponential fitting between the softening term in the original JC model with (cf. Figure 35a), and exponential fitting between the strain rate term and is used to compute constants and at the reference strain rate (cf. Figure 35b).  The modified JC model that was presented by Iturbe et al. [49] provided good predictions of the flow behavior of the tested alloy compared with the predicted stresses obtained using the original JC model for the tested alloy. Still, precise predictions were not achieved (cf. Figure 36). The modified JC model that was presented by Iturbe et al. [49] provided good predictions of the flow behavior of the tested alloy compared with the predicted stresses obtained using the original JC model for the tested alloy. Still, precise predictions were not achieved (cf. Figure 36).

Tao JC-Based Modification
Tao et al. [63] presented a modification of the original JC model for the prediction of the flow behavior of a Ti−6Al−4V tube during warm bending at different strain rate and temperature ranges. In their modification, the strain rate term of the original JC model is replaced by a strain rate term that was previously introduced by Kang et al. [69], while the softening term is modified considering constant in the original JC model as a function of temperature. The modified JC that was presented by Tao et al. [63] can be expressed as = ( + )(1 + 1 ln( • * ) + 2 (ln( • * )) 2 ) (1 − * + * + * 2 + * 3 + * 4 ) where , , , 1 , 2 , , , , , and are material constants. Constants , , and are determined as explained for the original JC model in Section 3.1. Constants 1 and 2 are determined using quadratic fitting by plotting /( + ) − 1 vs. ln( • * ) at the reference temperature (cf. Figure 37a). Thus, the parameters , , , , and are determined to have different values at different strain values when plotted for ln[1 − (+ )] vs. ln ( * ) (cf. Figure 37b), and, finally, they are expressed as fourth-order polynomials fitted by plotting the obtained values vs. strain (cf. Figure 37c).

Tao JC-Based Modification
Tao et al. [63] presented a modification of the original JC model for the prediction of the flow behavior of a Ti−6Al−4V tube during warm bending at different strain rate and temperature ranges. In their modification, the strain rate term of the original JC model is replaced by a strain rate term that was previously introduced by Kang et al. [69], while the softening term is modified considering constant m in the original JC model as a function of temperature. The modified JC that was presented by Tao et al. [63] can be expressed as σ = (A + Bε n ) 1 + C 1 ln(ε · * ) + C 2 (ln(ε · * )) 2 1 − T * a+bT * +cT * 2 +dT * 3 +eT * 4 (50) where A, B, n, C 1 , C 2 , a, b, c, d, and e are material constants. Constants A, B, and n are determined as explained for the original JC model in Section 3.1. Constants C 1 and C 2 are determined using quadratic fitting by plotting σ/(A + Bε n ) − 1 vs. ln(ε · * ) at the reference temperature (cf. Figure 37a). Thus, the parameters a, b, c, d, and d are determined to have different values at different strain values when plotted for ln[1 − (+Bε n )] vs. ln (T * ) (cf. Figure 37b), and, finally, they are expressed as fourth-order polynomials fitted by plotting the obtained values vs. strain (cf. Figure 37c). The modified JC model that was presented by Iturbe et al. [49] provided good predictions of the flow behavior of the tested alloy compared with the predicted stresses obtained using the original JC model for the tested alloy. Still, precise predictions were not achieved (cf. Figure 36).

Tao JC-Based Modification
Tao et al. [63] presented a modification of the original JC model for the prediction of the flow behavior of a Ti−6Al−4V tube during warm bending at different strain rate and temperature ranges. In their modification, the strain rate term of the original JC model is replaced by a strain rate term that was previously introduced by Kang et al. [69], while the softening term is modified considering constant in the original JC model as a function of temperature. The modified JC that was presented by Tao et al. [63] can be expressed as = ( + )(1 + 1 ln( • * ) + 2 (ln( • * )) 2 ) (1 − * + * + * 2 + * 3 + * 4 ) where , , , 1 , 2 , , , , , and are material constants. Constants , , and are determined as explained for the original JC model in Section 3.1. Constants 1 and 2 are determined using quadratic fitting by plotting /( + ) − 1 vs. ln( • * ) at the reference temperature (cf. Figure 37a). Thus, the parameters , , , , and are determined to have different values at different strain values when plotted for ln[1 − (+ )] vs. ln ( * ) (cf. Figure 37b), and, finally, they are expressed as fourth-order polynomials fitted by plotting the obtained values vs. strain (cf. Figure 37c).  Good agreement can be seen between the experimental stresses and predicted stresses for the prediction of the flow behavior of the Ti-6Al-4V alloy obtained using the modified JC model that was introduced by Tao et al. [63] (cf. Figure 38a). The modified JC model is verified by using new experimental data that were not included in the determination of the constants. Very good predictions are also obtained, which enhances the use of the coupling effect of strain and temperature in this modification (cf. Figure 38b). Good agreement can be seen between the experimental stresses and predicted stresses for the prediction of the flow behavior of the Ti−6Al−4V alloy obtained using the modified JC model that was introduced by Tao et al. [63] (cf. Figure 38a). The modified JC model is verified by using new experimental data that were not included in the determination of the constants. Very good predictions are also obtained, which enhances the use of the coupling effect of strain and temperature in this modification (cf. Figure  38b).

He JC-Based Modification
He et al. [86] introduced a modification of the original JC model to predict the flow behavior of a 10%Cr steel alloy at elevated temperatures and different strain rates. In their modification, strain-hardening terms are modified, and the coupling effect between the strain and strain rate is presented. The modified JC model that was introduced by He et al. [86] can be written as where 1 , 1 , 1 , 2 , 3 , 1 , and 2 are material constants. At the reference strain rate and reference temperature, Equation (51) reduces to = 1 1 ; by taking logarithms for both sides and plotting ln vs. ln , constant 1 can be determined from the intercept, while constant 1 can be determined from the slope. At the reference temperature, Equation (51) lowers to: To obtain constants 1 , 2 , and 3 , a new parameter, , is introduced, in which = 1 + 2 + 3 2 ; then, after performing some rearrangements, Equation (51) can be written as / 1 = 1 + 1 ln( • * ); after that, by plotting / 1 vs. ln( • * ), different values can be obtained from the slope ( 1 ) at different strain values (cf. Figure 39a). Finally, the obtained values of can be plotted vs. its corresponding values, and quadratic fitting can be implemented to obtain the values of 1 , 2 , and 3 (cf. Figure 39b). Constants 1 and 2 can be determined in the same way that was explained in Section 3.2.5. Figure 38. Experimental stresses (markers) compared to predicted stresses (solid lines) using the modified JC model that was presented by Tao et al. [63] for Ti-6Al-4V alloy for (a) the tested range of temperature and strain rates and (b) new experiments with different strain rates and temperatures for verification.

He JC-Based Modification
He et al. [86] introduced a modification of the original JC model to predict the flow behavior of a 10%Cr steel alloy at elevated temperatures and different strain rates. In their modification, strain-hardening terms are modified, and the coupling effect between the strain and strain rate is presented. The modified JC model that was introduced by He et al. [86] can be written as where A 1 , n 1 , b 1 , b 2 , b 3 , λ 1 , and λ 2 are material constants. At the reference strain rate and reference temperature, Equation (51) reduces to σ = A 1 ε n 1 ; by taking logarithms for both sides and plotting ln σ vs. ln ε, constant A 1 can be determined from the intercept, while constant n 1 can be determined from the slope. At the reference temperature, Equation (51) lowers to: To obtain constants b 1 , b 2 , and b 3 , a new parameter, D, is introduced, in which D = b 1 + b 2 ε + b 3 ε 2 ; then, after performing some rearrangements, Equation (51) can be written as σ/ε n 1 = A 1 + A 1 D ln(ε · * ); after that, by plotting σ/ε n 1 vs. ln(ε · * ), different values can be obtained from the slope (A 1 D) at different strain values (cf. Figure 39a). Finally, the obtained values of D can be plotted vs. its corresponding ε values, and quadratic fitting can be implemented to obtain the values of b 1 , b 2 , and b 3 (cf. Figure 39b). Constants λ 1 and λ 2 can be determined in the same way that was explained in Section 3.2.5.
Comparisons between the experimental stresses and predicted stresses for the 10%Cr steel alloy obtained using both the JC model and the modified JC model that was introduced by He et al. [86] are shown in Figure 40a and 40b, respectively. The modified JC model achieved a better agreement with the experimental stresses than the original JC model; however, not the predictions were not very precise. Comparisons between the experimental stresses and predicted stresses for the 10%Cr steel alloy obtained using both the JC model and the modified JC model that was introduced by He et al. [86] are shown in Figure 40a and 40b, respectively. The modified JC model achieved a better agreement with the experimental stresses than the original JC model; however, not the predictions were not very precise.

Hou X. JC-Based Modification
Hou X. et al. [87] modified the original JC model for the prediction of the flow behavior of Ti-6Al-4V at elevated temperatures and different strain rates. The strainhardening term is correlated with temperature in their modification, while both the strain rate and softening terms are kept as in the original JC model. The modified JC model that was presented by Hou X. et al. [87] can be written as = ( + (1 + 1 ln ( ⁄ )) ) (1 + ln • * )(1 − * 2 ) where , , 1 , , , and 2 are material constants. At the reference strain rate, after introducing a new parameter, = (1 + 1 ln ( ⁄ )), considering is the yield stress  Comparisons between the experimental stresses and predicted stresses for the 10%Cr steel alloy obtained using both the JC model and the modified JC model that was introduced by He et al. [86] are shown in Figure 40a and 40b, respectively. The modified JC model achieved a better agreement with the experimental stresses than the original JC model; however, not the predictions were not very precise.

Hou X. JC-Based Modification
Hou X. et al. [87] modified the original JC model for the prediction of the flow behavior of Ti-6Al-4V at elevated temperatures and different strain rates. The strainhardening term is correlated with temperature in their modification, while both the strain rate and softening terms are kept as in the original JC model. The modified JC model that was presented by Hou X. et al. [87] can be written as where , , 1 , , , and 2 are material constants. At the reference strain rate, after introducing a new parameter, = (1 + 1 ln ( ⁄ )), considering is the yield stress

Hou X. JC-Based Modification
Hou X. et al. [87] modified the original JC model for the prediction of the flow behavior of Ti-6Al-4V at elevated temperatures and different strain rates. The strain-hardening term is correlated with temperature in their modification, while both the strain rate and softening terms are kept as in the original JC model. The modified JC model that was presented by Hou X. et al. [87] can be written as where A, B, m 1 , n, C, and m 2 are material constants. At the reference strain rate, after introducing a new parameter, Q = B(1 + m 1 ln( T /T r )), considering A is the yield stress and m 2 = m, which are determined based on the original JC model, Equation (53) reduces to After performing some rearrangements, Equation (54) can be written as Taking logarithms for both sides and plotting ln[σ/(1 − T * m ) − A] vs. ln ε, constant n is the slope, and Q can be determined from the intercept. Different values of Q are obtained at different values of T; then, constants B and m 1 can be determined by fitting the obtained data of Q with T/T r . Constants C and m 2 can be determined in the same way that was explained for the original JC model (see Section 3.1).
A comparison between the experimental stresses and predicted stresses for the Ti-6Al-4V alloy under hot deformation using both the original JC model (cf. Figure 41a) and the modified JC model that was presented by Hou X. et al. [87] (cf. Figure 41b) is shown in Figure 41. As can be seen, the modified JC model gave a more accurate prediction (RE = 4.19%) than the original JC model (RE = 10.43%), which might be a result for taking the interaction of strain and temperature into account. and 2 = , which are determined based on the original JC model, Equation (53) reduces to After performing some rearrangements, Equation (54) can be written as Taking logarithms for both sides and plotting ln[ /(1 − * ) − ] vs. ln , constant is the slope, and can be determined from the intercept. Different values of are obtained at different values of ; then, constants and 1 can be determined by fitting the obtained data of with / . Constants and 2 can be determined in the same way that was explained for the original JC model (see Section 3.1).
A comparison between the experimental stresses and predicted stresses for the Ti-6Al-4V alloy under hot deformation using both the original JC model (cf. Figure 41a) and the modified JC model that was presented by Hou X. et al. [87] (cf. Figure 41b) is shown in Figure 41. As can be seen, the modified JC model gave a more accurate prediction (RE = 4.19%) than the original JC model (RE = 10.43%), which might be a result for taking the interaction of strain and temperature into account. Promoppatum et al. [187] studied the effect of constitutive models for the prediction of the mechanical behavior of the laser powder bed fusion of Ti-6Al-4V alloy using finite element simulations and the modified JC model that was introduced by Hou X. et al. [87], along with other models. The results showed that temperature-dependent strain hardening has a minimal effect on the predicted stress and strain fields. On the other hand, the rate-dependent term has a large effect on them. In contrary findings, Kugalur-Palanisamy et al. [188] reported the significance of the temperature-dependent strain hardening consequence on the deformation behavior of Ti-6Al-4V alloy through the machining process when using the modified JC model that was introduced by Hou X. et al. [87], along with others, in finite element simulations of the orthogonal cutting process.

Zhang JC-Based Modification
Zhang et al. [89] modified the softening term in the original JC to predict the flow behavior of AZ31 magnesium alloy at high temperatures and high strain rates. The modification is based on the dependence of the softening parameter on temperature, which can be expressed as Promoppatum et al. [187] studied the effect of constitutive models for the prediction of the mechanical behavior of the laser powder bed fusion of Ti-6Al-4V alloy using finite element simulations and the modified JC model that was introduced by Hou X. et al. [87], along with other models. The results showed that temperature-dependent strain hardening has a minimal effect on the predicted stress and strain fields. On the other hand, the rate-dependent term has a large effect on them. In contrary findings, Kugalur-Palanisamy et al. [188] reported the significance of the temperature-dependent strain hardening consequence on the deformation behavior of Ti-6Al-4V alloy through the machining process when using the modified JC model that was introduced by Hou X. et al. [87], along with others, in finite element simulations of the orthogonal cutting process.

Zhang JC-Based Modification
Zhang et al. [89] modified the softening term in the original JC to predict the flow behavior of AZ31 magnesium alloy at high temperatures and high strain rates. The modification is based on the dependence of the softening parameter on temperature, which can be expressed as where A, B, n, and C are material constants that are determined in the same way as explained in Section 3.1. In their modification, m(T) is determined as e m 0 +m 1 T+m 2 T 2 . After performing some rearrangements, at different values of strain rate, Equation (56) can be written as 1 − σ (A + Bε n )(1 + C ln(ε · * )) = T * m(T) By taking logarithms for both sides, ln[1 − σ(A + Bε n )(1 + C ln(ε · * ))] can be plotted vs. ln T * , in which m(T) is the slope of the equation. At different temperature and strain rate values, different values for m(T) can be obtained. Constants m 0 , m 1 , and m 2 are obtained by quadratically fitting m(T) vs. T using those values (cf. Figure 42a). 1 − ( + )(1 + ln( • * )) = * ( ) By taking logarithms for both sides, ln[1 − ( + )(1 + ln( • * ))] can be plotted vs. ln * , in which ( ) is the slope of the equation. At different temperature and strain rate values, different values for ( ) can be obtained. Constants 0 , 1 , and 2 are obtained by quadratically fitting ( ) vs.
using those values (cf. Figure 42a). A comparison between experimental stresses and predicted stresses obtained using the modified JC model that was introduced by Zhang et al. [89] for the prediction of the flow behavior of AZ31 magnesium alloy is shown in Figure 42b. Good agreement is obtained, which might be attributed to the dependence of the softening parameter on temperature.

Niu JC-Based Modification
Niu et al. [90] introduced a modification of the original JC model for the accurate prediction of the flow behavior of A356 alloy at elevated temperatures and different strain rates. In their modification, the strain-hardening rate is modified to be interconnected with strain. Furthermore, the strain rate parameter is established to be interrelated with the strain and strain rate. In addition, the softening parameter is considered a function of the strain, strain rate, and temperature. The modified JC model that was presented by Niu et al. [90] can be written as where , , , and are material constants. Constant is a function of and • , and constant is a function of , • , and ( − ). Constants and can be expressed as = 0 + 1 + 2 • + 3 ( − ) + 4 3 + 5 •3 + 6 ( − ) 3 At the reference strain rate and reference temperature, Equation (58) lowers to A comparison between experimental stresses and predicted stresses obtained using the modified JC model that was introduced by Zhang et al. [89] for the prediction of the flow behavior of AZ31 magnesium alloy is shown in Figure 42b. Good agreement is obtained, which might be attributed to the dependence of the softening parameter on temperature.

Niu JC-Based Modification
Niu et al. [90] introduced a modification of the original JC model for the accurate prediction of the flow behavior of A356 alloy at elevated temperatures and different strain rates. In their modification, the strain-hardening rate is modified to be interconnected with strain. Furthermore, the strain rate parameter is established to be interrelated with the strain and strain rate. In addition, the softening parameter is considered a function of the strain, strain rate, and temperature. The modified JC model that was presented by Niu et al. [90] can be written as where A, B, C, and A are material constants. Constant E is a function of ε and ε · , and constant F is a function of ε, ε · , and (T − T r ). Constants E and F can be expressed as At the reference strain rate and reference temperature, Equation (58) lowers to By plotting σ vs. ε, constants A, B, C, and D can be obtained using regression analysis. After performing some rearrangements, at the reference temperature, Equation (58) Consequently, constants E 0 , E 1 , E 2 , E 3 , E 4 , and E 5 can be determined using regression analysis. At different values of strain rates, and taking logarithms after performing some rearrangements, Equation (58) can be written as Accordingly, constants F 0 , F 1 , F 2 , F 3 , F 4 , F 5 , F 6 , and F 7 can be determined using regression analysis.
Very accurate predictions of the flow stress for the tested alloy are obtained when comparing predicted stresses obtained using the modified JC model that was introduced by Niu et al. [90], with experimental stresses with AARE = 1.46% (cf. Figure 43b). On the other hand, the original JC model provided inaccurate predictions, with AARE = 26.31% (cf. Figure 43a). The coupling effect of strain, strain rate, and temperature might be a reason for the accurate predictions using the modified JC model.
By plotting vs. , constants , , , and can be obtained using regression analysis. After performing some rearrangements, at the reference temperature, Equation (58) reduces to Consequently, constants 0 , 1 , 2 , 3 , 4 , and 5 can be determined using regression analysis. At different values of strain rates, and taking logarithms after performing some rearrangements, Equation (58) can be written as Accordingly, constants 0 , 1 , 2 , 3 , 4 , 5 , 6 , and 7 can be determined using regression analysis.
Very accurate predictions of the flow stress for the tested alloy are obtained when comparing predicted stresses obtained using the modified JC model that was introduced by Niu et al. [90], with experimental stresses with AARE = 1.46% (cf. Figure 43b). On the other hand, the original JC model provided inaccurate predictions, with AARE = 26.31% (cf. Figure 43a). The coupling effect of strain, strain rate, and temperature might be a reason for the accurate predictions using the modified JC model.

Chakrabarty JC-Based Modification
Chakrabarty et al. [91] presented a modification of the original JC model to predict the flow behavior of copper at very high strain rates. The modified JC model that was presented by Chakrabarty et al. [91] can be expressed as where constants , , , and m are determined in the same way that was explained for the original JC model (see Section 3.1). At the reference temperature, after performing some rearrangements, Equation (64) can be expressed as Chakrabarty et al. [91] determined both constants, . and , using a fitting toolbox in MATLAB that is based on the non-linear least squares method (cf. Figure 44a). The

Chakrabarty JC-Based Modification
Chakrabarty et al. [91] presented a modification of the original JC model to predict the flow behavior of copper at very high strain rates. The modified JC model that was presented by Chakrabarty et al. [91] can be expressed as where constants A, B, n, and m are determined in the same way that was explained for the original JC model (see Section 3.1). At the reference temperature, after performing some rearrangements, Equation (64) can be expressed as Chakrabarty et al. [91] determined both constants, ε . c and D, using a fitting toolbox in MATLAB that is based on the non-linear least squares method (cf. Figure 44a). The modified JC model provided good predictions of the flow stress of copper when compared with the original JC model (cf. Figure 44b).  Figure 44b).

Li JC-Based Modification
Li et al. [92] modified the JC model to predict the flow behavior of Roma Plastilina No. 1 (RP # 1) clay, which is used as a backing material in ballistic tests, as where 1 , , 2 , , 1 , 2 , 1 , 2 , and 3 are material constants. At the reference strain rate and reference temperature, Equation (66) reduces to By plotting vs. , constants 1 , , 2 , and can be determined using curvefitting based on the non-linear least squares method incorporated in MATLAB (cf. Figure  45a). At the reference temperature, constants 1 and 2 can be determined using the same method that was explained in Section 3.2.21. Taking the logarithm, after performing some rearrangements, at different strain rate values, Equation (66) where (ln( • * )) is the slope of the equation. At different temperatures and different strain rates, different values of (ln( • * )) can be obtained, in which constants 1 , 2 , and 3 can be obtained by quadratically fitting (ln( • * )) vs. ln( • * ) using those values. The modified JC model that was introduced by Li et al. [92] succeeded in providing good predictions of the flow behavior of the tested alloy at different strain rates and different temperatures, with R values ranging from 0.9896 to 0.9995 and AARE values ranging from 15.74% to 28.93% (cf. Figure 45b).

Li JC-Based Modification
Li et al. [92] modified the JC model to predict the flow behavior of Roma Plastilina No. 1 (RP # 1) clay, which is used as a backing material in ballistic tests, as where A 1 , n, A 2 , K, C 1 , C 2 , m 1 , m 2 , and m 3 are material constants. At the reference strain rate and reference temperature, Equation (66) reduces to By plotting σ vs. ε, constants A 1 , n, A 2 , and K can be determined using curve-fitting based on the non-linear least squares method incorporated in MATLAB (cf. Figure 45a). At the reference temperature, constants C 1 and C 2 can be determined using the same method that was explained in Section 3.2.21. Taking the logarithm, after performing some rearrangements, at different strain rate values, Equation (66) can be written as where m(ln(ε · * )) is the slope of the equation. At different temperatures and different strain rates, different values of m(ln(ε · * )) can be obtained, in which constants m 1 , m 2 , and m 3 can be obtained by quadratically fitting m(ln(ε · * )) vs. ln(ε · * ) using those values. The modified JC model that was introduced by Li et al. [92] succeeded in providing good predictions of the flow behavior of the tested alloy at different strain rates and different temperatures, with R values ranging from 0.9896 to 0.9995 and AARE values ranging from 15.74% to 28.93% (cf. Figure 45b).

Qian JC-Based Modification
Qian et al. [93] presented a modification of the JC model for the prediction of a CuCrZr alloy at high strain rates and high temperatures. The coupling effect between the strain and both the strain rate and temperature is taken into consideration in this modification. The modified JC model that was presented by Qian et al. [93] can be expressed as where , , , 0 , 1 , and 2 are material constants. Qian et al. [93] implemented the strain rate parameter ( , ln( • * )), which was first presented by Abd El-Aty et al. [189] in their modification, which can be expressed as ( , ln( • * )) = 0 + 1 + 2 2 + 3 ln( • * ) + 4 ln( • * ) + 5 (ln( • * )) 2 Constants , , and are determined as explained in the original JC model (see Section 3.1). At the reference temperature, Equation (69) Constant 0 can be determined from the intercept of the equation, while 1 and 2 can be determined from the slope, as explained in Section 3.2.24. The modified JC model that was introduced by Qian et al. [93] improved the accuracy of the prediction of the flow behavior of the CuCrZr alloy, with AARE values ranging from 0.59% to 3.22% (cf. Figure 46b) compared with the predictions obtained using the original JC model, with AARE values range from 7.41% to 10.22% (cf. Figure 46a). Taking the interconnection between the strain, strain rate, and temperature might be a reason for the obtained accuracy of the predictions.

Qian JC-Based Modification
Qian et al. [93] presented a modification of the JC model for the prediction of a CuCrZr alloy at high strain rates and high temperatures. The coupling effect between the strain and both the strain rate and temperature is taken into consideration in this modification. The modified JC model that was presented by Qian et al. [93] can be expressed as where A, B, n, m 0 , m 1 , and m 2 are material constants. Qian et al. [93] implemented the strain rate parameter C(ε, ln(ε · * )), which was first presented by Abd El-Aty et al. [189] in their modification, which can be expressed as C(ε, ln(ε · * )) = C 0 + C 1 ε + C 2 ε 2 + C 3 ε ln(ε · * ) + C 4 ln(ε · * ) + C 5 (ln(ε · * )) 2 (70) Constants A, B, and n are determined as explained in the original JC model (see Section 3.1). At the reference temperature, Equation (69) lowers to σ = (A + Bε n )(1 + C(ε, ln(ε · * )) ln(ε · * )) (71) Constants C 0 , C 1 , C 2 , C 3 , C 4 , and C 5 can be determined by plotting σ/(A + Bε n ) − 1 vs. ln(ε · * ) using regression analysis. Taking the logarithm after performing some rearrangements and at different strain rate values, Equation (69) can be expressed as Constant m 0 can be determined from the intercept of the equation, while m 1 and m 2 can be determined from the slope, as explained in Section 3.2.24.
The modified JC model that was introduced by Qian et al. [93] improved the accuracy of the prediction of the flow behavior of the CuCrZr alloy, with AARE values ranging from 0.59% to 3.22% (cf. Figure 46b) compared with the predictions obtained using the original JC model, with AARE values range from 7.41% to 10.22% (cf. Figure 46a). Taking the interconnection between the strain, strain rate, and temperature might be a reason for the obtained accuracy of the predictions.

Liu JC-Based Modification
Liu et al. [94] modified both the strain-hardening and strain rate terms in the original JC model to predict the flow behavior of SWRH82B steel alloy at high temperatures and different strain rates. The modified JC that was presented by Liu et al. [94] can be expressed as where , , and are material constants. Parameter ( , ln( • * )) is determined as ( , ln( • * )) = exp( 0 + 1 + 2 2 + 3 ln( • * ) + 4 ln( • * ) + 5 (ln( • * )) 2 ) At the reference strain rate and reference temperature, Equation (73) lowers to By plotting vs. , constants and can be determined using regression analysis. Taking logarithms for both sides after performing some rearrangements and at the reference temperature, Equation (73) where constants 0 , 1 , 2 , 3 , 4 , and 5 can be determined using regression analysis. The modified JC model that was presented by Liu et al. [94] showed good predictability (cf. Figure 47b) when compared with the original JC model (cf. Figure 47a) for the prediction of the flow behavior of SWRH82B steel alloy with two different specimen configurations.

Liu JC-Based Modification
Liu et al. [94] modified both the strain-hardening and strain rate terms in the original JC model to predict the flow behavior of SWRH82B steel alloy at high temperatures and different strain rates. The modified JC that was presented by Liu et al. [94] can be expressed as where A, B, and m are material constants. Parameter C(ε, ln(ε · * )) is determined as C(ε, ln(ε · * )) = exp C 0 + C 1 ε + C 2 ε 2 + C 3 ε ln(ε · * ) + C 4 ln(ε · * ) + C 5 (ln(ε · * )) 2 (74) At the reference strain rate and reference temperature, Equation (73) lowers to By plotting σ vs. ε, constants A and B can be determined using regression analysis. Taking logarithms for both sides after performing some rearrangements and at the reference temperature, Equation (73) can be written as (76) where constants C 0 , C 1 , C 2 , C 3 , C 4 , and C 5 can be determined using regression analysis. The modified JC model that was presented by Liu et al. [94] showed good predictability (cf. Figure 47b) when compared with the original JC model (cf. Figure 47a) for the prediction of the flow behavior of SWRH82B steel alloy with two different specimen configurations.

Yu JC-Based Modification
Yu et al. [60] proposed a modification of the original JC model to predict the flow behavior of TA23 titanium alloy under hot deformation. In their modification, the softening term is modified to take the coupling effect between the strain rate and temperature into account. The modified JC model that was presented by Yu et al. [60] can be written as where , 1 , 2 , 1 , 1 , 2 , and 3 are material constants. Constants , 1 , 2 , and 1 are determined in the same way that was explained in Section 3.2.5. Taking the logarithm after performing some rearrangements, Equation (77) can be introduced as ln [ ( + 1 + 2 2 )(1 + 1 ln • * ) ] = ( 1 + 2 ln • * + 3 (ln( • * )) 2 )( − ) By introducing a new parameter, = 1 + 2 ln • * + 3 (ln( • * )) 2 , different values of (the slope of Equation (78)) can be determined at different strain rate and temperature values, which can be fitted as quadratic functions in ln • * ; then, constants 1 , 2 , and 3 can be determined from the fitting (cf. Figure 48a).
Compared with the original JC model, the modified JC model that was presented by Yu et al. [60] prediction a good prediction of the flow behavior of the TA23 titanium alloy under hot deformation, with an RE value of 3.28% compared with 5.32% for original JC model (cf. Figure 48b). The interrelation between the strain rate and temperature, which has been taken into account, might be a reason for the obtained agreement.

Yu JC-Based Modification
Yu et al. [60] proposed a modification of the original JC model to predict the flow behavior of TA23 titanium alloy under hot deformation. In their modification, the softening term is modified to take the coupling effect between the strain rate and temperature into account. The modified JC model that was presented by Yu et al. [60] can be written as where A, B 1 , B 2 , C 1 , λ 1 , λ 2 , and λ 3 are material constants. Constants A, B 1 , B 2 , and C 1 are determined in the same way that was explained in Section 3.2.5. Taking the logarithm after performing some rearrangements, Equation (77) can be introduced as By introducing a new parameter, λ = λ 1 + λ 2 ln ε · * + λ 3 (ln(ε · * )) 2 , different values of λ (the slope of Equation (78)) can be determined at different strain rate and temperature values, which can be fitted as quadratic functions in ln ε · * ; then, constants λ 1 , λ 2 , and λ 3 can be determined from the fitting (cf. Figure 48a).
Compared with the original JC model, the modified JC model that was presented by Yu et al. [60] prediction a good prediction of the flow behavior of the TA23 titanium alloy under hot deformation, with an RE value of 3.28% compared with 5.32% for original JC model (cf. Figure 48b). The interrelation between the strain rate and temperature, which has been taken into account, might be a reason for the obtained agreement.

Wang (4) JC-Based Modification
Wang et al. [95] modified the softening term in the original JC model to predict the flow behavior of metallic materials at a wide range of strain rates and temperatures. The dependency of temperature on the elastic modulus and Poisson's ratio is taken into consideration in the modified JC model. The modified JC model that was introduced by Wang et al. [95] can be expressed as where T 0 is defined as room temperature; parameters E T and E 0 are the elastic modulus at the tested temperature and room temperature, respectively; and µ T and µ T 0 are Poisson's ratio at the tested temperature and room temperature, respectively. Plots of the elastic modulus vs. temperature for the Ti-6Al-4V and tantalum alloys are shown in Figures 49a and 49b, respectively. The modification contains only four constants, A, B, n, and C; however, they must be determined at the same time. Otherwise, regression analysis will not help. Wang et al. [95] presented a multi-objective technique along with the Latin hypercube sampling method, Spearman rank correlation analysis, and a modern genetic algorithm to determine the four constants at once.

Wang (4) JC-Based Modification
Wang et al. [95] modified the softening term in the original JC model to predict the flow behavior of metallic materials at a wide range of strain rates and temperatures. The dependency of temperature on the elastic modulus and Poisson's ratio is taken into consideration in the modified JC model. The modified JC model that was introduced by Wang et al. [95] can be expressed as where 0 is defined as room temperature; parameters and 0 are the elastic modulus at the tested temperature and room temperature, respectively; and and 0 are Poisson's ratio at the tested temperature and room temperature, respectively. Plots of the elastic modulus vs. temperature for the Ti-6Al-4V and tantalum alloys are shown in Figure 49a and Figure 49b, respectively. The modification contains only four constants, , , , and ; however, they must be determined at the same time. Otherwise, regression analysis will not help. Wang et al. [95] presented a multi-objective technique along with the Latin hypercube sampling method, Spearman rank correlation analysis, and a modern genetic algorithm to determine the four constants at once.

Wang (4) JC-Based Modification
Wang et al. [95] modified the softening term in the original JC model to predict the flow behavior of metallic materials at a wide range of strain rates and temperatures. The dependency of temperature on the elastic modulus and Poisson's ratio is taken into consideration in the modified JC model. The modified JC model that was introduced by Wang et al. [95] can be expressed as where 0 is defined as room temperature; parameters and 0 are the elastic modulus at the tested temperature and room temperature, respectively; and and 0 are Poisson's ratio at the tested temperature and room temperature, respectively. Plots of the elastic modulus vs. temperature for the Ti-6Al-4V and tantalum alloys are shown in Figure 49a and Figure 49b, respectively. The modification contains only four constants, , , , and ; however, they must be determined at the same time. Otherwise, regression analysis will not help. Wang et al. [95] presented a multi-objective technique along with the Latin hypercube sampling method, Spearman rank correlation analysis, and a modern genetic algorithm to determine the four constants at once.  Comparisons between the experimental stresses and predicted stresses obtained by the modified JC model that was introduced by Wang et al. [95] for the prediction of the flow behavior of Ti-6Al-4V and tantalum are shown in Figures 50a and 50b, respectively. As can be seen, the modified model fits the experimental data for Ti-6Al-4V and tantalum very well.
Comparisons between the experimental stresses and predicted stresses obtained by the modified JC model that was introduced by Wang et al. [95] for the prediction of the flow behavior of Ti-6Al-4V and tantalum are shown in Figure 50a and Figure 50b, respectively. As can be seen, the modified model fits the experimental data for Ti-6Al-4V and tantalum very well.

Shokry (2) JC-Based Modification
Shokry et al. [96] introduced an improved generic modification of the original JC model for flow behavior predictions for different element-based alloys at elevated temperatures and a wide range of strain rates. The coupling effect between the strain and strain rate, as well as between the strain, strain rate, and temperature, is taken into consideration in this modification. The improved generic modification of the original JC model that was introduced by Shokry et al. [96] can be written as The four constants of constitute the strain-hardening term, and the nine constants of constitute the strain rate term, while the twenty-seventh constant of constitutes the softening term. All forty constants are determined using regression analysis in MATLAB. At the reference strain rate and reference temperature, Equation (80) reduces to After expansion, Equation (81) extends to four terms with four constants that are constituted with strain and are determined by utilizing regression analysis. At the reference temperature, after performing some rearrangements, Equation (80)

Shokry (2) JC-Based Modification
Shokry et al. [96] introduced an improved generic modification of the original JC model for flow behavior predictions for different element-based alloys at elevated temperatures and a wide range of strain rates. The coupling effect between the strain and strain rate, as well as between the strain, strain rate, and temperature, is taken into consideration in this modification. The improved generic modification of the original JC model that was introduced by Shokry et al. [96] can be written as The four constants of A i constitute the strain-hardening term, and the nine constants of C ij constitute the strain rate term, while the twenty-seventh constant of m ijk constitutes the softening term. All forty constants are determined using regression analysis in MATLAB. At the reference strain rate and reference temperature, Equation (80) reduces to After expansion, Equation (81) extends to four terms with four constants that are constituted with strain and are determined by utilizing regression analysis. At the reference temperature, after performing some rearrangements, Equation (80) is simplified to Nine constants are constituted with the strain and strain rate and are determined utilizing regression analysis. At different strain rate values, after performing some rearrangements, Equation (80) is expressed as The right term provides 27 constants of m ijk after an expansion that is constituted by strain, strain rate, and temperature and is determined using regression analysis.
Precise predictions of the flow behavior of nickel-based (U720LI) and aluminumbased (AA7020) alloys using the improved generic modification of the original JC that was presented by Shokry et al. [96] are achieved (cf. Figure 51c1,c2) compared with experimental stresses, as well as predicted stresses obtained using both original JC model (cf. Figure 51a1,a2) and the modified JC model that was presented by Lin et al. [73] (cf. Figure 51b1,b2). Among the employed models, the improved generic modification provides the highest R, with values of 0.994 ± 0.013, and the lowest AARE, with values of 1.95 ± 1.08%, for the six tested alloys. The very accurate predictions of the flow behavior are due to the coupling effect between the strain, strain rate, and temperature.
Nine constants are constituted with the strain and strain rate and are determined utilizing regression analysis. At different strain rate values, after performing some rearrangements, Equation (80) The right term provides 27 constants of after an expansion that is constituted by strain, strain rate, and temperature and is determined using regression analysis.
Precise predictions of the flow behavior of nickel-based (U720LI) and aluminumbased (AA7020) alloys using the improved generic modification of the original JC that was presented by Shokry et al. [96] are achieved (cf. Figure 51c1,c2) compared with experimental stresses, as well as predicted stresses obtained using both original JC model (cf. Figure 51a1,a2) and the modified JC model that was presented by Lin et al. [73] (cf. Figure 51b1,b2). Among the employed models, the improved generic modification provides the highest R, with values of 0.994 ± 0.013, and the lowest AARE, with values of 1.95 ± 1.08%, for the six tested alloys. The very accurate predictions of the flow behavior are due to the coupling effect between the strain, strain rate, and temperature. Figure 51. Comparison between experimental stresses (markers) and predicted stresses (solid lines) for U720LI (first raw) and AA7020 (second raw) alloys obtained using (a1,a2) the original JC model, (b1,b2) the modified JC model that was presented by Lin et al. [73], and (c1,c2) the improved generic modified JC model that was presented by Shokry et al. [96].

Priest JC-Based Modification
Priest et al. [97] introduced a modification of the original JC model for the prediction of machining simulations in C45 steel at different temperatures and different strain rates. Both strain rate and softening terms are modified in the presented modification to take the coupling effect between the temperature and strain rate into account. The modified JC that was introduced by Priest et al. [97] can be written as ) (84) Figure 51. Comparison between experimental stresses (markers) and predicted stresses (solid lines) for U720LI (first raw) and AA7020 (second raw) alloys obtained using (a1,a2) the original JC model, (b1,b2) the modified JC model that was presented by Lin et al. [73], and (c1,c2) the improved generic modified JC model that was presented by Shokry et al. [96].

Priest JC-Based Modification
Priest et al. [97] introduced a modification of the original JC model for the prediction of machining simulations in C45 steel at different temperatures and different strain rates. Both strain rate and softening terms are modified in the presented modification to take the coupling effect between the temperature and strain rate into account. The modified JC that was introduced by Priest et al. [97] can be written as where A, B, n, a c , b c , c c , a m , b m , and c m are material constants. Constants A, B, and n are determined in the same way that was explained with the original JC model (see Section 3.1). At the reference temperature, after performing some rearrangements, Equation (84) can be written in the form of Using linear regression, at different temperature and strain rate values, different values for the slope of Equation (85), a c exp(b c T) + c c , are obtained, which can be exponentially fitted with temperature to obtain constants a c , b c , and c c (cf. Figure 52a). At the reference strain rate, and by taking logarithm after performing some rearrangements, Equation (84) can be written as Using linear regression, at different temperature and strain rate values, different values for the slope of Equation (85), exp( ) + , are obtained, which can be exponentially fitted with temperature to obtain constants , , and (cf. Figure 52a). At the reference strain rate, and by taking logarithm after performing some rearrangements, Equation (84) can be written as At different temperatures, different values for the slope of Equation (86), exp (− (( − )/ ) 2 ), are obtained, which can be fitted using a gaussian fitting with temperature to obtain constants , , and (cf. Figure 52b). A comparison between the experimental stresses and predicted stresses for the flow behavior of C45 steel alloy under hot deformations obtained using both the original JC model and the modified JC model that was introduced by Priest et al. [97] is shown in Figure 53. The modified model (cf. Figure 53b) provided better predictions than those obtained using the original JC model (cf. Figure 53a).  Figure 52b).
A comparison between the experimental stresses and predicted stresses for the flow behavior of C45 steel alloy under hot deformations obtained using both the original JC model and the modified JC model that was introduced by Priest et al. [97] is shown in Figure 53. The modified model (cf. Figure 53b) provided better predictions than those obtained using the original JC model (cf. Figure 53a).

Discussions and Summary
The flow behavior of metals and alloys is highly affected by hot working conditions, strain, strain rate, and temperature. As temperature decreases and strain rate increases, strain hardening increases. This is due to the emergence of crystal defects, mainly dislocations, strain-induced stages, or twin boundaries through plastic deformation [190][191][192]. Both the strain strength coefficient and strain-hardening index rise very quickly in the initial strain, followed by a decrease and steady value at high strain rates [14,193]. In general, an increasing strain rate is followed by an increase in flow stress; this is due to the restricted time that is required for dynamic recovery and growth associated with the

Discussions and Summary
The flow behavior of metals and alloys is highly affected by hot working conditions, strain, strain rate, and temperature. As temperature decreases and strain rate increases, strain hardening increases. This is due to the emergence of crystal defects, mainly dislocations, strain-induced stages, or twin boundaries through plastic deformation [190][191][192].
Both the strain strength coefficient and strain-hardening index rise very quickly in the initial strain, followed by a decrease and steady value at high strain rates [14,193]. In general, an increasing strain rate is followed by an increase in flow stress; this is due to the restricted time that is required for dynamic recovery and growth associated with the nucleation of dynamic recrystallization. In contrast, plenty of time leads to a slow rate of dynamic recovery as well as dynamic recrystallization at low strain rates and high temperatures [84,134].
Modeling the flow behavior of metals and alloys at a wide range of temperatures and strain rates is essential to optimize hot working conditions as well as predict the mechanical behavior of metals and alloys with applications under severe conditions such as dynamic loadings and elevated temperatures. The Johnson-Cook model [47] is one of the most important phenomenological models and is widely used to predict the flow behavior of materials at different strain rates and different temperatures. The model contains three independent terms: (i) strain hardening, (ii) strain rate, and (iii) thermal softening. One of the most powerful sides of the JC model is that it has only five constants, which constitute the relationship between flow stress and strain. However, a lot of criticism has been reported vs. the accuracy of predictions of flow behavior using the JC model, particularly with complex non-linear flow stress behavior. As a matter of fact, the JC model does not take the coupling effect of the strain, strain rate, and temperature into consideration in the variation of material parameters, which might be an acceptable reason for the inaccurate predictions of the JC model. Consequently, many modifications of the JC model have been introduced to precisely predict the flow behavior of metals and alloys.
The JC model starts by predicting flow stress using the strain-hardening term, followed by multiplying strain rate and softening terms. The JC model employs the Ludwik equation [194] to represent the strain-hardening term, σ = A + Bε n , in which A is yield stress, and B and n are the strain-hardening strength and strain-hardening exponent, in which the strain-hardening exponent considers the balance between softening and work hardening. The JC model achieves good agreement for strain hardening between the experimental and predicted stresses using the Ludwik equation for different alloys, such as nickel [48], iron [30], aluminum [51,52], and titanium-based [60] alloys. In order to describe the saturation of strain hardening, as well as the continuance of strain hardening with the progress in strain, Shin and Kim [75], replaced the Ludwik equation with the Voce hardening [180] equation in a modified JC model, which is represented by σ = [A + B{1 − exp(−Cε)}], in which A represents yield stress, constants B and C introduce saturated strain hardening and the strain-hardening index respectively. Good agreements between the experimental and predicted stresses for copper can be obtained. Similar findings were obtained using Voce hardening as a modification term in the original JC model that was presented by Wang et al. [77] for 30Cr2Ni4MoV rotor steel alloy. To overcome the difficulty of finding the real value of the yield stress, the yield stress constant, A, is replaced by beak stress, σ 0 , using the Zener-Hollomon parameter, which is provided by σ 0 = (1/α) ln (Z/A) 1/m + (Z/A) 2/m + 1 1/2 , where Z, α, A, and m are defined in Section 3.2.10. This modification of the strain-hardening term is introduced with good predictions of the flow stress of Al-Zn-Mg-Cu alloy [78] and 7075 Al alloy [79]. A different replacement for the strain-hardening term of the original JC model was presented by Lin et al. [73], in which the stress-strain data at both the reference strain rate and reference temperature are fitted using quadratic function as A + B 1 ε + B 2 ε 2 with material constants A, B 1 , and B 2 . A similar replacement with a quadratic polynomial was implemented as in [60,81] for TA23 titanium alloy and titanium matrix composite, and with polynomial fitting using a three-order function as in [84,85,96] for alloy 800H, FeCr alloy, and different element-based alloys.
The strain rate term has received a lot of modifications in the original JC model, which can be implemented when studying hot deformation as well as dynamic loadings. Rule and Jones [68] modified the strain rate term in the original JC model to mimic the quick rise in stress due to dynamic loadings, with good predictions of flow stress for different alloys such as 7075-T6 aluminum and high-strength steel alloys, in which the term C 4 (1/(C 5 − ln ε · * ) − 1/C 5 ) is added to the strain rate term in the original JC model, where C 4 and C 5 are material constants. Another modification of the prediction of the dynamic behavior of a vehicle body under a very high strain rate was presented by Kang et al. [69], in which the strain rate term is expressed as a quadratic function of the logarithm of the strain rate by adding the term C 2 (ln ε · * ) 2 to the strain rate term in the original JC model, and C 2 is a material constant, which gave provided predictions for steel plate material. Good predictions can also be obtained when using the same strain rate modification for Ti-6Al-4V tubes [63] and ballistic Roma Plastilina No. 1 clay [92]. Another modification was introduced by considering parameter C in the strain rate term of the original JC model depending on the strain, with a linear relationship [77,84] for the predictions of the flow behavior of 30Cr2Ni4MoV rotor steel and alloy 800H and a quadratic relationship [86] to predict the flow behavior of 10% Cr steel alloy. More interpretations of the strain rate term as quadratic functions in both the strain and strain rate are presented in [90,93,94,96], with good agreements between the experimental and predicted stresses. Couque et al. [71] obtained good predictions of the flow behavior of nickel when an exponential term is added to the strain rate in the original JC model as E(ε · /ε · 1 ) k , in which E and k are material constants. A sine wave representation of the strain rate term is introduced to modify the original JC model in [48,85,88] for the prediction of the flow behavior of Inconel 718, FeCr, and SnSbCu alloys, accompanied by good predictions.
The thermal-softening term in the original JC model has also received a lot of modifications. The early modified softening term was introduced by Meyers et al. [67], in which the softening term in the original JC model is replaced by an exponential term e −λ(T−T r ) , and λ is a material constant. Lin et al. [73] introduced a modification of the softening term that takes the coupling effect between the strain rate and temperature into account, considering the parameter λ = λ 1 + λ 2 ln ε · * , in which λ 1 and λ 2 are two constants that constitute the linear relationship between λ and ln ε · * . This modification has been widely used, with good predictions of flow behavior for different alloys, such as 28CrMnMoV steel [170], Hastelloy C-276 alloy [171], 40CrNi steel [172], and 2297 Al-Li alloy [174]. It is also employed in other modifications of the original JC model, as in [77], with 30Cr2Ni4MoV rotor steel alloy and [86] 10%Cr steel alloy. The same modification, but with a quadratic relationship instead of a linear one, provided λ = λ 1 + λ 2 ln ε · * + λ 3 (ln(ε · * )) 2 for the prediction of the flow behavior of TA23 titanium alloy, which was utilized in [60]. Another modification of the softening term was introduced by Hou Q. Y. et al. [74] to overcome the problem of predicting flow behavior at a lower temperature than the reference temperature. In this modification, T * in the original JC model is replaced by λ e T/T m − e T r /T m / e − e T r /T m , which provides good predictions for Mg-10Gd-2Y-0.5Zr alloy. The parameter λ is replaced by λ ε in another modification that was presented by Perez et al. [178] so that λ has a constant value considering the change in strain. Lin et al. [78] replaced constant m in the original JC model with a linear function in the strain rate as 0.71291 + 0.04391ε · , with good agreement with the prediction of the flow behavior of a Zn-Mg-Cu alloy with experimental stresses. In the same way, constant m is replaced by a linear relationship but in the strain as m 1+ m 2 ε for predictions of flow behavior of alloy 800H [84] and CuCrZr alloy [93], providing good predictions. By using the fourth polynomial equation in temperature, the m constant is modified as a + bT * + cT * 2 + dT * 3 + eT * 4 in another modification that was presented by Tao et al. [63] for the prediction of the flow behavior of Ti−6Al−4V, in which constants a, b, c, d, and e are determined using polynomial fitting. The strain rate is coupled with the quadratic function of temperature in another modification for the softening term in the original JC model, which was introduced by Song et al. [81] for the prediction of the flow behavior of a titanium matrix composite and expressed as exp λ 1 T * + λ 2 T * 2 ln ε · * , providing good predictions. Interactions between the effects of the strain, strain rate, and temperature are taken into consideration in other modifications, in which constant λ, which was introduced in [67], is replaced by a third-order multiplier of the strain, strain rate, and temperature, with eight selected terms in [90] for the prediction of A356 alloy and with twenty-seven terms in [96] for the prediction of four different element-based alloys, providing precise predictions.

Conclusions and Future Directions
In this review article, the flow behavior of metals and alloys at a wide range of strain rates and different temperatures is studied. Thus, the constitutive model of the original JC model, as well as more than thirty modified JC-based models, are critically reviewed and commented on. In addition, the methods and techniques that are used to determine model constants are presented and explained. Finally, a summary of modifications based on the three terms of the original JC model is presented.
The combination of multiple factors and their interactions, which may affect the flow stress response, is a much more efficient research methodology for empirical modeling. The Johnson-Cook model is a strong phenomenological model that has been extensively used for predictions of the flow behavior of metals and alloys. It has been implemented in finite element software packages, which enhances its importance in performing simulation analysis of hot working processes; optimizing strain, strain rate, and temperature; and simulating real applications in severe conditions. The findings of this review can be summarized as follows:

•
The Johnson-Cook model has been widely used to predict the flow behavior of metals and alloys at a wide range of temperatures and strain rates. In this regard, modified JC-based models were introduced for accurate predictions of the flow behavior for metals and alloys with complex non-linear behavior in which the JC model fails to precisely predict the flow behavior. • Lin et al. [73] provided one of the most important modifications of the original JC model, which has been used for the prediction of the flow behavior of different metals and alloys. The modification provides good predictability of flow behavior when compared with experimental stresses in many different alloys.

•
The improved generic modification for the JC model that was introduced by Shokry et al. [96] can be considered one of the most promising modifications, in which the interaction between the strain, strain rate, and temperature are taken into account. Accurate predictions are obtained using the improved generic JC model when implemented in six different element-based alloys. One of the limitations of this modification is the large number of constants it has: forty constants. However, this problem is not a concern nowadays due to the rapid improvements in computers and the software used for the determination of the constants. • Comparing predicted stresses using the original JC model and the modified JC-based models for the same types of alloys might be a future direction to precisely assess and evaluate the predictability of the JC model and modified JC-based models. • Another future direction is considering the inverse analysis methods and the techniques that are based on non-linear least squares methods to minimize the mean square error between experimental and predicted values for the accurate determination of model constants.

•
Coupling the original JC model and the modified JC-based models with other models such as the Zerilli-Armstrong and Arrhenius models might be another future direction.