Pareto Optimization of a Half Car Passive Suspension Model Using a Novel Multiobjective Heat Transfer Search Algorithm

Most of the modern multiobjective optimization algorithms are based on the search technique of genetic algorithms; however the search techniques of other recently developed metaheuristics are emerging topics among researchers. This paper proposes a novel multiobjective optimization algorithm named multiobjective heat transfer search (MOHTS) algorithm, which is based on the search technique of heat transfer search (HTS) algorithm. MOHTS employs the elitist nondominated sorting and crowding distance approach of an elitist based nondominated sorting genetic algorithm-II (NSGA-II) for obtaining different nondomination levels and to preserve the diversity among the optimal set of solutions, respectively.The capability in yielding a Pareto front as close as possible to the true Pareto front ofMOHTS has been tested on themultiobjective optimization problem of the vehicle suspension design, which has a set of five second-order linear ordinary differential equations. Half car passive ride model with two different sets of five objectives is employed for optimizing the suspension parameters usingMOHTS and NSGA-II.The optimization studies demonstrate that MOHTS achieves the better nondominated Pareto front with the widespread (diveresed) set of optimal solutions as compared to NSGA-II, and further the comparison of the extreme points of the obtained Pareto front reveals the dominance of MOHTS over NSGA-II, multiobjective uniform diversity genetic algorithm (MUGA), and combined PSO-GA based MOEA.


Introduction
In recent time, the multiobjective evolutionary algorithms (MOEAs) have gained enormous attention in solving the engineering optimization problems with more than one objective.The multi/many objective optimization problems (MOOPs) differ from their single objective counterparts (SOOPs) in terms of both their problem definitions/statements and methods to solve such problems.MOOPs have the conflicting objectives to be optimized simultaneously and solving these yields a set of optimal or trade-off or Pareto solutions (Pareto front points), whereas SOOPs have a single objective at a given time and the solution is usually a single optimal point.The classical methods such as calculus based methods, gradient based methods, and elimination and interpolation methods and the like for solving SOOPs and methods, for instance, weighted sum method and  constraint method for solving MOOPs are limited to the problems with simple objective functions and constraints.This is due to their nature to get stuck to a suboptimal solution, their convergence dependency on initial guess, and their unsuitability in solving a large variety of optimization problems [1].Alternatively evolutionary algorithms (EAs), or metaheuristics, have had a remarkable success in finding the global optimum of complex problems nearly in all the disciplines of the knowledge.This success of EAs and their nature of using a population of solutions had led the researchers to employing the search techniques of EAs to optimize the MOOPs.Such algorithms to solve MOOPs are referred to as multiobjective evolutionary algorithms (MOEAs) when they use EAs as their basic search techniques and are referred to as simply multiobjective optimization algorithms (MOOAs) when they use any metaheuristics in general.Primarily all the methods to solve MOOPs have two goals to attain.The first is to find the Pareto front solutions as close as possible to the optimal Pareto front and the second is to maintain diversity among the optimal set of solutions [1].
One of the oldest attempts to have employed the EA to form a population based multiobjective evolutionary algorithm was by Schaffer [2] and the method is known as vector evaluated genetic algorithm (VEGA).VEGA does the selection for each objective separately [3], but its incapability in finding all the Pareto front solutions had limited its applications to very few real world optimization problems.These demerits of VEGA were overcome by algorithms which were inspired by the work of Goldberg and Holland [4], such as multiobjective genetic algorithm (MOGA) by Knowles and Corne [5], nondominated sorting genetic algorithm (NSGA) by Srinivas and Deb [6], and niched Pareto genetic algorithm (NPGA) by Horn et al. [7].These new algorithms introduced the concept of Pareto optimality into the selection process which awarded them a great success in application to various disciplines of engineering as thoroughly described in [8].Consequently the researchers had started experimenting various ways in assigning fitness values to the populations and in maintaining the diversity of optimal points which led to the development of new state-of-the-art MOEAs, for example, strength Pareto evolutionary algorithm (SPEA2 and SPEA) by Zhou et al. [9] and by Zitzler and Thiele [3], Pareto archived evolution strategy (PAES) by Knowles and Corne [5], Pareto envelope-based selection algorithm (PESA and PESA-II) by Corne et al. [10] and Horn et al. [7], an elitist based nondominated sorting genetic algorithm-II (NSGA-II) by Deb et al. [11], and a multiobjective evolutionary algorithm based on decomposition (MOEA/D) by Zhang and Li [12] and its improved versions described in [9].In addition, swarm intelligence based search algorithms for multiobjective optimization [13][14][15] have also been developed and applied to a variety of MOOPs.Several more recently developed multiobjective algorithms [16][17][18][19] based on the search techniques other than the ones used in genetic algorithm [4] or EAs and particle swarm optimization [20] have also shown a great potential in obtaining the Pareto front closest to the true Pareto front.
In this study a nonlinear MOOP of optimizing the vehicle suspension system design having conflicting objectives is optimized.The goal is to optimize the design of passive suspension system of the passenger road vehicle by using a 2-dimensional pitch-heave ride model, also known as half car model.This half car ride model had been used by many researchers as an MOOP to test the performance of their proposed MOOAs, for instance, multiobjective uniform diversity genetic algorithm (MUGA) by Jamali et al. [21] and hybrid PSO and GA algorithm for MOOPs by Mahmoodabadi et al. [22].The following paragraph illustrates the MOOP of half car ride model.
Vehicle suspension system has essentially three objectives, good ride comfort, good road holding or handling, and limited suspension deflection or suspension rattle space.These objectives are conflicting with one another and the designers developing such products are always seeking an optimized passive suspension system fulfilling largely all three conflicting objectives simultaneously.In the present study the multiobjective optimization problem is solved by employing a lumped mass-spring-damper vibration model for the half car vehicle suspension system and then applying MOOAs to get the desired set of optimal design solutions, that is, trade-off or Pareto optimum design points.The employed 2D half car model has five degrees of freedom with one lumped sprung and two lumped unsprung masses (front and rear tires) and thus has a set of five second-order, linear, ordinary differential equations of motion.
In the present paper a multiobjective heat transfer search (MOHTS) algorithm for solving MOOPs is proposed and has been applied to the MOOP of optimizing the suspension parameters of the half car model with driver's seat having total five degrees of freedom as used in [21][22][23] to evaluate the potential of the proposed MOHTS algorithm.Moreover, optimization results of MOHTS are compared with the ones obtained from the popular MOEA named NSGA-II and with the other algorithms which were developed and applied to solve the multiobjective half car optimization problem.The proposed algorithm (MOHTS) employs a very recently developed metaheuristic algorithm named heat transfer search algorithm [25] as a base search technic.Additionally, for sorting of the population an elitist nondominated sorting approach and for maintaining the diversity among multiple optimal solutions of each nondominated front the crowding distance approach of NSGA-II are adopted.
The second section briefly introduces the basic concepts and definitions pertaining to multiobjective optimization.The third section describes the overview of the search technic of heat transfer search (HTS) algorithm, elitist based nondominated sorting and crowding distance approach of NSGA-II, and then explains the architecture and working of the proposed MOHTS algorithm.The fourth section of the paper illustrates two different half car multiobjective optimization problems which are solved numerically, elucidating the potential of the proposed MOHTS algorithm.Firstly the half car MOOP with five conflicting objective functions as used in [21][22][23] is employed in numerical study 1 to test and compare the performance of MOHTS with NSGA-II, MUGA [21], and with the combined PSO and GA based MOEA [22]; secondly another half car MOOP with a different set of five more realistic conflicting objective functions from the ones used in numerical study 1 is formulated and tested on the platforms of MOHTS and NSGA-II.

Terminologies in Multiobjective Optimization
A multiobjective optimization problem can be stated as follows [1,26].
Here,  is the vector of design variables known as design vector with its bounds as shown in (4);  forms an dimensional Cartesian space called design space  created by representing each variable with one coordinate axis. is the vector of objectives or decisions known as objective vector which is to be minimized meaning all the objectives are to be minimized simultaneously.Unlike the single objective optimization problem (SOOP), MOOP creates a multidimensional objective space in the present case it is an dimensional objective space with Cartesian each axis as one objective function.Equations ( 2) and (3) show the inequality and equality constraints, respectively, which must be satisfied by the design vector in order to create the feasible solutions.The set of all feasible solutions is known as the feasible region .In general, there exists no solution vector  which can minimize all the  objective functions simultaneously; hence a concept called Pareto optimum solution which is used in solving MOOPs is now introduced.

Pareto Dominance and Optimality.
For the multiobjective minimization problem as stated above, the set of optimum solutions (Pareto optimal design vectors) includes those vectors which cannot minimize any member of the corresponding objective vector without deteriorating another member [3].To illustrate the Pareto optimality mathematically, consider two design vectors , V ∈ ; then vector  is said to be dominant over vector V (denoted by  ≺ V) iff ∀ ∈ {1, 2, . . ., },   () ≤   (V) ∧ ∃ ∈ {1, 2, . . ., }   () <   (V).In a given set of solutions, the design vectors that are nondominated by any other design vector of the same set are known as nondominated regarding that set.Furthermore, the design vectors which are nondominated over the entire feasible region  are known as Pareto optimal solutions and these solutions collectively form a set called Pareto optimal set (PS) or Pareto optimal front (PF).

Solving MOOPs by Extending EAs and Metaheuristics.
Evolutionary or population based algorithms for solving SOOPs yield a set of solutions known as populations instead of a single solution after one iteration.This nature of metaheuristics in obtaining multiple solutions in each iteration makes it suitable to solve MOOPs.As discussed in the first section, MOOPs have primarily two goals: the first is to obtain the Pareto front as close as possible to the optimal Pareto front and the second is to preserve the diversity between the Pareto optimal solutions.The first aim can be achieved by choosing an appropriate fitness assignment scheme that prefers nondominated solutions, and the second can be accomplished by using an appropriate strategy that preserves the diversity among the solutions of each nondominated front.For instance, NSGA [6] used nondominated sorting as a fitness assignment procedure in which each individual from the population is compared with the others to find its nondominancy and thus to obtain the first Pareto front this will be followed by sorting all the nondominated individuals from the population and this will be repeated until the entire population is sorted in different fronts.In addition, for maintaining the diversity among the solutions in a front, a front wise sharing function method was used in NSGA which calculates the Euclidean distance for each solution in the Pareto front from another solution in the same front.The performance of NSGA was further improved in terms of converging the solution more close to the Pareto front by employing an elitist based nondominated sorting, in maintaining a widespread of solutions by applying crowding distance approach instead of sharing function method, and in reducing the computation time.This improved NSGA was named NSGA-II [11].The proposed algorithm (MOHTS) adopts the fitness assignment and diversity preserving mechanisms of NSGA-II over the search method of HTS for solving MOOPs.

Multiobjective Heat Transfer Search Algorithm
This section explains the working of the proposed MOHTS algorithm.The different elements of the MOHTS algorithm are basic search technic which is HTS and nondominated sorting method and diversity preserving crowding distance approach of NSGA-II.These elements which constitute the proposed algorithm are briefly explained first followed by the MOHTS procedure section described by the pseudocode for better understanding and software implementation for the readers and users of this algorithm.

Basic Heat Transfer Search Algorithm.
A very recently developed metaheuristic algorithm called heat transfer search (HTS) algorithm [25] which is based on the thermodynamic mechanism of the heat transfer has been employed as a basic search strategy in the proposed MOHTS algorithm.The search agents in the algorithm are the molecules which interact with one another and with the surrounding for gaining the thermal equilibrium.Unlike other metaheuristic algorithms, the HTS has been developed by putting more efforts in setting a well-balanced trade-off between the intensification (exploitation) and diversification (exploration).The algorithm is divided into three phases as conduction, convection, and radiation; and in the course of entire search process each phase is performed with an equal probability.All the three phases intensify the search space in initial iterations and later diversify during remaining iterations.The following sections briefly explain all three phases.
Else (  ) < (  ); then Else, for the remaining populations (number of iterations  > maximum number of iterations/CDF), if (  ) < (  ) then Else (  ) < (  ); then Here,  and  indicate randomly picked solutions and vary from 1 to  (population size); and  indicates the design variable index ranging from 1 to V (total number of design variables), which is also selected randomly, and CDF is the conduction factor.Moreover,  and   are probability value for conduction phase (this value is between 0 and 0.3333) and random number between 0 and 1, respectively.For the first half of the total generations (CDF should be taken as 2), a function evaluation comparison is made between two randomly selected solutions; the inferior of the two is replaced by ( 5) and (6).For the latter half the inferior solution is replaced by (7) and (8).Thus the switching from ( 5) and ( 6) to ( 7) and ( 8) obtains both the intensification and diversification for maintaining the algorithm's potential to seek an optimum solution.

Convection Phase.
In this phase, a thermal equilibrium is attained as system's mean temperature ( mean ) interacts with the surrounding temperature (  ).The latter temperature is assumed to be the best solution and the solution is updated as follows: If current iteration number  ≤ maximum number of iterations/CVF (initial populations) then Else, for the remaining population, It should be noted that, along with each solution   , each design variable is also updated during this phase.Here  is taken in the range from 0.6666 to 1, and CVF (convection factor) is taken as 10.That is to say, for the first onetenth of the iterations, the algorithm explores the feasible design space and for the remaining iterations it exploits the same, but in actual practice this condition does not take place as exploration (diversification) and exploitation (intensification) are taken care of by the last terms of ( 9) and (10).

Radiation Phase.
In this phase, system (solution) communicates either with the surroundings which is the best solution or within the system which is some other solution to attain the thermal equilibrium.The solutions are improved as described below.
If number of iterations  ≤ maximum number of iterations/RF (initial populations) and if (  ) < (  ) then Else (  ) < (  ); then Else, for the remaining populations, if (  ) < (  ) then Else (  ) < (  ); then Here  and  are randomly selected solutions for which their fitness values are compared with each other and that determines how each solution is updated using ( 11), ( 12), (13), and (14).The range of  in this phase is from 0.3333 to 0.6666, and radiation factor is taken to be 2.
All three phases of the HTS algorithm are employed as basic search technic in the proposed multiobjective heat transfer search (MOHTS) algorithm.Notice that the proposed MOHTS algorithm uses the values of the constants CDF, CVF, and RF as used in [25].Readers are advised to refer to [25] for clear justification behind using these particular values.

Elitist Nondominated Sorting and Crowding Distance
Approach.Elitist nondominated sorting method and diversity preserving crowding distance approach of NSGA-II are introduced in the proposed MOHTS algorithm for sorting of the population in different nondomination levels with computed crowded distance.Firstly for each solution obtained from the basic search method (HTS) or from initially generated random population  0 , all the objectives from the objective vector  are evaluated.In addition, a domination count   defined as number of solutions dominating the solution  Initialize   ,   , CDF, CVF, RF and Rand Generate  0 randomly in S and evaluate F for the generated  0 Sort the  0 based on elitist non dominated sort method and find the NDR and fronts Compute the CD for each front For  = 1 to   // MOHTS main loop starts For  = 1 to   If Rand ≤ 0.3333 // Conduction phase Get   from equations ( 5)-( 8) Else if 0.3333 < Rand ≤ 0.6666 // Radiation phase Get   from equations ( 11)-( 14) Else // Convection phase Get   from equations ( 9)-( 10 and   which is a set of solutions dominated by solution  are calculated.Secondly, all the solutions  are assigned a domination count zero and are put in first nondominated level also known as Pareto front (PF) and their nondomination rank (NDR  ) is set to 1. Thirdly, for each solution  with   = 0, each member  of the set   is visited and its domination count   is reduced by one.While reducing   count if it falls to zero the corresponding solution  is put in second nondomination level and NDR  is set to 2. The procedure is repeated for each member of second nondomination level to obtain the third nondomination level, and subsequently the procedure should be repeated until the whole population is sorted into different nondomination levels.
In crowding distance approach for maintaining diversity among the obtained solutions firstly the population is sorted according to value of each objective function in ascending order.An infinite crowding distance is then assigned to the boundary solutions,  = 1 and  = , of each objective.Here  is the total number of solutions in a particular nondominated set.The boundary solutions are the minimum ( = 1) and maximum ( = ) function values.Except the boundary solutions, all the other solutions of the sorted population ( = 2 to  − 1) for each objective  ( = 1, 2, . . ., ) are assigned the crowding distance (   ) as In (15) the right hand side term is the difference in values of objective function  for two neighbouring solutions ( + 1 and  − 1) of solution .Now each solution  is assigned two entities, nondomination rank NDR  and crowding distance CD  .A crowded comparison operator (≺  ) is used as follows to compare two solutions ( and )  ≺  , if (NDR  < NDR  ) or ((NDR  = NDR  ) and (CD  > CD  )).That is to say, between two solutions, the one with the lower nondomination rank is preferred and if both the solutions have the same nondomination rank then the one with the higher crowding distance is preferred.

MOHTS Procedure.
The procedure of the proposed MOHTS algorithm has been shown in Pseudocode 1. Firstly the parameters such as population size ( pop ), termination criteria, here being the maximum number of generations ( gen ), conduction, convection, and radiation factors (CDF, CVF, and RF, resp.) are initialized.Secondly a random parent population   in feasible region  is generated and each objective function of the objective vector  for   is evaluated.Next, elitist based nondominated sorting and crowding distance computation as explained in earlier section is applied on   .Thirdly HTS algorithm is employed to create the offspring population   , which is then merged with   to form the merged parent-offspring population   .This   is sorted based on elitism nondomination, and based on the computed values of NDR  and CD  the best  pop solutions are updated to form a new parent population.This process is repeated till the maximum number of generations (iterations) is reached.It should be noted that the same algorithm can also be used with the termination criteria set on the basis of number of function evaluations.Since nondominated sorting and crowding distance assignment of MOHTS area is adopted from NSGA-II, the computational complexity of MOHTS is also the same as NSGA-II which is (( pop ) 2 ), where  is the total number of objective functions and  pop is the population size.
Figure 1: Half car passive ride model with five degrees of freedom, adopted from [23].

Multiobjective Optimization of Half Car Ride Model
A half car model of the vehicle's passive suspension system with five degrees of freedom has been adopted from the literature by Bouazara [23] and is shown in Figure 1.The lumped mass model consists of driver's seat mass   , a sprung mass   , and two unsprung masses as front and rear tires  1 and  2 , respectively.The elastic and damping properties of the driver's seat are represented as stiffness and damping coefficients   and   , respectively; similarly elastic properties of the front and rear tires are depicted as stiffness coefficients  1 and  2 , respectively.The damping properties of the tires are assumed to be minimal and thus are neglected.Sprung mass   along with   rests on the front and rear suspensions.The stiffness and damping coefficients of front suspension are  1 and  1 , respectively, and of rear suspension are  2 and  2 , respectively.In addition to the four heave (vertical linear) degrees of freedom of four lumped masses, a rotational degree of freedom (pitch motion) has been considered as a fifth degree of freedom of the system.In Figure 1,   , ,  1 , and  2 are mass moment of inertia of sprung mass, position of the driver's seat, and positions of front and rear tires with respect to the centre of mass of the vehicle, respectively.In addition, the heave or vertical motion of the vehicle is assigned by the vertical displacements  with subscripts showing displacement of the associate lumped masses.The rotational motion (pitch) with respect to the axis perpendicular to the vehicle's plane and passing through the centre of mass of the vehicle is shown by .In Bouazara [23] model, the variables  1 and  2 are the damping coefficients for active suspension system.Since the present study has adopted this model as a passive half car model, these coefficients are neglected.
Applying the Newton-Euler equations to the multibody spring-mass-damper model leads to the equations of motion of the half car model which are given by The equations of motion given by ( 16) are a set of five, second-order ordinary differential equations.Here,   ,   ,  1 ,  2 , and  are the vertical linear displacements of the driver's seat, the gravitational centre of the sprung mass, front tire, and rear tire, and the angular displacement of the sprung mass with respect to its centre of mass, respectively.In addition,  1 ,  2 ,  1 , and  2 depict the vertical displacements of front and rear end of the sprung mass and the road disturbances on the front and rear tires, respectively.
Practically, the road disturbances are random in nature and thus for simulating the mathematical models they are treated as Gaussian isotropic random field [27], stationary  and ergodic stochastic process, or are modeled as single obstacles such as speed bumps and potholes [24,28,29].Single obstacle is a simple way to model the road excitation; the double bump excitation as used by Bouazara [23] has been used in the present study and is as shown in Figure 2.
The equations of motion of the half car were modeled in Simulink/MATLAB and were solved using fourth-order Runge-Kutta numerical method (RK4) using MATLAB.
The multiobjective optimization of the half car model using the proposed MOHTS and NSGA-II is presented in two numerical studies.The first study adopts the optimization problem from [21,22] with a total of five objectives and seven design variables, whereas the second numerical study considers the more realistic approach in choosing the objective functions for optimizing the passive half car model and the same design variables as used in study 1.The seven design variables taken for the optimization studies are stiffness and damping coefficients of driver's seat (  ,   ), front suspension ( 1 ,  1 ), and rear suspension ( 2 ,  2 ) and the position of the driver's seat from the centre of mass of the vehicle ().The rest of all the parameters of the half car model as shown in Figure 1 are considered as vehicle constants and are the same for both the numerical studies and are given in Table 1.Furthermore, the bounds on the design variables are also kept the same for both studies and are as specified in Table 2.The overall procedure followed for optimizing and simulating the half car ride model is shown in Figure 3.

Numerical Study 1.
The multiobjective optimization problem's statement for study 1 is given as follows.
with no constraints, and the bounds on the design variables are as presented in Table 2.The primary objectives of the suspension system are explained in Section 1 and are developed as measured performance criteria of the suspension system as given in (17).The first objective in (17) is that maximum ride comfort can be obtained by minimizing the absolute vertical linear acceleration of the driver's seat.The fourth and fifth objectives in (17) ensure minimum suspension rattle space and can be achieved by minimizing the front and rear suspension deflections, respectively.Notice that second and third objectives in (17) do not fulfill any primary objective of the suspension system.These are only chosen for study 1 merely to be consistent in the comparison of the proposed MOHTS algorithm with other algorithms which published the results based on the objective function vector as given by (17).More importantly all five selected objectives take the absolute values of their respective quantities and are conflicting with each other; hence they suite the multiobjective optimization quite effectively.Firstly, four 2-objective optimization studies, as taken by Jamali et al. [21], are executed with the pairs of objectives as using MOHTS and NSGA-II.For both MOHTS and NSGA-II, the population size is selected as 80, and the total number of generations is set to be 240.Furthermore for NSGA-II the crossover and mutation probabilities were taken as 0.9 and 0.1, respectively.The obtained Pareto front for the 2objective optimization studies is presented in Figures 4-7   Seat acceleration (m/s 2 )   given by NSGA-II.In addition, points   and   refer to the intermediate trade-off points given by MOHTS and NSGA-II, respectively.These trade-off points are different from the extreme points but show significant design information which would be discussed individually later in this section.

MOHTS NSGA-II
As can be seen from Figure 4, overall MOHTS has produced a Pareto front with the better diversity among the optimal solutions and with the best nondominated set of solutions compared to that of NSGA-II.Extreme left tradeoff point of MOHTS  1 has lower seat acceleration than its equivalent point  1 of NSGA-II but at the cost of higher front tire velocity.Similarly extreme bottom trade-off point  1 of MOHTS has lower front tire velocity than its match point  1 of NSGA-II, though in expense of higher seat acceleration.In addition to the extreme points on the Pareto front, the intermediate trade-off points are equally important.For instance, the optimal point of MOHTS  1 has lower seat acceleration than any optimal point of NSGA-II and lower front tire velocities than many optimal points of NSGA-II.Even more  1 registers 48.5% lower seat acceleration than  1 and 24.65% lower front tire velocity than  1 .Likewise,  1 of NSGA-II has 44.34% lower seat acceleration than  1 and 4.12% lower front tire velocity than  1 .The points  1 and  1 would not have been obtained without using Pareto based multiobjective optimization algorithms; thus such algorithms are able to help the decision maker by giving alternative design trade-off points.Similar to points  1 and  1 , improvements in other intermediate trade-off points in comparison with their respective Pareto fronts' extreme trade-off points are demonstrated in Table 3.It should be noted that the chosen intermediate trade-off points have the highest average percentage improvement (as compared to the rest of the intermediate trade-off points) in both the objectives with respect to the extreme worst points of the same objective.The dashes in Table 3 show that the chosen intermediate trade-off point has no relevance to the selected objective function and hence no improvement or deterioration can be measured.
Furthermore, as can be observed in Figures 5, 6, and 7, MOHTS also yields better nondominated set of Pareto optimal solutions with better and wider spread of the optimal solutions than NSGA-II in the rest of the three 2-objective optimization studies.
The comparison of the extreme trade-off points given by the proposed algorithm MOHTS, NSGA-II, MUGA [21], and a combined PSO-GA based MOEA [22] is presented in Tables 4-7.The points    and    are the extreme left and bottom points given by MUGA, and points    and    are the extreme left and bottom points given by combined PSO-GA based MOEA.As can be observed from the last two columns of the tables, MOHTS either outperforms or performs similar to NSGA-II but clearly outperforms the previous work by Bouazara [23], Jamali et al. [21], and Mahmoodabadi et al. [22].Since Bouazara [23] decomposed the multiobjective optimization into single objective optimization by using the weighted coefficients, his study yielded a single optimum solution.
The response of the seat acceleration, front suspension deflection, and rear suspension deflection of the Pareto  optimal design points obtained from MOHTS and NSGA-II and that of the optimal point found by Bouazara [23] are as shown in Figures 8-10, respectively.Figure 8 shows the seat acceleration response of the half car ride model for the design points  3 of MOHTS,  3 of NSGA-II, and  of the work by Bouazara [23].Figure 9 depicts the front suspension deflection response of the system for the design points which yield the minimum front suspension deflection, and these points are  3 of MOHTS,  3 of NSGA-II, and  of the work by Bouazara [23].Similarly rear suspension deflection response of the system for the optimal points  4 of MOHTS,  4 of NSGA-II, and  of the work by Bouazara [23] is demonstrated in Figure 10.Lastly, 5-objective optimization study, which seeks the Pareto front for simultaneous optimization of all five objectives of (17), with the same algorithm parameters of 2objective optimization study results, is investigated.Since 5objective optimization study deals with the five-dimensional objective space and hence for better understanding of the obtained Pareto front, we use the idea proposed by [30] of using biobjective cross sections (slices) of the Pareto front.Such slices are known as decision maps and they give a clear picture of trade-offs between all the objectives.Such  4) in terms of the values of objective functions and their corresponding design variables between present work (MOHTS), NSGA-II, and previous work by Bouazara [23], Jamali et al. [21], and Mahmoodabadi et al. [22].6: The comparison of the extreme trade-off points (as shown in Figure 5) in terms of the values of objective functions and their corresponding design variables between present work (MOHTS), NSGA-II, and previous work by Bouazara [23], Jamali et al. [21], and Mahmoodabadi et al. [22].7: The comparison of the extreme trade-off points (as shown in Figure 6) in terms of the values of objective functions and their corresponding design variables between present work (MOHTS), NSGA-II, and previous work by Bouazara [23], Jamali et al. [21], and Mahmoodabadi et al. [22].crossover probabilities for both the algorithms as indicated earlier in this section.

Numerical Study 2.
In contrast to the first numerical study, we choose more realistic approach in developing the objective functions for optimizing the design of a half car model which is considered with the same set of design variables, vehicle constants, and bounds on the design variables as chosen in study 1.As stated earlier, vehicle suspension system's three opposing performance indices or objectives are used to evaluate the suspension system's dynamic response.A better method to estimate the ride discomfort is to measure the root mean square (RMS) of the vertical acceleration felt by the passenger/driver.In this study RMS of the seat acceleration is chosen as the first objective.Equivalently road holding is determined through the RMS of tire deflection, and suspension rattle space is evaluated by measuring the RMS of suspension deflection.The RMS of front and rear suspension deflection and RMS of front and rear tire deflections are selected as other four objectives to form the 5-objective optimization problem in this study.The optimization problem is stated as follows.Find the design vector,  = {     1  1  2  2 }  , which minimizes with no constraints.Unlike the previous study, the Pareto fronts obtained from 2-objective studies are superimposed over the corresponding Pareto fronts yielded from the 5-objective studies to shorten the length of the paper.We again use the same slicing  8.It is clear from the table that MOHTS is superior to NSGA-II in minimizing the objectives of 2-objective optimization studies and the minimum values of each objective are shown in bold letters.
As can be observed from Figures 19 and 23, the Pareto front for the 5-objective optimization is quite similar for both MOHTS and NSGA-II and this is also true for the rest all the planes, thus, showing that the potential of the proposed MOHTS algorithm is as good as NSGA-II in finding the Pareto front solutions of the multiobjective optimization problem.Moreover, when 2-objective optimization Pareto fronts from MOHTS are compared for all the shown pairs      of objectives with NSGA-II, it appears that the proposed MOHTS algorithm dominates NSGA-II in obtaining the better nondominated Pareto front with the widespread (diveresed) set of optimal solutions.The reason behind the success of MOHTS in yielding a better Pareto front can be attributed to the search mechanism of HTS, as it takes a good care in maintaining the diversified set of solutions in the search space with its better exploration capabilities, which in turn allows the nondominated sorting and crowding distance approach mechanisms to produce a widespread Pareto front.Furthermore, as exploitation capability has also been wisely implemented in HTS algorithm it helps the dominance sorting and diversity preserving mechanisms in producing a better nondominated set.In contrast to the first numerical study, the average computational time for MOHTS and NSGA-II in this study is 1545 and 1633 seconds, respectively.The reason behind this conflict can be attributed to the uncertainties in the numerical solution of the half car model's dynamic equations of motion.

Conclusions
In this paper, a new multiobjective optimization algorithm called multiobjective heat transfer search (MOHTS) has been proposed which works on the search technic of recently developed metaheuristic HTS and on the elitist nondominated sorting for obtaining different nondomination levels and crowding distance approach to preserve the diversity among the optimal set of solutions of NSGA-II.The proposed algorithm has been applied to the MOOP of the vehicle dynamics to optimize the passive suspension parameters of the vehicle by employing half car lumped mass-springdamper model.Since the chosen MOOP represents a complex optimization problem because of five conflicting objectives and ordinary differential governing equations of the system, it is assumed that it can test the real potential of the proposed MOHTS algorithm.Two optimization studies are presented with a different set of optimization objectives and the nature of the objective functions is higher order linear ordinary differential equations in both studies.Since the proposed multiobjective optimization algorithm (MOHTS) has been developed on the platform of NSGA-II for sorting the population based on the nondominancy and preserving diversity among the obtained set of optimal solutions, its performance at all the stages is compared with NSGA-II.The optimization studies reveal that in 2-objective optimization MOHTS is able to yield the better nondominated Pareto front with the widespread (diveresed) set of optimal solutions as compared to NSGA-II.Furthermore from the analysis of extreme trade-off points obtained from the 2-objective optimization numerical study 1 it was observed that MOHTS produces a better Pareto front of the optimal solutions than multiobjective uniform diversity genetic algorithm (MUGA) and combined PSO-GA based MOEA.In 5-objective optimization studies MOHTS has produced a better and equivalent Pareto front as compared to NSGA-II in numerical studies 1 and 2, respectively.Hence, the performance of the proposed multiobjective optimization algorithm MOHTS was tested on the optimization problem with the linear differential equations of higher order and was found to give better Pareto front as compared to NSGA-II, MUGA, and combined PSO-GA based MOEA.

Figure 3 :
Figure 3: Optimization and simulation procedure followed in optimization of the half car ride model.

Figure 4 :
Figure 4: Pareto front of seat acceleration and front tire velocity in 2-objective optimization using MOHTS and NSGA-II.

Figure 5 :
Figure 5: Pareto front of seat acceleration and rear tire velocity in 2-objective optimization using MOHTS and NSGA-II.

Figure 6 :
Figure 6: Pareto front of seat acceleration and front suspension deflection in 2-objective optimization using MOHTS and NSGA-II.

Figure 7 :
Figure 7: Pareto front of seat acceleration and rear suspension deflection in 2-objective optimization using MOHTS and NSGA-II.

Figure 8 :Figure 9 :
Figure 8: Seat acceleration response of the half car model for different optimal design points.

Figure 16 :Figure 17 :
Figure 16: Pareto front of seat acceleration and front tire velocity in 5-objective optimization and 2-objective optimization using NSGA-II for numerical study 1.

Figure 18 : 5 RMSFigure 19 :
Figure 18: Pareto front of seat acceleration and front tire velocity in 5-objective optimization and 2-objective optimization using NSGA-II for numerical study 1.

8 RMSFigure 22 :
Figure 22: Pareto front of seat acceleration RMS and rear tire deflection RMS in 5-objective optimization and 2-objective optimization using MOHTS for numerical study 2.

5 RMS of seat acceleration (m/s 2 )Figure 23 :Figure 24 :
Figure 23: Pareto front of seat acceleration RMS and front suspension deflection RMS in 5-objective optimization and 2-objective optimization using NSGA-II for numerical study 2.

Figure 25 :
Figure 25: Pareto front of seat acceleration RMS and front tire deflection RMS in 5-objective optimization and 2-objective optimization using NSGA-II for numerical study 2.

Figure 26 :
Figure 26: Pareto front of seat acceleration RMS and rear suspension deflection RMS in 5-objective optimization and 2-objective optimization using NSGA-II for numerical study 2.
) End if End for Merge  0 and   to create   , i.e.   = P 0 ∪   For   apply the elitist non dominated sort method and find the NDR of the fronts Apply the crowding distance for the sorted set of solutions for each objective Based on NDR and CD sort   in descending order Replace  0 with   for first   members of End for // MOHTS main loop endsPseudocode 1: Pseudocode of MOHTS algorithm.

Table 3 :
Design information of intermediate trade-off points.

Table 4 :
[22]comparison of the extreme trade-off points (as shown in Figure3) in terms of the values of objective functions and their corresponding design variables between present work (MOHTS), NSGA-II, and previous work by Bouazara[23]Jamali et al.[21], and Mahmoodabadi et al.[22].

Table 5 :
The comparison of the extreme trade-off points (as shown in Figure

Table 8 :
The comparison of the extreme trade-off points (as shown in Figures18-25) in terms of the values of objective functions and their corresponding design variables between present work (MOHTS) and NSGA-II.