# Vibration Control of a Semiactive Vehicle Suspension System Based on Extended State Observer Techniques

Authors: Ze Zhang; Hamid Reza Karimi; Hai Huang; Kjell G. Robbersmyr
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

A feedback control method based on an extended state observer (ESO) method is implemented to vibration reduction in a typical semiactive suspension (SAS) system using a magnetorheological (MR) damper as actuator. By considering the dynamic equations of the SAS system and the MR damper model, an active disturbance rejection control (ADRC) is designed based on the ESO. Numerical simulation and real-time experiments are carried out with similar vibration disturbances. Both the simulation and experimental results illustrate the effectiveness of the proposed controller in vibration suppression for a SAS system.


---

## Body

## 1. Introduction

Vibration is a common and unpredicted phenomenon for dynamic and static bodies. Vibrations are detrimental to comfort in many places, and most of the applications in civil, mechanical, and electrical engineering are easily affected by undesirable vibration or even the complete system will lose functionality. Those vibrations are harmful and must be eliminated or reduced. Vibration reduction can be achieved in many different ways. From the difference of isolation components, vibration suppression systems can be divided into three categories: passive, semiactive, and active.In the automobile industry, the vibration suppression system is introduced for decades. In a wheeled vehicle the mechanical system of springs and shock absorbers connect the wheels and axles to the chassis form vibration suspension system. The suspension system can provide stiffening and damping when the vehicle is running on an irregular road surface to isolate the vehicle body and ensure the comfort of the passengers. Many efforts were made to make the suspension system work in an optimal condition by optimizing the parameters of the suspension system passively; however, they have limitation of frequency range. Active suspensions were also introduced in [1–7], but the system is complex and in the need of much more power supply because of the added actuators. Semiactive concept combines the advantages of both active and passive suspensions and results in good performance with less complexity and power acquirement [8].One of the semiactive suspension systems currently used for vibration isolation is equipped with magnetorheological (MR) damper which creates braking torque by changing the viscosity of the MR fluid inside the brake [9]. MR fluid has magnetically sensitive rheological properties. Varying the magnetic field strength by changing the input current has the effect of changing the viscosity of the MR field and this leads to the changing of the damper torque output [10]. Thus, the output torque of the MR damper can be controlled by changing the input current. For a MR damper and spring suspension system, the input current must be controlled properly in order to suppress vibration. Usually a feedback control with the information of the body position is introduced, such as PID control [11], neural network control [12], backstepping control [13, 14], fuzzy logic control [15], LQG [16], and H∞ control [17].In this paper, a practical feedback control solution for controlling the MR damper based on extended state observer (ESO), a part of active disturbance rejection control (ADRC) technology [18, 19], is applied. The ADRC has been successfully employed in many mechanical and electronic systems [20–24]. A simulation is done based on the dynamic model of the semiactive suspension (SAS) system and ADRC technology. The experimental task is also carried out to measure the performance of the controller in the real system.


## 2. Semiactive Suspension System

The semiactive suspension (SAS) system studied here for vibration suppression is developed by the Polish Company Inteco Limited [25]. This SAS system can be used to analyze the vertical dynamics of the car wheel. As shown in Figure 1, the SAS laboratory model simulated a quarter of a wheeled vehicle, and it consists of an upper beam which represents the car body, a wheel, rotational MR damper, and a spring. It is driven by a DC motor with gear coupled to an eccentric small wheel. The suspended car wheel rolls due to the eccentric wheel rotation and oscillates up and down due to the small wheel eccentricity. The MR damper incorporated in SAS acts as an interface between sensors (encoders), control algorithms, and mechanical structure of the suspension, using the external damper coil current to adjust the damping. The torque generated by the MR damper depends on the rotary velocity of the damper and the magnetic field strength.Figure 1
Semiactive suspension (SAS) system.Figure2 depicts the geometrical view of the SAS system.Figure 2
Geometrical diagram of SAS.The dynamic model of the SAS is described by the following differential equations, in which the detailed definition of the anglesα
(
·
) and the distances r
(
·
) are referred to in the nomenclature.The dynamic equation of the upper beam is(1)
J
2
d
2
α
2
d
t
2
=
T
2
.T
2 is the total torque added to the upper beam. It consists of the friction torque of the upper beam at the torsional joint at the MR damper point, the moment dual to gravity, the torque caused by the connecting spring, and the MR damper torque. The T
2 is given as the following equation:
(2)
T
2
=
k
2
d
α
2
d
t
-
M
2
cos
⁡
α
2
+
r
2
k
s
l
s
+
(
d
α
1
d
t
-
d
α
2
d
t
)
M
MR
(
i
)
with
(3)
l
s
=
l
0
s
-
(
r
2
cos
⁡
α
2
-
r
1
cos
⁡
α
1
)
2
+
(
r
2
sin
α
2
-
r
1
sin
α
1
)
2
.The dynamic equation of the lower beam is(4)
J
1
d
2
α
1
d
t
2
=
T
1
.T
1 is the total torque added to the lower beam. It consists of the friction torque of the lower beam at the torsional joint at the MR damper point, the moment dual to gravity, the torque caused by the connecting spring, the elastic torque caused by the moving tire, the friction torque caused by the moving tire, and the MR damper torque. The T
1 is given as the following equation:
(5)
T
1
=
k
1
d
α
1
d
t
-
M
1
cos
⁡
(
β
1
-
α
1
)
-
r
1
k
s
l
s
+
k
g
R
cos
⁡
(
β
-
α
1
)
×
(
l
0
g
+
R
sin
(
β
-
α
1
)
+
r
-
D
x
+
u
kin
)
+
f
g
(
d
(
D
x
-
u
kin
)
d
t
-
d
a
1
d
t
R
cos
⁡
(
β
-
α
1
)
)
+
(
d
α
1
d
t
-
d
α
2
d
t
)
M
MR
(
i
)
.To emphasize the MR damper torque, we rewrite (1) and (4) depending on (2) and (5) and define
(6)
T
MR
=
(
d
α
1
d
t
-
d
α
2
d
t
)
M
MR
(
i
)
.We can give out the dynamic equation of the upper and lower beam as the following:(7)
J
2
d
2
α
2
d
t
2
-
k
2
d
α
2
d
t
+
M
2
cos
⁡
α
2
-
r
2
k
s
l
s
=
T
MR
,
(8)
J
1
d
2
α
1
d
t
2
-
k
1
d
α
1
d
t
+
M
1
cos
⁡
(
β
1
-
α
1
)
+
r
1
k
s
l
s
-
k
g
R
cos
⁡
(
β
-
α
1
)
×
(
l
0
g
+
R
sin
(
β
-
α
1
)
+
r
-
D
x
+
u
kin
)
-
f
g
(
d
(
D
x
-
u
kin
)
d
t
-
d
a
1
d
t
R
cos
⁡
(
β
-
α
1
)
)
=
T
MR
.


## 3. Feedback Control with Extended State Observer (ESO)

To suppress the vibration of the vehicle and make passengers feel comfortable mean that the aim of the control methodology is to decease the changing ofα
2 and try to make the system stable at certain degree. Therefore, in a feedback control method we mainly focus on the vibration of the body, so that we only use the dynamic of body movement.Equation (7) can be expressed as
(9)
J
2
α
¨
2
-
k
2
α
˙
2
+
M
2
cos
⁡
α
2
-
r
2
k
s
l
s
=
T
MR
(
i
)
.Definex
1
=
α
2 and x
2
=
x
˙
1
=
α
˙
2; then, we can rewrite (9) in the following:
(10)
x
˙
1
=
x
2
,
x
˙
2
=
f
(
x
1
,
x
2
,
d
(
t
)
)
+
b
u
,
y
=
x
1
,
where the ESO compensation coefficient b
=
1
/
J
2; control torque produced by MR damper u
=
T
MR
(
i
); d
(
t
)
=
r
2
k
s
l
s
/
J
2 refers to the external disturbance. f
(
x
1
,
x
2
,
d
(
t
)
)
=
f
(
α
2
,
α
˙
2
,
d
(
t
)
), or simply donated as f, is a function of both the internal dynamics and the external disturbance and
(11)
f
=
k
2
J
2
α
˙
2
-
M
2
J
2
cos
⁡
α
2
+
r
2
k
s
l
s
J
2
.The output of the system isy.The basic idea of the control methodology here is to obtain the estimatedf
(
or
f
^
) in the real-time by an extended state observer (ESO) and to actively compensate for it in the control law [18, 19].Now we construct the ESO [18, 19], in the form of
(12)
e
=
z
1
-
y
,
fe
=
fal
(
e
,
0.5
,
δ
)
,
fe
1
=
fal
(
e
,
0.25
,
δ
)
z
˙
1
=
z
2
-
β
01
e
,
z
˙
2
=
z
3
+
b
u
-
β
02
fe
z
˙
3
=
-
β
03
fe
1
with
(13)
fal
(
e
,
α
,
δ
)
=
{
e
δ
1
-
α
,
|
e
|
≤
δ
|
e
|
α
sign
⁡
(
e
)
|
e
|
>
δ
.There are many ways to select the observer gainsβ
01, β
02, and β
03 for a particular problem. In [19], the author proposes a group of gains when the sampling time h of the system (simulation or real experiment) is chosen; the gains can be given as follows:
(14)
β
01
=
1
h
,
β
02
=
1
1.6
h
1.5
,
β
03
=
1
8.6
h
2.2
.The inputs to ESO are the system outputy and the control signal u, and the output of the ESO gives the important information F
(
t
)
=
f
(
x
1
,
x
2
,
d
(
t
)
). We can set the control law as (
u
0
-
F
(
t
)
)
/
b to reduce the plant in (10) to a cascade integral form of
(15)
x
˙
1
=
x
2
,
x
˙
2
=
u
0
,
y
=
x
1
.Then, a PD controller is implemented to makeu
0 a function of the tracking error and its derivative, and we can control the output of system y to a near constant value to achieve the vibration suppression of the body. The whole process is shown in Figure 3.Figure 3
Vibration suppression feedback control of SAS based on ESO.


## 4. MR Damper Model

The MR damper is a type of semiactive damper where the viscosity of the MR fluid is controlled by varying the amount of the input current and hence changes the output torque. The output of the controller is torque and a MR damper model is needed to calculate the current both in simulation and experiment. To describe the hysteretic behavior of MR dampers, several models have been proposed in the literature, including the Bingham model, Dahl model, the Bouc-Wen model, Lurge model, and neural networks [10, 26]. In this paper Dahl model is selected to model the MR damper [27].The governing equations of the MR damper are(16)
T
MR
=
K
x
(
i
)
θ
˙
+
K
y
(
i
)
z
,
z
˙
=
ε
(
θ
˙
-
|
θ
˙
|
z
)
,
where θ is the angle, K
x is the damping coefficients which depend linearly on the current i, and z is the hysteretic variable. Parameters K
y and ε control the shape of the hysteresis curve. K
x and K
y are described by the following equations:
(17)
K
x
=
K
a
+
K
b
i
,
K
y
=
K
1
+
K
2
i
.Damper torque given from the Dahl model is as follows:(18)
T
MR
=
(
K
a
+
K
b
i
)
θ
˙
+
(
K
1
+
K
2
i
)
z
.Therefore, the control current that should be sent to the MR damper is(19)
i
=
T
MR
-
K
a
θ
˙
-
K
1
z
K
b
θ
˙
+
K
2
z
.We test and measure the MR damper parameters in the Dahl model from the SAS laboratory model and get them as the following:(20)
K
1
=
5
,
K
2
=
1.5
,
K
a
=
0.001
,
2222222
K
b
=
0.001
,
ε
=
5
.


## 5. Numerical Results

The system equations and feedback control method based on ESO are implemented in Matlab/Simulink to evaluate the system performance. In the simulation, we chose simulation environment as fix-step, with the ode3 solver and sampling timeh
=
0.01 s. The ESO gains can be calculated by (14), and the results are β
01
=
100, β
02
=
625, and β
03
=
2920.8. The proportional gain and the differential gain of the PD controller are chosen as P
=
1 and D
=
5 for the purpose of making the loop stable. Different excitation inputs, step, continuous sinusoidal, and random sequence are taken to show the control results. Figures 4, 5, 6, and 7 illustrate that with the feedback control method based on ESO, the body is stabilized more quickly and the vibration is suppressed. The transmissibility of excitement input (mm) to body angle output (degree) is calculated based on random simulation data. The frequency range is between 0 and 5 Hz, which is proper for the real excitation system. The simulation results in Figure 7 show that the vibration of the body is well suppressed at the frequency range by feedback control based on ESO.Figure 4
Step excitement input simulation results.Figure 5
Continuous sinusoidal excitement input simulation results.Figure 6
Random sequence excitement input simulation results.Figure 7
Transmissibility of excitement input (mm) to body angle output (degree) of simulation.


## 6. Experiment with Laboratory SAS System

The SAS has been designed to operate with an external PC-based digital controller. The control computer communicates with the level sensors, valves, and pump by a dedicated I/O board and the power interface. The I/O board is controlled by the real-time software which operates in Matlab/Simulink/RTW rapid prototyping environment. RTW is short for real-time workshop.From the MR damper model we can see that a velocity of angle is needed to calculate the current; however, the hardware I/O board just returns angle information. In the real signal processing, it is difficult to get the differential information of a certain signal. So here we introduce the tracking differentiator (TD for short) method to calculate velocity information from the angle. The tracking differentiator can be found in [18, 28]. The form of TD here is nonlinear because of its good property of tracking differential information of original signals which are covered by sensor white noise [28, 29]. The discrete signals v
1 and v
2 of TD are the traces of the reference input v following an optimal path:
(21)
v
1
(
k
+
1
)
=
v
1
(
k
)
+
h
v
2
(
k
)
v
2
(
k
+
1
)
=
v
2
(
k
)
+
h
fhan
(
v
1
(
k
)
-
v
(
k
)
,
v
2
(
k
)
,
r
0
,
h
0
)
,
where h is the sampling step, r is the tracking velocity factor for signal tracking, and h
0 is the filtering factor to reject stochastic noise. The update function fhan
(
p
1
,
p
2
,
r
,
H
) is of the form
(22)
fhan
(
p
1
,
p
2
,
r
,
H
)
=
{
-
r
0
sign
⁡
(
a
)
,
|
a
|
>
d
-
r
0
a
d
,
|
a
|
≤
d
with
(23)
sign
⁡
(
a
)
=
{
-
1
,
a
<
0
0
,
a
=
0
1
,
a
>
0
,
a
=
{
p
1
+
c
-
d
2
sign
⁡
(
γ
)
,
|
γ
|
>
d
2
p
2
+
γ
d
,
|
γ
|
≤
d
2
,
d
=
r
0
H
,
γ
=
r
0
p
1
+
d
p
2
,
c
=
d
2
+
8
|
γ
|
.The smoothed velocityx
2, which is the velocity of MR damper angle α
2, is derived from
(24)
x
1
(
k
+
1
)
=
x
1
(
k
)
+
h
x
2
(
k
)
x
2
(
k
+
1
)
=
x
2
(
k
)
+
h
fhan
(
x
1
(
k
)
-
α
2
(
k
)
,
x
2
(
k
)
,
r
0
,
h
0
)
.In the experiment the parameters are chosen as the following:r
0
=
5, h
=
0.01, and h
0
=
0.1. Figure 8 shows the results of MR damper angular velocity when input excitement is continuous sinusoidal.Figure 8
Tracking results of MR damper angular velocity with a continuous sinusoidal excitement.Feedback control method based on the ESO has been taken experiment on the SAS laboratory system, and the MATLAB/Simulink program environment was used to realize the control method. With the tracking differentiator, the MR damper model, PD control, ESO compensator, excitation input control, and the SAS plant, the control structure diagram of the experiment is shown in Figure9. And the control program diagram for experiment with SAS hardware interface in the Matlab/Simulink is shown in Figure 10.Figure 9
The control structure diagram of the experiment.Figure 10
Real-time experimental Simulink model of the SAS with feedback control base on ESO.The similar excitation inputs as in the numerical simulation are taken in the experiment in order to compare with the former simulation results. The vibrations are introduced by the eccentric wheel, so the excitements are always from 0 to 6 mm. The vibration response of the body angle is obtained with and without feedback controller based on ESO. Figures11, 12, 13, and 14 illustrate the experiment results.Figure 11
Step excitement input experimental resultsFigure 12
Continuous sinusoidal excitement input experimental results.Figure 13
Random excitement input experimental results.Figure 14
Transmissibility of excitement input (mm) to body angle output (degree) of experiment.According to the experimental results, as in the simulation, in the cases of step, continuous sinusoidal, and random sequence excitements, we can see that with the feedback control based on ESO the vibration of the upper beam (the body) is suppressed and the body is stabilized much faster than without control. The transmissibility of excitement input (mm) to body angle output (degree) is calculated based on random experiment data. The frequency range is between 0 and 5 Hz, which is proper for the real excitation system. The experiment results in Figure14 show that the vibration of the body is well suppressed at the frequency range by feedback control based on ESO. There was an input current limitation of the MR damper in the SAS system, so the control results in experiments did not achieve the fine ability as in simulation.The vibration response of a freely released beam after suppression is depicted in Figure15. With the controller based on ESO, the body is stabilized much faster.Figure 15
Experimental results of the body angle for releasing the beam after suppression.


## 7. Conclusion

The SAS system represents a quarter wheeled vehicle model, and it uses a MR damper to semiactively reduce the vibration. The output torque of the MR damper can be changed by the input of the current, and we use feedback controller to implant vibration suppression control to the SAS system based on the function of MR damper. The controller is based on ESO, the main part of ADRC technology.The proposed controller is tested in numerical simulation as well as real-time experiments with various vibration input signals. The simulation and experimental results both show the good quality of the controller. The controller depends on less information of system model and can observe inner and outer disturbances in real-time to counterweigh them, so it has potential to be applied to a whole vehicle model with a group of four independent controllers base on ESOs in the future research.



---

*Source: 248297-2014-05-20.xml*