Optimum suspension design for nonlinear half vehicle model using particle swarm optimization (PSO) algorithm
Ahmet Yildiz^{1}
^{1}Bursa Uludag University, Engineering Faculty, Automotive Engineering Department, Bursa, Turkey
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 27, 2019, p. 4348.
https://doi.org/10.21595/vp.2019.21012
Received 10 September 2019; accepted 17 September 2019; published 30 September 2019
JVE Conferences
This paper is considered with a nonlinear suspension design for half vehicle model by using particle swarm optimization technique. To analyze the ride comfort, a fivedegree of freedom system is built, and it is integrated with the Particle Swarm Optimization (PSO) for optimizing the vehicle vibrations. A multiobjective function is proposed as the sum of the minimum seat and vehicle body acceleration, the minimum suspension deflection and the constraints and the design variables of the optimization problem are selected as the spring and damping coefficients of the front and rear suspension and the nonlinear spring and the linear damping coefficients of the seat. The simulations are carried out and the results are compared with the nonoptimized values. It is demonstrated that the vehicle vibration is decreased significantly with the help of the optimum values of the suspension parameters.
Keywords: vehicle vibrations, nonlinear, seat, suspension, optimization, PSO.
1. Introduction
The vehicle vibration is one of the most important issues for the vehicle design process and many researches has been performed on this important topic in order to obtain high ride comfort. The vibration is mainly caused by road profile it should be compensated by using a welldesigned suspension system. The suspension systems are generally divided into three main categories: Passive suspension, semiactive suspension and active suspension system [13]. Among them, the most common one is the passive suspension system because of its simplicity, efficiency and low cost [4]. In this type of suspension system, the spring and damper have the predefined characteristic, while in the semiactive and active suspension system, the reaction forces are controlled. On the other hand, by adding a nonlinear spring into passive suspension, one can provide a more proper vibration response for different road profiles.
The vehicle vibration is analyzed by using three main mathematical models: Quarter vehicle model (QVM) [5, 6], half vehicle model (HVM) [7] and full vehicle model (FVM) [8]. The simplest of these is the QVM. This model supposes that the vehicle moves only in the vertical direction while HVM has four degrees of freedom (DOF) and it is used to determine the characteristic of the pitch motion as well. If the roll behavior of the vehicle needs to be analyzed, the most complex model, FVM, should be used, which has seven DOF. Furthermore, the seat and human model can be added to these models as well [9].
The optimization of vehicle vibration is very important to achieve better ride comfort and many tools are available for solving such complex engineering problems. The optimization problem can be solved by using a deterministic or a stochastic method that solves the objective function with systematical random inputs. One of the most common stochastic optimization methods is the Particle Swarm Optimization (PSO) technique and it has been applied for many different engineering problems as well as vehicle vibration.
In this paper, a half vehicle model is considered with nonlinear seat suspension and five equations of motion are derived and then they are integrated with Particle Swarm Optimization (PSO) algorithm. The design variables of the optimization problem are the spring and damping coefficients of the front and rear suspension, and the nonlinear spring and the linear damping coefficients of the seat. The objective function of the optimization problem is defined as the sum of the minimum seat and vehicle body acceleration and the minimum suspension deflection. The simulations are performed, and the results are compared with the nonoptimized simulations. It has been found that much lower vibration amplitudes can be achieved by using optimum suspension parameters.
2. Mathematical model of 5 DOF vehicle
In this chapter, the mathematical model of the vehicle is derived to characterize the vehicle vibration with nonlinear seat suspension. As shown in Fig. 1, the half vehicle model with nonlinear seat spring system has five DOF: The vertical motion of the front (${x}_{wf}$) and rear wheel (${x}_{wr}$), the vertical (${x}_{p}$) and pitch motion of the vehicle body ($\theta $) and the vertical motion of the seat (${x}_{p}$).
Fig. 1. A schematic diagram of a half vehicle model
In the mathematical model, the seat spring force ${F}_{sp}$ is considered as a cubic nonlinear term:
where the ${k}_{p}$ and ${k}_{pn}$ are the spring coefficients of the seat suspension. Considering NewtonEuler approach, the equation of motions of HVM are derived as below:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+{c}_{sf}({\dot{x}}_{wf}{\dot{x}}_{v}+{L}_{f}\dot{\theta})=0,$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+{c}_{sr}({\dot{x}}_{wr}{\dot{x}}_{v}{L}_{r}\dot{\theta})=0,$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{k}_{sf}({x}_{wf}{x}_{v}+{L}_{f}\theta ){c}_{sf}({\dot{x}}_{wf}{\dot{x}}_{v}+{L}_{f}\dot{\theta})$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{k}_{sr}({x}_{wr}{x}_{v}{L}_{r}\theta ){c}_{sr}({\dot{x}}_{wr}{\dot{x}}_{v}{L}_{r}\dot{\theta})=0,$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{k}_{pn}({x}_{v}{x}_{p}{L}_{p}\theta {)}^{3}{L}_{p}{c}_{p}({\dot{x}}_{v}{\dot{x}}_{p}{L}_{p}\dot{\theta}){L}_{p}$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{k}_{sr}({x}_{wr}{x}_{v}{L}_{r}\theta ){L}_{r}{c}_{sr}({\dot{x}}_{wr}{\dot{x}}_{v}{L}_{r}\dot{\theta}){L}_{r}=0,$
where the ${m}_{wf}$, ${m}_{wr}$ and ${m}_{v}$ denote the mass of front wheel, rear wheel and vehicle body while ${k}_{wf}$,${k}_{wr}$ and ${c}_{wf}$, ${c}_{wr}$ are the stiffness and damper coefficients of the front and rear wheels respectively. The parameters ${k}_{sf}$,${k}_{sr}$ and ${c}_{sf}$, ${c}_{sr}$ represent the stiffness and damper coefficients of the front and rear suspension system, which are the design variables. Also, the $\theta $ and ${I}_{v}$ indicate the pith angle and the mass moment of the inertia of the vehicle body while ${L}_{r}$, ${L}_{f}$ and ${L}_{p}$ are the distance from the rear suspension, front suspension and passenger axes to the center of gravity respectively. The subscript ${w}_{f}$ and ${w}_{r}$ denote the front and rear wheel.
The input road profile is represented with $u$ and it can be defined as a sinus function and it is given in Fig. 2 by using following equations:
where $\omega $ and ${V}_{v}$ represent the angular frequency and the vehicle velocity which are equal to $\omega =\left(2\pi {V}_{v}\right)/\lambda $ and $\xi =({L}_{f}+{L}_{r})/{V}_{v}$. The amplitude of the road profile is $h=$ 0.05 m, and the period of it $\lambda =$ 11.1 m.
Table 1. Half vehicle parameters with the range of design variables
Parameter

Value/range

Parameter

Value/range

${m}_{wf}$ [kg]

60

${c}_{wf}$ [Ns/m]

10

${m}_{wr}$ [kg]

50

${c}_{wr}$ [Ns/m]

10

${m}_{v}$ [kg]

700

${c}_{sf}$ [Ns/m]

1001500

${m}_{p}$ [kg]

95

${c}_{sr}$ [Ns/m]

1001500

${k}_{wf}$ [kN/m]

190

${c}_{p}$ [Ns/m]

1001500

${k}_{wr}$ [kN/m]

190

${I}_{v}$ [kgm^{2}]

1400

${k}_{sf}$ [kN/m]

1555

${L}_{f}$ [m]

1.2

${k}_{sr}$ [kN/m]

1555

${L}_{r}$ [m]

1.7

${k}_{p}$ [kN/m]

25105

${L}_{p}$ [m]

0.4

${k}_{pn}$ [kN/m]

250010500

${V}_{v}$ [m/s]

22.2

3. PSO optimization algorithms and its procedures
Particle Swarm Optimization (PSO) method is proposed by [10] firstly, which is a stochastic populationbased optimization algorithm. The PSO algorithm is inspired by the group behavior of birds or fishes and it has a huge usage in the engineering applications [11]. A random group of particles is used to search the extreme point of the objective function and they follow the nearest ones in order to the find the best value of the target by means of the reproduced velocities until the finalize criteria is satisfied. About the PSO pseudocodes more details are given as follows:
Step1. Set parameters of PSO and generate the position of particles randomly according to the constrains and give the initial velocities to zero.
Step2. Calculate the objective function and find the particle best position and global best of swarm.
Step3. Update each particle position and velocity with the help of the equations (a) and (b). Restrict the velocity and position of particle in determined ranges to avoid exceeding the limits:
where $i=$1, 2,…, $M$, and $M$ is the swarm size, $j=$1, 2,…, $N$, and $N$ is the variable size, $k=$1, 2,…, ${k}_{max}$, the maximum number of iterations ${k}_{max}=$ 50 is, ${v}_{i,j}$ is the particle velocity, $X$ is the position of the particle, ${X}_{PB}$ and ${X}_{GB}$ are the best position of particles and global best of swarm, and the coefficients ${w}_{max}=$ 0.9 and ${w}_{min}=$ 0.4. ${r}_{1}$ and ${r}_{2}$ are random numbers between (0, 1) while the coefficient of cognitive and social acceleration factors are ${c}_{1}=$ 1.3 and ${c}_{2}=$ 2.8 respectively.
Step4. Update the particle best position and global best of swarm by evaluating the fitness of each particle.
Step5. Finalize the iterations if the requested conditions are achieved. If is it not satisfied, return the step 3.
In order to achieve better ride comfort of the vehicle, there are three main characteristic requirements from the system: Low amplitude of the seat acceleration to improve the ride comfort and low tire deflection for better roadholding ability and low suspension working space. Thus, the objective function is determined as the sum of the maximum acceleration of the seat and the sum of the squares of the constraints. There are three constrain functions regarding the real working conditions of the vehicle suspension system. These are the maximum acceleration of the seat and the vehicle body. Therefore, the objective function (${F}_{obj}$) of the optimization problem and the three constrains (${g}_{1}$,${g}_{2}$ and ${g}_{3}$) are given as below:
${g}_{1}={\ddot{x}}_{p}\left(t\right)4.5\mathrm{}\mathrm{m}/{\mathrm{s}}^{2}\le 0,$
${g}_{2}=\left{x}_{wf}\left(t\right){u}_{wf}\left(t\right)\right0.06\mathrm{}\mathrm{m}\le 0,$
${g}_{3}=\left{x}_{wr}\left(t\right){u}_{wr}\left(t\right)\right0.06\mathrm{}\mathrm{m}\le 0.$
According to the penalty function approach, the square of the constrains are added to the objective function. Design parameters of the optimization problem are selected as follow and the upper and lower limits of them are presented in Table 1:
Since the minimum amplitude of vehicle vibration is requested for better ride comfort, the cost function given in Eq. (12) is minimized by using the PSO algorithm. This minimization is carried out by solving the five equations of motions given Eqs. (26) using fourth order RungeKutta methods. The input of the differential equations is the road excitations given in Eqs. (7, 8). The values of the parameters of the simulations are given in Table 1.
4. Optimization results
The optimization results of the design variables obtained from the simulations are given in Table 2. Furthermore, the performances of the PSO are compared with unoptimized simulations as illustrated in the Fig. 2. The unoptimized simulation is performed with the following initial values of the suspension: ${k}_{sf}={k}_{sr}=$ 35 kN/m, ${c}_{sf}={c}_{sr}={c}_{p}=$ 750 Ns/m, ${k}_{p}=$ 55 kN/m, ${k}_{pn}=$ 5500 kN/m.
It is possible to see from the figure that the maximum amplitude of the seat and vehicle body acceleration are decreased significantly by using the optimization algorithm. The maximum seat acceleration of unoptimized simulation is 6.3965, while the maximum of optimized seat acceleration is equal to 3.1751. Thus, the maximum amplitude of the seat acceleration was reduced by about 50 percent. Similarly, the vehicle body acceleration is decreased by about 73 %. In Fig. 3, the displacements of the suspension system are illustrated to visualize the optimization effect. The results show that the optimum parameters can provide small amplitudes of displacements for the seat and the vehicle body. The constrains are also in the requested range, which are equal to ${g}_{1}=$ –0.9729 m/s^{2}, ${g}_{2}=$–0.0083 m and ${g}_{3}=$ –0.0076 m.
Table 2. The optimum design variables
Design variables

Optimization results

${k}_{sf}$ [Nm]

15000

${k}_{sr}$ [Nm]

15000

${k}_{p}$ [Nm]

105000

${k}_{pn}$ [Nm]

2500000

${c}_{sf}$ [Ns/m]

1216.62

${c}_{sr}$ [Ns/m]

99.99

${c}_{p}$ [Ns/m]

170.34

Fig. 2. Optimized and unoptimized simulation results: a) seat and b) vehicle body accelerations
Fig. 3. Simulation results: a) unoptimized and b) optimized displacements by using PSO algorithm
5. Conclusions
In this paper, five degrees of freedom of half vehicle with nonlinear seat suspension model is handled to analyze and minimize the vehicle vibration for better ride comfort. A multiobjective function approach is used by using PSO algorithm. The design variables of the optimization problems are the stiffness and damping coefficients of the front and rear suspension, the nonlinear and linear stiffness of the seat suspension and the damping coefficients of the seat. The fitness function is proposed as the sum of the minimum seat and vehicle body acceleration with constraints of optimization problems. The results of this paper show that the vibration amplitude can be decreased significantly thanks to the optimum design variables of the vehicle suspension by using PSO optimization technique. In future work, it is planning to validate these results by experiments and also to couple with different controller for the active suspension system.
References
 Gillespie T. D. Fundamentals of Vehicle Dynamics. SAE International, Warrendale, 1996. [CrossRef]
 Jazar R. N. Vehicle Dynamics: Theory and Application. Springer International, New York, 2017. [Publisher]
 Cao D., Song X., Ahmadian M. Editors’ perspectives: road vehicle suspension design, dynamics, and control. Vehicle System Dynamics, Vol. 49, Issues 12, 2011, p. 328. [Publisher]
 Alkhatiba R., Nakhaie Jazarb G., Golnaraghi M. F. Optimal design of passive linear suspension using genetic algorithm. Journal of Sound and Vibration, Vol. 275, Issues 35, 2004, p. 665691. [Publisher]
 Kuznetsov A., Mammadov M., Sultan I., Hajilarov E. Optimization of a quartercar suspension model coupled with the driver biomechanical effects. Journal of Sound and Vibration, Vol. 330, Issue 12, 2011, p. 29372946. [Publisher]
 Gundogdu O. Optimal seat and suspension design for a quarter car with driver model using genetic algorithms, International Journal of Industrial Ergonomics, Vol. 37, Issue 4, 2007, p. 327332. [Publisher]
 Nariman Zadeh Salehpour N.M., Jamali A., Haghgoo E. Pareto optimization of a fivedegree of freedom vehicle vibration model using a multiobjective uniformdiversity genetic algorithm (MUGA). Engineering Applications of Artificial Intelligence, Vol. 23, Issue 4, 2010, p. 543551. [Publisher]
 Guclu R. Fuzzy logic control of seat vibrations of a nonlinear full vehicle model. Nonlinear Dynamics, Vol. 40, 2005, p. 2134. [Publisher]
 Ning D., Sun S., Li H., Du H., Li W. Active control of an innovative seat suspension system with acceleration measurement based friction estimation. Journal of Sound and Vibration, Vol. 384, Issue 11, 2016, p. 2844. [Publisher]
 Eberhart R. C., Kennedy J. A new optimizer using particle swarm theory. Proceedings of the 6th International Symposium on Micro Machine and Human Science, Nagoya, Japan, 1995. [Publisher]
 Yildiz A. R. A novel particle swarm optimization approach for product design and manufacturing. The International Journal of Advanced Manufacturing Technology, Vol. 40, Issues 56, 2009, p. 617628. [Publisher]