Multi-Objective Optimization of VBHF in Deep Drawing Based on the Improved QO-Jaya Algorithm

Blank holder force (BHF) is a crucial parameter in deep drawing, having close relation with the forming quality of sheet metal. However, there are different BHFs maintaining the best forming effect in different stages of deep drawing. The variable blank holder force (VBHF) varying with the drawing stage can overcome this problem at an extent. The optimization of VBHF is to determine the optimal BHF in every deep drawing stage. In this paper, a new heuristic optimization algorithm named Jaya is introduced to solve the optimization efficiently. An improved “Quasi-opposi-tional” strategy is added to Jaya algorithm for improving population diversity. Meanwhile, an innovated stop criterion is added for better convergence. Firstly, the quality evaluation criteria for wrinkling and tearing are built. Secondly, the Kriging models are developed to approximate and quantify the relation between VBHF and forming defects under random sampling. Finally, the optimization models are established and solved by the improved QO-Jaya algorithm. A VBHF optimization example of component with complicated shape and thin wall is studied to prove the effectiveness of the improved Jaya algorithm. The optimization results are compared with that obtained by other algorithms based on the TOPSIS method.


Introduction
Deep drawing is a forming process that apply complex external force to the sheet through the punch, making the sheet material flow to the ideal direction.Refs.[1] and [2] review many parameters affecting the forming quality in the deep drawing process.For example, Zheng et al. [3] and Su et al. [4] studied the formability and performance of plates with different material properties.Modi et al. [5] and Srirat et al. [6] focused on the BHF to improve the forming quality.Lela et al. [7] determined the mathematical model of the friction coefficient in deep drawing.Karupannasamy et al. [8] and Gong et al. [9] explored influences of lubrication conditions on deep drawing process.Lin et al. [10] established the model related the tools shape parameters to the drawing stressstrain.Among these parameters, the BHF that the binder applies to the sheet is one of the most vital parameters.The friction resistance generated by the BHF can increase the tensile stress in the sheet, controlling the material flow.However, the BHF is a double-edged sword.An inappropriate BHF can result to forming defects.For example, a high BHF may cause sheet tearing, while a low BHF may cause sheet wrinkling.It is necessary to find a method that balance the BHF to maintain the best forming quality.
Compared with the BHF, the VBHF is an effective method that can significantly improve the forming quality.VBHF means the BHF varying differently in different deep drawing steps.The VBHF can reduce or even eliminate the forming defects such as wrinkling, tearing and springback.The division of deep drawing steps are generally based on the shape of forming parts and material fluidity.Kitayama et al. [11][12][13] studied how to determine the VBHF trajectory.The VBHF trajectory has reference value for the experiment of determining the optimal VBHF value at an extent.However, the VBHF trajectories are not the same for deep drawing parts with different shapes.It is difficult to determine the appropriate VBHF only through continuous experiments.A high-efficiency optimization method based on advanced models and algorithms need be adopted.
The traditional optimization of the sheet deep drawing process is developed by the software simulation.Ablat et al. [14] took a review on numerical simulation of sheet metal forming.Kim and Hong [15] adopted finite element (FE) methods and Singh et al. [16] used LS-DYNA to optimize design of deep drawing process.Modi et al. [5] developed programmable logic controller and data acquisition system to determine the variable BHF path, predicting formability by FE simulations.A single simulation always costs a lot of time, even though the modern workstation.Meanwhile, after the samples collected, the establishment of optimization model is also a vital part.The methods based on mechanics theories have complex formula, various parameters and strict prerequisites.It had better be used as a mechanism analysis of forming defects rather than a large-scale solution.To improve the solution effectivity, the surrogate models are proposed to establish the mathematical expression of optimization.Manoochehri and Kolahan [17] and Tian et al. [18] adopted the artificial neural network (ANN) to build relation between process parameters and forming defects, Feng et al. [19] and Xie et al. [20] applied the Kriging model to replace time-consuming optimization solving by simulation test.Kitayama et al. [21][22][23] conducted sequential approximate optimization of process parameters with the radial basis function (RBF) model or the support vector machine (SVM) model.These models can approximate high-precision expression by fitting or interpolation with few samples.On the one hand, the surrogate models improve the optimization speed.On the other hand, the solution of surrogate models also needs high efficiency.
With the development of high-effective algorithms, more and more researchers combine the optimization algorithms with the mathematical simulation software and apply it in the solution of sheet deep drawing.The optimization algorithms mainly divide into deterministic and heuristic algorithms.Lin et al. [24] research that the deterministic algorithms (e.g., linear programming, dynamic programming and integer programming) obtain a global or an approximately global optimum by analyzing properties of the problem.Actually, the deterministic algorithms are not practical in engineering for it too complex to solve non-convex or large-scale optimization problems.Heuristic algorithms use the empirical rules to choose effective methods instead of seeking the answers systematically and in definite steps.It is more flexible, efficient and universal than deterministic algorithms, especially when the optimization problem is expressed in surrogate models.Choosing the heuristic algorithms as solver, the forming defects are set as optimization objectives and the VBHF as design variables.Algorithms solve objectives to find optimal design variables and the simulation software verifies whether the VBHF is excellent out of the simulation results.For instance, Manoochehri and Kolahan [17] employed simulated annealing algorithm to optimize the process parameters.Tian et al. [18] used ant colony algorithm, and Li and Wang [25] used differential evolutionary algorithm to solve VBHF.As the most popular algorithm in solving multi-objective problems, NSGA-II algorithm is applied in the VBHF optimization by Feng et al. [26].
According to summary of Rao [27], the population based heuristic algorithms include evolutionary algorithms (EA) and swarm intelligence (SI) algorithms.Almost all heuristic algorithms own common control parameters such as population size, iterations.Furthermore, different algorithms need their own specific control parameters.For instance, the specific control parameters of genetic algorithm are mutation probability, crossover probability, etc.; simulated annealing algorithm uses annealing start and end temperature, cooling speed, etc.These specific control parameters relate to the algorithm's performance significantly.The inappropriate adjustment of algorithm-specific parameters has the risk of increasing the solution time or falling into the local optimum.To avoid the risk brought by algorithm-specific parameters, Rao et al. [27,28] introduced the teachinglearning-based optimization (TLBO) algorithm and the Jaya algorithm without any algorithm-specific parameters.They only need the common controlling parameters such as population size, generation number, random number.Compared with the TLBO algorithm's two-stage update, the Jaya algorithm just has one update stage.It solves optimization problems more simply with the similar precision.Therefore, we introduce the Jaya algorithm and improve it more suitable as the VBHF optimization solver in deep drawing.
After the comprehensive consideration, a new VBHF optimization method based on the improved Jaya algorithm is proposed.First, we focus on the VBHF optimization of complex shape thin-walled parts.Tearing and wrinkling are the common defects of this type of parts.Thus, we choose tearing and wrinkling as multiobjectives.The evaluation criterions of wrinkling and tearing are established as constrains for limiting defects overlarge.Second, the BHF sample points are selected by Latin Hypercube Sampling (LHS).The Kriging model have advantages of good adaptability, smoothness and minimum estimation variance.We use Kriging surrogate model to approximate the mathematical relation between the VBHF and two defects.Third, we adopt the Jaya algorithm as optimization solver.The Jaya algorithm is added quasi-oppositional strategy and innovated stop criterion for enhancing the convergence speed.Finally, the method we proposed is applied in a case study of complex shape thin-walled parts.The Pareto solution are scored and sorted by decision algorithm TOPSIS.The improved Jaya algorithm are compared with the original Jaya algorithm and NSGA-II algorithm, which prove the efficiency of the QO-Jaya algorithm.The main contributions of this paper are as follows.
1) The optimization model considering wrinkling and tearing is established by LHS and Kriging surrogate method.The operation is solved by a novel algorithm named Jaya that is unique in without algorithm-specific parameters and only having common control parameters, thus superior in saving solution times and universality of the algorithm.2) "Quasi-oppositional" strategy is added to Jaya algorithm for improving population diversity.The strategy of opposition based learning is used to generate a population opposite to the current population.It further diversifies the population and accelerates the convergence rate of Jaya algorithm.3) To ensure the optimization process converges reasonably, a stop criterion based on the "spread" measurement for each iteration is proposed to evaluate the Jaya algorithm optimization situation.The measurement index "spread" is defined to reflect the convergence of the solutions objectively.
The rest of the paper is organized as follow: In Section 2, the evaluation criteria of tearing and wrinkling are put forward and the optimization problem is described preliminary.In Section 3, the Kriging surrogate model provide optimization problem the mathematical expression, including objectives and constrains.Section 4 introduces the improved QO-Jaya algorithm and the improvement scheme.Section 5 is a case study about VBHF optimization on a complex shape thin-wall part.The conclusion and future research are presented in Section 6.

Defect Evaluation Criteria and VHBF Optimization Problem Framework
The main defects in deep drawing process are the wrinkling of the flange, and the tearing at the tangent of the bottom fillet and the straight wall.Wrinkling can be solved by applying BHF.However, with the increase of BHF, tearing at the tangent of the bottom fillet and straight wall are likely to occur.Generally, the wrinkling and tearing are considered at the same time in practical applications.

Tearing Assessment
The tearing of sheet is related to material strength limit.
The tensile stress exceeds the strength limit of the material during the deep drawing.The dangerous section of the drawing part is where the lower end of the cylinder wall meets the outer corner.The excessive strain at the dangerous section causes the wall thickness excessively thin, which leads to tearing.The degree of tearing can be reflected by the maximum thinning rate.Shown in the study by Feng et al. [19,26], the maximum thinning rate is defined as follows: where η max is the maximum thinning rate, t 0 is the initial sheet thickness, and t min is the thinnest sheet thickness of the thin-walled member after deep drawing.In order to prevent tearing during the deep drawing process, it is usually stipulated that the maximum thinning rate shall not exceed the critical value.Therefore, the evaluation criterion of tearing is where η cr is the critical thinning rate.It usually set as 25% according to experience.The maximum thinning rate has been widely used in actual production for its intuition and easy understanding.

Wrinkling Assessment
Wrinkling is a major harmful phenomenon in the deep drawing process.Slight wrinkling affects the forming accuracy and surface smoothness of parts.Severe wrinkling prevents the sheet material from flowing into the gap between the punch and die, resulting in the sheet tearing and become waste.Wrinkling is a phenomenon of plastic deformation instability.There are three reasons for wrinkling: the pressure bar instability caused by the pressure stress, the uneven stress in the deformation zone, and the shear stress effect.Flange wrinkling and sidewall wrinkling are two forms of wrinkling.Flange wrinkling usually reaches the maximum in the (1) initial stage of deep drawing.Sidewall wrinkling is usually accompanied by flange wrinkles.Thus, the flange part is mainly considered in the calculation of the critical wrinkling BHF.
The evaluation indicators of wrinkling are the maximum thickness strain and the wrinkling height.The maximum thickness strain can reflect the wrinkling trend of the sheet during the deep drawing process, and the wrinkling height reflects the wrinkling severity of the drawn part truly during deep drawing.
When the sheet is wrinkled, the distance L between the die and the blank is greater than the initial thickness t 0 .Therefore, the wrinkle evaluation criterion can be expressed by the distance between the die and the blank where L is the distance between the die and the blank, t 0 is the initial thickness of thin-walled components, B is a safety factor.According to engineering experience, B always takes 1.2.
Springback is another common forming defect in addition to wrinkling and tearing.In this paper, we omitted the consideration of springback because it seldom appears in the complex shape thin-wall part.

Model of the VHBF Optimization Problem
The optimization of VBHF includes the determination of design variables, objectives and constraints.The design variables are BHFs in different deep drawing stage, that is, in the different punch stroke.In Figure 1, the total punch stroke is partitioned into n steps S = [S 1 , S 2 , . . ., S n ] .Correspondingly, the BHFs takes n different values . Every F VBH is a design variable in the optimization.The divided steps are not changeless.It is mainly determined by the characteristics of the hydraulic press itself and the drawing parts shape.
After the VBHF curve established, the objectives and constrains can be determined.The objectives are that minimize the evaluation of tearing and wrinkling.There are three constrains.One is the range of BHF, and another are the evaluation of wrinkling and tearing that not exceed the limitation settled in Section 2.1 and Section 2.2.The description of optimization problem is as follows: (3) After the model of optimization problem built initially, the specific expression of problem is needed for solving.There are two main methods to build a solution model.One is the method based on the software simulation, that is, use finite element analysis model to select BHF and experiment continuously.The other one is using the surrogate model to approximate the evaluation of tearing and wrinkling.The simulation is time-consuming and requires high computer hardware.However, it is easier to use surrogate model and the accuracy of the surrogate model is guaranteed.In this paper, the Kriging surrogate model is selected for modeling.The simulation software, Dynaform, is applied to verify the effectiveness of the optimums.

Mathematical Expression Based on Kriging
Model for VBHF Optimization in Deep Drawing

Introduction of Kriging Surrogate Model
The Kriging surrogate model has the characteristics of smoothness and minimum estimation variance.It can fit problems with high nonlinearity ideally.It is widely used to fit low-order or high-order nonlinear problems, concluded by study of Zhao et al. [29].The Kriging surrogate model is essentially an interpolation method (4) The step partition of the punch stroke and the VBHF based on statistical theory.Its basic principle is to construct an approximate model of the objectives through known samples and their responses.And then use the constructed model to predict unknown samples and their responses.The Kriging surrogate model consists of a global model and a local random deviation function, that is, a combination of a parametric model and a non-parametric random process.The Kriging surrogate model is more flexible than a single parametric model.And it also overcomes the limitation of non-parametric models in processing high-dimensional data.
The expression form of Kriging surrogate model is where T is a known regres- sion model, generally a polynomial model, including polynomials of zero-order, first-order, second-order and so on; β = [β 1 , β 2 , ..., β q ] T is the corresponding regression coefficients to be estimated; z(x) is a stationary stochastic process, and ) is a correlation function about parameter θ, which generally takes Gaussian correlation function Fβ provides a global approximation of the Kriging model and z(x) provides an approximation of the local variation of the model.Kriging surrogate model is used to realize the nonlinear fitting between the input data (VBHF) and the output data (the response value of each VBHF).In Section 3.2, the VBHF optimization model constructed by Kriging method will be introduced in detail.

Kriging Surrogate Model of VBHF Optimization
In this section, the optimal Latin Hypercube Design (LHD) will be used to select sample points.The Kriging model will be constructed based on these sample points.Then the relationship will be established between VBHF and forming defects in the deep drawing process.The whole steps of establishing and solving the VBHF optimization problem is shown in Figure 2.
McKey et al. [30] proposed the Latin Hypercube Design (LHD).The LHD have been widely used in sampling in various large-scale design spaces.In the VBHF optimization, the sample array where k is the number of samples, V i (i = 1,2,...,k) is the ith sample point in VBHF design space, and Y i is the response corresponding to the sample V i obtained by the Dynaform simulation software. (5) In view of the fact that the sample points taken in the symmetric LHD are more uniform, it is also more conducive to the establishment of subsequent approximation models.In this paper, the symmetric LHD is used for the adoption of the design space.Considering the 6-dimensional design space is hard for visualization, we take a 2-dimensional design space to extract 30 sample points as an example to illustrate the symmetric LHD sampling scheme, as shown in Figure 3.
Based on the VBHF optimization model represented by Eq. ( 4), we convert the objectives and constrains with solvable mathematical form.The approximate models of wrinkling and tearing are constructed based on Kriging surrogate model, as follows: where fW (F VBH ) and fT (F VBH ) are the approxi- mate expresses of wrinkling and tearing respectively, β W F W (F VBH ) and β T F T (F VBH ) are the regression mod- els of wrinkling and tearing respectively, z W (F VBH ) and z T (F VBH ) are the random deviation in the regression models of wrinkling and tearing respectively.

Multi-objective Quasi-oppositional Based Jaya Algorithm
Jaya algorithm is a relatively new swarm intelligence algorithm, which was proposed first by Rao [27] in 2016.The optimization principle of Jaya algorithm is that the solution should move towards the best solution but avoid the worst one.The most significant feature of Jaya algorithm is which only has common control parameters without algorithm-specific control parameters.This feature makes the algorithm more easily in initial parameters setup, thereby reducing the time cost and ensuring the stable performance.Rao et al. [31] and Zhang et al. [32] proved the effectiveness of algorithm through solving various constrained and unconstrained engineering optimization problems or benchmark problems.And the computational results reveal that the Jaya algorithm ( 7) is superior to or comparable with other optimization algorithms.
After the Jaya algorithm presenting, there have been various improvement and strategies added in it for better performance.In order to handle multiple objectives simultaneously, Rao et al. [33] proposed the multi-objective Jaya (MO-Jaya) algorithm.Some changes are aimed to further diversify the population and increase the convergence speed of the Jaya algorithm.For example, Rao and More [34] proposed the self-adaptive Jaya algorithm that determined the population size automatically to improve the population diversity; Rao and Rai [35] advanced the quasi-oppositional based Jaya algorithm, which generated a population opposite to the current population to keep the randomness of Jaya algorithm.Meanwhile, it is also important for algorithm to balance the global exploration and local exploitation.In researches of Farah and Belazi [36] and Yu et al. [37], the chaotic Jaya algorithm use the chaotic sequence generated by the chaotic graph rather than random numbers to balance the exploration and exploitation.Rao and Saroj [38] added the elitist stratagem in Jaya algorithm.It replaces the worst solution(s) with the elitist one(s) to effectively control Jaya algorithm for the transition from global exploration to local exploitation.Besides, there are some other tiny improvements, such as adopting a neighborhood search strategy to enhance the population diversity, linearly decreasing inertia weights to enhance development efficiency, or a combination of the above strategies.The summary of the improvement strategies on Jaya algorithm is shown in Figure 4, which is divided according to the specific improvement objects.
Nowadays, the Jaya algorithm is applied in continuous/discrete, high-dimensional/low-dimensional, or linearity/nonlinearity optimization problems.The professional fields of application include electronic and electrical engineering, energy engineering, mechanical engineering, management science, computer science and so on.It has not been applied to drawing forming process optimization yet.In this paper, we will first use

Introduction of the Multi-objectives Jaya Algorithm
Like the common heuristic algorithm, Jaya algorithm has the common control parameters such as the population size, generation number, etc.The initial population is generated randomly within the ranges of variables.However, it is different from other heuristic algorithms that Jaya algorithm only have one update function.It is more quickly in updating population.In every generation, the population updates based on Eq. ( 8): where G is the max generation number; g means the gth generation.Pop is the population size; p means the pth individual.n is the number of variables; j means the jth variable.In Eq. ( 7), W ′ g,p,j is the updated new individual, W g,p,j is the old individual, r 1,g,j and r 2,g,j are the random number in [0,1].W g,best,j and W g,worst,j are respectively the best and worst solutions in the gth generation, the jth variable.(W g,best,j − |W g,p,j |) means the tendency of the solu- tion towards the best solution, while (W g,worst,j − |W g,p,j |) means the tendency of the solution far from the worst solution.In the function, the random number act as scaling factor to ensure good exploration performance.The best and worst solutions guarantee the positive direction of the update solution.At the end of the iteration, all acceptable solutions will be retained and become the inputs for the next generation.
Since the VBHF optimization we researched is a multi-objective problem, the Jaya algorithm needs some strategies added to convert to the multi-objectives Jaya (8) W ′ g,p,j = W g,p,j + r 1,g,j (W g,best,j − |W g,p,j |) − r 2,g,j (W g,worst,j − |W g,p,j |), g = 1, 2, . . ., G; p = 1, 2, . . ., Pop; j = 1, 2, . . ., n, algorithm.In the multi-objectives Jaya algorithm, the concepts of constraint-dominance sorting, non-dominance sorting and crowding distance computation are used to determine the rank of solutions.These principles are the keys for conducting the searches of Pareto frontier.The sorting principles are explained by Figure 5.The flowchart of the multi-objectives Jaya algorithm is shown in Figure 6.
The detailed solving steps of multi-objective Jaya algorithm is described as follows.
Step 1: Pop solutions are generated randomly as an initial population.The initial population is sorted based on the principles of constraint-dominance and non-dominance.
Step 2: First, the constraint-dominance is used to determine the superiority between solutions initially.Then the non-dominance and crowding distance are conducted to further determine the priority of the solution.The solution with a higher rank is better than the other one.If the solutions have the same rank, the solution with a higher crowding distance is considered better than the other one.Step 3: Choose the solution with the highest rank (rank = 1) as the best solution and the solution with the lowest rank is the worst solution.Then the solution of next generation can be updated according to Eq. (7).
Step 4: After the solutions updated, the new solutions combine with the old solution as 2Pop solutions.Reorder these solutions according to the principles of constraintdominance, non-dominance sorting and crowding distance computation.Then select Pop solutions as the new population based on the new sorting.
Step 5: Turn to Step 3 to update generation until satisfy the stop criterion.
The concept of constraint dominance guarantees that feasible solutions have a higher rank than infeasible solutions.Among the feasible solutions, the superior solution (non-dominated solution) ranks higher than the dominated solution.Among the infeasible solutions, a higher rank is assigned to the solution with less overall constraint conflict.The use of crowding distance for nondominated sorting ensures that the solution is selected from the sparse area of the search space.

Introduction of the Multi-objectives Jaya Algorithm 4.2.1 Quasi-oppositional Based Strategy
QO-Jaya algorithm means Quasi-oppositional based Jaya algorithm, which is added a concept of opposition based learning in common Jaya algorithm.A population opposite to the current population is generated to further diversify the population and accelerate the convergence rate of Jaya algorithm.We improve the quasi-oppositional population generation strategy for better practicality in programming.The opposite population is generate based on Eqs. ( 9)-( 11): where W L j and W U j are the lower and upper bounds of the variables range.W g,p,j is the current population and W q g,p,j is the opposite population to the current population.

Improved Stop Criterion
In most optimization algorithm, the stop criterion is usually the generation number reaches a max generation number or the cumulative individuals reach the limitation.Both the stop criteria are set rigidly by cumulative sum of counting parameters without reflecting the solutions condition.In this case, some solving progress that solutions have converged may continue searching meaninglessly.This greatly reduces the solution speed and waste resources on unnecessary computation.
In this paper, a measurement index 'spread' is defined to reflect the solutions convergence objectively.We propose a stop criterion based on spread: when the spread does not change much, the spread of current generation is smaller than the average of previous spread, and this condition is maintained for a certain number of times, the algorithm will stop.The spread is calculated by Eq. ( 12).We set the population size to 50, 100, and 500 for testing, which proves the stop criterion effective.where k is the number of objective functions, σ is the standard deviation of the crowding distance measure of points that are on Pareto front with finite distance.μ is the norm of the difference between the minimum objective on the Pareto front of the current generation and that of the previous generation.As there are two objectives in this study, μ is the sum over the two norms of objective functions shown as Eq. ( 13).W g,best and W g−1,best indicate the parameters for the minimum objective on the Pareto front of the current generation and the previous generation, respectively.
The spread-based stop criterion can make the convergence measurement of the design variables more objective, for it is not affected by parameters other than the number of iterations.The spread not only reflects the change of the congestion degree of the Pareto front, but also reflects the movement of the Pareto front.It is also found from experiment that many Pareto fronts may jump and escape from the best solution obtained before after satisfying the stop criterion once.This indicates that the Pareto front is not stable finally yet.Therefore, in order to maintain more stable and reliable Pareto optimal solutions, the algorithm will finally stop after reaching the stop criterion 5 times.Write the program in the algorithm to count the total number of the times that accumulatively satisfying the stop criterion.The flowchart of the improved QO-Jaya algorithm is shown in Figure 7.

Case Study
The complex outer shell thin-walled parts have many sharp inflection points.It is easy to cause wrinkling, tearing and other problems in the deep drawing process.Figure 8 shows the complex outer shell thin-walled parts that we optimized in case study.In this section, we (12) conduct the VBHF optimization of a specific complex outer thin-walled parts.

Parameters of the Complex Outer Shell Thin-walled Part
The material is 304L stainless steel.The simulation analysis model of the complex outer shell thin-walled part is shown in Figure 9.
The entire drawing stroke is divided into 6 phases, each phase adopts a different BHF.According to experience, the range of each BHF is [500 kN, 2000 kN], taking wrinkling and tearing as objectives.The constrains are all set • Tearing objective function

• Tearing objective function
The kriging model is the best unbiased estimate based on the interpolation method.It completely passes through the sample points.Figure 10 shows the mean square error of Kriging model.
The improved QO-Jaya algorithm is used to solve the established optimization model.Set the max individual number as 37500, and the population size is 100.The 100 initialized individuals.

Iteration Process
After initializing the 100 individuals of a population, generate the quasi-opposition population of the initialized population.Combine the two populations to sort and ( 14) selected the first 100 individuals as a new population.Then update the current population, sort the new and old individuals to the best 100 individuals.Update the quasiopposition population and repeat the iteration as before.
The algorithm will stop when the principle satisfied in 5 times.The principle is that the spread has not changed much and the final spread is less than the average of the most recent spreads.

Optimization Results
The solutions are scored and sorted by the decision algorithm TOPSIS.The solution with highest score is shown in Table 1.It can be seen that the risk of wrinkling has reduced 5.1% and that of tearing has reduced 42.6%.It can be concluded that the VBHF optimization improves the forming quality of complex outer shell thin-walled part in deep drawing process.For better representation of the optimization effect, we use Dynaform to simulate the optimized VBHF, as shown in Figures 13, 14.The Kriging surrogate model is adopted to express the optimization in mathematical form.Then the multiobjectives QO-Jaya algorithm is introduced in solving the optimization problem.A numerical example of the complex outer shell thin-walled part is conducted to verify the effectiveness of the method.We use TOPSIS to select the best solutions that each algorithm generates, and put them in Dynaform for simulation and visual comparison.
We adopt Jaya algorithm as solver for it without algorithm-specific parameters and only having common control parameters, which save solving times and ensure the universality of the algorithm.The improved multi-objective QO-Jaya algorithm in this paper includes three parts of contribution.First, the concepts of constraint-dominance sorting, non-dominance sorting and crowding distance computation are used to determine the rank of solutions.Second, the strategy of opposition-based learning is used to generate a population opposite to the current population.It further diversifies the population and accelerates the convergence rate of Jaya algorithm.Thirdly, the measurement index 'spread' is defined to reflect the solutions convergence objectively.A new stop criterion is raised based on the spread.The improved multi-objective QO-Jaya algorithm (improved QO-Jaya algorithm in short) is compared with the common Jaya algorithm, NSGA-II algorithm and so on.It verifies the good performance of improved QO-Jaya algorithm in solving optimization.
In future, we can study the VBHF coupling with other parameters in process optimization of deep drawing.Consider whether the parameters have uncertainty in optimization.The solving algorithm can be focused on how to balance the exploitation and exploration.We have considered the hybrid algorithm initially to improve the problem.

Figure 2
Figure 2 Framework of VBHF optimization using proposed method in deep drawing

Figure 5
Figure 5The sorting principles of constraint-dominance sorting, non-dominance sorting and crowding distance computation

Figure 6
Figure 6 Flowchart of the multi-objectives Jaya algorithm r 3 is the random number in [0, 1]. a means the midpoint of the variables interval and b means the mirror point of the variables.

Figure 7
Figure 7 Flowchart of the improved QO-Jaya algorithm

Figure 8 Figure 9 Figure 10
Figure 8The complex outer shell thin-walled parts

Figure 11
Figure 11The comparison between spread and the mean of spread

Figure 12
Figure 12The changes of spread

Figure 13 Figure 14 Figure 15
Figure13 The comparison of tearing between complex outer liner thin-walled parts before and after optimization: (a) Before optimization, (b) After optimization

Figure 16 Figure 17
Figure 16 Comparison of solution scores between QO-Jaya algorithm and the NSGA-II algorithm

Table 1
Comparison between the original constant BHF and the optimized VBHF