30.06.2013 Views

smart technologies for safety engineering

smart technologies for safety engineering

smart technologies for safety engineering

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

128 Smart Technologies <strong>for</strong> Safety Engineering<br />

problems exists in which such a model is not known or does not exist. In this case, in order to<br />

estimate on-line the load acting on a mechanical system, the unknown input technique of Hou<br />

and Müller [17] can be used. In their approach, the only required in<strong>for</strong>mation is the model<br />

of the system (mass and stiffness matrices) and measurements of a part of the state vector<br />

(displacements or velocities).<br />

Consider a dynamical system described by Equation (11). Additionally assume that m < p,<br />

which means that the number of sensors exceeds the number of <strong>for</strong>ces to be identified.<br />

The first step of the present method is to decompose the system of Equation (11) into two<br />

subsystems, so that the unknown input u(t) acts on one of them only. By the singular value<br />

decomposition (SVD), the matrix B can be restated in the following <strong>for</strong>m:<br />

B = UΣV T<br />

where U ∈ R n×n and V ∈ R m×m are unitary matrices, which satisfy the following condition:<br />

Matrix Σ is a rectangular matrix (n × m) of the <strong>for</strong>m<br />

Σ =<br />

(14)<br />

UU T = I, VV T = I (15)<br />

<br />

Σ0<br />

=<br />

0<br />

⎡<br />

σ1<br />

⎢<br />

0<br />

⎢ .<br />

⎢ .<br />

⎢ .<br />

⎢ 0<br />

⎢ 0<br />

⎢ .<br />

⎣ .<br />

0<br />

σ2<br />

.<br />

0<br />

0<br />

.<br />

···<br />

···<br />

. ..<br />

···<br />

···<br />

. ..<br />

⎤<br />

0<br />

0 ⎥<br />

.<br />

⎥<br />

. ⎥<br />

σm ⎥<br />

0 ⎥<br />

. ⎥<br />

. ⎦<br />

0 0 ··· 0<br />

where matrix Σ0 has singular values of B on the main diagonal. By substituting Equation (14)<br />

into (11),<br />

˙x(t) = Ax(t) + UΣV T u(t)<br />

Next, multiplying on the left by U T and using the orthogonality (Equation (15)),<br />

Linear trans<strong>for</strong>mations of the variables,<br />

yields the two following coupled subsystems:<br />

U T ˙x(t) = U T Ax(t) + ΣV T u(t) (16)<br />

x := U¯x, u := Vū (17)<br />

<br />

˙¯x1(t) = Ā11 ¯x1(t) + Ā12 ¯x2(t) + Σ0ū(t)<br />

˙¯x2(t) = Ā21 ¯x1(t) + Ā22 ¯x2(t)<br />

(18)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!