Abstract
Nonisothermal forging is an efficient plastic forming method for titanium alloys, but at the same time, it can produce large and uneven residual stress, which seriously affects the service life of components. In order to quantitatively analyze the influence of forging process parameters on the residual stress of Ti-6Al-4V alloy forgings, a numerical model was first established and optimized in combination with experiments. Then, the effects of deformation temperature, deformation degree, and deformation speed on the residual stress of forgings were analyzed by orthogonal test, and the optimal combination of forging process parameters was obtained. Finally, the multiple regression analysis was employed to propose multivariate regression models for the prediction of the average equivalent residual stress. Results show that the prediction model can be used for predicting the residual stress of Ti-6Al-4V alloy forgings with a higher reliability.
1. Introduction
Titanium and its alloys are often selected for use in high-value structural components in a number of applications due to their high specific strength, relatively low density, excellent corrosion resistance, and fatigue resistance [1–3]. Aerospace industry accounts for the majority of titanium consumption in the world, and more than 80% of titanium is used to produce a variety of alloys, among which Ti-6Al-4V is the most widespread variant [4]. In order to ensure the mechanical properties and dimensional stability of related parts, nonisothermal forging is usually used to process them [5]. Nonisothermal forging is a deformation method with high strain rate, which can greatly improve the mechanical properties of titanium alloys. However, residual stress often occurs in the forging process due to thermal conductivity, high deformation resistance, and narrow forging temperature range of Ti-6Al-4V alloy [6, 7]. In the process of machining, residual stress will cause the forging to deform and even crack, which reduces the qualified rate [8, 9]. In the service process, residual stress will affect the mechanical properties of components and reduce their service life [10, 11]. Therefore, it is necessary to reduce the residual stress of Ti-6Al-4V alloy forgings in the forging process.
At present, many researchers have conducted a lot of studies on reducing residual stress. Araghchi et al. [12] studied the effect of different concentrations of quenchant on the residual stress of the alloy. The results showed that quenching in 15% polymer solution and aging at 190°C for 12 h could reduce the residual stress, but the mechanical properties of the alloy decreased. Robinson et al. [13] analyzed the influence of quenching temperature on residual stress, and the results showed that quenching had little effect on residual stress when the quenching temperature was low, but the residual stress of the alloy could be reduced when the water temperature was above 100°C. Ya Bo Dong et al. [14] analyzed the effect of heat treatment on the residual stress of forged plate and improved the balance of mechanical properties and residual stress of forged plate by changing the heat treatment process parameters. It can be seen that most researchers reduce the residual stress of forgings by adjusting the heat treatment process parameters. However, as far as the present technology is concerned, it is impossible to completely eliminate the residual stress of forgings by heat treatment [15]. In addition, it was found that the qualified rate of a certain type of blade is increased by 20%–30% after heat treatment and the adjustment of cutting process parameters, but there are still a lot of scrapped products due to deformation, and the total qualified rate is only 70%. After the analysis of the scrapped products, it was found that forging process parameters have a great influence on the residual stress of forgings. In addition, Karunathilaka et al. [16] analyzed the influence of lubrication and forging load on residual stress and found that there is a positive relationship between the forging load and residual stress. Ameli et al. [17] investigated the effects of process parameters such as workpiece and die geometries, percentage of deformation, and workpiece motions on residual stresses. Ye Zhang et al. [18] analyzed the influences of different process parameters on residual stress during turbine disc hot die forging by numerical simulation. They all found that the key forging process parameters have a considerable influence on the residual stress of forgings [19]. Therefore, it is of great practical significance to reduce the residual stress and improve the distribution of residual stress in forgings by optimizing the forging process parameters.
In this paper, the residual stress of Ti-6Al-4V alloy forgings under different forging processes (deformation temperature, deformation degree, and deformation speed) were analyzed based on the finite element method, and the optimal combination of forging process parameters was obtained through the orthogonal test. Then, a mathematical model for predicting the average equivalent residual stress was established by regression analysis of the simulation data, which can provide effective guidance for the actual forging process.
2. Calculation Model
2.1. Establishment of the Simulation Model
In the process of numerical simulation, the blank is defined as an elastic-plastic body because of the elastoplastic transformation of materials while the dies are defined as a rigid body. In addition, due to the temperature difference between the blank and the dies, as well as between the blank and the air, the simulation control parameters such as the heat conduction between the blank and the dies, the heat convection between the blank and the air, and the heat radiation of the blank should be considered when establishing the simulation model. The heat transfer of forgings is shown in Figure 1 [20].

The geometric models of the upper and lower dies and the forging blank were established using three-dimensional drawing software, as shown in Figure 2(a). The size of the forging blank is Φ60 mm × 100 mm. Considering the lengthwise symmetrical characteristic and the complexity of the elastic-plastic finite element analysis, one fourth of the forging billet was taken for analysis, and the finite element model is shown in Figure 2(b).

(a)

(b)
The accuracy of the simulation results is affected by the control parameters. Therefore, in order to obtain more accurate control parameters, this paper combines physical experiments and numerical simulation to analyze the residual stress of Ti–6Al–4V alloy forgings after nonisothermal forging [21]. Firstly, the forging experiment and residual stress detection test were carried out to obtain the residual stress of the forging. Two sets of different forging process parameters were selected (one set is deformation temperature 975°C, deformation degree 50%, and deformation speed 10 mm/s, and the other set is deformation temperature 925°C, deformation degree 50%, and deformation speed 20 mm/s). Five experimental samples were taken under each set of process parameters, and the average value of residual stress was taken as the final result. The specific experimental process is shown in Figure 3 [20]. Then, the corresponding numerical simulation was carried out. Based on the forging residual stress measured in physical experiment, the error between experimental results and simulation results was analyzed. Finally, the simulation control parameters were adjusted until the residual stress values of the simulation are as close as possible to the experimental measurements. The final simulation control parameters are shown in Table 1.

2.2. Results’ Analysis
In order to intuitively show the distribution of residual stress of Ti-6Al-4V alloy forgings under nonisothermal forging, a set of forging process parameters (deformation temperature 925°C, deformation degree 50%, and deformation speed 100 mm/s) was selected for simulation analysis according to the engineering practice, and the result is shown in Figure 4.

In order to further reveal the distribution characteristics of residual stress, 11 tracking points with a spacing of 5 mm were selected from the upper surface to the lower surface along the central axis of the forging. The distribution of the residual stress at the tracking points along the XYZ direction is shown in Figure 5. It can be seen from Figure 5 that the distribution characteristics of the residual stress in the X and Y directions of the forging are basically similar, which are symmetrical but uneven about the center of the axis. In the direction of the central axis, the residual stress on the surface of the forging is tensile stress while the center is compressive stress. From the surface to the core, the tensile stress is transformed into compressive stress. The residual stress in Z direction is less than that in X and Y directions, and the distribution is more uniform.

3. Optimization of Forging Process Parameters Based on Orthogonal Test
3.1. Orthogonal Test Factors
In nonisothermal forging, it is necessary to select the appropriate forging process parameters for analysis to correctly guide the actual processing. According to the characteristics of Ti-6Al-4V alloy, such as transformation temperature (985°C), forging method, allowable forging deformation degree, and speed, the deformation temperature, deformation degree, and deformation speed were determined as the key process parameters. Through further research, it was found that these parameters have a great influence on the residual stress of forgings [20].
3.2. Orthogonal Test Schemes
Taking deformation temperature, deformation degree, and deformation speed as test factors, five levels are set for each factor, as shown in Table 2.
In accordance with the L25 (53) orthogonal table, 25 sets of orthogonal test schemes are established, as shown in Table 3, where T is deformation temperature, D is deformation degree, S is deformation speed, and σ is the average equivalent residual stress.
3.3. Orthogonal Test Results
Range analysis is usually used to analyze the results of an orthogonal test [22]. It can be described as shown in the following equation [23]:where Ej is the value of a certain factor, Ki is the average of each factor, and R is the range and used to estimate design variable sensitivities [24]. In general, the greater the Ki value, the higher the indicator values under this level, and if the value of R is larger, the factor value would be more influential [25]. The range analysis of the orthogonal test results was calculated according to equations (1) and (2), as shown in Table 4
It can be clearly found from Table 4 that the primary and secondary relationship of factors affecting the residual stress of the forging is deformation temperature, deformation degree, and deformation speed. And, their contribution rates to the residual stress of forgings are 56%, 28%, and 16%, respectively. Comparing the minimum values of K1∼K5 under each factor, the optimal combination of process parameters with the smallest residual stress can be determined as deformation temperature 1025 °C, deformation degree 30%, and deformation speed 1000 mm/s.
3.4. Numerical Simulation of Optimized Process Parameters
The forging process of Ti-6Al-4V alloy was simulated under the optimal process parameters of deformation temperature 1025 °C, deformation degree 30%, and deformation speed 1000 mm/s. The results are shown in Figure 6. It can be seen from Figure 6 that the residual stress of the forgings under the optimal combination of process parameters basically belongs to normal distribution, without obvious stress concentration, and the distribution is relatively uniform. This is because when the deformation temperature (1025 °C) far exceeds the transformation temperature of Ti-6Al-4V alloy (985°C), the forging is completed in the β single-phase region. The plastic forming ability of the forgings is significantly improved and the deformation resistance is greatly reduced, which greatly improves the deformation uniformity. Furthermore, when the deformation degree is 30%, the main softening mechanism of Ti-6Al-4V alloy forgings is dynamic recovery. The residual stress of the forgings decreases sharply during the recovery stage. In addition, the deformation temperature is high and the degree of deformation is small at this time, which reduces the heat converted and weakens the thermal effect of metal deformation, thereby greatly reducing the residual stress of the forging and improving the distribution of the residual stress.

(a)

(b)

(c)

(d)
However, when the forging temperature is above the β-transition temperature, there are a lot of coarse β-grains in the forgings. The coarser the grain, the more uneven the plastic deformation and the greater the internal stress concentration. Moreover, the coarse grains reduce the chance of crisscross between grains, which will be conducive to the propagation and development of cracks, and the strength and toughness will become worse. So, the mechanical properties will be reduced. Therefore, for the parts that require high-dimensional stability but not too much mechanical properties, this combination of forging parameters can be used as an effective guide.
4. Development of Prediction Model
4.1. Regression Model and Calculation
In order to further quantitatively analyze the relationship between the key forging process parameters and the average equivalent residual stress, regression analysis was carried out on the simulation data to find the mapping relationship between them, so as to determine a mathematical model that can effectively predict the average equivalent residual stress of the forgings and provide effective guidance for the actual forging process [26].
Based on Table 2 and the analysis results of Table 4, the change curve of the average equivalent residual stress under each factor is shown in Figure 7.

It can be seen from Figure 7 that when the deformation temperature exceeds the transformation temperature (985°C) of Ti-6Al-4V alloy, the average equivalent residual stress of forgings decreases sharply. However, there is a linear relationship between the deformation temperature and the average equivalent residual stress in the two different temperature ranges above and below the transformation temperature. In addition, when the deformation degree is between 40% and 70%, there is a linear relationship between the deformation degree and the average equivalent residual stress. And, when the deformation velocity is between 100 mm/s and 1000 mm/s, the relationship between the deformation velocity and the average equivalent residual stress is basically linear.
Based on the above characteristics, the deformation temperature was divided into two intervals. Then, regression model I was established by taking deformation temperature, deformation degree, and deformation velocity as independent variables and average equivalent residual stress as dependent variable. At the same time, considering the difference in the order of magnitude of the selected parameters, the logarithmic transformation was carried out to establish regression model II.
Set the average equivalent residual stress as σ (MPA), deformation temperature as t (°C), deformation degree as d (%), and deformation speed as s (mm/s), among which t∈[925, 980]∪[ 990, 1025 ], d∈[40, 70], and s∈[100, 1000].
The two regression models were as follows:
Fifty sets of sample data were obtained using the modified finite element model of residual stress of Ti-6Al-4V alloy forgings under nonisothermal forging. Among them, 25 sets of sample data with deformation temperature between 925°C and 980°C are shown in Table 5, and the other 25 sets of sample data with deformation temperature between 990°C and 1025°C are shown in Table 6.
Based on the sample data in Tables 5 and 6, the mathematical analysis software was used to analyze and solve the two regression models. The results are shown in Figures 8 and 9.

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)
The regression models obtained after calculation are as follows:where t∈[925, 980]∪[990, 1025], d∈[ 40, 70], and s∈[100, 1000].
By comparing Figures 8 and 9, it can be seen that the normalized residual of regression model II is closer to the normal distribution and has higher accuracy. Therefore, regression model II was selected as the prediction model of the average equivalent residual stress.
4.2. Model Validation
Since the simulation model has been modified by experiments, it can be used to verify the accuracy of the prediction model. Another 5 sets of forging process parameters were selected as test samples for simulation and model prediction, as shown in Table 7. The comparison result is shown in Figure 10, from which it can be seen that the simulation data are very close to the predicted data and the error is small. This indicates that the prediction model can accurately reflect the mapping relationship between the forging process parameters and the average equivalent residual stress of Ti-6Al-4V alloy forgings. It also indicates that the prediction model can effectively guide the prediction of the residual stress of Ti-6Al-4V alloy forgings in practical engineering.

5. Conclusions
By using the modified finite element model, the primary and secondary relationship of key factors affecting the average equivalent residual stress of Ti-6Al-4V alloy forgings was analyzed and the optimal process parameters with minimum residual stress were obtained. Then, the multiple linear regression analysis was carried out based on the simulation data, and the mathematical model which can effectively predict the average equivalent residual stress was determined. The significant conclusions have been drawn below:(1)In nonisothermal forging, the deformation temperature is the main factor affecting the residual stress of the forging, followed by the degree of deformation, and the deformation speed has the least influence on the residual stress of the forging.(2)The optimal combination of process parameters for the minimum residual stress of forgings without considering mechanical properties is deformation temperature 1025°C, deformation degree 30%, and deformation speed 1000 mm/s. For the parts that require high-dimensional stability but not too much mechanical properties, this optimal combination of process parameters has certain reference significance.(3)In the regression analysis, according to the distribution characteristics of the residual stress of Ti-6Al-4V alloy forgings in the upper and lower regions of the transformation point (985°C), the regression model was divided into two temperature regions, which can accurately reflect the mapping relationship between the key process parameters and the average equivalent residual stress.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors gratefully acknowledge the financial support of this study by the National Natural Science Foundation of China (Grant no. 51775427).