PSO optimum control strategy of 7 degrees of freedom semi-active suspensions

Jubin Qiao1 , Yeonsung Choi2 , Fuqiang Yang3

1, 2, 3School of Electrical and Information Engineering, Kunsan National University, Gunsansi, Jeollabukdo, Korea

1School of Mechanical and Electrical Engineering, Binzhou University, Binzhou, 256600, China

3School of Data and Computer Science, Shandong Women’s University, 250300, Jinan, China

2Corresponding author

Journal of Mechanical Engineering, Automation and Control Systems, Vol. 2, Issue 2, 2021, p. 98-108.
Received 31 July 2021; received in revised form 3 November 2021; accepted 20 November 2021; published 15 December 2021

Copyright © 2021 Jubin Qiao, 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.
Creative Commons License
Editor's Pick
Table of Contents Download PDF References
Cite this article
Views 20
Reads 13
Downloads 51
CrossRef Citations 0

The performance of vehicle body vibration and ride comfort of active or semi-active suspension with proper control is better than that with passive suspension. It is important to use simple, reliable, effective and low-cost optimal control methods to control the vehicle suspension. An important issue in the optimal control of vehicle semi-active suspension is to determine the weighting coefficient reasonably. This paper established a whole-car model of semi-active suspension systems with 7 degrees of freedom in Matlab, built optimum control system with system function, and then optimized the weight of control system by Particle Swarm Optimization (PSO) [1]. The results show that under different road input, the seven-degree-of-freedom control model for the whole vehicle controlled by particle swarm optimization algorithm can obtain better control effect, and effectively improve the comprehensive performance of the semi-active suspension system, both the vehicle's ride comfort and handling stability.

Keywords: semi-active suspension, particle swarm optimization, optimum control strategy.

1. Introduction

The suspension system is an important part of the vehicle, and its performance plays a decisive role in the ride comfort, operation, and stability. The passive vehicle suspension system was firstly presented by Olley in 1930s, which was generally composed of the stiffness-damping system. However, the parameters of passive vehicle suspension system could not be adjusted after they were determined. Normally, it cannot gain the satisfying performance in all of those three indexes or must rely on a complex structure [2], [3] Semi-active suspension can adjust suspension parameters according to the control strategy to achieve satisfying performance in different condition of the vehicle and the road [4]–[6]. And compared with active suspension, semi-active suspension has a simple structure and low cost advantages, which has been widely popularized in the market.

In practical vehicle suspension systems, the parameter uncertainties are usually unavoidable. To solve the problem, researchers have studied different modern control strategies. In addition to the traditional PID control strategy, they also adopt control strategies such as SDC (Skyhook Damper Control) [7], LQG (Linear Quadratic Gaussian) control [8], [9], fuzzy control [10], [11], SMC (Sliding Mode Control) [12], H∞ control [13], [14], adaptive control [15], [16], predictive control [17], [18], neural network control [19], [20], nonlinear intelligent control [21], and compound control [22], to reduce the vibration amplitude of vehicle body and improve the vehicle ride smoothness, stability, and ride comfort.

In the previous studies of active suspension control, PID control is often used as it has a simple principle and strong adaptability. However, it is not very effective in suspension systems with uncertain parameters. The above-mentioned traditional classic control methods have the characteristics of simple structure, reliable performance, and they are easy to implement in industry. However, under certain conditions, the system will have some problems such as delay and jitter, which will cause the system to be unstable. The system adopting the neural network control method needs to provide many labeled samples and consumes a lot of computing power. Moreover, these control algorithms are only applied to a quarter of the car model, and are not used in the simulation application of the 7-degree-of-freedom vehicle control model.

To solve the problem of the ride comfort and stability of the active suspension control system during the learning process, the passive suspension is still retained in the design of the suspension system. In order to make full use of the advantages of passive and active suspension force control and its practical application, a force control strategy for vehicle semi-active suspension system based on PSO optimization is proposed. The focuses of semi-active vehicle research are the development of adjustable damping shock absorbers, the study of control theory and methods, and the computer realization of operating systems. The control methods and controller design based on modern control theory are the cores of the semi-active suspension system research. The control strategy directly determines the performance of semi-active suspension. In the design of vehicle semi-active suspension, the optimal controller with different objective functions can be determined according to different performance requirements. Therefore, it’s important to determine the weighting coefficient reasonably for optimal control of vehicle semi-active suspension, because the selection of weighting coefficient in objective function decides the control effect of the control system itself.

The PSO algorithm is a new optimized computing technology derived from the study of bird predation behavior. It belongs to the category of stochastic global optimization technology, which is fast, effective and robust. In order to ensure the accuracy of the analysis results, this paper takes the automobile as the research object, establishes a seven-degree-of-freedom vehicle dynamics model in MATLAB, and designs the particle swarm optimization controller and different working conditions according to different working conditions. Under the vibration response of driving, the particle swarm algorithm is used to calculate the weight coefficient of each evaluation index in the optimal controller under different working conditions to meet the requirements of ride comfort and handling stability under different road grades.

2. Establishment of suspension system model

2.1. Pavement model establishment

Based on the known power spectral density, the method of filtering white noise [23] is used to construct the time domain model of the road surface as follows:

z ˙ g t = - 2 π f 0 z g t + 2 π G q n 0 u W t .

In the formula, zgt is pavement displacement, f0 is lower cut-off frequency, Gqn0 is pavement roughness coefficient, u is vehicle speed, Wt is filtered white noise. The generated model of filtered white noise random road spectrum established by Simulink tool kit is shown in Fig. 1.

Fig. 1. Random filtering white noise random rode spectrum generation model

Random filtering white noise random rode spectrum generation model

In Matlab, the white noise power Wh is 20 dB [24], vehicle speed is 60 km/h, using three pavement grades of A, B and C, the road roughness coefficient Gqn0 is 1.6×10-5 m3/cycle, 6.4×10-5 m3/cycle, 2.56×10-4 m3/cycle, and the lower cut-off frequency is 0.068 Hz. The solver is fixed-step ode4 (Runge-Kutta) solver [25], and generated A, B and C grade road spectrum are shown in Fig. 2.

Fig. 2. Integrated white noise method generated grade A, Band C road spectrum

Integrated white noise method generated grade A, Band C road spectrum

2.2. Seven degrees of freedom vehicle model establishment

In order to simplify the dynamics modeling process, following the reasonable assumptions, a dynamic model of a semi-active vehicle suspension system with 7 degrees of freedom for a certain car is established as shown in Fig. 3 [26].

(1) With the exception of tires, suspension springs and shock absorber, the rest of the components are simplified to rigid bodies and no structural deformation occurs.

(2) The suspension system consists of a linear spring and a viscous damper. The wheel is simplified as a linear spring and it is always in contact with the ground.

(3) The body pitch angle is small, and the sprung mass center of mass is vertically moved relative to the road surface excitation.

(4) Ignore the coupling effects of the front and rear axles of the vehicle, that is, the suspension mass distribution coefficient is 1.

(5) Ignore the effects of non-linear factors such as tire vertical stiffness and suspension stiffness.

Through the analysis of the motion of the car, the dynamic equation of the system is as follows:

m 1 z ¨ 1 = K f l 1 q 1 - z 1 + K f l 2 z f l - z 1 - U a 1 , m 2 z ¨ 2 = K r l 1 q 2 - z 2 + K r l 2 z r l - z 2 - U a 2 , m 3 z ¨ 3 = K f r 1 q 3 - z 3 + K f r 2 z f r - z 3 - U a 3 , m 4 z ¨ 4 = K r r 1 q 4 - z 4 + K r r 2 z r r - z 4 - U a 4 , m b z ¨ b = K f l 2 ( z 1 - z f l ) + K r l 2 ( z 2 - z r l ) + K f r 2 ( z 3 - z f r )             + K r r 2 z 4 - z r r + U a 1 + U a 2 + U a 3 + U a 4 , I p θ ¨ b = [ K r l 2 ( z 2 - z r l ) + U a 2 + K r r 2 ( z 4 - z r r ) + U a 4 ] × b             - K f l 2 z 1 - z f l + U a 1 + K f r 2 z 3 - z f r + U a 3 × a , J r φ ¨ b = [ K f l 2 ( z 1 - z f l ) + U a 1 + K r l 2 ( z 2 - z r l ) + U a 2 ] × d             - K f r 2 z 3 - z f r + U a 3 + K r r 2 z 4 - z r r + U a 4 × c ,

where Ua1 is the force of the left front wheel actuator, Ua2 is the force of the left rear wheel actuator, Ua3 is the force of the right front wheel actuator, Ua4 is the force of the right rear wheel actuator.

Among them, zfl=-zb+aθb+dφ,zrl=-zb-bθb+dφ,zfr=-zb+aθb-cφ,zrr=-zb-bθb-cφ.

In the formula, m1 is the left front wheel quality, m2 is the left rear wheel mass, m3 is the right front wheel quality, m4 is the right rear wheel mass, mb is the body quality, Ip is the body pitch moment of inertia, Jr is the body roll moment of inertia, Kfl1 is the left front wheel rigidity, Krl1 is the left rear wheel tire rigidity, Kfr1 is the right front wheel stiffness, Krr1 the right rear wheel rigidity, Kfl2 the left front wheel suspension stiffness, Krl2 is the left rear wheel suspension stiffness, Kfr2 is the right front wheel suspension stiffness, Krr2 is the right rear wheel suspension stiffness, Cfl is the left front wheel suspension damping, Crl is the left rear wheel suspension damping, Cfr is the right front wheel suspension damping, Crr is the right rear wheel suspension damping, a is the front axle to centroid distance, b is the rear axle to centroid distance, c is the right rear wheel center to centroid distance, d is the left rear wheel center to centroid distance.

For the seven-degrees of freedom semi-active suspension dynamic model of the car, taking the system’s state variables X=[z˙1,z˙2,z˙3,z˙4,z˙b,θ˙b,φ˙b,z1,z2,z3,z4,zb,θb,φb], control vectors Uc=[Ua1,Ua2,Ua3,Ua4], road disturbances W=[zg1,zg2,zg3,zg4], U=[W;Uc], output variables Y=[z˙1,z˙2,z˙3,z˙4,z˙b,θ˙b,φ˙b,z1,z2,z3,z4,zb,θb,φb], and establishing the system’s state equation: X˙=AlqrX+BlqrU,Y=ClqrX+DlqrU.

Fig. 3. Seven degrees of freedom vehicle model

Seven degrees of freedom vehicle model

3. Designing of optimal particle swarm optimization controller

3.1. Particle swarm optimization

The PSO algorithm can be understood as an iterative optimization tool. The system initializes a set of random solutions and iterative searches for the optimal value. In the PSO algorithm, each particle’s position Xj represents a candidate solution to the optimization problem. The particles fly at a certain speed in the search space. This speed will dynamically adjust its flight direction and speed based on its own flight experience and the flight experience of its companions. All particles have a fitness value determined by the fitness function according to their current position (current position Xj) and the best position found so far (the individual extremum Pj). The particle searches in the solution space according to the current optimal particle, and finally finds the optimal position.

In the particle swarm optimization algorithm, it is firstly necessary to initialize the current position of the particle swarm and determine the objective function of the particle swarm optimization control. Then it needs to update the position, velocity, inertia weight and acceleration factor of the particle swarm according to certain rules, and update the iteration by recalculating the function of the target value particle and comparing the size of the particle [27]. Each variable in the particle is represented by a real number. The iterative steps of particle swarm optimization algorithm are as follows:

(1) There are n particles in the D-dimensional search space. During random initialization, the position and velocity of the j-th particle are represented by Xi=(X1,j,X2,j,Xd,j) and Vi=(V1,j,V2,j,Vd,j) respectively;

(2) After evaluating the fitness of each particle and calculating the fitness value of the objective function, then store the current particle position and fitness value in each particle Pi=(P1,j,P2,j,Pd,j).

(3) Then calculate the individual extreme value. For each particle, compare the fitness value of the current position Xi with the best position of the fitness of its own experience Pi, if it is better than Pl, its position will replace Pl;

(4) For the global extremum calculation, the fitness value of the local optimal position Pl is compared with the fitness value of the best position Pg=(P1,j,P2,j,,Pd,j) experienced in the group. If it is better than Pg, the position will be replaced;

(5) Particle parameter update.

Update the inertia weight coefficient: w(jstep)=wmin+wmax-wminiter(iter-jstep), where, jstep – current number of iterations, wmax – maximum inertia factor, iter – maximum evolutionary number, wmin – minimum inertia factor.

Update accelerating factor: c1(jstep)=(c1e-c1j)jstepiter+c1j, c2(jstep)=(c2e-c2j)jstepiter+c2j, where, c1j– the initial value of c1, c2j – the initial value of c2, c1e– the final value of c1, and c2e – the final value of c2.

(6) Particle speed update.

Update the speed of particles:

V i , j j s t e p + 1 = ω V i , j j s e p + c 1 r 1 P i , l - X i , j j s t e p + c 2 r 2 P i , g - V i , j j s t e p ,
i = 1,2 , , d ,           j = 1,2 , , n .

Determine if the new speed value is within the range: Vi,j(jstep+1)Vmax, i=1,2,,d, j=1,2,,n.

If vi,j is greater than vmax , use vmax  instead of vi,j;

(7) Particle position update.

Substituting the updated velocity value into the formula:

X i , j j s t e p + 1 = X i , j j s t e p + V i , j j s t e p + 1 ,         i = 1,2 , , d ,           j = 1,2 , , n ,

the new particle position value xi,j is obtained;

(8) Judge whether the termination condition is satisfied, and if it is satisfied, the process goes to step (8); otherwise, it goes to step (2) to continue the search calculation;

(9) Stop searching when the stop condition is met, output Pg, end the seek operation.

3.2. Optimal controller design

The goal of the optimal control of semi-active suspension is to make the car get higher smoothness and control stability. If the vertical acceleration Z value of the sprung center of mass is smaller, it means that the ride comfort of the vehicle is better. The evaluation index of vehicle handling stability mainly depends on sprung mass pitch angular acceleration, sprung mass lateral acceleration, suspension dynamic deflection, four-wheel dynamic load, etc. The smaller their value, the better the performance. Taking into account the complexity of the model and a large number of observations, the linear quadratic function of the semi-active suspension is defined as follows [28]:

J = P 1 x ¨ b 2 + P 2 θ ¨ b 2 + P 3 φ ¨ b 2 + P 4 ( z f l - z 1 ) 2 + P 5 ( z r l - z 2 ) 2 + P 6 ( z f r - z 2 ) 2             + P 7 ( z r r - z 4 ) 2 + P 8 ( z 1 - z g 1 ) 2 + P 9 ( z 2 - z g 2 ) 2 + P 10 ( z 3 - z g 3 ) 2 + P 11 ( z 4 - z g 4 ) 2 ,

where P1, P2, P3, P4, P5, P6, P7, P8, P9, P10, P11 are respectively the body acceleration, the pitch angle acceleration, the lateral inclination acceleration, the left front suspension dynamic stroke, the left rear wheel suspension dynamic stroke, the right front suspension dynamic stroke, the right rear wheel suspension dynamic stroke, the left front wheel tire change, left rear wheel tires deformation, right front tires deformation, the weighting coefficient of the deformation of the rear wheel tire, and written in the form of a matrix:

J = X T Q X + U T R U + 2 X T N U ,

where Q and R respectively represent the weighted matrices of the state variables and the control variables. The different values of Q and R allow the different weighting coefficients to be added to the different performance indexes in the function. For different performance requirements, a different optimal control feedback gain matrix can be obtained by changing the weighted coefficient matrix of the state variable and the control variable. And the controller can be designed to meet the requirements of different performance. In order to control the weighting coefficient matrix better, this paper proposes a particle swarm optimization algorithm to optimize the weighting coefficient of the optimal controller.

3.3. Designing of optimal controller for particle swarm

The semi-active suspension simulation model of Particle swarm optimal control is constructed as shown in Fig. 4.

Fig. 4. The semi-active suspension simulation model of Particle swarm optimal control

The semi-active suspension simulation model of Particle swarm optimal control

The operation block diagram of the optimal controller of particle swarm is shown in Fig. 5.

Fig. 5. Simulation model of PSO-based best control strategy of semi-active suspensions

Simulation model of PSO-based best control strategy of semi-active suspensions

4. The optimal result subgroup of the optimal controller

According to the evaluation indexes corresponding to the vehicles on different roads, the fitness function corresponding to different road conditions is established, and the particle swarm algorithm is used to determine the weighting coefficient of each evaluation index. The relevant parameters of the vehicle model are shown in Table 1. These parameters, together with the parameters associated with the road, are programmed with Matlab and stored in the initial_parameters.m folder.

Table 1. Model parameters of vehicle suspension system.

The parameter
The numerical
The parameter
The numerical
m 1 / kg
C f l / (kN / ms-1)
m 2 / kg
C r l / (kN / ms-1)
m 3 / kg
C f r / (kN / ms-1)
m 4 / kg
C r r / (kN / ms-1)
m b / kg
K f l 1 / ( kN.m-1)
K f l 2 / (kN.m-1)
K r l 1 / ( kN.m-1)
K r l 2 / (kN.m-1)
K f r 1 / ( kN.m-1)
K f r 2 / (kN.m-1)
K r r 1 / ( kN.m-1)
K r r 2 / (kN.m-1)
d / m
I p / (kgm2)
c / m
J r / (kgm2)
a / m
b / m

The test vehicle traverses A-C class roads at a speed of 60 km/h. Set the solver step size to ode4 (Runge-Kutta), and then run the program. The following Table 2 shows the optimized controller parameters.

Table 2. Weighting coefficients under different levels of pavement

The weighted coefficient
The road level
Class A pavement
Class B pavement
Class C pavement
P 1
P 2
P 3
P 4
P 5
P 6
P 7
P 8
P 9
P 10
P 11

5. Conclusions

In order to verify the accuracy of the vehicle model and the rationality and effectiveness of the control strategy, the semi-active suspension system was simulated and analyzed under different road conditions on the loop simulation platform. The test vehicle traverses A, B, and C roads at a constant speed of 60 km/h. The solver is seated to a fixed step length ode15 s (stiff/NDF) [29], and the simulation work step length is 0.02 s. The entire simulation time is 60 s, of which 0-20 s to A-class road, B-class road 20-40 s, and C-class road 40-60 s [30]. Through the particle swarm optimization algorithm of the semi-active suspension system, the body acceleration, pitch acceleration and roll acceleration are obtained. Their simulation results are shown in Figs. 6-7.

Fig. 6. The characteristics of the body vertical acceleration in time domain

The characteristics of the body vertical acceleration in time domain

Fig. 7. The characteristics of the pitch acceleration in time domain

The characteristics of the pitch acceleration in time domain

The evaluation index of car ride comfort is mainly the three axial weighted accelerations from the seat support points. More importantly, the human body is mostly affected by vertical vibration caused by uneven road surfaces. Therefore, in vehicle control, the vertical vibration at the center of mass of the vehicle body is extremely important. It can be seen in Fig. 6, compared with passive suspension, based on particle swarm optimal control of semi-active suspension system of vehicle body acceleration value in three grade pavement are much smaller than the passive suspension system. The pitch and roll movement of the vehicle also affect the comfort of the human body. It can be seen in Fig. 7 that the pitch acceleration has been controlled to a certain extent, and the control effect is still obvious. With the deterioration of traffic conditions, the pitch acceleration value of the C-level road is still controlled within a certain range. It can be seen from Fig. 8 that the lateral inclination acceleration has little effect on the change of the road slope, the acceleration values are all well attenuated, and the control effect is the most ideal. The simulation results show that the semi-active suspension controlled by the particle swarm optimization algorithm has a good control effect in the above three important degrees of freedom compared with the passive suspension. That means that the semi-active suspension controlled by the particle swarm optimization algorithm has a better attenuation effect, reducing the impact of the human body in the vertical vibration and roll and pitch movement of the vehicle body.

Fig. 8. Characteristic diagram of lateral inclination acceleration in time domain

Characteristic diagram of lateral inclination acceleration in time domain

The simulation results of vehicle suspension system domain under random time road show that semi-active suspension based on particle swarm optimization control can effectively reduce the vertical vibration of vehicle acceleration on different levels of road surface. This shows that the particle swarm optimal controller can adjust the posture of the car body according to changes in traffic conditions. When the road conditions are poor, it can ensure driving safety while improving comfort and smoothness.


  1. J. Kennedy and R. C. Eberhart, “Particle swarm optimization,” in Proceedings of IEEE International Conference on Neural Networks, pp. 760–766, 2011. [Search CrossRef]
  2. M. Z. Q. Chen, K. Wang, Y. Zou, and G. Chen, “Realization of three-port spring networks with inerter for effective mechanical control,” IEEE Transactions on Automatic Control, Vol. 60, No. 10, pp. 2722–2727, Oct. 2015, [Publisher]
  3. K. Wang, M. Z. Q. Chen, and Y. Hu, “Synthesis of biquadratic impedances with at most four passive elements,” Journal of the Franklin Institute, Vol. 351, No. 3, pp. 1251–1267, 2014. [Publisher]
  4. X. Ma, P. K. Wong, J. Zhao, J.-H. Zhong, H. Ying, and X. Xu, “Design and testing of a nonlinear model predictive controller for ride height control of automotive semi-Active air suspension systems,” IEEE Access, Vol. 6, pp. 63777–63793, 2018, [Publisher]
  5. X. Tang, H. Du, S. Sun, D. Ning, Z. Xing, and W. Li, “Takagi-Sugeno fuzzy control for semi-active vehicle suspension with a magnetorheological damper and experimental validation,” IEEE/ASME Transactions on Mechatronics, Vol. 22, No. 1, pp. 291–300, Feb. 2017, [Publisher]
  6. M. Z. Q. Chen, Y. Hu, C. Li, and G. Chen, “Application of semi-Active inerter in semi-active suspensions via force tracking,” Journal of Vibration and Acoustics, Vol. 138, No. 4, pp. 1–11, Aug. 2016, [Publisher]
  7. J. L. Yao, S. Taheri, and S. M. Tian, “A novel semi-active suspension design based on decoupling skyhook control,” Journal of Vibroengineering, Vol. 16, No. 3, pp. 1318–1325, 2014. [Search CrossRef]
  8. B. Lan and F. Yu, “Design and simulation analysis of LQG controller of active suspension,” Transactions of the Chinese Society of Agricultural Machinery, Vol. 35, No. 1, pp. 13–18, 2004. [Search CrossRef]
  9. Y. W. Yu, C. C. Zhou, and L. L. Zhao, “Design of LQG controller for vehicle active suspension system based on alternate iteration,” Journal of Shandong University, Vol. 47, No. 4, pp. 50–58, 2017. [Search CrossRef]
  10. R. W. Wang, X. P. Meng, and D. H. Shi, “Fuzzy control of vehicle ISD semi-active suspension,” Transactions of the Chinese Society of Agricultural Machinery, Vol. 44, No. 12, pp. 1–5, 2013. [Search CrossRef]
  11. D. Y. Wang and H. Wang, “Control method of vehicle semi active suspensions based on variable universe fuzzy control,” Chinese Journal of Mechanical Engineering, Vol. 28, No. 3, pp. 366–372, 2017. [Search CrossRef]
  12. J.-L. Yao, W.-K. Shi, J.-Q. Zheng, and H.-P. Zhou, “Development of a sliding mode controller for semi-active vehicle suspensions,” Journal of Vibration and Control, Vol. 19, No. 8, pp. 1152–1160, Jun. 2013, [Publisher]
  13. H. P. Du and N. Zhang, “Constrained H∞ control of active suspension for a half-car model with a. time delay in control,” Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, Vol. 222, No. 5, pp. 665–684, 2008. [Publisher]
  14. W. Sun, H. Gao, and O. Kaynak, “Finite frequency H∞ control for vehicle active suspension systems,” IEEE Transactions on Control Systems Technology, Vol. 19, No. 2, pp. 416–422, Mar. 2011, [Publisher]
  15. W. Sun, Z. Zhao, and H. Gao, “Saturated adaptive robust control for active suspension systems,” IEEE Transactions on Industrial Electronics, Vol. 60, No. 9, pp. 3889–3896, Sep. 2013, [Publisher]
  16. F. Zhao, S. S. Ge, F. Tu, Y. Qin, and M. Dong, “Adaptive neural network control for active suspension system with actuator saturation,” IET Control Theory and Applications, Vol. 10, No. 14, pp. 1696–1705, Sep. 2016, [Publisher]
  17. M. Canale, M. Milanese, and C. Novara, “Semi-active suspension control using “fast” model-predictive techniques,” IEEE Transactions on Control Systems Technology, Vol. 14, No. 6, pp. 1034–1046, Nov. 2006, [Publisher]
  18. D. Wang, D. Zhao, M. Gong, and B. Yang, “Research on robust model predictive control for electro-hydraulic servo active suspension systems,” IEEE Access, Vol. 6, pp. 3231–3240, 2018, [Publisher]
  19. S. H. Zhu, B. Z. Lv, and H. Wang, “Neural network control method of automotive semi-active air suspension,” Journal of Traffic and Transportation Engineering, Vol. 6, No. 4, pp. 66–70, 2006. [Search CrossRef]
  20. H. Li, Y. H. Feng, and L. Su, “Vehicle active suspension vibration control based on robust neural network,” Chinese Journal of Construction Machinery, Vol. 15, No. 4, pp. 324–337, 2017. [Search CrossRef]
  21. Y. Qin, J. J. Rath, C. Hu, C. Sentouh, and R. Wang, “Adaptive nonlinear active suspension control based on a robust road classifier with a modified super-twisting algorithm,” Nonlinear Dynamics, Vol. 97, No. 4, pp. 2425–2442, Sep. 2019, [Publisher]
  22. S. K. Sharma, U. Saini, and A. Kumar, “Semi-active control to reduce lateral vibration of passenger rail vehicle using disturbance rejection and continuous state damper controllers,” Journal of Vibration Engineering and Technologies, Vol. 7, No. 2, pp. 117–129, Apr. 2019, [Publisher]
  23. Fan Yu, “Vehicle Dynamics and Control,” Ph.D. Thesis, China Machine Press, Beijing, 2010. [Search CrossRef]
  24. Zhang Yufen, Long Jinlian, Li Jing, and Lu Jiaxuan, “Research on LQG of active suspension based on immune particle swarm optimization.,” Computer Engineering and Applications, Vol. 54, No. 6, pp. 252–256, 2018. [Search CrossRef]
  25. O. P. Agrawal and A. A. Shabana, “Application of deformable-body mean axis to flexible multibody system dynamics,” Computer Methods in Applied Mechanics and Engineering, Vol. 56, No. 2, pp. 217–245, Jun. 1986, [Publisher]
  26. Jun Yao and Songhui Ma, “Modeling and Simulation of Simulink,” Xidian University Press, Xian, 2009. [Search CrossRef]
  27. Andries Engelbrecht P., Fundamentals of Computational Swarm Intelligence. John Wiley and Sons, 2011. [Search CrossRef]
  28. X. Song and M. Ahmadian, “Characterization of semi-active control system dynamics with magneto-rheological suspensions,” Journal of Vibration and Control, Vol. 16, No. 10, pp. 1439–1463, Sep. 2010, [Publisher]
  29. Zhihong Guo, Konghui Guo, and Xiaolin Song, “Study of semi-active suspension control strategy based on an identification model,” Journal of Hunan University (Natural Sciences), Vol. 37, No. 12, 2010. [Search CrossRef]
  30. Xiaopeng Wang, Jianjun Liu, and Long Wu, “A simulation research on seven DOF semi-active full vehicle suspension,” (in Chinese), Journal of Hunan University of Technology, Vol. 30, No. 6, 2016. [Search CrossRef]