Abstract

This paper introduces a novel method to control the operation of the active suspension system. In this research, a quarter-dynamic model is used to simulate the vehicle’s vibrations. Besides, the sliding mode-PID-integrated algorithm with five state variables is proposed to be used. This is a completely original and novel algorithm. The process of establishing the control algorithm is clearly described. The simulation is performed by the MATLAB software. The results of the paper have shown the advantages of the sliding mode-PID algorithm used in this research. Accordingly, the displacement and acceleration of the sprung mass were significantly reduced when this algorithm was used. The maximum and average values of the displacement of the sprung mass are only 1.31% and 1.29%, respectively, compared with the situation of the vehicle using a passive suspension. Similarly, the value of the acceleration is 6.98% and 2.94%, respectively. In addition, the phenomenon of “chattering” has also been significantly reduced when using this controller. In the future, some more complex algorithms can be proposed.

1. Introduction

The vehicle’s vibrations are generated by stimuli from the road surface. These vibrations can affect passengers and cargo on the vehicle. The suspension is used to quell the vehicle’s vibrations quickly. The suspension system divides the vehicle into two separate parts, including the sprung mass ms and the un-sprung mass mu [1]. According to [2], the entire vehicle compartment and the sub-assemblies located above the suspension system are called the sprung mass. The remainder is called the un-sprung mass. The excitations will be transmitted from the pavement to the un-sprung mass. It will then continue to be transferred onto the sprung mass through the suspension system. Therefore, the vehicle’s vibration can be reduced thanks to the suspension system, which is equipped with the vehicle.

The conventional suspension system used on a vehicle is called a mechanical suspension system (passive suspension system). According to Kashem, a passive suspension system consists of a spring, a shock absorber, and a lever arm [3]. Springs are used to regulate the oscillation. Besides, dampers help to quell the vibrations of the vehicle. For a passive suspension system, the stiffness of the springs and dampers cannot be changed. As a result, the vehicle’s smoothness may be affected when the vehicle is on the road. To overcome this situation, several solutions have been proposed to change the stiffness of the suspension system. In [4], Anaya-Martinez et al. suggested the use of a semi-active suspension system. This system uses electromagnetic dampers to replace conventional dampers. According to Sha et al., when an electric current is applied to the shock absorber, the metal particles inside will be tightly arranged. Therefore, the circulation rate of the liquid flow inside will change. This means that the damping stiffness will also change [5]. In addition, changing the spring’s stiffness has also been used on some vehicles today [6]. In [7], Sun et al. introduced the air suspension system, which can adjust the stiffness of the air spring. According to Nguyen, the stiffness of the air spring is adjusted based on the air pressure applied to the inside of the balloon. This change is continuous and fully automatic. As a result, the vehicle’s smoothness can be further improved [8]. Another method that has been proposed by Fazeli et al. is to use an active suspension system [9]. The active suspension system has a hydraulic actuator that is placed between the sprung mass and the un-sprung mass. According to Lee et al., the hydraulic actuator operates based on the opening and closing of the servo valves in the system. The operation of the servo valve depends on the supplied voltage signal [10]. The active suspension system can generate an impact force on both the sprung mass and the un-sprung mass. Therefore, the vibration of the vehicle can be improved most effectively.

As for the active suspension system, its most important issue is the control algorithms. The PID controller for active suspension was introduced by Moaaz and Ghazaly [11]. According to them, this controller is determined based on three parameters, including KP, KI, and KD. These parameters can be determined by the self-tuning method. In [12], this method was pointed out by Emam. This is a simple method; however, its accuracy is still not good. Besides, a fuzzy control algorithm is also used to adjust the parameters of the PID controller. This was shown in the study by Ding et al. [13]. Similarly, this method is also mentioned in the paper by Mahmoodabadi and Nejadkourki [14]. The PID controller is only used to control one object. This is the SISO method. Normally, only the value of displacement or acceleration of the sprung mass is considered for this algorithm. It is suitable for a quarter-dynamic model. If more complex models are considered, such as half models or spatial models, MIMO control methods should be used. According to Anh, the LQR controller can control many parameters of the control model [15]. This algorithm is established based on the matrix of the control model [16]. The choice of the coefficients of the matrix is very important. To limit the noise of the system, a Gaussian filter is integrated into this controller [17]. Therefore, the name “LQG controller” was born [18]. The performance of this controller is quite high [19, 20]. However, this is still only a linear control method. In many situations, the system’s response is still not stable. Therefore, nonlinear control methods are proposed for using.

Many nonlinear control algorithms have been used for the active suspension system. In [21], Wang et al. introduced the concept of the SMC. According to Chu et al., the SMC helps the active suspension system operate more stably. From there, the phenomenon of rollover can be improved [22]. According to Nguyen, the SMC needs to establish a suitable sliding surface [23]. Controlled objects will move along the sliding surface to move toward a stable position [24]. If the hydraulic actuator is considered in this model, it needs to be linearized appropriately. This has been shown in the study of Bai and Guo [25]. However, these authors only give linear equations without explanation. Therefore, the specific demonstration of the linearization of the actuator for the SM algorithm has been shown in Nguyen’s paper [26]. The input signals of the controller can be obtained through a sensor or a camera that is built into the vehicle [27]. According to Azizi and Mobki, the “chattering” phenomenon still exists when using this controller [28]. This phenomenon can be improved by using a hybrid fuzzy SMC; this has been demonstrated by Shaer et al. in [29]. Because the SM control algorithm is a nonlinear method, therefore, the controller design process is extremely complicated [3033]. Overall, this controller provides stable performance in a wide range of oscillating situations [3437]. In addition, robust and adaptive control algorithms are also used for the active suspension system. In [38], Yao et al. used the H∞ robust algorithm to improve the efficiency of the suspension system to reduce the roll angle of the vehicle. This algorithm is determined based on the matrix of oscillating states [2]. The system always ensures stability when using this algorithm [3942]. In some situations, the vehicle’s suspension system needs to oscillate continuously to respond well to changes in external stimuli. The adaptive control method is suitable for these situations [43, 44]. Besides, a few nonlinear control algorithms have also been used. Their efficiency is very high [4547].

In recent years, intelligent control methods have also been studied and widely applied to automotive mechatronic systems. In [48], Haddar et al. used the i-PD algorithm for the active suspension system. In their paper, the MFC method is used. Besides, the PSO algorithms with high performance are also used for the half model to optimize the vehicle’s vibration [49]. In many situations, it is difficult to accurately determine the vehicle’s vibrational state. Therefore, the use of the fuzzy control algorithm is completely appropriate [50]. According to Khan et al., the membership functions and fuzzy control rules are built based on the experience of the authors [51]. An intermediate state can exist using this algorithm [5255]. In addition, the use of ANN control algorithms also brings high efficiency to the system [56, 57]. The fundamental oscillation states are stored with the concept of “big data.” Based on this data, the controller can self-learn and make judgments in random situations [5860].

Based on the statements introduced above, this paper focuses on active suspension system control to improve the vehicle’s ride comfort and stability. Most of the previous studies used only linear control algorithms to control the active suspension system. Some studies have used the SMC nonlinear control algorithm. However, the “chattering” phenomenon still occurs when using this algorithm. Therefore, this research proposes a method using two integrated control algorithms SMC-PID to limit the “chattering” phenomenon. Besides, the hydraulic actuator is also linearized to establish the five variable models. This is considered the new contribution of the paper. The controller’s states that are established in this paper are assumed to be available. For conditions that are not available, the interval observation technique may be used. This is a unique technique [6163]. The content of this paper consists of four sections. Suspension system issues and literature review are introduced in the first section. In the second section, the vehicle dynamics model and the SMC-PID control algorithm are established. Next, the simulation process will be introduced in the third section. This process is performed in the MATLAB environment. And in the last section, some suggestions are made.

2. Materials and Methods

To simplify the simulation of the vehicle’s vibration, a quarter-dynamic model is used in this paper (Figure 1). This model includes the sprung mass ms and the un-sprung mass mu. The link between them is the suspension system, which includes the linear spring and damper. To enhance ride comfort, a hydraulic actuator is fitted alongside the suspension system. Thus, the passive suspension system becomes the active suspension system.

The quarter-dynamics model has two degrees of freedom corresponding to the two directions of motion, including the vertical displacement of the sprung mass zs and the vertical displacement of the un-sprung mass zu. Using mechanical theory, the equations describing the oscillations of the vehicle are given as follows:where

K: Stiffness of the spring

C: Stiffness of the damper

KT: Stiffness of the tire.

To reduce the vehicle’s vibration, the hydraulic actuator must generate the impact force FA. This force acts on both the sprung mass and the un-sprung mass. The impact force FA is generated based on the pressure change between the two chambers inside the cylinder of the actuator. The value of this force is proportional to the cross-sectional area of the piston (Sp) and the pressure difference (ΔP).

The actuator operates on the voltage signal supplied by the controller. When voltage is applied to the hydraulic actuator, the internal valves move. This displacement changes the fluid pressure inside the cylinder. The relationship between the voltage signal u(t), servo valve displacement xsv, and pressure difference is a very complex multivariable function. It even depends on the displacement of the suspension system xs. This relationship is shown through equations

If the original equations of the hydraulic actuator are used, the design of the SM controller is very difficult. According to Nguyen, the impact force FA generated by the actuator can be linearized according to equation [26]

According to [24], the acceleration of the sprung mass can be approximated as follows:

Let state variables

Taking derivatives of state variables,

Let e1(t) be the error between the output signal and the set point:

Take the derivative of both sides of equation (9):

Take the second derivative of equation (9):

Take the derivative of both sides of (11), and this is the third derivative of the error e1(t):

Continuing to derive both sides of (12), the fourth-order differential of the error signal is obtained:

Because the oscillator has five state variables, it is necessary to take the fifth derivative of the output signal y(t):

If r(t) is the disturbance, this value can be ignored. Equation (14) is rewritten as follows:where

According to the Lyapunov theory, the controlled object will slide along the surface to return to a steady state (Figure 2). Therefore, the slip surface needs to be determined based on the stability condition of the system. The slip surface is a function that depends on the signal of error e1(t):

Once the slip surface has been determined, the output signal of the SM controller can be easily determined as equation

The output signal of the SM controller is the desired threshold of the subsequent designed PID controller (Figure 3). Let e2(t) be the error signal of the PID controller:

According to [64], the PID controller includes three stages: a proportional stage, an integral stage, and a derivative stage. For each stage, there is a corresponding coefficient that needs to be selected: KP, KI, and KD. The current signal generated by the PID controller is shown in equation

After the controller has been designed, the simulation process is performed. This process happens in the MATLAB environment with specific parameters.

3. Results and Discussion

3.1. Simulation Situations

The results of the research are obtained based on the simulation process. The vehicle parameters used for the calculation are shown in Table 1.

There are three simulation cases used in this paper. In the first two cases, sine pavement excitation is used. However, the excitation amplitudes of the oscillations are different. In the third case, the excitation of the pavement is pulse. In each case, the vehicle will be simulated through four situations, including the vehicle using the passive suspension system (none), the active suspension system controlled by a PID controller (PID), the active suspension system controlled by the SM controller (SMC), and the active suspension system controlled by the SMC-PID controller (SMC-PID). The results of the simulation process are analyzed as shown below.

3.2. Results

Case 1. In the first case, the excitation amplitude of the pavement reaches 50 (mm). The displacement of the sprung mass is shown in Figure 4. Accordingly, the maximum value of displacement can be up to 61.96 (mm) if the vehicle uses only the passive suspension system. The displacement of the sprung mass has been markedly reduced when the hydraulic actuator is equipped. Its maximum values are only 23.35 (mm), 7.73 (mm), and 0.81 (mm), respectively. Besides the maximum value, the average value should also be considered. The average value of oscillations calculated according to the RMS criteria of the four situations above reached 32.08 (mm), 12.99 (mm), 3.86 (mm), and 0.42 (mm), respectively. This result shows that the SMC-PID controller helps to reduce displacement of the sprung mass. The vehicle body is almost unchanged during the oscillation.
The change in the acceleration of the sprung mass in the first case is shown in Figure 5. For a conventional mechanical suspension system, the value of the acceleration is the maximum. Its maximum value can be up to 1.29 (m/s2). After that, this value fluctuates more steadily, with a maximum amplitude of about 0.50 (m/s2). The average value of the acceleration corresponding to this situation is also the highest, reaching 0.34 (m/s2).
If the active suspension system is used, the vehicle’s vibration can be further improved. According to the results shown in the graph of Figure 5, the maximum value of acceleration when the active suspension system is equipped is 1.04 (m/s2), 0.76 (m/s2), and 0.09 (m/s2), which corresponds to three situations, including PID, SM, and SMC-PID. The average values calculated according to the RMS in the above order are 0.14 (m/s2), 0.13 (m/s2), and 0.01 (m/s2). According to this result, the average value of acceleration when the PID controller and the SMC controller are used is equivalent. In the case of a vehicle using a PID controller, the acceleration takes the form of a pulse in the first phase. From the following phases onward, the acceleration is periodic with time. In contrast, the vehicle’s acceleration when using the SMC controller continuously fluctuates. This phenomenon is called “chattering.” When using the SMC controller, it is difficult to avoid this phenomenon. However, the SMC-PID integrated controller can solve this problem better. Both the maximum and average values of acceleration when using the SMC-PID controller are very small. In addition, the phenomenon of “chattering” has been significantly improved. It only appears for a short time when the vehicle starts to oscillate with a rather small amplitude. As a result, the vehicle’s smoothness can be enhanced.
For an active suspension system, to reduce displacement and acceleration of the sprung mass, the hydraulic actuator will operate continuously. The impact force generated by the actuator will act on both the sprung mass and the un-sprung mass. The goal of the system is to ensure the stability and comfort of the sprung mass. As a result, the un-sprung mass is likely to fluctuate more. The vibrations of the un-sprung mass influence the overall vibrations of the vehicle. However, this effect is very small. If the oscillation of the un-sprung mass is small, the vehicle’s stability and smoothness will be further improved. The changes in displacement and acceleration of the un-sprung mass are shown in Figures 6 and 7.
According to previous results, the displacement of the un-sprung mass is quite large when the vehicle uses the active suspension system [24, 26]. If the displacement of the sprung mass is smaller, the displacement of the un-sprung mass will increase. The same goes for acceleration. However, the SMC-PID controller maintains steady oscillation for the un-sprung mass. Displacement of un-sprung mass in the case of vehicles with an active suspension system controlled by the SMC-PID method is always closely tracked in the case of vehicles using a conventional passive suspension system. Even though this value is slightly smaller, in the first phase, the acceleration of the un-sprung mass is still quite large. This value rapidly decreases and changes cyclically in subsequent periods. The impact of the first stage is very small. The vibration of the whole vehicle is still guaranteed in the allowed conditions.

Case 2. In the second case, the amplitude of the pavement excitation will be increased. Therefore, the vibration of the vehicle will be greater. This can demonstrate the effectiveness of the active suspension system. In this case, the pavement bump has a magnitude of up to 100 (mm). As a result, the displacement of the sprung mass has increased for the four situations, to 123.92 (mm), 46.70 (mm), 15.60 (mm), and 1.62 (mm), respectively. The average value of displacement also changed accordingly, reaching 78.46 (mm), 31.95 (mm), 9.38 (mm), and 1.01 (mm). The difference in value between the vehicle using the passive suspension system and the vehicle using the active suspension system controlled by the SMC-PID controller can be up to 76.49 times and 77.68 times (Figure 8). This difference is huge. It helps to regulate the vibrational state of the vehicle.
The acceleration of the sprung mass has a great influence on the stability and comfort of the vehicle when traveling on the road. If the value of acceleration is too large, smoothness will no longer be guaranteed. The active suspension system effectively improves this problem. This is demonstrated in Figure 9. According to this result, the acceleration of the sprung mass starts to stabilize in the second phase corresponding to the situation of the vehicle using the passive suspension system and the active suspension system controlled by the PID or SMC-PID algorithm. As for the conventional SMC algorithm, the “chattering” phenomenon still occurs during the whole oscillation process. In particular, the effect of this phenomenon is greatest at the first stage. The maximum values of the acceleration corresponding to four simulated situations are 2.58 (m/s2), 2.08 (m/s2), 1.51 (m/s2), and 0.18 (m/s2), respectively. Similarly, their average values have been specifically calculated to be around 0.81 (m/s2), 0.34 (m/s2), 0.31 (m/s2), and 0.03 (m/s2). The obtained results show that the acceleration of the sprung mass when using the SMC-PID algorithm is only 6.98% and 3.70%, respectively, compared to the situation of the vehicle using the passive suspension system. Again, the efficiency of the control algorithm established in this paper is demonstrated in detail.
Like in the first case, the displacement of the sprung mass in the second case also tended to be closely tracked in all four situations examined (Figure 10). This difference is quite small; it is not more than 2%.
The acceleration of the un-sprung mass has also changed more than in the first case. In subsequent phases of the oscillation, the values of the accelerations are always tracked to each other (Figure 11). In general, the oscillation in the first phase when the vehicle uses the SMC-PID controller does not greatly affect the oscillation of the whole vehicle.

Case 3. For some controllers, even though pavement excitation has ceased, the controller continues to operate. This will cause vehicle instability even though the bumpy road no longer exists. According to Figure 12, the displacement of the sprung mass is greatest when the vehicle uses the passive suspension system. This value can go up to as close as 150.91 (mm). Meanwhile, the displacement of the vehicle body when the SMC-PID algorithm is used is only 2.51 (mm). For the simple SMC controller, the actuator remains active even though the excitation from the pavement has ceased. Even its oscillation is larger than the situation of using a PID controller.
The acceleration of the sprung mass also varies continuously with the excitation conditions of the pavement (Figure 13). The active suspension system controlled by the SMC-PID controller still ensures the stability and smoothness of the vehicle. However, the “chattering” phenomenon persisted for a short period.
The results of the simulation process are summarized in Tables 25.

4. Conclusions

Bumps on the road can cause vibrations in the vehicle when traveling on the road. The suspension system is used to quench these oscillations. To improve the stability of the vehicle, an active suspension system is proposed to replace the conventional passive suspension system. This suspension system has a hydraulic actuator, which generates force to act on the sprung mass and un-sprung mass. As a result, the vehicle’s smoothness and comfort can be improved. There are several methods used to control the suspension system. Linear control methods such as PID and LQR cannot provide stable performance under many complex conditions. Therefore, the nonlinear control algorithm SMC is used to control the operation of the system. However, the phenomenon of “chattering” still occurs for the SMC algorithm.

This research established a quarter-dynamics model to simulate the vibrations of the vehicle. Besides, the integrated control method SMC-PID is proposed to be used. The simulation is done in the MATLAB environment. There are three types of pavement excitation used for system vibration survey purposes. The outputs of the research include displacements and accelerations of the sprung mass and un-sprung mass. The results of the paper have shown the superiority of the SMC-PID integrated algorithm. Accordingly, the displacement and acceleration of the sprung mass were significantly reduced when this algorithm was used. Compared with other control methods mentioned in the paper, such as PID and SMC, this algorithm provides higher efficiency in all conditions. Besides, the phenomenon of “chattering” has also been significantly reduced. In the future, an experimental process can be conducted to verify the effectiveness of this controller.

Nomenclature

FA:Force of the hydraulic actuator, N
FC:Force of the damper, N
FK:Force of the spring, N
FKT:Force of the tire, N
i(t):Output signal of the PID controller
r(t):Roughness on the road, m
u(t):Output signal of the sliding mode controller
zs:Displacement of the sprung mass, m
zu:Displacement of the un-sprung mass, m

Abbreviation

ANN:Artificial neural network
i-PD:Intelligent proportional-derivative
LQG:Linear quadratic Gaussian
LQR:Linear quadratic regulator
MFC:Model-free control
MIMO:Multi-input, multi-output
PID:Proportional-integral-derivative
PSO:Particle swarm optimization
RMS:Root mean square
SISO:Single input, single output
SMC:Sliding mode controller.

Data Availability

The data used to support this research are included within this paper.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.