Damage detection and identification of parameter matrices using residual force vector
HeeChang Eun^{1} , SuYong Park^{2}
^{1, 2}Kangwon National University, Department of Architectural Engineering, Chuncheon, Korea
^{1}Corresponding author
Journal of Measurements in Engineering, Vol. 2, Issue 1, 2014, p. 17.
Received 20 January 2014; received in revised form 23 January 2014; accepted 26 January 2014; published 30 March 2014
JVE Conferences
Beginning with incomplete mode shape measurement data, this study presents analytical equations to predict the actual stiffness and mass matrices. The measured modal data, including the measurement, manufacturing and modeling errors, should be updated for subsequent analysis. In this study, the incomplete mode shape data are expanded to a full set of degreesoffreedom (DOFs) based on the generalized inverse method and the concept of residual force vector. The corrected parameter matrices are straightforwardly derived using the estimated mode shape data and the pseudo inverse method. The validity of the proposed method is evaluated based on the number of measured modes in an application, and its limitations are investigated.
Keywords: residual force vector, eigenvalue function, update, constraint, measurement, data expansion.
1. Introduction
Modal data provide extremely useful information that can assist in the design of virtually any structure. The visualization of mode shapes is included in the design process, and the development of a modal model is useful for simulation and design studies. Measured and analytical data differ because of measurement noise and model inadequacies. It is desirable to correct the parameter matrices to provide satisfactory measurement data for subsequent analysis.
The parameter matrices can be corrected using the measured frequency response function (FRF) or mode shape data. The FRFbased parameter matrix update methods can utilize available measurement data in a certain frequency range and estimate the damping matrix [14].
The number of transducers for measurement is less than the number of DOFs in a finite element model, requiring the expansion of measured mode shapes. It is expected that a gap exists in the eigenvalues before and after updating and it should be reduced. There have been numerous attempts to expand the measured data [59]. The expansion of mode shape data requires the modification of the eigenvalue function. The residual force method is one modification method. The residual force is obtained by substituting the experimentally measured modal frequencies and mode shape, along with the analytical mass and stiffness matrices from the analytical model structure, into the eigenvalue function. Yang and Liu [10] presented a mode shape expansion method based on the best achievable eigenvector concept to solve the incomplete measurement problem. Zhao and Zhang [11], Yun et al. [12] and Gafka and Zimmerman [13] introduced a mode shape expansion method to estimate motion at all DOFs. Friswell and Mottershead [14] provided a comprehensive overview that illustrates many of the different techniques and factors involved in updating a finite element model. Many researchers have presented analytical methods to correct both analytical and actual stiffness matrices to satisfy the dynamic equation and the orthogonality constraints [1527].
This work presents analytical methods to expand incomplete measurement data to a full set of DOFs using the generalized inverse method provided by Eun [28] and the residual force vector. Using the estimated modal data set, the residual force vector, fundamental linear algebra and the MoorePenrose inverse, the variation in parameter matrices is predicted. The validity of the proposed method is evaluated based on the number of measured modes in an application, and its limitations are investigated.
2. Expansion of incomplete mode shape data
Considering modeling and measurement errors, the parameters of the finite element model should be corrected for proper subsequent analysis. However, it is not practical to measure the responses at a full set of DOFs of the system. The incomplete measurement data should be expanded, and the estimated full data utilized in predicting the parameter matrices. The mode shape expansion to satisfy the measured mode shape data is derived from the modification of the eigenvalue function.
The dynamic behavior of a structure, which is assumed linear and approximately discretized for $n$ DOFs without damping and external excitations, can be described by the following equation of motion:
where ${\mathbf{M}}_{a}$ and ${\mathbf{K}}_{a}$ denote the $n\times n$ analytical mass and stiffness matrices, respectively.
Assuming that the analytical parameter matrices ${\mathbf{M}}_{a}$ and ${\mathbf{K}}_{a}$ are exact in modeling, they should be corrected to eliminate the inconsistency between the analytical and experimental results. The expansion of incomplete data is derived using the generalized inverse method.
By substituting $\mathbf{u}={\mathbf{\phi}}_{a}{e}^{j\mathrm{\Omega}t}$, where $j=\sqrt{1}$ and ${\mathbf{\phi}}_{a}$ is the analytical mode shape vector, into Eq. (1), the eigenvalue equation for the undamped system is derived as follows:
Assume that the first $r$ ($r<n$) modal frequencies and the corresponding mode shape vectors at $s$ ($s<n$) measurement positions are provided for data expansion and parameter identification. The actual mode shape vector ${\mathbf{\phi}}_{x}$ is composed of unmeasured and measured sets. Assuming that the first $r$ mode shape sets at $s$ different locations are measured, the constraints at each mode are expressed by the $s$ measured mode shape data as follows:
where ${\mathbf{A}}_{i}$ is a $s\times n$ Boolean matrix defining the measurement locations at the $i$th mode and ${\mathbf{b}}_{i}$ is the corresponding $s\times 1$ measurement data vector.
The updated mode shape ${\mathbf{\phi}}_{i,x}$ is expanded to a full data set to satisfy the constraint equation of Eq. (3). The coefficient matrix $({\mathbf{K}}_{a}{\mathrm{\Omega}}_{i}^{2}{\mathbf{M}}_{a})$ in Eq. (2) is a rankdeficiency matrix. By inserting the frequencies ${\mathrm{\Omega}}_{i,x}$ ($i=1,2,\dots ,r$), $r<n$ and corresponding mode shape vector ${\mathbf{\phi}}_{i,x}$ into Eqn. (2), it is modified as the eigenvalue function of the full rank coefficient matrix. The matrix is expressed by the summation of the analytical mode shape ${\mathbf{\phi}}_{i}$ and the variation in the mode shape ${\u2206\mathbf{\phi}}_{i}$ as follows:
where ${\mathbf{\phi}}_{i}$ represents the mode shape vector obtained from the analytical eigenvalue equation of Eq. (2), and $\mathrm{\Delta}{\mathbf{\phi}}_{i}$ denotes the variation in mode shapes to be updated with the actual measurements. Based on the generalized inverse method, the variation in the mode shape can be derived.
The variation in mode shape is derived using the measured mode shape data based on the generalized inverse method. The variation is expressed as:
where ${\mathbf{R}}_{i}={\mathbf{K}}_{a}{\mathrm{\Omega}}_{i,x}^{2}{\mathbf{M}}_{a}$.
Substituting Eq. (5) into Eq. (4), the $i$th estimated mode shapes for a full set of DOFs are simply obtained. The estimated mode shape data are utilized to identify the parameter matrices.
3. Prediction of parameter matrices
We derive the mathematical forms of the corrected stiffness and mass matrices based on the estimated mode shapes. Based on the expanded mode shapes and the measured frequencies, the eigenvalue equation for the undamped system is expressed as:
where ${\mathrm{\Omega}}_{i,x}$ and ${\mathbf{\phi}}_{i,x}$ denote the $i$th measured frequency and corresponding expanded mode shape vector. Expressing the actual parameter matrices $\mathbf{M}$ and $\mathbf{K}$ by $\mathbf{M}={\mathbf{M}}_{a}+\mathrm{\Delta}\mathbf{M}$ and $\mathbf{K}={\mathbf{K}}_{a}+\mathrm{\Delta}\mathbf{K}$ and inserting them into Eq. (6), results in the following:
where ${\mathbf{R}}_{i}=\left({\mathbf{K}}_{a}{\mathrm{\Omega}}_{i,x}^{2}{\mathbf{M}}_{a}\right){\mathbf{\phi}}_{i,x}$, and $\u2206\mathbf{M}$ and $\u2206\mathbf{K}$ are the variation in the mass and stiffness matrices, respectively. Expanding Eq. (7) to the first s modes results in the following:
Solving Eq. (8) with respect to the parameter variation matrices results in the following:
Consequently, the variations in parameter matrices can be simply and explicitly obtained under the assumption that the mode shapes are exact. Eq. (9) exhibits that the exactness of the parameter matrices depends on the precision of estimated mode shape data. In the following application, the validity of the proposed method is investigated, and its limitations are examined.
4. Application
The derived mode shape expansion and correction approaches are considered in a plane truss structure in Fig. 1. The nodal points and the members are numbered in the figure. Each corresponding pair of nodal displacement components $({u}_{i},{v}_{i})$ is expressed by a set of forces $\left({H}_{i},{V}_{i}\right)\text{.}$ The structure has nine DOFs, except for the boundary conditions: ${\left[\begin{array}{ccccccccc}{u}_{2}& {v}_{2}& {u}_{3}& {v}_{3}& {u}_{4}& {v}_{4}& {u}_{5}& {v}_{5}& {u}_{6}\end{array}\right]}^{T}$. Each member has an elastic modulus of 200 GPa, a crosssectional area of $\text{2.5\xd7}{\text{10}}^{\text{3}}{\text{m}}^{\text{2}}$ and unit mass per unit volume of $\text{7.860}{\text{kg/m}}^{\text{3}}$.
For this numerical experiment, the actual stiffnesses of the truss structure are assumed as follows:
where ${k}_{i}^{\text{'}}$ denotes the $i$th analytical stiffness. The stiffness varies from 1.5 % to 37.5 % of the theoretical stiffness. The mass is assumed as invariant.
The measurement was performed at five different DOFs corresponding to ${\left[\begin{array}{ccccc}{v}_{2}& {v}_{3}& {u}_{4}& {u}_{5}& {v}_{5}\end{array}\right]}^{T}\text{.}$ This numerical example considered the following three different measurement cases with varying natural frequencies and corresponding mode shape data: a) Case 1: only the first natural frequency and mode, b) Case 2: the first two natural frequencies and corresponding modes, and c) Case 3: the first three natural frequencies and corresponding modes. Using and expanding the measured natural frequencies and corresponding mode shape data for each case, the stiffness and mass matrices were estimated utilizing Eq. (9). Using the predicted parameter matrices, the eigenvalues were reevaluated and compared as shown in Table 1. The eigenvalues measured in each case were identical with those predicted using the estimated parameter matrices. The eigenvalue estimates were as follows:
a) Case 1: eigenvalue of $\text{0.0680\xd7}{\text{10}}^{\text{6}}{\text{(rad./sec.)}}^{\text{2}}$;
b) Case 2: eigenvalues of $\text{0.0680\xd7}{\text{10}}^{\text{6}}{\text{}\text{(rad./sec.)}}^{\text{2}}$ and $\text{0.2633\xd7}{\text{10}}^{\text{6}}{\text{}\text{(rad./sec.)}}^{\text{2}}$;
c) Case 3: eigenvalues of $\text{0.0680\xd7}{\text{10}}^{\text{6}}{\text{}\text{(rad./sec.)}}^{\text{2}}$, $\text{0.2633\xd7}{\text{10}}^{\text{6}}{\text{(rad./sec.)}}^{\text{2}}$ and
$\text{0.6189\xd7}{\text{10}}^{\text{6}}{\text{(rad./sec.)}}^{\text{2}}$.
However, the other eigenvalues, except for the measured modes, were not consistent with the actual values. This indicates that accurate parameter matrices cannot be obtained from incomplete measurements.
Fig. 1. A planar truss structure
Table 1. Eigenvalues obtained analytically and experimentally (unit: $\text{\xd7}{\text{10}}^{\text{6}}{\left(\text{rad./sec.}\right)}^{2}$)
Mode number

1st

2nd

3rd

4th

5th

Before updating

0.0634

0.2422

0.5851

1.7548

2.8366

Actual truss

0.0680

0.2633

0.6189

1.9585

2.9804

Numerical results considering only the first mode

0.0680

0.2422

0.5851

1.7582

2.8315

Numerical results considering the first two modes

0.0680

0.2633

0.5870

1.7598

2.8200

Numerical results considering the first three modes

0.0680

0.2633

0.6189

1.7617

2.8158

Fig. 2 represents the difference between the estimated mode shapes and the actual mode shapes for each case. The mode shape data in this work were transformed as normalized data. The first mode shape for each case is estimated with a small margin of error as depicted in Fig. 2(a). However, the second and third modes exhibit greater inconsistencies than the first mode. It is recognized that this inconsistency is a result of the incomplete mode shape measurement data and the inaccurate prediction of parameter matrices.
Fig. 3 represents the inconsistency in the satisfaction of actual parameter matrices in the three cases. This numerical example was designed to estimate the stiffness matrix as established the above. However, the plots indicate that the mode shape calculations produced greater change in the mass matrix than in the stiffness matrix. A small change in the stiffness matrix is observed, contrary to expectations. The change of the analytical stiffness matrix is demonstrated by the variation in the mass matrix, yielding the prescribed eigenvalues and corresponding mode shapes. In all cases, the numerical values in the predicted parameter matrices are estimated with a similar margin of error.
Fig. 2. Difference between the estimated and actual normalized mode shape: a) Case 1, b) Case 2, c) Case 3
a)
b)
c)
This application demonstrates that it is not easy to accurately estimate the mode shape and the parameter matrices if insufficient measurement data are provided. Although the natural frequencies and corresponding mode shapes are obtained by the proposed method, the change in the analytical stiffness matrix is mostly lost in the variation of the mass matrix.
From the application, it is understood that it’s not easy to accurately estimate the mode shape as well as the parameter matrices if the insufficient measurement data are given. Although the natural frequencies and corresponding mode shapes are obtained by the proposed method, the change of the analytical stiffness matrix is mostly pervaded into the variation of the mass matrix.
5. Conclusions
In this study, the expansion of incomplete mode shape measurement data and the estimation of the parameter matrices were derived. The first mode shape is estimated with a small margin of error regardless of the number of measured mode shapes. The second and third modes exhibit greater inconsistencies than the first mode. The method does not provide the accurate mode shapes for a full set of DOFs and parameter matrices because of the incomplete measurement of mode shapes. The change of the analytical stiffness matrix is mostly lost in the variation of the mass matrix by the proposed method.
Fig. 3. Inconsistency between estimated and actual parameter matrices: a) stiffness variation in Case 1, b) mass variation in Case 1, c) stiffness variation in Case 2, d) mass variation in Case 2, e) stiffness variation in Case 3, f) mass variation in Case 3
a)
b)
c)
d)
e)
f)
Acknowledgements
This work was supported by the National Research Foundation of Korea (NRF) Grant funded by the Korea Government (MEST) (No. 20110012164).
References
 Mottershead J. E., Stanway R. Identification of structural vibration parameters by using a frequency domain filter. Journal of Sound and Vibration, Vol. 109, Issue 3, 1986, p. 495506. [Search CrossRef]
 Friswell M. I., Penny J. E. T. Updating model parameters from frequency domain data via reduced order models. Mechanical Systems and Signal Processing, Vol. 4, Issue 5, 1990, p. 377391. [Search CrossRef]
 Lee U., Shin J. A frequency response functionbased structural damage identification method. Computers and Structures, Vol. 80, Issue 2, 2002, p. 117132. [Search CrossRef]
 Rahmatalla S., Eun H. C., Lee E. T. Damage detection from the variation of parameter matrices estimated by incomplete FRF data. Smart Structures and Systems, Vol. 9, 2012, p. 5570. [Search CrossRef]
 O’Callahan J. A procedure for an improved reduced system (IRS) model. Proceedings of the 7th IMAC, Las Vegas, Nevada, USA, 1989, p. 1721. [Search CrossRef]
 Kenneth F. A., Francois H. Dynamic mode shape expansion using mass orthogonality constraints. Proceedings of the 18th IMAC, Texas, USA, 2000, p. 14961502. [Search CrossRef]
 Ewins D. J. Adjustment or updating of models. Sādhanā, Vol. 25, 2000, p. 235245. [Search CrossRef]
 Avitabile P., O’Callahan J. Dynamic expansion of FRFs for the full FRF matrix. Proceedings of the 19th IMAC, Orlando, Florida, USA, 2001. [Search CrossRef]
 Avitabile P., O’Callahan J. Frequency response function expansion for unmeasured translation and rotation dofs for impedance modeling applications. Mechanical Systems and Signal Processing, Vol. 17, Issue 4, 2003, p. 723745. [Search CrossRef]
 Yang Q. W., Liu J. K. Structural damage identification based on residual force vector. Journal of Sound and Vibration, Vol. 305, Issue 12, 2007, p. 298307. [Search CrossRef]
 Zhao J., Zhang L. A method for structural damage identification using residual force vector and mode shape expansion. International Conference on Multimedia Technology, Hangzhou, China, 2011, p. 945949. [Search CrossRef]
 Yun G. J., Ogorzalek K. A., Dyke S. J., Song W. A parameter subset selection method using residual force vector for detecting multiple damage locations. Structural Control Health Monitoring, Vol. 17, Issue 1, 2010, p. 4867. [Search CrossRef]
 Gafka G. K., Zimmerman D. C. Structural damage detection via least squares dynamic residual force minimization with quadratic measurement error inequality constraint. Proceedings of the 14th International Modal Analysis Conference, 1996, p. 12781284. [Search CrossRef]
 Friswell M. I., Mottershead J. E. Finite element model updating in structural dynamics. Solid Mechanics and its Applications, Kluwer Academic Publishers Group, 1995, p. 1286. [Search CrossRef]
 Baruch M., BarItzhack I. Y. Optimal weighted orthogonalization of measured modes. AIAA Journal, Vol. 17, Issue 8, 1979, p. 927928. [Search CrossRef]
 Baruch M. Optimal correction of mass and stiffness matrices using measured modes. AIAA Journal, Vol. 20, Issue 11, 1982, p. 16231626. [Search CrossRef]
 Berman A. Mass matrix correction using an incomplete set of measured modes. AIAA Journal, Vol. 17, Issue 10, 1979, p. 11471148. [Search CrossRef]
 Berman A., Nagy E. J. Improvement of a large analytical model using test data. AIAA Journal, Vol. 21, Issue 8, 1983, p. 11681173. [Search CrossRef]
 Kabe A. M. Stiffness matrix adjustment using mode data. AIAA Journal, Vol. 23, Issue 9, 1985, p. 14311436. [Search CrossRef]
 Caeser B., Pete J. Direct update of dynamic mathematical models from modal test data. AIAA Journal, Vol. 25, Issue 11, 1987, p. 14941499. [Search CrossRef]
 Zimmerman D. C., Widengren M. Correcting finite element models using a symmetric eigenstructure assignment technique. AIAA Journal, Vol. 28, Issue 9, 1990, p. 16701676. [Search CrossRef]
 Wei F. S. Stiffness matrix correction from incomplete test data. AIAA Journal, Vol. 18, Issue 10, 1980, p. 12741275. [Search CrossRef]
 Wei F. S. Analytical dynamic model improvement using vibration test data. AIAA Journal, Vol. 28, Issue 1, 1990, p. 175177. [Search CrossRef]
 Wei F. S. Mass and stiffness interaction effects in analytical model modification. AIAA Journal, Vol. 28, Issue 9, 1990, p. 16861688. [Search CrossRef]
 Lee E. T., Eun H. C. Update of corrected stiffness and mass matrices based on measured dynamic modal data. Applied Mathematical Modelling, Vol. 33, Issue 5, 2009, p. 22742281. [Search CrossRef]
 Lee E. T., Eun H. C. Correction of stiffness and mass matrices utilizing simulated measured modal data. Applied Mathematical Modelling, Vol. 33, Issue 6, 2009, p. 27232729. [Search CrossRef]
 Lee E. T., Rahmatalla S., Eun H. C. Estimation of parameter matrices based on measured data. Applied Mathematical Modelling, Vol. 35, Issue 10, 2011, p. 48164823. [Search CrossRef]
 Eun H. C., Lee E. T., Chung H. S. On the static analysis of constrained structural systems. Canadian Journal of Civil Engineering, Vol. 31, Issue 6, 2004, p. 11191122. [Search CrossRef]