Recent Advances in Hybrid Dynamical SystemsView this Special Issue
LMI-Based Sliding Mode Observers for Incipient Faults Detection in Nonlinear System
This paper presents a diagnosis scheme based on a linear matrix inequality (LMI) algorithm for incipient faults in a nonlinear system class with unknown input disturbances. First, the nonlinear system is transformed into two subsystems, one of which is unrelated to the disturbances. Second, for the subsystem that is free from disturbances, a Luenberger observer is constructed; a sliding mode observer is then constructed for the subsystem which is subjected to disturbances, so that the effect of the unknown input disturbances is eliminated. Together, the entire system achieves both robustness to disturbances and sensitivity to incipient faults. Finally, the effectiveness and feasibility of the proposed method are verified through a numerical example using a single-link robotic arm.
An electronic system is structurally complex, involving a number of electronic components. It is therefore difficult to accurately analyze the relationship between the input and output for each component and the cause of any faults within, as each fault may exhibit diverse manifestations. Incipient faults are important in electronic systems, and some common manifestations include zero drift, reduced precision, delayed response, and equipment aging. If these potential incipient faults cannot be detected in time, the long-term stable operation of the devices will be affected and, more importantly, they will cause a decline in production capacity and an increase in the cost of production, or even accidents. Therefore, the study of incipient faults in electronic devices has practical significance.
The sliding mode variable structure control algorithm has the advantages of simplicity, high robustness, and high reliability and is extensively used in motion control. The sliding mode observer has attracted attention for fault diagnosis since the sliding mode variable structure control enables the system to be robust against unknown input disturbances, especially against external disturbances [1, 2]. Observer technology is one of the most important methods for fault diagnosis and detection, estimation and fault-tolerant control [3–5]. A robust fault diagnosis method for the sliding mode observer has been proposed . This scheme adopts the Walcoot-Zak sliding mode observer, which makes system state error robust against disturbances, and is maintained on the sliding mode surface. Much progress has been made in measuring the deviation in the trajectory of the motion of the system by using a sliding mode surface [5, 7, 8]. However, the current fault diagnosis based on sliding mode observers is mostly used for mutant or intermittent faults that have significant numerical changes. There have been few reports on slow faults with a small initial value. This is because the sliding mode control is essentially a type of continuous nonlinear control, and buffeting is inevitable for the system output, which causes the incipient fault signals to be submerged in the buffeting signals for a long period of time after their generation, thus they cannot be detected.
Chen and Chowdhury  have proposed a fault detection scheme for actuator’s incipient faults in linear systems. Building on the work of Chen and Chowdhury , this paper presents an LMI-based fault diagnosis scheme, which is designed for incipient fault detection in nonlinear systems with unknown input disturbances. The adoption of the LMI algorithm relaxes the selection criteria for key parameters in the design, making it easier to obtain these parameters.
2. System Descriptions
For a nonlinear systems with disturbances, where is unmeasurable state vector; is measurable input vector; is the measurable output vector; is the known nonlinear function; is unknown bounded nonlinear function, representing the unknown input disturbances and modeling errors, which are collectively referred to as input disturbances; is known disturbance distribution matrix. , , and are known matrices, where ,??,??,??.
Assumption 2.1. is observable.
Assumption 2.2. is full column rank and .
Assumption 2.3. is a bounded disturbance such that , where is a known function greater than 0.
3. Design of Fault Diagnosis Scheme
3.1. Coordinate Transformation
Coordinate transformation is used, under certain geometric conditions, to decouple the unknown input disturbance and faults . Assumption 2.2 implies that there are two transformation matrixes, and . The system (2.1) can be decomposed into the following two subsystems (the details of the transformation process are given elsewhere ): where , , , , , , ,??,??,??.
The transformation matrix is
3.2. Observer Design
Assumption 3.1. and are observable.
Assumption 3.1 implies that there are matrixes and , which enable and to be stable matrices:
There are also two Lyapunov equations: where ,??,??, and are all symmetric positive definite (SPD) matrices.
Assumption 3.2. The function satisfies the Lipschitz condition: where , and is the Lipschitz constant.
A Luenberger observer is constructed for subsystem (3.1): where superscript “” indicates estimate value.
A sliding mode variable structure observer is constructed for subsystem (3.2): where is the input signal of sliding mode variable structure, expressed as where is the matrix to be designed; is the scalar function to be designed, ; is a positive constant of small value.
Assumption 3.3. The matrices and have to be chosen such that .
Define ,??,?? as the state estimation errors and ,?? as the output estimation errors.
Lemma 3.4 (see ). If satisfies the Lipschitz condition, then there will be an SPD matrix that satisfies the following equation: where and is the Lipschitz constant.
Lemma 3.5 (Schur Complement ). For a given symmetric matrix , where . The following three conditions are equivalent:(1);(2),??;(3),??.
Proof. Consider the following Lyapunov function:
The derivative of the Lyapunov function with respect to time is
From Lemma 3.4 we find that
According to Lemma 3.5, when the following inequality is satisfied, then .
Since , , , the above LMI (3.19) can be transform into inequality (3.13), which is listed in Theorem 3.6.
This completes the proof.
The main advantage of the method presented in this paper is that the complex equations for obtaining parameters given by Chen and Chowdhury  are transformed into inequality (3.13), a standard LMI. Then the solution becomes much easier when the LMI toolbox is used in MATLAB.
3.3. Fault Diagnosis
When a fault occurs, the systems (3.1) and (3.2) are transformed into the following equations: where is a nonzero matrix, , , is a known fault distribution matrix. is an unknown bounded nonlinear function, which represents incipient faults of the system.
It can be concluded from Theorem 3.6 that, in spite of the unknown input disturbance, and still converge to zero field when there is no fault, or deviate from zero when the faults occur. In (3.22), which is an explicit of both faults and the unknown input disturbance, a sliding mode observer is designed for weakening the impact resulting from the disturbance. However, the inherent chattering phenomenon of slide mode variable structure control makes it difficult to distinguish the early incipient faults signals from the chattering signals until the faults develop to be serious enough. Thus, is unsuitable to be the faults detecting residual. By contrast, the in (3.21) is designed with no sliding mode observer in order to avoid the impact from chattering; therefore, can be used as the residual for incipient fault detecting. The simulation results in Section 4 will verify the above statements.
4. Simulation Study
Consider a nonlinear system, a single-link robotic arm with a revolute elastic joint rotating in a vertical plane whose motion equations are  where and are the link displacement and the rotor displacement, respectively. The link inertia , the motor rotor inertia , the elastic constant , the link mass , the gravity constant , the center of mass , and the viscous friction coefficients , ??are all positive constant parameters. The control is the torque delivered by the motor.
The state variables are chosen as ,??,??,??.
When handling different objects, the load was carried by the manipulator changes. In addition, the friction coefficient of the joint and other parameters also varies over time. All these factors are uniformly classified as unknown input disturbances, denoted by . Suppose a fault occurs to the manipulator, then the single-link robotic arm model with unknown input disturbances and faults can be expressed as the following fourth-order nonlinear state equation:
The selected manipulator parameters are ?Nm/rad, , ?Nm/(rad/s), ?Nm2, ?Nm2, ?kg, , ?m.
Corresponding to system (2.1), the parameter matrices for each equation are, respectively,
The selected transformation matrices and are, respectively,
Using the transformation matrices described above, system (4.2) is decomposed into the following two subsystems:
After the coordinate transformation, a reduced-order subsystem, as shown in (4.5), is decomposed from the system (4.2). The subsystem only contains fault ; it does not contain the unknown input disturbances . The impact of on subsystem (4.5) is delivered by the state variables and . Using the above observer, the impact of on the system is eliminated, and the unknown input disturbance and the fault are decoupled.
It is assumed that the incipient fault begins at time instant 30?s, and , when s. The input signal is and the simulation parameters are , .
With the above simulation parameters, we use three kinds of faults to verify the effectiveness of the proposed method. In the first case, an incipient fault is selected as , and a high frequency disturbance is used. As shown in Figure 1, the output estimation error from the sliding mode observer vibrates at a high frequency after the fault occurs, though its amplitude is small. This indicates that it is insensitive to the fault. However, the output estimation error increases rapidly after the fault occurs. Thus, the error is defined as the fault diagnosis residual , and a fault alarm can be sent immediately after a time of 30?s.
|(c) at 29.5–30.5 s|
Another low frequency disturbance is used to verify the proposed method. The simulation results are shown in Figure 2. It can be seen that output estimation error from the Luenberger observer is sensitive to the fault, but that from the sliding mode observer is still insensitive.
|(c) at 29.5–30.5 s|
In the second case, a low-frequency sinusoidal signal is selected to illustrate that the fault detection is sensitive to incipient faults, that is, . A high frequency disturbance and a low frequency disturbance are used, respectively. The associated simulation results in Figures 3 and 4 verify that the proposed approach can be applied to detect incipient fault rapidly.
In the third case, . A high frequency disturbance and a low frequency disturbance are used, respectively. The associated simulations are shown in Figures 5 and 6. Through the above verifications, we can see that the proposed method can effectively detect incipient faults.
This paper has introduced a fault diagnosis scheme which combines a sliding mode observer and a Luenberger observer for nonlinear systems. The proposed method makes full use of the complementarity between the sliding mode observer and the Luenberger observer and transforms the system, using a coordinate transformation, into two subsystems. The input estimation error from the Luenberger observer, which is constructed for the unrelated subsystem with unknown input disturbances, is defined as the fault detection residual in incipient fault diagnosis. A sliding mode observer is constructed to eliminate the impact of faults and disturbances in the other subsystem. An LMI-based approach is employed to design the observer, which makes it easy to obtain its parameters. The simulation results verify the validity of the proposed method.
This work is supported by the Natural Science Foundation of China (61273157 and 61104024), Hunan Province Education Department (12A040), and Construct Program of the Key Discipline in Hunan Province.
W. T. Chen and M. Saif, “Actuator fault isolation and estimation for uncertain nonlinear systems,” in Proceedings of the IEEE International Conference on Systems, Man and Cybernetics, vol. 3, pp. 2560–2565, Waikoloa, Hawaii, USA, October 2005.View at: Google Scholar
K. Zhang, B. Jiang, and V. Cocquempot, “Adaptive observer-based fast fault estimation,” International Journal of Control, Automation and Systems, vol. 6, no. 3, pp. 320–326, 2008.View at: Google Scholar
R. Sreedhar, B. Ferandez, and G. Y. Masada, “Robust fault detection in nonlinear systems using sliding mode observers,” in Proceedings of the IEEE Conference on Control Applications, pp. 715–721, Vancouver, Canada, September 1993.View at: Google Scholar
X. G. Yan and C. Edwards, “Robust sliding mode observer-based actuator fault detection and isolation for a class of nonlinear systems,” in Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference (CDC-ECC'05), pp. 987–992, Seville, Spain, December 2005.View at: Publisher Site | Google Scholar
X. G. Yan, S. K. Spurgeon, and Y. Orlov, “Sliding mode observer based control for a class of nonlinear time delay systems with delayed uncertainties,” in Proceedings of the 50th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC'11), pp. 7611–7616, Orlando, Fla, USA, December 2011.View at: Google Scholar
L. Yu, Robust Control-Linear Matrix Inequality Processing Method, TsingHua University Press, Beijing, China, 2002.