Research Article  Open Access
Multiobjective Optimization of Precision Forging Process Parameters Based on Response Surface Method
Abstract
In order to control the precision forging forming quality and improve the service life of die, a multiobjective optimization method for process parameters design was presented by applying Latin hypercube design method and response surface model approach. Meanwhile the deformation homogeneity and material damage of forging parts were proposed for evaluating the forming quality. The forming load of die was proposed for evaluating the service life of die. Then as a case of study, the radial precision forging for a hollow shaft with variable cross section and wall thickness was carried out. The 3D rigidplastic finite element (FE) model of the hollow shaft radial precision forging was established. The multiobjective optimization forecast model was established by adopting finite element results and response surface methodology. Nondominated sorting genetic algorithmII (NSGAII) was adopted to obtain the Paretooptimal solutions. A compromise solution was selected from the Pareto solutions by using the mapping method. In the finite element study on the forming quality of forging parts and the service life of dies by multiobjective optimization process parameters, the feasibility of the multiobjective optimization method presented by this work was verified.
1. Introduction
Precision forging technology is an advanced technology in recent years; it has important significance for improving mechanics properties of forming and energy saving. However, the forming process is an unsteady large deformation problem, so holes with cracks and other defects often appear; these defects will directly influence the properties and service life of forging parts.
In order to obtain good performance forging parts, recent studies have employed the forging parts optimum by universally adopting two approaches: sensitivity analysis and response surface methodology. In the sensitivity analysis method, Zhao et al. [1, 2] proposed the relationship between preform tool shape and the difference between actual final forging shape and the desired final forging shape objective. And the multiobjective optimum preform tool shape is achieved in the case of axisymmetric disk upsetting problem. Zhao et al. [3] considered the deformation uniformity and deforming force of the final forging, and then the preform die shape of an Hshaped forging process is designed. Acharjee and Zabaras [4] presented that the continuum sensitivity method is extended to three dimensions to accurately compute sensitivity fields. And the cylindrical billet upsetting preform optimization design was carried out. Guan et al. [5] presented the microstructure optimal design at the aim of obtaining the netshape parts and the uniformity of the grain distribution in the final parts. Finally, the forming process of a typical Hshaped axisymmetric disk is illustrated as an example and it shows good results. Though the sensitivity analysis optimization theory could solve problems of the forging part quality, but the sensitivity derivation is very difficult for complex metal plastic forming process especially for threedimensional plastic forming process.
In the response surface optimization method, Yang et al. [6] proposed an approach that combines the finite element method and the response surface method to optimize the preform shapes in order to improve the deformation homogeneity in aerospace forgings. With a typical aero engine disk as an example, the proposed method is verified by achieving an optimal combination of design variables. Tang et al. [7] employed the neural network approximation model for response surface construction in order to overcome the limitation of quadratic polynomial model in solving nonlinear problems. Optimum was achieved by using pattern search optimization method to search response surface describing relationship between preform shape and die cavity fill ratio. Hu et al. [8] proposed adaptive response surface method with particle swarm optimization intelligent sampling for optimization of initial blank shape and blank hold force in sheet forming process. The results show that developed method is able to produce remarkable metamodel for highly nonlinear problems with multiparameter. Oudjene et al. [9] presented a response surface methodology, based on moving leastsquare approximation and adaptive moving region of interest for shape optimization problem. The geometries of both the punch and the die are optimized to improve the joints resistance to tensile loading. There are deficiencies about these studies, mainly through optimization of the shape of die to achieve the purpose of optimization forming and less flash. But in fact the quality of forged not only depends on the level of die design but also depends on the forging process parameters.
The present study is aimed at developing an approach to optimizing the process parameters with the combination of RSM and FEM to improve the precision forging parts quality and the service life of die. And the multiobjective optimization of precision forging process parameters is researched. The multiobjective optimization method based on leastsquare method and Latin hypercube design is proposed. Furthermore, the details of optimizing process parameters are presented with a variable cross section and thickness drive shaft radial precision forging as an example.
2. Theory of Multiobjective Optimization
Multiobjective optimization problem also is called multicriteria optimization problem, which is composed of a set of objective and associated constraint functions. Its mathematical model is shown in the following [10]:where , , and , respectively, are the objective function and inequality and equality constrained function. , , and are number of corresponding functions. is the decision vector.
In most cases, each objective function probably conflicts, which does not make the multiobjective optimization problems produce one global optimum solution. However, such solutions exist: one or several objective functions impossible further optimize while the other objective function not deterioration, such a solution called Pareto optimal [11].
Definition 1. If is nonPareto optimal which is a point in the search space, there is no (in the search space feasibility domain) which makes , .
Definition 2. The set which is made up of all nonParetooptimal solutions is called Paretooptimal set of multiobjective optimization problems.
Kalyanmoy Deb’s nondominated sorting genetic algorithmII (NSGAII) is an improved multiobjective evolutionary algorithm based on NSGA; it can keep the diversity without specifying any additional parameters. It can obtain the true Paretooptimal solutions by the elitepreserving compared to NSGA. And it is widely used to obtain the set of Paretooptimal solutions to handle constrained multiobjective optimization problems [12]. So this work adopts the NSGAII to solve the multiobjective optimization problem.
3. Latin Hypercube Experiment Design
The experimental design is an important link of building response surface model; the sample points selected should be as much as possible to reflect the entire design space. The prediction precision of response surface model depends on the sample point selection, so the sample point selection is very important.
Latin hypercube sampling technique is a kind of constraint random generating uniform sample points experimental design and sampling method, which is a type of design specifically proposed for simulation experiment and a common design method of research on multiple factor experimental design. It has sample memory function and can avoid repeated extraction sample point which has appeared and make distribution at the boundary of sample points involved in sampling [13]. Because each factor can be uniformly applied at each level, therefore, Latin hypercube sampling becomes an effective sample reduction technology. The experimental points produce according towhere is a random number of . One has , , where is number of variables and is number of experiments. is independent random arrangement of .
4. Multiobjective Optimization Function
4.1. Function of Deformation Homogeneity
The forging deformation more homogeneity, local deformation difference smaller, the tissue more homogeneity. Accordingly, internal tissue stress smaller, the crack opportunities caused by the uneven deformation of tissue smaller, the fatigue life of forging parts has increased [14]. When forging parts with uneven deformation, tensile stress in larger areas vulnerable hole or crack, which results in the decrease of forging parts serving time or internal quality of products completely unqualified.
Equivalent strain is an important index to measure the homogeneity forging deformation. In this work, deformation homogeneity evaluation function of forging parts by variance of all unit equivalent strain of work pieces and whole final forging parts average equivalent strain. Considerwhere is the average equivalent strain of final forging parts , is the maximum equivalent strain of final forging parts, is the actual equivalent strain of final forging parts, is the total element of final forging parts, and is function of deformation homogeneity.
4.2. Function of Material Damage
Material damage value at a point is the integral quantity about strain, which reflects the material fracture tendency. When material damage reaches the critical damage value, the material will break. And less than the critical damage value, reflecting the material fracture trend [15]. The smaller the forging internal damage value, the smaller the chance of cracking in use of process. Thus forging products serving relatively prolonged, the fatigue life is relatively increased.
According to different fracture criterion, there are various definitions of forming damage. This work adopts Cockroft & Latham model to describe the damage characteristics, the mathematical model as shown in (4). Considering minimum internal damage after forging deformation, adopts maximum damage of forging all unit for measuring forging deformation damage index. So the function of forming damage is shown in (5). One haswhere is the equivalent stress, is the maximum stress, when , , when , , is equivalent strain material fracture, is equivalent strain, is value of forming damage, is maximum damage of forging parts, and is function of material damage.
4.3. Function of Forming Load
The forming load directly determines the equipment tonnage. If the forming loads are lower, smaller tonnage equipment can complete the product manufacturing. In addition, smaller forming load can also reduce equipment wear and impact fatigue [16] and so can prolong the service life of the equipment and die, thus a corresponding reduction the production cost. The bigger forming load will increase wear between the forged parts and hammer. For a long time, it will cause the cavity size of hammer inconsistent the forged parts accuracy and size requirement, that the forged parts substandard quality.
The practice found that the real influence of die life is the peak load stroke curve of process of precision forging. Therefore this work adopts peak of load stroke curve as the forming load objective function:where is the maximum load of forging processes and is the forming load objective function.
5. Case Study
5.1. Establishment of Finite Element Model
Numerical simulation of radial forging billet dimensions is 41 mm × 5 mm × 580 mm, its material is alloy steel 34MnB5, the hammer and collet materials are Cr12MoVA, and the geometric model of hammer is shown in Figure 1. Figure 2 describes the FE model of billet forming by radial forging, the hammer and billet use fournode tetrahedron element, and collet uses eightnode hexahedron element. Coulomb friction model adopts hammer and billet, and the friction force , where is the friction coefficient and is yield stress of billet material.
5.2. Establishment of Approximation Model
In this work, the DEFORM3D command flow mode called by using a special program, realize the numerical simulation in the process of radial precision forging process parameters are modified and generation of the database file. The Latin hypercube experiment design is shown in Table 1. Since the variation range of process parameters is extremely different, data are encoded. The midpoint of the design parameter is 0, transformed into interval. The scope of process parameters distribution is as follows. Collet axial speed is mm/min. Billet rotating speed is r/min. Hammer radial feed rate is [1.4, 2.2] mm. Friction coefficient between hammer and billet is . Deformation homogeneity function is . Material damage function is . Forming load function is .

Response surface methodology (RSM) can be used to build the approximation model and multiobjective optimization. The secondorder response surface model, which has a good balance between accuracy and cost, is applied in this work. The final quadratic models of response equations are presented in the following:
The analysis of variance (ANOVA) statistical analysis techniques are used to check the fitness of response surface approximation model. The major statistical parameters used for evaluating model fitness are and adjusting (). Generally speaking, the larger the values of and , the better the fitness. The accuracy evaluation indexes of response surface model are shown in Table 2. The data shows that the complex correlation coefficient and the adjusting complex correlation coefficient are close to 1.

The experiment was carried out by comparing forging loads between the actual and response surface model results. The experiment was conducted at GFM SKK10 cold forging machine, and the forging loads were measured in the radial forging processing. Table 3 lists the forging load from the response surface model results comparison with the experimental results. It can be found that the results of response surface model are in agreement with the experimental results. That means the response surface prediction model has high precision and can accurately reflect the mathematical relation between objective function and design parameter and can replace the actual physical model for multiobjective optimization calculation.

5.3. Multiobjective Optimization and Analysis
In this study, the NSGAII was adopted in the optimization process and the optimization parameters are set as follows. The population size is 200, the maximum number of generations is 200, the mutation probability is 0.3, and crossover probability is 0.7.
Figure 3 shows all Paretooptimal solutions obtained by NSGAII in the multiobjective optimization. It can be concluded that NSGAII could provide a variety of alternative design solutions. Among the Paretooptimal solutions, it is clearly found that none of the solutions is absolutely better than any other: any one of them is an acceptable solution. Figures 4–6 show the Paretooptimal solutions vary with deformation homogeneity, material damage, and forming load. The results show that there are collisions between objective and , and . Reducing will usually increase and vice versa. Reducing will increase and vice versa. There are no collisions between objective and . Reducing will usually decrease and vice versa.
In the current work, after the Paretooptimal solution sets are obtained in the multiobjective optimization, the satisfaction function is defined as (8) by the mapping method to find the compromise solution. The smaller is, the better the solution will be. The option 91 in Figure 3 is selected as satisfactory solution based on the satisfaction function defined in (8), shown in Table 4. It shows the initial and optimal value of each design parameter, the deformation homogeneity function increases by 26.92% to 0.297, and the forming load function increases by 3.83% to , while the material damage function decreases by 12.76% to 0.294 compared with the original design parameters:where are the maximum and minimum values of among the Paretooptimal sets, are the maximum and minimum values of among the Paretooptimal sets, and are the maximum and minimum values of among the Paretooptimal sets.

Figure 7 is the equivalent strain contours of forging parts (forging by optimal process parameters). It shows that the deformation homogeneity function increases by 23.93% to 0.290, and the error of approximation model and finite element model is 2.41%. Figure 8 is the material damage contours of forging parts (forging by optimal process parameters). It shows that the maximum damage of forging parts decreases by 15.43% to 0.285, and the error of approximation model and finite element model is 3.16%. Figure 9 is the forming load of die (forging by optimal process parameters). It shows that the maximum damage of forging parts increases by 0.77% to , and the error of approximation model and finite element model is 3.04%. That means the response surface approximation model has high prediction precision.
6. Conclusions
In order to obtain good performance precision forging parts, a multiobjective optimization method for precision forging process parameters is proposed based on the combination of numerical optimization techniques and finite element simulation technology.(1)Based on numerical simulation, a method of character precision forging forming parts quality and improving the service life of die is proposed. The proposed method can cut down the number of experiments by using the Latin hypercube experimental method to build response surface approximation model.(2)Mathematical model of the multiobjective optimization was established based on the secondorder response surface method, and NSGAII was used to obtain the Paretooptimal solutions. It is clearly found that there are collisions between deformation homogeneity, material damage, and forming load. The satisfaction function is defined by the mapping method, which is used to obtain the optimal solution.(3)A variable wall thickness shaft radial precision forging was conducted as case study. In order to solve the conflict between the evaluation function of forming quality and the service life of die, the multiobjective optimization was carried out. After the optimization, the deformation homogeneity and forming load were increased, while the material damage was reduced compared to original process parameters. The feasibility of the proposed multiobjective optimization method was verified by numerical simulation.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This project was supported by the Graduate Innovation Fund Project of Wuhan University of Science and Technology (JCX0018) and the Foundation for Sci & Tech Commission Research Project of Chongqing (CSTC2014YYKFB0187).
References
 G. Q. Zhao, E. Wright, and R. V. Grandhi, “Preform die shape design in metal forming using an optimization method,” International Journal for Numerical Methods in Engineering, vol. 40, no. 7, pp. 1213–1230, 1997. View at: Publisher Site  Google Scholar
 G. Q. Zhao, X. W. Ma, X. H. Zhao, and R. V. Grandhi, “Studies on optimization of metal forming processes using sensitivity analysis methods,” Journal of Materials Processing Technology, vol. 147, no. 2, pp. 217–228, 2004. View at: Publisher Site  Google Scholar
 X. H. Zhao, J. F. Li, X. H. Huang, G. Q. Zhao, and G. C. Wang, “Optimal preform die shape design through controlling deformation uniformity and deforming force in metal forging,” Journal of Mechanical Engineering, vol. 45, no. 5, pp. 193–197, 2009. View at: Publisher Site  Google Scholar
 S. Acharjee and N. Zabaras, “The continuum sensitivity method for the computational design of threedimensional deformation processes,” Computer Methods in Applied Mechanics and Engineering, vol. 195, no. 4849, pp. 6822–6842, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 J. Guan, G.C. Wang, and G.Q. Zhao, “The microstructure optimal design using sensitivity analysis methods in forging process,” Journal of Plasticity Engineering, vol. 14, no. 6, pp. 6–10, 2007. View at: Google Scholar
 Y. H. Yang, D. Liu, Z. Y. He, and Z. J. Luo, “Optimization of preform shapes by RSM and FEM to improve deformation homogeneity in aerospace forgings,” Chinese Journal of Aeronautics, vol. 23, no. 2, pp. 260–267, 2010. View at: Google Scholar
 Y.C. Tang, X.H. Zhou, and J. Chen, “Preform tool shape optimization and redesign based on neural network response surface methodology,” Finite Elements in Analysis and Design, vol. 44, no. 8, pp. 462–471, 2008. View at: Publisher Site  Google Scholar
 W. Hu, L. G. Yao, and Z. Z. Hua, “Optimization of sheet metal forming processes by adaptive response surface based on intelligent sampling method,” Journal of Materials Processing Technology, vol. 197, no. 1–3, pp. 77–88, 2008. View at: Publisher Site  Google Scholar
 M. Oudjene, L. BenAyed, A. Delamézière, and J.L. Batoz, “Shape optimization of clinching tools using the response surface methodology with Moving LeastSquare approximation,” Journal of Materials Processing Technology, vol. 209, no. 1, pp. 289–296, 2009. View at: Publisher Site  Google Scholar
 A. Sadollah, H. Eskandar, and J. H. Kim, “Water cycle algorithm for solving constrained multiobjective optimization problems,” Applied Soft Computing, vol. 27, pp. 279–298, 2015. View at: Publisher Site  Google Scholar
 C. Bazgan, F. Jamain, and D. Vanderpooten, “Approximate Pareto sets of minimal size for multiobjective optimization problems,” Operations Research Letters, vol. 43, no. 1, pp. 1–6, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 S. Chatterjee, K. Abhishek, S. S. Mahapatra, S. Datta, and R. K. Yadav, “NSGAII approach of optimization to study the effects of drilling parameters in AISI304 stainless steel,” Procedia Engineering, vol. 97, pp. 78–84, 2014. View at: Publisher Site  Google Scholar
 M. Ai, K. Li, S. Liu, and D. K. J. Lin, “Balanced incomplete Latin square designs,” Journal of Statistical Planning and Inference, vol. 143, no. 9, pp. 1575–1582, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 X. H. Zhao, G. Q. Zhao, G. C. Wang, and T. H. Wang, “Preform die shape design for uniformity of deformation in forging based on preform sensitivity analysis,” Journal of Materials Processing Technology, vol. 128, no. 1–3, pp. 25–32, 2002. View at: Publisher Site  Google Scholar
 X.W. Chen, J. Wang, J. Chen, and Y. Xue, “Technological parameter optimization of gear billet hot forging process with damage minimization,” Journal of Shanghai Jiaotong University, vol. 39, no. 7, pp. 1070–1072, 2005. View at: Google Scholar
 T. Kroiß, U. Engel, and M. Merklein, “Comprehensive approach for process modeling and optimization in cold forging considering interactions between process, tool and press,” Journal of Materials Processing Technology, vol. 213, no. 7, pp. 1118–1127, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Fayuan Zhu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.