Numerical study for single and multiple damage detection and localization in beam-like structures using BAT algorithm

This paper presents a new damage detection and localization technique based on the changes in vibration parameters using BAT and Particle Swarm Optimization algorithm. The finite element method is used to apply damage at specific element(s) of the considered beams. The damage is represented by a reduction in Young’s modulus and the identification of damage is formulated as an optimization problem using objective function based on Modal Scale Factor and changes in natural frequencies. A procedure for detecting and locating damage of beam-like structures based on BAT algorithm is used. This approach presents a method that can be used to detect the single and multiple-damage positions and the rate of damage in structural elements with high accuracy after the first iteration. The results obtained using BAT algorithm are compared to those obtained using Particle Swarm Optimization Algorithm. By taking noise into account in the damage detection and localization problem, it is shown that our approach based on BAT algorithm can detect the damage locations with high accuracy.


Introduction
In the recent years, a structural damage detection became a field that has attracted a lot of interest in the scientific community.There have been many studies and researches intending to find a reliable method to detect and localize damage in structural elements using optimization methods.BAT algorithm (BA) is recently introduced as an effective optimization algorithm and was successfully applied to several problems [1][2][3][4][5][6].Identifying structural damage using measured vibration data is an inverse approach in mathematics.The normal damage detection methods produce an objective function, which is defined in terms of the differences between the data identified by modal vibration tests and those calculated from analytical or finite element model [7,8].The genetic algorithm (GA) was used to produce a sufficiently optimized amplitude characteristic filter to extract damage information from the strain energy mode shapes [9].A finite element (FE) model was used to generate training data set with the known location.The changes in the natural frequencies together with FE models to find the damage was proposed [10].An uncertainty approach for damage detection was suggested [11].The coordinate modal assurance criterion (COMAC), which makes use of the numerical and experimental vibration modes, was used to determine the magnitude and position of damage [12].To validate the method, two structures were chosen and modelled using FEA.In order to detect damage throughout the whole structure, especially some large, complicated structures, a methodology called global vibration-based structural damage detection [13,14] has been proposed.The damage location and severity can be determined by finding the differences of structural dynamic characteristics between the intact and damaged structures [15].A few natural frequencies were used and a methodology to non-destructively locate and estimate the size of damage in structures was proposed [16].A detailed revision of the methods for damage detection can be found in [17].Damage detection and estimation in structures based on incomplete modal data of a damaged structure can be done via a pattern search algorithm.The objective function based on the condensed mass and stiffness Matrices was formulated [18].GAs were employed by many researchers for structural damage detection and localization based on model updating methods [19].A damage detection and localization technique in beam-like structure using MAC (modal assurance criterion), COMAC (coordinate modal assurance criterion) and LFCR (Local Frequency Change Ratio) was introduced [20].A genetic algorithm (GA) was used to detect damage in beam structure using three objective functions based on frequencies and eigenvectors [21].In order to allow a comparison between different techniques, the following classification is proposed [22]: Detection, Localization, Assessment and Consequence.
In the present work, a new damage identification method is used by applying BAT algorithm to detect and locate single-and multiple-damage in beam-like structures.The objective functions used in the optimization process are based on vibration data of the structure.A finite element analysis is performed on simply supported beams discretized into 10 elements.Through the results obtained, we noticed that this method can accurately identify the damage, even though noise is present.This paper is organized as follows.In Section 2, methodological approach to the damage detection and localization problems is detailed.In Section 3, we present a numerical model of the beam structure.Implementation of BAT and PSO algorithm for damage detection and localization is given in Section 4. Results and discussion are shown in Section 5.

Methodological approach to the damage detection and localization using BAT algorithm
The optimization method refers to solving a problem in the best possible way.In most cases, the goal is to minimize or maximize an objective function based vibration analysis by choosing values for its variables within an acceptable range.The damage detection and localization procedures, which we propose here in, are carried out according to the following steps: 1) Expressing the damage by reducing Young's modulus using FEM and calculating the natural frequencies and mode shapes (real damage).
2) Quantifying damage using BAT algorithm form the optimization problem.3) Solving the inverse problem after providing frequencies and mode shapes (real damage) to the objective function in order to calculate damage locations and severity using BAT algorithm and FEM (estimated damage).
The Methodological approach to the damage detection and localization problems is illustrated by a flow chart as shown in Fig. 1.

Numerical model of the beam structure
We consider a simply supported unidirectionally [0 0 ] (AS4/3501-6) graphite epoxy beam.The finite element model is discretized in 10 elements.Each node of the finite element has three degrees of freedom; normal displacement along the axis, a rotation around the axis and the longitudinal displacement .Since the beam is macroscopically considered homogeneous, i.e. there is only one layer, the shear correction coefficient is the same as for isotropic beam [23].The material properties and beam dimensions of simply supported unidirectionally [0 0 ] (AS4/3501-6) graphite epoxy beam [24] discretized into 10 elements are given in Table 1.The FEA model and the beam structure is shown in Fig. 2.

BAT algorithm
BAT algorithm is a recently introduced nature inspired algorithm.It is introduced by Yang and He [25] and inspired by behaviors of bats.The mammals, known as bats, have wings and a very good feature of generating echolocation.Yang and He [26] have used these characteristics of bats for the development of their BA.This algorithm is based on three important rules, which are given as: 1) For sensing distance, bat uses its echolocation capacity.It also uses echolocation to differentiate between food and prey and background barriers even in the darkness.
2) Bats used to fly randomly with some characteristics like a velocity, fixed frequency and loudness to search for a prey.
3) It also features the variations in the loudness from a large to minimum loudness.
In simulations, we use virtual bats naturally.We have to define the rules how their positions and velocities in a -dimensional search space are updated.The new solutions and velocities at time step are given by: where ∈ [0, 1] is a random vector drawn from a uniform distribution.Here is the current global best location.
In the local search part, once a solution is selected among the current best solutions, a new solution for each bat is generated locally using random walk: where ǫ ∈ [−1, 1] is a random number, while = 〈 〉 is the average loudness of all the bats at this time step.
The Algorithm parameters and Pseudo-code of BAT Algorithm (BA) are illustrated in Table 2.

Particle swarm optimization (PSO)
The basic algorithm of the PSO works on a population of possible solutions called swarm, themselves called particles.These particles are randomly placed in the search space represented the damage location and level based on objective function (frequencies and mode shape).At each iteration, the particles move, taking into account their best position (egoistic movement) but also the best position of its vicinity [26].
In fact, the new speed is calculated from the following formula: where: and are the velocities of the particle for iterations and + 1. "Best particle‰ is the best position of the particle, "best vicinity‰ is the best position of its vicinity at iteration , "position particle‰ is the position of the particle at iteration ., , are constants, is randomly generated at each iteration and, in general, = .
We can then determine the next position of the particle with the speed that we just calculated, i.e.: where is the position of the particle at iteration .The principle of the algorithm can be more easily seen through the flowchart given in Fig. 3.

The objective function of the optimization problem
The objective function is based on changes in natural frequencies and modal scale factor (MSF).

Changes in natural frequencies
Natural frequencies or eigenvalues are used as diagnostic parameters in structural assessment procedures using vibration monitoring.One great advantage of using only eigenvalues in the damage assessment of structures is that they are cheaply acquired and the approach can provide an inexpensive structural assessment technique.The objective function to be minimized is defined as follows [27]: where -mode number ( = 1, 2, 3,…, ), -test natural frequencies and -calculated natural frequencies.
The are the natural frequencies, which are applied to our damage detection system as inputs.An objective function equal to zero indicates an exact match between the measured and calculated frequencies.

Modal scale factor (MSF)
The modal scale factor is used as a diagnostic parameter in structural assessment procedures using vibration monitoring [28]: where -calculated mode shape and -test mode shape The objective function based on Modal Scale Factor and Changes in Natural Frequencies to be minimized in this algorithm is given by:

Results and discussion
In this study, instead of experimental measurements, numerically generated measurements, were used for a simply supported beam structure with 10 elements and different single-and multiple-damage scenarios.In the first scenario D1, damage is near the hinged support end, positioned in element 2 with damage rate 30 %.In the second scenario D2, damage is near the roller support in element 9 with damage rate 10 %.While in the third scenario D3, damage is spread over three elements; namely element 2 with damage rate 60 %, element 5 with damage rate 40 % and element 8 with damage rate 20 %.The objective function based on Modal Scale Factor and Changes in Natural Frequencies is minimized in this unconstrained optimization problem using BAT and PSO algorithms.

Single damage scenarios
In this case, we created single damage in element 2 with damage rate 30 % (D1) and element 8 with damage rate 10 % (D2), as shown in Fig. 3, by finite element method and calculated the frequencies and mode shapes of the damaged structure.These parameters are used, as experimental data, in our approach for detecting and locating damage.The results obtained from optimization are shown in Figs. 4 and 5.

Multiple damage scenario
We created three damages in Elements 2, 5 and 8 (damage scenario D3) with different percentage 60 %, 40 % and 20 %, respectively using finite element method and we calculated the frequencies and mode shapes of the damaged structure using BAT and PSO algorithms.These parameters are used, as experimental data, in our approach for detecting three simulated damages.The results obtained from optimization are shown in Fig. 6.
Figs. 3-6 show that BAT algorithm is able to detect the location and the magnitude of damaged elements in all scenarios using complete data within the first few iterations (between 5-20 iterations) with high accuracy.We conclude that BAT algorithm is more efficient than PSO algorithm.

Damage detection with noise
In order to investigate the effects of noise on the results of our approach (Eq.( 10)), the th noise response (noisy), is simulated as [28]:

E l e m e n t 1 E l e m e n t 2 E l e m e n t 3 E l e m e n t 4 E l e m e n t 5 E l e m e n t 6 E l e m e n t 7 E l e m e n t 8 E l e m e n t 9
E l e m e n t 10 where is the noise level and is a random number in the interval [−1, 1].In this study σ is considered as 5 % and 10 %.The results of the application of the algorithms to damage detection and localization for beam taking into account noise are plotted in Figs.7 and 8, and Tables 3 and 4.  with noise level of 5 % and 10 %.
The Optimization results for single element damaged beam with noise level of 5% and 10% using BAT algorithm are shown in Fig. 7 and Table 3.

Multiple damage scenario
Optimization results for multiple damage beam with noise level of 5 % and 10 % using BAT algorithm are shown in Fig. 8 and Table 4.
For multiple damage scenario with noise, the stiffness of elements 2, 5 and 8 were reduced by 60 %, 40 % and 20 %, respectively with noise level of 5 % and 10 %.
From the results, we can observe that, when noise is taken into account in the damage detection problem, our approach based on BAT algorithm can detect the damage locations with good accuracy. Khatir

Conclusions
In this paper, an inverse problem method is proposed in order to detect and estimate damages in beam-like structures using PSO algorithm and new powerful algorithm called BAT algorithm.Several damage scenarios are considered to find the effect of damage severity and number of damaged elements.The BAT algorithm and PSO evaluates the location and severity of damage by minimizing an objective function based on the natural frequencies and mode shapes for undamaged and damaged tested beam with and without noise and considering different damage scenarios.The results showed high accuracy of the method using BAT algorithm, even when multiple-damage is present, within the first few iterations.When noise is taken into account in the damage detection problem, our approach based on BAT algorithm can detect the damage locations correctly with high accuracy.

Fig. 1 .
Fig. 1.Methodological approach to the damage detection and localization using BAT algorithm

1 . 5 .
Damage scenario D1 The results of damage scenario D1 are shown in Figs. 4 and 5. a) D2 using BAT algorithm b) D2 using PSO algorithm Fig. Damage results for single element damaged beam using BAT and PSO algorithms (Element 9, Damage rate 10 % (D2))

Fig. 7 .
Single element damaged beam with different noise level

Fig. 8 .
Multiple damage beam with different noise level 2818.REPUBLISHED PAPER.NUMERICAL STUDY FOR SINGLE AND MULTIPLE DAMAGE DETECTION AND LOCALIZATION IN BEAM-LIKE STRUCTURES USING BAT ALGORITHM.SAMIR KHATIR, IDIR BELAIDI, ROGER SERRA, MAGD ABDEL WAHAB, TAWFIQ KHATIR . REPUBLISHED PAPER.NUMERICAL STUDY FOR SINGLE AND MULTIPLE DAMAGE DETECTION AND LOCALIZATION IN BEAM-LIKE STRUCTURES USING BAT ALGORITHM.SAMIR KHATIR, IDIR BELAIDI, ROGER SERRA, MAGD ABDEL WAHAB, TAWFIQ KHATIR

Table 3 .
Single element damage with noise 5 % and 10 %

.1. Single damage scenario
For single damage scenario with noise, the stiffness of element 8 was reduced by 10 percent Samir, has conducted the research work presented in this paper in the framework of his Ph.D. research degree.Belaidi Idir has directed and supervised this research work, contributed in the conception, design and proposed methods of the work during his Ph.D. studies at LEMI in Boumerdes University, Algeria.Roger Serra has supervised this research work, contributed in the vibration analysis of the work and directed the first author during his research stay at LMR Laboratory INSA Centre Val de Loire, France.Magd Abdel Wahab, has supervised this research work, contributed in the conception and design of the work and directed the first author during his research stay at Ghent University, Belgium.Khatir Tawfiq has supervised this research work, contributed in the dynamic analysis of the work and Matlab programming.