Differential Evolution: An Inverse Approach for Crack Detection

This paper presents a damage detection technique combining analytical and experimental investigations on a cantilever aluminium alloy beam with a transverse surface crack. Firstly, the first three natural frequencies were determined using analytical methods based on strain energy release rate. Secondly, an experimental method was adopted to validate the theoretical findings. The damage location and severity assessment is the third stage and is formulated as a constrained optimisation problem and solved using the proposed differential evolution (DE) algorithm based on the measured and calculated first three natural frequencies as inputs. Numerical simulation studies indicate that the proposed method is robust and can be used effectively in structural health monitoring (SHM) applications.


Introduction
A crack is a potential source of catastrophic failure in structures.Extensive investigations by researchers have been done to develop structural integrity monitoring techniques.Vibration measurement and analysis being an effective and convenient way to detect cracks in structures is mostly being used for development of various such techniques.
Several nondestructive techniques (NDT) are available for local damage detection [1] using experimental methods like radiography, the magnetic field method, the acoustic method, and so forth.However, for health monitoring of critical and complex structures, in such experimental methods which require prior knowledge of the damage vicinity, accurate predictions may not be suitable.This has led to the development of quantitative global damage detection methods which are based on modal analysis [2,3].Researchers [4] also argue that in view of prohibitive costs and efforts involved in predicting damage to a high level accuracy a better idea is to roughly locate damage in the structure and then use standard NDT methods for closer analysis of the damaged area.
Recently a lot of work has been done using modal analysis to detect, locate, and predict crack severity to a greater accuracy level.Dimarogonas [2] in 1996 provided a comprehensive review of vibration based mode shape analysis followed by some recent reviews [4][5][6][7].The damage detection problem can be defined as a nonlinear inverse problem [3].In conventional model-based detection methods, the minimization of an objective function is defined in terms of the differences between the vibration data obtained by modal testing and those computed from the analytical model.These conventional optimization methods are gradient based and usually lead to a local minimum only.For a more accurate and reliable solution, global optimization techniques are developed such as genetic algorithms (GAs) [8,9], artificial neural networks (ANNs) [10,11], fuzzy logic [12,13], particle swarm optimization (PSO) [14][15][16][17][18][19], and some hybrid methods like genetic fuzzy systems [20].All these approaches have their strengths and weaknesses.However, there is a growing preference for the soft computing approaches based on dfferential evolution (DE) [21,22] due to faster and accurate predictions of global optimization to the objective function.A differential quadrature method combined with an evolutionary optimization algorithm has been proposed for crack detection in cylindrical shell structures.The circumferential crack, which is assumed to be open, is modeled by the extended rotational spring.For acceptability of this method an experimental analysis is also carried out [23].
A comparison has been made between the coupled local minimizers (CLM) method and the differential evolution (DE) algorithm to perform FE model updating for the damage detection in a cracked beam.CLM method is a gradient-based method with multiple local optimization runs whereas the DE algorithm is a direct search approach which uses a population of solution vectors collecting the design parameters [24].
In the current work, a systematic procedure has been developed to calculate the natural frequencies and mode shapes of a cracked cantilever beam with a transverse crack.The process endeavours to study the influence of a crack on natural frequencies and mode shapes.For different relative crack depths and locations, natural frequencies and the mode shapes of the cracked cantilever beam specimen have been found out using a theoretical method.
The rest of this paper is organized as follows.Section 2 introduces the theoretical analysis of a cracked beam structure.The experimental results are presented in Section 3, followed by Section 4 which includes the proposed DE algorithm.The results and discussions as well as conclusions are discussed in Sections 5 and 6, respectively.

Computation of Flexibility Matrix of a Damaged Beam
Subjected to Complex Loading.A beam with cracks has smaller stiffness than that of a normal beam.This decreased local stiffness can be formulated as a matrix.The dimension of the matrix would depend on the degrees of freedom in the problem.Figure 1 shows a cantilever beam of width  and height , having a transverse surface crack of depth  1 .The beam experiences combined longitudinal and transverse motion due to the axial force  1 and bending moment  2 .
Here we consider two degrees of freedom, leading to a 2 × 2 local stiffness matrix.
The relationship between strain energy release rate () and stress intensity factors (  ,  = 1 to 2) at the crack section is given by Tada et al. [25] as where   = /(1 − ] 2 ), for plane strain condition,   = , for plane stress condition,  1 = stress intensity factor for opening mode I due to load  1 ,  2 = stress intensity factor for opening mode I due to load  2 .
From earlier studies [25], the values of stress intensity factors are where the experimentally determined functions  1 and  2 are expressed as follows: The strain energy release rate (also called strain energy density function) at the crack location is defined as () =   /( × ), where ( × ) is the newly created surface area of the crack.
So the strain energy release (  ) due to the crack of depth  1 is calculated as Then from Castigliano's theorem, the additional displacement along the force   is From ( 1) and ( 3), thus we have Advances in Acoustics and Vibration 3 The flexibility influence coefficient   will be by definition Substituting ( 1) in ( 6), we have Using  = (/),  = /, the following expressions are found out  = , and, when  = 0,  = 0;  =  1 ,  =  1 / =  1 .
From the above condition, (9) can be written as where C  = flexibility influence coefficient in  direction (-direction or -direction) due to the load in  direction ( 1 or  2 ).
Calculating  11 ,  12 (=  21 ), and  22 is as follows: The local stiffness matrix can be obtained using the inverse of compliance matrix Converting the influence coefficient into dimensionless form we get

Governing Equations for Vibration Mode of the Cracked
Beam.The cantilever beam as mentioned in Section 2.1 is being considered for free vibration analysis.A cantilever beam of length ", " width ", " and depth "" with a crack of depth " 1 " at a distance "  " from the fixed end is considered as shown in Figure 1.Considering  1 (, ) and  2 (, ) as the amplitudes of longitudinal vibration for the sections before and after the crack position and  1 (, ) and  2 (, ) as the amplitudes of bending vibration for the same section as shown in Figure 2, the free vibration of an Euler-Bernoulli beam of a constant rectangular cross-section is given by the following differential equations: longitudinal vibration lateral vibration The normal functions for the cracked beam in nondimensional form for both the longitudinal and bending vibrations in steady state can be defined as where ( = 1, 12) are the constants to be determined, from boundary conditions.The boundary conditions of the cantilever beam considered in the present study are At the cracked section Also at the cracked section, we have Multiplying both sides of the above equation by / 11  12 we get Similarly, Multiplying both sides of the aforementaioned equation by / 2  22  21 we get where The normal functions ((16a)-(16d)) along with the boundary conditions as mentioned above ((18a) to ( 23)) yield the characteristic equation of the system as where  is a 12 × 12 matrix as given below, whose determinant is a function of natural circular frequency (), the relative location of the crack (), and the local stiffness matrix (), which in turn is a function of the relative crack depth where

Experimental Results
To verify the integrity of the proposed crack detection method and to find out the errors associated with the modeling and measurements, several experiments also have been conducted in the laboratory.Figure 3 shows a schematic diagram of the experimental setup and its description.A cracked cantilever beam has been rigidly clamped to the concrete foundation base.The geometry and material properties of the beam are presented in Table 1.The free end of the beam is excited with a vibration exciter.The vibration exciter is excited by the signal from the function generator.The signal is amplified by a power amplifier before being fed to the vibration exciter.The amplitude of vibration of the uncracked and cracked cantilever beam is taken by the accelerometer and is fed to the vibration indicator for vibration analysis.
The vibration signatures are analyzed graphically by vibration indicators.
Several tests are conducted using the experimental setup on aluminum beam specimens (800 × 50 × 6 mm) with a transverse crack for determining the natural frequencies and mode shapes for different crack locations (i.e., 200 mm, 400 mm, and 600 mm from the clamped end) and crack  depths varying from 1 mm to 5 mm by a step of 1 mm.The cracks were prepared by fine saw cuts perpendicular to the longitudinal axis.This ensures that the crack remains open during the vibrations.At each step the first three bending natural frequencies of the cracked beams were measured.Table 1 gives the corresponding bending natural frequencies of the intact and cracked beams.These specimens are set to vibrate under the 1st, 2nd, and 3rd modes of vibrations and the corresponding amplitudes are recorded in the vibration indicator.Experimental results of frequencies of transverse vibration at various locations along the length of the beam are recorded by positioning the vibration pickup and tuning the vibration generator at the corresponding resonant frequencies.DE was employed to detect cracks utilizing the results from the experimental study.A comparison was made between the experimentally measured natural frequencies of the damaged beam and the ones obtained through the cracked beam model using the objective function as discussed in Section 4.4.
The theoretical results are better than the experimental ones, because of measurement errors.The proposed methods Advances in Acoustics and Vibration have been applied to nine damage cases obtained by combining three different crack positions and three different crack depths.
Then theoretical and experimental results are compared by using the differential evolution method.

Differential Evolution
The differential evolution (DE) algorithm is a population based evolutionary algorithm developed by nondifferentiable continuous space functions [26].Like the genetic algorithm, the optimization process in differential evolution is subject to three basic operations: mutation, crossover, and selection.DE, at the start, randomly initializes a population of size NP solution vectors, which is evolved over  generations to reach an optimal solution.The population size NP is kept constant during the minimization process.Each individual   of the population is a vector that contains a set of  optimization parameters as the problem decision variables and can be expressed as In the crack detection problem the search space is bounded by the parameters such as crack location and depth in the cantilever beam.The minimum and maximum parameter bounds are  min = { 1 min , . . .,   min } and  max = { 1 max , . . .,   max }, respectively.The initial value of the th parameter in the th individual at the iteration  = 0 can be expressed as where rand represents a uniformly distributed random variable within the range [0, 1].

Mutation Operation.
For each target vector  ,,  = 1, 2, . . ., NP, a mutant vector is generated according to the following mutation strategy that is used in the present work: Where  1 ,  2 ∈ {1, 2, . . ., } are mutually different integer numbers and are also different from the index .
The indices   1 ,   2 ,   3 ,   4 , and   5 are mutually exclusive integers randomly generated within the range [1,NP], which are also different from the index .These indices are randomly generated once for each mutant vector.The scaling factor  is a positive control parameter for scaling the difference vector.

Crossover Operation.
The crossover operation is introduced in the DE algorithm, in order to increase the diversity of the vectors.The crossover operation is carried out by randomly exchanging between the original vectors of the population  , and those of the mutant population  ,+1 to obtain the trial vectors  ,+1 = ( where  rand is a randomly chosen integer in the range [1, ].

Selection.
In order to decide if a vector   may be the element of the new population of generation +1, each vector  ,+1 is compared to the previous corresponding target vector  , .If vector  ,+1 yields a smaller objective function value than  , , then  ,+1 is set to  ,+1 ; otherwise, the old value  , is retained.The selection operation can be expressed as with  = 1, 2, . . ., NP.
The aforementioned steps are repeated generation after generation until some specific termination criteria are satisfied.The algorithmic description of DE is summarized in Table 1.

Objective Function Based on Vibration
Data.Damage in a structure makes changes in vibrational parameters such as natural frequencies and mode shapes.In this current study natural frequencies are taken as the damage indicator as they are easier to measure than mode shapes.The objective function chosen for damage estimation is a minimization optimization problem.The objective function based on natural frequencies can be expressed as where    and    are the measured natural frequencies and natural frequencies obtained from a theoretical spring model based on the assumed crack location and depth, respectively.The  is the number of natural frequencies used to evaluate the objective function and   is a weighting factor, whose values are considered to be .In the present work, the first three natural frequencies are considered to determine the objective function.
Stopping Rule.The iterative procedure is terminated when a predefined number of generations or a computational error of 10 −4 is reached.

Pseudo code for DE
(1) Generate a population of solution vectors.
(2) Evaluate the best member of the population  best, .

Results and Discussions
The theoretical analysis and proposed DE algorithm were implemented using MATLAB 7.0 [27].Based on the results obtained from the numerical, experimental, and proposed DE model, the following discussions can be made.A crack in a beam structure causes a change in the stiffness of the beam, which is in turn a function of crack location and  crack depth.As discussed in Section 2, the dimensionless compliance matrix is used to evaluate the stiffness matrix.The dimensionless compliances ( 11 ,  12 =  21 ,  22 ) increasing with the increase in relative crack depth are as shown in Figure 4.
For crack location (e.g., 100 mm) and relative crack depths ( 1 =  1 /) (e.g., 0.1 ), the first three mode shapes are presented graphically in Figure 5, where it is observed that there are reasonable changes in mode shapes due to the presence of a crack in the beam.
The numerical results for the relative amplitude of transverse vibration at different locations of aluminium alloy 2014- 4 cracked specimens for the first three modes are obtained using the theoretical model as per (16a) to (16d) with the help of computer programming.These results for cracked and uncracked specimens are presented in Figure 5 for comparison.It can be seen that the deviations are more prominent with higher modes of vibrations.Table 1 presents the first three natural frequencies obtained from the theoretical model, which is subsequently used in the DE algorithm for crack identification.Figure 6 shows the convergence trend of the fitness value of the objective function with the number of iterations.Table 2 shows the first three natural frequencies of the cracked cantilever beam at different crack locations and depths.The theoretical natural frequencies of the test specimens are presented in Table 3.It is found that the estimation error is reduced by increasing the crack depths for the same locations (Table 5).Also by increasing the distance of the crack from the fixed end of the beam, the estimation error has been found to be decreasing for similar crack depth.Tables 4  and 5 represent the validation of the proposed DE model with the results obtained through the theoretical and experimental models, respectively.

Conclusions
Based on the results and discussions thereof, the following conclusions can be made.The mode shapes and bending frequencies of the cracked elastic structures are strongly influenced by the crack location and its intensity.Though significant changes in mode shapes are observed in the vicinity of crack location, these deviations in mode shapes cannot be used as a measuring tool in the prediction of crack location and its intensity.In the present work, the proposed DE algorithm is found to be an efficient method for damage quantification in terms of crack location and crack depth by minimizing the error between measured and predicted frequencies.It was also observed that the error associated with the prediction is less in the theoretical model as compared to the experimental model.The present work can be implemented for damage assessment in different structures.

Figure 3 :
Figure 3: Schematic block diagram of experimental setup for the cantilever beam.
Figures 6(a)-6(c) show the convergence trend for the theoretical model and Figure 6(d) is for the experimental study.As it can be seen from Figures 6(a)-6(c), the of iterations for convergence increases with the increase in crack location for the same crack depth.Figure 6(d) shows that the fitness value does not converge for the predefined number of iterations.

5. 2 .
Simulation Results.Different crack conditions have been taken to evaluate the performance of DE.Simulation results of 6 test points have been presented in this paper.The error is calculated using the following formula: %error =         predicted value − actual value actual value         × 100.(34)

Table 2 :
Measured natural frequencies for single crack cantilever beam.

Table 3 :
Calculated frequencies using theoretical analysis.

Table 4 :
Crack detection results of the theoretical study by applying DE.

Table 5 :
Crack detection results of the experimental study by applying DE.

Table 4
represents the validation of the proposed DE model with the results obtained as per the theoretical model.The percentage deviations in the results obtained by the experimental study and DE model are presented in Table5.