Study on Parameter Optimization Design of Drum Brake Based on Hybrid Cellular Multiobjective Genetic Algorithm

In consideration of the significant role the brake plays in ensuring the fast and safe running of vehicles, and since the present parameter optimization design models of brake are far from the practical application, this paper proposes a multiobjective optimization model of drum brake, aiming at maximizing the braking efficiency and minimizing the volume and temperature rise of drum brake. As the commonly used optimization algorithms are of some deficiency, we present a differential evolution cellularmultiobjective genetic algorithm DECell by introducing differential evolution strategy into the canonical cellular genetic algorithm for tackling this problem. For DECell, the gained Pareto front could be as close as possible to the exact Pareto front, and also the diversity of nondominated individuals could be better maintained. The experiments on the test functions reveal that DECell is of good performance in solving high-dimension nonlinear multiobjective problems. And the results of optimizing the new brake model indicate that DECell obviously outperforms the compared popular algorithm NSGA-II concerning the number of obtained brake design parameter sets, the speed, and stability for finding them.


Introduction
The safe reliability of the vehicles is attracting more and more attentions with the sharp increase in the vehicles.Since the improvement of the brake design and manufacturing is of great significance in ensuring the brake quality thus the vehicle's safe running, the brake parameter optimization design which is a key part in the whole brake design procedure is becoming increasingly important.Therefore, the research in this field is becoming a hot topic.
In general, in the brake parameters optimizing procedure, the first step is to set up parameter optimization model and the next step is to calculate the model with an optimization algorithm.A lot of researches have been done in this field 1-13 , most of which are done by Chinese scholars.Their researches mainly focus on two aspects: the establishment of the optimization model 1, 3, 6, 7 and the design, selection, and application of the algorithms for solving the model 2, 4, 9, 10 .Although plenty of achievements have been gained in the existing studies, they still need to be improved to some extent, especially in the following three aspects.First, the complexity of the model is not satisfactory.In general, there are a lot of brake performance indicators, such as brake temperature rise, brake drum volume, and braking efficiency.However, in existing models, only one or two of them are selected as optimization objectives, which is far from engineering practice.Second, the multiobjective optimization problems are often solved by being transferred into single-objective problems with weighting or separate calculating method, which ignores the fact that the solution of a multiobjective optimization problem is a set.Third, the performance of the adopted algorithms cannot meet the requirements.The frequently used algorithms in the existing studies are the MATLAB optimization toolbox, basic genetic algorithms, basic ant colony algorithm, and so forth.However, these algorithms are easy to be trapped in local convergence when solving complex problems, for which the gained results are undesirable.
Based on foregoing statement, this paper establishes a multiobjective optimization design model of drum brake, which aims at maximizing the braking efficiency and minimizing the volume and temperature rise of the drum brake.It is expected that this model could meet the requirements of practical engineering design better.Besides, in order to deal with the model more efficiently, a differential evolution cellular multiobjective genetic algorithm DECell is proposed by introducing a differential evolution strategy into the cellular genetic algorithm in which an individual may only interact with individuals belonging to its neighborhood.The proposed algorithm is then applied to solve the brake model.
The rest of the paper is organized as follows.The next section shows the details of the multiobjective optimization model of drum brake.Section 3 presents the differential cellular multiobjective genetic algorithm DECell .Section 4 tests the performance of DECell using several representative test functions.The optimization procedures of design parameters of new brake model are presented in Section 5.And in Section 6, we study a case of designing brake using DECell.This paper concludes with a summary in Section 7.

Multiobjective Optimization Design Model of Drum Brake
Drum brake is one of the most commonly used brakes in vehicle design; it can be categorized into leading-and trailing-shoe brake, two-leading-shoe brake, two-trailing-shoe brake, and servo brake concerning the arrangement of the brake shoes.The optimization design object in this paper is the leading-and trailing-shoe brake, which is shown in Figure 1.
Figure 1 shows the structural parameters and force diagram of the leading-and trailing-drum brake only left part of the structure is presented .With the effect of braking force P , the two brake shoes at both sides rub the drum brake to generate frictional resisting torque and thus to brake.In Figure 1, F 1 is the pressing force of the brake shoe, while R f is the action radius of the frictional force.δ represents the angle between x axis and the active line of F 1 .The main structural parameters are the starting angle of brake friction plate θ 0 , wrap angle of friction plate θ , the distance from braking force P to the brake center a , the central location of the brake shoe pin c, c , and so forth.

Design Variables
One has X T x 1 , x 2 , . . ., x 11 θ 0 , θ, b 1 , R, c , c, a, b 2 , e, d, p , where θ 0 is the starting wrap angle of the friction plate, θ is the wrap angle of the friction plate, b 1 is the friction plate width, R is the brake drum radius, c and c are the central sites of the brake shoe pin, a means the distance from the braking force to the brake center, b 2 is the brake width, e is the thickness of the brake drum, d is the brake cylinder diameter, and p is the braking oil pressure.

Objective Functions
The optimization target is the comprehensive performance of the brake.And the subobjectives are maximizing the braking efficiency factor and minimizing the volume and temperature rise of the drum brake.The subobjectives are presented as follows.

Braking Efficiency Factor
One has the following 14 : where M μ is the braking output torque, P is the input force, and R f represents the action radius of the frictional force.Therefore, efficiency factor of the brake with tight loose shoe is as follows: where

2.3
μ is the friction factor between the brake drum and the friction plate.
With the same braking force, the bigger the braking efficiency is, the bigger the braking torque is, and the better the brake performs, hence the more efficient the brake is.Therefore, the enhancement of braking efficiency is of great significance to the safe running of vehicle.On this basis, the braking efficiency factor is selected as the optimization objective, and the objective function is where K 1 represents efficiency factor of the brake with tight shoe and K 2 represents efficiency factor of the brake with loose shoe see 2.2 .

The Volume of the Drum Brake
With the same function and efficiency, the smaller the brake is, the fewer materials and the less space it will take up.Hence, the volume of the brake drum should be as small as possible, and it is taken as the optimization objective as follows:

3 Temperature Rise of Drum Brake
It can be regarded that the heat energy generated through braking is absorbed by the front brake and rear brake due to the fast braking process.The heat energy is then distributed to the front brake and rear brake in accordance with their distribution coefficients.The brake drum is installed at the rear wheel, so the rear brake temperature can be calculated by the following formula: where m a is the total mass of the vehicle kg ; m d is the total mass of each drum brake kg ; m h is the total mass of the heating metal materials linking with the drum brake; C d is the specific heat capacity J/ kg•K of material of the drum brake, while C h is the specific heat capacity J/ kg•K of the heating metal materials linking the drum brake; v a is the initial velocity m/s of the brake drum; β is the braking force distribution coefficient.For the convenience of calculation, the brake mass m d and the total mass of the metal materials linking the drum brake m h are selected through the following formula: where ρ d and ρ c are the density of the brake drum and the density of the brake shoe, respectively.Let ρ d ρ c 7570 kg/m 3 ; e and Δ are the thickness of brake drum and thickness of brake shoe, respectively, where Δ 10 mm.

Constraint Conditions 1 Self-Locking Constraint
The possibility of the self-locking for the brake shoe should be taken into consideration when designing the drum brake.To avoid self-locking, the following constraint should be satisfied: 2

Constraint of the Maximum Pressure on the Liners
The maximum pressure on the frictional liners should be less than the defined pressure value: 3

Constraint of the Uniformity of Pressure Distribution
The nonuniform distribution of pressure would cause uneven wear of the frictional liners, which would change the contact area, location, and pressure characteristic of the friction interface.Then the braking efficiency and stability is greatly lowered.The constraint is as follows: 4

Constraint of the Wearing Characteristic of Frictional Liners
During emergency braking, the specific energy dissipating rate of the frictional liners should be less than the defined rate: where E u is the upper limit of the specific energy dissipating rate of rear wheels, while the general upper limit of the specific energy dissipating rate of the brake drum is 1.9 W/mm 2 ; m a is the total mass of the vehicle; R and b 1 are measured in millimeter and braking time t v a /j; v a is the initial velocity m/s ; j is the braking deceleration m/s 2 , and in general j 0.6g, where g is the gravity acceleration. 5

Constraints of the Friction Area of Brake Lining
While braking, the frictional characteristic would be better if the positive pressure and energy load of the brake shoe lining area get lower.In cars, the general value range of the friction area of drum brake lining should obey the following constraints: where A p b 1 Rθ mm 2 .

Heating Load Constraint
To calculate the heating load of the brake, the temperature rise of brake from the initial speed of 30 km/h to full stop should be less than the defined value Δt 15 K 15 : 7

Constraints of the Radius of Brake Drum
Since the diameter of car brake drum is usually 64% ∼ 74% of rim diameter 10 , thus, the ratio of car brake drum radius R against rim radius r d should be in the range: where r d is the rim radius mm .
Mathematical Problems in Engineering 7 8

Constraints of the Liner Width
When choosing the liner width, the pressure in a unit area should not be oversized since too large width would cause nonuniform contact.In general, the ratio of liner width against brake drum radius should meet the following constrains: 2.15

Differential Cellular Multiobjective Genetic Algorithm
In order to solve the above-mentioned multiobjective optimization model much more efficiently, we introduce differential evolution DE strategy into the cellular genetic algorithm and thus form a new algorithm called differential cellular multiobjective genetic algorithm DECell , which is able to minimize the distance from gained Pareto front set to exact optimal front while maintaining good uniformity and distribution range.

Differential Evolution Strategy
DE strategy is well suited for solving multiobjective problems of multidimension and nonlinear as the main feature characterizing it is performing global parallel direct search by using the information of the distance and direction among the individuals in the population, and what is more, it can be easily implemented 16, 17 .The detailed operation process states as follows.
As for each individual x i,D , i ∈ {0, 1, . . ., N − 1}, i ∈ Z, N means the population size, and the definition of the operation is: where D is the number of the decision variables; x r1,D , x r2,D , and x r3,D are the three selected parents; r1, r2, and r3 are the index positions of the parents; and F means the scaling factor.The previous formula can be explained as follows: select two parent vectors x r2,D and x r3,D randomly; get the difference value of the two vectors; add the zoomed difference to the third parent vector x r1,D , and then the individual vector after crossover is v i,D v i,1 , . . ., v i,j , . . ., v i,D ; as for each component x i,j , which is in the parent vector x i,D x i,1 , . . ., x i,j , . . ., x i,D locating in the ith index position, generate randomly a number rand j ∈ 0, 1 ; make a comparison between rand j and the crossover factor CR or between j and K to judge whether to replace x i,j with v i,j ; then get the new individual u i,D u i,1 , . . ., u i,j , . . ., u i,D , where x i,j , otherwise.

DECell Algorithm Principle
Cellular genetic algorithms cGAs are a kind of decentralized genetic algorithms GAs in which the population is structured in such way that the interaction of the individuals composing it is limited to a certain subset of individuals in the immediate vicinity of their location.Usually, the population is arranged in a 2D toroidal mesh in cGAs, and individuals are allowed to interact only with their nearby neighbors see Figure 2 in the dotted box .This endows the cGA with useful properties for the optimization and also facilitates a parallel implementation because of its inherent parallel design 18 .
Based on cGA, the detailed principle of DECell is presented in Figure 2. First, the population is structured in a two-dimension annular gird.Since cGAs are based on the concept of neighborhood, each individual namely, current individual has a set of close individuals or neighbors according to some regular structure like Moore structure, von Neumann structure .In this paper, the neighborhood structure uses Moore structure, where current individual has eight nearby neighbors see Figure 2 .We select two better ones from the neighbors of current individual based on fitness.The three individuals containing the two selected individuals and current one compose the parent individuals, on which crossover operation is then performed by using DE strategy.An offspring is generated after mutation operation.If the offspring dominates the current individual, or both the offspring and the current individual are nondominated, but the offspring is with larger crowding distance than the current individual, then it will replace the current individual.Exploration and diffusion of solutions to the rest of the population happen in a smooth way, due to the continuous overlap existing among the different neighborhoods, which make high-quality solutions propagate over the population.
In the genetic operation process for each cellular individual, the nondominated individuals are stored in an external archive at the same time, and all the individuals in an external archive are ranked according to their crowding distance.If the number of the superior individuals exceeds the specified capacity of the archive, the individuals with smaller crowding distance will be removed.At the end of each generation, some of the individuals in the external archive are chosen randomly to replace the individuals of the same number in the original population so as to update the population constantly.By this method, the nondominated individuals in the external archive are able to obtain a reasonably good approximation to the Pareto optimal front when maintaining the diversity.The main pseudocode of DECell is in Pseudocode 1.

DECell Performance Test and Results Analysis
To assess how competitive the DECell is, a set of unconstraint, multivariable, and multiobjective test functions i.e., DTLZ family 19 and two constraint multiobjective test functions i.e., the Viennet2 and Tanka 20, 21 are adopted.DECell, NSGAII 22 , MOCell 23 , PESA2 24 , and PAES 25 are used simultaneously to compute the above-mentioned test functions, respectively.The performances of the five algorithms are then compared and analyzed under the same assessment system.

Performance Assessment Indicators of Algorithms
This paper uses three assessment indicators: generational distance 26 , spread 23 , and hypervolume 27 .

Generational Distance
Generational distance is used to measure the distance between the obtained Pareto front and the Pareto optimal front.The computational formula is as follows: where n is the number of the obtained Pareto front and d i is the Euclidean distance between the ith solution and the nearest solution in the Pareto optimal front in the objective space.It is obvious that all the obtained solutions are exactly the Pareto optimal solutions if GD 0. Therefore, the smaller the value of GD is, the better it converges to the Pareto optimal solutions.

Spread
The spread measurement proposed by Deb is an indicator that measures the distribution and spread of the obtained Pareto front.The calculation formula is as follows: 3 Hypervolume Hypervolume is adopted to calculate the volume covered by the individuals of the obtained Pareto front.It is defined as where Q is the amount of the obtained Pareto front; as for each individual i in the Pareto front, v i is the hypercube formed by the reference point W 0, .., 0 and the individual i.Larger values of HV indicate that the obtained Pareto front individuals can be broadly distributed in the objective space.

Computational Results and Analysis
The parameters settings of DECell, NSGA-II, MOCell, PESA2, and PAES are as follows: all the algorithms adopt real coding and polynomial mutation.NSGAII, MOCell, PESA2, and PAES employ SBX 28 crossover operator, while DECell adopts DE crossover operator.Let F 0.1  and CR 0.5.The population size is 100, the maximum evaluation generation is 20000, and crossover probability is 0.9, while mutation probability is 1/len, where "len" means the variable length.These five algorithms run 20 times independently for every test function, respectively.Quantitative analysis on the performances of these five algorithms is presented in Tables 1, 2, and 3, in which the values in the dark grey blanks are the optimal values, while those in light grey blanks are the second optimal values.We consider first the generational distance GD .As the test results of the computation on the benchmarks in Table 1 show that, DECell can get the best indicator values in six of the eight multiobjective problems, MOCell and PESA2 can only obtain the best values in one of the eight test functions, respectively.Moreover, NSGA-II and PAES cannot get any best GD value in any of the problems.Therefore, a conclusion can be made that in general DECell outperforms the other four algorithms in terms of convergence when solving highdimensional multiobjective optimization problems MOPs mentioned previously.
Regarding the spread value Table 2 , the results reveal that DECell outperforms the other four algorithms in terms of the diversity of the nondominated solutions along with the obtained Pareto front since it obtains the best lowest spread value in five of the eight multiobjective functions, along with two suboptimal values.The Δ indicators suggest that PAES performs the worst, as it cannot get best value in any problem, only with one suboptimal value.The values of hypervolume indicator HV in Table 3 indicate that DECell gets better results than the other algorithms.DECell obtains the best values the highest values in six of the eight test functions, while NSGA-II gets the best values in two of the eight functions, which denotes that the Pareto front individuals obtained by DECell are able to be distributed more broadly in the objective space.
To compare the performance of algorithms more clearly, we present Figures 3, 4, 5, 6, and 7 the Pareto fronts gained by DECell, MOCell, NSGA-II, PESA2, and PAES when solving DTLZ2, where the grid areas are the Pareto optimal front surface and the dots represent the Pareto front obtained by each algorithm.It is clear that DECell outperforms the other four algorithms concerning distribution and convergence when solving DTLZ2.The PEAS2 has the worst distribution, while MOCell has the worst convergence.And as for the coverage of the gained solutions, DECell outperforms the other ones, which is followed by MOCell and NSGA-II, and next is PAES, while PESA2 is the worst, which matches the results shown in Tables 1, 2, and 3.
Another analysis method is needed to further analyze the influence of the evaluation number on the algorithm performance.The simple way is to observe the evolution of three given indicators during the whole execution of the algorithms from which we can tell what is happening in every generation.Figure 8 shows the convergence progress of GD, HV, and spread when solving DTLZ1 problem with different evaluation number.In this figure, with evaluation number increasing, we can see that the GD values obtained by all the algorithms have converged to the lowest or best value and the HV value obtained by DECell is larger or better than those obtained by the other algorithms during evolution.We also observe that fluctuation in spread value is unstable, but the fluctuation of DECell is smaller than that of the other algorithms, and the spread value obtained by DECell is lower or better than those by the other algorithms with the progress of evolution as well.
A more detailed analysis can be performed if we show the results using boxplots, which constitutes a useful way of depicting groups of numerical data.Figure 9 presents the boxplot distribution of GD, spread, and HV values obtained by the five algorithms in ZTLZ3.We can see that GD and spread values obtained by DECell have been lower or better than the ones obtained by the other algorithms though HV values obtained by these algorithms vary slightly in fact, HV values have no confidence interval .It can be concluded that in general DECell performs better than the other algorithms when solving complicated multiobjective problems, while NSGA-II takes the second place.Though DECell is outstanding in accuracy, it is worse than some other MOEAs in efficiency.That is to say, DECell costs more computational time than NSGA-II, PESA2 for the problems, because DECell is a kind of decentralized cellular genetic algorithms, where individual interacts only with its nearby neighbors and diffusion of solutions to the population happens slowly in a smooth fashion.

Parameter Optimization Design of Brake Based on DECell
The computational procedures of the brake model by using DECell are as follows. 1

Population Initialization
Randomly generate initial population in the constraint range of optimization parameters by real coding.The optimization parameters are proc crossoverInDE r1, r2, r3 // P x1 and P x2 are the corresponding initial position index and terminal position index of the parameters Rand j randInt P x1 , P x2 ; // perform differential crossover operation for j P x1 , j ⇐ P x2 , j j 1 if rand 0,1 < CR or j Rand j v r3 j ,t x r3 j ,t F × x rl j ,t − x r2 j ,t ; //t means current functions evaluations else v r3 j ,t x r3 j ,t ; end if end for return v r3,t ; end proc crossoverInDE

Mutation
This operation aims to avoid local convergence.

Termination Condition
If the maximum evaluation generations are reached, output the results or go back to step 2 .

Case Study
For a vehicle, some of its parameters and its leading-and trailing-shoe brake are shown in Table 4.
To search optimal solutions of the multiobjective problems, Table 5 presents the boundary of design variables.
To evaluate the competitiveness of this algorithm, the most commonly used multiobjective genetic algorithm NSGA-II is also adopted to compute the brake model besides DECell.The parameters settings are as follows: population size M 100, archive capacity of DECell N 100, feedback number C 25, crossover probability p c 0.9, and mutation probability p m 1/var, where "var" is the number of the decision variables.

Qualitative Analysis
After each evolution, the superior solutions gained by NSGA-II are put in the final population, while those obtained by DECell are stored in the external archive.Since the population size and the archive capacity are both 100, each algorithm obtains 100 superior solutions after computation.However, there are still some dominated solutions existing in the 100 individuals, which means that there would be 100 or less non-dominated ones after each computation.
Let the maximum evaluation generations be G 25000.Two algorithms are operated under the same conditions; they both get 100 nondominated individuals, which are shown in Figure 8.
Even though we cannot easily and specifically tell which algorithm performs better in terms of the three objectives from Figure 10, the overall situation shows that the braking efficiency factor obtained by DECell is obviously bigger than that got by NSGA-II.Since the bigger the braking efficiency factor is, the better the braking performance is, it is clear that DECell outperforms NSGA-II for the objective of braking efficiency.In order to make clear the relationship between temperature rise and the volume of the brake drum, the twodimensional graph of the two factors is shown separately in Figure 11.
In Figure 11, the Pareto fronts in upper left are of higher temperature rise and smaller brake drum volume, which means that the smaller the brake drum is, the higher the brake temperature rise will be, that is, the heat dispersion is worse.As for the Pareto fronts in lower  right, they have lower brake temperature rise and larger brake drum volume, which indicates that bigger brake drum volume accompanies smaller brake temperature rise; that is, the heat dispersion is much better.In Figure 11, as for the two factors, the Pareto front individuals gained by DECell distribute more uniform than those of NSGA-II; moreover, it is of better continuity.What is more, more Pareto front individuals obtained by DECell are distributed in the compromise position.It can be concluded that DECell outperforms NSGA-II in terms of distribution and convergence.
Figure 12 presents the relationship between braking efficiency factor and the volume of the brake drum.In the figure, Pareto front solutions obtained by the two algorithms distribute disorderly in 2D multiobjective space, which implies that the efficiency factor is irrelative with the volume of drum brake.That is to say, there is no promoting or conflicting relationship between them.

Convergence Rate Analysis
To further assess the algorithms performance in searching Pareto optimal solutions, different maximum evolution generations, which are got from 1000 to 10000 with the interval 1000, are adopted.For every maximum evolution generation, the aforementioned two algorithms are applied to compute the model of brake for 20 times, respectively, and the averages of the number of the nondominated individuals gained by each algorithm during the 20 times computations are plotted in Figure 13.
Figure 13 shows that DECell can always find more nondominated individuals than NSGA-II with the same maximum evolution generations and the former algorithm outperforms the latter in terms of space searching ability, optimal solution saving ability, and exploring ability.

Convergence Accuracy Analysis
To further analyze the computational stability and accuracy of the algorithm, the maximum evolution generation is set as G 25000 and 20 independent computations on the brake optimization design model are performed by NSGA-II and DECell, respectively.Each algorithm gets 2000 Pareto solutions i.e., 2000 sets of objective function values are obtained , which are then statistically analyzed in terms of three indicators: the mean, standard deviation and optimal value of all the values gained for each objective.The results are presented in Table 6.
Each data in the grey blanks presented in Table 6 is the superior value of the two data in the same type.If the brake is of large braking efficiency factor, small volume, and low temperature rise, it is of good performance.
As is shown in Table 6, DECell gets better mean efficiency factor and temperature rise values, except that the average volume of drum brake optimized by DECell is a little larger than that of NSGA-II the volume is still acceptable .The results denote that DECell outperforms NSGA-II in finding much better design parameters of brake.Meanwhile, as for the deviation of three indicators, two of them obtained by DECell are lower than those obtained by NSGA-II, which means that DECell is of low volatility and good stability.As for each objective, all the optimal values found by DECell are better than those got by NSGA-II.It indicates that DECell can converge to the much better Pareto solutions and it is of higher accuracy.
Therefore, it can be concluded that differential cellular multiobjective genetic algorithm performs better than NSGA-II in terms of solution space searching ability and speed for finding much more and better solutions when solving the brake optimization design model.Moreover, the solutions could be kept effectively.What is more, this algorithm is of good stability and high accuracy, in which much more good design parameters of brake could be easily found for designers to choose.

Conclusions
This paper established a multiobjective optimization design model of drum brake with the goals of maximizing the efficiency factor of braking, minimizing the volume of drum brake, and minimizing the temperature rise of brake, in order to better meet the requirements of engineering practice.
We present a new optimization algorithm by introducing the differential strategy into the canonical cellular genetic algorithm to solve the model effectively, which is called differential evolution cellular multiobjective genetic algorithm DECell .In DECell, all the population individuals are structured in a bidimensional grid, and each individual may only interact with its neighborhood.Besides, three parent individuals are selected for each individual to take part in the crossover operation together.With this method, the obtained Pareto front would be as approximate as possible to the Pareto optimal front.Besides, their distribution uniformity and coverage are well maintained.The evaluation on test function indicates that the proposed algorithm outperforms the other four popular multiobjective algorithms i.e., NSGAII, MOCell, PESA2, and PAES greatly when solving high-dimensional nonlinear multiobjective problems.The experimental computation of DECell and NSGA-II on parameter optimization design of drum brake suggests that the DECell can find much more and better solutions stably.It is an effective algorithm that can be applied to solve the drum brake parameters optimization and other complicated engineering problems.

Figure 1 :
Figure 1: Structural parameters and force diagram.

2 where
d i is the Euclidean distance between every two adjacent individuals in the obtained Pareto front, d is the mean of these distances, d f and d l are the Euclidean distances between the boundary individuals of the obtained Pareto front and the boundary individuals of the Pareto optimal front, and n is the number of the individuals in the obtained Pareto front.A zero value of this indicator means an ideal uniform distribution.Therefore, the smaller this indicator's value is, the more uniform it distributes.

Figure 8 :
Figure 8: Evolution of the three indicators for different algorithms in ZDT1.

Figure 9 :
Figure 9: Boxplot representation of GD, spread, and HV values of the fronts computed by five algorithms in ZTLZ3.

PSEUDOCODE 2 :
Crossover operation pseudocode based on DE strategy.

Figure 10 :
Figure 10: Distribution of the Pareto front.

Figure 11 :
Figure 11: Pareto front set schematic diagram between volume and temperature rise.

Figure 12 :
Figure 12: Pareto front set schematic diagram between volume and efficiency factor.

Figure 13 :
Figure 13: The number of the nondominated individuals obtained by the two algorithms.
Evaluate the offspring individuals.// The superior offspring replaces the present individual.insert position individual , offspring, population ; // Put the non-dominated individuals into the external archive.

Table 1 :
Average avg and standard deviation std of the GD indicator of the Pareto front obtained by five algorithms for different problems.

Table 3 :
Average avg and standard deviation std of the HV indicator of the Pareto front obtained by five algorithms for different problems.

Table 5 :
Domain of decision variables.