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
. 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 (
), medium (
) and high (
) 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).
where
is the order of the partial derivative,
and
are indices of the strain rate and the temperature in the stress matrix, respectively (
Table 3). For example, when
= 1 and
= 1,
represents the logarithmic stress with a strain rate of 0.001 s
−1 and a temperature of 900 °C. The maximum values of
and
, 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
= 1 and
= 1.
Equation (7) can be replaced by Equation (9) when = 5 and = 5.
Equation (8) can be replaced by Equation (10) when
= 1.
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 (
), medium (
) and high (
) 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 (
) 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 (
) and third (
) 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 (
) 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 (
) partial derivative is increased. (3) The second (
) 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 (
) 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
. According to the binary Taylor expansion formula,
is expanded at
to obtain Equation (11).
where
m is the number of terms,
is a combination operator and
is m-order partial derivative operator. Equation (11) shows that function
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,
and
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
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).
where
,
and
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).
where
is the combined number of temperatures and strain rates when the strain is fixed (
w = 5 × 5 in this study), and
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).
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:
where
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).
where
is the combined number of temperatures and strain rates when the strain is fixed (
w = 5 × 5 in this study), and
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).
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).