Damage detection in a fixed-fixed beam using natural frequency changes

Muhammad Wasif Khan1 , Naveed Akmal Din2 , Rizwan Ul Haq3

1, 2, 3National University of Sciences and Technology, Islamabad, Pakistan

1Corresponding author

Vibroengineering PROCEDIA, Vol. 30, 2020, p. 38-43. https://doi.org/10.21595/vp.2019.21081
Received 3 October 2019; accepted 11 October 2019; published 2 April 2020

Copyright © 2020 Muhammad Wasif Khan, et al. 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.
Creative Commons License
Table of Contents Download PDF References
Cite this article
Views 72
Reads 32
Downloads 433
CrossRef Citations 0

This paper presents a damage detection technique which uses change in Frequency Response Function (FRF) and Modal Strain Energy (MSE) for damage detection in beam and truss like structures. It’s a robust approach for Structural Health Monitoring (SHM) using natural frequency for structural damage assessment. This technique couples FRF with Iterative Modal Strain Energy (IMSE) method to accurately localize and quantify the damage with only few frequencies, overcoming the limitation of mode shape based damage detection methods which requires mode shapes at full coordinates but damage is either over/underestimated. In the current approach, a damage detection algorithm is developed and validated by conducting numerical studies for a Fixed-Fixed beam, both noise-free and noise-effected cases are simulated. The numerical studies reveal that proposed algorithm is capable of identifying, localizing and quantifying the damage accurately at reduce computational cost.

Keywords: damage, FRF, modal strain energy, detection, estimation, iterative.

1. Introduction

SHM finds its application in mechanical, civil and aerospace Engineering. Damage in any structure is result of deterioration, degradation or any natural disaster in a structure during its service life. Visual inspection methods for damage assessment are dependent on instrumentation and have limited capability to determine the damage extent and severity. Vibration based Damage Identification (VBDIT) techniques have significant importance in SHM problems [1-4]. Many structural damage detection techniques are proposed in the recent decades which use the measured mode shapes and natural frequencies in the damaged structure to localize and estimate the damage. Stubbs et al. first proposed Damage Index (DI) method which uses MSE for structural damage detection [5]. Modal Strain Energy Decomposition (MSED) method is an improved form of DI method to improve the localization ability of SDI method [6, 7]. Shi et al. introduced Modal Strain Energy Change (MSEC) method using sensitivity of Modal Strain Energy Change Ratio (MSECR) to apply and validate this approach on a real life 2D structure in laboratory [8, 9]. DI & MSEC methods require paired modes of intact and damaged structures, as correct modal identification is practically difficult in some cases. Considering this limitation, Cross Modal Strain Energy (CMSE) method is developed by hu et al. to quantify the damage accurately [10]. Hwang et al. used few measured FRF for damage identification in a cantilever beam [11]. Zenzen et al. coupled FRF with BAT algorithm, minimizing the objective function based on intact and damaged structure natural frequency [12]. Khatir et al. proposed a damage detection technique using change in natural frequency and Modal Scale Factor (MSF) to minimize an objective function through PSO and BAT algorithm [13]. Pan and Zhang et al. Formulated a Noise Response Rate (NRR) in an inverse problem to select the specific mode shape for damage detection based on mode shape sensitivity to noise. This method improves the accuracy for prediction of size and location of damage [14]. An Improved modal strain energy method is formulated and validated on same Fixed-Fixed beam with noise up to 5 % provides better convergence rate than MSEC method. It was recommended to extend this study to 3D complex structures [15]. Wange et al. proposed IMSE method and validated on a cantilever beam using experimentally measured frequencies [16]. In this research work a damage detection algorithm is proposed, coupling FRF with IMSE method for structural damage detection in a Fixed-Fixed beam.

2. Methodology

2.1. Frequency response function (FRF)

Vibration based damage identification techniques (VBDIT) use change in natural frequency and mode shaps for damage detection in structure. Basic dynamic response equation for n DOF system can be expressed as Eq. (1):

M a ¨ t +   C a ˙ t + K a t c = f t .

For an external force and displacement ft=Fωejωt and at=aωejωt, a damped free condition FRF is expressed as:

H ω = - ω 2 M + K - 1 .

The analytical and measured FRFs is presented as [Hω] and [Hω*], where:

K * = H + ω 2 M .

It is assumed that mass of structure remains constant and stiffness changes:

K = [ H ] - 1 - H - 1 * .

When multiplied by H*, Eq. (4) gives:

H * K = H * H - 1 - H - 1 * .

Based of analytical and measured FRF, beta (β) will be calculated as:

β 1 ,   i = ( [ H ] 1 n * ) *   H - 1 : , i - I 1 , i .

2.2. Iterative modal strain energy (IMSE) method

Damage severity will be estimated for damage location identified form FRF, since M*=M and Global stiffness matrix will be written as linear combination of local stiffness matrix for each element:

[ K ] * = K + n = 1 N d   α n K l n ,

where Nd is the total number of damaged elements while αnand ln show the damage severity coefficient and the damaged element. Damage severity equation will be as follows:

n = 1 N d α n Φ i t K l n Φ i * = λ i * λ i - 1 .

Structural MSE and elemental MSE changes due to damage which is applied here as an input for estimation of damage extent:

C i = Φ i t K Φ i * ,
C n , i = α n Φ i t K l n Φ i * .

Here Φi, Φi*, λi and λi* represent the mode shapes and natural frequencies for intact and damaged structure. Using the equation Eqs. (9, 10), Eq. (8) can be written as follow:

n = 1 N d α n C n , i = b i ,
b i = λ i * λ i - 1 .

For m Eq. (11) can be simplified as:

C α = b ,

where C, α and b present elemental MSE, damage severity coefficient and change in natural frequency. For mNd, Least square solution method will be used to calculate Damage severity coefficient:

α = ( C T C - 1 )   C T b .

Eq. (14) requires mode shapes for damaged structure at full coordinates which is difficult to obtain, here initially zero damage is assumed for Fixed-Fixed beam and IMSE method is applied to quantify the damage using measured natural frequencies:

Φ i * K * , M = Φ i * K + n = 1 N d   α n K l n   .

Eqs. (14, 15) will be used for each Iteration of IMSE Method. IMSE method consists of four steps for damage quantification.

Step 1: Initialize the solution with α0= 0, calculate Φi*(0)K*,M where K*=K.

Step 2: Solve for α using Φi*(0), first iteration for IMSE completes here.

Step 3: Compute Φi*(k-1) from αk-1, and estimate αk using Φi*(k-1), where k= 2, 3….

Step 4: if αk-αk-1<r damage severity is estimated, otherwise move to step 3, where r= 0.0001.

2.3. Damage detection algorithm

A damage detection algorithm is proposed based on FRF & IMSE method is shown in Fig. 1.

Fig. 1. Damage detection algorithm based on natural frequencies

 Damage detection algorithm based on natural frequencies

2.4. Noise-effect

In actual measured modal parameters differ by simulated FE parameters due to noise. In order to simulate this effect on dynamic response of structure, here k j* and kj  are stiffnesses with and without Gussian Noise for damaged elements localized. Here γj present the random number with standard deviation of 1 and of mean 0, n shows the percentage of noise:

k   j * = k j   1 + n γ j .

3. Numerical simulations

A Fixed-Fixed Beam of length 2 m with elastic modulus of 210e09 N/m2, area of 0.0014 m2 and density of 7850 Kg/m3 is considered for eigen analysis consisting of 12 Euler Bernoulli beam elements and 13 nodes with 2 DOF per node is shown in Fig. 2. Three Damage cases are simulated and first three natural frequencies are shown in Table 1.

Fig. 2. FE model of a fixed-fixed beam

 FE model of a fixed-fixed beam
 FE model of a fixed-fixed beam

Table 1. First three natural frequencies (Hz) for beam

Damage location pe (element)
Damage severity (percentage)
Natural frequency (Hz)
35 %
3rd and 8th
25 % and 30 %
3rd ,6th and 9th
40 % each

4. Results and discussions

4.1. Damage localization using FRF

FRF for each damage case indicates presence of damage in the structure. Damage localization indicator β1,i uses Eq. (6) for damage location at element and its corresponding DOF.

Fig. 3. Damage indicator (β) for DCI: a) based on element, b) based on DOF

 Damage indicator (β) for DCI: a) based on element, b) based on DOF


 Damage indicator (β) for DCI: a) based on element, b) based on DOF


4.2. Damage quantification using IMSE method

IMSE method utilizes the natural frequencies of Intact and damaged structures. For a multiple damage case it requires only first three natural frequencies for damage quantification.

Fig. 4. Damage indicator (β) for DCII: a) based on element, b) based on DOF

 Damage indicator (β) for DCII: a) based on element, b) based on DOF


 Damage indicator (β) for DCII: a) based on element, b) based on DOF


Fig. 5. Damage indicator (β) for DCIII: a) based on element, b) based on DO

 Damage indicator (β) for DCIII: a) based on element, b) based on DO


 Damage indicator (β) for DCIII: a) based on element, b) based on DO


For beam structure single and multiple damage cases FRFs indicate the presence of damage in all three damage cases. For DCI, β indicator shows a higher value at 6th element and its corresponding DOFs. Similarly, for DC II and DC III, there exists a higher value of β at 3rd, 6th, 8th and 9th Element. This damage localization is used as an input for damage quantification in the damaged elements. For DCI, IMSE Method uses only first natural frequency of damaged structure. The damage severity estimator α is estimated within 10 iterations with a tolerance of 0.0001, similarly DCII and DCIII require first two and three natural frequencies to estimate the damage severity within 9 and 11 iterations.

Fig. 6. Damage severity estimator (α) for DCI, DCII and DCIII

 Damage severity estimator (α) for DCI, DCII and DCIII
 Damage severity estimator (α) for DCI, DCII and DCIII
 Damage severity estimator (α) for DCI, DCII and DCIII

5. Conclusions

A damage detection algorithm is proposed which uses FRF as Damage indicator and uses its characteristics for damage localization. Damage severity is estimated by IMSE method which requires minimum input in form of few measured frequencies from damaged structure. The proposed approach is better than conventional mode shape based methods using only natural frequency of damaged structure as an input for damage detection. It requires only one sensor in impact hammer based modal testing setup, proving to be a cost effective approach with reduced computational cost and higher convergence rate. This method may be tested on a real structure in laboratory using actually measured natural frequencies.


  1. Montalvão D., Maia N. M. M., Ribeiro A. M. R. A review of vibration-based structural health monitoring with special emphasis on composite materials. The Shock and Vibration Digest, Vol. 38, 2006, p. 295-324. [Publisher]
  2. Yan Y. J., Cheng L., Wu Z. Y., et al. Development in vibration-based structural damage detection technique. Mechanical Systems and Signal Processing, Vol. 21, 2007, p. 2198-2211. [Publisher]
  3. Fan W., Qiao P. Z. Vibration-based damage identification methods: a review and comparative study. Structural Health Monitoring, Vol. 10, 2011, p. 83-111. [Publisher]
  4. Seo J., Hu J. W., Lee J. Summary review of structural health monitoring applications for highway bridges. Journal of Performance of Constructed Facilities, Vol. 30, Issue 4, 2016, https://doi.org/10.1061/(ASCE)CF.1943-5509.0000824. [Publisher]
  5. Stubbs N., Kim J., Topole K. An efficient and robust algorithm for damage localization in offshore platforms. Proceedings of the ASCE Structures Congress, 1992. [Search CrossRef]
  6. Yang H. Z., Li H. J., Wang S. Q. Damage localization of offshore platforms under ambient excitation. China Ocean Engineering, Vol. 17, 2003, p. 495-504. [Search CrossRef]
  7. Li H. J., Yang H. Z., Hu S. J. Modal strain energy decomposition method for damage localization in 3D frame structures. Journal of Engineering Mechanics, Vol. 132, 2006, p. 941-951. [Publisher]
  8. Shi Z. Y., Law S. S. Structural damage localization from modal strain energy change. Journal of Sound and Vibration, Vol. 218, 1998, p. 825-844. [Publisher]
  9. Shi Z. Y., Law S. S., Zhang L. M. Structural damage detection from modal strain energy change. Journal of Engineering Mechanics, Vol. 126, 2000, p. 1216-1223. [Publisher]
  10. Hu S. L. J., Li H. J., Wang S. Q. Cross-model cross-mode method for model updating. Mechanical Systems and Signal Processing, Vol. 21, 2007, p. 1690-1703. [Publisher]
  11. Hwang H. Y., Kim C. Damage detection in structures using a few frequency response measurements. Journal of Sound and Vibration, Vol. 270, 2004, p. 1-14. [Publisher]
  12. Roumaissa Zenzena, Idir Belaidia, Samir Khatir, Magd Abdel Wahab A damage identification technique for beam-like and truss structures based on FRF and bat algorithm. Comptes Rendus Mécanique, Vol. 346, Issue 12, 2018, p. 1253-1266. [Publisher]
  13. Samir Khatir, Idir Belaidi, Roger Serra, Magd Abdel Wahab, Khatir Tawfiq Numerical study for single and multiple damage detection and localization in beam-like structures using BAT algorithm. Journal of Vibroengineering, Vol. 18, Issue 1, 2016, p. 202-213. [Search CrossRef]
  14. Pan Jingwen, Zhang Zhifang, Wu Jiurong, Ramakrishnan Karthik Ram, Singh Hemant. A novel method of vibration modes selection for improving accuracy of frequency-based damage detection. Composites Part B: Engineering, Vol. 159, 2019, p. 437-446. [Publisher]
  15. Li Yingchao, Wang Shuqing, Zhang Min, Zheng Chunmei An improved modal strain energy method for damage detection in offshore platform structures. Journal of Marine Science and Application, Vol. 15, Issue 2, 2016, p. 182-192. [Publisher]
  16. Shuqing Wang Iterative modal strain energy method for damage severity estimation using frequency measurements. Journal of the International Association for structural control and Monitoring, Vol. 20, Issue 2, 2013, p. 230-240. [Search CrossRef]