Parameters of desired prescribed performance and applied prescribed performance.

## Abstract

This work proposes two novel adaptive fuzzy controllers and applies them to vibration control of a vehicle seat suspension system subjected to severe road profiles. The first adaptive controller is designed by considering prescribed performance of the sliding surface and combined with adaptation laws so that robust stability is guaranteed in the presence of external disturbances. As for the second adaptive controller, both the H-infinity controller and sliding mode controller are combined using inversely fuzzified values of the fuzzy model. In order to evaluate control performances of the proposed two adaptive controllers, a semi-active vehicle suspension system installed with a magneto-rheological (MR) damper is adopted. After determining control gains, two controllers are applied to the system and vibration control performances such as displacement at the driver’s position are evaluated and presented in time domain. In this work, to demonstrate the control robustness two severe road profiles of regular bump and random step wave are imposed as external disturbances. It is shown that both adaptive controllers can enhance ride comfort of the driver by reducing the displacement and acceleration at the seat position. This excellent performance is achieved from each benefit of each adaptive controller; accurate tracking performance of the first controller and fast convergence time of the second controller.

### Keywords

- adaptive fuzzy control
- sliding mode control
- H-infinity control
- prescribed performance of the sliding surface
- vibration control
- seat suspension system

## 1. Introduction

Nowadays, modern control-based technical devices such as robotics, assistive machines and home appliances are popularly used to improve the level of human being’s life. In these devices, control algorithm is one of the most important components which brings comfortable requirements to the consumer. The development of control algorithms in recent years is abundantly being undertaken from the aspect of classical control to salient characteristics of intelligent control. The classical control methods are frequently combined with modern control technique to resolve parameter uncertainties and disturbances those are existed in most of control devices. A controller which is formulated using more than two different control schemes is called “a hybrid controller” or “composite controller” [1, 2]. Among many candidates of the hybrid controller, the type of hybrid adaptive controller is the most popular since its structure is relatively simple and its control performance is very robust against the uncertainties or/and external disturbances. A hybrid adaptive control with fuzzy model and wavelet neural networks was presented in [1, 3] in which the sliding mode control was used to connect the parameters of the fuzzy model and the neural networks. This method is the typical model to develop the adaptive control in the last few years. Besides of uncertain nonlinear system, the problem of unknown input nonlinearity such as dead-zone or backlash-like hysteresis was also studied through the hybrid adaptive control [4]. It has been also shown that the neural works can be designed for a good performance of the hybrid adaptive control to deal with the uncertain system [5]. A hybrid adaptive controller possessing the robustness against input and parameter uncertainties was studied using the sliding mode controller associated with the fuzzy model [6, 7]. When a hybrid adaptive controller is formulated, in general the adaptation laws are simultaneously calculated. Furthermore, the back-stepping method was integrated with the fuzzy mode to achieve high performance of the hybrid adaptive controller [8].

As mentioned earlier, both the fuzzy model and the neural networks model are frequently used for the formulation of high performance of a hybrid adaptive controller [9]. Recently, a modified type of the fuzzy model called interval type 2 was combined with the back-stepping method to design of a hybrid adaptive control [10, 11]. It is remarked that the fixed fuzzy model always provides a safe choice in design of a hybrid adaptive control. However, this choice may cause a large error in finding the final values. To resolve this problem, an adaptive interval type 2 fuzzy neural network was developed on the basis of the online technique which can strengthen the flexibility of design parameters against the uncertainties [12]. Besides the above, there are many approaches to formulate new hybrid adaptive controllers such as output feedback control approach to take account for unknown hysteresis [13]. From the aspect of experimental implementation of hybrid adaptive controllers, several dynamic systems featuring magneto-rheological (MR) mount and MR damper are adopted for vibration control [2, 14–18]. Most of hybrid adaptive controllers used in these experimental realizations have been formulated by combining the models of interval type 2 fuzzy and interval type 2 fuzzy neural networks, and the control techniques of H-infinity control and sliding mode control. The advantage of using the interval type 2 fuzzy model is its flexibility in which optimized fuzzy values can be achieved unlike the classical fuzzy rule with the fixed value [19]. In order to improve the fuzzy model, clustering method [20] and data-driven for fuzzy rules [21] were also introduced.

As a subsequent work to develop a new hybrid adaptive controller, in this work two different new hybrid adaptive controllers are developed and their control performances are evaluated by investigation on vibration control of a semi-active seat suspension system installed with MR damper. The first hybrid adaptive controller is designed by combing online interval type 2 fuzzy neural networks model and prescribed performance of the sliding surface associated with adaptation laws to guarantee robust stability (HAC-PP in short). The second hybrid adaptive controller is formulated by combining inversely fuzzified value with H-infinity control to minimize computational cost algorithm (HAC-IFV in short).The stability of both adaptive controllers are rigorously proved based on the Lyapunov stability and appropriate control gains are determined to evaluate vibration control performance. It is shown that both proposed adaptive controllers are very effective and robust for controlling unwanted vibrations or excitations from the road profiles. These are validated by presenting control results showing significant reduction of both the displacement and acceleration at the seat position subjected to external excitations.

## 2. Formulation of HAC-PP

As mentioned in Introduction, the online interval type 2 fuzzy neural networks (OIT2FNN in short) model is used to formulate two adaptive controllers. The rule base of OIT2FNN can be expressed as follows [22].

where, * m*is the number of rules, and

In the above,

As a problem formulation, consider a single-input and single-output (SISO) nonlinear system governed by the following equation:

where ** f**(

**) ∈**x

*and*R

^{n}

**(**g

**) ∈**x

*are two unknown non-linear function vectors,*R

^{n}

*(*u

*) ∈*t

R

^{1}is control function,

**(**d

**) ∈**t

*is an external disturbance vector, |*R

^{n}

**(**d

**)| ≤**t

**where**δd

**∈**δd

*is upper bound of*R

^{n}

**(**d

**),**t

*is defined as follows:*s

_{s}

where, ** K** = [

*,*k

_{n}

k

_{n − 11},

k

_{n − 2}, …,

k

_{1}] is defined as the coefficients such that all of the roots of the polynomial

*+*σ

^{n}

k

_{n − 1}

σ

^{n − 1}+

k

_{n − 2}

σ

^{n − 2}+ … +

k

_{1}are in the open left-half complex plane. The sliding surface (4) is rewritten using the state variables as follows:

A new vector

where,

The tracking error is defined as * e* =

x

_{1}−

*with the desired states of*x

_{d}

*. Then the error performance function is defined as follows [23]:*x

_{d}

where, * l* > 0, 0 < |

*(0)| <*e

*(0),*λ

λ

_{∞}> 0,

λ

_{∞}<

*(0) then*λ

*> 0 and*λ

_{t}

*(*λ

*) tend to*t

λ

_{∞}exponentially. In order to guarantee fast convergence of tracking error, and obtain a certain convergence accuracy, the tracking error is set as follows:

In the above, the prescribed error performance function * S*(

*) found as follows:*φ

The function * S*(

*) must satisfy the following conditions.*φ

(i) * S*(

*) is smooth continuous and monotone increasing function*φ

(iii) lim_{φ → + ∞}* S*(

*) = 1 and lim*φ

_{φ → − ∞}

*(*S

*) = − 1*φ

From the above conditions (10), the function * S*(

*) can be determined as follows:*φ

Then using Eq. (8), the tracking error is obtained by

Hence, the tracking error can be summarized as * Ξ* = {

*∈*e

*: |*R

*(*e

*)| <*t

*∀*λ

*≥ 0 and*t

*(*e

*) <*t

λ

_{∞}for

*→ ∞}. On the other hand, the inverse function of (11) is expressed as:*t

Hence, the derivatives of Eq. (13) are obtained as:

where,

In order to realize * φ* → 0, the second sliding surface is defined as follows:

where * c* > 0. The derivative of Eq. (16) is obtained as:

_{s}

The lumped uncertainty of system is defined as:

where * γ* =

_{f}

*(*f

*) −*x

f

^{∗}(

*),*x

*=*γ

_{g}

*(*g

*) −*x

g

^{∗}(

*). Using Eqs. (17) and (18), the derivative Eq. (17) is rewritten as:*x

Based on Eq. (2), the relationship between Eq. (19) and OIT2FNN is expressed as follows:

where

* Δθ* = {

_{f}

*∈*θ

_{f}

*, ‖*R

^{n}

*‖ ≤*θ

_{f}

*},*Θ

_{f}

*= {*Δθ

_{g}

*∈*θ

_{g}

*, ‖*R

^{n}

*‖ ≤*θ

_{g}

*},*Θ

_{g}

*= {*Δx

*∈*x

*, ‖*R

^{n}

*‖ ≤*x

*}*Θ

_{x}

Now, an equivalent control is determined from Eq. (20) based on the assumption

The equivalent control _{1} cannot control the system because it cannot compensate the error from the fuzzy approximation. To guarantee the robustness and stability in control, a robust control part _{2} should be introduced as follows:

Then, the total control * u*of the system is determined as follows:

The control _{2} is the combination of two sliding surfaces * s*and

_{s}

*. The value*σ

_{s}

*is the adaptive parameter where its boundary is given by*Γ

*= {*ΔΓ

*∈*Γ

*, ‖*R

*‖ ≤*Γ

*,*Θ

_{Γ}

*≤*σ

_{s}Γξ

_{fz}

*}, and*ρ

*is constant boundary. The matrix*Θ

_{Γ}

**=**P

P

*≥ 0 in which its result is a solution of Riccati-like equation given by*

^{T}

where, * ρ* ≥

*,*σ

_{s}Γξ

_{z}

*is the prescribed attenuation level,*ρ

**=**Q

Q

^{T}≥ 0,

*is consequent membership value of the OIT2FNN. When the value*ξ

_{z}

*=*ρ

*, the Riccati-like equation is rewritten as:*σ

_{s}Γξ

_{z}

Now, Eq. (20) can be analyzed as follows:

where

Now, the stability of the proposed adaptive control system can be solidly proved with Eqs. (21)–(23) and adaptation laws as follows:

In order to make a proof, in this work the following Lyapunov function candidate is proposed.

The derivative of Eq. (29) is then obtained by

Substituting Eq. (27) into Eq. (30), Eq. (30) is rewritten as follows:

It is noted that Eq. (24) is used in finding Eq. (31). Substituting Eq. (28) into Eq. (31), the following is achieved.

Eq. (32) cannot use for conclusion of stability. Hence, it will be integrated from * t* = 0 to

*=*t

*, we have:*T

where, * L*(

_{v}

*) is always positive, so Eq. (33) is determined as:*T

From Eqs. (32) and (34), the stability is guaranteed.

From the boundedness of the parameters

In here, ℵ_{f}, ℵ_{g}, _{1}, _{2} are the choosing parameters. Hence, the adjusted adaptation laws are redefined as follows:

In the above, _{1}, _{2} and _{3} are choosing parameters related boundaries of * f*(

*),*x

*(*g

*) and*x

*. It is noted here that in order to utilize the states of the system, the Luenberger observer [24] has been used in this work. Figure 1 presents a flow chart of the HAC-PP showing the combination of each controller and the prescribed performance.*Γ

## 3. Formulation of HAC-IFV

As a first step to design the controller, consider the system (3) rewritten by

where, the function _{0}(** x**) and

g

_{0}(

**) are the functions of**x

**(**f

**) and**x

**(**g

**) which are determined as:**x

In the above, ** δf**and

**are two positive vectors. It is noted that**δg

**=**D

**+**δf

δg

*(*u

*) +*t

**(**d

**) denotes the uncertain disturbance and**t

**= [0, 0, …,**D

D

_{0}]

*. In order to formulate the controller, the following assumption is made: There exists a constant*

^{T}

*∈*g

_{m}

ℜ

^{+}to satisfy |

*(*g

*)| >*x

*. Without loss of generality, it is assumed that the equation*g

_{m}

*(*g

*) >*x

*. The error between a desired output*g

_{m}

*and the measured output*x

_{d}

*is*x

*=*e

*−*x

_{d}

*. Hence, the error vector is defined by*x

*can be written as*s

_{s}

*(*s

*,*x

*) =*t

K

^{T}

**, and its derivative is found as**E

*is determined by:*u

Assuming the disturbance of ** D** ≈ 0, then Eq. (39) can be rewritten as:

The relationship of Eq. (40) and OIT2FNN is expressed by

where, * f*(

_{00}

*) and*x

*(*g

_{00}

*) are the fuzzified functions of*x

*(*f

*) and*x

*(*g

*), respectively. The derivative of*x

**is expressed through Eqs. (40) and (41) as follows:**E

Define the minimum approximation error due to fuzzy approximation as follows.

Substituting functions of * f*(

_{00}

*),*x

*(*g

_{00}

*) and (43) into Eq. (42) yields the following equation.*x

Let _{1} established without the minimum approximation error * w*is defined as follows:

where, * γ*and

_{f}

*, respectively. The control*γ

_{g}

u

_{1}cannot use for control the system because of the error from the fuzzy approximation. To deal with this problem, a new robust compensator based on the inversely fuzzified value is suggested as follows:

where, * Γ*is a constant, and

**=**P

P

^{T}≥ 0 is the solution of the following Riccati-like equation.

where, * ρ*is the prescribed attenuation level,

**=**Q

Q

^{T}≥ 0,

*is consequent membership value of the OIT2FNN. When the value*ξ

_{z}

u

_{1}is the main controller to retain the stability of the system. From Eqs. (45) and (46), the final fuzzy control of the system is determined as follows:

Now, substituting Eq. (48) into (44) yields he following.

where,

The derivative of Eq. (50), and then substituting Eq. (25) into the derivative, the result is obtained as follows:

From Eq. (51), adaptation laws are established as follows:

Applying Eqs. (52) and (53), Eq. (51) can be written as follows:

where,

Now, the integration of (54) from * t* = 0 to

*=*t

*yields the following equation.*T

The value of * V*(

*) ≥ 0, and thus Eq. (55) is rewritten as follows:*T

where, _{1}, _{2} are the choosing parameters. Hence, the adjusted adaptation laws are redefined as follows:

Figure 2 presents a flow chart of the HAC-IFV showing the combination process of each controller with the adaptation laws.

## 4. Application to seat suspension system

### 4.1. Control results of the HAC-PP

In order to implement two adaptive controllers, principal parameters of the seat suspension and MR damper as shown in Figure 3 are given in [25]. And two different road profiles of random step wave road and regular bump road are adopted to emulate severe external disturbances as shown in Figure 4. The first excitation is collected from the real road, and the second excitation is used same as in [25]. The process of simulation is expressed as follows: The proposed control will be simulated following an objective trajectory, which is control of [25]. Then, the outputs of the proposed control and the objective will be used for calculating error. This error will be checked by desired prescribed performance. It is remarked that the desired prescribed performance is different from the applied prescribed performance which is shown in Eqs. (7)–(13). The parameters of both the desired and the applied prescribed performance are listed in Table 1. The damping force of the MR damper is designed 1000 N (±5%) at 2 A. The fuzzy model is established based on the online model with the centroid vector as shown in [25]. It is noted that two main variables for the fuzzy models are displacement and acceleration. The fuzzy models include 6 clusters, and then the outputs of fuzzy rules become also 6. The sigma value for Gaussian function of the fuzzy model is chosen as 0.4 [22, 25], and this value is not changed through the simulation. The values of the sliding surface [_{1}, _{2}] are chose by [1, 20] for both random step wave road and regular bump road. The constant value * Γ*of the Riccati-like equation is chosen by 10 for both roads. The constant

*is 500 and 5000 for regular bump road and the random step wave road, respectively. In addition, the matrix*c

_{s}

**of the Riccati-like equation is chosen as**Q

**= [−2 0; 0 − 2]. The constants**Q

μ

_{1},

μ

_{2},

μ

_{3}of adaptation laws are chosen as 10 for two road profiles. The values of

δ

_{1},

δ

_{2},

δ

_{3}are chosen by 0.1. In this simulation, the initial states for the dynamic states are used as

Parameter | Desired prescribed performance | Applied prescribed performance |
---|---|---|

Initial value (0) | 0.5 | 0.5 |

Infinity value _{∞} | 0.001 | 0.001 |

Exponential value | 1 | 0.00047 |

Figures 5–8 present control responses of the HAC-PP. It is clearly observed from Figures 5 and 6 that the initial excitation has been significantly reduced by activation the proposed adaptive controller in terms of both displacement and acceleration. In addition, it is seen that the proposed control well tracks the objective trajectory which directly indicates high performance of the prescribed performance of the sliding surface. Figure 7 presents the error of performance of the proposed adaptive controller which is always less than the boundary of the prescribed performance. These results mean that the application of the prescribed performance in design of the hybrid adaptive controller can improve the quality of control with high robustness against severe excitations.

### 4.2. Control results of the HAC-IFV

In simulation of the HAC-IFV, the values of the sliding surface [_{1}, _{2}] are chosen by [1, 1.10^{−5}]. The constant value * Γ*of the Riccati-like equation is chosen by 40, 10 for the regular bump road, the random step wave road, respectively. The constants

α

_{1},

α

_{2}of adaptation laws are chosen as 10 for all road profiles. The values of

*,*ε

_{f}

*of the expanded adaptation laws are chosen by 10 and the values of*ε

_{g}

δ

_{1},

δ

_{2}are chosen by 0.05. In this simulation, the initial states for the dynamic states are used as

k

_{1},

k

_{2}] are chosen as [1, 1.5] for random regular bump and [1, 5] for random step wave bump.

Figures 8–10 present control responses of the HAC-IFV. As similar to the HAC-PP, the initial excitations were remarkably reduced by applying the proposed controller. The displacements at the seat and driver positions are reduced resulting in the improvement of the ride comfort. In order to demonstrate a salient benefit of the proposed controller, its control response is compared obtained from the controller proposed in [17, 25]. It is clearly identified that the convergence time of the displacement of the proposed controller is 2 seconds for both excitations, while that is 15 seconds for the random step wave excitation, 6 seconds for regular bump excitation in [17, 25]. In Figure 8, the sliding surfaces of three controllers are shown. It is observed that the proposed control obtains stable motion much faster than the comparative controls at 0.1 second. It is noted here that the better control responses of the proposed controller comes from the inversely fuzzified values in given Eqs. (46)–(48). In Eq. (48), the independent of the inversely fuzzified value helps the controller to increase its robustness. This new exploration is the outstanding property of the proposed controller in the severe operation environment subjected to strong and random disturbances.

## 5. Concluding remarks

In this study, two new adaptive controllers were formulated and their effectiveness was validated by applying them to vibration control of a semi-active vehicle seat suspension system featuring MR damper. The first adaptive controller includes two sliding mode controls: one for initial states of the system and the other for prescribed performance associated with the parameters of the modified Riccati-like equation. By doing this way, the tracking performance is enhanced resulting in the improved control responses. The second adaptive controller was formulated on the basis of the inversely fuzzified value with the H-infinity control to minimize computational cost algorithm. Hence, by doing this way, the convergence time can be reduced resulting in high stability of the system subjected to severe external disturbances. It has been sown that the proposed two adaptive controllers can significantly reduce the excitation from the road profiles at both the seat and driver positions. In reality, this can enhance the ride comfort of the driver. Especially, the HAC-PP provides good tracking performance with the error in range of the defined boundary and the HAC-IFV can reduce the convergence time compared with two comparative adaptive controllers. It is finally remarked that the development of a new hybrid adaptive controller needs to be connected with desired control performances to appropriately select each control scheme.

## Acknowledgments

This research was funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 107.01-2017.28. The financial support is gratefully acknowledged.

## Declaration of conflicting interest

The authors declare that there is no conflict of interest.