Abstract

The present trend in industries is to improve the techniques currently used in design and manufacture of products in order to meet the challenges of the competitive market. The crucial task nowadays is to find the optimal design and machining parameters so as to minimize the production costs. Design optimization involves more numbers of design variables with multiple and conflicting objectives, subjected to complex nonlinear constraints. The complexity of optimal design of machine elements creates the requirement for increasingly effective algorithms. Solving a nonlinear multiobjective optimization problem requires significant computing effort. From the literature it is evident that metaheuristic algorithms are performing better in dealing with multiobjective optimization. In this paper, we extend the recently developed parameter adaptive harmony search algorithm to solve multiobjective design optimization problems using the weighted sum approach. To determine the best weightage set for this analysis, a performance index based on least average error is used to determine the index of each weightage set. The proposed approach is applied to solve a biobjective design optimization of disc brake problem and a newly formulated biobjective design optimization of helical spring problem. The results reveal that the proposed approach is performing better than other algorithms.

1. Introduction

Engineering design optimization has recently received lot of attention from designers so as to produce better designs. It saves time, cost, and energy involved. Engineering design optimization often deals with many design objectives under nonlinear and complex constraints. Also the design problem is subjected to constraints limited by cost, weight, and material properties like strength, design specifications, and availability of resources. The design objectives are often conflicting in nature and hence finding the true Pareto optimal front is difficult. Conventionally, numerous mathematical methods such as linear, nonlinear, dynamic, and geometric programming have been developed in the past to solve engineering design optimization problems. But these methods have many drawbacks and no single method is found suitable for solving all types of engineering design optimization problems.

Even for a single objective design optimization problem involving large number of nonlinear design variables, arriving at a global best solution is not an easy task. The shortcoming of mathematical methods gives way for metaheuristics algorithm to solve engineering design optimization problems. Metaheuristics algorithms have proven to be powerful in solving this kind of optimization since they combine rules and randomness to mimic natural phenomena such as biological systems (Genetic Algorithm), animal behavior (ant algorithm, tabu search), and physical annealing process (Simulated Annealing). GA is a popular metaheuristics algorithm which has been broadly applied to solve various design optimization problems and proven to be successful in finding the global best than traditional methods. These include the design optimization of machine elements such as gears [1], gearbox [2], journal bearing [3], magnetic thrust bearing [4], rolling element bearing [5], and automotive wheel bearing unit [6]. Swarm intelligence is another important concept in many recent metaheuristics algorithms such as particle swarm optimization, firefly algorithm, artificial bee colony algorithm, and cuckoo search.

Particle swarm optimization mimics the group behavior of birds and fish schools [7]. Sabarinath et al. [8] have applied PSO algorithm for solving optimal design of belt pulley system. Yang [9] developed firefly algorithm that imitates the flashing behavior of tropic firefly swarms. Yang and Deb [10] also proposed cuckoo search algorithm based on the behavior of cuckoo birds. Recently Geem et al. [11] proposed a new harmony search (HS) that draws its motivation not from a biological or physical process like most other metaheuristic optimization algorithms, but from an arty one, that is, the improvisation process of musicians in the hunt for a magnificent harmony. Geem et al. [11] explained the similarity between musical performance and optimization that explains the strength and language of HS. The attempt to find the harmony in music is similar to finding the optimality in an optimization process and the musician’s improvisations are similar to local and global search schemes in optimization techniques. HS algorithm uses random search based on harmony memory considering rate (HMCR) and the pitch adjusting rate (PAR) instead of gradient search. Even though the basic algorithm is successful in solving variety of problems still there is a growing diversity of modified HS algorithms that try to find improved performance. Recently Kumar et al. [12] had proposed a new version of HS algorithm, namely, parameter adaptive harmony search (PAHS) to solve standard benchmark functions in single objective optimization successfully. This PAHS algorithm is extended to solve multiobjective design optimization problems by weighted sum approach. In this weighted sum PAHS approach, a performance index based on least average error [13] is used to evaluate the performance of each weightage set.

2. Literature Review

2.1. Previous Work on Disc Brake Optimization

The proposed multiobjective disc brake optimization problem was introduced by Osyczka and Kundu [21]. The authors used the modified distance method in genetic algorithm to solve the disc brake problem and compared their results with that of a plain stochastic method. Ray and Liew [22] used a swarm metaphor approach in which a new optimization algorithm based on behavioural concepts similar to real swarm was proposed to solve the same problem. Yıldız et al. [23] used hybrid robust genetic algorithm combining Taguchi’s method and genetic algorithm. The combination of genetic algorithm with robust parameter design through a smaller population of individuals resulted in a solution that lead to better parameter values for design optimization problems. L16 orthogonal arrays were considered to design the experiments for this problem. The optimal levels of the design parameters were found using ANOVA with respect to the effects of parameters on the objectives and constraints. Yıldız [24] used hybrid method combining immune algorithm with a hill climbing local search algorithm for solving complex real-world optimization problems. The results of the proposed hybrid approach for multiobjective disc brake problem were compared with the previous solutions reported in literature.

Yang and Deb [25] used multiobjective cuckoo search (MOCS) algorithm for solving this disc brake problem by considering two objective functions. The authors extended the original cuckoo search for single objective optimization by Yang and Deb for multiobjective optimization by modifying the first and third rule of three idealized rules of original cuckoo search. The same author [26] used multiobjective firefly algorithm (MOFA) for solving this disc brake problem. By extending the basic ideas of FA, Yang et al. developed multiobjective firefly algorithm (MOFA). Yang et al. [27] successfully extended a flower algorithm for single objective optimization to solve multiobjective design problems. The author solved the biobjective disc brake problem using multiobjective flower pollination algorithm (MOFPA). Reynoso-Meza et al. [28] used the evaluation of design concepts and the analysis of multiple Pareto fronts in multicriteria decision-making using level diagrams. They addressed multiobjective design optimization problem of disc brake by considering the friction surfaces as 4 and 6 to obtain Pareto fronts.

2.2. Previous Work on Helical Compression Spring Optimization

In single objective spring design optimization problem, two different cases were proposed with objective function of volume minimization as case I and weight minimization as case II. This single objective design problem had been solved using many optimization algorithms as two different cases. In this section we will look into a few reviews on case I of a compression spring design problem since the proposed multiobjective helical compression spring optimization problem is the conversion of well-known standard benchmark problem of single objective design optimization, that is, volume minimization of helical compression spring (case I). For both the cases, the design variables are common. The three design variables are the wire diameter , the mean coil diameter , and the number of active coils . But the data type of design variables, objective function, and constraints of these two cases are different. Sandgren [15] used integer programming to solve this problem. Chen and Tsao [16] used simple genetic algorithm to minimize the volume of the helical compression spring. Wu and Chow [17] used metagenetic parameter in GA to solve the same problem. Another improved version of GA called geneAS was used by Deb and Goyal [29] to address this problem. The same problem was solved by discrete version of PSO by Kennedy and Eberhart [20] and by using differential evolution algorithm by Lampinen and Zelinka [19]. Guo et al. [18] used swarm intelligence to optimize the design of helical compression spring. Datta and Figueira [30] used real-integer-discrete-coded PSO for solving this case I problem. He et al. [31] used an improved version of PSO to solve this case I spring problem. Deb et al. [32] proposed NSGA II algorithm for solving some biobjective design optimization of mechanical components. He used case I spring problem and converted it to a biobjective problem by adding one more objective of minimizing stress induced in the spring. In this work, we have converted the case I spring problem into a biobjective problem by adding one more objective of maximizing the strain energy stored in the spring. In all works related with case I spring problem, the previous researchers used FPS unit system. But in this work, all data are converted into SI units and the results of this proposed multiobjective spring problem are presented in SI units.

The main objectives of this work are (i) multiobjective design optimization of the disc brake and compression spring using the weighted sum approach of parameter adaptive harmony search algorithm and (ii) to demonstrate the effectiveness of this proposed algorithm for multiobjective optimization of machine components. So far, the weighted sum approach of parameter adaptive harmony search algorithm has not been tried or experimented for the multiobjective optimization of machine components. In this paper, ability of the algorithm is demonstrated using a performance index based on least average error. The optimization results obtained by using this proposed algorithm are compared with those obtained by previous researchers using other methods.

2.3. Previous Works on HS Algorithm

HS algorithm is widely used in solving various optimization problems such as pipe network design [33], design of coffer dam drainage pipes [34], water distribution networks [35], satellite heat pipe design [36], design of steel frame [37], power economic load dispatch [38], economic power dispatch [39], power flow problem [40], vehicle routing [41], orienteering [42], robot application [43], and data clustering [44]. Mahdavi et al. [45] proposed improved HS algorithm by changing the parameters of HS algorithm dynamically. HS algorithm has been implemented by many researchers for the past few years to solve optimization problems in many fields of engineering, science, and technology. To increase the performance of HS algorithm, several improvements were done periodically in the past. From the literature, it is observed that mainly two improvements had been considered. First improvement is in terms of tuning the parameters of HS algorithm and the other one is in terms of hybridizing the components of HS algorithm with other metaheuristic algorithms. In this work, we are concentrating only with the first improvement.

Mahdavi et al. [45] first attempted an improvement in HS algorithm by proposing a new variant of HS known as improved HS (IHS) algorithm. They changed the variants such as pitch adjusting rate (PAR) and bandwidth (BW) dynamically with generations. Linear increase in PAR and exponential decrease in BW were allowed between prespecified minimum and maximum range.

Kong et al. [46] proposed an adaptive HS (AHS) algorithm to adjust PAR and BW. In this approach, PAR was changed dynamically in response to objective function values while BW was tuned for each variable. Omran and Mahdavi [47] proposed a global best HS algorithm inspired by PSO algorithm to improve the performance of HS algorithm. The difficulty in finding the lower and upper bounds of BW was solved by using global best particle which is the fittest particle in the swarm in terms of objective function values than other particles. A self-adaptive mechanism for selecting BW was proposed by Das et al. [48], known as explorative HS algorithm. BW was recomputed for each iteration while the other parameters were kept fixed. The other improvements in HS variants were reported in the literature [14, 41, 44, 4855]. In all the above cases, only PAR and BW values were changed whereas HMCR was kept fixed. This fixed value of HMCR is the key factor which prevents to get global optimal solution. In order to overcome this difficulty, a parameter adaptive harmony search (PAHS) algorithm was recently proposed by Kumar et al. [12]. A detailed description of the proposed parameter adaptive harmony search (PAHS) algorithm can be seen from [12]. However, for the sake of easiness a brief introduction and a step by step computational procedure for implementing this algorithm are given in the following section. A flow chart showing the Steps 1 to 5 of the proposed PAHS approach is given in Figure 1.

3. Parameter Adaptive Harmony Search (PAHS) Algorithm

The harmony search algorithm (HS) is one of the most recent metaheuristic optimization algorithms, analogous to music improvisation process where the musicians improvise the pitches of their instruments to obtain better harmony. Mahdavi et al. [45] further improved HS in the form of improved harmony search algorithm (IHS) wherein the parameters are changed dynamically. Harmony memory consideration rate (HMCR) and pitch adjustment rate (PAR) are the two important parameters that largely affect the performance of HS algorithm. While searching for global best solution, we make use of HMCR, which in turns requires proper tuning, failing which the solution may stuck with local optima. PAR plays a vital role as far as the local search is concerned. During improvisation process, the fixed value of these parameters may lead to a delay in convergence because of the absence of proper balance of global and local search capabilities. The only solution for this issue may be obtained by changing both HMCR and PAR dynamically. The above issue is addressed in the newly proposed algorithm named as parameter adaptive harmony search algorithm (PAHS) by Kumar et al. [12]. In PAHS, a dynamic change in PAR and HMCR values were proposed consequently by modifying improvisation step of IHS. The value of HMCR is initially kept small in order to explore the entire search space to a larger extent and later on the value of HMCR is increased to confine the search space HM (local search) only. Similarly, the value of PAR is initially kept high in order to modify the solutions either stored in HM or from feasible range. Finally PAR value is reduced to obtain better solutions. Based on the modifications stated above, (1) in linear form for HMCR and (2) in exponential form for PAR had been proposed as given below: where HMCR(gn) is the harmony memory consideration rate for generation gn and and are the minimum and maximum values of harmony memory consideration rate, respectively. Similarly and are the minimum and maximum values of pitch adjustment rate, respectively. NI represents the maximum number of improvisations and gn represents generation number. The bandwidth is given as where BW(gn) is the bandwidth for generation gn and and are the minimum and maximum values of bandwidth, respectively.

They have considered four different cases wherein besides the linear change, two parameters have been allowed exponential change also. PAHS’s variants were applied on fifteen standard benchmark test functions. They have compared the results of PAHS with the other versions of HS as proposed by Geem et al. [11], Mahdavi et al. [45], and Kong et al. [46]. The authors also applied this PAHS to data clustering problem. From their results it is evident that the PAHS-3, that is, linear change in HMCR and exponential change in PAR, has outperformed the results of HS, IHS, AHS, and other variants of HS. PAHS-3 also provided better results in both dimensional and noisy environments. Hence, in this work, PAHS-3 is used to solve multiobjective design optimization problems using weighted sum approach.

3.1. Computational Procedure of PAHS

Step 1 (initialization). In this step, the optimization problem and the algorithmic parameters are initialized. In this paper, the minimization of objective function is considered for the two design optimization scenarios.
The optimization problem is stated as minimize such that is the decision variable set. is the number of decision variables. and are the upper and lower bounds for decision variable , respectively. The parameters of the PAHS algorithm are harmony memory size (HMS), range for harmony memory consideration rate (, ), range for pitch adjustment rate (, ), range for distance bandwidth (, ), and number of improvisation (NI).

Step 2 (initialize harmony memory (HM)). Solution vectors are randomly generated for the size of HM, that is, (HMS) in order to fill the harmony memory. Next sorting of harmony memory is done based on the values of the objective function .

Step 3 (generation of new harmony). Using the procedure given below, a new harmony vector is generatedFor each do/ Proposed equation (1)// Proposed equation (2)//(3)/if then/ memory consideration /begin where if then/ pitch adjustment /begin, endifelse/ random selection /endifdone

Step 4 (update harmony memory). Evaluate the newly generated harmony vector in terms of objective function value. It is compared with that of the initial harmony memory. The better harmony vector is added in the harmony memory and the worst harmony vector is dropped from the harmony memory. This step is the deciding stage of the algorithm whether a new harmony vector is to be included in the harmony memory or not.

Step 5 (check for the termination criterion). If the maximum number of improvisation step (termination criterion) is reached, computation is terminated and the best harmony vector is returned by the algorithm. Otherwise, repeat Steps 3 and 4 till the best harmony vector is obtained.
For solving the two design optimization problems of machine elements (disc brake and compression coil spring) explained in Section 4 using weighted sum approach of PAHS, the harmony memory size (HMS) is set to 20, and the maximum generations (set as stopping condition) are set to 200 generations. The harmony memory consideration rate (HMCR) is varied in the range of 0.7 to 0.99 while the pitch adjustment rate (PAR) is varied in the range of 0.01 to 0.99. The bandwidth (BW) is varied in the range of 0.001 to 1/20(UB-LB). Static penalty method is applied for handling the constraints. A total of 4000 fitness function evaluations were made with this optimization approach in each run. To study the performance of proposed parameter adaptive harmony search algorithm, constrained optimization of input parameters to minimize the objective functions was solved and the best results obtained through the mentioned optimization approach in 40 trials were compared with those reported in the literature.

3.2. Weight Based Multiobjective Optimization Approach

The easiest way of implementing the multiobjective optimization is done through weighted sum method. In this approach, the parameter adaptive harmony search algorithm is modified to facilitate multiobjective design optimization. The implementation of weighted sum approach by assigning suitable weight values to represent the preference in choice is explained by Marler and Arora in multiobjective optimization (MOO) [56]. The implementation of weighted sum approach to artificial bee colony algorithm is addressed in Hemamalini and Simon [57] where the multiobjective optimization is treated as composite objective function. The composite objective function is expressed as follows: where and are the assigned weights. The relationship between the variables and is given by where is chosen in the range of [0-1].

Min max normalization is performed for the objective functions so as to assess the fitness of the composite objective function. The normalization is done for scaling the objective functions within a specified range. The following normalization formula is used [58]: where represents the fitness to be normalized and fitoverall represents the overall fitness. is taken as 0.999 in order to avoid zeroes during normalization process.

3.3. Selection of Weightage Value and Performance Index

An analysis is performed by changing the weight values from 1 to 0 with a step size of 0.1 so as to assign the proper weightage values in (4). It is done to see the importance of each weightage set () towards the objective function. A performance index based on the least average error (LAE) formulated by Naidu et al. [13] is used to evaluate the performance of each weightage set. The LAE is evaluated by where represents fitness 1 value and represents fitness 2 value at weightage set (where ). and represent the overall fitness 1 and fitness 2 values. The performance index based on LAE () is expressed as

4. Problem Formulation

4.1. Multiobjective Disc Brake Optimization Problem

The multiobjective disc brake optimization problem was solved by Osyczka and Kundu [21] using plain stochastic method and genetic algorithms for optimization of disc brake problem. They have shown that genetic algorithm was giving better results compared with that of plain stochastic method. The objectives of the problem are to minimize the mass of the brake and to minimize the stopping time. The disc brake optimization model has four variables (as shown in Figure 2) that are(1), inner radius of the discs, in ,(2), outer radius of the discs, in ,(3)F, engaging force, in ,(4)n, number of the friction surfaces (integer) = .

The objective functions and constraints of the disc brake design optimization model provided by Osyczka and Kundu [21] are defined as follows.

Objective FunctionsMass of the brake: Stopping time:

The constraints are as follows.(1)Side constraints , , , .(2)Geometric constraints(i)Minimum distance between radii: (ii)Maximum length of the brake: (3)Behaviour constraints(i)Pressure constraint (ii)Temperature constraint (iii)Generated torque constraint

4.2. Compression Coil Spring Design: Case I

This problem involves the design of a coil spring, which is to support a constant axial compressive load. It is required to minimize the wire volume of the spring so that it can support a given load without failure. At the same time, it is desired to store maximum strain energy possible. So these two objectives are conflicting in nature. As shown in Figure 3, there are three design variables: the real-valued outside diameter of the spring (), the integer-valued number of spring coils (), and the discrete-valued spring wire diameter ().

The objective, in this problem, is to find the optimal combination of the above three design parameters that provides minimum weight and maximum storage of energy subjected to the constraints.

Accordingly, the biobjective optimization problem is formulated as below:

To minimize the weight of the spring, that is, and to maximize the strain energy stored, that is, where The supplied numerical data for the problem are given in SI units in Table 1 along with the allowable discrete values of the spring wire diameter () in Table 2.

The present work aims to minimize the volume of the spring and maximize the energy stored. Since the objectives are conflicting in nature, modification of second objective, that is, strain energy (SE), is done to get it converted for minimization.

The objective functions are given below:

5. Results and Discussion

The multiobjective PAHS algorithm has been applied for constrained design optimization of input parameters of disc brake and helical compression spring problems. The optimization process has been implemented using Matlab 2009 to run on a PC compatible with Pentium IV, a 3.2 GHz processor, and 2 GB of RAM (Random Access Memory).

5.1. Disc Brake Problem

The weightage set is varied in the range of 0 to 1 with a step size of 0.1. It is evident from the results that single objective optimization is achieved, when the weights are defined at their extreme values. The optimal design parameters of disc brake at different weights are shown in Table 3. When and , minimizing the mass of the brake is the effective objective function and it can be observed that the mass of the brake is minimum but the stopping time value is large. In the other extreme, when and , minimizing the stopping time of the brake is the effective objective function in which the stopping time is minimum but the mass of the brake is large. These results are in accordance with the literature discussed above. By varying the weightage set with a step size of 0.1, we can generate 11 possible combinations. In order to find the best combination of weightage set, a performance index based on least average error (LAE) is used as explained in Section 3.2. The purpose of this performance index is to find the least possible error for both objective functions. The first part of (5) is to find the lowest error ratio for the first objective function, that is, mass of the brake. This is done by subtracting the minimum error value (). Every objective function value of mass of the brake in Table 3 will be subtracted from this particular value and then the resulting value will be divided by the minimum error value () to obtain the ratio. A negative answer would not appear through this method since we are subtracting the least value. The same principle is applied to the second objective function, that is, stopping time of the brake. The results for both the calculations are then averaged to obtain a single index value. Based on the highest value observed in Table 3, the weightage set of and is chosen as the best combination. Figure 4 shows the plot of performance index based on LAE () calculated from the average results of 40 independent trials, each trial having 200 generations for the three different weight sets. It is evident from the plot that the weightage set of and shows better performance with higher values. The best combination is marked by bold font in Table 3.

The statistical measures like mean, best, worst, and standard deviation are calculated from 40 independent trials. Statistical performances are given in Table 4.

The results of the proposed PAHS approach are compared with the previous solutions reported in literature. Yıldız [24] used Hybrid Immune-Hill Climbing algorithm for solving this problem. Osyczka and Kundu [21] used distance method in GA to solve this problem and compared their results with that of the plain stochastic method. Sabarinath et al. [59] used fast and elitist NSGA II algorithm for solving this disc brake problem. Table 5 compares the results obtained by various methods including the proposed approach for solving the biobjective disc brake problem. It is evident from Table 5 that the results provided by the proposed PAHS approach are better than the results of Yıldız [24] and Osyczka and Kundu [21] with less function evaluations. Yıldız [24] also reported the extreme points (i.e., minimum value points for each separate criterion) obtained from three methods in their paper. The extreme points reported for NSGA II algorithm is the same as that of the proposed PAHS algorithm. Comparison of the extreme points obtained by PAHS against other four methods is given in Table 5.

5.2. Compression Coil Spring Problem

The weightage set is varied in the range of 0 to 1 with a step size of 0.1. It is evident from the results that single objective optimization is achieved, when the weights are defined at their extreme values. The optimal design parameters of spring at different weights are shown in Table 6. When and , minimizing the volume of the spring is the effective objective function and it can be observed that the volume of the spring is minimum. In the other extreme, when and , maximizing strain energy stored is the effective objective function in which the strain energy stored is maximum. Here also the concept of performance index based on least average error (LAE) as explained in disc brake problem is applied. Based on the highest value observed in Table 6, the weightage set of and is chosen as the best combination. Figure 5 shows the plot of performance index based on LAE () calculated from the average results of 40 independent trials, each trial having 200 generations for the three different weight sets. It is evident from the plot that the weightage set of and shows better performance with higher values. The best combination is marked by bold font in Table 6.

The statistical measures like mean, best, worst, and standard deviation are calculated from 40 independent trials. Statistical performances are given in Table 7.

Since we are proposing a new biobjective design optimization for spring problem by converting the standard single objective benchmark problem, we are unable to compare the results as in the case of disc brake. But when and , minimizing the volume of the spring is the effective objective function and thus it becomes single objective benchmark problem. The results obtained for single objective, that is, volume minimization for the proposed PAHS approach, are compared with that of the results reported in literature in Table 8. Since the values for comparison are available in inches, the result of the proposed PAHS algorithm is reported in inches as given in Table 8.

For validating our results for biobjective spring problem, the same problem is solved using NSGA II algorithm with the following values of the parameters of NSGA-II algorithm.

Variable type = real variable, population size = 50, crossover probability = 0.8, real-parameter mutation probability = 0.01, real-parameter SBX parameter = 10, real-parameter mutation parameter = 100, and total number of generations = 100. The results are presented in Table 9. It is evident from Table 9 that the results provided by the proposed PAHS approach are slightly better than the results of NSGA II algorithm.

6. Conclusion

In this paper, a recently proposed parameter adaptive harmony search algorithm based on the combinations of linear and exponential changes has been applied successfully to solve multiobjective design optimization problems using weighted sum approach. In this weighted sum approach, a performance index based on least average error (LAE) was used to find the performance of each weightage set. Based on the higher LAE value, optimum weightage selection was determined. First, a standard benchmark problem of biobjective disc brake optimization has been solved and the results are compared with that of a plain stochastic method, simple genetic algorithm, HIHC, and NSGA II. Based on the analysis carried out, the weightage set of and is determined as the best selection for disc brake optimization problem. It is evident from the results that this approach performs better than other algorithms. Hence this approach is extended to solve a newly formulated biobjective design optimization of spring problem and the results are presented. The optimal weightage set of and is determined as the best selection for spring optimization problem. The result of spring problem is compared with the results of NSGA II algorithm and the results are better than NSGA II. This weighted sum approach based on LAE can be extended to any metaheuristic algorithm for implementing multiobjective optimization. Also Pareto optimal solutions can be found by using adaptive weighted sum method as a future extension of the present problem.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.