HYDIAG is a software developed in Matlab by the DISCO team at LAAS-CNRS. It is currently a software designed to simulate, diagnose and prognose hybrid systems using model-based techniques. An extension to active diagnosis is also provided. This paper aims at presenting the native HYDIAG tool, and its different extensions to prognosis and active diagnosis. Some results on an academic example are given.
HYDIAG is a software developed in Matlab, with Simulink. The development of this software was initiated in the DISCO team with contributions about diagnosis on hybrid systems [1]. It has undergone many changes and is currently a software designed to simulate, diagnose and prognose hybrid systems using model-based techniques [2; 3; 4]. An extension to active diagnosis has been also realized [5; 6]. This article aims at presenting the native HyDiag tool and its different extensions to prognosis and active diagnosis.
Section 2 recalls the hybrid formalism used by HYDIAG. Section 3 presents the native HYDIAG tool that simulates and diagnoses hybrid systems. Section 4 explains how HYDIAG has been extended in HYDIAGPRO to prognose and diagnose hybrid systems. Section 5 presents the extension to active diagnosis. Experimental results of HYDIAG and its extension HYDIAGPRO are finally presented in Section 6.
HYDIAG deals with hybrid systems defined in a monolithic way. Such a system must be modeled by a hybrid automaton [7]. Formally, a hybrid automaton is defined as a tupleS = (ζ, Q, Σ, T, C, (q0, ζ0)) where: • ζ is a finite set of continuous variables that comprises input variables u(t) ∈ Rnu , state variables x(t) ∈ Rnx , and output variables y(t) ∈ Rny . • Q is a finite set of discrete system states. • Σ is a finite set of events. • T ⊆ Q × Σ → Q is the partial transition function between states. • C = Sq∈Q Cq is the set of system constraints linking continuous variables.
Each state q ∈ Q represents a behavioural mode that is characterized by a set of constraints Cq that model the linear continuous dynamics (defined by their representations in the state space as a set of differential and algebraic equations). A behavioural mode can be nominal or faulty (anticipated faults). The unknown mode can be added to model all the non anticipated faulty situations. The discrete part of the hybrid automaton is given by M = (Q, Σ, T, q0), which is called the underlying discrete event system (DES). Σ is the set of events that correspond to discrete control inputs, autonomous mode changes and fault occurrences. The occurrence of an anticipated fault is modelled by a discrete event fi ∈ Σf ⊆ Σuo, where Σuo ⊆ Σ is the set of unobservable events. Σo ⊆ Σ is the set of observable events. Transitions of T model the instantaneous changes of behavioural modes. The continuous behaviour of the hybrid system is modelled by the so called underlying multimode system Ξ = (ζ, Q, C, ζ0). The set of directly measured variables is denoted by ζOBS ⊆ ζ.
An example of a hybrid system modeled by a hybrid automaton is shown in Figure 1. Each mode qi is characterized by state matrices Ai, Bi, Ci and Di.
Hybrid system u
q1 q2 C1 xY11((nn+)=1C)=1xA11(xn1)(+nD)+1uB(1nu)(n) σ12 C2 xY22((nn+)=1C)=2xA22(xn2)(+nD)+2uB(2nu)(n) σ21 y σ1i
qi Ci xYii((nn+)=1C)=ixAi(inx)i(+nD)+iuB(un()n) σ
…
The method developed in [1] for diagnosing faults on-line in hybrid systems can be seen as interlinking a standard diagnosis method for continuous systems, namely the parity space method, and a standard diagnosis method for DES, namely the diagnoser method [8].
HYDIAG allows the user to edit the modes of a hybrid automaton S as illustrated in Figure 1. To model the system, the user must first provide in the Graphical User Interface of the HYDIAG software the following information: the number of modes, the number of discrete events that can be observable or unobservable, and the sampling period used for the underlying multimode system (defined by the set of state matrices of the state space representation of each mode).
There are optional parameters that are helpful to initialize the mode matrices automatically before editing them: the number of entries for the continuous dynamics, the number of outputs for continuous dynamics, the dimensions of each matrix A. The number of entries (resp. outputs) must be the same for all the modes.
The simulator of the edited model has no restrictions on the number of modes or the order of the continuous dynamics, it is generically designed. Online computations are performed using Matlab / Simulink. Results provided by Matlab can be reused if a special need arises. Figure 2 shows an overview of the software interface.
HYDIAG automatically computes the analytical redundancy
relations (ARRs) by using the parity space approach [9].
Details of this computation can be found in [
The idea of HYDIAG is to capture both the continuous dynamics and the discrete dynamics within the same mathematical object. To do so, the discrete part of the hybrid system M = (Q, Σ, T, q0) is enriched with specific observable events that are generated from continuous information. The resulting automaton is called the Behaviour Automaton (BA) of the hybrid system. HYDIAG then builds the diagnoser of the Behaviour Automaton (see [8]) by using the DIADES1 software also developed within the DISCO team at LAAS-CNRS (see an example of diagnoser in Figure 7).
Given the built hybrid diagnoser, HYDIAG then loads a set of timed observations produced by the system and it provides at each observation time an update of the diagnosis 1http://homepages.laas.fr/ypencole/DiaDes/ of the system by triggering the current transition of the hybrid diagnoser that matches the current observation. It is possible to define in HYDIAG a simulation scenario for the modeled system with a duration and a time sample defined by the user. 3.2
The general architecture of HYDIAG and its two extensions (see the next sections for their description) is presented on Figure 3. Ellipses represent the objects handled by the software, rectangles with rounded edges depict HYDIAG functions and rectangles with straight edges correspond to external DIADES packages. The behaviour automaton is at the heart of the architecture as HYDIAG and both its extensions rely on it to perform diagnosis, active diagnosis and prognosis.
ActDiades Model display Enriched hybrid model
Specialized Active diagnosers
Additional Signature event ARRs computation
AND/OR Graph
AO* Algorithm
Behaviour
Automaton display Behaviour Automaton
Diades Prognoser prognosis prognosis
ActHyDiag Conditional
plan Conditional plan display
HyDiag Diagnoser display Diagnoser diagnosis
Diagnosis display diagnosis
Prognosis display
HyDiagPro
HYDIAG has been extended in order to provide a prognosis functionality to the software [4]. The prognosis function computes (1) the fault probability of the system in each behavioural mode, (2) the future fault sequence that will lead to the system failure, (3) the Remaining Useful Life (RUL) of the system.
In HYDIAGPRO, the initial hybrid model is enriched
by adding for each behavioural mode a set of aging laws:
S+ = (ζ, Q, Σ, T, C, F , (q0, ζ0)) where F = {F q, q ∈ Q}
and F q is a set of aging laws one for each anticipated fault
f ∈ Σf in mode q. The aging modeling framework that
is adopted in HYDIAGPRO is based on the Weibull
probabilistic model [
The prognoser uses the aging laws in S+ to predict fault occurrences (see Figure 3). The prognoser uses the current diagnosis result to update on-line these aging laws (the parameters γjq) according to the operation time in each behavioural mode. For each new result of diagnosis, the prognosis function computes the most likely sequence of dated faults that leads to the system failure. From this sequence is estimated the system RUL [4]. 5
The second extension of HYDIAG provides an active diagnosis functionality to the software (see Figure 3). The inputs are the same as for HYDIAG but an additional file indicates the events of S that are actions, as well as their respective cost. Based on the behaviour automaton, we compute a set of specialised active diagnosers (one per fault): such a diagnoser is able to predict, based on the behaviour automaton, whether a fault can be diagnosed with certainty by applying an action plan from a given ambiguous situation [6]. From these diagnosers, we also extract a planning domain as a AND/OR graph.
At runtime, when HYDIAG is diagnosing, the diagnosis might be ambiguous. An active diagnosis session can be launched as soon as a specialised active diagnoser can analyse that the current faulty situation is discriminable by applying some actions. If the active diagnosis session is launched, an AO∗ algorithm starts and computes a conditional plan from the AND-OR graph that optimises an action cost criterion. It is important to note that in the case of a system with continuous dynamics, only discrete actions are contained in the active diagnosis plan issued by ACTHYDIAG. In particular, it is assumed that if it is necessary to guide the system towards a value on continuous variables, the synthesis of control laws must be performed elsewhere. 6
h
hmax h2
h1
HYDIAGPRO has been tested on a water tank system (Figure 4) composed of one tank with two hydraulic pumps (P1, P2). Water flows through a valve at the bottom of the tank depending on the system control. Three sensors (h1, h2, hmax) detect the water level and allow to set the control of the pumps (on/off). It is assumed that the pumps may fail only if they are on. The discrete model of water tank and the controls of pumps are given in Figure 5. Discrete events in Σ = {h1, h2s, h2i, hmax, f1, f2} allow the system to switch into different modes. Observable events are Σo = {h1, h2s, h2i, hmax}. Two faults that correspond to the pump failures are anticipated Σf = {f1, f2} and are not observable.The Weibull parameter values of aging models F = {F qi } are reported in Table 1.
The underlying continuous behaviour of every discrete mode qi for i ∈ {1..8} is represented by the same state pump mode
Pump1 Pump 2
Time (h) Time (h) Time (h)
Left hand side of Figure 6 shows the diagnoser belief state just before and after the fault f1 occurrence. Results are consistent with the scenario: before 3500h, the belief states of the diagnoser are always tagged with a nominal diagnosis. After 3500h, all the states are tagged with f1.
Middle of Figure 6 illustrates the predicted date of fault occurrence (df1 and df2 ). At the beginning of the process, the prognosis result is: Π0 = ({f1, 4120}, {f2, 5105}). It can be noted that the predicted dates df1 and df2 of f1 and f2 globally increase. Indeed, the system oscillates between stressful modes and less stressful modes. To make it simple, we can consider that in some modes, the system does not degrade, so the predicted dates of f1 and f2 are postponed. Before 3500h, the predicted date of f1 is lower than the one of f2. After 3500h, the predicted date of f2 is updated, knowing that the system is in a degraded mode. Finally, the prognosis result is Π3501 = ({f2, 5541}). Figure 6 shows the evolution of the RUL of the system. At t = 3501, as the fault f2 is estimated to occur at t = 5541, the system RUL at t = 3501 is 5541 − 3501 = 2040h. 7
HYDIAG is a software developed in Matlab, with Simulink, by the DISCO team, at LAAS-CNRS. This tool has been extended into HYDIAGPRO to simulate, diagnose and prognose hybrid systems using model-based techniques. Some results on an academic example are exposed in the paper. An extension to active diagnosis is also presented. The active diagnosis algorithm is currently tested on a concrete industrial case. HYDIAG and its user manual will be soon available on the LAAS website.