Damage localization and quantification of composite beam structures using residual force and optimization

Structural Health Monitoring (SHM) and impact monitoring of composite structures have become important research topics in the recent year. In this research, a non-destructive vibration-based damage detection method is formulated using Genetic Algorithm (GA) and compared with classical method. The robustness and reliability of the capability to locate and to estimate the severity of damage, based on changes in dynamic characteristics of a structure, is investigated. The objective function for the damage identification problem is established by using the residual force method (FRM). Numerical experiments using finite element analysis are performed on composite beams with different damage scenarios in order to clarify the validity of the developed technique. The comparison between estimated and real damage illustrates the efficiency of the algorithm in damage detection. The results show that the present approach is correct and efficient for detecting structural local damages in composite beam structures.


Introduction
The non-destructive testing evaluation and condition assessment of aging infrastructure and mechanical engineering have become important research topics and received much attention for structural damage detection.There is always a need to develop and implement accurate methods of damage detection in structures using intelligent artificial methods.The frequencies and mode shapes are the most popular modal parameters used in the damage detection techniques in recent years [1][2][3][4][5][6].The basic idea of these techniques is that the change in modal parameters can be used to detect the damage in a structure.If damage exists in a structure, an unbalance error resulting from the substitution of a refined finite element model and the measured modal data into the structure eigenvalue equation, which is called the residual modal force, takes place.This residual modal force can be used as an indicator of damage detection.In this paper, to quantitatively identify the extent and location of the damage, we used the residual modal force as an objective function.
Damage detection and localization in thin plates based on vibration analysis using BAT algorithm in beam like and complex structures by Khatir et al. [7].Sheinman [8] proposed a modified residual force vector to detect damage in beam-like structures.The identification of damage was formulated as an optimization problem using three objective functions; a) the change of natural frequencies, b) Modal Assurance Criterion (MAC) and c) Modal Assurance Criterion natural frequency.As introduced in [9], the residual forces were used to detect damage in beam and truss structures based on mass and stiffness matrices.The Coordinate Modal Assurance Criterion (COMAC), which made use of the numerical and experimental vibration modes, was used to determine the magnitude and position of damage by Iturrioz, et al. [10].The application of the change in modal curvatures to detect damage in the concrete bridge was introduced by Wahab and De Roeck [11].The application of ultrasonic guided waves generated by piezoelectric smart transducers has become one of the widely-used techniques in structural health monitoring.In references [12,13], the authors used this concept for damage identification in a shear structure and a tapered and non homogeneous beam using all mode shapes that were considered in the objective function formulation.In reference [4], the authors presented a new approach of inverse damage detection and localization based on model reduction.The problem was formulated as an inverse problem, where an optimization algorithm was used to minimize the cost function expressed as the normalized difference between a frequency vector of the tested structure and its numerical model.The damage indicator based on vibration data are used to detect and locate defects in stratified beam structures [14].A novel wave number analysis approaches were developed and discussed to show how they could be used to investigate Lamb wave interactions with delaminated plies [15].
The residual force method is effective in damage localization using modal data [16].The authors of reference [16] used residual modal force vectors for damage localization.The quantitative of multidamage monitoring method for large-scale complex composite was presented in [17].In reference [18], the authors proposed a new approach to capture and process the full matrix of all transmit-receive time-domain signals from the array.The damage localization approach based on the residual forces method to locate damage.The subspace rotation damage identification algorithm was introduced by Kahl, and Sirkis [19] using residual forces.The residual force vector was used in many methods in the literature or detecting and locating damage in structures [20,21].
In the present paper, in the first section, we develop a damage indicator residual force coupled with FEM.In the second section, we quantify damage using FEM coupled with optimization method, i.e.GA.The residual modal force is used as an objective function to compare the measured data with the ones calculated by GA.

The residual forces method
Testing structural damage detection and localization methods based on the residual force vector [22] is studied in this paper.The damage index of the th element is expressed as the change of the rigidity of a finite element, i.e.: where [ ] and [ ] are the th element of the elementary matrix of the damaged and undamaged structure, respectively.In Eq. (1), ∆[ ] represents the variation of stiffness, is in the interval [0,1] and indicates a loss of rigidity of th element, i.e. = 0 for undamaged and = 1 for damaged element.We consider that the mass matrix of the damaged structure is not affected by damage, and rigidity matrix of the damaged element change as given below: where = 1, … , .The modal residual force vector can be written as: Eq. ( 4) can be written in a matrix form as: where the coefficient of the matrix [ ] is: where is the actually th mode node force vector of the th element in globe coordinates.The modal residual force vector can be written as: Eq. ( 5) can be rewritten as: where is the number of modes, while is the number of elements.The solution of the system of equations in Eq. ( 8) allows us to determine the values of the damage indicators:

Finite element method
We consider in this article, a unidirectional composite beam.We model the structure using SI12 finite elements [23].Each node of a finite element has three degrees of freedom, i.e. displacement normal to the beam, longitudinal displacement and rotation around the -axis, as shown in Fig. 1.Therefore, the total number of degrees of freedom of the finite element SI12 is 12.

Fig. 1.
Unidirectional CFRP [18] The coefficients of stiffness of the beam using the stiffness are given by the following relationships: where : Shear correction factor, : width of beam and is bottom coordinate of the th layer: where : Young's modulus of the th layer in the direction, : transverse shear modulus of the th layer.
The generalized mass densities , and are given by the following equations: where : width of the beam, : coordinate of the th layer.The elementary stiffness matrix [ ] can be written as: With: The elementary mass matrix [ ] is defined by the following relationship: where is the matrix of shape functions and is the generalized matrix densities.The material properties and beam dimensions of CFRP (Carbon Fiber Reinforced Polymer) composite beam [18] are given in Table 1.In this paper, we proposed several damage scenarios as summarized in Table 2, in which damaged element(s) and damage rate are given for each damage scenario.

Evaluation of damage indicators
The natural frequencies of the undamaged and damaged structures for all damage scenarios are shown in Table 3.Values of the damage indicator , for damage scenarios 1 to 5, are shown and plotted in Figs. 2 to 6, respectively.For damage scenarions 1 and 2, it can be seen from Figs. 2 and 3 that the value of the indicator is significant for the damaged element and therefore this method allowed us to detect and localized the damage for a single damage scenario.For the multiple damage scenarios, Figs. 4 to 6, it can be seen that the damage elements are correctly identified using the damage indicator.Therefore, this technique is also suitable for multiple damage cases.

Genetic algorithm (GA)
In this section, we present the optimization method, Genetic Algorithm (GA), which will be used along with the residual forces to quantify damage.The genetic algorithm used in this paper is the Standard genetic algorithm, which is based on Binary Encoding.The residual forces are used in the objective function to compare the data between measured and calculated by GA.Genetic Algorithm is a global probabilistic search algorithm inspired by Darwin's survival-ofthe fittest theory originally developed by Holland [25][26][27] combines algorithms to solve optimization problems using the principles of evolution.In this optimization method, information about a problem posed, such as variable parameters, is coded into a genetic string known as an individual (chromosome).In this study, we have two chromosomes that represent damage location and severity.Each of these individuals has an associated fitness value, which is usually determined by the objective function using data of residual forces to be minimized.This optimization method has been shown to be able to solve the optimization problem through mutation, crossover and selection operation applied to individuals in the population.All numerical studies presented in this paper are implementing in MATLAB.In applying GA, the following steps are programmed: 1) Initializing population of individuals, created in real encoding as a random generation.Each individual has two chromosomes corresponding to the location and level of damage.
2) Evaluating each individual by introducing the proposed parameters (damage element and level) into FEM using genetic algorithm.
3) Evaluating populations according to their fitness value and then ranking them.A proportion for breeding a new generation will be selected.
4) Selecting the best populations.Flowchart of a basic genetic algorithm is given in Fig. 1.Through series of identification tests, the following genetic parameters were chosen based on the accuracy of results: The number of generation used between (50-300) depend number of damaged element, population size = 300, crossover rate = 0.8 i.e. 800 individuals were selected for crossover, the mutation rate = 0.01.The mutation was used to avoid the convergence of the solution toward local optimums by creating diversity.

Results and discussion
In case of using GA, the residual force method (FRM) is used as an objective function combined with finite element models.The same damage scenarios as presented in Table 2 are used herein.The results are shown in Figs. 8 to 12, in which three plotted for each damage scenario are presented; a) Fitness b) Damage indicator c) Convergence of damaged elements.For the single damage scenario, D1 and D2, it can be seen, in Figs. 8 and 9, that the algorithm is converged after few iterations, i.e. before reaching 10 iterations, and the damage location and severity are correctly identified with higher precision.Similarly, for multiple damage scenarios, scenarios D3 to D5, in Figs. 10 to12, all damage locations and severities are correctly found.However, as the number of damage elements increases, the number of iterations increases and the convergence becomes slower.In general, we conclude that the results show that GA algorithm is able to quantify damage in most scenarios using data of residual force method within the first few iterations with high accuracy.However, for damage scenario D5, more than 130 iterations are required in order to reach convergence, as it can be seen in Fig. 12.

Conclusion
In this article, a method for inverse problem is proposed for damage detection and severity in CFRP composite beam structures.The proposed technique is based on Genetic Algorithm (GA), residual force method (RFM) and finite element method (FEM).The objective function is based on using the calculated data of residual forces and compare them to the measured ones using GA.Numerical experiments using FEM are performed on composite beams with different damage scenarios in order to clarify the validity of the developed technique.The comparison between estimated and real damage illustrates the efficiency of the algorithm in damage detection.The results show that the present approach is correct and efficient for detecting structural local damages in case of single and multiple damage scenarios.

5 )
Crossover of individuals to produce the population of the next generation.6) Mutation of a specified level of the resulting population.7) Replacing of the old population by new one and coming back to the step 2.

10 . 11 .
a) Fitness convergence = [0 0 0.25 0 0 0 0 0.25 0 0 ].b) Convergence of damaged elements c) Damage identification Fig. Convergence and damage identification for damage scenario 3 a) Fitness convergence = [0 0.10 0 0.20 0 0.30 0 0.40 0 0 ].b) Convergence of damaged elements c) Damage identification Fig. Convergence and damage identification for damage scenario 4 a) Fitness convergence = [0.050.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 ]. b) Convergence of damaged elements c) Damage identification Fig. 12. Convergence and damage identification for damage scenario 5 A. Behtani has carried out the research work related to the residual forces and has written large part of the manuscript.A. Bouazzouni has supervised the first author for the research on residual forces method.S. Khatir has carried out the finite element analysis and the genetic algorithm optimization, and has written part of the manuscript.S. Tiachacht has contributed in the implementation of the genetic algorithm method.Y.-L.Zhou has supervised S. Khatir for the implementation of the finite element analysis.M. Abdel Wahab has supervised all the work, contributed in the research concept and revised the manuscript.

Amar
Behtani is an Assistant Professor at Mouloud Mammeri University of Tizi-Ouzou.Algeria.He is a research fellow in the Department of Mechanical Engineering Laboratory of Mechanics, Structure and Energetics (LMSE).His research interests include damage detection, steel structure, aeronautical and space engineering, delamination, composites structure and vibration analysis.Amar Bouazzouni is a Professor of structure dynamics at Mouloud Mammeri University of Tizi-Ouzou.Algeria.He is a research fellow in Department of Mechanical Engineering Laboratory of Mechanics, Structure and Energetics (LMSE).His research based on damage identification, dynamics structures, modal analysis, analysis of mechanical structures, and identification of damage of mechanical structures.Khatir Samir is a Ph.D. Student at Ghent University.Received his Master's degree in 2011 and a first Ph.D. in 2016 from University of Boumerdes in Algeria.His research interests include damage detection, steel structure, aeronautical and space engineering, composites structure, vibration analysis, model reduced, 3D complex structure optimization methods and vibration experimental.He has published more than 20 scientific papers in solid mechanics and dynamics of structures.Samir Tiachacht is an Assistant Professor at Mouloud Mammeri University of Tizi-Ouzou.Algeria.He is a research fellow in Department of Mechanical Engineering Laboratory of Mechanics, Structure and Energetics (LMSE).His research interests include damage detection, steel structure, complex structure, composites structure and vibration analysis.Yun-Lai Zhou is a research fellow in Department of Civil and Environmental Engineering at National University of Singapore, Singapore.He received his B.Sc., 2010, in theoretical and applied mechanics, from Northwestern Polytecnic University, China; M.Sc, 2011 and Ph.D. in 2015, both in civil engineering from Technical University of Madrid, Spain.Dr. Zhou has published more than 10 scientific papers in dynamic mechanics, and his research interests include structural health monitoring, system identification, and fracture mechanics.Magd Abdel Wahab is a full time Professor of applied mechanics in the Faculty of Engineering and Architecture at Ghent University, Belgium and an adjunct Professor of computational mechanics at Ton DucThang University, Vietnam.He received his B.Sc., 1988, in civil engineering and his M.Sc., 1991, in structural mechanics, both from Cairo University.Prof. Wahab completed his Ph.D. in fracture mechanics in 1995 at KU Leuven, Belgium.He was awarded the degree of Doctor of Science from the University of Surrey in 2008.He has published more than 290 scientific papers and technical reports in solid mechanics and dynamics of structures.His research interests include fracture mechanics, damage mechanics, fatigue of materials, durability, and dynamics and vibration of structures.

Table 1 .
Material properties and beam dimensions

Table 3 .
Natural frequencies (Hz) of the undamaged and damaged beam structures Mode 1 Mode 2 Mode 3 Mode 4 Mode 5