# Comparison of Damage Detection Methods Depending on FRFs within Specified Frequency Ranges

Authors: Yong-Su Kim; Hee-Chang Eun
Journal: Advances in Materials Science and Engineering
 (2017)

---

## Abstract

Structural damage can be detected using frequency response function (FRF) measured by an impact and the corresponding responses. The change in the mechanical properties of dynamic system for damage detection can seldom be estimated using FRF data extracted from a very limited frequency range. Proper orthogonal modes (POMs) from the FRFs extracted in given frequency ranges and their modified forms can be utilized as damage indices to detect damage. The POM-based damage detection methods must be sensitive to the selected FRFs. This work compares the effectiveness of the damage detection approaches taking the POMs estimated by the FRFs within five different frequency ranges including resonance frequency and antiresonance frequency. It is shown from a numerical example that the POMs extracted from the FRFs within antiresonance frequency ranges provide more explicit information on the damage locations than the ones within resonance frequency ranges.


---

## Body

## 1. Introduction

Structural damage is detected based on the variation in dynamic responses due to the local change of physical properties at damage region. There have been many attempts to provide the accurate damage detection methods related to structural health monitoring. Many indices like mode shapes, mode shape curvature, frequency response function (FRF) curvature, modal strain energy, and so forth to evaluate the structural performance have been utilized [1, 2].One of the measurement data types is FRF data set. The FRF can be estimated by impact hammer test and the energy is propagated from the impact point. The FRF data measured from experimental work are utilized to estimate the characteristics of dynamic system. It indicates that the measured FRF data set can be used as an index to detect damage. Phani and Woodhouse [3] proposed the FRF curvature method based on only the measured data without the need for any modal identification. Lee and Kim [4] proposed a structural damage detection method to determine both location and magnitude of damage from perturbation equations of FRF data. Fritzen [5] evaluated a FRF-based model updating algorithm using experimentally collected data and presented a protocol for measurement selection and a regularization technique. Rahmatalla et al. [6] presented a feasible method for structural vibration-based health monitoring to reduce the dimension of the initial FRF data and to employ artificial neural network. Sampaio et al. [7] provided a structural damage identification technique using changes of the FRF to be related to the changes of the stiffness and mass through damage sensitivity equations. Wang et al. [8] developed a methodology by coupled response monitoring through the emergence of peaks on a FRF plot to determine the damage existence. Garcia-Palencia et al. [9] exhibited that the curvatures of FRF at frequencies other than natural frequencies can be used for identifying both the existence of damage and the location of damage. Bandara et al. [10] presented a damage identification method using the changes in the distribution of the compliance of the structure due to damage. Bandara et al. [11] investigated the advantages and limitations to use vibration-based damage detection methods from the measurements of mode shapes and FRFs. Nozarian and Esfandiari [12] identified parameter matrices using FRF data measured at specific positions and constraints and provided a damage detection method from their variations.Proper orthogonal decomposition (POD) analysis captures most of the kinetic energy in the least number of modes possible. Its application is similar to that of Fourier analysis, except that it normally requires far less modes to represent the system within a desired level. POMs extracted from the FRF data in a prescribed frequency range are utilized as an index to recognize the existence of damage. The POD method has been widely applied in various fields of engineering and science. Ramanamurthy et al. [13] related the POM to normal modes of vibration in systems for lightly modally damped systems. They exhibited that the POMs represent the principal axes of inertia formed by the distribution of data on the modal coordinate curve. Using statistical process control technique, Mondal et al. [14] investigated the applications of principal component analysis for damage diagnosis. Choi et al. [15] provided a damage detection algorithm based on the POD technique to filter out the influence of operational/environmental variation in the dynamic response. De Medeiros et al. [16] applied a POD statistical method to the field of damage detection for continuous static monitoring systems. Feeny and Kappagantu [17] presented a damage detection method based on fuzzy c-means clustering algorithm and measured FRF data reduced by principal component analysis. De Boe and Golinval [18] developed a POMs-based damage identification approach using acceleration dynamic responses of shear-type buildings. Shane and Jha [19] presented a neural networks-based damage detection method using FRF data. The method reduces the dimension of the initial FRF data and employs an artificial neural network for the actual damage.The parameter matrices as well as damage region of dynamic systems can be predicted using FRF data sets within a specific frequency range rather than at a specific frequency. The FRF data sets are transformed to the POM to represent the principal axes of inertia formed by the distribution of data on the modal coordinate curve. The POM can be changed depending on the FRF data sets extracted from full sets of FRFs. Beginning with the FRFs measured from the dynamic finite element model, this work investigates the sensitivity of the damage detection method using the POMs to be estimated from the FRFs corresponding to five different frequency ranges including resonance frequency and antiresonance frequency. A numerical example compares the sensitivity of the damage detection method depending on the extracted FRFs.


## 2. FRF and POM

The dynamic behavior of a structure, which is assumed to be linear and approximately discretized forn DOFs, can be described by the equations of motion:(1)Mu¨+Cu˙+Ku=Ft,whereM,K, andC denote the n×n analytical mass, stiffness, and damping matrices, respectively, u=u1u2⋯unT, andFt is the n×1 load excitation vector.Frequency response functionHΩ is defined as the ratio of the complex spectrum of the response to the complex spectrum of the excitation:(2)HΩ=UΩFΩ,where UΩ and FΩ denote the complex spectrum of the response and the complex spectrum of the excitation, respectively. The magnitude of the FRF HΩ denotes the ratio of UΩ/FΩ.For the case of a displacement response at stationp and a disturbing force at station q, the numerical frequency response can be constructed as(3)HpqΩ=∑i=1nϕipϕiqωi2-Ω2+2jξiωiΩ,where ϕip denotes the pth element of the vector ϕi and ωi and ξi denote the circular natural frequency and the damping ratio, respectively, for the ith mode, j=-1.The FRF data,H^pq, contaminated by external noise rarely provide the accurate information on the dynamic system. In this study, they can be described by adding a series of random numbers to the calculated FRF response data expressed by (4)H^pq=Hpq1+ασpq,where α denotes the relative magnitude of the error, and σpq is a random number variant in the range -1,1.The measured FRFs can be reduced by the POD and are transformed to the POD to extract extremal data set. The POD technique is effective method because basis elements are formed in an optimal way. The POD basis collects snapshots. The FRFs of a system are generated by forcing the system.HpqΩ expresses an FRF on a finite number of points in space and a finite frequency interval Ωmin≤Ω≤Ωmax. The FRFs at n dofs are sampled m frequencies and the data are arranged in a snapshot matrix H^:(5)H^=H^1qΩ1H^1qΩ2⋯H^1qΩmH^2qΩ1H^2qΩ2⋯H^2qΩm⋮⋮⋱⋮H^nqΩ1H^nqΩ2⋯H^nqΩmn×m.Let them×m autocovariance matrixCbe defined as(6)C=HmeaTHmea,where m is the number of the FRF data sets extracted. MatrixC is a Hermitian positive semidefinite matrix to possess a complete set of orthogonal eigenvectors with corresponding nonnegative real eigenvalues. Solving the eigenvalue problem of (6) at the core of the POD method, (6) satisfies(7)Cφk=λkφk,k=1,…,m,where the eigenvalues reflect the energies in different POMs and are arranged in descending order as follows:(8)λ1≥λ2≥⋯≥λm≥0,where the eigenvalues λk are the POVs, and the corresponding eigenvector φk of the extreme value problem is associated with a POV λk. The greatest POV is the optimal vector to characterize the snapshots.The POMs may be used as a basis for the decomposition ofHmea. The POM associated with the greatest POV is the optimal vector. If the eigenvalues are normalized to unit magnitude, then they represent the relative energy captured by the corresponding POM. The eigenvalue reflects the relative kinetic energy associated with the corresponding mode. The energy is defined as the sum of the POVs. The POMs are written as(9)ψk=H^φkH^φk,k=1,…,m,where H^ denotes the m×m FRF matrix in (5) and φk represents the m×1 eigenvector corresponding to the eigenvalue λk. The POMs ψ are arranged as follows:(10)ψ=ψ1ψ2⋯ψm.The slope of the POM data corresponding to two adjacent nodes from Figure1 can be defined as(11)θj=ψj+11-ψj1h,j=1,2,…,n-1,where h is the distance between two adjacent nodes. The sum of the squared slope of the POM at two adjacent nodes is defined as(12)ηj=θj2+θj+12,j=1,2,…,n-2.Figure 1
Slope of POM.Damage introduced into the flexural beam to be modeled as finite elements also leads to local changes in the shape of POM curvature obtained by POMs. The curvature at each locationj on the structure, ψj′′, is numerically obtained by a central difference approximation:(13)ψj′′=ψj-1-2ψj+ψj+1h2,j=2,3,…,n-1,where ψj′′ is the second derivatives at the jth node and h is the distance between two successive nodes. The damage is evaluated by the POM curvature extracted from the measured FRFs without the baseline data.


## 3. Numerical Experiment

A numerical experiment was performed in detecting the damage in a finite element model of a three-span continuous beam shown in Figure2 to examine the validity of the proposed method depending on the selected FRF data sets. The nodal points and the elements are numbered as shown in the figure. The continuous beam with a length of 9 m is modeled using 100 beam elements. The beam has an elastic modulus of 2.0×105 MPa and a unit mass of 7,860kg/m3. The beam’s gross cross section is 500 mm × 200 mm, and its damage section is established as 90% of the moment of inertia of the undamaged beam. The damping matrix is assumed as a Rayleigh damping to be expressed by the stiffness matrix and a proportionality constant of 0.0001. This application considers the damage detection of a beam with multiple damage at elements  and . The FRF data are numerically simulated by the responses at all nodes due to the impact at node 18. The measurement data in (4) include 1% external noise.Figure 2
A three-span continuous beam structure model.The damage leads to the local stiffness deterioration and the variation of the natural frequencies. Its existence can be recognized through the comparison with the natural frequencies at the undamaged state. The first three natural frequencies of the damaged dynamic system exist at 24.67 Hz, 40.89 Hz, and 69.69 Hz, respectively. Figure3(a) represents the mode shapes corresponding to the first, second, and third natural frequencies, respectively. The damage can be detected by the comparison of the mode shapes at both states. However, it is not easy to collect the accurate mode shapes corresponding to the resonance frequencies due to external noise. Thus this work utilizes a set of FRFs and they are transformed to the POM to take extremal data set. Figure 3(b) displays the FRF receptance curve of the responses at all nodes due to a single impact at node 18. It is shown that the first three resonance frequencies coincide with the natural frequencies. The resonance frequencies correspond to the peaks in the plots and the antiresonance frequencies correspond to the valleys.Figure 3
Dynamic characteristics of the system: (a) mode shapes and (b) FRF curve.
(a)
(b)This numerical study investigates the sensitivity of the damage detection depending on the FRFs extracted in the neighborhood of two resonance frequency ranges at points (B) and (D) in Figure3(b), two antiresonance frequency ranges at points (C) and (E), and point (A) less than the first resonance frequency in Figure 3(b). The FRFs were calculated at the intervals of 0.02 Hz for this numerical experiment. Eleven FRF data sets simulated in the frequency range of ±0.1Hz from the resonance and antiresonance frequencies are extracted to establish the POMs as shown in Table 1.Table 1
Frequency ranges of the extracted FRFs.
Location
(A)
(B)
(C)
(D)
(E)
Selected frequencies
10 Hz± 0.1 Hz
24.67 Hz± 0.1 Hz
34.25 Hz± 0.1 Hz
40.89 Hz± 0.1 Hz
50.93 Hz± 0.1 Hz
Less than first resonance
Resonance
Antiresonance
Resonance
AntiresonanceFigure4 represents the normalized POM curve of beams 1, 2, and 3 corresponding to the first POV depending on the noise-free FRFs in the selected five frequency ranges shown in Figure 3(b). Figure 4(a) displays the POM curves of beam 1 depending on the FRFs to be chosen. It is observed that the POM plots corresponding to the first resonance frequency (B) and its neighboring frequency (A), and the second frequency (D) and its neighboring frequency (E) are very similar except the curve (C) corresponding to the antiresonance frequency. The valley of the curve corresponding to (C) corresponds to the impact location. Figure 4(b) represents the POM curves of beam 2. It is shown that the peak of the curves moves little by little to the right with the increase of the frequency region to extract the FRF data except the POM curve corresponding to region (E). Figure 4(c) exhibits the POM curves of beam 3. The peak of the curves moves little by little to the right with the increase of the frequency region to extract the FRF data. It is found that the plots in Figure 4 do not provide the information on the damage in the beam.Figure 4
POM curves using noise-free FRFs within the five different frequency ranges: (a) beam 1, (b) beam 2, and (c) beam 3.
(a)
(b)
(c)Figure5 exhibits the POM curvature curves approximated numerically using the central difference method of (11). The damage locates at the region to display the abrupt variation of the curvature curve. The impact at node 18 of beam 1 leads to the abrupt change of the POM curvature curve as shown in Figure 5(a). Curve (E) in Figure 4(b) looks like the second mode shape to take the absolute. The abrupt change in the neighborhood of the inflection point is observed. Except those variations, the local change in the neighborhood of elements  and  is observed. It is observed that the noise-free POM curvature curves taken from the FRF data of the first four regions shown in Table 1 properly indicate the damage regions.Figure 5
POM curvature curves using noise-free FRFs: (a) beam 1, (b) beam 2, and (c) beam 3.
(a)
(b)
(c)Figures6–10 display the POM curvature curves estimated from the FRF data corresponding to the five different frequency regions. The simulated FRF data contain 1% external noise in (4). The plots exhibit the large and small variations of the curvature due to the presence of the external noise so that the damage cannot be detected. The abrupt variation nearby the impact location at node 18 of beam 1 is observed in the plots. Figure 6 represents the plots corresponding to region (B) including the first resonance frequency. It is shown that the POM curve in Figure 6(a) corresponds to first mode shape. The other curves indicate explicitly the damage at the 37th element only. The other damage cannot be recognized because of the noise contained in the measurement data. It is understood that the POM extracted in the first resonance frequency range is sensitive to the external noise.Figure 6
POM curvature curves using FRFs in the region (A) containing 1% external noise: (a) beam 1, (b) beam 2, and (c) beam 3.
(a)
(b)
(c)Figure 7
POM curvature curves using FRFs in the region (B) containing 1% external noise: (a) beam 1, (b) beam 2, and (c) beam 3.
(a)
(b)
(c)Figure 8
POM curvature curves using FRFs in the region (C) containing 1% external noise: (a) beam 1, (b) beam 2, and (c) beam 3.
(a)
(b)
(c)Figure 9
POM curvature curves using FRFs in the region (D) containing 1% external noise: (a) beam 1, (b) beam 2, (c) beam 3.
(a)
(b)
(c)Figure 10
POM curvature curves using FRFs in the region (E) containing 1% external noise: (a) beam 1, (b) beam 2, and (c) beam 3.
(a)
(b)
(c)Figure7 represents the curves corresponding to region (C) including the second antiresonance frequency. It is found in Figure 7(a) that the POM curve after the impact node 33 does rarely respond and cannot obtain the damage information. The other plots in Figure 7 give the damage information at the 22nd element only and the other damage location cannot be recognized. It is observed that the POM extracted in this region is also sensitive to the external noise.Figure8 represents the plots corresponding to region (D) including the second resonance frequency. It is shown that the POM curve in Figure 8(a) corresponds to the absolute curve of the second mode shape. It is obtained from  Figures 8(b)-8(c) that the damage at the only 22nd element can be explicitly detected and the other damage should be recognized through close examination. The abrupt change of the plots in the center of the span is due to the result that the POM curve takes the absolute curve as shown in Figure 8(a). It is known that the POM extracted in this region is also sensitive to the external noise.Figure9 exhibits the plots corresponding to region (E) including the third antiresonance frequency. It is analyzed that the POM curve in Figure 9(a) cannot obtain any information on the damage location. It is shown that the curves in  Figures 9(b)-9(c) lead to the abrupt variation at the damage locations except impact node 33 and the damage can be explicitly detected. From this numerical experiment, the damage detection method can be explicitly and widely carried out by the three damage indices of the POM corresponding to the first POV to extract the FRFs in the first and third antiresonance frequencies despite the external noise. The POM-based approaches using the FRFs extracted in the other ranges should be sensitive to the external noise and can be utilized by close examination of the POM curve.


## 4. Conclusions

This work compared the damage detection methods utilizing three different damage indices depending on the extracted FRFs within five different frequency ranges including resonance frequency and antiresonance frequency. It is shown from a numerical example that the damage detection method can be explicitly and widely carried out by the three damage indices of the POM corresponding to the first POV to extract the FRFs in the first and third antiresonance frequencies despite the external noise. The POM-based approaches using the FRFs extracted in the other ranges should be sensitive to the external noise and can be utilized by close examination of the POM curve.



---

*Source: 5821835-2017-08-29.xml*