1. Introduction
Fault-Tolerant Control (FTC) systems are methods developed to deal with possible faults caused by sensors and/or actuators in terms of reducing the impacts of these faults and keeping the performance of the system acceptable. These approaches have the ability to maintain a performance close to the desirable performance while preserving stability conditions in the presence of several types of faults caused by sensors and/or actuators. Over the past several decades, several FTC approaches have been introduced and investigated in the literature, aiming to reduce the faults’ impacts on the total stability and performance of systems [
1,
2,
3].
Generally, there are two major classes of FTC methods: active (AFTC) and passive (PFTC). In the PFTC approach, the FTC system cannot respond in real-time to all fault events since the controller’s structure and its parameters are formerly set and trained to only tolerate a specific set of faults [
4]. As a result, a PFTC approach cannot adjust its parameters to deal with the occurrence of faults, in real-time, that are not within the predefined set. On the other hand, the AFTC methods are mainly developed to deal with the occurrence of faults by accommodating faults in such a way that the controller updates its parameters to reduce the faults’ impacts on the control system [
5,
6]. As a result, the AFTC adjusts its structure or parameters based on the estimated faults obtained by the fault detection and identification (FDI) block, and the proposed AFTC design manipulates the estimated faults to maintain better performance in terms of stability, robustness, and fault tolerability [
7].
Linear Parameter-Varying (LPV) modeling approaches have received considerable attention over the past three decades in modeling and control areas. This is due to their efficiency in dealing with tractable mathematical descriptions for relatively complex nonlinear systems. One of the main types of LPV systems is the polytopic form, which allows the system to be described as a form of a convex combination of sub-models, and these sub-models are defined by the vertices of a convex polytope. As a result, the sub-models are combined to form the LPV model of the system by convex weighing functions, resulting in the global model of the system. Many nonlinear systems have been successfully represented by LPV models [
8,
9,
10].
Recently, many researchers have been attracted by the theory of solving nonlinear systems subject to faults based on the LPV representations, and many research papers have been published in the last few years. Li and Zhang [
11] investigated robust H-1 filtering for singular LPV systems subject to time-varying delays, while the states of the system were estimated but without any FDI purposes. Hamdi et al. [
12] proposed an approach to deal with polytopic models with unknown inputs and used a proportional integral observer for fault detection and isolation. Wang et al. [
13] proposed and investigated an adaptive fault diagnosis observer-based approach for regular LTI systems to deal with actuator faults. Although this proposed approach can be easily modified to deal with LPV systems, unfortunately, the methodology itself can only detect and estimate constant actuator faults. Lira et al. [
14] proposed a technique for LPV systems in which an observer-based approach is synthesized to estimate polymer electrolyte membrane (PEM) system states while varying parameters were planned with state variables. The proposed approach’s performance was also tested for several common fault scenarios that may appear during normal PEM system operation. Yadigar et al. [
15] proposed an approach based on a virtual actuator for fault-tolerant control of linear systems subject to time-varying additive actuator faults as well as an external disturbance. The proposed approach assumes that the structure of the nominal controller is known. Even though this approach was used to synthesize a controller for LTI systems, it can easily be modified to accommodate LPV systems. In [
16,
17], a H-infinity observer-based approach is designed to estimate states and sensor and/or actuator faults. Based on the synthesized observer, an output feedback fault-tolerant controller is established to guarantee the stability of a satellite attitude system. Recently, a novel observer-based approach for robust stable hybrid fault-tolerant predictive control was proposed by Zahaf et al. to investigate actuator faults in systems [
18]. Recently, the virtual actuator and virtual sensor-based techniques have attracted more attention in dealing with different types of faults for LPV systems. The authors in [
19] proposed an LPV virtual sensor design for FTC nonlinear systems, in which the nonlinear system is described as a LPV representation. In this proposed approach, the scheduling parameters for the controller, plant, and reconfiguration blocks are the same. In another approach proposed in [
20], the input and the output matrices of the LPV model are set to have dependent parameters while different virtual actuator (VA) and virtual sensor (VS) parameter designs are allowed depending on whether the faults are total or partial. In addition, several practitioners have considered the implementation of a bank of reconfiguration blocks [
1] technique in which each fault scenario is considered in the design. The proposed approach, based on predictive control, has shown good performance in handling actuator faults and maintaining robust trajectory tracking.
Generally, model predictive control (MPC) is considered to have been one of the most successful advanced control algorithms in many industrial applications over the past few decades [
21,
22,
23]. However, the FTC problem has not been fully investigated for model predictive control, and the FTC–MPC topic has recently become attractive to many researchers, with several algorithms being proposed to deal with actuator faults. Zou et al. [
24] proposed a model predictive fault-tolerant control (MPFTC) methodology optimized and adjusted by a Genetic algorithm. The proposed methodology was investigated for batch processes under disturbances and partial actuator faults. Zarch et al. [
25] proposed an actuator fault-tolerance controller based on nonlinear model predictive control systems. In this approach, the viability theory is used to synthesize the controller. Shi et al. [
26] proposed a robust constrained model predictive fault-tolerant control methodology for systems subject to uncertainties, unknown disturbances, interval time-varying delays, and partial actuator failures. Similarly, Sheikhbahaei et al. [
27] proposed a new methodology for explicit model predictive control of linear discrete-time systems under linear constraints and subject to unknown disturbances, uncertainties, and actuator faults. Zhang et al. [
28] proposed and presented a state-space model predictive fault-tolerant control methodology for batch processes subject to unknown disturbances and partial actuator faults. It is worth mentioning that most of the proposed approaches based on MPC algorithms for nonlinear discrete-time systems are based on the linearization of systems around operating points, which makes these approaches less accurate compared with nonlinear methodologies. In addition, it can be clearly observed that FTC-based MPC problems have not been fully investigated for LPV discrete-time systems subject to input constraints and sensor/actuator faults.
In this paper, an approach based on a virtual observer is proposed for discrete-time LPV systems subject to sensor faults. The proposed LPV model is used to simplify the complexity of the nonlinear system for the fault-tolerant predictive control (FTPC) and, at the same time, to maintain high accuracy for the controlling procedure. The proposed virtual observer is similar in theory to the observer presented in [
16,
17] but incorporated with the MPC approach. As a result, an extended state-space LPV model is used to construct a virtual observer, then the formulated virtual observer is used to propose a real observer. Knowing that the real observer is based on the parameters of the virtual observer, both proposed observers will be used to formulate the fault-tolerant predictive control strategy. The proposed methodology has the ability to reduce the influences of faults in the dynamics of the system and ensure the robust stability of the closed-loop system. In addition, sufficient stability conditions are proposed, in terms of linear matrix inequality (LMI) constraints, to ensure the robust stability of the overall closed-loop system containing system states, and the estimation of error dynamics. Moreover, the gains of both observers and the fault-tolerant controller are computed by solving the linear matrix inequalities (LMIs). This paper is structured as follows: the formulation of the LPV as well as fault representations are presented in
Section 2.
Section 3 presents both virtual and real observers to estimate the states and sensor faults at the same time. Based on the LPV model, an augmented robust observer fault-tolerant predictive control strategy is presented in
Section 4. The simulation results are illustrated in
Section 5. Finally,
Section 6 concludes the paper.
Notation. The superscripts and are used to denote the transposition and inverse of matrix , respectively. The symbols describe the controller gains matrices while the terms represent the observer gain matrices. denotes the th dimensional Euclidean space. and are the zero matrix and the identity matrix, respectively. In addition, the symbol represents the Euclidean norm of vector The term RMS(·) represents the root mean square function. The symbol * is used to describe the symmetric part in equations. The term Co{·} denotes a convex hull, where its element “·” refers to the vertex of the polyhedron. In addition, the matrices , , , , and are symmetric positive definite matrices. The expression (or ) denotes that the matrix is positive semi-definite (or positive definite, respectively). Finally, the symbols are the weighting matrices.
2. Convex Polytopic Model
Consider the following polytopic discrete-time model given by:
where
is the system state,
is the system output vector, and
is the system input vector. The system matrices can be written in the polytopic form as
,
, and
, where
and
.
The parameters are bounded in a predefined compact set , i.e., .
Consequently, the possible values of
are held within a polytope convex
of
vertices
. The images of the matrix
for each vertex
correspond to a set
, while the components of the set
are known as the extreme of the polytope [
15]. The system in Equation (1) can be represented by a convex interpolation of the vertices
and each vertex
of the polytope, which is obvious because of the linear dependence on
.
Note that the control action is set in order to satisfy the constraints below:
where
is the maximum value (upper limit) of the control input.
Since the proposed polytopic discrete-time system is subject to both sensor faults and external disturbances, the real output signal measured by the sensor may be altered when a sensor fault occurs. The output of the system is then modeled as follows:
where
is the vector of additive sensor faults.
Throughout this work, the following Lemmas and assumptions are needed to derive the main results.
Lemma 1 (Schur complements lemma [
29])
. Given any real matrices ,
, and with and , then Assumption 1. Sensor faults and the external disturbance are bounded, and there exists a positive scalar such that .
Remark 1. Note that Assumption 1 ensures that the increments of faults and the disturbances between two sampling time instants are bounded.
In this paper, only the sensor faults are considered and the model in Equation (1) can be rewritten as:
where
,
, and
denote the vector of the faulty control input, the vector of the faulty system state, and the faulty measurement output vector, respectively. Note that
is assumed to be an additive sensor fault.
3. Design of Observer-Based Control Law
The main objective of this work is to propose a virtual observer to improve the estimations and provide better accuracy for both fault observation and the system’s states.
First, an augmented system, which is proposed based on the system in Equation (4) and exposed to an external disturbance, is proposed for designing an observer that can estimate system state vector
and sensor faults
. The augmented system is defined as:
where
Multiplying
to both sides of (4), and by using fact (5), i.e.,
, where
and
, we have:
Let us define the following virtual observer:
where
is the estimation of the system state
and sensor fault signals
, and
is the observer gains, which will be obtained later.
To compute these gains, we will consider an augmented system with estimation error
of the observer in Equation (7) as
. Subtracting Equation (7) from Equation (6), the error dynamic equation is obtained as
Remark 2. The virtual observer cannot be realized since it includes unreachable information, i.e.,. Consequently, a real observer based on the virtual observer’s information is synthesized by eliminating the unknown terms.
Therefore, we define the auxiliary variable
. The real observer for the system in Equation (4) can be obtained by:
where a new auxiliary variable,
is introduced along with the matrices
,
, and
, and the gain matrices
are the observer matrices to be computed later as follows:
where
is a time-varying observer gain matrix that depends on
, which is defined as
4. Robust Fault-Tolerant Predictive Control
Next, a fault-tolerant predictive controller based on an observer is designed for LPV systems subject to disturbances and sensor faults. As the initial part of this work, an observer was constructed to estimate faulty signals of the system subject to an external disturbance. Next, based on the observer information a fault-tolerant predictive controller is proposed to guarantee stability and maintain control performance.
Next, the control law in Equation (10) is considered for compensating for the impacts of the faults in the closed-loop system and that rejects the disturbance’s influence.
where
is a time-varying gain matrix that depends on
, and can be defined as
.
Replacing Equation (10) in Equation (4), we obtain:
To obtain the predictive controller, an optimization problem that minimizes the following worst-case quadratic objective function in an infinite horizon is considered [
30]:
Subject to (6) and
, where
where
and
are known positive-definite weighting matrices.
Next, sufficient stability conditions, in terms of LMIs, are derived to ensure the stability of the overall closed-loop system in Equation (5).
Theorem 1. Consider the faulty LPV system in Equation (5). The state-feedback controller given by Equation (4) will robustly stabilize the system in Equation (6) subjected to external disturbances and sensor faults if there exist symmetric positive definite matrices,,,,and, and a positive scalarthat satisfies the following convex optimization problem:Subject towhere and .
Additionally, the control gain can be calculated by and the gain of the state observer is obtained by , .
Proof of Theorem 1. First, the following Lyapunov function is defined to obtain the stability conditions:
where
with
, and
with
. Furthermore,
and
and
For any
, suppose
satisfies the following stability constraint:
As it is assumed that summation is up to
, i.e.,
,
, summing from
to
produces:
By defining , an upper bound on the performance index is obtained as .
Hence, the first inequality of Equation (20) holds.
Next, we show that the second inequality of Equation (15) holds.
Since the inequality (18) implies that
strictly decreases as
goes to
and
from (18), we have
by the Schur complement, we obtain
Substituting and into the above inequality, and by applying the congruence transformation to the resulting inequality with Diag [1, , ], we conclude that Equation (15) holds.
Consider the input constraint in Equation (2)
Using the Schur complement, we obtain:
Substituting and followed by multiplying the right by and the left by , we obtain (16).
Now, (17) implies that:
where
and the matrices
and
are positives. Furthermore, since
, we have
By taking into account Equations (32) and (33), we can state that (17) implies
where
with
, and
with
. Furthermore,
and
, and by substituting
and
into Equation (29), pre-multiplying and post-multiplying by Diag [
], and using the Schur complement, we see that this is equivalent to:
Again, applying the Schur complement to Equation (35):
Then, multiplying the resulting inequality from the left by
and from the right by
and taking into account Equations (2) and (3), we have:
This will conclude the proof of Theorem 1. □
Next, a closed-loop system is implemented as a polytopic observer-based fault-tolerant predictive controller for LPV systems subject to sensor faults.
5. Results and Discussion
In this section, the proposed methodology is considered for a quarter-vehicle active suspension system and illustrated in
Figure 1 [
31,
32]. The quarter-vehicle suspension model parameters used in this study are displayed in
Table 1.
The mathematical model of the quarter-vehicle active suspension system is a nonlinear system described by the following nonlinear differential equations [
27]:
where
zus(
t) and
zs(
t) represent the vertical displacements of the unsprung mass and the vertical displacements of the sprung mass, respectively. The vector
zo(
t) represents the road displacements. The active force (control signal) of the suspension system is given by
u(
t). The sprung mass, which is also known as the car chassis, is given by
ms. The mass of the wheel assembly, also known as the unsprung mass, is given by
mus. The parameters
ks and
cs are the spring coefficient and the damping coefficient of the suspension system, respectively. In addition, the parameters
kus and
cus are the compressibility and the damping of the pneumatic tire, respectively. Finally,
Ns is the nonlinear stiffness property of the primary suspension. In this study, more attention was paid to the vertical behavior of the tire, which is useful in the design and the analysis of the comfort of the passengers and the road handling of the vehicle. The tire can be modeled by a simple spring
kus that is in parallel with a damper
cus. In addition, the stiffness of the spring is considered constant while its damping ratio is negligible.
The nonlinear system described in Equation (32) can be modified to have the following state-space representation:
Note that the tire deflection is the state variable . The rest of the state variables are the unsprung mass velocity , the suspension stroke , and the sprung mass velocity . The disturbance input that affects the system is represented as .
Note that the system matrix A can be expressed in polytopic form as:
where
are the vertex coefficients with a square signal
and
,
. The values of the vertex matrices
are
and
. In this work, we assume that the suspension deflection cannot exceed a maximum value, and can be constrained by the mechanical structure
. Note that the variable
is given by the expression
. As a result,
and
; that is,
and
, and
is a predefined constant.
Based on Equation (39), the faulty discrete state-space representation in Equation (38) can be expressed as a polytopic LPV system:
where the subscript
refers to the corresponding vertex of the polytope, and
is the vertex coefficient subject to
and
.
Considering the LPV model as described in Equation (1) with
N = 2 vertices, the nonlinear system in Equation (31) is described by Equation (35) with the following matrices:
Next, the MPC parameters are selected as follows. The weight matrices and of the cost function, in Equation (12), are set to and , respectively, where is the identity matrix with a proper dimension. A constraint is set on the active control force to limit the power of the hydraulic actuator. Thus, a hard constraint is imposed on the active suspension such that , and umax = 1500. In addition, the initial state values of the LPV system are set as follows: and for the augmented vector that includes initial fault . To test the robustness of the proposed method, the external disturbance, which is considered as a bump road profile, is expressed as . In the vehicle industry, the sensors are frequently subjected to strong vibrations that may add disturbances to the readings and cause sensor failure. Unfortunately, sensors with high vibration resistance are relatively expensive. In this study, we assumed that faults can be modeled in terms of vibration amplitude and frequency. So, the faults of the sensors injected into the closed-loop system are given by . Note that suspension deflection is bounded and should not exceed a maximum value defined by the mechanical structure. The suspension deflection was bounded as , where the maximum suspension deflection value is 0.04 m.
In this work, a LMI control toolbox in the MATLAB environment was used to solve the optimization problem in Equations (30)–(34). The results of the system states and their estimates and the fault with its estimate are illustrated in
Figure 2,
Figure 3,
Figure 4,
Figure 5 and
Figure 6, respectively.
The results indicate that the proposed control methodology can robustly improve the ride comfort compared with passive FTC methods. The tire deflection and sprung mass acceleration responses together with their estimated values are plotted in
Figure 2 and
Figure 3, respectively. The results in
Figure 2 show that the tire deflection response has been considerably reduced despite the fact that the system is subject to real-time sensor faults and disturbances. The same comment can be made about the sprung mass acceleration results shown in
Figure 3. In addition, the vertical acceleration of the vehicle’s vibration, presented in
Figure 3, is rapidly reduced by the fault-tolerant predictive control strategy while the actual value has been quickly and accurately estimated. As a result, the system may achieve lower values of body acceleration over time for the proposed active controller, which is better than the passive suspension systems proposed in [
21] and, as a consequence, this will ensure better ride comfort.
Figure 4 illustrates the suspension deflection, which is bounded by a maximum value of 0.04 m, where the system achieves lower values over time for the closed-loop systems. In addition, the proposed methodology has reasonably estimated the suspension deflection values, and the values have never exceeded the maximum value, which is better than the results presented in [
21].
Figure 5 illustrates the control signal
u(
t), i.e., the actual force applied by the actuator. Note that the proposed robust predictive methodology considers the failures at the sensor and reacts by increasing the force needed. As a result, the proposed robust controller has improved the vehicle’s ride comfort while satisfying all hard constraints, including both suspension deflection and actuator saturation constraints, even though the system is subject to an actuator fault.
Figure 6 and
Figure 7 present the fault estimation results, which show that the stability of the closed-loop system is guaranteed under external disturbances and sensor faults. The estimation results show that the proposed methodology has excellent estimation performance and efficiently tolerates sensor faults.
Table 2 summarises the RMS comparison between the proposed robust approach and the method discussed in [
21]. The comparison describes the RMS errors of the vehicle’s dynamic responses in terms of the suspension stroke, body acceleration, and tire dynamic load under random road excitation.
Table 2 presents the maximum values reached in body acceleration, suspension deflection, and tire deflection with the estimated RMS errors. The results suggest that the proposed robust predictive controller has better performance than the passive suspension system discussed in [
21] when the system is subject to sensor faults and external disturbances. In addition, the proposed FTPC structure has the ability to improve the ride quality and road-holding performance under different types of active suspension constraints and in the presense of external disturbances and sensor failures.
Note that the illustrated performance of the proposed fault predictive controller is due to the following:
The proposed predictive control structure is well adapted for LPV models such as LPV representations for automotive suspension systems. In addition, the external disturbances and faults in the output channel are already considered in the predictive control design procedures, which means that the controller is designed to meet both the stability of the vertical motion of the car’s body and the passenger’s isolation from vibrations caused by rough roads.
The proposed predictive controller has the ability to anticipate any errors caused by external disturbances and sensor failures, which eventually helps to reduce the impact of these errors while upper limits are imposed on the dynamic force of the motor. Note that the upper limits are used to avoid discomfort to passengers and damage to the vehicle’s components.