A two-step approach for damage detection in beam based on influence line and bird mating optimizer

Zhong Rong Lu2 , Dong Tan1

2, 1Department of Applied Mechanics, Sun Yat-sen University, Guangzhou, Guangdong Province, P. R. China

2Corresponding author

Vibroengineering PROCEDIA, Vol. 14, 2017, p. 102-107. https://doi.org/10.21595/vp.2017.19211
Received 20 September 2017; accepted 27 September 2017; published 21 October 2017

Copyright © 2017 - JVE International Ltd.


This paper presents a two-step approach for structural damage identification in beam structure using the influence line and bird mating optimizer (BMO). Local damage is simulated as the reduction of the elemental Young’s modulus and mass of beam element. The technique for damage localization based on influence line and its derivatives before and after damage for beam structure was outlined. An objective function comprised of dynamic acceleration is utilized for BMO algorithm. The dynamic response data under external force is calculated by Newmark integration method. Numerical examples of a simply supported beam was investigated. Effect of measurement noise is studied. Studies in the paper indicate that the proposed method is efficient and robust for identifying damages in beam structures.

Keywords: damage identification, influence line residue, bird mating optimizer, dynamic responses.

1. Introduction

It is well-known that structures are suffering from all kinds of damage problems, thus a reliable structural damage identification method is needed.

Vibration based damage identification basically comprised of methods using data from frequency domain and time domain. In the last few decades, many damage identification methods have been developed using the vibration data in frequency domain [1-4]. Moreover, with the advantage of massive vibration data, methods in time domain [5-8] emerged serve as a more reliable but complicated aspect in damage identification.

In mathematical aspect, damage identification problem can be regarded as an optimization problem that aims to find the optimal solution of an objective. And with the fast development of computation, more attention was paid in algorithm these years. Neural network was among the first to be adopted to identified damage [9, 10], and therefore others heuristic algorithm such as genetic algorithm (GA) [11], particle swarm optimization (PSO) [12] and artificial bee colony (ABC) [13] are proved to the available in damage identification.

A new heuristic algorithm named Bird Mating Optimizer (BMO) was proposed by Askarzadeh [14, 15] recently, which imitates the mating behavior of birds. The algorithm was simulated and compared with GA, PSO and GSO [14] and had better accuracy.

In this paper, a two-step method for damage identification is proposed, damage is quantified as both stiffness and mass reduction in different levels. The disadvantage of heuristic algorithm is the probability of local optimum, especially with large amount of identified parameters existed. Here the displacement influence line will be utilized to determine the damage location so as to reduce the number of parameters. Then, the damage is quantified in the second step using bird mating optimizer. An objective function is established by minimizing the discrepancies between the simulated ‘measured’ dynamic acceleration responses and calculated ones. A simulation of simply supported beam is studied to show the promising results of the method.

2. Damage identification methodology

2.1. Damage model construction

The finite element method is utilized to calculate the responses of the damage beams. Without loss of generality, we quantify the local damage as the reduction of both stiffness and mass parameters in damage element(s). For a structure with n elements, we preset n-dimension vectors α and β, whose elements αi and βi(i=1, 2,…, n) are between 0 and 1, stand for the damage parameters of stiffness and mass of the ith element, respectively:

K d = i = 1 n α i k i e ,         Μ d = i = 1 n β i m i e ,

where kie and mie represent the elemental stiffness and mass matrices for element i, and Kd, Md as global stiffness and mass matrices. The dynamic responses of the damaged structure can be obtained by Eq. (2) using direct integration method:

M d ¨ + C d ˙ + K d = F t ,

where C is the damping matrix of the system and Rayleigh damping is adopted that Cd=η1Md+η2Kd, η1 and η2 are two constants determined by two given damping ratios that corresponding to the first two modal frequencies of the structure.

2.2. Localization of the damage by influence line residue

The localization of damage is totally based on the static displacement of the beam when a external force exerted on it. In a beam structure with length l, if loaded with force P, the static deflection f(xs) of point xs can be obtained by Eq. (3):

f x s = 0 l M P η M ¯ x s , η E I η d η ,

where MP(η) is the bending moment of beam when only external force P exerted on the beam, and M¯(xs,η) is the bending moment when a unit force exerts at the measuring point xs, respectively, f is the deflection function of the beam.

If we assume that there is damage existing between a and b with stiffness damage index α, then the residue between damaged and intact situation is:

Δ f x s = a b M p η M ¯ x s , η E I 1 α - 1 d η .

If we set the external force P as a single moving load whose location is xp, then the displacement influence line residue for point xs is formed:

Δ f x p , x s = a b M p x p , η M ¯ x s , η E I 1 α - 1 d η .

The value of Mp(xp,η) within [a,b] will change with xp moving from a to b. In mechanics aspect, the change is produce directly to the change of shear force. Taking the partial derivative of Mp(xp,x) with respective to xp of Eq. (5) we get Eq. (6):

Δ f x p x p , x s = a b F s p x p , η M ¯ x s , η E I 1 α - 1 d η ,

where Fsp(xp,η) represents the shear force of the beam caused by moving load. A sudden change of Δfxp(xs,xp) will occur within [a,b], which can be explained as: when a moving force passes through a damaged area, the shear force within the area would change from positive to negative as is illustrated in Fig. 1.

Fig. 1. The change of shear force Δfxp(xs,xp)

For a structure with multiple damages, the values of Δfxp(xs,xp) will be the linear sum of derivative of influence line residue (DILR) of every damage individual. For discretized model with n elements and elemental length is l/n, the discrete influence line residue is:

D I L R x p = Δ f x p x p , x s = i = 1 N D F s p x p , i M ¯ x s , i l n E I 1 α i - 1 ,

where ND is the number of damages. One can further specify the change of DILR by taking the derivative of DILR(DDILR) to xp again. In the discretized model, the second-order central difference quotient is used to calculate the DDILR as shown in Eq. (8), where the damage parts will present as non-zero DDILR values:

D D I L R n i = D I L R n i - D I L R n i - 1 l e n g = Δ f n i + 1 - 2 Δ f n i + Δ f n i - 1 l e n g 2 .

Note that in numerical study, a damage is assumed as an elemental damage that associates with two nodes, on the other hand Δfxp is nodal value, so a single damage will be presented as two consecutive changes of corresponding nodes on the curve.

2.3. Objective function for identification assessment

To get results of higher precision especially of structure with lots of elements, the acceleration responses was to made up the objective function. For node n the acceleration responses vector with m sampled time points can be expressed as:

a c c n = a n t 1 a n t 2 a n t m T .

By measuring the deviance of acceleration responses between actual and calculated in the form of Modal Assurance Criterion (MAC):

a c c M A C n = ( { a c c n C } T { a c c n M } ) 2 { a c c n C } 2 { a c c n M } 2 ,

where the superscript C and M denote that the responses originated from calculation and measurement respectively. When the calculation in accordance with the measured data, the accMAC reaches 1, so the objective function can be:

f = j = 1 N A w ϕ j 2 1 - ( { a c c j C } T { a c c j M } ) 2 { a c c j C } 2 { a c c j M } 2 ,

where wϕ is the weighting factor of the acceleration data and NA is the number of the measuring points.

3. Numerical simulation

The identification procedure is shown in. A simply supported beam made of aluminum is the numerical example of the study, with length l= 10 m and cross-sectional width b= 0.6 m and height h= 0.4 m, the beam is discretized into 12 Euler-Bernoulli elements, Young’s modulus is E= 6.9 GPa and mass density ρ= 2700 kg/m3. There are three designed damages, elements 3, 8 and 9 with Young’s modulus reduction by 10 %, 20 % and 15 %, mass reduction by 5 %, 10 % and 5 % respectively. The damping model is Rayleigh damping and the two Rayleigh coefficients are both assumed to be 0.01.

To maintain practicality, 3 % random noise presented by Eq. (15) is added to the displacement values during localization:

f m = f 1 + μ 2 r a n d - 1 ,

where μ is the noise level, two sensors are set in nodes 4 and 8, and the localization is processed three times to reduce the affection of noise. The DDILR are shown in Fig. 2, nodes 3, 4, 8, 9, 10 show stable changes in the chart of sensor set in node 4, and nodes 8, 9 and 10 have constant trend in the chart of the other sensor in node 8. Thus, the corresponding elements 3, 8 and 9 are selected as suspected elements.

Fig. 2. DDILR of the simply supported beam

For acceleration responses a, a Gaussian distributed random noise Noise with zero mean and unit standard deviation is added as Eq. (16), where std() stands for the standard deviation of the acceleration responses in time history. The weighting factors wϕ of the accMAC are all set as one:

a = a c a l + E p N o i s e s t d a c a l .

The acceleration responses of the structure of a constant sinusoidal force acting at node 7 in the global z direction with P(t)=103sin(10πt)N is calculated by Newmark-beta method, the time increment is 0.01 second and the time duration for the response calculation is 6.0 seconds and 600 time steps in total. The arrangement of acceleration measurements are two nodes locating at nodes 5, 8.

The identified results of noise level 5 % and 10 % are shown in Fig. 3, with comparison to the result of noise free in condition. In noise free condition the identification yielded good results with maximum relative error 0.9667 % by elemental stiffness of the 3rd element, in low noise level the result is still satisfactory with maximum deviation less than 5 %, and in high noise level the maximum relative error is 5.2526 % which is not too great compared to the result of 5 %.

In the best logarithmic fitness values chart show in Fig. 4, it can be observed that with higher noise level, the searching process reach convergence faster but on the other hand the precision dropped.

Fig. 3. Comparison of identification results under different noise levels

Fig. 4. Comparison of logarithmic best fitness under different noise levels

4. Conclusions

By making use of the influence line residue and acceleration data, a two-step damage detection method for beam structures is proposed in this study. A damage model combined with stiffness and mass reduction is applied. The influence line data of displacement is used to reduce the dimension of the identified parameters, objective function based on acceleration and the form of MAC (Modal Assurance Criterion) is established, then BMO is used to minimizing the discrepancies between the simulated measured data and the data from damaged structure. A simply supported beam numerical simulation reflect the effectiveness of the proposed method, results under three different noise levels demonstrated that the proposed method is insensitive to measurement noise, which also manifested the robustness of the method.


  1. Cawley P., Adams R. D. The location of defects in structures from measurements of natural frequencies. The Journal of Strain Analysis for Engineering Design, Vol. 14, Issue 2, 1979, p. 49-57.
  2. Pandey A. K., Biswas M., Samman M. M. Damage detection from changes in curvature mode shapes. Journal of Sound Vibration, Vol. 145, Issue 2, 1991, p. 321-332.
  3. Lee J. Identification of multiple cracks in a beam using natural frequencies. Journal of Sound and Vibration, Vol. 320, Issue 3, 2009, p. 482-490.
  4. Zhong S., Oyadiji S. O. Detection of cracks in simply-supported beams by continuous wavelet transform of reconstructed modal data. Computers and Structures, Vol. 89, Issues 1-2, 2011, p. 127-148.
  5. Majumder L., Manohar C. S. A time-domain approach for damage detection in beam structures using vibration data with a moving oscillator as an excitation source. Journal of Sound and Vibration, Vol. 268, Issue 4, 2003, p. 699-716.
  6. Zhu X. Q., Law S. S. Wavelet-based crack identification of bridge beam from operational deflection time history. International Journal of Solids and Structures, Vol. 43, Issues 7-8, 2006, p. 2299-2317.
  7. Law S. S., Lu Z. R. Crack identification in beam from dynamic responses. Journal of Sound and Vibration, Vol. 285, Issues 4-5, 2005, p. 967-987.
  8. Lu X. B., Liu J. K., Lu Z. R. A two-step approach for crack identification in beam. Journal of Sound and Vibration, Vol. 332, Issue 2, 2013, p. 282-293.
  9. Wu X., Ghaboussi J., Jr J. H. G. Use of neural networks in detection of structural damage. British Journal of Surgery, Vol. 81, Issue 11, 1992, p. 578-581.
  10. Rhim J., Lee S. W. A neural network approach for damage detection and identification of structures. Computational Mechanics, Vol. 16, Issue 6, 1995, p. 437-443.
  11. Maity D., Tripathy R. R. Damage assessment of structures from changes in natural frequencies using genetic algorithm. Structural Engineering and Mechanics, Vol. 19, Issue 1, 2005, p. 21-42.
  12. Seyedpoor S. M. A two stage method for structural damage detection using a modal strain energy based index and particle swarm optimization. International Journal of Non-Linear Mechanics, Vol. 47, Issue 1, 2012, p. 1-8.
  13. Xu H. J., Ding Z. H., Lu Z. R., Liu J. K. Structural damage detection based on chaotic artificial bee colony algorithm. Structural Engineering and Mechanics, Vol. 55, Issue 6, 2015, p. 1223-1239.
  14. Askarzadeh A. Bird mating optimizer: an optimization algorithm inspired by bird mating strategies. Communications in Nonlinear Science and Numerical Simulation, Vol. 19, Issue 4, 2014, p. 1213-1228.
  15. Askarzadeh A. Parameter estimation of fuel cell polarization curve using BMO algorithm. International Journal of Hydrogen Energy, Vol. 38, Issue 35, 2013, p. 15405-15413.