Open access peer-reviewed chapter

Robust Diagnosis by Observer Using the Bond Graph Approach

Written By

Abderrahmène Sallami

Submitted: January 26th, 2018 Reviewed: May 22nd, 2018 Published: December 12th, 2018

DOI: 10.5772/intechopen.79046

Chapter metrics overview

869 Chapter Downloads

View Full Metrics


In this chapter, we have proposed a Luenberger observer linear systems diagnostic technique using the bond graph. Indeed, an observer can reconstruct or estimate the current state of a real system using the available measurements, without prior knowledge of the initial conditions. In addition, it allows to estimate the nonmeasurable states of a system. The design of the observer is carried out using graphical methods from the structural properties of the model bond graph becomes simple and practical to build. We presented the bond graph approach for the construction of a full-order observer and proposed a new BG-based observer diagnostic method. Subsequently, we presented the uncertain parameter systems modeled by the bond graph approach, and we also proposed a new method for diagnosing systems with uncertain parameters by Luenberger observer. In the last part of this chapter, we developed and proposed an observer bench diagnosis technique (BG-DOS/BG-GOS) to detect and locate defects.


  • robust diagnosis
  • bond graph
  • Luenberger observer
  • DC motor
  • generation of robust residues

1. Introduction

Nowadays, the engineering sciences rely heavily on the estimation of the state of the systems. Indeed, the complete knowledge of the state of a system is often necessary for the elaboration of a control law or the setting up of a strategy of monitoring or diagnosis. In practice, the state of a system is not always available and the input and output signals are the only quantities accessible by measurement. The most widespread solution to this problem consists in coupling to the system another auxiliary system, called estimator or observer of state. The observer provides an estimate of the state of the system based on its model and measurements of its inputs and outputs. The observer conventionally used, within the framework of linear systems, is said to have proportional gain or Luenberger. In recent years, observer diagnosis is used to estimate shareholder defects and sensor defects [1].

In this chapter, we show how the leap graph model can be used for modeling, simulation and residual determination by Luenberger observers for diagnosis system.


2. Observer diagnosis

2.1. Observational diagnosis using the analytical model

The principle of diagnosis consists in estimating, by appropriate techniques, all the components of the state vector or, more generally, the output of the process, using the error of estimation as a residue [2].

This operation is carried out by means of a proportional observer. The functional diagram of such a method is given in Figure 1.

Figure 1.

Observational diagnosis.

The residue equations are defined as follows:

State estimation residue:


Exit estimate residue:


2.1.1. Residual in the case of normal operation

We consider a linear system whose equations of states are defined as follows:


The proportional state observer is defined by equations:


The state estimation or state reconstruction residual is by definition:


For a proportional observer, the evolution of the residue r(t) is:


The Laplace transform of Eq. (6) is written:


The output reconstruction residue is defined by:


Hence using the previous result, Eq. (8) is written:


2.1.2. Residual with a sensor failure

We consider again the same linear system in which the measurements made are subjected to a defect which we denote by fc(t) of unknown amplitude and appearing at an unknown time (Figure 2):

Figure 2.

Proportional observer with sensor fault.

The equations of states become:


In order to diagnose the state of operation of the system, the states reconstruct or must be sensitive to this bias and must highlight this sensor defect, isolate it and possibly quantify it. The residue r(t) is:


The equation is written using the Laplace transformation:


To see the error of reconstruction of output according to the biases of the sensors, we do not take into account the initial conditions of or:


2.1.3. Residual with an actuator failure

In the case where the control applied to the system is misinterpreted by the actuator, the latter introduces a fault which we denote by fa(t)(Figure 3).

Figure 3.

Proportional observer with actuator fault.

The equation describing the faulty system is:


For zero initial conditions, the expression of the state reconstruction residue is of the form:


The output reconstruction residue is defined by equation:


2.1.4. Residual with sensor and actuator failures

In this case, we have the addition of two defects at the system level:

  • Fault at sensor fc(t);

  • Fault at actuator fa(t).

The equations of the system are then written as follows:


Applying the principle of superposition, since the system is linear, we can determine the output reconstruction residue which is the sum of the residuals caused by the sensor alone and the actuator alone. The expression of this residue is as follows:


This residue is sensitive to sensor and actuator defects.

2.2. Observational diagnosis using the bond graph model

2.2.1. Construction of a Luenberger observer based on model bond graph

To construct the observers, one must check the observability of the system. From a Bond Graph perspective, proposed by [3], a system modeled by leap graph is structurally observable, if the following two conditions are met:

  • First condition: There is at least one causal path linking a sensor to each dynamic element I or C in the integral causality when the hop graph is the preferred integral causal model.

  • Second condition: All elements I or C admitting a derivative causality when we put the leap graph model in derivative causality, and dualize the sensors.

Figure 4 presents, respectively, the proportional observer using the bond graph for elements I and C [4].

Figure 4.

Construction of a proportional observer case of: (a) element I and (b) element C.

2.2.2. Proportional observer diagnosis using the bond graph model

The observer provides an estimate of the state of the system based on its model and measurements of its inputs and outputs. The observer conventionally used in linear systems is said to have Proportional (P) or Luenberger gain [5].

Consider a continuous system described by the equation of state using the hop graph represented in Figure 5:

Figure 5.

Schema of the continuous system described by bond graph.

The state observer of proportional type is represented by Figure 6. The equation of state is of the following form:


Figure 6.

Structure of a Luenberger observer using graph graph.

2.2.3. Generation of residues

The generation of residues from a proportional observer using the bond graph model is summarized by the following steps:

Verify that the bond graph model of the system is structurally observable, if so, then continue the next steps;

Construction of the observer using the graph;

The symbolic expression of residue is deduced from equation:


After calculation, the residue is in the form:


2.3. Determination of the gain of the Luenberger observer

2.3.1. Determination of the gain of the Luenberger observer by the analytical method

To guarantee the asymptotic convergence toward zero of the estimation error, choose the gain K to stabilize the matrix (A-KC).

The gain K can be determined by two methods:

  • Placement of the poles: This method makes it possible to have a stable observer and allows playing on the criterion of rapidity of the convergence of the observer by choosing eigenvalues greater or less in absolute values.

  • LYAPUNOV criterion: This method guarantees the stability of the observer but gives no idea about the speed of convergence of the observer since it does not allow to choose the eigenvalues.

2.3.2. Determination of the gain of the Luenberger observer by the bond graph method

First, we will compute the coefficients of the characteristic polynomial of the linear system:


The calculation of the gain of the observer by leap graph is based on the Rahmani theorem, [6].

Theorem:The value of each coefficient (ai) of the characteristic polynomial PA isequal to the total gain of families of causation cycles of order i of the model graph graph:

The gain of each affected family of causality cycles must be multiplied by (−1) d if the family is constituted by disjoint causality cycles.

2.4. Example

Consider the electrical system of a DC motor and its model bond graph given in Figure 7. On this system, we will detect and locate faults at the flow sensors (speed sensor Df1 and current sensor Df2).

Figure 7.

(a): Direct current motor, (b): Bond model of the DC motor.

2.4.1. State equation

The equation of state of the DC motor is as follows:


With R: armature resistance, L: armature inductance, m: torque coefficient, J: moment of inertia, b: coefficient of friction.

R = 1 Ω; L = 5 mH; b = 10−4 Nm/rd. S−1; J = 10−3 Kg.m2, m = 0.2 Nm/A.


2.4.2. Determination of the gain of an observer Luenberger by the analytical method

It is desired that the poles have the following values s1 = −160 and s1 = −75, then the characteristic polynomial:

Let us calculate the polynomials characteristic of the observer P1(A-K1C1)and P2(A-K2C2)

  • For the characteristic polynomial P1(A-K1C1)

  • For the characteristic polynomial P2(A-K2C2)


2.4.3. Determination the gain of a Luenberger observer by bond graph model

Using the theorem of Rahmani [6]:

  • The coefficient a1is equal to the total gain of families of first order causality cycles in the jump graph model of Figure 8.

  • The coefficient a2is equal to the total gain of second-order families of the causation cycles in the jump graph model of Figure 8.

Figure 8.

Schema of the LFT principle.

To calculate the gains k11and k12, it is necessary to use Table 1 that indicates the causality cycles in the observer’s bond graph model.

  • For the BGO1 model

Table 1.

Causation cycles in the bond graph model of observers BGO1 and BGO2.

The coefficient a1 is equal to the sum of the gains of the first order:

a1 = G1 + G2 + G3


So k11 = 0.1745.

To calculate k12, it is necessary to determine the coefficient of a2(a2equal to the sum of the gains of the second order):

  • For the BGO2 model

The coefficient a1 is equal to the sum of the gains of the first order:

a1 = G1 + G2 + G4


So k22 = 0.0349

To calculate k21, it is necessary to determine the coefficient of a2(a2equal to the sum of the gains of the second order):


Table 2 represents the total earnings of i-order causal cycle families in the first and second observer graph (BGO1 and BGO2).

Causal of order iFamily of cycles
Total gain in BGO1Total gain in BGO2
1200.1 + 200 k11200.1 + 1000 k22
28020 + 20 k11 + 4000 k128020 + 40,000 k21 + 200,000 k22

Table 2.

Total gains of families of causal cycles in the bond graph model of observers BGO1 and BGO2.


3. Robust diagnosis using bond graph model

3.1. Introduction

The diagnosis of industrial systems with uncertain parameters has been studied by several researchers in recent years. Indeed, the method most used is the method of the form linear fractional transformations (LFT), this method offers several advantages such as the flexibility of the diagnosis point of view that allows to model all the uncertainties. But unfortunately, the transition to the LFT form is not always possible (e.g., the nonlinear state models) because the separation of the nominal part from its uncertain part is very difficult (even impossible). Djeziri [7] opted for bond graph modeling using the LFT form method, using a single leap graph model in LFT form generates residuals and adaptive thresholds of normal operation with perfect separation.

3.2. Building a bond graph

Two methods are proposed by Dauphin-Tanguy [8] and Sié Kam [9] to construct parametric uncertainties by BG. The first is to represent uncertainty on a leap graph element as another element of the same type, causally related to the nominal element or the rest of the model. These uncertainties are kept in derived causality when the model is in integral preferential causality so as not to modify the order of the model. The second method is linear fractional transformation (LFT) introduced on mathematical models by Redheffer [10].

3.3. Representation LFT

Linear fractional transformations (LFTs) are widely used in the modeling of uncertain systems. The universality of LFT is due to the fact that any rational expression can be written in this form according to Alazard [11]. This form of representation is widely used for the synthesis of the control laws of uncertain systems using the principle of μ-analysis. It consists of separating the nominal part of a model from its uncertain part as illustrated in Figure 8.

The nominal values are grouped in an augmented matrix denoted M, supposed to be proper, and the uncertainties whatever their type (structured and unstructured parametric uncertainties, modeling uncertainties, measurement noises, etc.) are combined in a matrix Δ of structure diagonal. In the linear case, this standard form leads to a state representation of the form (Figure 8):

  • x Rn: the state vector of the system;

  • uRm: the vector combining the control inputs of the system;

  • yRp: the vector grouping the measured outputs of the system;

  • wRland zRl: group respectively the inputs and the auxiliary outputs;

  • n, m, l and p are positive integers.

The matrices (A, B1, B2, C1, C2, D11, D12, D21, D22) are matrices of suitable dimensions.

3.4. Formatting assumptions LFT

Forming LFT requires that the model be clean and observable [12]. The bond graph methodology allows by causal manipulations to check these properties directly on the bond graph model.

Property 1: A leap graph is unique if and only if it contains no derivative causal dynamic component when it is a preferential integral causality, and reciprocally.

Property 2: A leap graph is structurally observable if and only if the following conditions are met:

  • On the bond graph model in integral causality, there exists a causal path between all dynamic elements I and C in integral causality and a detector De or Df;

  • All dynamic elements I and C admit causality derived on the leap graph model in preferential derived causality. If dynamic elements I or C remain in integral causality, the dualization of detectors De and Df must make it possible to put them into derivative causality.

3.5. Modeling of BG elements by LFT

The modeling of linear systems with uncertain parameters has been developed in Djeziri’s thesis [13], the modeling of uncertain graph leap elements (R, I, C, TF and GY) has been determined. We will therefore limit ourselves in this section to show the two methods of modeling uncertain BG elements, as well as the advantages of BG-LFT for robust diagnosis.

3.5.1. BG element with a multiplicative uncertainty

The introduction of a multiplicative uncertainty on, for example, the element R in causality resistance gives Eq. (25):

  • Rn: The nominal value of the R element;

  • aR: The multiplicative uncertainty on the parameter;

  • eRand fR: represent the force and the flux in the element R, respectively;

  • enand einc: represent respectively the effort provided by the nominal parameter and the effort introduced by the multiplicative uncertainty.

The bond graph model equivalent to the mathematical model of Eq. (24) is given in Figure 9.

Figure 9.

(a): Model BG-LFT of an element R causality resistance with multiplicative uncertainty, (b): Model BG-LFT of an element R with conductance causation with multiplicative uncertainty.

3.5.2. Construction of a BG-LFT model

The complete BG-LFT model is represented by the diagram in Figure 10.

Figure 10.

Representation in form BG-LFT.

3.5.3. Rugged diagnosis by ARRs

The generation of robust analytical redundancy relationships (ARRs) from a clean, observable and overdetermined leap graph is summarized by the following steps:

Step 1: Checking the state of the coupling on the model graph deterministic preferential derived causality; if the system is overdetermined, then proceed with the following steps;

Step 2: The leap graph template is set to LFT;

Step 3: The symbolic expression of the RRA is deduced from the equations at the junctions. This first form will be expressed by:

For junction0:bifinc+Sf+wi̇=0E25
For junction1:bieinc+Se+wi̇=0E26

With Sḟthe sum of the flux sources linked to the junction 0, the Sėsum of the force sources linked to the junction 1, b = ± 1 according to whether the half-arrow enters or leaves the junction, einand finare the unknown variables .

Is the wi̇sum of the modulated inputs corresponding to the uncertainties on the elements related to the junction.

Step 4: Unknown variables are eliminated by traversing causal paths between detectors or sources and unknown variables;

Step 5: After removing the unknown variables, the uncertain RRAs are in the form:


TFn and GYn are, respectively, the nominal values of the modules of the elements TF and GY, also Rn, Cn and In are the nominal values of the elements R, C and I.

3.5.4. Robust diagnosis by Luenberger observer using bond graph

The model BG-LFT of the system with the observer of Luenberger is represented by the diagram in Figure 11.

Figure 11.

Representation in BG-LFT form of the system with observer from Luenberger.

3.5.5. Location of defects per observer bench

A single residue allows the detection of a fault. However, the location of a defect requires a set of structured residues. These residues must be designed to be sensitive to certain defects and insensitive to others. The use of observer banks constructed from only part of the inputs and/or outputs of the system makes it possible to respond to this problem. The symptoms generated are compared with fault signatures.

There are two observer banks:

  • The Dedicated Observer Scheme (DOS): the ith observer is controlled by the ith output (input) and all inputs (outputs). The other outputs (inputs) are considered unknown.

  • The Generalized Observer Scheme (GOS) structure: the ith observer is controlled by all outputs (inputs) except the ith and all inputs (outputs).

For a given observer bank, the signature of the various defects is defined in the signature table. The rows and columns of this table correspond, respectively, to defects and symptoms. The cells in the table are filled with binary values. A zero (0) means that the symptom is not sensitive to the defect. One (1) means that the symptom is sensitive to this defect.

Figures 12 and 13, respectively, present the principle of detection and localization of faults by observers (DOS) and (GOS).

Figure 12.

Sensor faults with DOS structure.

Figure 13.

Sensor faults with GOS structure.

Tables 3 and 4 illustrate the signatures of the various defects as a function of the residuals.


Table 3.

Signature of sensor faults with DOS structure.


Table 4.

Signature of sensor faults with GOS structure.

3.6. Example

Consider the same example (DC motor) given in Figure 7.

3.6.1. Proportional observer diagnosis using the bond graph model

The application of leap graph for the diagnosis of industrial systems is mainly justified by the fact that the model can be fine-tuned by adding or removing jump graph elements (graphically) according to simplifying assumptions. It is therefore also desirable that for the state estimation, the observer design is carried out using graphic methods and taking advantage of the properties.

3.6.2. Requirements

For observer construction, we need to check the following conditions:

The first condition is satisfied: When the bond graph model of the DC motor is put into a preferred integral causality, it is clear that there exists a causal path linking the sensors Df1and Df2to each dynamic element I in the integral causality (Figure 14(a)). The second condition is satisfied: when we put the jump engine model of the DC motor in derived causality, all the elements I admit a derivative causality, and the sensors Df1and Df2are dualized (Figure 14(b)).

Figure 14.

(a): Bond model of the DC motor in integral causality, (b): bond model of the DC motor in derived causality.

3.6.3. Calculation of residues

Figure 15 represents a construction of the proportional observer of the system using the bond graph approach.

Figure 15.

Proportional observer using bond graph of DC motor monitor.

We simulated the system with the 20sim software. Figure 16 shows the evolution of the real and estimated state variables.

Figure 16.

Output variables, (a): output variablesy1andŷ1, (b): output variablesy2andŷ2.

From model BG of Figure 15, the residues r1(t) and r2(t) can be deduced.

Residue r1: f2f2̂=Df2f2̂=0.



Residue r2: f6f6̂=Df1f6̂=0

We proceed in the same way, we end:


Figure 17 shows the convergence of residuals toward zero.

Figure 17.

(a): Residuer1(t), (b): Residuer2(t)in the case of normal operation.

3.6.4. Génération des résidus avec défauts capteurs

The sensors Df1and Df2are affected by defects (FC1and FC2), then:

Residue r1: Df1+FC1f6̂=0


Residue r2: Df2+FS2f2̂=0


Eqs. (30) and (31) show that residues r1(t)and r2(t)are sensitive to sensor defects. Figure 18 confirms the sensitivity of the residuals to the defects of the sensors Df1and Df2.

Figure 18.

(a): Residuer1(t)with defects of the sensorsDf1andDf2,(b): Residuer2(t)with defects of the sensorsDf1andDf2.

3.6.5. Robust diagnosis by observer from Luenberger

Figure 19 shows the Luenberger observer of the DC motor using the BG-LFT model.

Figure 19.

Luenberger observer of the DC motor using the BG-LFT model.

From the BG-LFT model (Figure 19), residues R1(t)and R2(t)can be deduced.

Residue R1: f2f2̂=0


Eq. (32) is composed of two parts: the first part corresponds to the evolution of the normal residue (r1n) and the second part represents the evolution related to the uncertainty of the parameters (d1).


Residue R2: f6f6̂=0


Eq. (33) is composed of two parts: the first part corresponds to the evolution of the normal residue (r2n) and the second part represents the evolution related to the uncertainty of the parameters (d2).


3.6.6. Detection and localization of sensor faults

Figure 20 shows the evolution of residues based on the BG model using the BG-DOS and BG-GOS structures, the residue r1(t)is sensitive to the defects occurring on the Df1sensor while the residue r2(t)is sensitive to faults on the Df2sensor.

Figure 20.

(a): Residuer1(t)andr2(t)for the DOS structure, (b): residuer1(t)andr2(t)for the DOS structure.

Table 5(a) and (b) represents the binary signatures of the defects for the deduced DOS and GOS structures to perfectly isolate the defects.

Structure DOSStructure GOS

Table 5.

Signatures of faults. (a) Structure DOS and (b) structure GOS.


4. Conclusion

In this chapter, we proposed a technique for the diagnosis of linear systems by Luenberger observer using the leap graph. We presented the leap graph approach for constructing a full-order observer and proposed a new BG-based observer method. Subsequently, we presented the systems with uncertain parameters modeled by the leap graph approach and we also proposed a new method of diagnosis of systems with uncertain parameters by Luenberger observer.

In the last part of this chapter, we developed and proposed an observer-based diagnostic technique (BG-DOS / BG-GOS) to detect and locate faults.


  1. 1. Khedher A, Othman KB, Benrejeb M, Maquin D. Adaptive observer for fault estimation in nonlinear systems described by a Takagi-Sugeno model. In: 18th Mediterrranean Conference on Control and Automation, Marrakech, Morocco; 2010. pp. 23-25
  2. 2. Maquin D, Ragot J. Diagnostic Des Systèmes Linéaires. Hermes ed. Paris, France: Hermes; 2000
  3. 3. Sueur C, Dauphin-Tanguy G. Structural Controllability and Observability of linear Systems Represented by Bond Graphs. Journal of the Franklin Institute. 1989;326:869-883
  4. 4. Pichardo-Almarza C, Rahmani A, Dauphin-Tanguy G, Delgado M. Observateur Proportionnel-Intégral pour des Systèmes Linéaires Modélisés par Bond Graphs. In: Actes des Journées Doctorales d’Automatique; 2003
  5. 5. Luenberger DG. An introduction to observers. IEEE Transactions on Automatic Control. 1971;16:596-602
  6. 6. Rahmani A, Sueur C, Dauphin-Tanguy G. Pole assignment for systems modelled by bond graph. Journal of the Franklin Institute. 1994;331(3):299-312
  7. 7. Djeziri MA, Ould Bouamama B, Merzouki R. Modelling and robust FDI of steam generator using uncertain bond graph model. Journal of Process Control. 2009;19:149-162
  8. 8. Dauphin-Tanguy G, Sié Kam C. How to Model Parameter Uncertainies in a Bond Graph Framework. Erlangen, Germany: ESS’99; 1999. pp. 121-125
  9. 9. Sié Kam C, Dauphin-Tanguy G. Bond graph models of structured parameter uncertainties. Journal of the Franklin Institute. 2005;342:379-399
  10. 10. Redheffer R. On a certain linear fractional transformation. Journal of Mathematics and Physics Banner. 1960;39:269-286. l'aide d'éléments caractérisés». 8ème colloque Electronique dePuissance du Futur EPF 2000, Lille; 2000
  11. 11. Alazard D, Cumer C, Apkarian P, Gauvrit M, Fereres G. Robustesse et Commande Optimale. Cépadues ed. 1999. ISBN: 2.85428.516.6
  12. 12. Sié Kam C. Les bond graphs pour la modélisation des systèmes linéaires incertains [Thèse de doctorat]. USTLille1-ECLille; Décembre 2001. N°d’ordre 3065, 2001
  13. 13. Djeziri MA. Diagnostic des systèmes incertains par l’approche bond graph [thèse de doctorat]. École Centrale de Lille; 2007

Written By

Abderrahmène Sallami

Submitted: January 26th, 2018 Reviewed: May 22nd, 2018 Published: December 12th, 2018