Open access peer-reviewed chapter

Fault Detection and Isolation of Nonlinear Systems with Generalized Hamiltonian Representation

By Luis Humberto Rodriguez-Alfaro, Efrain Alcorta-Garcia, Cornelio Posadas-Castillo and David Alejandro Diaz-Romero

Submitted: May 13th 2016Reviewed: February 23rd 2017Published: May 31st 2017

DOI: 10.5772/68084

Downloaded: 1016


The problem of fault diagnosis in a class of nonlinear system is considered. Systems that can be written in the so‐called Generalized Hamiltonian Representation (which is equivalent to an Euler‐Lagrange representation) are studied, and a model‐based observer approach for this class of systems is developed. The main advantage of the proposed approach is the facility to design the required observers, which take advantage of the system structure given by the Hamitonian representation. In order to show the proposed schema, a model of a permanent magnet synchronous machine is revised and the fault diagnosis schema presented. Simulation results confirm the effectivity of the proposed schema.


  • fault diagnosis
  • Hamiltonian systems
  • nonlinear systems
  • observers
  • fault isolation

1. Introduction

Safety operation and reliability of industrial processes are highly prized by the contemporary society. A key to achieve safety and reliability in industrial processes is through the use of diagnosis and fault‐tolerant control algorithms. Note that a fault is understood as a change of a parameter out of the tolerance limits. Physical systems are liable to potentially harmful fault events, which could cause a negative effect on the system functionality, as well as under‐performance. Faults can be originated by diverse reasons, for example, natural wear caused by common use, aging, use under stress conditions and so on. The importance of detecting and isolating the fault occurrence in a system lies in the possibility to reduce the maintenance and/or dead‐time for repairing on a production line.

There exist a lot of results related to fault diagnosis for linear systems, as it can be seen in the literature, for example, in Refs. [16], among others. For the case of nonlinear systems, there are also some available solutions based on diverse model structures, see, for example, Refs. [7, 8]. Fault diagnosis in nonlinear systems has been considered in Ref. [9], where the solution is based on a geometric approach, and the conditions of existence are not easily satisfied. Other approaches consider Lipschitz‐type nonlinear systems together with an observer‐based method [10, 11]. Ref. [12] is related to the problem of fault estimation for a class of switched nonlinear systems of neutral type, where the problem formulated as an Hfiltering is solved using a switched observer‐based fault estimator. In Ref. [13], the fault diagnosis is made for a class of bilinear systems considering only the case of faults on the actuator.

In Ref. [14], an unknown input observer (UIO) for a class of nonlinear state‐affine systems for fault diagnosis is proposed. By using sum‐of‐squares (SOS) theory and Lie geometry as the main tools, the rank constraint in the traditional UIO approach is relaxed and the design procedure simplified, especially for the case of nonlinear polynomial systems. In Ref. [15, 16], an approach to fault detection and isolation for the class of nonlinear systems with linear parameter varying (LPV) systems is shown. A different idea is to use a energy index in the diagnosis process, as in Ref. [17]. In Ref. [18], an algorithm for the diagnostics of nonlinear systems is presented where the solution is based on the estimation of the system parameters using the nonlinear response. The use of a bank of high‐order sliding mode observers has been proposed in Ref. [19].

From the above discussion, it is clear that even if some approaches are available to settle the fault diagnosis problem, in general there is no systematic way to design it (a model‐based or an observer‐based approach), because of the difficulty to design an observer for nonlinear systems even if the system is known. The available solutions consider a specific class of nonlinear systems, but each of these class of systems is more related to some mathematical (or system) properties and not necessarily to a wide class of systems from a practical point of view. Systems in Hamiltonian representation form represent a wide range of physical systems considering the relationship between Euler‐Lagrange and Hamiltonian systems [2022].

In this chapter, a solution to the problem of fault detection and isolation applying the observer‐based residual generation method is proposed. The class of nonlinear systems considered includes all systems, which admit a generalized Hamiltonian representation. The proposed solution begins with a mathematical nonlinear model of a system with faults. A nonlinear decoupling is applied to the faulty system in order to obtain a set of subsystems with sensibility to a particular fault or group of faults. Then, each subsystem is represented in a generalized Hamiltonian form, for which, a nonlinear observer is designed. Using the nonlinear observer, the residual generator is designed for each subsystem. One contribution of this work is the systematic way for residual generator design (an observer‐based approach with weak design requirements). Note that the observer‐based approach is guaranteed because of the Hamiltonian representation. Fault detection and isolation follow from the residual analysis. The structure of the Hamiltonian system representation is exploited to guarantee the residual existence for each subsystem. The approach is then applied to the model of a permanent magnet synchronous machine with additive faults. The faults are detected and isolated conveniently, showing the effectiveness of the proposed approach.


2. Generalized Hamiltonian representation of a system

Consider a nonlinear system described in general form as follows:


where xnis the state vector, uris the input vector, ymrepresents the output vector, and the function f:n×rnassociates with each value of xand ua corresponding ndimensional vector.

A special class of Generalized Hamiltonian representation is defined by Sira‐Ramierez and Cruz‐Hernandez [23],


where xndenotes the state vector, uris the input vector, Gn×ris a constant matrix, F(x)ndenotes a vector that contains the nonlinearities, ymdenotes the output vector, and Cm×nis a constant output matrix. Some nonlinear systems such as these described by Eq. (1) can be represented by Eq. (2) if satisfies the following conditions: There exists a smooth energy function H(x)that is positive definite in nand described by:


the column gradient vector denoted by H(x)x=Mxcan be obtained using Eq. (3), where Mn×nmust be a symmetric matrix constant and positive definitive, J(x)n×nmust be satisfied for all xn, and Sn×nis a constant symmetric matrix,

J(x)=JT(x), S=ST.E4

These conditions allow that a wide set of nonlinear systems can be brought to a generalized Hamiltonian representation, such as electromechanical systems, electric systems, mechanical systems, etc.

In the generalized Hamiltonian representation, the additive faults can be represented as in Eq. (5), where these appear as additional inputs (unknown inputs).


where N(Δf)nand Q(Δf)mrepresent the additive faults of the system.

3. Fault detection and isolation

In a general sense, if the fault diagnosis consists in the detection of a fault, then it is called fault detection(FD), and similarly, if the fault diagnosis consists in the detection and isolation of a fault, then it is called fault detection and isolation (FDI). The fault detection consists in determine the occurrence of faults in the functional units of the process, which leads to undesired behavior of the system, and the fault isolation consists in to classify the detected faults. The observer based fault diagnosis technique is a scheme of the model‐based fault diagnosis approach. In this technique, the idea is to replace the process model by an observer which estimates the fault‐free process outputs. The difference between the measured process variables and the estimated process variables defines the residual. The fault effect is contained in the measured process variables. Thus, a residual signal includes the fault effect. Ideally, if the residual is different from zero then a fault has occurred, otherwise the process is fault free. The residual generation allows to know the occurrence of faults, and the residual evaluation is necessary to extract the fault information. Figure 1 shows a common diagnosis scheme.

Figure 1.

Fault diagnosis scheme.

In this contribution, a fault detection and isolation approach to nonlinear systems that admit a generalized Hamiltonian representation is considered. The proposed approach follows the classical procedure of fault diagnosis: First, a fault decoupling in order to get subsystems with sensibility to a specific fault is developed. Second, an observer‐based residual generator for each subsystem is designed. Third, a residual analysis is performed to determine which functional unit has failed.

Figure 2 shows the proposed fault detection and isolation scheme, where Cnis the nominal control, ΣHis a system in Hamiltonian representation, and the diagnostic block contains the observer and the residual generator.

Figure 2.

Diagnostic scheme.

The fault decoupling consists in to define a transformation over the system in order to get a subsystem with sensibility to a fault or a set of faults, and this subsystem is coupled with a fault and decoupled from the rest of possible faults. There are some works on the analysis and synthesis of these transformations, see, for example, Refs. [9, 24, 25].

For the case of the generalized Hamiltonian representation with faults, Eq. (5) considers the following nonlinear transformation


It is required that


and the transformation T(x)be selected in such a way that the resulting transformed system has the desired fault sensibility, that is, suppose N(Δf)=[n1(Δf) n2(Δf)  nl(Δf)]where N¯(Δf)represents the columns associated with the faults that requires no to affect a specific subsystem and N¯¯(Δf)are the columns related to the faults that are required to affect the subsystem. With T(x)xN¯(Δf)=0and T(x)xN¯¯(Δf)0. In Ref. [24], it can be found details about the existence of this transformation.

Assumption 1.Consider the systemEq. (1)in generalized Hamiltonian representation with faults as inEq. (2)as well as the nonlinear transformationT(x)satisfying decoupling requirements. Also, the transformed system (decoupled) can be represented in the Hamiltonian form given byEq. (2).

For some examples, at least, the assumption is satisfied and consequently, a systematic way to fault isolation is obtained. At the moment, we do not have a result on the characterization of the class of systems for which the assumption is satisfied. The resulting decoupled system is represented in a Hamiltonian form.


where ζnζdenotes the state vector, J¯(ζ)nζ×nζ, S¯nζ×nζ, urζis the input vector, G¯nζ×rζis a constant matrix, F¯(ζ)nζdenotes a vector that contains the nonlinearities, yζmζdenotes the output vector, and C¯mζ×nζis a constant output matrix.

After a subsystem has been determined, the next step is to design an observer for each subsystem. From the decoupled subsystem in the generalized Hamiltonian representation Eq. (9), an observer can be designed as follows [26]:


where Kn×mis the observer gain, ζ^nis the estimated state, y^ζmis the estimated output calculated in terms of ζ^, H(ζ^)ζ^=Mζ^is the gradient vector with Mn×nas a symmetric positive definite matrix.

For this observer, the conditions design is described in the following Theorem:

Theorem 1.The statexof the nonlinear system in the generalized Hamiltonian representationEq. (9)can be globally, exponentially, asymptotically estimated by the observerEq. (10), if the pair (C,S) is observable or at least detectable and the matrix


is negative definite. WithΠ=12[MF(r)x+(F(r)x)TM]andρis a vector such thatρ(x,ζ^).


The proof of Theorem 1 is fully defined and explained in Ref. [26]. Then, for the decoupled system, a residual generator is defined as follows

Theorem 2.For the decoupled nominal system (Eq. (9)withN¯¯=0). The system


is a directional residual generator if the pair(C¯, S¯)is observable or at least detectable and the matrix


is negative definite. WithΠ=12[MF¯(r)ζ+(F¯(r)ζ)TM]andρis a vector such thatρ(ζ,ζ^).



The proof of Theorem 2 is a consequence of the proof of Theorem 1.


4. Application example

In this section, the results to apply in the permanent magnet synchronous machine (PMSM) the proposed approach for fault detection and isolation are presented. The closed loop system is used in the fault diagnosis analysis where any specific control law is used.

The PMSM mathematical model in the stationary reference frame dq0(direct‐quadrature‐zero axes) is taken from Ref. [27] and is described by:


where x=[id iq ω]T, Gm=diag[1L1L1Jm], um=diag[uduqτL]and


where Bis the viscous friction coefficient, Ris the stator resistance, Lis the inductance, Φis the flux linkage, Pis the pole pairs, idand iqare the electric currents on the direct and quadrature axis, respectively, udand uqare the voltages on the direct and quadrature axes, respectively, ωis the rotor speed, Jmis the rotor inertia, and τLis the load torque.

In the fault diagnosis analysis, it is considered that the system is operating in nominal conditions, which implies that the system is in closed loop with any controller. In this case, a back‐stepping nonlinear control [22] is used in the PMSM.

In order to obtain the Hamiltonian representation Eq. (2) of the PMSM described by Eq. (15), a Hamiltonian energy function is defined as follows:


with a gradient vector

H(x)x=Mx  M=diag[1L2Jm3],E18

where Mis a symmetric, positive definite and constant matrix so that the Hamiltonian representation of the PMSM is as follows:


where x=[id iq ω]T, u=[ud uq τL]T,

J(x)=[000003PΦ2JmL03PΦ2JmL0],S=[RL000RL20003B2Jm2], C=[10001/L00032Jm], F(x)=[PωiqPωid0], G=[1L0001L0001Jm].E105

Solving the Hamiltonian representation Eq. (19) for each of the state equations, the same model described by Eq. (15) is obtained, so that the Hamiltonian representation is correct since it fulfills the conditions Eq. (4).

An intermittent connection, signal lost or signal offset are some of the sensor faults prone to occur in electrical machines [28], the control objective is affected mainly by first and second faults. The nominal value of the load torque is known, an unknown change in this parameter is considered as an additive fault. The PMSM may occur faults on elements such as sensors, actuators and components. The following additive faults are considered in this contribution: Δωis a fault in the speed sensor, fais a fault in the control input, and ΔτLis an unknown change in the load torque.

When these faults are considered, the Hamiltonian representation of the PMSM is as follows:


where x, S, J(x), F(x), Gand uare the same as in the nominal case when there are no faults and

N(Df)=[fafaΔτL], Q(Δf)=[00Δω].E106

Once defined the mathematical model of the PMSM with faults, the fault decoupling is done for each fault presented in the system. From this fault, decoupling analysis is obtained subsystems with sensibility to a particular fault and without sensibility for the rest.

Subsystem sensitive to the control input faultfa:For this subsystem, decoupling the output y3is not used to avoid the sensor fault effect. Considering the first two equations and the outputs y1and y2of the faulty system Eq. (21) a subsystem sensitive to the actuator fault is obtained, as follows:


solving Eq. (23) for x3


now replacing x3in Eq. (22)


multiplying Eq. (26) by (ϕL+x1)and solving for ϕLx˙1+x1x˙1+x2x˙2=υ˙1,


Eqs. (27) and (28) are the subsystem 1with sensitivity to the control input fault fa, where x1and x2are quantities available in measurable outputs y1and y2, respectively.

Subsystem sensitive to the load torque faultΔτL:once more the output y3is not used to avoid sensitivity to the sensor fault. Subtracting Eq. (22) to Eq. (23)


if a new state υ˙2=x˙1x˙2is defined, Eq. (29) becomes Eq. (30), this equation and the third equation of Eq. (21) define the subsystem 2,


where x2is available in the measurable output y2.

Subsystem sensitive to the sensor faultΔω:since this subsystem must be sensitive to the sensor fault, the output y3is used. Using the transformed state υ˙3=x˙1x˙2, the subsystem 3 is obtained with sensibility to the sensor fault:


where x3and x2are quantities available in the measurable outputs y1and y2, respectively.

Once decoupled subsystems were obtained, for the residual generator design an observer for each one for each of the decoupled subsystem is designed.

For decoupled subsystems sensitive to faand Δω, the observer design using the proposed approach in Ref. [26] coincides with a Luenberger observer [29, 30], but, however, this does not apply for decoupled subsystems sensitive to ΔτL.

The observer design and the residual generator for the decoupled subsystem sensitive to ΔτLare presented. The decoupled subsystem sensitive to ΔτLcan be expressed as follows:


Which can be written in the form Eq. (9) with


where xa=[υ2 x3]Twhose gradient vector is defined as follows

H(xa)xa=[RLυ2, BJmx3]TM=diag[RL, BJm],E37

and with

J(xa)=[0JmPϕ2BLJmPϕ2BL0], S=[1JmPϕ2BLJmPϕ2BL1], u=[uduqΔτL], F(xa)=[Px3(v2+2y2)3Pϕ2Jmy2τLJm], G=[1L001], C=[LRJmB].E107

For this case, the pair (C,S)is observable, and thus, there exists a matrix Kthat satisfies both the requirements of Theorem 1 for the observer design and the requirements of Theorem 2 for the residual generator design, and thus the observer is as follows


and the directional residual generator for the decoupled subsystem sensitive to ΔτLis given by


For decoupled subsystem sensitive to control input faEqs. (27) and (28), the observer and its directional residual generator are as follows:

υ^˙1=RLυ1R2L(y22+y12)+ud(y1L+ϕL2)+ uqy2L+L1(υ1υ^1),E40

Finally, for the decoupled subsystem sensitive to the sensor fault ΔωEq. (33), the observer and its directional residual generator are as follows:


To summarize the fault sensitivity results of each residual (associated to each subsystem), see Table 1.

Residual (subsystem)Fault
Actuator fault (fa)
Sensor fault (Δω)
Change of charge (ΔτL)

Table 1.

Fault incidence table.

Where the fault affecting the residual is indicated with , and the symbol means that there is no connection between the fault and the corresponding residual.

As can be appreciated from Table 1, there is a one‐to‐one relationship between faults and residuals so that perfect decoupling has been attached. One nice thing of perfect decoupling is that the occurrence of faults can be detected and isolated without problems.

The following results were obtained by computer simulation. Table 2 shows the considered faults. About 10%of the nominal value of each variable is the fault magnitude considered, where 34.62is the nominal value of the control input on the stationary reference frame dq0, 100 rad/secis the nominal value of the angular speed, and 1.4 Nmis the nominal value of the load torque.

CaseFaultFault interval(sec.)Magnitude
2fa[2,2 . 5]3 . 462
3Δω[3,3 . 5]10
4ΔτL[4,4 . 5]0 . 14
5fa, Δω, ΔτL[3,3 . 5], [4,4 . 5], [2,2 . 5]3 . 462, 10, 0.14

Table 2.

Fault cases.

Table 3 shows the PMSM parameters, which were taken from Ref. [22].

ParameterNumerical value
Resistance (R)1.6 Ω
Rotor inertia (Jm)76.5 × 10−6 kg m2
Viscous friction coefficient (B)4 × 10−6 Nm/rev/min
Flux linkage (Φ)0.29 Nm/A
Inductance (L)9.4 H
Load torque (τL)1.4 Nm
Pole pairs (P)1

Table 3.

PMSM parameters.

Figure 3 shows the evolution of PMSM states in the time (the time scale is given in seconds), that is, the current in the direct axis id(t), the current in the quadrature axis iq(t)and the angular velocity ω(t)in nominal conditions (without faults). Actually, this figure represents the response of the PMSM with nominal parameters.

Figure 3.

Nominal states.

The evolution of the residuals when an actuator fault faoccurs is depicted in Figure 4, where both residuals 2and 3are equal to zero since these are insensitive to the fault fa, while residual 1is different from zero indicating the sensitivity to control input fault fa. Note that the magnitude of the two first residuals represents deviations between nominal and measurement currents with respect to the time (time is given in seconds). The third residual represents the deviation of the nominal and measured angular velocity of the PMSM rotor.

Figure 4.

Residual sensitive to the actuator faultfa.

Figure 5 shows the residuals evolution when the sensor fault occurs, where residuals 1and 3are zero at all time due to its insensitivity to this fault, while residual 2differs from zero due to its sensitiveness to this fault.

Figure 5.

Residual sensitive to sensor faultΔω.

Figure 6 shows the residuals evolution when the load torque fault occurs, where residuals 1and 2are null and residual 3is different from zero, indicating the sensitivity to load torque fault. Note that the fault magnitude is of 0.14; however, the residual becomes a value around 3, that is, the effect of the fault is not directly the magnitude of this. Extra work is required in the design of the observer‐based residual in order to get at the residual a more approximated value of the fault magnitude.

Figure 6.

Residual sensitive toΔTL.

Figure 7 shows the case when all three faults occur, even if not at the same time. The effect of the faults is manifested in the correct residual. It means that the problem of fault isolation in multiple faults can be carried out effectively.

Figure 7.

Residuals occurring simultaneously atfa,ΔωandΔTL.

As a final note, it can be appreciated a minimum transient at the beginning of all residuals signals, and this transient does not affect the fault detection and isolation process.

5. Concluding remarks

Fault diagnosis for a wide class of nonlinear systems, the class of systems that admit a Hamiltonian representation, has been considered. An observer‐based solution with weak existence conditions for the fault diagnosis has been proposed, and this approach allows the detection and isolation of additive faults.

The proposed procedure follows the traditional way, namely: First, a decoupling methodology is applied to systems with Hamiltonian representation in order to obtain subsystems that preserve the Hamiltonian structure. Observer‐based residual generators are designed for each subsystem so that each residual generator is sensible to a fault (or to a specific group of faults). The residual has the property of remain close to zero (or under a threshold value) if no fault is present in the system and non zero (or greater than a threshold value) when a fault affects the system. The proposed approach solves the fault isolation problem, and it permits a systematic design of the required residual generators. In contrast with other methodologies, for systems with Hamiltonian representation, an easy way to design an observer has been introduced. In addition, a wide set of nonlinear systems can be represented in the Hamiltonian structure, making the proposed solution widely applicable.

The proposed methodology has been applied to a synchronous machine, showing that, using the proposed approach, it is possible to detect and isolate additive faults in scenarios such as a fault in the control input, a change in the load torque as well as a fault in the angular velocity sensor.

Future work includes the study of multiplicative fault type.

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Luis Humberto Rodriguez-Alfaro, Efrain Alcorta-Garcia, Cornelio Posadas-Castillo and David Alejandro Diaz-Romero (May 31st 2017). Fault Detection and Isolation of Nonlinear Systems with Generalized Hamiltonian Representation, Fault Diagnosis and Detection, Mustafa Demetgul and Muhammet Ünal, IntechOpen, DOI: 10.5772/68084. Available from:

chapter statistics

1016total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Process Monitoring Using Data-Based Fault Detection Techniques: Comparative Studies

By Mohammed Ziyan Sheriff, Chiranjivi Botre, Majdi Mansouri, Hazem Nounou, Mohamed Nounou and Mohammad Nazmul Karim

Related Book

First chapter

Sliding Mode Control and Fuzzy Sliding Mode Control for DC-DC Converters

By Kamel Ben Saad, Abdelaziz Sahbani and Mohamed Benrejeb

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us