An Enhanced Lemaitre Model and Fracture Map for Cr5 Alloy Steel during High-Temperature Forming Process

To effectively control and predict crack defects in the high-temperature forming process of Cr5 alloy steel, based on the traditional Lemaitre damage model, a new high-temperature damage model of Cr5 alloy steel was proposed which considered the change of material elastic modulus with temperature, the influence of material hydrostatic pressure as well as temperature and strain rate on material damage. Because Cr5 alloy steels are usually forged at high temperatures, tensile testing is an important method to study the damage behaviour of materials. Through the high-temperature tensile test and elastic modulus measurement test of the Cr5 alloy steel, the stress–strain curves and the relationship curves of the elastic modulus value with the temperature of Cr5 alloy steel under different temperatures and strain rates were obtained. A new high-temperature damage model of Cr5 alloy steel was built by introducing the Zener–Hollomon coefficient considering the influence of temperature and strain rate. The established high-temperature damage model was embedded in Forge® finite element software through the program’s secondary development method to numerically simulate the experimental process of Cr5 alloy steel. Comparing the difference between the displacement–load curves of the numerical simulation and the actual test of the tensile process of the experimental samples, the correlation coefficient R2 is 0.987 and the difference between the experimental value and the simulated value of the tensile sample elongation at break is 1.28%. The accuracy of the high-temperature damage model of Cr5 alloy steel established in this paper was verified. Finally, the high-temperature damage map of Cr5 alloy steel was constructed to analyse the variation law of various damage parameters with the temperature and strain rate of the high-temperature damage model of Cr5 alloy steel.


Introduction
As the key components of the rolling mill of the steel rolling production line, with the continuous development of the steel rolling industry, the requirements for the hardness and wear resistance of the rolls in industrial production are getting higher and higher. Cr5 alloy steel has been widely used in the production of modern rolls as a substitute for 60CrMnMo, 50CrNiMo, 70Cr3NiMo and Cr4 [1] due to its good hardness, wear resistance and resistance. Fatigue performance has always been one of the key research topics for roll steels in the iron and steel production industry. Many scholars from all over the world have conducted a lot of research into the mechanical properties and fatigue resistance of Cr5 alloy steel and have achieved certain results. Junkai Fan et al. [2] determined the forming properties of Cr5 alloy steel by studying the mechanical properties of Cr5 alloy steel. Qing Gu et al. [3] determined the mechanical properties of Cr5 alloy steel by studying the fatigue behaviour of Cr5 alloy steel. M. Surendran [4] and others determined the fatigue characteristics of Cr5 alloy steel by studying the hoop stress of Cr5 alloy steel. Many scholars conducted steel. Ashwani Verma et al. [23] simulated the tensile process of the dual-phase steel based on the Lemaitre damage model, and obtained the damage parameters of the dual-phase steel in the tensile state. Based on the Lemaitre damage model, the fatigue life of the alloy was studied and predicted by L.M. Araújo et al. [24]. Kazem Malekipour et al. [25] observed and simulated the damage evolution process of AISI 316L parts based on the Lemaitre damage model. The axial stretching process was simulated, and the damage evolution law of the thin-walled steel pipe was obtained by Sheng He et al. [26]. Many scholars have used the traditional Lemaitre damage model to accurately predict the damage behaviour of different forming processes. Tandon Puneet et al. [27] simulated the establishment of progressive sheet metal forming process based on the Lemaitre damage model and investigated the effectiveness of the Lemaitre damage model in the simulation of the ISF process. However, this model can only be used for the calculation and prediction of damage defects in the cold forming process of metal materials. Because this model does not affect the crack initiation and propagation process of temperature and strain rate during the high-temperature forming process, it cannot be used for the calculation and prediction of damage defects in the high-temperature forming process; On the other hand, this damage model lacks consideration of the effect of temperature-induced elastic modulus changes on material damage [28].
The traditional Lemaitre damage model does not consider the effect of temperature and strain rate on damage, and it is difficult to predict damage behaviour at high temperature. In this paper, considering the influence factors of temperature and strain rate, the hightemperature Lemaitre damage model of Cr5 alloy steel is improved and established, and its high temperature damage map is drawn to accurately predict the damage initiation of Cr5 alloy steel at high temperature.
In order to study the damage law of Cr5 alloy steel at high temperature, the hot tensile test and high temperature elastic modulus measurement test of Cr5 alloy steel were carried out. Based on the test data, a high-temperature Lemaitre damage model of Cr5 alloy steel considering the temperature and strain rate was established to predict its high-temperature damage behaviour. The high-temperature damage model of Cr5 steel was embedded in the finite element software FORGE ® to simulate the tensile test of the Cr5 alloy steel. The accuracy of the established high-temperature damage model was verified by comparing the displacement-loads and fracture lengths of the simulation and the test. In order to avoid the crack defect of the Cr5 roll during hot forming, the corresponding Lemaitre high-temperature damage model damage map was proposed.

Materials and Test Procedures
The material used in the test is Cr5 alloy steel, which is used in conventional industrial production and its chemical compositions are shown in Table 1. The equipment used in the test is a Gleeble-1500D thermal simulation testing machine, as shown in Figure 1. The machine has the function of the real-time calculation of strain, meaning that the strain ε equals ln(L/L0) in real time. In a certain time interval delta t, the strain changed delta ε, then under the linear control of the strain, a strain rate equal to (delta ε)/(delta t) and equal to a constant can be achieved. The central part of the tensile specimen is processed by means of arc transition to ensure that the tensile specimen can be smoothly broken from the central part [29]. The samples have a taper that gradually decreases from the end to the centre, and both ends are stretched with M10 standard thread. The specific shape and size of the tensile specimen are shown in Figure 2, which was prepared in accordance with ASTM E8/E8M-2016a.    Because 800 °C-1150 °C is the common forging temperature range of the Cr5 roll, and the strain rate of the material ranges between 0.01 s −1 and 5 s −1 in actual production, the temperatures used in the test are 800 °C, 850 °C, 900 °C, 950 °C, 1000 °C, 1050 °C, 1100 °C, and 1150 °C, and the real strain rates are set as 0.01 s −1 , 0.1 s −1 , 1 s −1 , 5 s −1 by Gleeble-1500D thermal simulation testing machine, and a tensile test is performed at each temperature and strain rate, and the total number of samples consumed are 32. Using thermocouple wire to measure the temperature of the tensile sample, the tensile specimen was first heated to the specified temperature at a speed of 10 °C/s, and then kept warm for 3 min to eliminate the temperature gradient. Finally, it remained tensile at a predetermined strain rate until fracture. Figure 3 shows the displacement-load curves of the Cr5 alloy steel obtained by tensile tests. Table 2 shows the fracture displacement of the Cr5 alloy steel at different temperatures and strain rates. It can be seen from the figure that at the same strain rate, with the continuous increase in temperature, the fracture displacement of Cr5 alloy steel gradually increases. When the rate is 1 s −1 , the fracture displacement of Cr5 alloy steel is 13.89 mm at 800 °C, and the fracture displacement of Cr5 alloy steel is 17.93 mm at 1150 °C. For example, when the temperature is 1000 °C, the fracture displacement of Cr5 alloy steel is 12.83 mm when the strain rate is 0.01 s −1 , while the fracture displacement of Cr5 alloy steel is 14.69 mm when the strain rate is 5 s −1 . At the same strain rate, with the continuous increase in temperature, the maximum load of Cr5 alloy steel during stretching gradually   Because 800 °C-1150 °C is the common forging temperature range of the Cr5 roll, and the strain rate of the material ranges between 0.01 s −1 and 5 s −1 in actual production, the temperatures used in the test are 800 °C, 850 °C, 900 °C, 950 °C, 1000 °C, 1050 °C, 1100 °C, and 1150 °C, and the real strain rates are set as 0.01 s −1 , 0.1 s −1 , 1 s −1 , 5 s −1 by Gleeble-1500D thermal simulation testing machine, and a tensile test is performed at each temperature and strain rate, and the total number of samples consumed are 32. Using thermocouple wire to measure the temperature of the tensile sample, the tensile specimen was first heated to the specified temperature at a speed of 10 °C/s, and then kept warm for 3 min to eliminate the temperature gradient. Finally, it remained tensile at a predetermined strain rate until fracture. Figure 3 shows the displacement-load curves of the Cr5 alloy steel obtained by tensile tests. Table 2 shows the fracture displacement of the Cr5 alloy steel at different temperatures and strain rates. It can be seen from the figure that at the same strain rate, with the continuous increase in temperature, the fracture displacement of Cr5 alloy steel gradually increases. When the rate is 1 s −1 , the fracture displacement of Cr5 alloy steel is 13.89 mm at 800 °C, and the fracture displacement of Cr5 alloy steel is 17.93 mm at 1150 °C. For example, when the temperature is 1000 °C, the fracture displacement of Cr5 alloy steel is 12.83 mm when the strain rate is 0.01 s −1 , while the fracture displacement of Cr5 alloy steel is 14.69 mm when the strain rate is 5 s −1 . At the same strain rate, with the continuous increase in temperature, the maximum load of Cr5 alloy steel during stretching gradually Because 800 • C-1150 • C is the common forging temperature range of the Cr5 roll, and the strain rate of the material ranges between 0.01 s −1 and 5 s −1 in actual production, the temperatures used in the test are 800 • C, 850 • C, 900 • C, 950 • C, 1000 • C, 1050 • C, 1100 • C, and 1150 • C, and the real strain rates are set as 0.01 s −1 , 0.1 s −1 , 1 s −1 , 5 s −1 by Gleeble-1500D thermal simulation testing machine, and a tensile test is performed at each temperature and strain rate, and the total number of samples consumed are 32. Using thermocouple wire to measure the temperature of the tensile sample, the tensile specimen was first heated to the specified temperature at a speed of 10 • C/s, and then kept warm for 3 min to eliminate the temperature gradient. Finally, it remained tensile at a predetermined strain rate until fracture. Figure 3 shows the displacement-load curves of the Cr5 alloy steel obtained by tensile tests. Table 2 shows the fracture displacement of the Cr5 alloy steel at different temperatures and strain rates. It can be seen from the figure that at the same strain rate, with the continuous increase in temperature, the fracture displacement of Cr5 alloy steel gradually increases. When the rate is 1 s −1 , the fracture displacement of Cr5 alloy steel is 13.89 mm at 800 • C, and the fracture displacement of Cr5 alloy steel is 17.93 mm at 1150 • C. For example, when the temperature is 1000 • C, the fracture displacement of Cr5 alloy steel is 12.83 mm when the strain rate is 0.01 s −1 , while the fracture displacement of Cr5 alloy steel is 14.69 mm when the strain rate is 5 s −1 . At the same strain rate, with the continuous increase in temperature, the maximum load of Cr5 alloy steel during stretching gradually decreases. For example, when the strain rate is 5 s −1 , the maximum load of Cr5 alloy steel when stretched at 800 • C is 1707.44 N and the maximum load of Cr5 alloy steel in tension is 571.53 N at 1150 • C. At the same temperature, with the gradual increase in strain rate, the maximum load of the Cr5 alloy steel during stretching gradually increases. Additionally, when the strain rate is 5 s −1 , the load is 250.52 N, and the maximum load when the Cr5 alloy steel is stretched is 571.53 N. The experimental results are consistent with the literature [30]. It can be obtained that the maximum load of Cr5 alloy steel material can be reached when plastic deformation occurs and decreases with the increase in temperature, and its maximum load increases with the increase in strain rate. Additionally, the effect of temperature on the maximum load of Cr5 alloy steel when plastic deformation occurs is greater than the effect of the strain rate. This is because, with the continuous increase in temperature, the Cr5 alloy steel material will recover and recrystallize at high-temperature, and the maximum load required for plastic deformation will gradually decrease due to the softening of the material. However, with the increase in strain rate, the maximum load required for the plastic deformation of Cr5 alloy steel will gradually increase with the increase in the strain rate due to the influence of the deformation degree per unit time [31].

Analysis of Test Results
the maximum load of the Cr5 alloy steel during stretching gradually increases. Ad ally, when the strain rate is 5 s −1 , the load is 250.52 N, and the maximum load wh Cr5 alloy steel is stretched is 571.53 N. The experimental results are consistent w literature [30]. It can be obtained that the maximum load of Cr5 alloy steel material reached when plastic deformation occurs and decreases with the increase in tempe and its maximum load increases with the increase in strain rate. Additionally, th of temperature on the maximum load of Cr5 alloy steel when plastic deformation is greater than the effect of the strain rate. This is because, with the continuous incr temperature, the Cr5 alloy steel material will recover and recrystallize at high-te ture, and the maximum load required for plastic deformation will gradually decrea to the softening of the material. However, with the increase in strain rate, the ma load required for the plastic deformation of Cr5 alloy steel will gradually increa the increase in the strain rate due to the influence of the deformation degree per un [31].
Therefore, in the actual production of Cr5 alloy steel products and the formula the hot forming process plan, in order to effectively control crack defects, the highe perature and a lowest strain rate possible should be selected.   Therefore, in the actual production of Cr5 alloy steel products and the formulation of the hot forming process plan, in order to effectively control crack defects, the highest temperature and a lowest strain rate possible should be selected.

Traditional Lemaitre Damage Model
Damage is a very important factor affecting the quality of metal products. However, studying the damage of a specific material is a very complicated process, which is closely related to the stress state of the material during plastic processing [32]. Therefore, it is of great significance to determine a more reliable material damage model to effectively predict and control the product crack defects. The Lemaitre damage model is a damage model that clearly expresses the damage of materials and considers more parameters and variables, and can be directly applied to guide the production and manufacture of products in industrial production. The traditional Lemaitre damage model considers the material to be damaged or destroyed when the absolute damage value D suffered by the material during processing is greater than or equal to the critical damage value D c of the material itself, and its mathematical formula is expressed as follows [33]: where dD represents the absolute damage value increment; f (ξ) represents the stress triaxiality function; ξ represents the stress triaxiality; E represents the elastic modulus of the material; S represents the anti-damage factor; ε represents the equivalent plastic strain; dε represents the equivalent plastic strain increment; v is the Poisson's ratio of the material; σ m represents the mean stress; σ represents the equivalent stress experienced by the material. However, in actual industrial production, the stress state experienced by materials is often not in a single direction, but in multiple axial stress states. Therefore, the expression of D c in the Lemaitre damage model under multiaxial stress states is as shown in Formula (3): where σ u and D 1c represent the stress limit value and the damage critical value of the material in the tensile test, respectively. The damage variable D represents the damage degree of the material in the unit. When the value of the damage variable D is 0, there is no damage inside the material. When the damage variable of the material reaches the critical damage value D c , this means that the material is damaged. In the Lemaitre damage model, the relative damage value D rel , that is, the ratio of the absolute damage value D to the critical damage value D c , is used to represent the degree of damage inside the material [34], and its specific expression is as follows: All four damage parameters D 1c , σ u , S and ε th in the Lemaitre damage model can be obtained through the tensile test of the material, and the damage parameters at different temperatures and different strain rates that can be obtained through the actual tensile test data according to Formulas (5)-(8) are calculated: where ε m is the tensile strain of the material; ε B is the fracture strain of the material; R m is the tensile strength of the material; R B is the breaking strength of the material.

Modification of Lemaitre Damage Model
Without considering the effect of temperature, the traditional Lemaitre damage model cannot be applied to the high-temperature forming process of the material, and the stress triaxiality term in the damage model is a square term, which is not able to distinguish between the influence of positive and negative hydrostatic stress on the material damage. Therefore, improvements and corrections are required. The Zener-Hollomon equation is a function related to temperature and strain rate, which can describe the rheological behaviour of materials at different temperatures and strain rates. Its specific expression is shown in Equation (9) [35]: where . ε represents the strain rate of the material; Q represents the deformation activation energy of the material; R is the gas constant; T is the deformation temperature of the material.
In previous experimental studies [36,37], it was found that compressive stress often plays a role in the healing of micro-defects and micro-cracks in the material in the process of plastic deformation, and there is a certain difference with the role of tensile stress in the process of material damage. The traditional Lemaitre damage model does not fully consider the different effects of tensile and compressive stress on the damage of the material. Thus, the damage model often cannot accurately describe the damage of the material. Therefore, the influence of the internal compressive stress of the material on the damage of the material (which will promote crack healing) can be removed through coupling, and the Lemaitre damage model of the Cr5 alloy steel with only tensile stress is obtained. The traditional Lemaitre damage model is modified to obtain the modified Lemaitre damage model [38].
The absolute damage value increment equation of the modified Lemaitre damage model is: In order to obtain the Lemaitre damage model of Cr5 alloy steel suitable for different temperatures and different strain rates, the Zener-Hollomon coefficient was coupled with the traditional Lemaitre damage model, and the damage parameter values of the Cr5 alloy steel at different temperatures and different strain rates were obtained. As such, a Lemaitre damage model of the Cr5 alloy steel which is suitable for different temperatures and different strain rates can be built, the Lemaitre damage model can be coupled using the Zener-Hollomon coefficient, and it can be used to manufacture Cr5 alloy steel products at different temperatures and strain rates in industrial production [39].
In summary, the coupled form of the damage parameters in the improved Lemaitre damage model based on the Zener-Hollomon coefficient are: As a model that can predict the damage of materials, the Lemaitre damage model is closely related to various parameters of the material. As a discipline that describes the damage of materials, damage mechanics are related to the properties of all aspects of materials. As a physical quantity describing the elastic properties of materials, the elastic modulus E is a very fundamental parameter variable in the Lemaitre damage model. The value of the elastic modulus E used in the calculation of the previous Lemaitre damage model is often a constant, and the elastic modulus E value of the material at different temperatures is not measured. The Gleeble-1500D thermal simulation testing machine used in this paper used C-Gauge sensor to accurately measure the elastic modulus E values of Cr5 alloy steel at different temperatures, and these values can be substituted into the Lemaitre damage model so as to obtain a more accurate Lemaitre damage model of the Cr5 alloy steel at different temperatures and different strains.

Measurement of Elastic Modulus E Value
As one of the essential performance parameters of engineering materials, the modulus of elasticity's value will fluctuate with different material composition and structure. Especially as the temperature of the material increases, the distance between atoms inside the material will gradually increase and then the bond between atoms will be weakened, resulting in a decreased E value, higher temperature, and greater changes in the value of the elastic modulus. The elastic modulus E values of the same material at room temperature and at high temperature are often quite different. Therefore, in the present paper, the study of the damage of a certain material at high-temperature was performed and the modulus of elasticity is a crucial factor.
The elastic modulus is a physical quantity that describes the elastic properties of materials. The value of the elastic modulus E of the material is quite different under different conditions. Especially under the influence of material temperature, it often changes several times or even dozens of times. It has a great influence on establishing an accurate high-temperature damage model. For Cr5 alloy steel, the elastic modulus value at room temperature is 206 GPa. However, with the temperature of the material rising to a higher temperature, its elastic modulus value sharply decreases. Therefore, this has great influence on the establishment of an accurate high-temperature Lemaitre damage model for Cr5 alloy steel [40]. Previous studies on the Lemaitre damage model have rarely involved attention to the change of material elastic modulus, the research results obtained are often quite different from the actual results and are of low accuracy in terms of material damage prediction. Therefore, in order to improve the accuracy of the high-temperature Lemaitre damage model of Cr5 alloy steel, it is necessary to measure the actual elastic modulus E value of Cr5 alloy steel at different temperatures.
In this study, the elastic modulus E values of Cr5 alloy steel at different temperatures were measured. The measurement tests were performed on a Gleeble-1500D thermal simulation testing machine; the elastic modulus of Cr5 alloy steel was measured at different temperatures for one time each; the samples were heated at 10 • C /s and kept for 3 min to reach the specified temperature; the elastic modulus was then measured at the strain rate of 0.01 s −1 ; and the elastic modulus E values of Cr5 alloy steel at different temperatures were measured with the help of the C-Gauge sensor of the DSI company in the United States. The specific measurement results are shown in Table 3, whilst the fitting curve of the elastic modulus E value of Cr5 alloy steel and temperature is shown in Figure 4. It can be seen from Figure 5 that the change curve of the elastic modulus E value of Cr5 alloy steel in the high-temperature stage with temperature is approximately a linear function. Therefore, it can be fitted using a linear function. The relationship curve of the elastic modulus of the Cr5 alloy steel at the high-temperature stage with temperature is obtained as follows: E = −0.15T + 176.89 (17) ulation testing machine; the elastic modulus of Cr5 alloy steel was measured at different temperatures for one time each; the samples were heated at 10 °C /s and kept for 3 min to reach the specified temperature; the elastic modulus was then measured at the strain rate of 0.01 s −1 ; and the elastic modulus E values of Cr5 alloy steel at different temperatures were measured with the help of the C-Gauge sensor of the DSI company in the United States. The specific measurement results are shown in Table 3, whilst the fitting curve of the elastic modulus E value of Cr5 alloy steel and temperature is shown in Figure 4.  It can be seen from Figure 5 that the change curve of the elastic modulus E value of Cr5 alloy steel in the high-temperature stage with temperature is approximately a linear function. Therefore, it can be fitted using a linear function. The relationship curve of the elastic modulus of the Cr5 alloy steel at the high-temperature stage with temperature is obtained as follows: From the analysis of test results, it can be obtained that the elastic modulus E value of Cr5 alloy steel is 55 GPa under the condition of 800 °C and it is only 3.5 GPa at 1150 °C. These are quite different from the elastic modulus of 206 GPa of Cr5 alloy steel at room temperature. Additionally, the fitting function of temperature and elastic modulus has a high linear relationship, and the elastic model at other temperatures can be calculated according to Equation (17). This shows that after the increase in the temperature of Cr5 alloy steel, the atomic spacing inside the material increases, the thermal motion of the molecules intensifies, and then the ability to resist external elastic deformation is weakened.  Figure 6 is a comparison diagram of the damage resistance factor S value of the Lemaitre damage model calculated by substituting the elastic modulus (206 GPa, red) of the Cr5 alloy steel at room temperature and the elastic modulus E value measured by the test (black). This figure shows that the damage resistance factor S value directly calculated by the constant elastic modulus of 206 GPa shows little difference at each temperature. Additionally, at the same temperature, the difference of the damage resistance factor S value under each strain rate is not obvious. However, the S value of the damage resistance factor calculated using the E value measured by the actual elastic modulus measurement test shows obvious differences at each temperature. Additionally, at the same temperature, the value of the damage resistance factor S at each strain rate is also quite different. This shows that the value of the elastic modulus E has great influence over obtaining an accurate material damage model. A more accurate Lemaitre damage model for the Cr5 alloy steel at different temperatures and strain rates can be obtained incorporating the E value measured by the actual elastic modulus measurement test.  From the analysis of test results, it can be obtained that the elastic modulus E value of Cr5 alloy steel is 55 GPa under the condition of 800 • C and it is only 3.5 GPa at 1150 • C. These are quite different from the elastic modulus of 206 GPa of Cr5 alloy steel at room temperature. Additionally, the fitting function of temperature and elastic modulus has a high linear relationship, and the elastic model at other temperatures can be calculated according to Equation (17). This shows that after the increase in the temperature of Cr5 alloy steel, the atomic spacing inside the material increases, the thermal motion of the molecules intensifies, and then the ability to resist external elastic deformation is weakened. Figure 6 is a comparison diagram of the damage resistance factor S value of the Lemaitre damage model calculated by substituting the elastic modulus (206 GPa, red) of the Cr5 alloy steel at room temperature and the elastic modulus E value measured by the test (black). This figure shows that the damage resistance factor S value directly calculated by the constant elastic modulus of 206 GPa shows little difference at each temperature. Additionally, at the same temperature, the difference of the damage resistance factor S value under each strain rate is not obvious. However, the S value of the damage resistance factor calculated using the E value measured by the actual elastic modulus measurement test shows obvious differences at each temperature. Additionally, at the same temperature, the value of the damage resistance factor S at each strain rate is also quite different. This shows that the value of the elastic modulus E has great influence over obtaining an accurate material damage model. A more accurate Lemaitre damage model for the Cr5 alloy steel at different temperatures and strain rates can be obtained incorporating the E value measured by the actual elastic modulus measurement test. calculated using the E value measured by the actual elastic modulus measurement test shows obvious differences at each temperature. Additionally, at the same temperature, the value of the damage resistance factor S at each strain rate is also quite different. This shows that the value of the elastic modulus E has great influence over obtaining an accurate material damage model. A more accurate Lemaitre damage model for the Cr5 alloy steel at different temperatures and strain rates can be obtained incorporating the E value measured by the actual elastic modulus measurement test.

Determination of Damage Model Parameters
The stress and strain data of the Cr5 alloy steel at different temperatures and strain rates can be obtained from the tensile tests. The Lemaitre damage parameters and Zener-Hollomon coefficient values of the Cr5 alloy steel at different temperatures and strain rates can be solved by Formulas (5)-(9). The calculation results are shown in Table 4.

Determination of Damage Model Parameters
The stress and strain data of the Cr5 alloy steel at different temperatures and strain rates can be obtained from the tensile tests. The Lemaitre damage parameters and Zener-Hollomon coefficient values of the Cr5 alloy steel at different temperatures and strain rates can be solved by Formulas (5)-(9). The calculation results are shown in Table 4.

Lemaitre High-Temperature Damage Model of Cr5 Alloy
Coupling the Lemaitre damage model with the Zener-Hollomon coefficient, a Lemaitre high-temperature damage model of Cr5 alloy steel was founded. The damage parameters in Table 4 were polynomial fitted with ln(Z) based on the least square method. The fitting parameters of the obtained Lemaitre damage model coupled with the Zener-Hollomon coefficient are shown in Table 5. The fitted curve is shown in Figure 7.  It can be seen from Figure 7c that the fitting of anti-damage factor S is greatly dispersed, which may be caused by the difference in the influence degree of temperature and strain rate on S. Finally, the high-temperature damage model of the Cr5 alloy steel Lemaitre can be obtained:  It can be seen from Figure 7c that the fitting of anti-damage factor S is greatly dispersed, which may be caused by the difference in the influence degree of temperature and strain rate on S. Finally, the high-temperature damage model of the Cr5 alloy steel Lemaitre can be obtained: The correlation coefficient R, R-squared R 2 and mean square error MSE are often used to measure the linear relationship between two variables. The larger the R and R 2 value, the smaller the MSE, and the better the fitting effect. They are calculated as follows: In the formula, X i is the value of the ith simulation; Y i is the ith test value; X is the average value of the simulation; Y is the average value of the test value. Substitute the experimental and simulation values of damage strain threshold ε th , damage resistance factor S, critical damage value D 1c , and ultimate stress value σ u into Equations (19)-(21) for calculation, and obtain the data in Table 6. It can be seen from Table 6 that the correlation coefficient R values obtained after fitting are all above 0.90. This shows that the established Cr5 alloy steel Lemaitre high-temperature damage model has high reliability. It can be seen from Figure 7 that the high-temperature Lemaitre damage model of the Cr5 alloy steel has a strong correlation with the actual test values [41]. This shows that the high-temperature Lemaitre damage model of Cr5 alloy steel can describe the damage characteristics of Cr5 alloy steel at different temperatures and strain rates in industrial production. By coupling and correcting the Lemaitre damage model with the Zener-Hollomon coefficient and the elastic modulus value E of the Cr5 alloy steel at different temperatures, it can provide an effective means of controlling the crack defects of Cr5 alloy steel roll products in the actual industry.

Finite Element Modelling of Tensile Process of Cr5 Alloy Steel
The reliability of the high-temperature damage model of the Cr5 alloy steel Lemaitre can effectively be verified by the finite element simulation technology. Using FORTRAN language, the Cr5 alloy high-temperature damage model was programmed, and it was embedded into Forge ® software through dynamic link library, and performed a numerical simulation of tensile test process of Cr5 alloy steel [42]. By comparing the simulation results with the test results, the accuracy of the high-temperature damage model of the Cr5 alloy steel Lemaitre can be verified.
Since only the temperature of the centre of the sample can reach the specified temperature during the actual tensile test, the temperature on both sides gradually decreases and the magnitude of the temperature decrease increases with the distance from the centre of the sample. During the tensile test, a thermocouple wire was used to measure the temperature of the sample and the sample was kept at this temperature. The temperature distribution of the tensile sample with the length of the sample is shown in Figure 8. The length scale takes the centre of the sample as the origin. The temperature of the simulated tensile specimen was set according to the actual temperature distribution. The finite element geometric model and temperature distribution are shown in Figure 9. The mesh size is set to 0.5, the temperature and strain rate refer to the experimental conditions, the Coulomb friction coefficient is 0.3, the heat exchange coefficient is 2000 W/m 2 ·K, and the emissivity is 0.88.

Verification of Lemaitre High-Temperature Damage Model of Cr5 Alloy Steel
By way of the observation and analysis of the digital simulation results of the tensile test, the displacement-load curves of the Cr5 alloy steel at various temperatures and strain rates are obtained, and the results obtained by digital simulation are compared with the actual tensile test results of the Cr5 alloy steel, as shown in Figure 10.

Verification of Lemaitre High-Temperature Damage Model of Cr5 Alloy Steel
By way of the observation and analysis of the digital simulation results of the tensile test, the displacement-load curves of the Cr5 alloy steel at various temperatures and strain rates are obtained, and the results obtained by digital simulation are compared with the actual tensile test results of the Cr5 alloy steel, as shown in Figure 10.

Verification of Lemaitre High-Temperature Damage Model of Cr5 Alloy Steel
By way of the observation and analysis of the digital simulation results of the tensile test, the displacement-load curves of the Cr5 alloy steel at various temperatures and strain rates are obtained, and the results obtained by digital simulation are compared with the actual tensile test results of the Cr5 alloy steel, as shown in Figure 10.

Verification of Lemaitre High-Temperature Damage Model of Cr5 Alloy Steel
By way of the observation and analysis of the digital simulation results of the te test, the displacement-load curves of the Cr5 alloy steel at various temperatures and s rates are obtained, and the results obtained by digital simulation are compared with actual tensile test results of the Cr5 alloy steel, as shown in Figure 10.   As can be seen from Figure 10, the displacement-load curve obtained by the ten simulation test of Cr5 alloy steel is compared with the actual tensile test, and at 800 850 °C, 900 °C and 950 °C, they basically overlap; at 1000 °C, 1050 °C, 1100 °C and °C, the deviation is small. By comparison, it can be found that the displacementcurve obtained by the tensile simulation test is basically consistent with the displacem load curve obtained by the actual tensile test. The correlation coefficient R 2 was use compare the deviation between the experimental and simulated displacementcurves. The simulation test data and the displacement-load data of the actual tensile are brought into Formula (20). It was calculated that R 2 = 0.9874. Figure 11 is a schematic diagram of the comparison of the gauge length and frac of the simulated tensile specimen of Cr5 alloy steel at 950 °C with a strain rate of 5 s −1 the actual tensile specimen.  As can be seen from Figure 10, the displacement-load curve obtained by the tensile simulation test of Cr5 alloy steel is compared with the actual tensile test, and at 800 • C, 850 • C, 900 • C and 950 • C, they basically overlap; at 1000 • C, 1050 • C, 1100 • C and 1150 • C, the deviation is small. By comparison, it can be found that the displacement-load curve obtained by the tensile simulation test is basically consistent with the displacement-load curve obtained by the actual tensile test. The correlation coefficient R 2 was used to compare the deviation between the experimental and simulated displacement-load curves. The simulation test data and the displacement-load data of the actual tensile test are brought into Formula (20). It was calculated that R 2 = 0.9874. Figure 11 is a schematic diagram of the comparison of the gauge length and fracture of the simulated tensile specimen of Cr5 alloy steel at 950 • C with a strain rate of 5 s −1 and the actual tensile specimen.
As can be seen from Figure 11, based on the high-temperature damage model of Cr5 alloy steel Lemaitre, under the conditions of 950 • C and strain rate of 5 s −1 , the fracture gauge length difference between the numerical simulation tensile specimen and the actual tensile specimen is 1.28%, the difference of the area reduction in the fracture surface is 1.1% and the difference of the area reduction at the centre between the fracture surface and the nearest large end face of the sample is 2.1%. The accuracy and reliability of the high-temperature damage model of Cr5 alloy steel Lemaitre are verified, and it shows that the high-temperature damage model is suitable for predicting the plastic damage of Cr5 alloy steel at different temperatures and strain rates. It can effectively prevent the occurrence of crack defects in the production of Cr5 rolls. load curve obtained by the actual tensile test. The correlation coefficient R 2 was used to compare the deviation between the experimental and simulated displacement-load curves. The simulation test data and the displacement-load data of the actual tensile test are brought into Formula (20). It was calculated that R 2 = 0.9874. Figure 11 is a schematic diagram of the comparison of the gauge length and fracture of the simulated tensile specimen of Cr5 alloy steel at 950 °C with a strain rate of 5 s −1 and the actual tensile specimen.  As can be seen from Figure 11, based on the high-temperature damage model of Cr5 alloy steel Lemaitre, under the conditions of 950 °C and strain rate of 5 s −1 , the fracture gauge length difference between the numerical simulation tensile specimen and the actual tensile specimen is 1.28%, the difference of the area reduction in the fracture surface is 1.1% and the difference of the area reduction at the centre between the fracture surface Figure 11. Schematic diagram of the comparison of gauge distance and tensile fracture between simulated tensile specimen and actual tensile specimen at 950 • C and 5 s −1 (unit: mm).

Damage Map of Lemaitre High-Temperature Damage Model
In order to have a more comprehensive understanding of the damage characteristics of Cr5 alloy steel and to more intuitively discern the relationship between the various damage parameter values of the Lemaitre damage model, the calculated damage parameter values of the Lemaitre high-temperature damage model at different temperatures and different strain rates of Cr5 alloy steel are plotted with their corresponding temperatures and strain rates. The damage map of the Lemaitre high-temperature damage model of the Cr5 alloy steel can be obtained, as shown in Figure 12.
As can be seen from Figure 12a, Cr5 alloy steel ultimate stress value σ u is inversely proportional to temperature and proportional to the strain rate. When the temperature is 800 • C and the strain rate is 0.01 s −1 , the ultimate stress value σ u is 150.95 MPa; when the strain rate is increased to 5 s −1 , σ u is 347.46 MPa; when the strain rate is 0.01 s −1 and the temperature is increased to 1150 • C, the ultimate stress value σ u is 43.28 MPa. It can be seen from the whole relationship diagram that the maximum value of the ultimate stress value σ u appears at 800 • C and 5 s −1 , and the ultimate stress value σ u is 347.46 Mpa; the minimum value of the ultimate stress value σ u appears at 1150 • C, 0.01 s −1 , and the ultimate stress value σ u at this time is 43.28 MPa.
It can be seen from Figure 12b that the damage resistance factor S of the Cr5 alloy steel is proportional to the temperature and inversely proportional to the strain rate. For example, when the strain rate is 0.01 s −1 and the temperature is 800 • C, the damage resistance factor S is 0.059 MPa, and the damage resistance factor S is 0.265 MPa when the strain rate is increased to 5 s −1 ; when the temperature is increased to 1150 • C, the damage resistance factor S is 0.134 MPa. It can be seen from the whole relationship diagram that the maximum value of the anti-damage factor S appears at 1150 • C, 5 s −1 , and the anti-damage factor S at this time is 0.870 Mpa; the minimum value of the anti-damage factor S appears at 800 • C/0.01 s −1 , and the damage resistance factor S is 0.059MPa.
It can be seen from Figure 12c that the critical damage value D 1c of the Cr5 alloy steel is basically inversely proportional to the temperature and proportional to the strain rate. For example, when the strain rate is 0.01 s −1 and the temperature is 800 • C, the critical damage D 1c is 0.919, and the critical damage D 1c is 0.995 when the strain rate is increased to 5 s −1 ; when the temperature is increased to 1150 • C, the critical damage D 1c at 1150 • C is 0.701. From the whole relationship diagram, it can be seen that the maximum value of critical damage D 1c occurs at 800 • C/5 s −1 and the critical damage D 1c at this time is 0.995; the minimum value of critical damage D 1c occurs at 1150 • C/0.01 s −1 , and the critical damage D 1c at 0.701. rence of crack defects in the production of Cr5 rolls.

Damage Map of Lemaitre High-Temperature Damage Model
In order to have a more comprehensive understanding of the damage characteristi of Cr5 alloy steel and to more intuitively discern the relationship between the variou damage parameter values of the Lemaitre damage model, the calculated damage param eter values of the Lemaitre high-temperature damage model at different temperatures an different strain rates of Cr5 alloy steel are plotted with their corresponding temperatur and strain rates. The damage map of the Lemaitre high-temperature damage model of th Cr5 alloy steel can be obtained, as shown in Figure 12. As can be seen from Figure 12a, Cr5 alloy steel ultimate stress value σu is inverse proportional to temperature and proportional to the strain rate. When the temperature 800 °C and the strain rate is 0.01 s −1 , the ultimate stress value σu is 150.95 MPa; when th strain rate is increased to 5 s −1 , σu is 347.46 MPa; when the strain rate is 0.01 s −1 and th temperature is increased to 1150 °C, the ultimate stress value σu is 43.28 MPa. It can b seen from the whole relationship diagram that the maximum value of the ultimate stre value σu appears at 800 °C and 5 s −1 , and the ultimate stress value σu is 347.46 Mpa; th minimum value of the ultimate stress value σu appears at 1150 °C, 0.01 s −1 , and the ultima stress value σu at this time is 43.28 MPa. It can be seen from Figure 12d that the damage strain threshold ε th of the Cr5 alloy steel is proportional to the strain rate and has little change under the influence of temperature. For example, when the strain rate is 0.01 s −1 and the temperature is 800 • C, the damage strain threshold ε th is 0.17, and the damage strain threshold ε th is 0.29 when the strain rate is increased to 5 s −1 ; when the strain rate is 0.01 s −1 and the temperature is 1150 • C, the damage strain threshold ε th is 0.13, and the damage strain threshold ε th is 0.28 when the strain rate is increased to 5 s −1 . When the strain rate is 1 s −1 , the damage strain threshold values ε th at 800-1150 • C are 0.22, 0.29, 0.27, 0.26, 0.26, 0.23, 0.23 and 0.22, and the change is not too obvious. It can be seen from the whole relationship diagram that the maximum damage strain threshold ε th appears at 850 • C/5 s −1 , the damage strain threshold ε th at this time is 0.33; and the minimum damage strain threshold ε th appears at 1150 • C/0.01 s −1 and the damage strain threshold ε th at this time is 0.13.
From the damage diagram of the Lemaitre high-temperature damage model of the Cr5 alloy steel, it can be seen that when the temperature is constant, as the strain rate gradually increases, the critical damage value D 1c is basically positively correlated with the damage strain threshold value ε th . However, when the strain rate is constant, the correlation between the two is lower with the gradual increase in temperature; when the strain rate is constant, the ultimate stress σ u and the damage resistance factor S are basically negatively correlated with the gradual increase in temperature; when the temperature is constant, there is a positive correlation between the ultimate stress σ u and the damage resistance factor S as the strain rate gradually increases.

Conclusions
In order to accurately predict the crack defects of the Cr5 alloy steel during hightemperature forming. In the present study, the thermal tensile test and elastic modulus measurement experiment of Cr5 alloy under the conditions of 800-1150 • C and 0.01-5 s −1 were carried out, and the high-temperature damage model and damage diagram of the Cr5 alloy steel Lemaitre were established.
The elastic modulus E value of the Cr5 alloy steel decreases with the increase in temperature in the temperature range of 800-1150 • C. This is because when the temperature rises, the atomic spacing inside the material increases, the thermal motion of molecules increases, and then the ability to resist external elastic deformation is weakened. The test results show that the decrease in elastic modulus at high temperature increases the damage resistance factor S. In order to more accurately describe the effect of temperature on the elastic modulus, the functional relationship between the elastic modulus of the Cr5 alloy steel and temperature is determined.
Considering the influence of the temperature and strain rate on the damage of Cr5 alloy steel, Zener-Hollomon coefficient was introduced, the parameters of Lemaitre hightemperature damage model were measured, and the Lemaitre damage model of the Cr5 alloy steel at a high-temperature stage was founded.
The established high-temperature damage model was embedded in a Forge ® finite element software through the program's secondary development method to carry out the numerical simulation calculation of the strength of the experimental samples. Comparing the difference between the simulation and the test displacement-load curve, the obtained correlation coefficient (R 2 ) was 0.987. The fracture gauge length of the tensile specimen differed by 1.28% and the fracture morphology was basically consistent. This shows that the high-temperature damage model of the Cr5 alloy steel Lemaitre established in this paper has high accuracy in predicting the damage behaviour of this material. Finally, the damage map of Cr5 alloy steel at the high-temperature stage was constructed. The relationship between the damage parameters of the high-temperature damage model of Cr5 alloy steel and the deformation temperature and strain rate were analysed. This provides an effective method for controlling crack defects in the Cr5 alloy steel during high-temperature plastic deformation. In this paper, a high temperature damage model is established for Cr5 alloy steel. In the future, the production process of the Cr5 large roller will be simulated and verified in real industry.
Author Contributions: Conceptualization, X.C.; original manuscript writing, L.G.; method theory, X.C.; software, B.Z. and R.B.; data processing and data analysis, L.G.; writing-review, R.B.; revision of the manuscript and final proofread manuscript, B.Z. and L.G. All authors have read and agreed to the published version of the manuscript.