Constitutive Modeling of Annealed OFHC with Wide Strain-Rate and Temperature Effects: Incorporating Dislocation Dynamics and Normalized Microstructural Size Evolution

The flow stress of face-centered cubic (FCC) metals exhibits a rapid increase near a strain rate of 104 s−1 under fixed-strain conditions. However, many existing constitutive models either fail to capture the mechanical characteristics of this plastic deformation or use piecewise strain-rate hardening models to describe this phenomenon. Unfortunately, these piecewise models may suffer from issues such as discontinuity of physical quantities and difficulties in determining segment markers, and struggle to reflect the underlying physical mechanisms that give rise to this mutation phenomenon. In light of this, this paper proposes that the abrupt change in flow stress sensitivity to strain rate in FCC metals can be attributed to microstructural evolution characteristics. To address this, a continuous semiempirical physical constitutive model for FCC metals is established based on the microstructural size evolution proposed by Molinari and Ravichandran and the dislocation motion slip mechanism. This model effectively describes the mutation behavior of strain-rate sensitivity under fixed strain, particularly evident in an annealed OFHC. The predicted results of the model across a wide range of strain rates (10−4–106 s−1) and temperatures (77–1096 K) demonstrate relative errors generally within ±10% of the experimental values. Furthermore, the model is compared with five other models, including the mechanical threshold stress (MTS), Nemat-Nasser–Li (NNL), Preston–Tonks–Wallace (PTW), Johnson–Cook (JC), and Molinari–Ravichandran (MR) models. A comprehensive illustration of errors reveals that the proposed model outperforms the other five models in describing the plastic deformation behavior of OFHC. The error results offer valuable insights for selecting appropriate models for engineering applications and provide significant contributions to the field.


Introduction
Face-centered cubic (FCC) metals, such as aluminum and copper, are widely utilized in a variety of civil and military applications [1].However, during processes involving high strain rates and large strains, such as impact, penetration, and deformation processing, these materials exhibit complex nonlinear dynamic mechanical behavior [2].This behavior encompasses a broad range of strain rates and temperatures, resulting in significant coupling between strain rate, temperature, and strain hardening [3,4].Moreover, it is influenced by the historical effects of strain rate and temperature, presenting a formidable challenge in accurately constructing plastic constitutive models.
To address these challenges, various dynamic plastic constitutive models have been proposed in recent decades.These models can be classified into two main categories [1]: phenomenological models and models incorporating microstructural evolution.Phenomenological models are purely mathematical and conform to physical laws while possessing simple forms.They require fewer parameters to be calibrated and can produce accurate results within the calibrated range of plastic deformation.The Johnson-Cook (JC) model [5] is one of the most widely used phenomenological constitutive models that consider strain hardening and rate-temperature effects.It is based on the overstress type of models introduced by Malvern [6], Perzyna [7], and Campbell [8] and has been modified.Other variations of the JC model [9][10][11][12] have also been proposed using experimental data.However, these models may deviate significantly from reality for plastic deformation beyond the calibrated range of materials.
Experimental studies by Follansbee and Kocks [13] have shown that certain FCC metals, such as annealed oxygen-free high conductivity copper (OFHC), exhibit a sharp increase in flow stress when the strain rate exceeds a critical value (~10 4 s −1 ).Despite the high strain rate, plastic deformation remains dominated by the thermal activation of dislocation motion.This remarkable rate sensitivity is primarily attributed to the microstructural evolution during deformation, which is not accounted for in phenomenological models.Microstructural evolution refers to the changes in the arrangement and characteristics of dislocations at the microscale [14].These changes can include a rapid increase in dislocation density and a significant decrease in characteristic lengths of dislocations (e.g., dislocation spacing, thermally activated area, cell size), leading to a rapid rise in flow stress [15,16].Recently, Lea's research [17,18] highlights the paramount significance of microstructure and composition as the primary determinants leading to substantial variations in the dynamic mechanical response of copper under high strain rates, despite minor variations in purity.As a result, it is imperative to consider the impact of microstructural evolution on a material's mechanical performance when modeling.
In response to the influence of microstructural evolution on material dynamics, several constitutive models have been proposed [15].These models postulate that the storage distance of dislocations restricts their range of motion, resulting in a limited displacement that dislocations can undergo in unit time.This limitation increases the possibility of interactions between dislocations and defects, such as solutes and vacancies, leading to the accumulation and pileup of more dislocations and enhancing the strain hardening effect of the material.Examples of such models include the Mecking-Kocks (MK) model [14], Zerilli-Armstrong (ZA) model [19], mechanical threshold stress (MTS) model [13], Nemat-Nasser-Li (NNL) model [20], Gao-Zhang (GZ-1, GZ-2) models [3,21], and Austin-McDowell models [22,23].Among these models, the MTS model [13] captures the mutation phenomenon of strain hardening sensitivity through a single "continuity" equation with the Voce law.However, this model requires various experimental data and phenomenologically based expressions to determine the thermal component of the material's behavior, which cannot be explicitly stated, leading to a complex constitutive form that can be challenging to manipulate.On the other hand, the Preston-Tonks-Wallace (PTW) model [24] incorporates an enhanced Voce strain hardening criterion to characterize the rate-sensitive mutation phenomenon in materials.It employs piecewise functions within the intermediate regime (10 4 s −1 to 10 8 s −1 ) to depict thermally activated dislocation glide mechanisms or phonon-drag mechanisms, which may result in a discontinuous strain-rate sensitivity of stress (m = ∂lnσ/∂ln .ε) with fixed strain.Similarly, the GZ-2 model [3] also employs piecewise functions to characterize the strain-rate sensitivity mutations, resulting in discontinuous strain-rate sensitivity coefficients.The abovementioned constitutive models provide valuable insights into the effects of microstructural evolution on the strain hardening behavior of FCC metals, contributing to the development of accurate predictive models.
In addition to phenomenological and microstructure-based models, semiempirical physical models have also been developed.These models combine experimental data with physical principles and employ empirical parameters or fitting coefficients to simulate the evolution of microstructure-related internal variables.One notable model in this regard is the Molinari-Ravichandran (MR) model [25], which successfully predicts the plastic flow behavior of annealed copper across a wide range of strain rates spanning seven orders of magnitude.The MR model introduces a single internal variable representing the effective microstructural length related to temperature and strain rate, offering valuable insights for metal microstructure design.Other extensions of the MR model are proposed by Durren-berger et al. [26,27], considering the flow stress as the sum of internal stress (long-range interaction) and thermally activated effective stress, providing simplicity and flexibility in describing thermoviscoplastic behavior under extreme loading conditions.However, their application under high-temperature working conditions still needs improvement.
Motivated by the cell characteristic size evolution mechanism in the MR model, this study develops a simplified semiempirical physical constitutive model with a single "continuity" equation to capture stress upturning.It utilizes a normalized cell size evolution equation and incorporates the thermal activation of dislocation motion to capture the plastic dynamic mechanical behavior of FCC metals across a wide range of strain rates and temperatures.In this paper, a comprehensive explanation of the theoretical modeling methods for the components of thermal and athermal stresses in the flow stress is presented in Section 2. Section 3 provides the constitutive parameters for annealed OFHC and the sources of experimental values used to validate the predictive capability of the models.Additionally, a comparative analysis is conducted between the proposed model and other mature models, including the MTS, NNL, PTW, JC, and MR models, to evaluate their performance in predicting experimental results (Section 4).Finally, the main conclusions are summarized in Section 5.
This study aims to contribute to the advancement of constitutive modeling by capturing the complex behavior of materials during plastic deformation.By incorporating the microstructural evolution and the thermal activation of dislocation motion, the proposed model provides valuable insights into the dynamic mechanical behavior of FCC metals under various strain rates and temperatures.

Constitutive Modelling
Plastic flow during the thermally activated motion of dislocations is primarily governed by the movement of dislocations, which will encounter obstacles in the form of short-range and long-range barriers during their motion.Long-range barriers are nontemperature dependent and are influenced by factors such as grain boundaries, far-field forest dislocations, and other microstructural elements [1,20].These barriers can impact dislocation motion through mechanisms such as externally applied stress or chemical potential differences.Conversely, the height of short-range barriers is closely related to temperature.At absolute zero temperature (0 K), the height of short-range barriers reaches its maximum value because there is no thermal energy available for dislocations to overcome them.However, the atomic vibration amplitude in the material increases with temperature rising, thereby reducing the height of the short-range barriers.This increase in thermal energy assists dislocations in overcoming the short-range barriers [1].Therefore, short-range barriers can be overcome through a thermal activation mechanism.In FCC metals, short-range barriers may include forest dislocations, point defects (such as vacancies and self-interstitials), solute atoms, impurities, and precipitates.The presence of these short-range barriers plays a significant role in plastic deformation and the mechanical properties of materials.Thus, based on the nature of barriers and the theory of MTS, the flow stress of a material can be expressed as [13] σ where σ ath represents an athermal component of flow stress (rate independent), which is influenced by material structure properties.σ th represents a thermal component of flow stress (rate dependent), which characterizes the resistance overcome by mobile dislocations through diffusion or slip, driven by thermal activation effects, leading to material plastic deformation.σth is the thermal-activation-related mechanical threshold stress at 0 K, and Φ . ε, T is the coefficient of influence related to strain rate and temperature (<1), indicating the role of thermal activation effects.

Thermal Stress
The thermal stress is closely related to the microstructural features of the material, as it utilizes the thermal vibrations of the crystal lattice to assist dislocations in overcoming weak barriers (such as forest dislocations and Peierls barriers).This facilitates the diffusion or glide of dislocations, resulting in material deformation.

Thermal Activation Effect in Fixed Structures
During the thermal activation stage, the time required for dislocations to overcome barriers is primarily determined by the waiting time of t w .In the mechanism of thermally activated dislocation motion, the average waiting time before encountering short-range barriers can be expressed in a modified form of the Arrhenius equation [22]: where ν 0 represents the vibration frequency (~10 11 s −1 ) of the dislocation before encountering barriers [28], k B is the Boltzmann constant, T is the temperature, and ∆G is the thermal activation free energy (Gibbs free energy) required for dislocations to overcome the barriers, which is a function associated with the shape of the barrier (e.g., rectangular, square, triangular, parabolic, exponential) as well as stress.The modified form indicates that the waiting time for dislocations in a continuous slip system (∆G = 0) is zero.Kocks and Ashby [28] proposed a general expression for activation energy: where G 0 is the reference free energy at 0 K (= g 0 µ 0 b 3 ), g 0 is the normalized free energy, µ 0 is the shear modulus (= µ(0K)), and b is the magnitude of the Burgers vector.p and q are a pair of parameters that characterize the shape of the barrier.q corresponds to the shape of the short-range barrier, where for a quasi-parabolic barrier, the range is 1 ≤ q ≤ 2. p describes the rate of decay of the field at longer distances, with a value of 0 < p ≤ 1.A value of p = 1 indicates no decay, while smaller values indicate faster decay of the activation energy at longer distances.
According to Orowan's law (γ = N m b∆l), the equivalent true strain rate can be expressed as follows: .
where ∆l represents the average distance traveled by dislocations between barriers, and N m denotes the mobile dislocation density.In the high-strain-rate regime, the mobile dislocation density dominates the total dislocation density, that is N m ≈ N tot [3,29].M is the Taylor factor, where σ = Mτ, ε = γ/M, τ, and γ represent the shear stress and strain, respectively.The reference strain rate is defined as .ε 0 = N m b∆lν 0 /M with a general order of 10 7 s −1 [13].
By substituting Equation (4) into Equation ( 3) and considering the relationship of Φ .ε, T = σ th / σth , the following expression is obtained: When the temperature exceeds the critical temperature T cr , the thermal activation stress term disappears, and the flow stress is solely determined by the athermal stress component, that is, σ = σ ath .The expression for the critical temperature is given by

Microstructural Evolution
Based on the analysis of microstructural evolution in deformation processing experiments, one of the key features of metal deformation is the variation of cell size with strain.Sevillano's experimental results [30] demonstrate that the cell size decreases with increasing equivalent plastic strain.For certain metals, the ratio of the average cell size δ to its saturated size δ s is a function of equivalent plastic strain, and the data points cluster around a curve that is independent of the metal type [31].The normalized cell size δ, defined as the ratio of the average cell size to its initial value δ 0 ( δ = δ/δ 0 ≤ 1), can be described by the following evolution equation: where δ r is a dimensionless coefficient that represents the rate of cell refinement, while δ s denotes the normalized saturated cell size (= δ s /δ 0 ).Both are functions related to strain rate and temperature.It is worth noting that this study is consistent with the theoretical framework proposed by Molinari and Ravichandran [25], which suggests that δ is a generalized characteristic length of the microstructure and does not have a direct relationship with geometric dimensions.Therefore, the use of the dimensionless size δ captures the essential trends in simulating microstructural evolution, where specific cell sizes of δ, δ 0 , and δ s become less important, and there is no need to precisely capture quantitative results regarding the microstructure.The integral Equation ( 7) yields the normalized cell size δ, denoted as The coefficients δ s and δ r can be expressed using the following empirical relationships [16]: where T r and .ε r represent the reference temperature and strain rate, respectively.a s and a r are non-negative constants, and ξ s , ν s , ξ r , and ν r are constant values ranging from 0 to 1, controlling the dependence of δ s and δ r on strain rate and temperature.δ s0 represents the maximum reference saturated cell size δ s0 in comparison with its initial size ( δ s ≤ δ s0 < 1), while δ r0 is the minimum reference refinement rate (δ r0 ≤ δ r ).The normalized saturated cell size δ s decreases with an increasing deformation rate and increases with an elevated temperature.Conversely, the cell refinement rate δ r accelerates with higher deformation rates and decelerates with an increasing temperature.The evolution of these two parameters indicates that under low-temperature and high-strain-rate conditions, the dynamic cell refinement rate is faster, resulting in smaller achievable saturated sizes.In Equation ( 9), a modification is made to the temperature term, suggesting that the evolution of cell size becomes insensitive to temperature when it is lower than T r .
Staker and Holt [32] proposed a simple relationship between the average cell size δ and the average total dislocation density N tot , which can be expressed as follows: where K c is a material constant, with a value of 16 for OFHC.
Therefore, the square root of the normalized dislocation density is inversely proportional to the normalized cell size δ: By substituting the saturation value into Equation ( 11), the normalized initial saturated cell size δ s0 = √ N 0 /N s0 can be obtained, where N s0 represents the minimum reference value (~10 14 m −2 [33]) of the saturated dislocation density.

Threshold Stress and Thermal Stress
In the theory of MTS, the threshold stress associated with thermal activation is related to microstructural evolution [14] and can be expressed in combination with Equation (11) as where α represents the intensity of dislocation interaction (<1), and µ(T) is the temperaturedependent shear modulus.σ0 = M αµ 0 b √ N 0 represents a constant threshold stress value associated with the initial microstructure state at 0 K.
By substituting Equation ( 8) into Equation ( 12) and combining it with Equation ( 5), the thermal activation stress can be expressed as here, µ is the normalized shear modulus (= µ(T)/µ 0 ).The thermal activation stress consists of two components: the structural evolution (threshold stress) part and the strain rate and temperature response of the material under a fixed structure.The structural evolution (normalized cell size) is related to strain rate and temperature.Therefore, the strain-rate sensitivity of the thermal activation stress can be divided into two parts: the rate sensitivity of material strain hardening caused by structural evolution and the instantaneous rate sensitivity of stress [15,34].Equation ( 13) can also be written in the standard Voce-Kocks format [35]: where σ ε=∞ th and σ ε=0 th represent the saturation and initial thermal activation stresses, respectively, with the following dependencies on strain rate and temperature:

Athermal Stress
The athermal stress is a result of the interaction between dislocations and long-range obstacles, known as the Hall-Petch strengthening effect [15].The long-range obstacles mainly include thermally independent grain boundaries and far-field forest dislocations.As the grain size decreases (~µm), the number of grain boundaries increases, or the density of far-field forest dislocations increases with increasing deformation.These long-range obstacles hinder dislocation motion and capture more dislocations, thereby enhancing the material strength [36,37].
The Hall-Petch strengthening effect caused by the grain size effect can be expressed as follows [37,38]: where σ 0 represents the friction stress required for the motion of dislocations, which is typically negligible for FCC metals.k s is the strengthening coefficient, and D 0 is the initial grain size.The density of far-field forest dislocations N f can be calculated using the rate-independent dislocation density evolution equation proposed by Klepaczko [15], where the effective forest dislocation density N f is represented as where M II and k a0 are constants representing the proliferation rate and annihilation rate of dislocations, respectively.As the strain increases, the effective dislocation density approaches a saturation value, which is The athermal stress can be expressed as [1,29] where B represents the stress at the saturation dislocation density at 0 K, denoted as B = α d µ 0 b N fs , and α d is a material constant (<1).The athermal stress indirectly achieves the thermal softening effect of the material through the relationship between the shear modulus (elastic field) and temperature.It realizes strain hardening of the material through the generation and storage of dislocations during deformation, while the deformation rate does not affect the evolution of nonthermal stress.

Flow Stress and Hardening Rate
Finally, combining Equations ( 13) and ( 18), the constitutive relationship for the flow stress of FCC metals is obtained as The strain hardening rate with respect to stress is given by where the first term represents the strain hardening rate induced by the evolution of long-range obstacles in the microstructure, such as grain size and grain boundaries.The hardening rate initially decreases inversely proportional to the athermal stress and follows a linear function with a constant slope k a0 at higher stresses.The second term accounts for the strain hardening rate due to thermal stress arising from the evolution of short-range obstacles, such as line defects represented by dislocation spacing.Under stable loading conditions, the hardening rate decreases linearly as the thermal stress increases, where δ r is a constant.
Observing the model's behavior under constant loading conditions (stable strain rate and temperature), the overall strain hardening rate exhibits a rapid decrease as the deformation intensifies due to the inverse proportional function, reflecting the hardening characteristics of FCC metals at stage III [14,31].Once the material reaches a certain level of strengthening, the overall strain hardening rate gradually decreases linearly until it approaches zero, with a slope approximately equal to −(k a0 + δ r ).This behavior is consistent with the hardening characteristics of stage IV [31,39].Since k a0 is a material constant, the sensitivity of the hardening rate to strain rate and temperature primarily depends on the microstructural refinement rate δ r .

Shear Modulus
The Nadal-Le Poac (NP) shear modulus model [40] replaces the empirical temperature dependence of the shear modulus in the SCG model with an equation based on the Lindemann melting theory.This modification enhances the robustness of the model when applied under high-temperature and high-pressure conditions.Moreover, the NP model mitigates the instantaneous decrease of the shear modulus at the melting point.The zero-pressure form of the NP shear modulus model can be expressed as where ρ m represents the material's standard density, T m denotes the melting point of the material under standard conditions, C is a material constant related to the Lindemann constant, m u signifies the atomic mass, Γ(T/T m ) controls its nonlinear behavior when T/T m ≈ 1, and ζ m 1 serves as a fitting parameter that governs the overall fit quality.

Adiabatic Temperature Rise
During the metal deformation process, it was experimentally found that the majority of the energy supplied to effect plastic deformation converted to heat, contributing to an increase in temperature [41,42].The conversion coefficient β generally falls within the range of 0.9 to 1.0 for high deformation rates (>1 s −1 ) in adiabatic processes [42,43].If adiabatic shear banding occurs under high strain rates, the heat generated during plastic deformation can lead to a rapid increase in temperature within the localized region of the adiabatic shear band.The equation connecting temperature raising and plastic work is presented in an integral format [44]: where T 0 is the initial temperature, and C p is the specific heat of the material at constant pressure.In the calculation of flow stress during adiabatic processes, it is necessary to update the temperature state variable at each incremental strain step.

Application
This section presents the material parameters used in the model for annealed OFHC and compares the predicted results with experimental data.The results cover the influence of strain rate on flow stress under fixed strains, as well as the deformation behavior of annealed OFHC across a wide range of strain rates (4 × 10 −4 to 6.4 × 10 5 s −1 ) at room temperature, and a broad temperature range (77 to 1096 K) at a strain rate of 4000 s −1 .

Parameter Calibration
The physical constants of annealed OFHC [43] and the parameters in the NP shear modulus model at zero hydrostatic pressure [40] are shown in Table 1.In the component of athermal stress, Meyers et al. [37] conducted a study and found that the value of k s for OFHC decreases linearly with increasing strain under low strain rates, ranging from 56 to 260 MPa √ µm.However, at high strain rates, the value remains approximately constant.This behavior is taken into account in the ZA model [19], where an intermediate value is employed.The coefficients for the athermal stress in the model are provided in Table 2.The material parameters for the thermal stress component are described in the following two tables.Table 3 provides the parameters that are independent of microstructural evolution, calculated σ ˆ0 as 50.2MPa based on the data.The parameters for the normalized cell size evolution model are proposed by Molinari and Ravichandran [25] and illustrated in Table 4.By utilizing the data in the table and Equation (11), the minimum value of the reference saturated dislocation density is approximately 7 × 10 14 m −2 , which consistent with the actual physical laws.

Simulation Results
The experimental values for a comparison with the model calculations can be found in Appendix A. In Figure 1, the relationships between flow stress, normalized cell size, and strain rate for OFHC at fixed strains and room temperature are depicted.It is noteworthy that the sensitivity of flow stress to strain rate undergoes a sudden change at 10 4 s −1 , which aligns with the evolution pattern of the normalized reciprocal cell size δ −1 .This observation suggests that the sensitivity of flow stress to strain rate is dominated by the microstructural evolution law.Specifically, when .ε < 10 4 s −1 , the normalized reciprocal cell size lines are horizontally parallel with fixed strains.This behavior arises from the fact that the saturation size and refinement rate, as controlled by Equation ( 9), remain approximately constant, indicating that the microstructure exhibits low sensitivity to changes in strain rate.Additionally, the sensitivity of flow stress to strain rate is influenced by the term Φ .ε, T .At extremely low strain rates ( .ε < 10 −3 s −1 ), the value of Φ .ε, T exhibits minimal variation with strain rate, implying that the flow stress is predominantly determined by the deformation level and temperature.
Figures 2 and 3 present the relationship between flow stress and strain within a wide range of strain rates at room temperature and a wide range of temperatures at a fixed strain rate for OFHC, respectively.The results demonstrate a strong agreement between the model's predicted values and the experimental data, which will be further analyzed in the subsequent section to quantify the error.In Figure 2, when the strain exceeds ε 1 (≈0.2), a distinct change in the surface occurs near 1 s −1 for a fixed strain.Within a narrow range of strain rates, an intriguing phenomenon emerges as the flow stress decreases with an increasing strain rate.This intriguing behavior can be attributed to the consideration of temperature increase induced by plastic deformation when .ε > 1 s −1 , where the thermal softening effect of the material surpasses the strengthening effect associated with strain rate.
cell size lines are horizontally parallel with fixed strains.This behavior arises from that the saturation size and refinement rate, as controlled by Equation ( 9), remain a imately constant, indicating that the microstructure exhibits low sensitivity to cha strain rate.Additionally, the sensitivity of flow stress to strain rate is influenced term ( )  Φ ε,T .At extremely low strain rates (  ε < 10 −3 s −1 ), the value of ( )  Φ ε,T e minimal variation with strain rate, implying that the flow stress is predominantly mined by the deformation level and temperature.Figures 2 and 3 present the relationship between flow stress and strain within range of strain rates at room temperature and a wide range of temperatures at strain rate for OFHC, respectively.The results demonstrate a strong agreement b the model's predicted values and the experimental data, which will be further an in the subsequent section to quantify the error.In Figure 2, when the strain exc (≈0.2), a distinct change in the surface occurs near 1 s −1 for a fixed strain.Within a range of strain rates, an intriguing phenomenon emerges as the flow stress decreas an increasing strain rate.This intriguing behavior can be attributed to the consid of temperature increase induced by plastic deformation when  ε > 1 s −1 , where the softening effect of the material surpasses the strengthening effect associated with stra 10 −2 10 0 10 2 10 4 10 6     A3, with blue indica stress and red indicating high stress.

Relative Error Analysis
Based on the experimental results of OFHC across a wide range of strain ra deformation temperatures, Biswajit Banerjee conducted a meticulous evaluation i search [47,48].He identified the hierarchical order of the model's evaluation error MTS, ZA, JC, and SCGL (Steinberg-Cochran-Guinan-Lund [49]).Notably, the model's microstructure evolution is predicated upon the MR model, while the va experiments predominantly stem from Nemat-Nasser and Li [20].To enhance the standing of the predictive accuracy of the constitutive models, the paper unde  A3, with blue indicating low stress and red indicating high stress.

Relative Error Analysis
Based on the experimental results of OFHC across a wide range of strain rates and deformation temperatures, Biswajit Banerjee conducted a meticulous evaluation in his research [47,48].He identified the hierarchical order of the model's evaluation errors: PTW, MTS, ZA, JC, and SCGL (Steinberg-Cochran-Guinan-Lund [49]).Notably, the present model's microstructure evolution is predicated upon the MR model, while the validation experiments predominantly stem from Nemat-Nasser and Li [20].To enhance the understanding of the predictive accuracy of the constitutive models, the paper undertook a thorough analysis of errors and evaluated the performance of not only the present model but also five alternative constitutive models: MTS, PTW, JC, NNL, and MR.This comprehensive assessment encompassed a total of six distinct models, providing valuable insights into their respective prediction capabilities.
The relative error values E i,j are computed using Equation ( 23): where σ pre i,j denotes the model-predicted flow stress and σ exp i,j represents the experimental flow stress.The subscripts 'i' and 'j' signify the relative for the i-th data point under the j-th deformation condition.
Figures 4 and 5 showcase the error curves representing the disparities between the predicted flow stress of OFHC using six distinct models and the corresponding experimental results.The order in which these models are presented is as follows:  The relative error curves of Figure 5a,b,f exhibit a high degree of similarity.These error curves demonstrate a consistent pattern, where the curves gradually shift upwards as the temperature increases, indicating that the predicted values are higher compared with the experimental results at high temperatures.Furthermore, for model (c) NNL, its predicted values are consistently higher than the experimental values at both high and low temperatures.This behavior can be attributed to the fact that, under high temperature, the solid-state diffusion rate in the crystal significantly increases, promoting grain growth.Larger grains may contain more structural defects or dislocations, resulting in a decrease in material strength.Since these models do not account for microstructural coarsening, the predicted results at high temperatures tend to be higher than the experimental ones.On the other hand, the PTW and JC models underestimate the material's flow stress at higher temperatures, and as the strain increases, their predictions deviate even further from the experimental results.From the perspective of the error curves, the MTS model exhibits the best predictive performance, while there is still room for improvement in the present model's predictions at higher temperatures (>0.7 Tm).As shown in Figure 4a, the present model exhibits relatively larger relative errors (at ±20%) in predicting results at room temperature and strain rates of 8000 s −1 and 1.3 × 10 4 s −1 .However, under other strain-rate conditions, the relative errors of the model are mainly within ±10%.Notably, the relative error curves of the other five models at the above two strain rates exhibit similar trends to that of the present model.Therefore, it is plausible that these significant discrepancies may stem from experimental errors at these two strain rates.Regarding extremely high strain rates (6.4 × 10 5 s −1 ), the present model, PTW, and MR demonstrate relatively accurate prediction results, while the MTS model exhibits slightly higher values compared with the present.Both the NNL and JC models significantly underestimate the experimental results at low strains, with relative errors approaching −40%.Nevertheless, the accuracy of NNL predictions notably improves as the strain increases, while the JC model remains approximately −20%.Moreover, under all other strain-rate conditions, the JC model overestimates material flow stress at low strains with relative errors exceeding +20%.Excluding the JC model, the other models exhibit satisfactory consistency with the experimental results across the range of 0.001-8500 s −1 .However, at an extremely low strain rate of 4 × 10 −4 s −1 , all models predict higher values than the experimental results.Generally speaking, the predictions of the present model outperform those of the other models.This discrepancy can be attributed to the fact that the dislocation slip velocity is slower at relatively low strain rates, resulting in prolonged interaction time between dislocations and internal defects or barriers within the material.Consequently, dislocation diffusion or dispersion mechanisms may become the dominant deformation mechanisms, leading to a decrease in material strength and causing the model to predict higher flow stress values compared with the experimental data [19,50].

Comprehensive Error Analysis
To provide a more comprehensive assessment of the predictive capabilities of the abovementioned six constitutive models for the flow stress of OFHC, three statistical average errors were employed to evaluate the model predictions.).The calculation formulas for these metrics are as follows: where j I represents the total number of data points with the deformation condition of index j, j E denotes the average value of the relative error curve j, and J refers to the total number of experimental data that satisfy the statistical criteria.For example, there are four sets of experimental data with ε ≤ 10 3 s −1 at room temperature, that is, J = 4.The statistical average errors for the predictions of the six models under different strain-rate gradients at room temperature are presented in Figure 6.It is observed that the majority of the models perform well in accurately capturing the flow stress behavior of the The relative error curves of Figure 5a,b,f exhibit a high degree of similarity.These error curves demonstrate a consistent pattern, where the curves gradually shift upwards as the temperature increases, indicating that the predicted values are higher compared with the experimental results at high temperatures.Furthermore, for model (c) NNL, its predicted values are consistently higher than the experimental values at both high and low temperatures.This behavior can be attributed to the fact that, under high temperature, the solid-state diffusion rate in the crystal significantly increases, promoting grain growth.Larger grains may contain more structural defects or dislocations, resulting in a decrease in material strength.Since these models do not account for microstructural coarsening, the predicted results at high temperatures tend to be higher than the experimental ones.On the other hand, the PTW and JC models underestimate the material's flow stress at higher temperatures, and as the strain increases, their predictions deviate even further from the experimental results.From the perspective of the error curves, the MTS model exhibits the best predictive performance, while there is still room for improvement in the present model's predictions at higher temperatures (>0.7 T m ).

Comprehensive Error Analysis
To provide a more comprehensive assessment of the predictive capabilities of the abovementioned six constitutive models for the flow stress of OFHC, three statistical average errors were employed to evaluate the model predictions.These include the mean absolute relative error of plastic work (MARE-P.WK, E AR.W ), the mean absolute relative error of flow stress (MARE-F.S, E AR.Y ), and the mean maximum absolute error of flow stress (MMAE-F.S, E MA.Y ).The calculation formulas for these metrics are as follows: where I j represents the total number of data points with the deformation condition of index j, E j denotes the average value of the relative error curve j, and J refers to the total number of experimental data that satisfy the statistical criteria.For example, there are four sets of experimental data with .ε ≤ 10 3 s −1 at room temperature, that is, J = 4.The statistical average errors for the predictions of the six models under different strain-rate gradients at room temperature are presented in Figure 6.It is observed that the majority of the models perform well in accurately capturing the flow stress behavior of the material at low strain rates, with the exception of the JC model, which exhibits a higher MMAE-F.S value.The strain-rate sensitivity of the flow stress undergoes a notable transition within the range of 10 3 s −1 -10 4 s −1 .Specifically, the present model and the MTS and NNL models exhibit relatively low statistical average errors of approximately 10% within this range, while the JC and MR models demonstrate larger errors exceeding >20%.This finding indicates that the abrupt change in flow stress sensitivity of OFHC is significantly influenced by microstructural evolution.Moreover, at extremely high strain rates, the values of MMAE-F.S for the MTS, NNL, and JC models exceed 20%, suggesting that these three models may not yield accurate results under such extreme high-strain-rate conditions.The statistical average errors for the predictions of the six models under different temperature gradients at a strain rate of 4000 s −1 are shown in Figure 7.It is observed that the present model and the MTS and NNL models consistently demonstrate relatively small statistical average errors across all temperature ranges (maximum value is approximately 10%).The PTW and MR models demonstrate statistical average errors below 10% at temperatures below 0.4 T m .However, as the temperature increases, the statistical average errors increase as well.Notably, the MR model exhibits statistical average errors exceeding 20% within the intermediate temperature range and even surpasses 60% at high temperatures.Therefore, the MR and PTW models are deemed unsuitable for high-temperature conditions.Moreover, throughout the entire temperature range, the JC model displays statistical average errors ranging from 15% to 30%.This indicates that the standard JC model is not suitable for simulations involving strain rates equal to or greater than 4000 s −1 .The JC model, being an empirical model that is based on experimental test data for parameter fitting, lacks a robust physical foundation.The information should only be utilized within the specified parameters for accurate results and should not be applied to conditions outside of those parameters.Otherwise, it may result in substantial errors and inaccurate predictions of material behavior.The statistical average errors in model predictions across all strain rates and temperature ranges are comprehensively evaluated and presented in Figure 8.The enclosed triangle area formed by the lines serves as a measure of the model's prediction accuracy, with a smaller area indicating higher precision.When considering different strain rates, the sizes of the triangles formed by the three statistical average errors follow this order: present ≈ PTW < MTS < MR < NNL < JC.Similarly, for varying temperatures, the triangle sizes adhere to this sequence: present ≈ NNL < MTS < PTW < JC < MR.From the results, it is evident that the proposed model demonstrates favorable simulation outcomes across a broad range of strain rates and temperatures.The statistical average errors in model predictions across all strain rates and temperature ranges are comprehensively evaluated and presented in Figure 8.The enclosed triangle area formed by the lines serves as a measure of the model's prediction accuracy, with a smaller area indicating higher precision.When considering different strain rates, the sizes of the triangles formed by the three statistical average errors follow this order: present ≈ PTW < MTS < MR < NNL < JC.Similarly, for varying temperatures, the triangle sizes adhere to this sequence: present ≈ NNL < MTS < PTW < JC < MR.From the results, it is evident that the proposed model demonstrates favorable simulation outcomes across a broad range of strain rates and temperatures.

Conclusions
This study develops a semiempirical physical constitutive model based on a normalized cell size evolution model, along with the dislocation density evolution mechanism and

Conclusions
This study develops a semiempirical physical constitutive model based on a normalized cell size evolution model, along with the dislocation density evolution mechanism and the impact of dislocation motion on plastic deformation.This single and continuous model, which captures the sensitivity of microstructure evolution to strain rates, is capable of effectively capturing the upturning in flow stress at strain rates exceeding the critical value (~10 4 s −1 ) at a fixed strain.Furthermore, it demonstrates remarkable precision across a wide range of strain rates (4 × 10 −4 -6.4 × 10 5 s −1 ) and temperatures (77-1096 K).Despite this, the model is not without its limitations.The high-temperature environment presents challenges in the nonlinear mechanical behavior of materials, such as grain growth and structural defect generation.The ramifications of these phenomena on the plastic deformation process and the impact of temperature on microstructural evolution warrant further examination.
Upon comparison of the predictive capabilities of the present model with those of the MTS, NNL, PTW, JC, and MR models, it was determined that the microstructure-based constitutive model, particularly the present model, offers a more accurate representation of material plastic deformation behavior under high strain rates.The proposed model surpasses others in terms of its predictive accuracy.The results of this investigation hold a significant value as they provide guidance for selecting appropriate models for engineering applications and augmenting our understanding of the rheological behavior of OFHC.

Figure 1 .
Figure 1.Relationship between flow stress, normalized cell size, and strain rate of anneale at fixed strain and room temperature.

Figure 1 .Figure 2 .
Figure 1.Relationship between flow stress, normalized cell size, and strain rate of annealed OFHC at fixed strain and room temperature.Materials 2023, 16, x FOR PEER REVIEW

Figure 2 .
Figure 2. Flow stress-strain curved surface of annealed OFHC over a wide strain rate range at room temperature, where the dots are experimental data referred in Table A2, with blue indicating low stress and red indicating high stress.

Figure 2 .
Figure 2. Flow stress-strain curved surface of annealed OFHC over a wide strain rate range temperature, where the dots are experimental data referred in TableA2, with blue indica stress and red indicating high stress.

Figure 3 .
Figure 3. Flow stress-strain curved surface of annealed OFHC over a wide temperature ra  ε = 4000 s −1 , where the dots are experimental data referred in TableA3, with blue indica

Figure 3 .
Figure 3. Flow stress-strain curved surface of annealed OFHC over a wide temperature range with .ε =4000 s −1 , where the dots are experimental data referred in TableA3, with blue indicating low stress and red indicating high stress.
(a) present, (b) MTS, (c) NNL, (d) PTW, (e) JC, and (d) MR.Models (a) to (d) encompass constitutive models that incorporate microstructure evolution, while model (e) is an empirical model and (f) is a semiempirical physical model.

Figure 4 .
Figure 4.The relative error curve between the predicted values of models and the experimental results over a wide range of strain rates at room temperature.The models are (a) Present, (b) MTS, (c) NNL, (d) PTW, (e) JC, and (f) MR.

Figure 4 .
Figure 4.The relative error curve between the predicted values of models and the experimental results over a wide range of strain rates at room temperature.The models are (a) Present, (b) MTS, (c) NNL, (d) PTW, (e) JC, and (f) MR.

Figure 5 .
Figure 5.The relative error curve between the predicted values of models and the experimental results over a wide range of temperatures at a strain rate of 4000 s −1 .The models are (a) Present, (b) MTS, (c) NNL, (d) PTW, (e) JC, and (f) MR.

Figure 5 .
Figure 5.The relative error curve between the predicted values of models and the experimental results over a wide range of temperatures at a strain rate of 4000 s −1 .The models are (a) Present, (b) MTS, (c) NNL, (d) PTW, (e) JC, and (f) MR.

Figure 8 .
Figure 8. Comprehensive evaluation of three statistical average errors in model prediction results: (a) in different strain rates and (b) in different temperatures.

Figure 8 .
Figure 8. Comprehensive evaluation of three statistical average errors in model prediction results: (a) in different strain rates and (b) in different temperatures.

Table 1 .
Parameters in the NP shear modulus model and temperature evolution model of OFHC.

Table 2 .
Material coefficients in the athermal stress component.

Table 3 .
Material parameters in the thermal stress component that are independent of microstructure.

Table 4 .
Parameters of the evolution equation for normalized cell size [25].

Table A3 .
Experimental results on the flow stress-strain relationship of OFHC over a wide range of temperatures at a fixed strain rate.