# Stiffness matrix modification with vibration test data by displacement feedback technique

## Yongxin Yuan1

1School of Mathematics and Statistics, Hubei Normal University, Huangshi, 435002, P. R. China

Journal of Vibroengineering, Vol. 19, Issue 7, 2017, p. 5149-5160. https://doi.org/10.21595/jve.2017.18407
Received 30 March 2017; received in revised form 20 May 2017; accepted 27 May 2017; published 15 November 2017

Copyright © 2017 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Views 90
Abstract.

A no spill-over method is developed which uses measured normal modes and natural frequencies to adjust a structural dynamics model in light of displacement feedback technique. By the method, the required displacement feedback gain matrix is determined, and thus the updated stiffness matrix which satisfies the characteristic equation is found in the Frobenius norm sense and the large number of unmeasured high-order modal data of the original model is preserved. The method directly identifies, without iteration, and the solution of this problem is of a compact expression. The numerical example shows that the modal measured data are better incorporated into the updated model.

Keywords: model updating, undamped vibration system, displacement feedback, modal measured data, optimal approximation.

#### 1. Introduction

After spatial discretization by using the finite element method, the equation of motion of a linear elastic time-invariant structure with $n$ degrees of freedom is given by:

(1)
${M}_{a}\stackrel{¨}{\mathbf{q}}\left(t\right)+{K}_{a}\mathbf{q}\left(t\right)=\mathbf{f}\left(t\right).$

In which ${M}_{a}\in {\mathbb{R}}^{n×n}$, analytical mass matrix, is symmetric positive definite and ${K}_{a}\in {\mathbb{R}}^{n×n}$, analytical stiffness matrix, is symmetric positive semidefinite. $\mathbf{q}\left(t\right)\in {\mathbb{R}}^{n×1}$ is the displacement vector and $\mathbf{f}\left(t\right)\in {\mathbb{R}}^{n×1}$ is the external force vector. Eq. (1) is known as the finite element (FE) analytical model of the structure. Let $\mathbf{q}\left(t\right)=\mathbf{x}{e}^{i\omega t}$ be a soluton of the homogeneous part of Eq. (1), then we can get the following eigenvalue-eigenvector equation:

(2)

where ${\lambda }_{j}={\omega }_{j}^{2}$. Let:

where:

Then Eq. (2) can be equivalently written as:

(3)
${M}_{a}X\mathrm{\Lambda }={K}_{a}X.$

Assume that the modal orthogonality relationship is satisfied:

(4)
$\left[\begin{array}{l}{X}_{1}^{Τ}\\ {X}_{2}^{Τ}\end{array}\right]{M}_{a}\left[{X}_{1},{X}_{2}\right]={I}_{n}.$

By substituting Eq. (4) into Eq. (3), we can get another orthogonality relation:

(5)
$\left[\begin{array}{l}{X}_{1}^{Τ}\\ {X}_{2}^{Τ}\end{array}\right]{K}_{a}\left[{X}_{1},{X}_{2}\right]=\mathrm{\Lambda }.$

When active control forces are exerted on an undamped vibration system, Eq. (1) now becomes:

(6)
${M}_{a}\stackrel{¨}{\mathbf{q}}\left(t\right)+{K}_{a}\mathbf{q}\left(t\right)=B\mathbf{u}\left(t\right),$

where $B\in {\mathbb{R}}^{n×m}$ is the full column rank control feedback matrix and $\mathbf{u}\left(t\right)\in {\mathbb{R}}^{m×1}$ is the control vector. In discussing the feedback control, we assume that the control vector $\mathbf{u}\left(t\right)$ is defined by the control law:

(7)
$\mathbf{u}\left(t\right)=-F\mathbf{q}\left(t\right),$

where $F\in {\mathbb{R}}^{m×n}$ is displacement feedback gain matrix. By substituting Eq. (7) into Eq. (6) yields the following closed-loop system:

(8)
${M}_{a}\stackrel{¨}{\mathbf{q}}\left(t\right)+\left({K}_{a}+BF\right)\mathbf{q}\left(t\right)=0.$

Separation of variables by $\mathbf{q}\left(t\right)=\mathbf{y}{e}^{i\omega t},$ which leads to the generalized eigenvalue problems of Eq. (8) as follows:

(9)

where ${\sigma }_{j}={\omega }_{j}^{2}$.

In engineering practice, accurate mathematical models are required in order to predict their dynamic characteristics accurately. However, current FE analysis cannot provide sufficiently accurate FE models, which are in good agreement with measured results. A vibration engineer then faces the problem of updating the existing FE model with minimal changes so that the updated FE model can better reflect the measured data from the physical structure being modelled . The updated model may then be considered a better predictions of the responses of the structure, and can be used with greater confidence for damage detection, health monitoring and structural control, and so on.

Model updating techniques are now extensively developed and studied for the structural systems. For undamped systems, to update coefficient matrices using vibration test data by means of direct matrix-updating methods has been considered by Thoren , Baruch and Bar-Itzhack , Berman and Nagy , Wei , Modak et al. , Yang and Chen , Yuan , Yuen , Yuan and Liu , Modak  and Sarmadi et al. . For damped structural systems, the theory and methods have been discussed by Friswell et al. , Kuo et al. , Bai , Lancaster , Yuan and Dai , Yuan and Liu  and Mao and Dai . Carvalho et al.  presented a numerical method for the stiffness matrix updating problem in an undamped model. by the IMDH (incomplete measured data handling) method, they overcame the difficulty of the incomplete measured data in an algorithmic way without using standard modal expansion or reduction techniques. The method is also capable of preserving the large number of eigendata of the FE model that are not affected by updating. Very more recently, Sehgal and Kumar  presented a review of structural dynamic model updating techniques. A number of direct and iterative techniques of model updating along with their applications to real life systems are reviewed and a number of future research directions have been highlighted which can be used for further advancements in the field of model updating.

It is of practical importance for updating an existent model that the newly measured parameters enter the system without altering other unrelated high-order vibration parameters. Such an updating is known as no spill-over. Assume that ${\mathrm{\Lambda }}_{1}\in {\mathbb{R}}^{p×p}$, ${X}_{1}\in {\mathbb{R}}^{n×p}$ are known for the first $p$ eigenvalues and associated eigenvectors of the original system. The remaining $n-p$ unknown eigenpairs ${\mathrm{\Lambda }}_{2}\in {\mathbb{R}}^{\left(n-p\right)×\left(n-p\right)}$ and ${X}_{2}\in {\mathbb{R}}^{n×\left(n-p\right)}$ remain unchanged. Some explanations for updating with no spill-over being required see the Ref. .

The aim of this paper is to modify the stiffness matrix by using the displacement feedback so that the modified second-order system Eq. (8) will contain a number of measured eigenvalues and eigenvectors. Mathematically, the problem of updating stiffness matrix with no spill-over by displacement feedback, therefore, can be stated as follows.

Problem 1. Let $B\in {\mathbb{R}}^{n×m}$ be a full column rank matrix. Assume that and ${Y}_{1}=\left[{\mathbf{y}}_{1},\cdots ,{\mathbf{y}}_{p}\right]\in {\mathbb{R}}^{n×p}$ are respectively the measured eigenvalue and eigenvector matrices, where $p\ll n.$ Find displacement feedback gain matrix $F\in {\mathbb{R}}^{m×n}$ such that:

(10)
${M}_{a}{X}_{2}{\mathrm{\Lambda }}_{2}=K{X}_{2},$
(11)
${M}_{a}{Y}_{1}{\mathrm{\Sigma }}_{1}=K{Y}_{1},$

where $K={K}_{a}+BF$ and $K={K}^{Τ}$. Generally, the solution of Problem 1 is not unique. Thus, we need to solve the following least-squares approximation problem.

Problem 2. Find $\stackrel{^}{F}\in {\mathbf{S}}_{E}$ such that:

(12)
$‖B\stackrel{^}{F}‖=\underset{F\in {\mathbf{S}}_{E}}{\mathrm{m}\mathrm{i}\mathrm{n}}‖BF‖,$

where ${\mathbf{S}}_{E}$ is the solution set of Problem 1. Once the solution $\stackrel{^}{F}$ of Problem 2 is obtained, the updated stiffness matrix can be expressed as:

(13)
$\stackrel{^}{K}={K}_{a}+B\stackrel{^}{F}.$

Many structural components are generally subjected to dynamic loadings in their working life. Very often these components may have to perform in severe dynamic environment where in the maximum damage results from the resonant vibration . Therefore, in order to avoid the undesired phenomena, one way is to use feedback control so that the unfavorable eigendata are replaced by some suitable ones [23-26]. The idea of using the eigenstructure assignment technique to solve the model updating problem has been considered by [27, 28]. The method can produce an updated FE model on damping and stiffness matrices that matches the measured modal data. More recently, Ouyang and Zhang  addressed passive structural modifications of mass-spring systems for partial assignment of natural frequencies, two solution methods were proposed to construct the required mass-normalised stiffness matrix, which satisfies the partial assignment requirement of natural frequencies and maintains the configuration of the original structure after modifications. Sen and Bhattacharya  adopted a control theory-based eigenstructure assignment technique to update the FE model of a linear time-invariant system. The proposed method uses state feedback to produce the gain matrix which in turn updates the existing system matrices through simultaneous assignment of eigenvalue and eigenvector pairs in the FE model generated system matrices. Richiedei and Trevisani  introduced a novel hybrid method for vibration control in lightly damped systems through the concurrent synthesis of passive structural modifications and active state feedback control gains. The passive modifications alter the set of eigenvectors that can be achieved through state feedback control and gives additional degrees of freedom in the controller synthesis, which overcoming the limitations of eigenstructure assignment through active control used alone.

It should be mentioned that the studies by Zhang et al. [33-35], Chu and Datta , Nichols and Kautsky , Datta et al.  and Lin and Wang  lead to a feedback design problem for a second-order control system. That consideration eventually results in either a full or a partial eigenstructure assignment for the eigenvalue problem. Nonetheless, these results cannot meet the basic requirement that the updated matrices should be symmetric. Our main contribution is to provide a new numerical method to solve the FE model updating problem using displacement feedback technique and the updated model has the following properties:

• The measured eigenvalues and eigenvectors will be embedded in the updated model.

• The updated stiffness matrix is also symmetric and positive semidefinite.

• The eigenvalues and eigenvectors corresponding to the unmeasured ones remain unchanged.

• The difference between the updated model and the original model is minimal.

The method directly identifies, without iteration, and works directly on the second-order system model. More importantly, the approach allows the control matrix to be specified beforehand and also leads naturally to a small norm solution of the feedback gain matrices. We believe that the method proposed should give considerable insight into the important model updating problem.

In what follows, in Section 2, by using the QR-decomposition and the singular value decomposition (SVD) of matrices, we provide a necessary and sufficient condition for the set ${\mathbf{S}}_{E}$ to be nonempty and construct the set ${\mathbf{S}}_{E}$ explicitly when it’s nonempty. In Section 3, when the set ${\mathbf{S}}_{E}$ is nonempty, we show that the solution of Problem 2 is unique and present the explicit expression of the unique solution $\stackrel{^}{F}$ of this problem. In Section 4, a numerical algorithm is proposed to determine the displacement feedback gain matrix and a numerical example is provided to demonstrate the effectiveness of the proposed method.

As usual, let ${\mathbb{R}}^{m×n}$ be the set of all $m×n$ real matrices and $\mathbb{S}{\mathbb{R}}^{n×n}$ the set of all $n×n$ symmetric matrices in ${\mathbb{R}}^{n×n}$. ${A}^{Τ}\text{,}$${A}^{+}$ and $‖A‖$ denote the transpose, the Moore-Penrose generalized inverse and the Frobenius norm of the matrix $A$, respectively. ${I}_{n}$ denotes the identity matrix of order $n$.

#### 2. The solution of Problem 1

In order to solve Problem 1, the following two lemmas are needed.

Lemma 1. If $\stackrel{~}{Y}\in {\mathbb{R}}^{m×l}$, $\stackrel{~}{Z}\in {\mathbb{R}}^{m×q}$, then $\stackrel{~}{Y}\stackrel{~}{N}=\stackrel{~}{Z}$ has a solution $\stackrel{~}{N}\in {\mathbb{R}}^{l×q}$ if and only if $\stackrel{~}{Y}{\stackrel{~}{Y}}^{+}\stackrel{~}{Z}=\stackrel{~}{Z}.$ In which case, the general solution of $\stackrel{~}{Y}\stackrel{~}{N}=\stackrel{~}{Z}$ can be expressed as $\stackrel{~}{N}={\stackrel{~}{Y}}^{+}\stackrel{~}{Z}+\left({I}_{l}-{\stackrel{~}{Y}}^{+}\stackrel{~}{Y}\right)\stackrel{~}{L}$, where $\stackrel{~}{L}\in {\mathbb{R}}^{l×q}$ is an arbitrary matrix .

Lemma 2. Suppose that $\stackrel{~}{A}$, $\stackrel{~}{B}\in {\mathbb{R}}^{n×p}$, then the matrix equation $\stackrel{~}{A}\stackrel{~}{Y}=\stackrel{~}{B}$ has a symmetric solution $\stackrel{~}{Y}\in \mathbb{S}{\mathbb{R}}^{p×p}$ if and only if :

in which case, the general symmetric solution is:

$\stackrel{~}{Y}={\stackrel{~}{A}}^{+}\stackrel{~}{B}+\left({I}_{p}-{\stackrel{~}{A}}^{+}\stackrel{~}{A}\right)\left({\stackrel{~}{A}}^{+}\stackrel{~}{B}{\right)}^{Τ}+\left({I}_{p}-{\stackrel{~}{A}}^{+}\stackrel{~}{A}\right)\stackrel{~}{J}\left({I}_{p}-{\stackrel{~}{A}}^{+}\stackrel{~}{A}\right),$

where $\stackrel{~}{J}\in \mathbb{S}{\mathbb{R}}^{p×p}$ is an arbitrary symmetric matrix.

Let:

(14)

Assume that the QR-decomposition of $B$ is:

(15)
$B=Q\left[\begin{array}{l}R\\ 0\end{array}\right],$

where is an $n×n$ orthogonal matrix (${Q}_{1}\in {\mathbb{R}}^{n×m}$) and $R$ is an $m×m$ nonsingular matrix. By Lemma 1 and Eq. (15), Eq. (14) with respect to $F$ is solvable if and only if:

(16)
${Q}_{2}^{Τ}\mathrm{\Delta }K=0,$

and thus the unique solution can be represented as:

(17)
$F={B}^{+}\mathrm{\Delta }K={R}^{-1}{Q}_{1}^{Τ}\mathrm{\Delta }K.$

By Lemma 2, Eq. (16) always has a symmetric solution $\mathrm{\Delta }K$ and the general symmetric solution is:

(18)
$\mathrm{\Delta }K={Q}_{1}S{Q}_{1}^{Τ},$

where $S\in \mathbb{S}{\mathbb{R}}^{m×m}$ is an arbitrary symmetric matrix. Substituting Eq. (18) into Eq. (17), we obtain:

(19)
$F={R}^{-1}S{Q}_{1}^{Τ}.$

By Eqs. (4) and (5), we get:

(20)

Using Eq. (20) and noting that ${X}_{2}$ is of full column rank, Eq. (10) is equivalent to:

(21)
${M}_{a}T=KP,$

where $T={M}_{a}^{-1}{K}_{a}{M}_{a}^{-1}-{X}_{1}{\mathrm{\Lambda }}_{1}{X}_{1}^{Τ}$, $P={M}_{a}^{-1}-{X}_{1}{X}_{1}^{Τ}.$ In practice, the matrices ${X}_{2}$ and ${\mathrm{\Lambda }}_{2}$ in Eq. (10) are usually unknown, we notice that the matrices ${X}_{2}$ and ${\mathrm{\Lambda }}_{2}$ don’t appear in Eq. (21) explicitly. It follows from Eq. (21) that:

${M}_{a}T=KP⇔\mathrm{\Delta }KP={M}_{a}T-{K}_{a}P.$

Observe that:

${M}_{a}T-{K}_{a}P={M}_{a}\left({M}_{a}^{-1}{K}_{a}{M}_{a}^{-1}-{X}_{1}{\mathrm{\Lambda }}_{1}{X}_{1}^{Τ}\right)-{K}_{a}\left({M}_{a}^{-1}-{X}_{1}{X}_{1}^{Τ}\right)$

That is, Eq. (10) is equivalent to:

(22)
$\mathrm{\Delta }KP=0.$

Substituting Eq. (18) into Eq. (22), we obtain:

${Q}_{1}S{Q}_{1}^{Τ}P=0,$

which implies that:

(23)
$S{Q}_{1}^{Τ}P=0.$

Assume that the singular value decomposition (SVD) of ${Q}_{1}^{Τ}P$ is:

(24)
${Q}_{1}^{Τ}P=U\left[\begin{array}{ll}\mathrm{\Theta }& 0\\ 0& 0\end{array}\right]{V}^{Τ},$

where , , $U=\left[{U}_{1},{U}_{2}\right]\in {\mathbb{R}}^{m×m}$ and $V=\left[{V}_{1},{V}_{2}\right]\in {\mathbb{R}}^{n×n}$ are orthogonal matrices, ${U}_{1}\in {\mathbb{R}}^{m×l}$, ${V}_{1}\in {\mathbb{R}}^{n×l}.$

Using Lemma 2 again, the general symmetric solution of Eq. (23) is:

(25)
$S={U}_{2}H{U}_{2}^{Τ},$

where $H\in \mathbb{S}{\mathbb{R}}^{\left(m-l\right)×\left(m-l\right)}$ is an arbitrary symmetric matrix. By substituting Eq. (25) into Eq. (18), we obtain the general symmetric solution of Eqs. (10) and (14) as:

(26)
$\mathrm{\Delta }K={Q}_{1}{U}_{2}H{U}_{2}^{Τ}{Q}_{1}^{Τ},$

where $H\in \mathbb{S}{\mathbb{R}}^{\left(m-l\right)×\left(m-l\right)}$ is an arbitrary symmetric matrix.

Now, to solve Problem 1 is equivalent to finding symmetric matrix $H$ such that:

(27)
${W}_{1}H{W}_{1}^{Τ}{Y}_{1}={M}_{a}{Y}_{1}{\mathrm{\Sigma }}_{1}-{K}_{a}{Y}_{1},$

where ${W}_{1}={Q}_{1}{U}_{2}\text{.}$ Observe that ${W}_{1}^{Τ}{W}_{1}={I}_{m-l}\text{,}$ if let ${W}_{2}=\left[{Q}_{1}{U}_{1},{Q}_{2}\right]\text{,}$ then ${W}_{2}^{Τ}{W}_{2}={I}_{n-m+l}$ and $W=\left[{W}_{1},{W}_{2}\right]\in {\mathbb{R}}^{n×n}$ be an othogonal matrix. Thus, the equation of Eq. (27) can be equivalently written as:

(28)
${W}_{2}^{Τ}\left({M}_{a}{Y}_{1}{\Sigma }_{1}-{K}_{a}{Y}_{1}\right)=0,$
(29)
$H{W}_{1}^{Τ}{Y}_{1}={W}_{1}^{Τ}\left({M}_{a}{Y}_{1}{\Sigma }_{1}-{K}_{a}{Y}_{1}\right).$

Let the SVD of ${W}_{1}^{Τ}{Y}_{1}$ be:

(30)
${W}_{1}^{Τ}{Y}_{1}=L\left[\begin{array}{ll}\mathrm{\Psi }& 0\\ 0& 0\end{array}\right]{J}^{Τ},$

where:

are orthogonal matrices, ${L}_{1}\in {\mathbb{R}}^{\left(m-l\right)×t}$, ${J}_{1}\in {\mathbb{R}}^{p×t}.$ By Lemma 2, Eq. (29) has a symmetric solution $H$ if and only if:

(31)

In which case, the general symmetric solution is:

(32)
$H={W}_{1}^{Τ}A\left({W}_{1}^{Τ}{Y}_{1}{\right)}^{+}+{\left({W}_{1}^{Τ}A\left({W}_{1}^{Τ}{Y}_{1}{\right)}^{+}\right)}^{Τ}\left({I}_{m-l}-{W}_{1}^{Τ}{Y}_{1}\left({W}_{1}^{Τ}{Y}_{1}{\right)}^{+}\right)+{L}_{2}G{L}_{2}^{Τ},$

where $A={M}_{a}{Y}_{1}{\mathrm{\Sigma }}_{1}-{K}_{a}{Y}_{1}$ and $G\in \mathbb{S}{\mathbb{R}}^{\left(m-l-t\right)×\left(m-l-t\right)}$ is an arbitrary symmetric matrix.

As a summary, we can get the following result.

Theorem 1. Let the QR-decomposition of $B$ be given by Eq. (15) and the SVD of the matrix ${Q}_{1}^{Τ}P$ be given by Eq. (24). Assume that ${W}_{1}={Q}_{1}{U}_{2},{W}_{2}=\left[{Q}_{1}{U}_{1},{Q}_{2}\right]$ and the SVD of the matrix ${W}_{1}^{Τ}{Y}_{1}$ is given by Eq. (30). If the conditions (28) and (31) hold, then Problem 1 is solvable and the solution set ${\mathbf{S}}_{E}$ of Problem 1 can be expressed as:

${\mathbf{S}}_{E}=\left\{F\in {\mathbb{R}}^{m×n}\left|F={R}^{-1}{U}_{2}H{U}_{2}^{Τ}{Q}_{1}^{Τ}\right\right\},$

where $P={M}_{a}-{X}_{1}{X}_{1}^{Τ}\text{,}$ and $H$ is given by Eq. (32) and $G\in \mathbb{S}{\mathbb{R}}^{\left(m-l-t\right)×\left(m-l-t\right)}$ is an arbitrary symmetric matrix.

#### 3. The solution of Problem 2

It has been shown in Section 2 that if the conditions (28) and (31) are satisfied, the solution set ${\mathbf{S}}_{E}$ is nonempty. Clearly, ${\mathbf{S}}_{E}$ is a closed convex subset of ${\mathbb{R}}^{m×n}$. It follows from the best approximation theorem  that there exists a unique solution $\stackrel{^}{F}$ in ${\mathbf{S}}_{E}$ such that Eq. (12) holds. Now, we will seek the unique solution $\stackrel{^}{F}$ in ${\mathbf{S}}_{E}$. For $F\in {\mathbf{S}}_{E}$, we can get:

$‖BF‖=‖{Q}_{1}{U}_{2}H{U}_{2}^{Τ}{Q}_{1}^{Τ}‖=‖{W}_{1}H{W}_{1}^{Τ}‖=‖{W}^{Τ}{W}_{1}H{W}_{1}^{Τ}W‖$

where $W=\left[{W}_{1},{W}_{2}\right].$ Thus, by Eq. (32) we can obtain:

${‖BF‖}^{2}={‖{H}_{0}+{L}_{2}G{L}_{2}^{Τ}‖}^{2}={‖{L}^{Τ}\left({H}_{0}+{L}_{2}G{L}_{2}^{Τ}\right)L‖}^{2}$

where:

(33)
${H}_{0}={W}_{1}^{Τ}A\left({W}_{1}^{Τ}{Y}_{1}{\right)}^{+}+{\left({W}_{1}^{Τ}A\left({W}_{1}^{Τ}{Y}_{1}{\right)}^{+}\right)}^{Τ}\left(I-{W}_{1}^{Τ}{Y}_{1}\left({W}_{1}^{Τ}{Y}_{1}{\right)}^{+}\right).$

Therefore, $‖BF‖=\mathrm{m}\mathrm{i}\mathrm{n}$ if and only if:

(34)
$G=-{L}_{2}^{Τ}{H}_{0}{L}_{2}.$

By substituting Eq. (34) into Eqs. (19) and (32), we obtain the following result.

Theorem 2. If the conditions (28) and (31) hold, then Problem 2 has a unique solution and it can be described as:

(35)
$\stackrel{^}{F}={R}^{-1}{U}_{2}\stackrel{^}{H}{U}_{2}^{Τ}{Q}_{1}^{Τ},$

where:

(36)
$\stackrel{^}{H}={H}_{0}-{L}_{2}{L}_{2}^{Τ}{H}_{0}{L}_{2}{L}_{2}^{Τ},$

and ${H}_{0}$ is given by Eq. (33).

#### 4. A numerical example

Based on Theorems 1 and 2 we can establish an algorithm for solving Problems 1 and 2 as follows.

Algorithm.

1) Input ${M}_{a}$, ${K}_{a}$, $B$, ${X}_{1}$, ${\mathrm{\Lambda }}_{1}$, ${\mathrm{\Sigma }}_{1}$, ${Y}_{1}.$

2) Compute the QR-decomposition of $B$ by Eq. (15).

3) Compute $P={M}_{a}^{-1}-{X}_{1}{X}_{1}^{Τ}$ and the SVD of the matrix ${Q}_{1}^{Τ}P$ by Eq. (24).

4) Compute ${W}_{1}={Q}_{1}{U}_{2}$ and ${W}_{2}=\left[{Q}_{1}{U}_{1},{Q}_{2}\right],$ and the SVD of the matrix ${W}_{1}^{Τ}{Y}_{1}$ by Eq. (30).

5) Compute $A={M}_{a}{Y}_{1}{\mathrm{\Sigma }}_{1}-{K}_{a}{Y}_{1}.$

6) If the conditions (28) and (31) hold, then continue, otherwise, go to 1).

7) Compute ${H}_{0}$ and $\stackrel{^}{H}$ by Eqs. (33) and (36), respectively.

8) Compute $\stackrel{^}{F}$ by Eq. (35).

9) Compute $\stackrel{^}{K}$ by Eq. (13).

Example. Consider a cantilever beam model (see, Fig. 1). The cross section of the beam is rectangular with length 2 m, width 60 mm and height 3 mm, respectively. The material of the cantilever beam is aluminum alloy with the modulus of elasticity = 71 GPa, Poisson ratio = 0.33 and mass density = 2.714×10-5 N/mm3$.$ The beam is discretised into 10 elements shown in Fig. 1.

Fig. 1. The model of a cantilever beam The mass matrix of the FE model is diagonal and the stiffness matrix is 3-diagonal, which are given by:

${K}_{a}=\left[\begin{array}{cccccccccc}28755& -14378& 0& 0& 0& 0& 0& 0& 0& 0\\ -14378& 28755& -14378& 0& 0& 0& 0& 0& 0& 0\\ 0& -14378& 28755& -14378& 0& 0& 0& 0& 0& 0\\ 0& 0& -14378& 28755& -14378& 0& 0& 0& 0& 0\\ 0& 0& 0& -14378& 28755& -14378& 0& 0& 0& 0\\ 0& 0& 0& 0& -14378& 28755& -14378& 0& 0& 0\\ 0& 0& 0& 0& 0& -14378& 28755& -14378& 0& 0\\ 0& 0& 0& 0& 0& 0& -14378& 28755& -14378& 0\\ 0& 0& 0& 0& 0& 0& 0& -14378& 28755& -14378\\ 0& 0& 0& 0& 0& 0& 0& 0& -14378& 14378\end{array}\right].$

${\mathrm{\Lambda }}_{1}$ and ${X}_{1}$ are given by:

The measured modal data are given by:

Let control feedback matrix $B$ be:

$B=\left[\begin{array}{cccc}9.0206e-017& -0.013541& -0.12815& 0.36264\\ 0.35388& -0.072726& -0.24211& 0.53765\\ -0.45575& -0.19181& -0.32024& 0.44423\\ 0.14096& -0.34151& -0.33384& 0.15324\\ -0.12794& -0.46385& -0.25796& -0.15566\\ 0.30127& -0.50566& -0.08754& -0.29933\\ 0.11106& -0.45106& 0.15068& -0.1971\\ -0.62835& -0.33348& 0.39956& 0.084173\\ 0.36456& -0.21942& 0.58834& 0.37206\\ 8.3267e-017& -0.08628& 0.32903& 0.24579\end{array}\right].$

It is easy to check that the conditions (28) and (31) hold:

$‖{W}_{2}^{Τ}\left({M}_{a}{Y}_{1}{\mathrm{\Sigma }}_{1}-{K}_{a}{Y}_{1}\right)‖=3.7989e-011,$
$‖{W}_{1}^{Τ}A\left({W}_{1}^{Τ}{Y}_{1}{\right)}^{+}{W}_{1}^{Τ}{Y}_{1}-{W}_{1}^{Τ}A‖=4.5297e-013,$
$‖{Y}_{1}^{Τ}{W}_{1}{W}_{1}^{Τ}A-{A}^{Τ}{W}_{1}{W}_{1}^{Τ}{Y}_{1}‖=3.005e-011.$

By the Algorithm, we can obtain the following unique solution of Problem 2:

${\stackrel{^}{F}}^{Τ}=\left[\begin{array}{cccc}-1.4248e-014& 67.459& -98.657& 106.18\\ -2.2092e-014& 88.246& -155.53& 158.09\\ -2.0327e-014& 42.007& -148.4& 130.64\\ -1.0395e-014& -54.103& -83.778& 40.537\\ 2.669e-015& -155.7& 8.6985& -64.427\\ 1.3126e-014& -216.78& 93.102& -133.09\\ 1.7578e-014& -215.7& 145.36& -139.3\\ 1.6446e-014& -166.45& 162.77& -94.79\\ 1.3131e-014& -108.89& 160.22& -39.389\\ 5.7334e-015& -41.968& 78.221& -7.521\end{array}\right].$

In which case, the optimal updated stiffness matrix can be figured out.

$\stackrel{^}{K}={K}_{a}+B\stackrel{^}{F}$

Fig. 2 indicates the absolute values of the stiffness discrepancy matrix obtained from the proposed model updating method. We define the residual as:

and the numerical results are shown in the Tables 1 and 2.

Fig. 2. The absolute values of the stiffness discrepancy matrix (SDM) using the proposed method Fig. 3. The frequencies of the measured and updated models Fig. 4. The frequencies of the analytical and updated models Table 1 shows that the measured modal data are embedded in the new model ${M}_{a}{Y}_{1}{\mathrm{\Sigma }}_{1}=\stackrel{^}{K}{Y}_{1}$ (Fig. 3 shows the comparison of frequencies of the measured model with the updated model) and Table 2 implies that the model is updated with no spill-over (Fig. 4 shows the comparison of frequencies of the analysis model with the updated model) and the updated stiffness matrix is also symmetric and positive definite.

Table 1. Numerical results

 Eigenpairs $\left({\sigma }_{1},{\mathbf{y}}_{1}\right)$ 1.6404e-011 $\left({\sigma }_{2},{\mathbf{y}}_{2}\right)$ 3.0154e-011 $\left({\sigma }_{3},{\mathbf{y}}_{3}\right)$ 2.2863e-011

Table 2. Numerical results

 Eigenpairs $\left({\lambda }_{4},{\mathbf{x}}_{4}\right)$ 7.0668e-011 $\left({\lambda }_{5},{\mathbf{x}}_{5}\right)$ 6.7161e-011 $\left({\lambda }_{6},{\mathbf{x}}_{6}\right)$ 3.3355e-011 $\left({\lambda }_{7},{\mathbf{x}}_{7}\right)$ 2.1904e-011 $\left({\lambda }_{8},{\mathbf{x}}_{8}\right)$ 2.7824e-011 $\left({\lambda }_{9},{\mathbf{x}}_{9}\right)$ 5.457e-011 $\left({\lambda }_{10},{\mathbf{x}}_{10}\right)$ 6.671e-011

#### 5. Conclusions

A no spill-over direct updating method for undamped vibration systems with vibration test data using displacement feedback technique has been presented. When the conditions (28) and (31) hold, the required displacement feedback gain matrix can be determined, and the optimal updated stiffness matrix which satisfies the characteristic equation can be achieved. The method is easy to implement, and allows the control matrix to be specified beforehand. Although the proposed method can guarantee that the updated matrix is symmetric and positive semidefinite, it failed to preserve the pattern of the FE stiffness matrix. How to maintain the physical connectivity of the updated matrix is worthy of further study.

1. Friswell M. I., Mottershead J. E. Finite Element Model Updating in Structural Dynamics. Klumer Academic Publishers, Dordrecht, 1995. [Publisher]
2. Thoren A. R. Derivation of mass and stiffness matrices from dynamic test data. AIAA Conference Paper, No. 72-346, 1972, 6 pages. [Publisher]
3. Baruch M., Bar Itzhack I.-Y. Optimal weighted orthogonalization of measured modes. AIAA Journal, Vol. 16, 1978, p. 346-351. [Search CrossRef]
4. Berman A., Nagy E. J. Improvement of a large analytical model using test data. AIAA Journal, Vol. 21, 1983, p. 1168-1173. [Search CrossRef]
5. Wei F. S. Mass and stiffness interaction effects in analytical model modification. AIAA Journal, Vol. 28, 1990, p. 1686-1688. [Search CrossRef]
6. Modak S. V., Kundra T. K., Nakra B. C. Model updating using constrained optimization. Mechanics Research Communications, Vol. 5, 2002, p. 543-551. [Search CrossRef]
7. Yang Y. B., Chen Y. J. A new direct method for updating structural models based on measured modal data. Engineering Structures, Vol. 31, 2009, p. 32-42. [Search CrossRef]
8. Yuan Y. A symmetric inverse eigenvalue problem in structural dynamic model updating. Applied Mathematics and Computation, Vol. 213, 2009, p. 516-521. [Search CrossRef]
9. Yuen K.-V. Updating large models for mechanical systems using incomplete modal measurement. Mechanical Systems and Signal Processing, Vol. 28, 2012, p. 297-308. [Search CrossRef]
10. Yuan Y., Liu H. An iterative updating method for undamped structural systems. Meccanica, Vol. 47, 2012, p. 699-706. [Search CrossRef]
11. Modak S. V. Direct matrix updating of vibroacoustic finite element models using modal test data. AIAA Journal, Vol. 52, 2014, p. 1386-1392. [Search CrossRef]
12. Sarmadi H., Karamodin A., Entezami A. A new iterative model updating technique based on least squares minimal residual method using measured modal data. Applied Mathematical Modelling, Vol. 40, 2016, p. 10323-10341. [Search CrossRef]
13. Friswell M. I., Inman D. J., Pilkey D. F. The direct updating of damping and stiffness matrices. AIAA Journal, Vol. 36, 1998, p. 491-493. [Search CrossRef]
14. Kuo Y.-C., Lin W.-W., Xu S.-F. New methods for finite element model updating problems. AIAA Journal, Vol. 44, 2006, p. 1310-1316. [Search CrossRef]
15. Bai Z. J. Constructing the physical parameters of a damped vibrating system from eigendata. Linear Algebra and its Applications, Vol. 428, 2008, p. 625-656. [Search CrossRef]
16. Lancaster P. Inverse spectral problems for semisimple damped vibrating systems. SIAM Journal on Matrix Analysis and Applications, Vol. 29, 2007, p. 279-301. [Search CrossRef]
17. Yuan Y., Dai H. On a class of inverse quadratic eigenvalue problem. Journal of Computational and Applied Mathematics, Vol. 235, 2011, p. 2662-2669. [Search CrossRef]
18. Yuan Y., Liu H. A gradient based iterative algorithm for solving structural dynamics model updating problems. Meccanica, Vol. 48, 2013, p. 2245-2253. [Search CrossRef]
19. Mao X., Dai H. A quadratic inverse eigenvalue problem in damped structural model updating. Applied Mathematical Modelling, Vol. 40, 2016, p. 6412-6423. [Search CrossRef]
20. Carvalho J., Datta B. N., Gupta A., Lagadapati M. A direct method for model updating with incomplete measured data and without spurious modes. Mechanical Systems and Signal Processing, Vol. 21, 2007, p. 2715-2731. [Search CrossRef]
21. Sehgal S., Kumar H. Structural dynamic model updating techniques: a state of the art review. Archives of Computational Methods in Engineering, Vol. 23, 2016, p. 515-533. [Search CrossRef]
22. Korabathina R., Koppanati M. S. Linear free vibration analysis of rectangular Mindlin plates using coupled displacement field method. Mathematical Models in Engineering, Vol. 2, 2016, p. 41-47. [Search CrossRef]
23. Chu E. K.-W., Datta B. N. Numerically robust pole assignment for second-order systems. International Journal of Control, Vol. 64, 1996, p. 1113-1127. [Search CrossRef]
24. Nichols N. K., Kautsky J. Robust eigenstructure assignment in quadratic matrix polynomials: Nonsingular case. SIAM Journal on Matrix Analysis and Applications, Vol. 23, 2001, p. 77-102. [Search CrossRef]
25. Datta B. N., Elhay S., Ram Y. M., Sarkissian D. R. Partial eigenstructure assignment for the quadratic pencil. Journal of Sound and Vibration, Vol. 230, 2000, p. 101-110. [Search CrossRef]
26. Lin W.-W., Wang J.-N. Partial pole assignment for the quadratic pencil by output feedback control with feedback designs. Numerical Linear Algebra with Applications, Vol. 12, 2005, p. 967-979. [Search CrossRef]
27. Zimmerman D., Widengren M. Correcting finite element models using a symmetric eigenstructure assignment technique. AIAA Journal, Vol. 28, 1990, p. 1670-1676. [Search CrossRef]
28. Minas C., Inman D. J. Correcting finite element models with measured modal results using eigenstructure assignment methods. Proceedings of the 4th IMAC Conference, Schenectady, NY, 1987, p. 583-587. [Search CrossRef]
29. Ouyang H., Zhang J. Passive modifications for partial assignment of natural frequencies of mass-spring systems. Mechanical Systems and Signal Processing, Vol. 50, Issue 51, 2015, p. 214-226. [Search CrossRef]
30. Sen S., Bhattacharya B. Non-iterative eigenstructure assignment technique for finite element model updating. Journal of Civil Structural Health Monitoring, Vol. 5, 2015, p. 365-375. [Search CrossRef]
31. Richiedei D., Trevisani A. Simultaneous active and passive control for eigenstructure assignment in lightly damped systems. Mechanical Systems and Signal Processing, Vol. 85, 2017, p. 556-566. [Search CrossRef]
32. Chu M. T., Datta B. N., Lin W. W., Xu S. F. Spillover phenomenon in quadratic model updating. AIAA Journal, Vol. 46, 2008, p. 420-428. [Search CrossRef]
33. Zhang J., Ouyang H., Yang J. Partial eigenstructure assignment for undamped vibration systems using acceleration and displacement feedback. Journal of Sound and Vibration, Vol. 333, 2014, p. 1-12. [Search CrossRef]
34. Zhang J., Ye J., Ouyang H. Static output feedback for partial eigenstructure assignment of undamped vibration systems. Mechanical Systems and Signal Processing, Vol. 68-69, 2016, p. 555-561. [Search CrossRef]
35. Zhang J., Yuan Y., Liu H. An approach to partial quadratic eigenvalue assignment of damped vibration systems using static output feedback. International Journal of Structural Stability and Dynamics, https://doi.org/10.1142/S0219455418500128. [Search CrossRef]
36. Ben Israel A., Greville T. N. E. Generalized Inverses. Theory and Applications. Second Ed., Springer, New York, 2003. [Search CrossRef]
37. Don F. J. H. On the symmetric solutions of a linear matrix equation. Linear Algebra and its Applications, Vol. 93, 1987, p. 1-7. [Search CrossRef]
38. Aubin J. P. Applied Functional Analysis. John Wiley & Sons, New York, 1979. [Search CrossRef]

#### Cited By

 Applied Sciences Meilong Chen, Huajiang Ouyang, Wanyou Li, Donghua Wang, Siyuan Liu 2019