Optimal control and parameters design for the fractional-order vehicle suspension system

Abstract

In this paper the optimal control and parameters design of fractional-order vehicle suspension system are researched, where the system is described by fractional-order differential equation. The linear quadratic optimal state regulator is designed based on optimal control theory, which is applied to get the optimal control force of the active fractional-order suspension system. A stiffness-damping system is added to the passive fractional-order suspension system. Based on the criteria, i.e. the force arising from the accessional stiffness-damping system should be as close as possible to the optimal control force of the active fractional-order suspension system, the parameters of the optimized passive fractional-order suspension system are obtained by least square algorithm. An Oustaloup filter algorithm is adopted to simulate the fractional-order derivatives. Then, the simulation models of the three kinds of fractional-order suspension systems are developed respectively. The simulation results indicate that the active and optimized passive fractional-order suspension systems both reduce the value of vehicle body vertical acceleration and improve the ride comfort compared with the passive fractional-order suspension system, whenever the vehicle is running on a sinusoidal surface or random surface.

Keywords

Fractional-order derivatives optimal control least square algorithm Oustaloup filter algorithm

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. Viscoelastic vibration isolation material is a major component of the vehicle suspension system, such as air springs, viscous dampers, magneto-rheological (MR) dampers. At present, air springs, and viscous dampers are widely applied to the design of suspension for many trucks, passenger cars and high-class cars. At the end of the 20th century, Lord company and Delphi company researched and developed different semi-active systems which contain MR dampers. Since then, MR dampers are also widely used in the design of vehicle suspension. Air springs, viscous dampers, and MR dampers are viscoelastic materials, so that it is not accurate to only use stiffness coefficients k or damping coefficients c to describe their dynamical properties. In 1936, Gemant firstly presented fractional-order derivative models of viscoelastic materials.¹ In 1991, Makris proposed a fractional-derivative Maxwell model of viscous dampers.² In 1995, Oustaloup presented the CRONE suspension which resulted from a traditional suspension model whose spring and damper were replaced by mechanical and hydropneumatic system defined by fractional-order force–displacement transmittance.³ In 2016, Chen used fractional-order differential equation to describe the dynamical behaviors of MR dampers.⁴ All these studies suggested that it is more reasonable to use fractional-order derivative model to describe the dynamical behaviors of viscoelastic materials than integer-order one. Accordingly, in order to study the vehicle suspension system conveniently and effectively, fractional-order derivative should be introduced into the model of vehicle suspension system. Nowadays, more and more attention has been paid on the fractional-order vehicle suspension system. Moreau and Altet researched the stability and control strategy of fractional-order vehicle suspension system.^5,6 Daou et al. emphasized the inertial effect and the nonlinear relations involved in the fractional-order vehicle suspension system.^7,8 Daou et al. also dealt with the effects of the structural uncertainties on the fractional-order vehicle suspension system.⁹

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. In other words, the parameters in the passive vehicle suspension system could not be changed with the external excitation, which limited the vehicle performance. The passive vehicle suspension system was most widely used because of its simple structure, reliable performance, and low cost. In recent years, along with the rapid development of vehicle technology, the requirements on ride comfort are getting higher and higher. Therefore, industrial countries had already begun to study the semi-active and active control systems, and the vehicle suspension system has experienced the changes from passive to semi-active and active forms.^10,11 Although there are so many significant works in the area of optimal control and parameters design of the traditional vehicle suspension system, there is little significant work on the problem of the optimal control and parameters design of the fractional-order vehicle suspension system.

In this paper, the optimal control and parameters design of the fractional-order vehicle suspension system are investigated. The primary aim of this study is to improve the ride comfort, which means that the vertical vibration acceleration of vehicle body is reduced while the suspension deformation and the dynamic load of tire are controlled in permissible ranges. The fractional-order derivatives in the system come from the viscoelasticity of isolation devices, and they could be solved by approximate analytical method¹² and Oustaloup filter algorithm respectively.¹³ The optimal control force of the active fractional-order vehicle suspension system is obtained based on optimal control theory, where the fractional-order derivatives in the system are replaced by the integer-order counterpart. A stiffness-damping system is added to the passive fractional-order vehicle suspension system, and the force arising from it should be as close as possible to the optimal control force of the active fractional-order suspension system. The passive fractional-order vehicle suspension system is optimized by using this method. Then, the parameters of the optimized passive fractional-order vehicle suspension system are obtained based on least square algorithm.

This paper is organized as follows. In “Dynamical model of active fractional-order vehicle suspension system” section, the dynamical model of active fractional-order vehicle suspension system is established. In “Processing fractional-order derivatives” section, the processing method of the fractional-order derivatives is presented. In “The solution of optimal control force” section, the optimal control force of the active fractional-order vehicle suspension system is obtained. In “Optimization of passive fractional-order vehicle suspension system” section, the parameter optimization of the passive fractional-order vehicle suspension system is discussed. In “Optimization example of fractional-order vehicle suspension system” section, the performances of the three kinds of fractional-order vehicle suspension systems are studied in time domain and frequency domain respectively. Finally, the conclusions are made in “Conclusions” section.

Dynamical model of active fractional-order vehicle suspension system

In this paper, a simplified quarter-vehicle model is selected as the research object, which is shown in Figure 1. According to Newtonian second law, the differential motion equation could be established as

{m_{1} {\overset{··}{z}}_{1} + k_{1} (z_{1} - z_{2}) + c_{1} ({\overset{\cdot}{z}}_{1} - {\overset{\cdot}{z}}_{2}) + {kD}_{t}^{p} (z_{1} - z_{2}) = u m_{2} {\overset{··}{z}}_{2} + k_{2} (z_{2} - q) + k_{1} (z_{2} - z_{1}) + c_{1} ({\overset{\cdot}{z}}_{2} - {\overset{\cdot}{z}}_{1}) + {kD}_{t}^{p} (z_{2} - z_{1}) = - u

(1)

where m₁, m₂, k₁, k₂, and c₁ are the sprung mass, the unsprung mass, suspension stiffness, tire stiffness, and suspension damping coefficients, respectively. Fractional-order derivative in this system comes from the viscoelasticity of isolation devices, such as the air spring or MR damper. The force–deformation relationship of viscoelastic materials in this system could be expressed as

F_{vm} = {kD}_{t}^{p} (z_{1} - z_{2})

, where

F_{vm}

is force,

z_{1} - z_{2}

is the suspension dynamic deformation, and

D_{t}^{p} (z_{1} - z_{2})

is the fractional-order derivative of

z_{1} - z_{2}

to t with the fractional coefficient

k (k > 0)

and the fractional order

p (0 \leq p \leq 1)

. u, q, z₁, and z₂ are the active control force, the road input, the displacement of vehicle body, and the tire, respectively.

Figure 1.

The 2-degree of freedom (DOF) model of quarter-vehicle body.

Processing fractional-order derivatives

The concept of fractional-order derivatives was firstly presented by Hospital and Leibniz in the late 1600s. After 300 years of development, several different definitions of fractional-order derivatives were presented by mathematicians.^14,15 Here Caputo’s definition is adopted

D_{t}^{p} x (t) = \frac{1}{Γ (1 - p)} \int_{0}^{t} \frac{\overset{\cdot}{x} (τ)}{(t - τ)^{p}} d τ

(2)

where

Γ (x)

is gamma function satisfying

Γ (x + 1) = x Γ (x)

At present, the computation methods of fractional-order derivatives are roughly divided into three kinds, i.e. approximate analytical method, numerical simulation, and filter algorithm.^16,17 In this paper, the approximate analytical method and filter algorithm are adopted.

Oustaloup filter algorithm of fractional-order derivatives

Under zero initial conditions, the Laplace transform of fractional-order derivatives is

L [D_{t}^{p} x (t)] = s^{p} X (s)

(3)

At present, the approximate filter algorithm of fractional-order derivatives includes three forms, i.e. continued fraction approximation,¹⁸ Charef approximation,¹⁹ and Oustaloup approximation.²⁰ In this paper, Oustaloup approximation is adopted. The principle of Oustaloup filter algorithm is that using higher-order transfer function approximately to replace the fractional-order operator $s^{p}$ . Assuming that the fitting frequency band is ( $ω_{b}$ , $ω_{h}$ ), and the filter order is N, Oustaloup filter can be obtained as following

s^{p} \approx ω_{h}^{p} Π_{i = 1}^{N} \frac{s + ω_{i}^{*}}{s + ω_{i}}

(4)

where

ω_{i}^{*} = ω_{b} (\sqrt{ω_{h} / ω_{b}})^{(2 i - 1 - p) / N}

(5)

ω_{i} = ω_{b} (\sqrt{ω_{h} / ω_{b}})^{(2 i - 1 + p) / N}

(6)

If Oustaloup filter order N is 3, 5, and 9, respectively, the effects of the filter approximation are shown in Figure 2 for $p = 0.5$ , $ω_{b} = 0.001 rad / s$ , $ω_{h} = 1000 rad / s$ . Comparing the filter approximation effects under different filter orders, it can be concluded that Oustaloup filter algorithm has a high accuracy in the selected frequency range when $N = 9$ .

Figure 2.

Comparisons of Bode diagrams with different filter orders.

Approximate analytical method of fractional-order derivatives

Substituting $x (t) = a \sin (ω t + ϕ)$ into equation (2), the following equation can be obtained

D_{t}^{p} x (t) = \frac{a ω}{Γ (1 - p)} \int_{0}^{t} \frac{\cos (ω τ + ϕ)}{(t - τ)^{p}} d τ

(7)

Let $t - τ = y$ , equation (7) can be rewritten as

D_{t}^{p} x (t) = \frac{a ω \cos ϕ}{Γ (1 - p)} \int_{0}^{t} \frac{\cos (ω t - ω y)}{y^{p}} d y - \frac{a ω \sin ϕ}{Γ (1 - p)} \int_{0}^{t} \frac{\sin (ω t - ω y)}{y^{p}} d y

(8)

Differentiating equation (8), the following result can be obtained

\frac{d {[D}_{t}^{p} x (t)]}{dt} = \frac{a ω \cos ϕ}{Γ (1 - p)} [- ω \sin (ω t) \int_{0}^{t} \frac{\cos (ω y)}{y^{p}} d y + \frac{1}{t^{p}} + ω \cos (ω t) \int_{0}^{t} \frac{\sin (ω y)}{y^{p}} d y] - \frac{a ω \sin ϕ}{Γ (1 - p)} [ω \cos (ω t) \int_{0}^{t} \frac{\cos (ω y)}{y^{p}} d y + ω \sin (ω t) \int_{0}^{t} \frac{\sin (ω y)}{y^{p}} d y]

(9)

Introducing two formulae²¹

\int_{0}^{t} \frac{\cos (ω y)}{y^{p}} d y = ω^{p - 1} {Γ (1 - p) \sin (\frac{p π}{2}) + \frac{\sin (ω t)}{(ω t)^{p}} + O [(ω t)^{- p - 1}]}

(10a)

\int_{0}^{t} \frac{\sin (ω y)}{y^{p}} d y = ω^{p - 1} {Γ (1 - p) \cos (\frac{p π}{2}) - \frac{\cos (ω t)}{(ω t)^{p}} + O [(ω t)^{- p - 1}]}

(10b)

where

O [(ω t)^{- p - 1}]

is high-order infinitesimal of

(ω t)^{- p - 1}

. And substituting equation (10) into equation (9), the following result can be obtained

\frac{d {[D}_{t}^{p} x (t)]}{dt} = - a ω^{p + 1} [\sin (\frac{p π}{2}) + \frac{O ((ω t)^{- p - 1})}{Γ (1 - p)}] \sin (ω t + φ) + a ω^{p + 1} [\cos (\frac{p π}{2}) + \frac{O ((ω t)^{- p - 1})}{Γ (1 - p)}] \cos (ω t + φ)

(11)

Here we assume that

D_{t}^{p} x (t) = bx (t) + c \overset{\cdot}{x} (t) + df (t)

(12)

Differentiating equation (12), the following equation can be established

\frac{d {[D}_{t}^{p} x (t)]}{dt} = ab ω \cos (ω t + ϕ) - ac ω^{2} \sin (ω t + ϕ) + d \overset{\cdot}{f} (t)

(13)

Because equation (11) and equation (13) are equal, one can obtain

b = ω^{p} \cos (\frac{p π}{2}) + \frac{O (t^{- p - 1})}{ω Γ (1 - p)}

(14a)

c = ω^{p - 1} \sin (\frac{p π}{2}) + \frac{O (t^{- p - 1})}{ω^{2} Γ (1 - p)}

(14b)

d = 0

(14c)

Therefore any trigonometric function meets the following equality

D_{t}^{p} x (t) = [ω^{p} \cos (\frac{p π}{2}) + \frac{O (t^{- p - 1})}{ω Γ (1 - p)}] x (t) + [ω^{p - 1} \sin (\frac{p π}{2}) + \frac{O (t^{- p - 1})}{ω^{2} Γ (1 - p)}] \overset{\cdot}{x} (t)

(15)

The solution of optimal control force

The fractional-order derivatives in the system can be replaced by integer-order counterpart, when the vehicle is running on a sine surface.²² In order to use optimal control theory to find the explicit form of the optimal control force, equation (1) should be rewritten as the state equation. Considering state vector as

X = [x_{1} x_{2} x_{3} x_{4}]^{T} = [z_{1} - z_{2} {\overset{\cdot}{z}}_{1} - {\overset{\cdot}{z}}_{2} z_{2} - q {\overset{\cdot}{z}}_{1}]^{T}

(16)

and the output vector

Y = [y_{1} y_{2} y_{3}]^{T} = [x_{1} x_{3} {\overset{\cdot}{x}}_{4}]^{T}

(17)

the system state equation can be obtained as

\overset{\cdot}{X} = AX + Bu + Lw

(18)

Y = CX + Du

(19)

where

A = (0100 - \frac{k_{1}}{m_{1}} - \frac{k_{1}}{m_{2}} - (\frac{k}{m_{1}} + \frac{k}{m_{2}}) ω^{p} \cos (\frac{p π}{2}) - \frac{c_{1}}{m_{1}} - \frac{c_{1}}{m_{2}} - (\frac{k}{m_{1}} + \frac{k}{m_{2}}) ω^{p - 1} \sin (\frac{p π}{2}) \frac{k_{2}}{m_{2}} 00 - 1 01 - \frac{k_{1}}{m_{1}} - \frac{k}{m_{1}} ω^{p} \cos (\frac{p π}{2}) - \frac{c_{1}}{m_{1}} - \frac{k}{m_{1}} ω^{p - 1} \sin (\frac{p π}{2}) 00) B = (0 \frac{1}{m_{1}} + \frac{1}{m_{2}} 0 \frac{1}{m_{1}}), L = (00 - 1 0), D = (00 \frac{1}{m_{1}}), w = \overset{\cdot}{q} C = (10000010 - \frac{k_{1}}{m_{1}} - \frac{k}{m_{1}} ω^{p} \cos (\frac{p π}{2}) - \frac{c_{1}}{m_{1}} - \frac{k}{m_{1}} ω^{p - 1} \sin (\frac{p π}{2}) 00)

The primary aim of this study is to find out an optimal control force which makes the active fractional-order vehicle suspension system achieve the optimal performance. The performance of vehicle suspension system can be quantitatively evaluated in three indexes, i.e. the vertical acceleration of vehicle body which represents the ride comfort, the suspension dynamic deformation which represents the body posture and the suspension structure, and the dynamic load of the tire which represents the tire grounding characteristics. But these three performance indexes are normally contradictory. Therefore, the target performance index J can be set as the mean integral value of the weighted square sum of the active force, the body vertical vibration acceleration, the dynamic deformation of the suspension, and tire

J = \frac{1}{2} \int_{0}^{\infty} [q_{1} y_{1}^{2} + q_{2} y_{2}^{2} + q_{3} y_{3}^{2} + {ru}^{2}] d t

(20)

where q₁, q₂, q₃, and r are weighting coefficients. Substituting equation (19) into equation (20), the following equation can be obtained

J = \frac{1}{2} \int_{0}^{\infty} [X^{T} C^{T} QCX + u^{T} Ru + 2 X^{T} Nu] d t

(21)

where

Q = (q_{1} 000 q_{2} 000 q_{3}), R = D^{T} QD + r, N = C^{T} QD

Based on optimal control theory, there is an optimal control force

u^{*} = - KX

(22)

where the feedback gain vector is

K = R^{- 1} (B^{T} P + N^{T}) = [k_{11} k_{12} k_{13} k_{14}]

The matrix P can be achieved from the Riccati matrix algebraic equation

PA + A^{T} P - (PB + N {) R}^{- 1} (B^{T} P + N^{T}) + C^{T} QC = 0

(23)

Optimization of passive fractional-order vehicle suspension system

A stiffness-damping system is added to the passive fractional-order vehicle suspension system, and the force arising from it should be as close as possible to the optimal control force of the active fractional-order vehicle suspension system. By using this method, the parameters of the optimized passive vehicle suspension system can be obtained. The force arising from the accessional devices is

u' = - k' (z_{1} - z_{2}) - c' ({\overset{\cdot}{z}}_{1} - {\overset{\cdot}{z}}_{2})

(24)

where

k'

and

c'

are the stiffness and damping coefficients of the accessional devices, respectively. The closer

f = | u^{*} - u' |

gets to zero, the better the performance of optimized passive vehicle suspension system will be. Substituting equation (22) and equation (24) into f, the following equation can be obtained

f = | (k' - k_{11}) z_{1} - (k' - k_{11} + k_{13}) z_{2} + (c' - k_{12} - k_{14}) {\overset{\cdot}{z}}_{1} - (c' - k_{12}) {\overset{\cdot}{z}}_{2} + k_{13} q |

(25)

Based on least square algorithm, the following function H can be constructed

H = λ_{1} (k' - k_{11})^{2} + λ_{2} (k' - k_{11} + k_{13})^{2} + λ_{3} (c' - k_{12} - k_{14})^{2} + λ_{4} (c' - k_{12})^{2} + λ_{5} k_{13}^{2}

(26)

where

λ_{i} > 0

(i = 1, 2, 3, 4, 5) are weighting coefficients. Taking the first-order partial derivative of equation (26), and let it equal to zero, one can obtain

{k' = k_{11} - \frac{λ_{2}}{λ_{1} + λ_{2}} k_{13} c' = k_{12} + \frac{λ_{3}}{λ_{3} + λ_{4}} k_{14}

(27)

λ_{i}

is generally adjusted relying on the designer’s own experience, which is subjective and blind. Therefore the approximation precision is low. In this paper, because of we choose the state vector X reasonably,

λ_{2} / (λ_{1} + λ_{2})

and

λ_{3} / (λ_{3} + λ_{4})

are both changed in

[0, 1]

, which reduces the subjective speculation and improves the approximation precision. By appropriately selecting

λ_{2} / (λ_{1} + λ_{2})

and

λ_{3} / (λ_{3} + λ_{4})

k'

and

c'

can be achieved from equation (27). If

k' > 0

c' > 0

, the corresponding spring or damper should be added to the system. On the contrary, the stiffness coefficient or the damping coefficient of the system should be reduced.

Optimization example of fractional-order vehicle suspension system

A set of suspension parameters is selected, which is shown in Table 1. The main factors affecting the ride comfort include the road input, the performance of tire, and suspension. Among them, the suspension performance has the biggest influence on the ride comfort. In this paper, the main purpose is to improve the ride comfort, and coordinate the relationship between the ride comfort, the driving safety, and the tire life. Through a lot of trials, the weighting coefficients can be selected as

q_{1} = 6.5 \times 10^{9}

q_{2} = 3 \times 10^{10}

q_{3} = 10^{6}

r = 1

λ_{1} = 1

λ_{2} = 20

λ_{3} = 1

, and

λ_{4} = 50

, respectively. When the vehicle is running on a sinusoidal surface, i.e. a sinusoidal wave with the frequency

ω = 30 rad / s

, the system feedback gain matrix K can be obtained by using MATLAB/LQR

K = [451 - 875 38202194]

Table 1.

The parameters of the fractional-order vehicle suspension system.

$m_{1} / kg$	240
$m_{2} / kg$	36
$k_{1} / (N \cdot m^{- 1})$	15,000
$k_{2} / (N \cdot m^{- 1})$	150,000
$c_{1} / (N \cdot s \cdot m^{- 1})$	1500
k	1000
p	0.5

The parameters of the accessional stiffness-damping system can be achieved from equation (27)

k' = - 3187, c' = - 832

Next, the simulation model of the optimized passive fractional-order vehicle suspension is developed, when the vehicle is running on a sinusoidal surface and a random surface, respectively. Its simulation results are compared with those of the original passive and the active systems in time domain and frequency domain, respectively. In the random case the main evaluation indexes of the vehicle suspension system are root mean square (RMS) of steady-state vertical vibration acceleration of the vehicle body, the suspension dynamic deformation and the dynamic load of the tire. The simulation results are shown in Figure 3 to Figure 6, when the vehicle is running on a sinusoidal surface whose amplitude is $0.005 m$ and the frequency changes in $[0 Hz, 30 Hz]$ . The following conclusion can be obtained from Figure 3 to Figure 6.

The force $u^{*}$ and $u'$ are very close except at the resonance frequencies. At the same time, the performances of the active and optimized passive suspension systems are also very close except at the resonance frequencies.

The active and optimized passive suspension systems both significantly reduce the vehicle body vertical acceleration, and ride comfort is improved significantly. In this regard, their effects are very well.

In regard of the dynamic deformation of suspension, the effects of the active and optimized passive suspension systems are not good compared with the passive suspension system. But the weakness is not so significant to affect the whole performance of the suspension system.

In regard of the dynamic load of the tire, the active and optimized passive suspension systems present good effect compared with the passive suspension system in $[0 Hz, 8 Hz]$ , but it is contradictory in $[8 Hz, 20 Hz]$ . The three different control strategies almost have the same effects, when the excited frequency is too large. The dynamic load of tire will not affect the driving safety, as long as it is changed in the allowed range.

Figure 3.

Root mean square (RMS) of the force $u'$ and $u^{*}$ .

The filtering white noise is used to simulate the random road input. The mathematical model in time domain is

\overset{\cdot}{q} (t) = 2 π η_{0} \sqrt{G_{φ} (η_{0}) v} w (t) - 2 π η_{1} vq (t)

(28)

where

η_{1}

is the lower spatial cut-off frequency, and it is usually taken as

0.011 / m

. v is the vehicle speed,

η_{0}

is the reference spatial frequency,

G_{q} (η_{0})

is the pavement roughness coefficient shown in Table 2, and

w (t)

is the filtering white noise with the noise power 0.5.

Table 2.

The pavement roughness coefficient.

The road level	A	B	C	D
$G_{q} (η_{0}) \times 10^{- 6} / m^{3}$	16	64	265	1024

First the aforementioned Oustaloup filter algorithm is adopted to simulate fractional-order derivatives. And then, the simulation model of optimized passive fractional-order vehicle suspension system is developed. The simulation results are compared with those of the active and passive suspension systems, and they are summarized in Table 3, where the vehicle is running on different road levels from A to D with

v = 20 m / s

. The B-class road model is shown in Figure 7. The time histories of the vertical accelerations of the vehicle body under three different control strategies are shown in Figure 8 to Figure 10, when the vehicle is running on B-class road with

v = 20 m / s

. From Figure 8, Figure 9, Figure 10, and Table 3, it can be seen that the active and optimized passive suspension systems could significantly reduce vertical acceleration of the vehicle body. In regard of the dynamic deformation of suspension and dynamic load of tire, the effects of the active and optimized passive suspension systems are not good compared with the passive suspension system. However, the weakness is not so significant to affect the whole performance of the suspension.

Table 3.

RMS of fractional-order vehicle suspension performance indicators.

Road level	Suspension type	Vertical acceleration, ( $m / s^{2}$ )	Suspension deformation, $10^{- 3} m$	Dynamic load of tire, $N$
A	Passive	0.39	2.2	152
	Active	0.267	2.4	185
	Optimized	0.313	3.1	183
B	Passive	0.78	4.4	304
	Active	0.535	4.8	369
	Optimized	0.625	6.2	366
C	Passive	1.57	8.7	606
	Active	1.07	9.6	743
	Optimized	1.25	12.5	731
D	Passive	3.12	17.5	1212
	Active	2.13	19.2	1485
	Optimized	2.5	25	1462

RMS: root mean square.

Figure 4.

Root mean square (RMS) of acceleration.

Figure 5.

Root mean square (RMS) of deformation of suspension.

Figure 6.

Root mean square (RMS) of dynamic load of tire.

Figure 7.

B-class road model.

Figure 8.

Passive fractional-order vehicle suspension.

Figure 9.

Optimized passive fractional-order vehicle suspension.

Figure 10.

Active fractional-order vehicle suspension.

Conclusions

In this paper, the optimal control force of the active fractional-order vehicle suspension system is obtained based on optimal control theory, where the fractional-order derivative in the system is replaced by the integer-order counterpart. A stiffness-damping is added to the passive fractional-order vehicle suspension system and the force arising from it should be as close as possible to the optimal control force of the active fractional-order vehicle suspension system. The passive fractional-order vehicle suspension system is optimized by using this method. Then, the parameters of the optimized passive fractional-order vehicle suspension system are obtained based on least square algorithm. Finally, we compared the response of three kinds of vehicle suspension systems. The simulation results indicate that the active and optimized passive suspension systems both reduce the values of the vertical acceleration of vehicle body and improve the ride comfort compared with the passive suspension system, whenever the vehicle is running on a sine surface or random surface. This paper provides a guideline for the optimal control and parameters design of fractional-order vehicle suspension system, and also provides a simple and easy method for solving the similar engineering problems.

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: National Natural Science Foundation of China (No. 11372198), the Cultivation plan for Innovation team and leading talent in Colleges and universities of Hebei Province (LJRC018), the Program for advanced talent in the universities of Hebei Province (GCC2014053), and the Program for advanced talent in Hebei Province (A201401001).

References

Gemant

. A method of analyzing experimental results obtained from elasto-viscous bodies. Physics 1936; 7: 311–317.

Makris

Constantinou

. Fractional-derivative Maxwell model for viscous dampers. J Struct Eng 1991; 177: 2708–2724.

Oustaloup

Moreau

Nouillant

. The CRONE suspension. Control Eng Pract 1996; 4: 1101–1108.

Chen

Wilson

, et al. Fractional modeling and analysis of coupled MR damping system. IEEE/CAA J Autom Sinica 2016; 3: 288–294.

Moreau

Altet

Oustaloup

. The CRONE suspension: management of the dilemma comfort-road holding. Nonlinear Dyn 2004; 38: 461–484.

Altet

Moreau

Moze

, et al. Principles and synthesis of hydractive CRONE suspension. Nonlinear Dyn 2004; 38(1–4): 435–459.

Daou

RAZ

Francis

Moreau

. Study of the inertial effect and the nonlinearities of the CRONE suspension based on the hydropneumatic technology. Nonlinear Dyn 2011; 63: 1–17.

Daou

RAZ

Moreau

Francis

. Effect of hydropneumatic components nonlinearities on the CRONE suspension. IEEE Trans Veh Technol 2012; 61: 466–474.

Daou

RAZ

Moreau

Francis

. Study of the effects of structural uncertainties on a fractional system of the first kind-application in vibration isolation with the CRONE suspension. Signal Image Video Process 2012; 6: 463–478.

10.

Shen

Liu

Yang

. Research on the application of optimal control theory on parameters optimization of vehicle suspension. J Low Freq Noise V A 2003; 22: 253–263.

11.

Poussot-Vassal

Spelta

Sename

, et al. Survey and performance evaluation on some automotive semi-active suspension control methods: a comparative study on a single-corner model. Annu Rev Control 2012; 36: 148–160.

12.

Shen

Yang

Xing

, et al. Primary resonance of Duffing oscillator with fractional-order derivative. Commun Nonlinear Sci Numer Simul 2012; 17: 3092–3100.

13.

Oustaloup

Levron

Mathieu

, et al. Frequency-band complex non-integer differentiator: characterization and synthesis. IEEE Trans Circuit Syst I 2000; 47: 25–39.

14.

Podludny

. Fractional differential equations, London: Academic Press, 1999.

15.

Machado

Kiryakova

Mainardi

. Recent history of fractional calculus. Commun Nonlinear Sci Numer Simul 2011; 16: 1140–1153.

16.

Atanackovic

Stankovic

. On a numerical scheme for solving differential equations of fractional order. Mech Res Commun 2008; 35: 429–438.

17.

Shakoor

Ricardo

Delfim

FMT

. Numerical approximations of fractional derivatives with applications. Asian J Control 2013; 15: 698–712.

18.

Podlubny

Petraš

Vinagre

, et al. Analogue realizations of fractional-order controllers. Nonlinear Dyn 2002; 29: 281–296.

19.

Charef

Sun

Tsao

, et al. Fractal system as represented by singularity function. IEEE Trans Automat Control 1992; 37: 1465–1470.

20.

Oustaloup

Levron

Mathieu

, et al. Frequency-band complex noninteger differentiator: characterization and synthesis. IEEE Trans Circuits Syst I 2000; 47: 25–39.

21.

Chen

Zhu

. Stochastic dynamics and fractional optimal control of quasi integrable Hamiltonian systems with fractional derivative damping. Fract Calc Appl Anal 2013; 05: 189–225.

22.

You

Shen

Yang

. Optimal design for fractional-order active isolation system. Adv Mech Eng 2015; 7(12): 1–11.