EN FR
EN FR


Section: Research Program

Subspace-based identification and detection

For reasons closely related to the vibrations monitoring applications, we have been investigating subspace-based methods, for both the identification and the monitoring of the eigenstructure (λ,φλ) of the state transition matrix F of a linear dynamical state-space system :

X k + 1 = F X k + V k + 1 Y k = H X k , (10)

namely the (λ,ϕλ) defined by :

det ( F - λ I ) = 0 , ( F - λ I ) φ λ = 0 , ϕ λ = Δ H φ λ (11)

The (canonical) parameter vector in that case is :

θ = Δ Λ vec Φ (12)

where Λ is the vector whose elements are the eigenvalues λ, Φ is the matrix whose columns are the ϕλ's, and vec is the column stacking operator.

Subspace-based methods is the generic name for linear systems identification algorithms based on either time domain measurements or output covariance matrices, in which different subspaces of Gaussian random vectors play a key role  [32] . A contribution of ours, minor but extremely fruitful, has been to write the output-only covariance-driven subspace identification method under a form that involves a parameter estimating function, from which we define a residual adapted to vibration monitoring [1] . This is explained next.

Covariance-driven subspace identification.

Let Ri=Δ𝐄YkYk-iT and:

p + 1 , q = Δ R 0 R 1 R q - 1 R 1 R 2 R q R p R p + 1 R p + q - 1 = Δ Hank R i (13)

be the output covariance and Hankel matrices, respectively; and: G=Δ𝐄XkYkT. Direct computations of the Ri's from the equations (10 ) lead to the well known key factorizations :

R i = H F i G p + 1 , q = 𝒪 p + 1 ( H , F ) 𝒞 q ( F , G ) (14)

where:

𝒪 p + 1 ( H , F ) = Δ H H F H F p and 𝒞 q ( F , G ) = Δ ( G F G F q - 1 G ) (15)

are the observability and controllability matrices, respectively. The observation matrix H is then found in the first block-row of the observability matrix 𝒪. The state-transition matrix F is obtained from the shift invariance property of 𝒪. The eigenstructure (λ,φλ) then results from (11 ).

Since the actual model order is generally not known, this procedure is run with increasing model orders.

Model parameter characterization.

Choosing the eigenvectors of matrix F as a basis for the state space of model (10 ) yields the following representation of the observability matrix:

𝒪 p + 1 ( θ ) = Φ Φ Δ Φ Δ p (16)

where Δ=Δ diag (Λ), and Λ and Φ are as in (12 ). Whether a nominal parameter θ0 fits a given output covariance sequence (Rj)j is characterized by [1] :

𝒪 p + 1 ( θ 0 ) and p + 1 , q have the same left kernel space. (17)

This property can be checked as follows. From the nominal θ0, compute 𝒪p+1(θ0) using (16 ), and perform e.g. a singular value decomposition (SVD) of 𝒪p+1(θ0) for extracting a matrix U such that:

U T U = I s and U T 𝒪 p + 1 ( θ 0 ) = 0 (18)

Matrix U is not unique (two such matrices relate through a post-multiplication with an orthonormal matrix), but can be regarded as a function of θ0. Then the characterization writes:

U ( θ 0 ) T p + 1 , q = 0 (19)

Residual associated with subspace identification.

Assume now that a reference θ0 and a new sample Y1,,YN are available. For checking whether the data agree with θ0, the idea is to compute the empirical Hankel matrix ^p+1,q:

^ p + 1 , q = Δ Hank R ^ i , R ^ i = Δ 1 / ( N - i ) k = i + 1 N Y k Y k - i T (20)

and to define the residual vector:

ζ N ( θ 0 ) = Δ N vec U ( θ 0 ) T ^ p + 1 , q (21)

Let θ be the actual parameter value for the system which generated the new data sample, and 𝐄θ be the expectation when the actual system parameter is θ. From (19 ), we know that ζN(θ0) has zero mean when no change occurs in θ, and nonzero mean if a change occurs. Thus ζN(θ0) plays the role of a residual.

It is our experience that this residual has highly interesting properties, both for damage detection [1] and localization [3] , and for flutter monitoring [8] .

Other uses of the key factorizations.

Factorization ( 3.5.1 ) is the key for a characterization of the canonical parameter vector θ in (12 ), and for deriving the residual. Factorization (14 ) is also the key for :

  • Proving consistency and robustness results [6] ;

  • Designing an extension of covariance-driven subspace identification algorithm adapted to the presence and fusion of non-simultaneously recorded multiple sensors setups [7] ;

  • Proving the consistency and robustness of this extension [9] ;

  • Designing various forms of input-output covariance-driven subspace identification algorithms adapted to the presence of both known inputs and unknown excitations [10] .

Research program

The research will first focus on the extension and implementation of current techniques as developed in I4S and IFSTTAR. Before doing any temperature rejection on large scale structures as planned, we need to develop good and accurate models of thermal fields. We also need to develop robust and efficient versions of our algorithms, mainly the subspace algorithms before envisioning linking them with physical models. Briefly, we need to mature our statistical toolset as well as our physical modeling before mixing them together later on.

Direct vibration modeling under temperature changes

This task builds upon what has been achieved in the CONSTRUCTIF project, where a simple formulation of the temperature effect has been exhibited, based on relatively simple assumptions. The next step is to generalize this modeling to a realistic large structure under complex thermal changes. Practically, temperature and resulting structural prestress and pre strains of thermal origin are not uniform and civil structures are complex. This leads to a fully 3D temperature field, not just a single value. Inertia effects also forbid a trivial prediction of the temperature based on current sensor outputs while ignoring past data. On the other side, the temperature is seen as a nuisance. That implies that any damage detection procedure has first to correct the temperature effect prior to any detection.

Modeling vibrations of structures under thermal prestress does and will play an important role in the static correction of kinematic measurements, in health monitoring methods based on vibration analysis as well as in durability and in the active or semi-active control of civil structures that by nature are operated under changing environmental conditions. As a matter of fact, using temperature and dynamic models the project aims at correcting the current vibration state from induced temperature effects, such that damage detection algorithms rely on a comparison of this thermally corrected current vibration state with a reference state computed or measured at a reference temperature. This approach is expected to cure damage detection algorithms from the environmental variations.

I4S will explore various ways of implementing this concept, notably within the FUI SIPRIS project.

Damage localization algorithms (in the case of localized damages such as cracks)

During the CONSTRUCTIF project, both feasibility and efficiency of some damage detection and localization algorithms were proved. Those methods are based on the tight coupling of statistical algorithms with finite element models. It has been shown that effective localization of some damaged elements was possible, and this was validated on a numerical simulated bridge deck model. Still, this approach has to be validated on real structures.

On the other side, new localization algorithms are currently investigated such as the one developed conjointly with University of Boston and tested within the framework of FP7 ISMS project. These algorithms will be implemented and tested on the PEGASE platform as well as all our toolset.

When possible, link with temperature rejection will be done along the lines of what has been achieved in the CONSTRUCTIF project.

Uncertainty quantification for system identification algorithms

Some emphasis will be put on expressing confidence intervals for system identification. It is a primary goal to take into account the uncertainty within the identification procedure, using either identification algorithms derivations or damage detection principles. Such algorithms are critical for both civil and aeronautical structures monitoring. It has been shown that confidence intervals for estimation parameters can theoretically be related to the damage detection techniques and should be computed as a function of the Fisher information matrix associated to the damage detection test. Based on those assumptions, it should be possible to obtain confidence intervals for a large class of estimates, from damping to finite elements models. Uncertainty considerations are also deeply investigated in collaboration with Dassault Aviation in Mellinger PhD thesis or with Northeastern University, Boston, within Gallegos PhD thesis.