Mathematical Problems in Engineering

Volume 2012 (2012), Article ID 734193, 23 pages

doi:10.1155/2012/734193

Research Article

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

Yi Zhang, Hu Zhang, and Chao Lu

Hubei Key Laboratory of Hydroelectric Machinery Design & Maintenance, China Three Gorges University, Yichang 44300, China

Received 19 July 2012; Revised 14 October 2012; Accepted 15 October 2012

Academic Editor: Jyh Horng Chou

Copyright © 2012 Yi Zhang 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.

Abstract.

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 cellular multiobjective 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.

1. 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.

2. 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: Structural parameters and force diagram.

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 , the two brake shoes at both sides rub the drum brake to generate frictional resisting torque and thus to brake. In Figure 1, is the pressing force of the brake shoe, while is the action radius of the frictional force. represents the angle between axis and the active line of . The main structural parameters are the starting angle of brake friction plate (), wrap angle of friction plate (), the distance from braking force to the brake center (), the central location of the brake shoe pin (), and so forth.

2.1. Design Variables

One has , where is the starting wrap angle of the friction plate, is the wrap angle of the friction plate, is the friction plate width, is the brake drum radius, and are the central sites of the brake shoe pin, means the distance from the braking force to the brake center, is the brake width, is the thickness of the brake drum, is the brake cylinder diameter, and is the braking oil pressure.

2.2. 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 is the braking output torque, is the input force, and represents the action radius of the frictional force. Therefore, efficiency factor of the brake with tight (loose) shoe is as follows:

where

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 represents efficiency factor of the brake with tight shoe and 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:

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 is the total mass of the vehicle (kg); is the total mass of each drum brake (kg); is the total mass of the heating metal materials linking with the drum brake; is the specific heat capacity (J/(kg·K)) of material of the drum brake, while is the specific heat capacity (J/(kg·K)) of the heating metal materials linking the drum brake; is the initial velocity (m/s) of the brake drum; is the braking force distribution coefficient. For the convenience of calculation, the brake mass and the total mass of the metal materials linking the drum brake are selected through the following formula:

where and are the density of the brake drum and the density of the brake shoe, respectively. Let  kg/m³; and are the thickness of brake drum and thickness of brake shoe, respectively, where  mm.

2.3. Constraint Conditions

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:

Constraint of the Maximum Pressure on the Liners

The maximum pressure on the frictional liners should be less than the defined pressure value:

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:

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 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²; is the total mass of the vehicle; and are measured in millimeter and braking time ; is the initial velocity (m/s); is the braking deceleration (m/s²), and in general g, where g is the gravity acceleration.

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 (mm²).

Heating Load Constraint

To calculate the heating load of the brake, the temperature rise of brake from the initial speed of km/h to full stop should be less than the defined value  K [15]:

Constraints of the Radius of Brake Drum

Since the diameter of car brake drum is usually of rim diameter [10], thus, the ratio of car brake drum radius against rim radius should be in the range:

where is the rim radius (mm).

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:

3. 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.

3.1. 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 , means the population size, and the definition of the operation is:

where is the number of the decision variables; , and are the three selected parents; , , and are the index positions of the parents; and means the scaling factor. The previous formula can be explained as follows: select two parent vectors and randomly; get the difference value of the two vectors; add the zoomed difference to the third parent vector , and then the individual vector after crossover is ; as for each component , which is in the parent vector locating in the th index position, generate randomly a number ; make a comparison between and the crossover factor or between and to judge whether to replace with ; then get the new individual , where :

3.2. 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].

Figure 2: Breeding loop in DECell.

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.

proc stepsUp()

population← creatPopulation(); //Create initial population.

archive← creatFront(); //Create a external archive for storing non-dominated solutions.

    while! StopCondition()//Termination condition.

        for  individuals←1  to  population.Size

            //Select the neighbors for the current individual (MOORE is adopted).

            neighborhood← selection(neighborhood);

            //Select two parent individuals from the neighbors.

            parent1← selection(neighborhood);

            parent2← selection(neighborhood);

            //Make sure the two parent individuals are different.

            while  parent1==parent2

                     parent2← selection(neighborhood);

            //Crossover operator with DE strategy and Pseudocode 2 presents

            the detailed operation.

            offspring← crossoverInDE(position(individual),position(parent1),

            position(parent2));

            offspring← Mutation(offspring); // Mutation.

            evalutionFitness(offspring); // 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.

            addToArchive(individual);

        population←replaceIndividuals(population, archive); //Feedback mechanism.

end_proc stepsUp;

Pseudocode 1: DECell pseudocode.

4. 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.

4.1. 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 is the number of the obtained Pareto front and is the Euclidean distance between the th 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:

where is the Euclidean distance between every two adjacent individuals in the obtained Pareto front, is the mean of these distances, and are the Euclidean distances between the boundary individuals of the obtained Pareto front and the boundary individuals of the Pareto optimal front, and 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.

Hypervolume

Hypervolume is adopted to calculate the volume covered by the individuals of the obtained Pareto front. It is defined as

where is the amount of the obtained Pareto front; as for each individual in the Pareto front, is the hypercube formed by the reference point and the individual . Larger values of HV indicate that the obtained Pareto front individuals can be broadly distributed in the objective space.

4.2. 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 and . The population size is 100, the maximum evaluation generation is 20000, and crossover probability is 0.9, while mutation probability is , where “” 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.

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

Table 2: Average (avg) and standard deviation (std) of the spread () 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.

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 high-dimensional 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.

Figure 3: DTLZ2 solved by DECell.

Figure 4: DTLZ2 solved by MOCell.

Figure 5: DTLZ2 solved by NSGA-II.

Figure 6: DTLZ2 solved by PESA2.

Figure 7: DTLZ2 solved by PAES.

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)

(b)

(c)

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

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).

(a)

(b)

(c)

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

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.

5. Parameter Optimization Design of Brake Based on DECell

The computational procedures of the brake model by using DECell are as follows.

Population Initialization

Randomly generate initial population in the constraint range of optimization parameters by real coding. The optimization parameters are .

Selection

The tournament method is used, which selects the parent individuals based on the individual rank and crowding distance. The detailed procedures are as follows: select the neighbors of the current individual; compare the ranks of the neighbors, and remain the individuals with lower ranks since the individuals with lower ranks performs better; if the neighbors are of the same rank, compare their crowding distances, and reserve the ones with bigger crowding distances. At last, two better individuals are selected from the neighbors.

Crossover

Crossover operation is performed to generate more new individuals and improve the space searching ability. DE strategy is introduced to the crossover operation, and the detailed procedures are explained in formulas (3.1) and (3.2). The parameters and are adopted in this operation, where the crossover coefficient (CR) equals 0.5 and the scaling factor equals 0.5. The detailed parameter setting is as per [29]. The code of the crossover flow is presented in Pseudocode 2.

proc crossoverInDE()

// and are the corresponding initial position index and terminal position

index of the parameters

;

// perform differential crossover operation

    for ()

                  if (rand(0,1) < CR or )

            ; //t means current functions evaluations

        else

            ;

        end if

    end for

    return ;

end_proc crossoverInDE

Pseudocode 2: Crossover operation pseudocode based on DE strategy.

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 .

6. Case Study

For a vehicle, some of its parameters and its leading- and trailing-shoe brake are shown in Table 4.

Table 4: Model parameters.

Parameters
The weight

of the vehicle (kg)
Distribution coefficient
Rim radius (mm)
Initial velocity

of brake (km/h)
Friction coefficient

of the friction plate
Specific heat capacity

,   (J/(kg·K))

Values
1560
0.77
177.8
30
0.35
482

To search optimal solutions of the multiobjective problems, Table 5 presents the boundary of design variables.

Table 5: Domain of decision variables.

Variable
 (degree)
 (degree)
 (mm)
  (mm)
  (mm)
  (mm)

Upper limit
10
90
36
120
120
85

Lower limit
45
120
75
180
160
120

Variable
 (mm)
 (mm)
 (mm)
 (mm)
 (MP)

Upper limit
60
90
5
15
6

Lower limit
130
120
16
35
12

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 , archive capacity of DECell , feedback number , crossover probability , and mutation probability , where “” is the number of the decision variables.

6.1. 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 . 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 two-dimensional graph of the two factors is shown separately in Figure 11.

Figure 10: Distribution of the Pareto front.

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

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.

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

6.2. Quantitative Analysis

6.2.1. 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: The number of the nondominated individuals obtained by the two algorithms.

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.

6.2.2. Convergence Accuracy Analysis

To further analyze the computational stability and accuracy of the algorithm, the maximum evolution generation is set as 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.

Table 6: Performance comparison.

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.

7. 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.

Acknowledgment

This work was supported by the National Science Foundation of China (51275274).

References

L. Cheng and W. M. Zhang, “Optimization design of mining trucks & caliper disc brake,” Mining Research and Development, vol. 28, no. 1, pp. 41–43, 2008.

J. Li, “Optimization design of tractor brake based on hybrid PSO,” Chinese Agricultural Mechanization, vol. 4, pp. 93–96, 2011.

Z. H. Li, X. L. Zhang, and L. C. Guo, “Optimal design for caliper disc brake,” Machine Design and Research, vol. 25, no. 2, pp. 83–85, 2009.

L. Chen, “Simulated annealing algorithm for function optimization and application in caliper-type brake,” Coal Mine Machinery, vol. 27, no. 10, pp. 38–41, 2006.

Q. B. Huang, M. G. Zhou, and Y. Wang, “Structural parameters optimization design of drum brake based on genetic algorithm,” Machinery, vol. 11, pp. 24–26, 2006.

P. Jiang and W. J. Huang, “Optimal design of disc brake based on MATLAB,” Mechanical Engineering & Automation, vol. 6, pp. 158–161, 2007.

J. Mi and Y. L. Wu, “Multi-objective optimization design of drum-fashioned brake technology,” Machinery Design & Manufacture, vol. 1, pp. 25–26, 2007.

J. Mi and Y. L. Wu, “The Application of multi-objective optimization method in brake design,” Journal of Mechanical Transmission, vol. 31, no. 5, pp. 71–73, 2007.

G. L. Qin, D. M. Wang, and X. H. Chen, “Optimization design of disk brake based on an improved genetic algorithm,” Modular Machine Tool & Automatic Manufacturing Technique, vol. 7, pp. 101–103, 2011.

H. Wang, Parametric Brake Optimization Based on Ant Colony Algorithm, Jilin University, Jilin, China, 2007.

J. Wu and W. J. Li, “Optimal design of disk brake based on improved swarm optimization,” Machinery Design & Manufacture, vol. 4, pp. 18–20, 2007.

R. H. Yang, “Multi-objective optimization design of auto drum-fashioned brake based on genetic algorithm of MATLAB,” Machine Tool & Hydraulics, vol. 23, pp. 91–93,97, 2011.

M.G. Zhou, Y. Chen, J. D. Zhou, and J. Yi, “Genetic algorithm optimal design and its implementation with MATLAB for drum brake,” Journal of Hubei University of Technology, vol. 24, no. 2, pp. 78–80, 2009.

W. Y. Wang, Automotive Design, China Machine Press, Beijing, China, 2004.

W. X. Liu, Structural Analysis & Design of Vehicle Brake System, Tsinghua University Press, Beijing, China, 2004.

Y. J. Shi, The Cooperative Co-Evolutionary Differential Evolution and Its Applications for Complex Layout Optimization, Dalian University of Technology, Dalian, China, 2006.

J. Durillo, A. J. Nebro, F. Luna, and E. Alba, “Solving three-objective optimization problems using a new hybrid cellular genetic algorithm,” in Proceedings of the 10th international conference on Parallel Problem Solving from Nature: PPSN X, pp. 661–670, 2008.

B. Dorronsoro and E. Alba, “A simple cellular genetic algorithm for continuous optimization,” in Proceedings of the IEEE Congress on Evolutionary Computation (CEC'06), pp. 2838–2844, IEEE Computer Society, Vancouver, Canada, July 2006.

Deb, K. L. Thiele, M. Laumanns, and E. Zitzler, “Scalable test problems for evolutionary multi-objective optimization,” Tech. Rep. 112, Computer Engineering and Networks Laboratory (TIK), Swiss Federal Institute of Technology (ETH), Zurich, Switzerland, 2001.

A. Vicini and D. Quagliarella, “Multipoint transonic airfoil design by means of a multiobjective genetic algorithm,” Tech. Rep. AIAA-97-0082, AIAA, Washington, DC, USA, 1997.

H. Tamaki, M. Mori, M. Araki, Y. Mishima, and H. Ogai, “Multi-criteria optimization by genetic algorithms: a case of scheduling in hot rollingprocess,” in Proceedings of the 3rd Conference of the Association of Asian-Pacific Operational Research Societies with IFORS (APORS'94), pp. 374–381, World Scientifc, Fukuoka, Japan, July 1994.

K. Deb, A. Pratap, S. Agarwal, and T. Meyarivan, “A fast and elitist multiobjective genetic algorithm: NSGA-II,” IEEE Transactions on Evolutionary Computation, vol. 6, no. 2, pp. 182–197, 2002.

A. J. Nebro, J. J. Durillo, F. Luna, B. Dorronsoro, and E. Alba, “MOCell: a cellular genetic algorithm for multiobjective optimization,” International Journal of Intelligent Systems, vol. 24, no. 7, pp. 726–746, 2009.

Corne, D. W. N. R. Jerram, J. D. Knowles, M. J. Oates, and J. Martin, “PESA-II: Region-based selection in evolutionary multi-objective optimization,” in Proceedings of the Genetic and Evolutionary Computation Conference (GECCO'01), pp. 283–290, Morgan Kaufmann, San Francisco, Calif, USA, July 2001.

J. D. Knowles and D. W. Corne, “Approximating the nondominated front using the Pareto archived evolution strategy,” Evolutionary Computation, vol. 8, no. 2, pp. 149–172, 2000.

D. A. van Veldhuizen and G. B. Lamomt, “On measuring multiobjective evolutionary algorithm performance,” in Proceedings of the Congress on Evolutionary Computation (GECCO'00), pp. 204–211, Morgan Kaufmann, Las Vegas, Nevada, USA, 2000.

E. Zitzler and L. Thiele, “Multiobjective evolutionary algorithms: a comparative case study and the strength Pareto approach,” IEEE Transactions on Evolutionary Computation, vol. 3, no. 4, pp. 257–271, 1999.

K. Deb and R. B. Agrawal, “Simulated binary crossover for continuous search space,” Complex Systems, vol. 2, pp. 115–148, 1995.

K. Zielinski, P. Weitkemper, R. Laur, and K. D. Kammeyer, “Parameter study for differential evolution using a power allocation problem including interference cancellation,” in Proceedings of the IEEE Congress on Evolutionary Computation (CEC'06), pp. 1857–1864, IEEE Computer Society, Vancouver, Canada, July 2006.