Next Article in Journal
How Retting Could Affect the Mechanical Behavior of Flax/Epoxy Biocomposite Materials?
Next Article in Special Issue
Mechanical Alloying Behavior and Thermal Stability of CoCrCuFeMnNix High-Entropy Alloy Powders Prepared via MA
Previous Article in Journal
Fatigue Experiment Study on Internal Force Redistribution in the Negative Moment Zone of Steel–Concrete Continuous Composite Box Beams
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Quadratic Constitutive Model Based on Partial Derivative and Taylor Series of Ti6242s Alloy and Predictability Analysis

1
Chongqing Key Laboratory of Advanced Mold Intelligent Manufacturing, College of Materials Science and Engineering, Chongqing University, Chongqing 400044, China
2
Chongqing Jiepin Technology Co., Ltd., Chongqing 400000, China
*
Author to whom correspondence should be addressed.
Materials 2023, 16(7), 2928; https://doi.org/10.3390/ma16072928
Submission received: 10 March 2023 / Revised: 29 March 2023 / Accepted: 31 March 2023 / Published: 6 April 2023
(This article belongs to the Special Issue Research on Forming and Serving Performance of Advanced Alloys)

Abstract

:
To solve the problem of insufficient predictability in the classical models for the Ti6242s alloy, a new constitutive model was proposed, based on the partial derivatives from experimental data and the Taylor series. Firstly, hot compression experiments on the Ti6242s alloy at different temperatures and different strain rates were carried out, and the Arrhenius model and Hensel–Spittel model were constructed. Secondly, the partial derivatives of logarithmic stress with respect to temperature and logarithmic strain rate at low, medium and high strain levels were analyzed. Thirdly, two new constitutive models with first- and second-order approximation were proposed to meet the requirements of high precision. In this new model, by analyzing the high-order differential data of experimental data and combining the Taylor series theory, the minimum number of terms that can accurately approximate the experimental rheological data was found, thereby achieving an accurate prediction of flow stress with minimal material parameters. In the new model, by analyzing the high-order differential of the experimental data and combining the theory of the Taylor series, the minimum number of terms that can accurately approximate the experimental rheological data was found, thereby achieving an accurate prediction of flow stress with minimal material parameters. Finally, the prediction accuracies for the classical model and the new model were compared, and the predictabilities for the classical models and the new model were proved by mathematical means. The results show that the prediction accuracies of the Arrhenius model and the Hensel–Spittel model are low in the single-phase region and high in the two-phase region. In addition, second-order approximation is required between the logarithmic stress and logarithmic strain rate, and first-order approximation is required between logarithmic stress and temperature to establish a high-precision model. The order of prediction accuracy of the four models from high to low is the quadratic model, Arrhenius model, linear model and HS model. The prediction accuracy of the quadratic model in all temperatures and strain rates had no significant difference, and was higher than the other models. The quadratic model can greatly improve prediction accuracy without significantly increasing the material parameters.

1. Introduction

The Ti6242s alloy, a typical high-temperature titanium alloy, is widely used in high-pressure compressor disks and blades in aeroengines; these parts are usually manufactured by a hot-working process [1,2], such as hot forging, hot extrusion, hot roll forging, etc. The most important purpose of hot-working processes is that metals are deformed into the desired shapes under certain thermomechanical processing conditions to obtain the appropriate microstructure and mechanical performance [3]. Numerical simulation technology is an important technical means for achieving the accurate control of shape and performance, and one of the cores of accurate numerical simulation is to build a high-precision constitutive model [4].
In past decades, many scholars have studied the constitutive model of titanium alloy. Three main categories of constitutive models are utilized to predict the flow behavior of metallic alloys, which are based on physical, phenomenological and artificial neural network constitutive models [5,6,7,8,9]. The most widely used constitutive models in numerical simulation software are the Arrhenius (AH) model [10,11], the Johnson-Cook (JC) model [12], the Hensel–Spittel (HS) model [13,14] and so on. Since these constitutive models were proposed, many subsequent studies have attempted to improve these classical models. For example, Kotkunde et al. [15] used four constitutive models (a modified JC model, modified AH model, modified Zerilli–Armstrong and the Rusinek–Klepaczko model) to analyze the accuracy of the rheological behavior of the Ti-6Al-4V alloy. The results showed that the modified AH model had high accuracy. Zhang et al. [16] proposed the strain-compensated AH and modified Norton–Hoff constitutive models for the Ti-6Al-4V alloy. The results showed that the Norton–Hoff model had higher accuracy. Xiao et al. [17] proposed a new constitutive model of the TiNiNb alloy, in which the effects of temperature and strain rate on stress were considered. The regression results of the experimental data showed that the predicted accuracy of the proposed constitutive model was very high. Luo et al. [18] analyzed the Kocks–Mecking constitutive model of the Ti-6Al-4V alloy, and the average relative error was about 6.13%. Cai et al. [19] conducted a constitutive analysis of the Ti-6Al-4V alloy by using the stress–strain data obtained from an isothermal hot compression experiment. Considering the independent effects of strain, strain rate and temperature, a modified Arrhenius constitutive model was proposed. The results showed that the modified parallel constitutive model based on multiple regression could better predict the flow stress of the Ti-6Al-4V alloy, and had good correlation and generalization ability. Ming et al. [20] compared the prediction errors of the modified AH model and the modified JC model for the Ti-6Al-4V alloy. The results showed that the modified AH model with strain compensation was more accurate than the modified JC model. Ga et al. [21] compared the prediction accuracy of five constitutive models for the flow data of the Ti-6Al-4V alloy. The results showed that the modified HS model had the highest accuracy and the JC model had the worst accuracy. Recently, an artificial neural network model was also introduced into the prediction of rheological data. For example, Ahmed et al. [22] compared the prediction accuracy between the AH model and an artificial neural network (ANN) model. The results showed that the predictability of the artificial neural network model was higher than that of the AH model for the Ti-2.5Al-1.8Mn alloy. Sun et al. [23] proposed a BP neural network constitutive model of the Ti40 alloy, and the results showed that this model had a high prediction accuracy for flow stress. Reddy et al. [24] established a BP neural network constitutive model to predict the flow stress of the Ti-6Al-4V alloy. This model was successfully trained by using experimental flow data in both the double phase α + β region and the single phase (β) region. This model seemed to have a higher prediction accuracy than the classical models. However, there are many obvious disadvantages for the neural network model, such as many offset and weight parameters, no gradient information (simple expression), over fitting, etc. In addition, compared with the traditional models, the calculation efficiency of the artificial neural network model is extremely low. For finite element simulation with a strict speed requirement, the artificial neural network model will obviously not be adopted by the mainstream commercial simulation software. Because the AH, JC and HS models have the characteristics of simple structure and high accuracy, they are currently the mainstream models.
However, it is difficult for classical models to accurately predict the flow stress in both the α + β phase and β phase regions due to the phase transition of the Ti6242s alloy within the hot forming temperature range. In addition, the second-order approximation is required between the logarithmic stress and logarithmic strain rate, and first-order approximation is required between logarithmic stress and temperature to establish a high-precision model. However, the classical constitutive models have only first-order accuracy; therefore, these modifications and enhancements conducted by the above researchers have limited improvement in the prediction accuracy, only within the range of the first-order accuracy.
To solve the problem of insufficient predictability in the classical constitutive models, a novel constitutive model of the Ti6242 alloy was established based on the partial derivative and Taylor series. The novel constitutive model can predict the high-temperature flow stress of the Ti6242s alloy in both the α + β phase and single β phase regions more precisely. Firstly, the isothermal hot compression tests were completed in the temperature range of 900~1100 °C and a strain rate range of 0.001~10 s–1 on a Gleeble-3500 simulator. Secondly, two new constitutive models based on the Taylor series were proposed without significantly increasing the material parameters. In this novel building method, the minimum number of terms that can accurately approximate the experimental rheological data was found by analyzing the high-order differentials of experimental data and combining the theory of the Taylor series, thereby achieving an accurate prediction of flow stress with minimal material parameters. Finally, the prediction accuracy of the proposed constitutive equation was evaluated by comparing with classical constitutive models. The results showed that the prediction accuracy of the classical models was low in the single-phase region and high in the two-phase region. The order of prediction accuracy of the four models from high to low is quadratic model, Arrhenius model, linear model and HS model. The prediction accuracy of the quadratic model in all temperatures and strain rates is higher than the other models and it can greatly improve the prediction accuracy without significantly increasing the material parameters.

2. Materials and Experiments

2.1. Material

The experimental material was the Ti6242s alloy and its nominal chemical composition is shown in Table 1. The main alloy elements of the experimental material were Al, Sn, Zr and Mo. It is a kind of near-α titanium alloy with a service temperature between 450 °C and 500 °C. This material is mainly used for medium and high pressure compressor discs and blades in aircraft engines [25]. It has excellent high temperature durability, weldability and processability.
Figure 1 shows the original forged microstructure of the material. It can be seen that the primary microstructure of the material displays the typical bimodal morphology which is mainly composed of massive equiaxed primary α phase (more than 50%), lamellar secondary α phase and a small amount of transformed phase β.
The hot processing or heat treatment of the titanium alloy is carried out in a two-phase zone or a single-phase zone, which has a certain impact on the structure and properties. Therefore, the phase transformation point of the titanium alloy is an important parameter in its processing specifications. For the Ti6242s alloy, the phase transform temperature from α to β is 1015 °C.

2.2. Experimental Details and Results

Generally, the range of the forming temperature of titanium alloy is located in a double phase region [26,27,28]. Therefore, in the isothermal hot compression, the temperatures located in both the double phase region (900, 950, 1000 °C) and single phase region (1050 and 1100 °C) were selected to build a complete database. Twenty-five experimental specimens were machined from the original metal. The size of the specimens was 8 mm in diameter and 12 mm in height. As shown in Table 2, 25 samples were heated to 900, 950, 1000, 1050 and 1100 °C at the rate of 5 °C/s and then subjected to isothermal compression on the Geeble-3500 thermal simulation testing machine. The compression deformation rates were 0.001, 0.01, 0.1, 1 and 10 s−1, respectively, and the compression amounts were 60% (the true strain was 0.916). To retain the microstructure, the specimens were immediately placed in water. To eliminate the friction, tantalum foil with a thickness of 0.1 mm was placed between die and specimen before hot compression.
Figure 2 shows the stress–strain curves at various strain rates and temperatures. The following rules can be obtained from Figure 2: (1) the flow stress decreases with the increase in temperature under the same strain rate. The higher the temperature is, the stronger the atomic activity and the stronger the dislocation motion ability are, and the speed of mutual cancellation between dislocations is accelerated. Meanwhile, with the increase in temperature, α→β phase transformation is triggered. The slip system of β phase is more than α phase. Then, this phase transformation causes deformation to occur more easily. In addition, the increase in temperature accelerates the nucleation and growth rate of recrystallization, and the softening effect is enhanced. (2) In the temperature range of both α + β phase region (900~1000 °C), a significant stress drop can be observed at all strain rates, and the flow softening phenomenon is more obvious when the temperature lowers. (3) In single β phase zone (1050~1100 °C), the curve will soften only when the strain rate is low, and the curve will decline slightly or remain basically unchanged after reaching the peak value at the high strain rate. The reason is that the temperature rises above the recrystallization termination temperature of the material; recrystallization stops in the grains, and the softening effect of dynamic recovery is not obvious; under the counteraction of deformation hardening, it shows relatively stable flow stress.

2.3. Strain Dispersion

In order to regress the constitutive model easily, the flow stress curve was divided into 10 equal parts according to the strain from 0.04 to 0.9. It is worth noting that the accuracy increases with more divided parts in a certain range. In order to eliminate the influence of the number of equal parts on the accuracy of the model, all subsequent models used the same number of equal parts, that is, the data in Table 3. The data in Table 3 can be obtained by linear interpolation using each rheological curve in Figure 2.

3. Classical Constitutive Model

The accurate construction of the prediction model of the flow curve is the premise of accurate numerical simulation. Before proposing a new phenomenological constitutive model, two classical constitutive models, namely the Arrhenius model and the Hensel–Spittel (HS) model, were analyzed.

3.1. Arrhenius Model

The Arrhenius constitutive model has been widely used due to its high accuracy and wide generalizability. The model was first proposed by Sellars and McTegart [29], and its specific form is presented in Equation (1).
{ ε ˙ = A σ β α exp ( Q R T ) ,   ( α σ 0.8 ) ε ˙ = A exp ( β σ ) exp ( Q R T ) , ( α σ 1.2 ) ε ˙ = A [ sinh ( α σ ) ] n exp ( Q R T ) ,   ( for   all )
where A, α, n and β are material parameters, Q is the activation energy of thermal deformation (J·mol−1), R is the universal gas constant (8.314 J·K−1·mol−1), T is the absolute temperature (K), ε ˙ is the strain rate (s−1), σ is flow stress (MPa). In Equation (1), the effect of strain on flow stress was not considered. In recent years, most scholars consider how to take the effect of strain into account [20]. Equation (2) can be obtained by logarithmic transform in Equation (1).
{ ln ε ˙ = ln A + β α ln σ Q R T , ( α σ 0.8 ) ln ε ˙ = ln A + β σ Q R T , ( α σ 1.2 ) ln ε ˙ = ln A + n ln sinh ( α σ ) Q R T , ( for   all )
The material parameters of ln A , α , n , Q at different strain levels can be obtained by multivariate nonlinear regression of equation ln ε ˙ = ln A + n ln sinh ( α σ ) Q / ( R T ) using flow stress data (shown in Figure 2). For example, the strain has been divided into 10 parts between 0.04 and 0.9. Then, the stress matrix (shown in Table 3), corresponding to different strain levels, was obtained by interpolation. The material parameters of ln A , α , n and Q corresponding to different strain levels can be obtained by multivariate nonlinear regression using the stress matrix and equation ln ε ˙ = ln A + n ln sinh ( α σ ) Q / ( R T ) . Then, the polynomial fittings of ln A - ε , α - ε , n - ε and Q - ε curves (shown in Figure 3) were performed, respectively, to obtain the constitutive equation with strain compensation.
As shown in Figure 3, the ten points in each subgraph correspond to strains of 0.04, 0.14, 0.23, …, 0.90. The fifth-degree polynomial is used to fit these data, and expression of each parameters is obtained. The degree of polynomial can be determined according to the regression accuracy, but too much will lead to overfitting.
{ ln A = 43.91 + 192.62 ε 969.19 ε 2 + 2137.32 ε 3 2236.13 ε 4 + 887.20 ε 5 α = 0.01 0.05 ε + 0.22 ε 2 0.45 ε 3 + 0.47 ε 4 0.19 ε 5 n = 4.61 3.33 ε 18.75 ε 2 + 46.38 ε 3 57.78 ε 4 + 26.69 ε 5 Q = 491161 + 1703450 ε 9111853 ε 2 + 20591837 ε 3 21622908 ε 4 + 8542954 ε 5
The Arrhenius constitutive equation of the Ti6242s alloy can be obtained by introducing Formula (3) into (1). In order to observe the prediction accuracy of the model more intuitively, the experimental data and the prediction data of the Arrhenius constitutive equation were plotted in the same figure.
As shown in Figure 4, the prediction accuracy of the Arrhenius constitutive equation is significantly different in different temperatures and strain rates. The prediction accuracy of the Arrhenius constitutive equation is low in the single-phase region (900, 950 and 100 °C) and high in the two-phase region (1050 and 1100 °C), where the prediction accuracy does not conform to this rule when the strain rate is 0.1 s−1 and the temperature is 1050 °C. An ideal constitutive model should have the same prediction accuracy at different temperatures and strain rates. Therefore, the standard Arrhenius constitutive equation needs to be improved when the prediction accuracy of the rheological properties of α+β titanium alloys is required to be high.

3.2. Hensel–Spittel Model

The Hensel–Spittel (HS) constitutive model, which has been used in famous commercial software (Forge NxT), is frequently used in hot forming simulation. Its form is simple and its parameters can be obtained easily [30,31]. The general form of the HS model is exhibited in Equation (4).
σ = Aexp ( m 1 T ) ε m 2 ε ˙ m 3 exp ( m 4 ε ) ( 1 + ε ) m 5 T exp ( m 6 ε ) ε ˙ m 7 T T m 8
where ε , σ , ε ˙ , T are strain, flow stress, strain rate and temperature, A and m 1 m 8 are material constants. Equation (5) can be obtained by taking the natural logarithm of Equation (4).
ln σ = ln A + m 1 T + m 2 ln ε + m 3 ln ε ˙ + m 4 / ε + m 5 T ln ( 1 + ε ) + m 6 ε + m 7 T ln ε ˙ + m 8 ln T
There is a linear relationship between ln σ , T , ln ε , ln ε ˙ , 1 / ε , T ln ( 1 + ε ) , ε , T ln ε ˙ , and ln T . Then, the solving problem of the material constants in Equation (5) is a typical multiple linear regression problem. Similarly, the strain has been divided into 10 parts between 0.04 and 0.9. Then, the flow stress at each strain level can be obtained by interpolating the original compressed data (shown in Figure 1). The stress, strain, temperature and strain rate corresponding to all of the strain levels should approximately satisfy Equation (5) (not absolutely satisfied due to regression error). All strains are brought into Equation (5) and a random error factor is introduced to obtain linear simultaneous Equation (6).
[ ln σ 1 ln σ 2 ln σ 3 ln σ n ] = [ 1 1 1 1   T 1 T 2 T 3 T n   ln ε 1 ln ε 2 ln ε 3 ln ε n   ln ε ˙ 1 ln ε ˙ 2 ln ε ˙ 3 ln ε ˙ n   1 / ε 1 1 / ε 2 1 / ε 3 1 / ε n   T 1 ln ( 1 + ε 1 ) T 2 ln ( 1 + ε 2 ) T 3 ln ( 1 + ε 3 ) T n ln ( 1 + ε n ) ε 1 ε 2 ε 3 ε n   T 1 ln ε ˙ 1 T 2 ln ε ˙ 2 T 3 ln ε ˙ 3 T n ln ε ˙ n   ln T 1 ln T 2 ln T 3 ln T n ] γ + μ
where n = 5 × 5 × 10 (5 temperatures, 5 strain rates and 10 strain levels), γ = [ ln A   m 1   m 2   m 3   m 4   m 5   m 6   m 7   m 8 ] , μ is the error vector with size of n × 1, which follows the normal distribution with the mean value of zero. Equation (6) is an overdetermined equation in which the number of variables is fewer than the number of equations. The least square method is a classical method used to solve this overdetermined equation. The γ can be solved by using function of “ γ = regress(Y, X)” in the software of MATLAB. The multiple linear regression results of the material parameters are listed in Table 4, and the comparisons of the HS model prediction flow stress and experimental flow stress are shown in Figure 5.
As shown in Figure 5, the prediction accuracy of the HS constitutive equation is low in the single-phase region (900, 950 and 100 °C) and high in the two-phase region (1050 and 1100 °C). This rule is the same as for the Arrhenius model. The prediction accuracy of the HS model is lower than that of the Arrhenius model, but the HS model has only 9 material parameters and the Arrhenius model has 24 material parameters. Generally, the more parameters of the model, the higher the prediction accuracy. So far, it is difficult to improve the prediction accuracy of either model. Based on this, a new constitutive model is urgently needed.

4. New Constitutive Model

4.1. Mathematical Principles

An idea can be drawn from the Arrhenius model; that is, when the strain is fixed, the core of building the constitutive model of materials is to build the functional relationship between stress, strain rate and temperature. The Arrhenius model is an implicit equation, which is complicated in practical application. In addition, the prediction accuracy of this model is also low for the Ti6242s alloy. In fact, the relationship between stress, strain rate and temperature can be defined as σ = g ( ε ˙ , T ) . However, when considering the large nonlinear relationship among stress, strain rate and temperature, this function is not used. The effects of strain rate and temperature on stress were analyzed before proposing a new constitutive model. The data for low ( ε = 0.04 ), medium ( ε = 0.23 ) and high ( ε = 0.90 ) strain levels in Table 2 were selected for analysis. In order to study the relationship between logarithmic stress, logarithmic strain rate and temperature, the partial derivatives of logarithmic stress with respect to temperature and logarithmic strain rate can be calculated by discrete formula, that is, the difference quotient replacing the derivative. The discrete formulas of derivatives include forward difference, backward difference and central difference. In order to improve the accuracy, forward and backward difference are used at the boundary, and central difference is used in the middle region. Then, the nth-order partial derivatives of logarithmic stress with respect to the temperature and logarithmic strain rate for the experimental flow stress curves can be calculated by Equation (7).
{ n ln σ T n | i , j = n 1 ln σ T n 1 | i , j + 1   n 1 ln σ T n 1 | i , j 1 T i , j + 1 T i , j 1 n ln σ ( ln ε ˙ ) n | i , j = n 1 ln σ ( ln ε ˙ ) n 1 | i + 1 , j   n 1 ln σ ( ln ε ˙ ) n 1 | i 1 , j ln ε ˙ i + 1 , j ln ε ˙ i 1 , j
where n is the order of the partial derivative, i and j are indices of the strain rate and the temperature in the stress matrix, respectively (Table 3). For example, when i = 1 and j = 1, ln σ i , j represents the logarithmic stress with a strain rate of 0.001 s−1 and a temperature of 900 °C. The maximum values of i and j , which are both five in this study, are the number of strain rates and temperatures corresponding to the experiment, respectively. Equation (7) can be replaced by Equation (8) when i = 1 and j = 1.
{ n ln σ T n | i , j = n 1 ln σ T n 1 | i , j + 1   n 1 ln σ T n 1 | i , j T i , j + 1 T i , j n ln σ ( ln ε ˙ ) n | i , j = n 1 ln σ ( ln ε ˙ ) n 1 | i + 1 , j   n 1 ln σ ( ln ε ˙ ) n 1 | i , j ln ε ˙ i + 1 , j ln ε ˙ i , j
Equation (7) can be replaced by Equation (9) when i = 5 and j = 5.
{ n ln σ T n | i , j = n 1 ln σ T n 1 | i , j   n 1 ln σ T n 1 | i , j 1 T i , j T i , j 1 n ln σ ( ln ε ˙ ) n | i , j = n 1 ln σ ( ln ε ˙ ) n 1 | i , j   n 1 ln σ ( ln ε ˙ ) n 1 | i 1 , j ln ε ˙ i , j ln ε ˙ i 1 , j
Equation (8) can be replaced by Equation (10) when n = 1.
{ ln σ T | i , j = ln σ i , j + 1 ln σ i , j 1 T i , j + 1 T i , j 1 ln σ ln ε ˙ | i , j = ln σ i + 1 , j ln σ i 1 , j ln ε ˙ i + 1 , j ln ε ˙ i 1 , j
In fact, the greater the number of the temperatures and the strain rates in the experiment, the higher the calculation accuracy of Equations (7)–(10) is. As shown in Table 5 and Table 6, the first, second and third partial derivatives of logarithmic stress with respect to the temperature and logarithmic strain rate at low ( ε = 0.04 ), medium ( ε = 0.23 ) and high ( ε = 0.90 ) strain levels were calculated, using Equations (7)–(10). The redder the color in the figure, the smaller the value represented.
There are several phenomena that can be found from Table 5: (1) with the increase in temperature, the logarithmic stress has a monotonic decreasing characteristic for all strain levels (Table 5 a, e and i). (2) The difference of the first ( ln σ / T ) partial derivative of logarithmic stress with respect to the temperature is small for a certain strain level (Table 5 b, f, and j). (3) The second ( 2 ln σ / T 2 ) and third ( 3 ln σ / T 3 ) partial derivatives of logarithmic stress with respect to the temperature is close to zero for all strain levels (Table 5 c, g, k, d, h and l). According to calculus theory and the above phenomena, a basic conclusion can be derived; that is, the linear model can construct the relationship between logarithmic stress and temperature with high accuracy.
Similarly, there are several phenomena that can be found from Table 6 as follows: (1) with the increase in logarithmic strain rate, the logarithmic stress has a monotonic decreasing characteristic for all strain levels (see a, e and i in Table 6). (2) The difference in the first ( ln σ / ln ε ˙ ) partial derivative of logarithmic stress with respect to the logarithmic strain rate is relatively large for all strain levels (Table 6 b, f and j), and with the increase in temperature and logarithmic strain rate, the first ( ln σ / ln ε ˙ ) partial derivative is increased. (3) The second ( 2 ln σ / ( ln ε ˙ ) 2 ) partial derivative of logarithmic stress with respect to the logarithmic strain rate is relatively large for all strain levels (Table 6 c, g and k), and the third ( 3 ln σ / ( ln ε ˙ ) 3 ) partial derivative of logarithmic stress with respect to the logarithmic strain rate is very small (Table 6 d, h and l). According to calculus theory and the above phenomena, a basic conclusion can be derived; that is, the quadratic model can construct the relationship between logarithmic stress and logarithmic strain rate with high accuracy.
The function of logarithmic stress, logarithmic strain rate and temperature has been defined as ln σ = f ( ln ε ˙ , T ) . According to the binary Taylor expansion formula, f ( ln ε ˙ , T ) is expanded at ( ln ε ˙ = 0 , T = 0 ) to obtain Equation (11).
ln σ = i = 0 m [ 1 m ! n = 0 m ( C m n ( ln ε ˙ ) n T m n ( m f ( ln ε ˙ ) n T m n | ( 0 , 0 ) ) ) ]
where m is the number of terms, C m n is a combination operator and m is m-order partial derivative operator. Equation (11) shows that function f can be approximated by polynomials. The more complex the nonlinear relationship is, the more items are needed to build the same precision model. It is a linear model when m = 1, it is a quadratic model when m = 2, and so on. The greater the m is, the higher the accuracy is. It can be predicted that if the rheological data of the materials are relatively complex (such as including single-phase and two-phase microstructures, especially for α+β titanium alloy), m should be appropriately increased. According to Table 5 and Table 6, 3 ln σ / T 3 and 3 ln σ / ln ε ˙ 3 of Ti6242s alloy are both close to zero. It can be asserted that the quadratic model can be used to build a high-precision constitutive model f ( ln ε ˙ , T ) of this material.

4.2. Linear Model

First, the accuracy of the first-order model (linear model) was analyzed. When m = 1, the constitutive Equation (11) is simplified to Equation (12).
ln σ = k 0 + k 1 T + k 2 ln ε ˙
where k 0 , k 1 and k 2 are material parameters, which can be obtained by using multiple linear regression based on Table 3. The linear regression expression of material parameters corresponding to each strain is shown in Equation (13).
[ ln σ 1 ln σ 2 ln σ 3 ln σ w ] = [ 1 1 1 1   T 1 T 2 T 3 T w   ln ε ˙ 1 ln ε ˙ 2 ln ε ˙ 3 ln ε ˙ w   ] [ k 0 k 1 k 2 ] + μ 1
where w is the combined number of temperatures and strain rates when the strain is fixed (w = 5 × 5 in this study), and μ 1 is an error variable and follows the normal distribution with the mean value of zero. Ten groups of material parameters can be obtained by linear regression for each strain datum, which are shown as points in Figure 6.
The expression of each material parameter and strain can be obtained by fitting the data points of linear regression with a fifth-order polynomial. Then, the expression of each material is shown in Equation (14).
{ k 0 = 10.7957 + 26.6315 ε 124.7369 ε 2 + 261.1311 ε 3 285.3326 ε 4 + 97.1624 ε 5 k 1 = 0.0060 0.0249 ε + 0.1154 ε 2 0.2428 ε 3 + 0.2519 ε 4 0.0916 ε 5 k 2 = 0.1634 + 0.1248 ε 0.2586 ε 2 + 0.1393 ε 3 + 0.2671 ε 4 0.2520 ε 5
The linear constitutive equation of the titanium alloy can be obtained by bringing Equation (14) into Equation (12). The comparisons of prediction flow stress and experimental flow stress for the linear model are shown in Figure 7.
This model has a similar phenomenon to the Arrhenius model and HS model, that is, the prediction accuracy is poor in the low temperature region (900, 950 and 1000 °C) and higher in the high temperature region (1050 and 1100 °C). The model has the following advantages: fewer material parameters, a simple form and easy solution of material parameters (multiple linear regression only). In order to further improve the accuracy, a high-order constitutive model was proposed.

4.3. Quadratic Model

The linear model has low prediction accuracy for low temperature regions. Therefore, a quadratic model was proposed, in which the m is equal to two. The equation of the quadratic model is as follows:
ln σ = a 0 + a 1 ln ε ˙ + a 2 T + a 3 ( ln ε ˙ ) 2 + a 4 T 2 + a 5 T ln ε ˙
where a 0 ~ a 5 are material parameters, which can be obtained by multiple linear regression, based on Table 3. The linear regression expression of the material parameters corresponding to each strain is shown in Equation (16).
[ ln σ 1 ln σ 2 ln σ 3 ln σ w ] = [ 1 1 1 1   ln ε ˙ 1 ln ε ˙ 2 ln ε ˙ 3 ln ε ˙ w   T 1 T 2 T 3 T w   ( ln ε ˙ 1 ) 2 ( ln ε ˙ 2 ) 2 ( ln ε ˙ 3 ) 2 ( ln ε ˙ w ) 2   T 1 2 T 2 2 T 3 2 T w 2   T 1 ln ε ˙ 1 T 2 ln ε ˙ 2 T 3 ln ε ˙ 3 T w ln ε ˙ w   ] [ a 0   a 1   a 2   a 3   a 4   a 5 ] ' + μ 2
where w is the combined number of temperatures and strain rates when the strain is fixed (w = 5 × 5 in this study), and μ 2 is an error variable and follows the normal distribution with the mean value of zero. Ten groups of material parameters can be obtained by linear regression for each strain data, which are shown as points in Figure 8.
The expression of each material parameter and strain can be obtained by fitting the data points of linear regression with a fifth-order polynomial. Then, the expression of each material is shown in Equation (17).
{ a 0 = 1.58977 + 243.78788 ε 1176.81606 ε 2 + 2530.31159 ε 3 2516.33047 ε 4 + 939.51775 ε 5 a 1 = 0.217555 0.222702 ε + 2.174965 ε 2 2.618603 ε 3 1.456050 ε 4 + 2.177008 ε 5 a 2 = 0.011670 0.461668 ε + 2.234179 ε 2 4.807239 ε 3 + 4.773571 ε 4 1.777977 ε 5 a 3 = 0.004285 0.022862 ε 0.075039 ε 2 0.146988 ε 3 + 0.141039 ε 4 0.053276 ε 5 a 4 = 0.000008 + 0.000218 ε 0.001061 ε 2 + 0.002284 ε 3 0.002263 ε 4 + 0.0008400 ε 5 a 5 = 0.000361 + 0.000242 ε 0.002088 ε 2 + 0.002081 ε 3 + 0.002373 ε 4 0.002674 ε 5
Similarly, we inserted (17) into (15) to obtain the constitutive equation of the quadratic model, and compared the predicted values of the model with the experimental values to give Figure 9.
Compared with the linear model (see Figure 7 and Figure 9), the prediction accuracy of the quadratic model was improved significantly. In addition, the prediction accuracy of the quadratic model is significantly higher than that of the HS model and the Arrhenius model. Contrary to the previous models, the prediction accuracy of the model is high in both the single-phase region (900, 950 and 1000 °C) and the two-phase region (1050 and 1100 °C).

5. Predictability Analysis

5.1. Prediction Accuracy

5.1.1. Overall Accuracy

In order to analyze the overall prediction accuracy, the accuracy of the classical models and the new models was quantified in terms of standard statistical parameters such as correlation coefficient (R), root mean square error (RMSE), sum of squares for error (SSE), and sum of absolute error (SAE). The formulas of these parameters are shown in Equations (18)–(21).
R = i = 1 N ( σ ^ i σ ^ ¯ ) ( σ i σ ¯ ) i = 1 N ( σ ^ i σ ^ ¯ ) 2 i = 1 N ( σ i σ ¯ ) 2
R M S E = i = 1 N ( σ ^ i σ i ) 2 N
S S E = i = 1 N ( σ ^ i σ i ) 2
S A E = i = 1 N | σ ^ i σ i |
where σ ^ i is the ith prediction value of the flow stress, σ i is the ith experiment value of flow stress, σ ^ ¯ is the mean value of the prediction flow stress, σ ¯ is the mean value of the experiment flow stress, and N is the number of comparison points (N = 5 × 5 × 10). The quantitative comparison of the prediction accuracy of the different models is listed in Table 7.
As can be seen from Table 7, the order of prediction accuracy from high to low is the quadratic model, Arrhenius model, HS model and linear model. The SSE of the quadratic model (1.0E3) is one order of magnitude smaller than the other models (1.0E5). The RMSE of the quadratic model is 1/2 of the Arrhenius model, 1/3 of the HS model and 1/4 of the linear model. The correlation coefficients of the quadratic model are significantly higher than the other models. The sum of absolute error of the quadratic model is 884.17 MPa, which is much less than the other models (1422.53 MPa for the Arrhenius model, 2097.60 MPa for the HS model and 2827.67 MPa for the linear model). Therefore, the overall accuracy of the quadratic model is significantly higher than the other models.
As can be seen from Figure 10, the data points of the quadratic model are distributed on both sides of the centerline and more concentrated than the other models. In addition, the square of the correlation coefficient between experimental stress and predicted stress in the quadratic model, Arrhenius model and Hensel–Spittel model are 0.9961, 0.988 and 0.9751, respectively. This also shows that the prediction accuracy of the quadratic model is higher than that of the other models.
As shown in Figure 11, the “” symbol represents the predictive value of the Arrhenius model, the “ ” symbol represents the predictive value of the linear model, the “×” represents the predictive value of the HS model, the “+” represents the predictive value of the quadratic model, and the solid line represents the experimental curve. Almost all the prediction points of the quadratic model are closer to the experimental curve than the other models. The precision of the quadratic model is significantly higher than that of the other models.

5.1.2. Local Accuracy

In order to analyze the prediction accuracy under different temperatures and strain rates, the correlation coefficient (R) and sum of absolute error (SAE) were calculated. The distribution of the correlation coefficient between predicted flow stress and experimental flow stress on temperature and strain rate is shown in Figure 12.
As shown in Figure 12, the difference from high to low in the correlation coefficient with temperature gives the ranking order as the HS model, Arrhenius model, linear model and quadratic model. The difference between the maximum value and the minimum value of the correlation coefficient for the quadratic model was 0.073, which was significantly less than that of the other models (0.2486 for HS model, 0.1228 for linear model and 0.1199 for Arrhenius model). The rules in Figure 12 show that the prediction accuracy of Arrhenius model and Hensel–Spittel model in the two phase region ( α + β ) and single phase ( α ) is significantly different. The prediction accuracy of the quadratic model remains high in both the single phase region and the two phase region, and the prediction accuracy of the quadratic model is significantly higher than the other models. The prediction accuracy of the linear model in the single phase region and the two phase region is also similar to the quadratic model, but the prediction accuracy is lower. In addition, the distribution of the correlation coefficient on the strain rate for all models also has the above rules, that is, the distribution of the correlation coefficient on strain rate for the Hensel–Spittel and Arrhenius models is very different. The correlation coefficient of the quadratic model at different strain rates and temperatures is large, and the difference is small.
The distribution of the sum of absolute error between predicted flow stress and experimental flow stress on temperature and strain rate is shown in Figure 13.
As shown in Figure 13, the maximum absolute error sums of the Arrhenius model, HS model, linear model and quadratic model are 237.4, 274.3, 727.2, and 95.53 MPa, respectively. This also shows that the accuracy of the quadratic model is significantly higher than that of the other models. In addition, the maximum and minimum values of the sum of absolute error for the Arrhenius model, HS model and linear model differ greatly, which indicates that the prediction accuracy of these three models for specific strain rate and temperature is poor. The distribution of the absolute error sum of the quadratic model is relatively uniform, which shows that the quadratic model has good prediction ability for all temperatures and strain rates.

5.2. Predictability of Model

Generally, linear and non-linear multiple regression can be used to obtain the global optimal material parameters for classical models and the linear model. However, the classical models and the linear model still produce large errors. Therefore, the reason for the large error could be that the classical models and the linear model have insufficient predictability. In this section, mathematical means are used to study the reasons for the insufficient precision of the classical models and the linear model.

5.2.1. Analytic Formula of Partial Differential

It can be seen from Table 5 and Table 6 that the third derivative of the logarithmic stress with respect to the logarithmic strain rate of the material is close to zero, and the second derivative of the logarithmic stress with respect to the temperature is close to zero. Therefore, the analytical partial derivatives of all models are derived first. The first-, second- and third-order partial derivatives of the logarithmic stress with respect to the logarithmic strain rate and temperature for all constitutive models can be solved by Equation (22).
{ ln σ ln ε ˙ = ε ˙ σ σ ε ˙ ln σ T = σ σ T 2 ln σ ( ln ε ˙ ) 2 = ε ˙ ( ln σ ln ε ˙ ) σ ε ˙ 2 ln σ T 2 = ( ln σ T ) σ T     3 ln σ ( ln ε ˙ ) 3 = ε ˙ ( 2 ln σ ( ln ε ˙ ) 2 ) σ ε ˙ 3 ln σ T 3 = ( 2 ln σ T 2 ) σ T      
σ = { 1 α ln ( ( 1 A ε ˙ exp ( Q R T ) ) 1 n + ( 1 A ε ˙ exp ( Q R T ) ) 2 n + 1 ) ,   for   Arrhenius   Model Aexp ( m 1 T ) ε m 2 ε ˙ m 3 exp ( m 4 ε ) ( 1 + ε ) m 5 T exp ( m 6 ε ) ε ˙ m 7 T T m 8 ,   for   HS   Model   exp ( k 0 + k 1 T + k 2 ln ε ˙ ) ,   for   Linear   Model   exp ( a 0 + a 1 ln ε ˙ + a 2 T + a 3 ( ln ε ˙ ) 2 + a 4 T 2 + a 5 T ln ε ˙ ) ,   for   Quadratic   Model
where A , α , n and Q are quintic polynomials of ε shown in Equation (3), ln A and m 1 m 8 are constants shown in Table 4. k 0 k 1 are quintic polynomials of ε shown in Equation (14) and a 0 a 5 are quintic polynomials of ε shown in Equation (17). When the strain value is fixed, these parameters are constants.

5.2.2. Results of Partial Differential

The first, second and third partial derivatives of logarithmic stress with respect to temperature for experimental data and different models are listed in Table 8 when strain is equal to 0.04 by using Equations (22) and (23). The redder the color in the figure, the smaller the value represented.
There are several rules that can be drawn from Table 8: (1) the second and third partial derivatives of all models are close to zero, which are in good agreement with the experimental data (Table 8 c, g, k, o, s, d, h, l, p and t). (2) The first partial derivatives for the Arrhenius model and the HS model have great regularity with temperature. The partial differential decreases with the increase in temperature for the Arrhenius model and the law is the opposite for the HS model (Table 8 f and j). (3) The first partial derivative for the linear model is not related to temperature or strain rate (Table 8 n). (4) The first partial derivative for the quadratic model is less dependent on temperature, but strongly dependent on strain rate (Table 8 r). (5) The first partial derivative for the experimental data is less dependent on temperature, but slightly dependent strain rate (Table 8 b). The above phenomena show that the experimental data can be approximated well only by one order accuracy, and all models meet the first-order accuracy. However, the quadratic model has second-order accuracy between logarithmic stress and temperature. Considering that all models have first-order accuracy for temperature, another error source is the strain rate factor.
The first, second and third partial derivatives of logarithmic stress with respect to logarithmic strain rate for experimental data and different models are listed in Table 9, when strain is equal to 0.04, by using Equations (22) and (23). The redder the color in the figure, the smaller the value represented. It can be seen from Table 9 that the second-order partial derivatives of logarithmic stress with respect to logarithmic strain rate for the Arrhenius model, the HS model and the linear model are close to zero. It is worth noting that the third partial differential of logarithmic stress with respect to logarithmic strain rate for the quadratic model is close to zero. This phenomenon is determined by the equation of the constitutive model and is independent of the regression error, because the partial derivatives are calculated by analytical formulas. However, the second derivative of logarithmic stress with respect to logarithmic strain rate for the experimental data is not close to zero. This proves that the Arrhenius model, the HS model and the linear model cannot make a second-order approximation to the experimental data at logarithmic strain rate, while the quadratic model can.
By comparing the first, second and third partial derivatives of all models, the following conclusions can be drawn: (1) the predictabilities of the classical models are only the first-order accuracy in the logarithmic stress, logarithmic strain rate and temperature space. (2) In essence, the predictabilities of the Arrhenius model, the HS model and the linear model are in the same order of magnitude, because they are linear approximations of the logarithmic stress, logarithmic strain rate and temperature space. (3) The predictability of the quadratic model is one level higher than that of the Arrhenius model, the HS model and the linear model, because it is a quadratic approximation. (4) The traditional modification methods for the Arrhenius and HS models cannot significantly improve the predictability of the models, because these model modifications are carried out within the linear level in the logarithmic stress, logarithmic strain rate and temperature space. (5) Due to the existence of the cross term of temperature and logarithmic strain rate ( T ln ε ˙ ) and the quadratic term of temperature ( T 2 ), the quadratic model can largely eliminate the predictability difference between the classical model and the linear model at different temperatures.

6. Conclusions

In this paper, hot compression tests under different strain rates and temperatures were carried out. A new constitutive model based on the Taylor series and partial derivatives of the experimental data is proposed without significantly increasing the material parameters. The predictabilities of the new model and the classical models (Arrhenius and HS models) were analyzed. The reason for the lack of predictability in the classical models is proved mathematically. The main conclusions are as follows:
(1) The prediction accuracies of the Arrhenius model and HS model are significantly different in different temperatures and with different strain rates. The prediction accuracies of the Arrhenius model and the HS model are low in the single-phase region (900, 950 and 100 °C) and high in the two-phase region (1050 and 1100 °C). The prediction accuracy of the quadratic model in all temperatures and with all strain rates has no significant difference, and is higher than the other models. The order of prediction accuracy from high to low is the quadratic model, the Arrhenius model, the HS model and the linear model. The SSE of the quadratic model (1.0 × 103) is one order of magnitude smaller than the other models (1.0 × 105). The RMSE of the quadratic model is 1/2 of the Arrhenius model, 1/3 of the HS model and 1/4 of the linear model. The correlation coefficients of the quadratic model are significantly high than the other models. The sum of absolute error of the quadratic model is 884.17 MPa which is much less than the other models (1422.53 MPa for the Arrhenius model, 2097.60 MPa for the HS model and 2827.67 MPa for the linear model). The square of the correlation coefficient between experimental stress and predicted stress in the quadratic model, Arrhenius model and Hensel–Spittel model is 0.9961, 0.988 and 0.9751, respectively. This also shows that the prediction accuracy of the quadratic model is higher than that of the other models.
(2) The second and third partial derivatives of logarithmic stress with respect to temperature for the experimental data is close to zero. The third partial derivatives of logarithmic stress with respect to logarithmic strain rate for the experimental data is close to zero. This phenomena indicates that second-order approximation is required between logarithmic stress and logarithmic strain rate, and first-order approximation is required between logarithmic stress and temperature to establish a high-precision model. The Arrhenius model, HS model and linear model meet the first-order approximation requirements between logarithmic stress and temperature, but do not meet the second-order approximation requirements between logarithmic stress and logarithmic strain rate. The quadratic model meets these two requirements and thus has higher prediction accuracy than the other models.
(3) The order of prediction accuracy for the four models from high to low is the quadratic model, Arrhenius model, linear model and HS model. The material parameters of the quadratic model can be solved only by multiple linear regression, while the Arrhenius model needs to use multiple nonlinear regression. Compared with the Arrhenius model, the prediction accuracy of the quadratic model can be improved significantly by adding only a few parameters, and the parameter solution is simpler. The prediction accuracy of the new model can be further improved by adding the number of items.

Author Contributions

J.Z. (Jiansheng Zhang) and G.X. put forward the relevant algorithms and ideas, and wrote the article. G.D. and Y.Z. carried out the experiments. J.Z. (Jie Zhou) guided and reviewed the whole process. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the National Natural Science Foundation of China (No. 52205337), General Project of Chongqing Natural Science Foundation (No.CSTB2022NSCQ-MSX0443) and Science and Technology Planning Project of Deyang (No.2021JBJZ008).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

All authors gave their consent for the manuscript to be published.

Data Availability Statement

All data generated or analyzed during this study are included in this published article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cui, Y.M.; Zheng, W.W.; Li, C.H.; Cao, G.H.; Wang, Y.D. Effectiveness of hot deformation and subsequent annealing for β grain refinement of Ti-5Al-5Mo-5V-1Cr-1Fe titanium alloy. Rare Met. 2021, 40, 3608–3615. [Google Scholar] [CrossRef]
  2. Li, C.M.; Huang, L.; Li, C.L.; Hui, S.X.; Yu, Y.; Zhao, M.J.; Guo, S.Q.; Li, J.J. Research progress on hot deformation behavior of high-strength β titanium alloy: Flow behavior and constitutive model. Rare Met. 2022, 41, 1434–1455. [Google Scholar] [CrossRef]
  3. Xu, W.X.; Zhao, Z.L.; Liu, N.; Huang, X.Y.; Li, S.; Wang, Y.B. Dynamic transformation of a near alpha high-temperature titanium alloy during hot deformation. IOP Cof. Ser. Mater. Sci. Eng. 2022, 1270, 012108. [Google Scholar] [CrossRef]
  4. Ye, B.; Li, L.X.; Liu, S.; Hu, S.D.; Liao, H.Q. A unified constitutive equation of a bainite steel during hot deformation. J. Mater. Eng. Perform. 2017, 25, 4581–4586. [Google Scholar] [CrossRef]
  5. Haghdadi, N.; Zarei-Hanzaki, A.; Abedi, H.R. The flow behavior modeling of cast A356 aluminum alloy at elevated temperatures considering the effect of strain. Mater. Sci. Eng. A 2012, 535, 252–257. [Google Scholar] [CrossRef]
  6. Marandi, A.; Zarei-Hanzaki, A.; Haghdadi, N.; Eskansdari, M. The prediction of hot deformation behavior in Fe-21Mn-2.5Si-1.5Al transformation-twinning induced plasticity steel. Mater. Sci. Eng. A 2012, 554, 72–78. [Google Scholar] [CrossRef]
  7. Zhang, H.; Wen, W.; Cui, H.; Xu, Y. A modified Zerilli-Armstrong model for alloy IC10 over a wide range of temperatures and strain rates. Mater. Sci. Eng. A 2009, 527, 328–333. [Google Scholar] [CrossRef]
  8. Voyiadjis, G.Z.; Almasri, A.H. A physically based constitutive model for fcc metals with applications to dynamic hardness. Mech. Mater. 2008, 40, 549–563. [Google Scholar] [CrossRef]
  9. Lin, Y.C.; Chen, X.M. A critical review of experimental results and constitutive descriptions for metals and alloys in hot working. Mater. Des. 2011, 32, 1733–1759. [Google Scholar] [CrossRef]
  10. Sellars, C.M.; Mctegart, W.J. On the mechanism of hot deformation. Acta Metall. 1996, 14, 1136–1138. [Google Scholar] [CrossRef]
  11. Zener, C.; Hollomon, J.H. Effect of strain rate upon plastic flow of steel. J. Appl. Phys. 1994, 15, 22–32. [Google Scholar] [CrossRef]
  12. Johnson, G.R.; Cook, W.H. A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures. Eng. Fract. Mech. 1983, 21, 541–548. [Google Scholar]
  13. Xu, Z.J.; Huang, F.L. Thermomechanical behavior and constitutive modeling of tungsten-based composite over wide temperature and strain rate ranges. Int. J. Plast. 2013, 40, 163–184. [Google Scholar] [CrossRef]
  14. Ying, L.; Gao, T.H.; Rong, H.; Han, X.; Hu, P.; Hou, W.B. On the thermal forming limit diagram (TFLD) with GTN mesoscopic damage model for AA7075 aluminum alloy: Numerical and experimental investigation. J. Alloys Compd. 2019, 802, 675–693. [Google Scholar] [CrossRef]
  15. Kotkunde, N.; Krushnamurthy, H.N.; Puranik, P.; Gupta, A.K. Microstructure study and constitutive modeling of Ti–6Al–4V alloy at elevated temperatures. Mater. Des. 2014, 54, 96–103. [Google Scholar] [CrossRef]
  16. Zhang, C.; Li, X.Q.; Li, D.S.; Jin, C.H.; Xiao, J.J. Modelization and comparison of Norton-Hoff and Arrhenius constitutive laws to predict hot tensile behavior of Ti–6Al–4V alloy. Trans. Nonferrous Met. Soc. China 2012, 22, s457–s464. [Google Scholar] [CrossRef]
  17. Xiao, M.; Li, F.; Wei, Z.; Yang, G.L. Constitutive equation for elevated temperature flow behavior of TiNiNb alloy based on orthogonal analysis. Mater. Des. 2012, 35, 184–193. [Google Scholar] [CrossRef]
  18. Luo, J.; Li, M.; Li, X.; Shi, Y.P. Constitutive model for high temperature deformation of titanium alloys using internal state variables. Mech. Mater. 2010, 42, 157–165. [Google Scholar] [CrossRef]
  19. Cai, J.; Shi, J.M.; Wang, K.S.; Li, F.G.; Wang, W.; Wang, Q.J.; Liu, Y.Y. A modified parallel constitutive model for elevated temperature flow behavior of Ti-6Al-4V alloy based on multiple regression. Int. J. Mater. Res. 2017, 108, 527–541. [Google Scholar] [CrossRef]
  20. Ming, H.; Dong, L.; Zhang, Z.; Lei, X.F.; Yang, R.; Sha, Y.H. Correction of Flow Curves and Constitutive Modelling of a Ti-6Al-4V Alloy. Metals 2018, 8, 256. [Google Scholar] [CrossRef] [Green Version]
  21. Gao, S.; Sang, Y.; Li, Q.; Sun, Y.; Wu, Y.; Wang, H. Constitutive modeling and microstructure research on the deformation mechanism of Ti-6Al-4V alloy under hot forming condition. J. Alloys Compd. 2022, 892, 162128. [Google Scholar] [CrossRef]
  22. Ahmed, M.; Anastasia, M.; Anton, K.; Theo, P.; Sergey, A.; James, K.; Vladimir, P. Modelling of the Superplastic Deformation of the Near-α Titanium Alloy (Ti-2.5Al-1.8Mn) Using Arrhenius-Type Constitutive Model and Artificial Neural Network. Metals 2017, 7, 568. [Google Scholar] [CrossRef] [Green Version]
  23. Sun, Y.; Zeng, W.D.; Zhao, Y.Q.; Zhang, X.M.; Shu, Y.; Zhou, Y.G. Modeling constitutive relationship of Ti40 alloy using artificial neural network. Mater. Des. 2011, 32, 1537–1541. [Google Scholar] [CrossRef]
  24. Reddy, N.S.; Lee, Y.H.; Park, C.H.; Lee, C.S. Prediction of flow stress in Ti–6Al–4V alloy with an equiaxed α + β microstructure by artificial neural networks. Mater. Sci. Eng. A 2008, 492, 276–282. [Google Scholar] [CrossRef]
  25. Prasad, Y.V.R.K.; Gegel, H.L.; Doraivelu, S.M.; Malas, J.C.; Morgan, J.T.; Lark, K.A.; Barker, D.R. Modeling of dynamic material behavior in hot deformation: Forging of Ti-6242. Metall. Trans. A 1984, 15, 1883–1892. [Google Scholar] [CrossRef]
  26. Ren, G.G.; Je, Y.T.; Zhan, Y.Z.; Shong, S.L. Effect of microstructure on deformation behavior of Ti-6Al-4V alloy during compressing process. Mater. Des. 2012, 36, 796–803. [Google Scholar] [CrossRef]
  27. Huang, S.H.; Zong, Y.Y.; Shan, D.B. Application of thermohydrogen pr°Cessing to Ti6Al4V alloy blade isothermal forging. Mater. Sci. Eng. A 2013, 561, 17–25. [Google Scholar] [CrossRef]
  28. Comley, P. Multi-rate Superplastic Forming of Fine Grain Ti-6Al-4V Titanium Alloy. J. Mater. Eng. Perform. 2007, 16, 150–154. [Google Scholar] [CrossRef]
  29. Richardson, G.J.; Sellars, C.M.; Tegart, W. Recrystallization during creep of nickel. Acta Metall. 1996, 14, 1225–1236. [Google Scholar] [CrossRef]
  30. Mehtedi, M.E.; Spigarelli, S.; Gabrielli, F.; Donati, L. Comparison Study of Constitutive Models in Predicting the Hot Deformation Behavior of AA6060 and AA6063 Aluminium Alloys. Mater. Today Proc. 2015, 2, 4732–4739. [Google Scholar] [CrossRef]
  31. Mehtedi, M.E.; Musharavati, F.; Spigarelli, S. Modelling of the flow behaviour of wrought aluminium alloys at elevated temperatures by a new constitutive equation. Mater. Des. 2014, 54, 869–873. [Google Scholar] [CrossRef]
Figure 1. The microstructure of the Ti6242s alloy.
Figure 1. The microstructure of the Ti6242s alloy.
Materials 16 02928 g001
Figure 2. Flow stress–strain curves of Ti6242s alloy during hot compression at strain rates of (a) 0.01 s−1; (b) 0.1 s−1; (c) 1 s−1; (d) 10 s−1; (e) 10 s−1.
Figure 2. Flow stress–strain curves of Ti6242s alloy during hot compression at strain rates of (a) 0.01 s−1; (b) 0.1 s−1; (c) 1 s−1; (d) 10 s−1; (e) 10 s−1.
Materials 16 02928 g002
Figure 3. Nonlinear regression values of material parameters at different strain levels and their polynomial fitting: (a) lnA-ε; (b) α-ε; (c) n-ε; (d) Q-ε.
Figure 3. Nonlinear regression values of material parameters at different strain levels and their polynomial fitting: (a) lnA-ε; (b) α-ε; (c) n-ε; (d) Q-ε.
Materials 16 02928 g003
Figure 4. Comparison of prediction flow stress and experimental flow stress for Arrhenius model: (a) 0.01 s−1; (b) 0.1 s−1; (c) 1 s−1; (d) 10 s−1; (e) 10 s−1.
Figure 4. Comparison of prediction flow stress and experimental flow stress for Arrhenius model: (a) 0.01 s−1; (b) 0.1 s−1; (c) 1 s−1; (d) 10 s−1; (e) 10 s−1.
Materials 16 02928 g004
Figure 5. Comparison of prediction flow stress and experimental flow stress for HS model: (a) 0.01 s−1; (b) 0.1 s−1; (c) 1 s−1; (d) 10 s−1; (e) 10 s−1.
Figure 5. Comparison of prediction flow stress and experimental flow stress for HS model: (a) 0.01 s−1; (b) 0.1 s−1; (c) 1 s−1; (d) 10 s−1; (e) 10 s−1.
Materials 16 02928 g005
Figure 6. Linear regression values of material parameters at different strain levels and their polynomial fitting: (a) Material parameter of k 0 ; (b) Material parameter of k 1 ; (c) Material parameter of k 2 .
Figure 6. Linear regression values of material parameters at different strain levels and their polynomial fitting: (a) Material parameter of k 0 ; (b) Material parameter of k 1 ; (c) Material parameter of k 2 .
Materials 16 02928 g006
Figure 7. Comparison of prediction flow stress and experimental flow stress for linear model: (a) 0.01 s−1; (b) 0.1 s−1; (c) 1 s−1; (d) 10 s−1; (e) 10 s−1.
Figure 7. Comparison of prediction flow stress and experimental flow stress for linear model: (a) 0.01 s−1; (b) 0.1 s−1; (c) 1 s−1; (d) 10 s−1; (e) 10 s−1.
Materials 16 02928 g007
Figure 8. Linear regression values of material parameters at different strain levels and their polynomial fitting for the quadratic model.
Figure 8. Linear regression values of material parameters at different strain levels and their polynomial fitting for the quadratic model.
Materials 16 02928 g008
Figure 9. Comparison of prediction flow stress and experimental flow stress for the quadratic model: (a) 0.01 s−1; (b) 0.1 s−1; (c) 1 s−1; (d) 10 s−1; (e) 10 s−1.
Figure 9. Comparison of prediction flow stress and experimental flow stress for the quadratic model: (a) 0.01 s−1; (b) 0.1 s−1; (c) 1 s−1; (d) 10 s−1; (e) 10 s−1.
Materials 16 02928 g009
Figure 10. Predictability of flow stress by (a) Arrhenius model; (b) HS model; (c) linear Model; (d) quadratic model.
Figure 10. Predictability of flow stress by (a) Arrhenius model; (b) HS model; (c) linear Model; (d) quadratic model.
Materials 16 02928 g010
Figure 11. Comparison of prediction flow stress and experimental flow stress for Arrhenius, HS, linear and quadratic models: (a) 0.001 s−1; (b) 0.01 s−1; (c) 0.1 s−1; (d) 1 s−1; (e) 10 s−1.
Figure 11. Comparison of prediction flow stress and experimental flow stress for Arrhenius, HS, linear and quadratic models: (a) 0.001 s−1; (b) 0.01 s−1; (c) 0.1 s−1; (d) 1 s−1; (e) 10 s−1.
Materials 16 02928 g011
Figure 12. The distribution of the correlation coefficient between predicted flow stress and experimental flow stress on temperature and strain rate: (a) Arrhenius model; (b) HS model; (c) linear model; (d) quadratic model.
Figure 12. The distribution of the correlation coefficient between predicted flow stress and experimental flow stress on temperature and strain rate: (a) Arrhenius model; (b) HS model; (c) linear model; (d) quadratic model.
Materials 16 02928 g012
Figure 13. The distribution of the sum of absolute error between predicted flow stress and experimental flow stress on temperature and strain rate: (a) Arrhenius constitutive model, (b) HS constitutive model, (c) linear model, and (d) quadratic model.
Figure 13. The distribution of the sum of absolute error between predicted flow stress and experimental flow stress on temperature and strain rate: (a) Arrhenius constitutive model, (b) HS constitutive model, (c) linear model, and (d) quadratic model.
Materials 16 02928 g013
Table 1. The normal chemical composition of Ti6242s alloy (wt, %).
Table 1. The normal chemical composition of Ti6242s alloy (wt, %).
AlSnZrMoSiFeONC
6.131.934.161.900.0900.0200.090.010.010
Table 2. Parameters used in isothermal compression test.
Table 2. Parameters used in isothermal compression test.
NumberHeating Temperature (T/°C)Strain Rate ( ε ˙ /s−1) Heating Rate (°C/s)Quenching MediumAmount of Deformation
1–59000.001, 0.01, 0.1, 1, 105water60%
6–10950
11–151000
16–201050
21–251100
Table 3. Stress matrix corresponding to different strain levels (The units of strain rate, temperature and stress are s−1, K and MPa, respectively).
Table 3. Stress matrix corresponding to different strain levels (The units of strain rate, temperature and stress are s−1, K and MPa, respectively).
Temperature Temperature
ε ε ˙ 11731223127313231373 ε ε ˙ 11731223127313231373
0.040.00188.2761.9433.4023.6815.630.140.00190.8660.4532.9922.9314.75
0.01107.1993.1055.7939.5327.930.01153.4789.1656.7339.2224.94
0.10180.23127.2988.2064.5641.710.10222.10135.2789.1267.5040.31
1216.19190.22126.7891.5971.241279.70192.35122.4894.1269.67
10287.08230.64162.64128.04104.7810347.07246.11169.92132.68105.52
ε ε ˙ 11731223127313231373 ε ε ˙ 11731223127313231373
0.230.00185.7956.5431.4921.9313.720.330.00180.8252.7429.5421.0312.68
0.01144.7183.4653.6938.2723.900.01136.6178.1751.7537.3423.37
0.10217.10133.4785.3265.7938.910.10205.28130.5282.2863.5037.47
1276.49190.86118.6891.5367.651262.37186.65117.1288.0065.71
10357.03239.83166.90128.26102.9210350.95233.02160.30122.5397.28
ε ε ˙ 11731223127313231373 ε ε ˙ 11731223127313231373
0.420.00175.5449.0827.1020.0712.220.520.00172.4445.6225.5919.6012.22
0.01130.2774.1650.5536.1823.110.01122.7770.1448.2134.6122.85
0.10195.80127.8879.3761.7635.640.10186.98125.0075.9360.1934.22
1250.85183.04115.5384.0662.751240.83179.25112.5080.5961.77
10336.61225.79155.49120.1494.0710316.86216.11151.21118.2892.16
ε ε ˙ 11731223127313231373 ε ε ˙ 11731223127313231373
0.610.00169.5643.2723.8818.7611.310.710.00167.5441.5122.1518.2710.59
0.01116.4366.8045.7733.8322.060.01110.9464.3945.0432.7921.38
0.10178.63122.6873.5958.5733.410.10171.54119.8571.9658.1733.00
1231.14175.09110.7579.2861.071222.91169.64109.7678.2560.79
10299.39210.62148.26117.0691.5510283.90206.85145.62116.9892.55
ε ε ˙ 11731223127313231373 ε ε ˙ 11731223127313231373
0.800.00166.3140.4421.2517.7210.460.900.00165.7039.6120.6117.4510.42
0.01106.6362.3242.6132.7820.840.01104.1261.0142.4533.3819.86
0.10167.04117.4272.3759.3633.000.10165.31118.2372.0058.8532.19
1216.97163.66108.3178.0161.281212.39160.23106.8477.4460.30
10270.08200.75143.65115.2993.9810260.68194.34140.07109.9692.97
Table 4. Multiple linear regression results of material parameters.
Table 4. Multiple linear regression results of material parameters.
ln A m 1 m 2 m 3 m 4 m 5 m 6 m 7 m 8
120.33770.0069−0.1544−0.2426−0.01350.0001−0.15810.0003−17.4000
Table 5. The first, second and third partial derivatives of logarithmic stress with respect to temperature at low ( ε = 0.04 ), medium ( ε = 0.23 ) and high ( ε = 0.90 ) strain levels for the experiments’ flow stress curves.
Table 5. The first, second and third partial derivatives of logarithmic stress with respect to temperature at low ( ε = 0.04 ), medium ( ε = 0.23 ) and high ( ε = 0.90 ) strain levels for the experiments’ flow stress curves.
Order of Partial Derivatives of Logarithmic Stress with Respect to Temperature
εln σ ln σ / T 2 ln σ / T 2 3 ln σ / T 3
0.04(a)(b)(c)(d)
Materials 16 02928 i001Materials 16 02928 i002Materials 16 02928 i003Materials 16 02928 i004
0.23(e)(f)(g)(h)
Materials 16 02928 i005Materials 16 02928 i006Materials 16 02928 i007Materials 16 02928 i008
0.9(i)(j)(k)(l)
Materials 16 02928 i009Materials 16 02928 i010Materials 16 02928 i011Materials 16 02928 i012
Table 6. The first, second and third partial derivatives of logarithmic stress with respect to the variable logarithmic strain rate at low ( ε = 0.04 ), medium ( ε = 0.23 ) and high ( ε = 0.90 ) strain levels for experiment flow stress curves.
Table 6. The first, second and third partial derivatives of logarithmic stress with respect to the variable logarithmic strain rate at low ( ε = 0.04 ), medium ( ε = 0.23 ) and high ( ε = 0.90 ) strain levels for experiment flow stress curves.
Order of Partial Derivative of Logarithmic Stress with Respect to the Variable Logarithmic Strain Rate
Strainln σ ln σ / ln ε ˙ 2 ln σ / ( ln ε ˙ ) 2 3 ln σ / ( ln ε ˙ ) 3
0.04(a)(b)(c)(d)
Materials 16 02928 i013Materials 16 02928 i014Materials 16 02928 i015Materials 16 02928 i016
0.23(e)(f)(g)(h)
Materials 16 02928 i017Materials 16 02928 i018Materials 16 02928 i019Materials 16 02928 i020
0.9(i)(j)(k)(l)
Materials 16 02928 i021Materials 16 02928 i022Materials 16 02928 i023Materials 16 02928 i024
Table 7. Quantitative comparison of predictability of different models.
Table 7. Quantitative comparison of predictability of different models.
IndexsArrhenius ModelHS ModelLinear ModelQuadratic Model
R0.99390.98740.98140.9980
RMSE8.43312.60916.8974.3737
SSE1764039430708105565
SAE (MPa)1422.522097.602827.67884.17
Table 8. The first, second and third partial derivative of logarithmic stress with respect to temperature for different models at ε = 0.04 .
Table 8. The first, second and third partial derivative of logarithmic stress with respect to temperature for different models at ε = 0.04 .
.Order of Partial Derivatives of Logarithmic Stress with Respect to Temperature
Modelln σ ln σ / T 2 ln σ / T 2 3 ln σ / T 3
Experimental data(a)(b)(c)(d)
Materials 16 02928 i025Materials 16 02928 i026Materials 16 02928 i027Materials 16 02928 i028
Arrhenius model(e)(f)(g)(h)
Materials 16 02928 i029Materials 16 02928 i030Materials 16 02928 i031Materials 16 02928 i032
HS Model(i)(j)(k)(l)
Materials 16 02928 i033Materials 16 02928 i034Materials 16 02928 i035Materials 16 02928 i036
Linear model(m)(n)(o)(p)
Materials 16 02928 i037Materials 16 02928 i038Materials 16 02928 i039Materials 16 02928 i040
Quadratic model(q)(r)(s)(t)
Materials 16 02928 i041Materials 16 02928 i042Materials 16 02928 i043Materials 16 02928 i044
Table 9. The first, second and third partial derivatives of logarithmic stress with respect to strain rate for different models at ε = 0.04 .
Table 9. The first, second and third partial derivatives of logarithmic stress with respect to strain rate for different models at ε = 0.04 .
.Order of Partial Derivatives of Logarithmic Stress with Respect to Temperature
Model ln σ ln σ / ln ε ˙ 2 ln σ / ( ln ε ˙ ) 2 3 ln σ / ( ln ε ˙ ) 3
Experimental data(a)(b)(c)(d)
Materials 16 02928 i045Materials 16 02928 i046Materials 16 02928 i047Materials 16 02928 i048
Arrhenius model(e)(f)(g)(h)
Materials 16 02928 i049Materials 16 02928 i050Materials 16 02928 i051Materials 16 02928 i052
HSModel(i)(j)(k)(l)
Materials 16 02928 i053Materials 16 02928 i054Materials 16 02928 i055Materials 16 02928 i056
Linear model(m)(n)(o)(p)
Materials 16 02928 i057Materials 16 02928 i058Materials 16 02928 i059Materials 16 02928 i060
Quadratic model(q)(n)(r)(s)
Materials 16 02928 i061Materials 16 02928 i062Materials 16 02928 i063Materials 16 02928 i064
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, J.; Xiao, G.; Deng, G.; Zhang, Y.; Zhou, J. The Quadratic Constitutive Model Based on Partial Derivative and Taylor Series of Ti6242s Alloy and Predictability Analysis. Materials 2023, 16, 2928. https://doi.org/10.3390/ma16072928

AMA Style

Zhang J, Xiao G, Deng G, Zhang Y, Zhou J. The Quadratic Constitutive Model Based on Partial Derivative and Taylor Series of Ti6242s Alloy and Predictability Analysis. Materials. 2023; 16(7):2928. https://doi.org/10.3390/ma16072928

Chicago/Turabian Style

Zhang, Jiansheng, Guiqian Xiao, Guoyong Deng, Yancheng Zhang, and Jie Zhou. 2023. "The Quadratic Constitutive Model Based on Partial Derivative and Taylor Series of Ti6242s Alloy and Predictability Analysis" Materials 16, no. 7: 2928. https://doi.org/10.3390/ma16072928

APA Style

Zhang, J., Xiao, G., Deng, G., Zhang, Y., & Zhou, J. (2023). The Quadratic Constitutive Model Based on Partial Derivative and Taylor Series of Ti6242s Alloy and Predictability Analysis. Materials, 16(7), 2928. https://doi.org/10.3390/ma16072928

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop