# Direct Index Method of Beam Damage Location Detection Based on Difference Theory of Strain Modal Shapes and the Genetic Algorithms Application

Authors: Bao Zhenming; Zhao Deyou; Lin Zhe; Ma Jun; Zhen Yaobao
Journal: Mathematical Problems in Engineering
 (2012)

---

## Abstract

Structural damage identification is to determine the structure health status and analyze the test results. The three key problems to be solved are as follows: the existence of damage in structure, to detect the damage location, and to confirm the damage degree or damage form. Damage generally changes the structure physical properties (i.e., stiffness, mass, and damping) corresponding with the modal characteristics of the structure (i.e., natural frequencies, modal shapes, and modal damping). The research results show that strain mode can be more sensitive and effective for local damage. The direct index method of damage location detectionISMSD is based on difference theory, without the modal parameter of the original structure. FEM numerical simulation to partial crack with different degree is done. The criteria of damage location detection can be obtained by strain mode difference curve through cubic spline interpolation. Also the genetic algorithm box in Matlab is used. It has been possible to identify the damage to a reasonable level of accuracy.


---

## Body

## 1. Introduction

Recently the development of new methods of evaluating the integrity of structures has been attached great importance. The vibrational characteristics of structures can be easily alternated once the damage occurs. In order to analyze the state of the structure integrity and the location of structural damage, the modal analysis with vibrational test data is widely employed [1], though the modal vectors were effectively ignored earlier [2, 3]. The reason is that the parameters of modal analysis depend only on the mechanical characteristics of the structure. Considering this, the characteristics of modal vibration, that is, the natural frequencies, mode shapes, and so on are considered to represent the state of structure [4].Although damage index based on strain modal shape is much more sensitive [5–8], in most cases, this index can only indicate the damage location both in undamaged and damaged state with modal data. In fact, these kinds of methods are difficult to employ, because the baseline under suddenness and geological disaster is difficult to master. The problem is solved by Gu et al. [9], who put forward the direct index method of damage location detection ISMSD based on difference theory without baseline modal parameters. The index ISMSD contains three elements: the effective distance ratio between two adjacent effective extreme points, the absolute difference value of the two effective extremums, and the maximum absolute value of the effective extremum.Then, the genetic algorithm was gradually applied to solve the problem of damage detection [10]. The method is demonstrated on a simulated beam example and an experimental plate example. Based on the genetic algorithm, Yi and Liu estimated the structural damage with measured dynamic data [11]. Some improved strategies such as multiparent crossover and adjustment of variables are introduced to detect the damage of fixed-end beams and continuous beams. And some satisfied results are achieved. Based on conventional modal analysis theory, Mares and Surace [12] coded a function of binary, and a case of truss-type structure is used to identify the structural damage. And the accurate identification of both the location and the degree of the damage had become possible to be achieved.Genetic algorithm appears to provide a robust search procedure for solving difficult problems. Due to the way the genetic algorithm explores, the region of interest avoids getting stuck at a particular local minimum and locates the global optimum. To formulate an objective function, with the parameters related to the physical properties and state of the structure, is the aim of the detection of damage in structures. When evaluated with the true parameters, the maximum value of the objective function is obtained. Genetic search algorithm is an optimization procedure which can be employed to determine the values of these parameters by following an iteration process, selecting parameters to maximize the objective function. When the optimization procedure is over, the state of the structure, that is, where and how it is damaged is known.The objectives of this paper are to develop and apply models to diagnose the damage location and damage degree. This paper is organized as follows: the first section provides a brief introduction to GA, the next section contains results and analyses including performance evaluation of the methodology, and finally the conclusions are presented.


## 2. Model Development

Modal test of continuum structure uses a series displacements of discrete measuring points to describe all order modes of mode functions. The mode results of all order are the ratio of normalization displacements of discrete measuring points. The strain distributing can be calculated by the difference of deformable displacements.The vibrational differential equation of straight Euler beam in transverse section is(2.1)∂2∂x2(EI∂2u∂x2)+ρA∂2u∂x2=0.The displacement containing all order modes of displacementϕr(x)(r=1,2,…) is
(2.2)u(x,t)=U(x)ejωt=∑r=1∞ϕr(x)Qrejωt.Curvature mode is the difference of displacement mode.The discrete curvature mode can be approximately calculated by center difference method as(2.3)ϕrk′′=ϕr(k-1)-2ϕrk+ϕr(k+1)dx2≈ϕr(k-1)-2ϕrk+ϕr(k+1)h2ϕrk is the rth order of displacement amplitude; k is the calculation point; h is the distance between two adjacent calculation points.Supposez(x) is the distance from the middle bending surface to any point, then the strain of this point in x direction is
(2.4)εz,x=∂u∂x≈∂2ϕrk∂x2z(x)=ϕrk′′z(x).The relationship between strain mode and curvature mode can be transmitted through formula (2.4). Strain mode shape is more sensitive than displacement mode shape at damage detection.The curve of strain mode shape has extremum and sharp variation at the damage location for simply supported beam, which is the common structural style. The difference curve of strain mode shape must be smoothed by cubic spline interpolation. Central difference method is a numerical method. The advantages of this method are symmetry, a certain accuracy, and easy to use. The format at the nonboundary nodes are central difference. The format at boundary nodes is eccentric difference.The beam is separated inton-1 units, namely, n nodes. The equidistance is h. Φrε(x) is the rth order strain mode function. The basic equidistance difference format is
(2.5)(dΦrεdx)i=Φrε(xi+1)-Φrε(xi-1)2h.The difference formulation at boundary node 1 is(2.6)(dΦrεdx)1=Φrε(x2)-Φrε(x1)h.The difference formulation at boundary noden is
(2.7)(dΦrεdx)n=Φrε(xn)-Φrε(xn-1)h.When difference curve fits, the 2nd order difference values of two extreme points are(2.8)(d2Φrεdx2)1=Φrε(x1)+Φrε(x3)-2Φrε(x2)2h2,(d2Φrεdx2)n=Φrε(xn-2)+Φrε(xn)-2Φrε(xn-1)2h2.Nearby the sharp variation section or the peak value of the mode shape, the extreme points at strain mode curve in damage are located. The damage location is at the curve’s sharp variation section.For actual engineering, Extreme point is unnecessarily just zero at the difference of ideal curve. Three reasons are explained as follows: firstly, the space principles of measuring point or unit classification in numerical simulation are not just at the damage location. Secondly, the accuracy problem of both the strain mode test and the computation is to be considered. Thirdly, the accuracy of difference computation is to be considered.Interpolation method is used to solve the question of the zero point of strain mode difference curve and confirm the extremums. The computational value at nodes is equal to the supposed function. Many function values are calculated and the function curves draw. Cubic spline interpolation is used in this paper, a widely used interpolation method.


## 3. Genetic Algorithm

Complex adaptive systems are extremely difficult to comprehend, when unexpected and unpredictable results occur. Holland succeeds to solve this problem by presenting genetic algorithm. In recent years, genetic algorithms are very popular procedure of robust research for solving maximizing or minimizing a given objective function often subject to some constraints [13–20].Genetic algorithm derives from the process of natural selection and evolution. Because of the inherent advantage of being able to process with a large population of designs and facilitating arrival at the globally optimal solution, the philosophy of “survival of the fittest” has been adopted. It is necessary to devise a general coding system for the representation of the design variables, namely, a directly analogy of the DNA structure of chromosomes. Most commonly, the design variables are coded by a bit-string which is a binary representation. It can only be coded as integers. The progress of genetic algorithm is in the same way of the natural evolution of a species: the fundamental concepts of reproduction, chromosomal crossover, and occasional mutation of genes. Fitness function, which determines possible solutions to the problem, is used to estimate the quality of the represented solution (chromosome). Crossover is a reproduction operation in GA, which is used to vary the programming from one generation to the next by exchanging genetic information between parent chromosomes.In applying the evolution theories to designing optimization, a number of candidate design variables either randomly or heuristically are created. Then, they evolve over generations to produce new designs which are “fitter.” The “fitness” of the designs is evaluated according to the objective function, a specific optimization problem. The point with no further improvement is the solution. This paper applies genetic algorithm to the problem of damage detection using strain mode.In the application to the damage detection in structures, the aim is to formulate an objective function in terms of damage degree of the structure. The objective function must be formulated in such a way that the maximum value is obtained when evaluated with the true parameters. Genetic algorithm can be employed to determine the values of these parameters by following an iteration process, selecting parameters to maximize the objective function. When the optimization procedure arrives at the solution, the values of the parameters indicate the state of the structure, that is, if, where and how it is damaged. The detailed steps are shown as follows.Step 1.
Initialize: the initial population is generated with binary coding, and each individual represents a initial solution.Step 2.
Individual evaluation: calculate the fitness value of each individual.Step 3.
Convergence judgment: implement Step4, if the fitness value is still not convergent or cannot reach the limit of iteration times, otherwise, terminate the calculation.Step 4.
Individual crossover, mutation, and selection:  return to Step2 after executing crossover (one-point crossover), mutation (discrete mutation), and selection (sampling randomly) process.Thus, the search process of genetic search can be seen in Figure1.Figure 1
Flowchart of the genetic search.


## 4. Example

In the genetic search procedure, selection, crossover, and mutation operators, each with a given probability, are applied to each current population to create the new generations. The results published by DeJong, Gerefenstette, and Schaffer indicate the most appropriate parameter setting: a population size of 20–30 individuals, a crossover rate of 0.6–0.95, and a mutation rate 0.01-0.02.Take the “damage rate 6.25%” as an example to apply the calculation. Figure2 shows that during the 100 evolutional generations, the result began to converge after the 80th iteration. Therefore, the final result is obtained: in the first order, the best value of ISMSD(2) is convergent to 0.0048.Figure 2
Calculation results.The mathematical model has solved the maximum of formula (4.5), in which x∈[1,41]. Selecting binary code, the swarm number is 10. The length of binary code is 19. Crossover probability is 0.95. The mutation probability is 0.01. After the 100 times of genetic iterative, the results were close to the data of Table 1.Table 1
Damage location indices.
Damage location
Order
Damage degree
6.25%
12.5%
25%
31.25%
Quarter span
First order
ISMSD(1)
0.0018
0.0049
0.0111
0.0159
ISMSD(2)
0.0048
0.0092
0.0197
0.0238
ISMSD(3)
0.0011
0.0013
0.0114
0.0167
ISMSD(4)
0.0016
0.0022
ISMSD(5)
0.0010
Second order
ISMSD(1)
0.0162
0.0205
0.0673
0.1003
ISMSD(2)
0.0202
0.0497
0.1054
0.1647
ISMSD(3)
0.0137
0.0162
0.0655
0.0997
ISMSD(4)
0.0127
0.0000078
0.0035
ISMSD(5)
0.0122
0.0105
ISMSD(6)
0.0154
Third order
ISMSD(1)
0.0787
0.0877
0.1005
0.1335
ISMSD(2)
0.0110
0.0661
0.1699
0.2660
ISMSD(3)
0.0457
0.0754
0.0885
0.1821
ISMSD(4)
0.0544
0.0539
0.0649
0.0644
ISMSD(5)
0.0701
0.0707
0.0753
Note: overstriking index and overstriking data correspond to damage location.The example was a numerically simulated simply beams (see Figure3) with a finite element model. The characteristics of the beam were as follows: length l=0.4 m, cross-sectional area A=0.0002 m2, in quarter span damage, a one fifth length of the beam width, modulus of elasticity E=211Gpa, density ρ=7850 kg/m3, Poisson radio μ=0.33, damage extent α1=6.25%h, α2=12.5%h, α3=25%h, α4=31.25%h, eight node solid 45 unit in Ansys, 25 grids at transverse section, 16 grids at vertical section, damage width w=0.002 m.Figure 3
FEM modal of local damage in simply supported beam.On the basis of numerical simulation results of the strain modes, the smooth strain mode difference curve can be draw based on central difference and cubic spline interpolation.The figures (from Figures4, 5, 6, and 7) were the smooth strain mode difference curves after cubic spline interpolation. The damage extents were 6.25%, 12.5%, 25%, and 31.25% of the height.When damage 6.25%, the first three strain modal difference curves.
(a)
First order
(b)
Second order
(c)
Third orderWhen damage 12.5%, the first three strain modal difference curves.
(a)
First order
(b)
Second order
(c)
Third orderWhen damage 25%, the first three strain modal difference curves.
(a)
First order
(b)
Second order
(c)
Third orderWhen damage 31.25%, the first three strain modal difference curves.
(a)
First order
(b)
Second order
(c)
Third orderWith the exception of supports, strain mode difference curves have sharp variation at the damage location. For different damage degrees, strain mode curve is slightly different, but the rule has some consistency.


### 4.1. Mathematical Model of Direct Index Method of Damage Location

#### 4.1.1. Direct Index Method of Damage Location

From the data change at sharp variation zone between two adjacent extremums, damage index was presented. The extents of the curve’s variation were reflected.


#### 4.1.2. Effective Distance Ratio

The distance between two adjacent effective extremums is effective distance. Effective distance ratio is the ratio of effective distance to span.Suppose: span l, effective extremum points:x0,x1,x2,…,xp, effective extremums: y0,y1,y2,…,yp, then the effective distance ratios are
(4.1)I1j=(xj-xj-1)l,j=(1,2,…,p).The smaller distance ratio of the two adjacent effective extremum is, the greater probability of damage is as follows:(4.2)I2j=|yj2-yj1|,j=(1,2,…,p).Let(4.3)I3j=I2jI1j,j=(1,2,…,p).


### 4.2. The Absolute Maximum of Effective Extremum:ISMSD

As strain mode difference curve of different orders was nonnormalized, the computational results had no comparison.In order to solve this problem, the results of formula (4.3), divided the absolute maximum of effective extremums of all orders, the data of different orders can be compared
(4.4)ISMSD(j)=I2jI1jmax(|y0|,|y1|,|y2|,…,|yq|).Namely,(4.5)ISMSD(j)=l|yj2-yj1|(xj-xj-1)max(|y0|,|y1|,|y2|,…,|yq|).Formula (4.5) is called the direct damage location index. The greater value it is, the more probability of damage occurs. Combined with the rules of difference curve, the damage location can be determined.Based on formula (4.5), the numerical computation results of damage location indices are listed in Table 1.From the illustrative examples used, the method using genetic algorithms is very useful to identify both the damage location and the damage degree with a reasonable accuracy. The striking feature of genetic algorithm is that the optimal solution can be determined even with perturbed data having noise. As a parallel search technique, the genetic algorithm makes this approach attractive for problems of large dimensions.


## 4.1. Mathematical Model of Direct Index Method of Damage Location

### 4.1.1. Direct Index Method of Damage Location

From the data change at sharp variation zone between two adjacent extremums, damage index was presented. The extents of the curve’s variation were reflected.


### 4.1.2. Effective Distance Ratio

The distance between two adjacent effective extremums is effective distance. Effective distance ratio is the ratio of effective distance to span.Suppose: span l, effective extremum points:x0,x1,x2,…,xp, effective extremums: y0,y1,y2,…,yp, then the effective distance ratios are
(4.1)I1j=(xj-xj-1)l,j=(1,2,…,p).The smaller distance ratio of the two adjacent effective extremum is, the greater probability of damage is as follows:(4.2)I2j=|yj2-yj1|,j=(1,2,…,p).Let(4.3)I3j=I2jI1j,j=(1,2,…,p).


## 4.1.1. Direct Index Method of Damage Location

From the data change at sharp variation zone between two adjacent extremums, damage index was presented. The extents of the curve’s variation were reflected.


## 4.1.2. Effective Distance Ratio

The distance between two adjacent effective extremums is effective distance. Effective distance ratio is the ratio of effective distance to span.Suppose: span l, effective extremum points:x0,x1,x2,…,xp, effective extremums: y0,y1,y2,…,yp, then the effective distance ratios are
(4.1)I1j=(xj-xj-1)l,j=(1,2,…,p).The smaller distance ratio of the two adjacent effective extremum is, the greater probability of damage is as follows:(4.2)I2j=|yj2-yj1|,j=(1,2,…,p).Let(4.3)I3j=I2jI1j,j=(1,2,…,p).


## 4.2. The Absolute Maximum of Effective Extremum:ISMSD

As strain mode difference curve of different orders was nonnormalized, the computational results had no comparison.In order to solve this problem, the results of formula (4.3), divided the absolute maximum of effective extremums of all orders, the data of different orders can be compared
(4.4)ISMSD(j)=I2jI1jmax(|y0|,|y1|,|y2|,…,|yq|).Namely,(4.5)ISMSD(j)=l|yj2-yj1|(xj-xj-1)max(|y0|,|y1|,|y2|,…,|yq|).Formula (4.5) is called the direct damage location index. The greater value it is, the more probability of damage occurs. Combined with the rules of difference curve, the damage location can be determined.Based on formula (4.5), the numerical computation results of damage location indices are listed in Table 1.From the illustrative examples used, the method using genetic algorithms is very useful to identify both the damage location and the damage degree with a reasonable accuracy. The striking feature of genetic algorithm is that the optimal solution can be determined even with perturbed data having noise. As a parallel search technique, the genetic algorithm makes this approach attractive for problems of large dimensions.


## 5. Conclusion

The damage indices of the first and second orders are obviously greater than others in the table, and the damage degrees are in accordance with it. The direct index method can accurately detect the damage location and damage degree, especially lower damage. Genetic algorithm is a powerful tool. From the example, accurate identification of both damage location and damage degree had been possible in the structure model. For a protracting beam, there is no macrocrack. The difference curve of strain modes at the damage location has sharp variation, but the numerical value probably is not zero.



---

*Source: 302482-2012-12-11.xml*