Material Behavior Description for a Large Range of Strain Rates from Low to High Temperatures: Application to High Strength Steel

Current needs in the design and optimization of complex protective structures lead to the development of more accurate numerical modelling of impact loadings. The aim of developing such a tool is to be able to predict the protection performance of structures using fewer experiments. Considering only the numerical approach, the most important issue to have a reliable simulation is to focus on the material behavior description in terms of constitutive relations and failure model for high strain rates, large field of temperatures and complex stress states. In this context, the present study deals with the dynamic thermo-mechanical behavior of a high strength steel (HSS) close to the Mars® 190 (Industeel France, Le Creusot, France). For the considered application, the material can undergo both quasi-static and dynamic loadings. Thus, the studied strain rate range is varying from 10−3–104 s−1. Due to the fast loading time, the local temperature increase during dynamic loading induces a thermal softening. The temperature sensitivity has been studied up to 473 K under quasi-static and dynamic conditions. Low temperature measurements (lower than the room temperature) are also reported in term of σ − ε | ε ˙ , T curves. Experimental results are then used to identify the parameters of several constitutive relations, such as the model developed initially by Johnson and Cook; Voyiadjis and Abed; and Rusinek and Klepaczko respectively termed Johnson–Cook (JC), Voyiadjis–Abed (VA), and Rusinek–Klepaczko (RK). Finally, comparisons between experimental results and model predictions are reported and compared.


Introduction
Mechanical response of steels can be very sensitive to the strain rate and the temperature, as shown by numerous studies such as [1][2][3]. The use of these materials for protective systems to ballistic impacts requires the knowledge of their mechanical behavior, due to the fact that these structures undergo a wide range of solicitation. Moreover, the development of such protection includes several material, leading to complex structures. In fact, the number of parameters to consider (thickness, arrangement and orientation of layers, type of materials, projectile shapes, etc.) is too significant to perform a parametric study using expensive and time-consuming experimental tests. Therefore, the design and optimization of such multi-layered structures is made possible through the use of numerical simulations.

Experimental Procedures
In order to determine the hardening, the strain rate, and temperature dependency of the studied material, compression tests were performed. The studied range of strain rate was varying from 10 −3 -10 4 s −1 . Regarding the temperature sensitivity, experiments at various temperatures were performed both in quasi-static and dynamic conditions. Regarding the quasi-static case, the strain rate was fixed at 0.001 s −1 through a strain rate control and the temperature was varying from room temperature up to 473 K. In dynamic conditions, the strain rate was in the order of 10 3 s −1 and the temperature was varying from 173 K up to 473 K with a step of 50 K. The compression samples used for the quasi-static experiments were cylindrical specimens with an initial diameter d 0 of 6 mm and an initial length l 0 of 6 mm. For the split Hopkinson pressure bar (SHPB) and direct impact experiments, compression specimens had an initial diameter d 0 of 8 mm and an initial length l 0 of 4 mm. This enables one to obtain a ratio of s 0 = l 0 /d 0 = 0.5, required to reduce the inertia effect [16][17][18]. For experiments under quasi-static loading or at various temperatures, each condition was performed three times. Each experiment at room temperature and high strain rate was considered as a unique condition, due to the difficulty of obtaining identical strain rates with accuracy.

Experiments at Room Temperature
Uniaxial quasi-static experiments were performed using an electromechanical press. The use of the internal load cell sensor and an extensometer allowed to measure the values of the force F(t) and displacement δ(t) respectively. Experiments at strain rates~500 and~2000 s −1 were done using a SHPB set-up. It is composed of a striker, an input and an output bar. All these parts have a diameter of 20.6 mm. The input and output bars are respectively 1900 mm and 1300 mm long, while the length of the striker is varying from 200-600 mm to change the loading time and the strain level induced into the specimen. From the gage measurement of the incident (ε i ), reflected (ε r ), and transmitted (ε t ) waves, the true stress and the true strain are calculated using the elastic waves propagation theory. The strain rate depending on time is defined as follows, Equation (1) whereε represents the nominal strain rate, C 0 the elastic wave velocity in the bars, and E bar and ρ bar the Young's modulus and the density of the bars, respectively. The nominal strain, Equation (2) is obtained using the previous equation Equation (1) as Moreover, for material such as high strength steel (HSS), punching effect can induce an overestimation of the strain level induced to the specimen. An analytic solution was developed by Safa and Gary [19] to take into account this effect, Figure 1. A comparison is reported using or not this correction. The variation of strain ∆ε is proportional to the stress level and therefore increases with strain for material with positive hardening. Thus, it is observed that the puncture effect is reduced to ∆ε max ≈ 0.004 corresponding to~2.6% for a strain level of ε max = 0.15. For this reason, the correction is applied during the test analyses. Concerning the true stress in dynamic compression, it is obtained using the transmitted wave ε t , Equation (3). This relation assumes the mechanical equilibrium during the test.
where S bar is the section of the bars and S 0 the section of the specimen. In general, the stress value is affected by the friction value µ between the bars and the specimen. Moreover, inertia may act on the stress in addition to friction. Due to these two effects, the stress level measured during the test must be corrected. Details of the related correction can be found in [16]. Highest strain rate experiments were performed using direct impact experiments [18,20]. The stress was obtained by the same procedure as the one used for SHPB and the strain was determined by optical measurement.
During dynamic experiments, a part of the plastic work is converted into heat and cannot be dissipated along the specimen, contrary to quasi-static loading. In the last case, the temperature is therefore assumed as constant and equal to the initial temperature imposed to the specimen. In dynamic conditions corresponding to adiabatic conditions, the temperature increases ∆T during the process of plastic deformation is defined in Equation (4) as where C p is the specific heat at constant pressure, and ρ s is the specimen density. The parameter β is the Taylor-Quinney coefficient and corresponds to the part of the plastic work converted into heat and assumed for metals equal to 0.9. The values used in this work are reported in Table 1. In addition, and as discussed previously, thermal softening appears at high strain rates for deformations larger than 0.1. Therefore, an experimental study was carried out to determine the temperature sensitivity of the studied HSS material. Experiments at high temperatures for a strain rate equal to 0.001 s −1 have been performed using a climatic chamber surrounding the sample and the compression platens. These platens were heated to an imposed temperature, and the specimen was then placed and heated during 30 min before the test. This procedure was used to avoid cooling the specimen by thermal conduction with the compression tools. The tests were performed with a constant strain rate controlled by an extensometer fixed on the compression platens. Each test was performed three times to ensure the repeatability of the experiments. High temperature experiments under dynamic conditions were carried out using a specific furnace heating the specimen and the bar ends. This furnace is similar to the one used by Rusinek et al. [21]. A description of this apparatus is given in Figure 2. If a hot specimen is placed in contact with the bars at room temperature, the thermal conduction induces a fast cooling of the specimen as described in [18]. To avoid this problem, the ends of the bars were heated in the furnace. All parts were heated for 30 min to reach a homogeneous temperature distribution within the specimen. Experiments at low temperature were also carried out in the dynamic regime, using liquid nitrogen mixed with ethanol to adjust the temperature. Ethanol was chosen because of its melting temperature close to 159 K. The bars were in contact with the mixture to decrease their temperature, as described in Figure 3. This study assumes that the evolution of the elastic properties is negligible for these temperatures [22]. Moreover, no bending of the bars was observed during the test. The specimen was not in contact with the mixture, but was cooled by conduction. An instrumented specimen with a thermocouple was used to determine the temperature and the waiting time required before dynamic compression testing. The temperature was set by adding either liquid nitrogen or ethanol. Moreover, a thick steel plate was placed at the bottom of the cooling box to increase the thermal inertia of the mixture. Once the temperature is stabilized, the instrumented specimen was removed and a sample was placed between the bars. The temperature recorded inside the instrumented specimen is shown in Figure 4. A stable temperature is reached by conduction after a cooling time between 30 and 40 s once the specimen is placed between the bars. Moreover, as the bars are cooled prior to the test, their thermal inertia allows one to maintain the temperature of the specimen for a duration long enough to perform the test. For example, for a target temperature of 173 K, the specimen is maintained between 172 K and 174 K for 35 s (from 30-65 s). For a target temperature of 273 K and 223 K, the temperature of the specimen is maintained (at ±1 K) for more than a minute.

Experimental Results
The true stress values for each strain rate for a given strain level are represented in Figure 5. A plastic strain value below 0.1 was chosen to reduce the adiabatic heating effect on the stress. Therefore, the thermal softening effect is assumed to be negligible for this strain level. These experiments highlight a stress increase with the strain rate. Moreover, it has to be noticed that the strain rate sensitivity is much higher under dynamic conditions in comparison with low strain rates. The strain rate sensitivity, defined by m = ∂ log σ/∂ logε, is equal to 0.0021 for quasi-static loading and reaches 0.058 for dynamic loading. Thus, a nonlinear strain rate sensitivity is observed for the studied material and in general for metals [23][24][25]. The true stress values for each temperature and an imposed plastic strain level are reported in Figure 6. It has to be noted that the considered temperatures are corresponding to the initial temperatures (at the beginning of the test). The stress decreases as the temperature increases. The temperature sensitivity, defined by ν = ∂ log σ/∂ log T, is equal to −0.21 in dynamic loading conditions and decrease to −0.17 in quasi-static conditions. Experiments have highlighted the effect of strain rate and temperature on the mechanical behavior of the studied HSS material.

Material Characterization of the Experimental Results
As discussed in the introduction, the theory of dislocation splitting the stress in two distinct parts is used to define the thermoviscoplastic behavior of materials. The equivalent stress for body cubic-centered (bcc) material is defined as τ = τ µ + τ * (with τ µ corresponding to the athermal stress and τ * the thermal stress), assuming an additive contribution of each part. Each component corresponds to a barrier preventing the dislocation motion. Two main categories are distinguished: long-range and short-range obstacles. The stress required to overcome the long-range barrier is represented by the athermal stress τ µ and considers the obstacles such as grain boundaries or other dislocations. Therefore, this stress is related to the dislocation density into the considered material. For bcc metals, τ µ is weakly affected by the strain rate or the temperature. The other component called the thermal stress τ * , is representing the stress required to overcome the short-range obstacles. These obstacles correspond to the easiness for a dislocation to move from one site to the next, at an atomic level. At this scale, the motion of the dislocations requires to overcome the Peierls-Nabarro stress [23]. This stress corresponds to that required for the dislocation to move from one equilibrium position to the next. The energy required to overcome this barrier is temperature and strain rate dependent. On one hand, the temperature increase reduces the required energy by increasing the amplitude of atom vibrations which in turn facilitates the jump of individual atoms. On the other hand, higher strain rate induces a higher dislocation velocity and then a reduced time for the dislocations to move from one position to another. Therefore, the contribution of the thermal energy to overcome this barrier is weaker. These mechanisms are called thermally activated mechanisms. At T = 0 K, there is no contribution of the thermal energy, and the stress corresponding to the thermally activated mechanisms is maximum. This stress is called the mechanical threshold stress. This theoretical change of the equivalent flow stress with the temperature is represented in Figure 7 [26]. Below a critical temperature T c , the flow stress corresponds to the addition of the thermal stress τ * and the athermal stress, τ µ . For temperature higher than T c , the thermal stress vanishes and the flow stress is equal to the athermal stress. Moreover, beyond a value of strain rate, dislocation motions become continuous. Their motion is then controlled by their interactions with phonons and electrons and corresponds to the phenomena of viscous drag. In these conditions, stress opposing to the dislocation motion is directly proportional to their velocity [27]. Nevertheless, the performed experiments did not reach strain rates high enough to observe such phenomenon.
The following section describes the constitutive relations studied in this work and the procedure used to define their constants based on previous experiments.

Constitutive Models
In this section, several models are compared with experiments to estimate their limits in term of strain rate and temperature. In this study, four models are analyzed. The first one is the well-known and widely used phenomenological Johnson-Cook model. The other models are defined as semi-physical since they are originally based on the dislocation theory discussed previously. The models described in this paper are the original versions proposed by Voyiadjis and Abed [15] and the Rusinek-Klepaczko models [1]. Moreover, two versions are proposed to update the model of Voyiadjis and Abed to describe in a more precise way the observed temperature and strain rate sensitivities.

Johnson-Cook Constitutive Model
The Johnson-Cook constitutive relation is a phenomenological model allowing one to describe the mechanical behavior of materials, Equation (5) [4]. It consists of splitting the different effects as the hardening, the temperature sensitivity, and the strain rate sensitivity in three multiplicative terms. Therefore, the first term represents the strain hardening, while the strain rate sensitivity and the thermal effect are introduced using respectively the second and third terms as shown below where A, B, C, n and m are the model parameters and T * is the homologous temperature. The parameters A, B and n represent respectively the yield stress, the modulus of plasticity and the hardening coefficient. The constant parameters are determined by fitting the true stress vs. true plastic strain curve for an assumed strain rateε 0 and temperature T 0 . The parameter C, representing the strain rate sensitivity, is determined by considering the stress level obtained for an imposed plastic strain level and different strain rates as shown in Figure 8a. The same procedure is applied to determine the parameter m, which represents the temperature sensitivity, as shown in Figure 8b. The values of the parameters for the studied HSS material are reported in Table 2. Considering this Johnson-Cook (JC) model, the parameters are easy to identify but its validity is limited. The term related to the strain rate assumes a linear sensitivity. This statement is valid for either the quasi-static or the dynamic observations, but does not fit the experimental data over the whole studied range. Moreover, a poor description of the strain rate sensitivity induces mistakes when instabilities are considered [28].

Voyiadjis-Abed Constitutive Model
The semi-physical model proposed by Voyiadjis and Abed [15] is based on the dislocations theory and the split of the equivalent stress into an athermal stress σ a and a thermal stress σ th . The authors developed four specific constitutive relations, for body-centered cubic (bcc), face-centered cubic (fcc), hexagonal close-packed (hcp), and stainless steel microstructure. This study deals with the version developed for bcc metals. For most of these metals, the hardening is weakly influenced by the strain rate and the temperature [14]. Taking into account this consideration, the athermal stress is defined only as a function of the plastic strain. On the contrary, the thermal stress does not depend on the equivalent plastic strain, but depends on the strain rate and the temperature. Therefore, the constitutive relation is defined as follows, in Equation (6) as where the previously discussed athermal stress σ a , is defined using the first part of the Johnson-Cook constitutive model, Equation (7) σ a ε p = Y a + Bε p n where Y a represents the yield stress, B is the modulus of plasticity, and n is the hardening coefficient. The thermal part σ th is defined by Equation (8). More details may be found in the work of Voyiadjis and Abed [15].
where σ represents the mechanical threshold stress, β 1 and β 2 are constants related to the dislocation density and p and q are constants related to the obstacles shape that dislocations have to overcome [29]. Thus, eight parameters have to be defined to model the behavior of the considered HSS material. This constitutive relation assumes that the contribution of thermally activated mechanisms vanishes beyond a certain temperature. This assumption is used to identify the first three parameters Y a , B, and n. As the stress observed for various temperatures do not evolve between 423 K and 473 K, it is assumed that for these temperatures, there is only a contribution of the athermal stress σ a . Therefore, the results obtained at 473 K are used to identify Y a , B, and n. The parameters p and q take the typical values as: p ∈ [0,1] and q ∈ [1,2]. In this work, 0.5 and 1.5 have been chosen for p and q, respectively as proposed by [29]. The mechanical threshold stress σ, then corresponds to the value added to the athermal stress σ a in order to obtain the flow stress at 0 K. Consequently, for a fixedε p , the term σ − Y a − Bε n p p is plotted as a function of T 1/q and approximated using a linear function, Figure 9. At T = 0 K, the following equality is reached σ − Y a − Bε n p p = σ p , allowing to determine σ p , as shown in Figure 9. In order to determine β 1 and β 2 , the term 1 − σ − Y a − Bε n p / σ p q versus ln ε p is plotted for a fixed temperature. The results are reported in Figure 10a. The linear approximation is supposed to be equal to β 1 T room − β 2 T room ln ε p with T room assumed as equal to 293 K. Nevertheless, the experimental data do not exhibit a linear behavior, especially under dynamic loading conditions. Therefore, all the domains defined in terms of strain rate may not be correctly modeled, as shown in Figure 10a. To obtain a better correlation with experimental observations a nonlinear approach is proposed, allowing one to model the non-linearity of the stress at high strain rates and to increase the strain rate sensitivity, as shown in Figure 10b. The thermal part, σ th is defined in this case using Equation (9) as The logarithm of the strain rate has been deleted as proposed in the modified JC model to avoid the stress linearity behavior with the strain rate [23][24][25]. However, using the second approach (Equation (9)) the quasi-static part is not well described. The solution to fit correctly the studied domain and to be in agreement with experiments is to combine the two previously detailed analytical relations Equations (8) and (9). This implies to define two domains corresponding to the quasi-static and dynamic regimes. Thus, the relation describing the strain rate sensitivity is composed of two parts. To ensure the continuity of the strain rate sensitivity definition between the two domains, it is necessary to impose the following conditions, Equations (10) and (11) β Static The results are reported in Figure 11. Thus, the strain rate transition is defined by only the two coefficients β 1 and β 2 . Using the HSS material, the transition strain rate to define the continuity of the strain rate sensitivity is equal to 403 s −1 . Figure 11. Determination of β 1 and β 2 using mixed approach and Equation (10) to describe the continuity condition. Therefore, coupling Equation (8) and Equation (9), it is possible to fit experimental results performed at room temperature for a large range of strain rates with an improved accuracy. It has to be noted that a strain rate transition in order of 10 2 s −1 is frequently observed for metals to describe the change between isothermal and adiabatic conditions. It induces a quick increase of the stress level with the strain rate [23][24][25]. Except for β 1 and β 2 , these three versions of the constitutive relations developed by Voyiadjis and Abed have the same parameters, as reported in Table 3.  The values for β 1 and β 2 for each approach are reported in Table 4.

Rusinek-Klepaczko Constitutive Model
The Rusinek-Klepaczko semi-physical model reported in [1] is initially based on the theory of dislocations and the model proposed by Klepaczko [30]. This study deals with the original version reported in [1], although several authors [31][32][33] have proposed an extended version of this model. It decomposes the equivalent stress in an internal stress σ µ , related to long-range obstacles and an effective stress σ * , related to short range obstacles, as in Equation (12) σ ε p ,ε p , T = E(T) E 0 σ µ ε p ,ε p , T + σ * ε p , T where E(T) and E 0 represent respectively the Young's modulus depending on the temperature and at T = 0 K. Although this value can vary, the evolution of the Young's modulus for steel [18] is assumed to be low enough (between 173 K and 473 K) in order to assume it as constant. The internal stress σ µ is expressed by Equation (13) σ µ ε p ,ε p , T = B ε p , T (ε + ε 0 ) n(ε p ,T) (13) where ε 0 is the strain level defining the yield stress, B and n are strain rate and temperature dependent coefficients and are defined as follows, in Equations (14) and (15) Finally, the effective stress σ * is given by Equation (16) σ where B 0 , ν, n 0 , D 2 , σ * 0 , D 1 , and m * are the model parameters,ε max andε min are the strain rate limits of the model and T m is the melting temperature. The Macaulay brackets are defined as: x = x if x > 0 and x = 0 otherwise. The internal stress σ µ mimics the evolution of the dislocation density (creation and annihilation) which is affected by the strain (hardening), the strain rate and the temperature [34]. On the contrary, the effective stress σ * describes the coupling and the reciprocity between the strain rate and the temperature. From a physical point of view, it represents the difficulty of a dislocation to overcome an obstacle. This additional stress increases with the strain rate but tends to vanish as the temperature increases [35]. Therefore, the model neglects the effective stress at low strain rates at room and at high temperatures.
The first step consists in defining D 1 , assuming that the effective stress σ * is equal to zero at low strain rate (10 −3 s −1 in the present case) and 293 K as suggested by Equation (17) Assuming that the effective stress is equal to zero at low strain rate, it allows to determine a first approximation of B(0.001,293) and n(0.001,293). As the effective stress σ * corresponds to an instantaneous effect and the internal stress σ µ to a history effect, the next step is to estimate the stress increase with strain rate for an imposed plastic strain level. Thus, Equation (18) allows one to determine the values of σ * 0 and m * .
Once the effective stress σ * is defined, the values of B ε p , T and n ε p , T for different sets of conditions can be identified. Finally, Equation (14) and (15) are used to determine B 0 , ν, n 0 , and D 2 . Parameter values determined for this model are given in Table 5.   Figure 12a for different strain rates. As the strain rate varies during dynamic experiments, an average value has been calculated. It can be seen that hardening is well represented by the model in quasi-static loading, emphasizing that the first term of Equation (5) provides a good description of the material behavior in these conditions. Nevertheless, in the dynamic regime, one can observe that this constitutive relation does not fit the experimental data. It overestimates the stress at the lowest strain rates and underestimates at the highest. The Johnson-Cook model assumes a linear strain rate sensitivity of the stress, which is in contrast with experiments ( Figure 8a). The comparison between the model and experiments is shown in Figure 12b for different temperatures. Experiments at various temperatures were performed at an average strain rate of about 1100 s −1 . This strain rate induces an overestimation of the flow stress due to the poor description of the strain rate sensitivity. Nevertheless, this model provides a correct description of the temperature sensitivity.

Comparison of Experiments with the Voyiadjis-Abed Model using Different Approaches
The comparison between constitutive relations and experiments are shown in Figures 13-15 for the original linear, the modified non-linear and the mixed approaches. The strain rate sensitivity is not be well represented for a dynamic loading in the original version, Figure 13a. This is due to the linear sensitivity assumed by this model, as expressed in Equation (8). The flow stress value described by this model at 7912 s −1 is~50 MPa higher than the description at 586 s −1 , although experiments showed an increase of~200 MPa. The non-linear modification brought to this model provides an overestimation of the strain rate sensitivity, as observed in Figure 14a. The flow stress at the highest reported strain rate (7912 s −1 ) is overestimated. In dynamic conditions, the mixed approach uses the same equation of the thermal part of the flow stress, σ th , defined by Equation (9). Therefore, the results obtained using the mixed approach are better than those obtained with previous versions, as shown in Figure 15a. The experimental values of the dynamic compression tests at various temperatures are well represented by these models, including at low temperature. The related comparisons are shown in Figures 13b, 14b and 15b. Regarding the temperature sensitivity at low strain rate, shown in Figures 13c, 14c and 15c, a correct description of the experiments is obtained at room temperature and for T 0 = 373 K. However, the flow stress is overestimated at 423 K for all of the approaches.

Comparison of Experiments with Rusinek-Klepaczko Model
A comparison, at different strain rates, between experiments and the model is shown in Figure 16a. The representation of the dynamic regime is similar to the one obtained with the modified Voyiadjis-Abed constitutive relation. Although the sensitivity in the quasi-static case is too low in comparison with experiment performed at 0.1 s −1 , this model is in good agreement with experimental measurements in most of the studied cases. Regarding experiments at various temperatures, this model provides a good description of the temperature sensitivity and the thermal softening (see Figure 16b). Figure 16c provides a comparison of the experimental results and the constitutive relation for the quasi-static compression tests at high temperature. The temperature sensitivity described by this model at 0.001 s −1 is lower than the one observed experimentally.

Comparison between Constitutive Relations and Limits
A comparison of the constitutive relations described previously is reported in Figure 17. It compares the macroscopic strain rate sensitivity described by the different models. It is observed that the Johnson-Cook model and the model defined originally by Voyiadjis and Abed do not allow to define the nonlinear strain rate sensitivity observed experimentally for metals as reported in [23][24][25]. To overcome this limitation, the model proposed by Rusinek and Klepaczko and the constitutive relation proposed by Voyiadjis and Abed (modified and using mixed approach) may be used. The three models are sensibly giving the same trends, Figure 17. The modified Voyiadjis-Abed model presents a strain rate sensitivity at ε p = 0.07 which is too low for low strain rates, and too high for high strain rates. The expression of the thermal part, σ th (Equation (9)) was chosen to represent the strong sensitivity observed for dynamic loadings. This explains the poor description obtained for quasi-static cases. Moreover, the parameters identification takes into account experiments at low strain rates and therefore increases the inaccuracies for dynamic loadings. The mixed approach used provides a compromise between the linear sensitivity observed at low strain rates and the nonlinear sensitivity at high strain rates. The Rusinek-Klepaczko constitutive relation provides a better representation of high strain rates experiments, but underestimates the sensitivity for quasi-static cases. It has to be noted that the values of Figure 17 represent the flow stress for fixed strain level. In order to compare the efficiency of the models, the mean error expressed in Equation (19) was calculated for several cases for each constitutive relation.

Mean error
The comparison of errors obtained for various strain rates at room temperature for each constitutive relation is shown in Figure 18a. For a strain rate value ofε = 0.001 s −1 , the Johnson-Cook constitutive relation provides a better description than the others models. This is explained by the fact that this loading conditions (0.001 s −1 , 293 K) was taken as a reference to identify the parameters of this model. The errors obtained are globally higher for high strain rates. The mixed approach of the Voyiadjis-Abed model and the Rusinek-Klepaczko constitutive relation present best agreement for dynamic cases. The error comparisons obtained for various temperatures at~1100 s −1 for each constitutive relation is shown in Figure 18b. Regarding the Johnson-Cook model, no comparison can be established between this model and experiments below the reference temperature used during the parameters identification (293 K). Other models provide a good description of the temperature sensitivity for dynamic loading. The gaps obtained for various temperatures at 0.001 s −1 for each constitutive relation is shown in Figure 18c. As shown in Figure 6, the thermal sensitivity of this material is strain rate dependent. The Johnson-Cook model assumes a constant thermal sensitivity regardless of the strain rate. Moreover, the Rusinek-Klepaczko underestimates this sensitivity at 0.001 s −1 . This explains higher errors observed in quasi-static conditions compared to the other constitutive relations. The other models provide good description of the temperature sensitivity for quasi-static loading. For every model, the highest error is observed at the highest studied temperature, in quasi-static conditions. As a consequence, when studying the case of impact loading, the best compromise can be found with the mixed approach of the Voyiadjis-Abed constitutive relation due to its ability to describe high strain rate phenomena without degrading low strain rate representation.

Conclusions
The mechanical behavior of a high strength steel was studied, with experiments performed over a wide range of strain rates varying from 10 −3 -10 4 s −1 and various temperatures starting from 173 K to 473 K. This experimental investigation allowed to highlight the different sensitivities of this high strength steel. The parameters of five different constitutive relations were identified and the results were compared to the experiments. The first one is the Johnson-Cook constitutive model. This model provides a good description of the temperature sensitivity, but only for temperatures above the reference temperature (293 K in this work). Furthermore, it assumes a linear strain rate sensitivity, which is not the case of the studied material, and cannot represent any coupling effect between the strain rate and the temperature. Nevertheless, among the considered relation, its parameters are the easiest to determine and it allows a correct description of the material mechanical behavior at conditions close to those chosen as references. In order to overcome these limitations, two semi-physical models are considered. The constitutive relation developed by Voyiadjis and Abed gives a better description of the thermal sensitivities over a wide range of strain rates. Nevertheless, the strain rate sensitivity is not well fitted, especially at high strain rates. A modified nonlinear version for dynamic cases has been suggested. Nevertheless, this modification degrades the representation of low strain rates. Therefore, a mixed approach using a transition between low and high strain rates has been suggested. This model provides a better description at high strain rates and is more consistent with the observed trend. Moreover, it offers a good correlation with tests at various temperatures. The last studied model has been developed by Rusinek and Klepaczko. This constitutive relation allows description of the sensitivities observed experimentally but presents higher mean errors than the mixed approach of the Voyiadjis and Abed model for temperature sensitivity description. It should be noted that the modified and the mixed versions of the Voyiadjis-Abed model have a maximal strain rate limit at~10 4 s −1 . From a physical point of view, other phenomena, such as viscous drag, become preponderant beyond this value and are not taken into account in these models. Regarding the Rusinek-Klepaczko model, an extension has been studied in [33] to introduce this effect.