An Optimized Damage Identification Method of Beam Using Wavelet and Neural Network

An optimized damage identification method of beam combined wavelet with neural network is presented in an attempt to improve the calculation iterative speed and accuracy damage identification. The mathematical model is developed to identify the structure damage based on the theory of finite elements and rotation modal parameters. The model is integrated with BP neural network optimization approach which utilizes the Genetic algorithm optimization method. The structural rotation modal parameters are performed with the continuous wavelet transform through the Mexico hat wavelet. The location of structure damage is identified by the maximum of wavelet coefficients. Then, the multi-scale wavelet coefficients modulus maxima are used as the inputs of the BP neural network, and through training and updating the optimal weight and threshold value to obtain the ideal output which is used to describe the degree of structural damage. The obtained results demonstrate the effectiveness of the proposed approach in simultaneously improving the structural damage identification precision including the damage locating and severity.


Introduction
In recent years, extensive engineering application of wavelet to structure damage identification has been reported, as it has been recognized that using such an approach yields improved damage detection effects [1]. Structure can be regarded as a mechanical system consisting of stiffness, mass and damping matrix.
Once the structure is damaged, structural parameters are also changed, which leads to the change of system frequency response function and modal parame-How to cite this paper: Miao, B.R., Wang, M.Y., Yang, S.W., Luo, Y.X. and Yang, C.J. (2020) An Optimized Damage Identification Method of Beam Using Wavelet and state is established in structural damage identification, which is widely used in fault diagnosis. The genetic algorithm is a kind of optimization method based on biological genetic and natural selection mechanism. Its advantages lie in that it does not depend on the specific domain of the problem, and does not require the continuity of the objective function and the type of the gradient information to the problem. For these reasons, NN optimization has found its ways in railway structural healthy monitoring applications.
The literature on structural damage identification used wavelet is quite extensive and the research activities in this field have focused on a wide variety of applications. As a new signal processing method, wavelet analysis has the characteristic of spatial localization. The wavelet transform of the signal at a certain point is determined by the local information in the low scale, which is determined by the local information, so it is able to analyze the singularity of the signal. Wavelet transform overcomes the problem that some characteristic signals are easily lost in the process of Fourier transform from time domain to frequency domain. In addition, it has more complete analysis functions in both the time domain and the frequency domain. It can not only perform signal vibration measurement and multi-scale analysis, but also perform signal vibration measurement and multi-scale analysis. By analyzing the time-varying properties of the system, the damage is identified, and the existence of the damage is shown directly in the wavelet transform. Fan et al. [1] proposed the wavelet method which is compared with the 2-D strain energy method (SEM) to demonstrate the proposed algorithm is superior with limited sensor data. Gokdag et al. [2] used a wavelet method based on the assumption that a damaged mode shape of a beam is approximately composed of an undamaged mode and other contributors such as measurement induced variations. Xu et al. [3] found two-dimensional directional Gaussian wavelet can characterize directional information about damage and described the operating deflection shapes can be used to locate the real-time dynamic characteristics of a plate. Le et al. [4] proposed a new method based on the similar formulation between the free decay response and the correlation function of the response of a linear mechanical system under ambient excitation.
The neural network method to identify structural damage has also made a lot of useful results. Mehrjoo et al. [5] applied a method for estimating the damage intensities of joints for truss bridge structures using a back-propagation based  [6] combined the back-propagation neural network with the computational mechanics to the detection and sizing of cracks. Yam et al. [7] presented an integrated method for damage detection of composite structures using their vibration responses, wavelet transform and artificial neural networks (ANN) identification. Lu et al. [8] used an innovative neural network-based pattern analysis tool to identify damage-sensitive frequency ranges autonomously and provide detailed information such as the damage type and severity. Adeli et al. [9] applied a new dynamic time-delay fuzzy wavelet neural network model for nonparametric identification of structures using the nonlinear autoregressive moving average with exogenous inputs approach. Liu et al. [10] performed a structure damage diagnosis method combining the wavelet packet decomposition, multi-sensor feature fusion theory and neural network pattern classification. Gao et al. [11] improved wavelet threshold algorithm to eliminate the noise for vibration signals. Sahoo et al. [12] used a hybrid neuro-genetic algorithm to automate the design of neural networks for different types of structures based on the fact that the damage has an important effect on the static as well as dynamic behavior of the structure. Vakil-Baghmisheh et al. [13] proposed a fault diagnosis method based on genetic algorithms (GAs) and a model of damaged (cracked) structure. Buezas et al. [14] presented methodology based on genetic algorithms for the crack detection on damaged structural elements. Meruane, et al. [15] proposed a hybrid real-coded Genetic Algorithm with damage penalization is implemented to locate and quantify structural damage. Hao, et al. [16] applied a genetic algorithm with real number encoding to identify the structural damage by minimizing the objective function, which directly compares the changes in the measurements before and after damage.
However, fewer studies have considered the application of wavelet method and NN optimization to structure damage identification location and damage severity quantitative problems. At present, most of the identification methods are based on natural frequency and modal parameters to identify structural damage, but there are still many shortcomings in the identification efficiency, noise resistance and quantitative analysis. The wavelet analysis method can effectively determine the location of damage using wavelet singularity detection theory. But it is still needed to improve the accuracy of calculation for quantitative determination of the degree of damage. The nonlinear mapping ability of a neural network can establish the mapping relationship between the input and output parameters in the structural damage detection. In this paper, an optimized damage identification method combined the wavelet transform with neural network optimization is presented to improve the identification effects of the damage location and degree.

Wavelet Singularity Identification
If the mother wavelet Mother wavelet satisfying the following admissibility condition: where ( ) ψ ω is the Fourier transform of ( ) where a and b are the scale and translation factor which reflect the frequency and time information of the signal.
( ) . When using the wavelet transform to detect the singularity of the function, the wavelet transform given by Equation (2) can be rewritten as a convolution product: The continuous wavelet transform of the function ( ) f t is: Equation (

BP Neural Network Technology
The artificial neural network has a strong ability of pattern classification and recognition, and has been widely used in structural damage detection. BP neural network is a multilayer feed-forward network with error back propagation. Its structure is simple and the algorithm is easy to implement, the system can solve learning problems in a multilayer network containing layer unit is connected with the right. The basic principle of the BP learning algorithm is the steepest descent method, the central idea is to adjust the weights to minimize the total error of the network [18]. The gradient search technique is used to minimize the error of the actual output and the expected output value of the network.
Taking the 3 layer BP network as an example, the network structure is shown in Figure 1.
 is the hidden layer node. The activation function is the Sigmoid function, therefore the output of the hidden layer and the output layer are respectively [19]:  θ are the threshold of the network; l s , j y the output of the hidden layer and the output layer, respectively.
The total error of n samples can be described as follow: where q E is the error of a sample, l t is described for the network expected output, q is the number of samples.
The output error of the network is a function of the weight of each layer, and the adjustment of the weight is proportional to the negative gradient of the error, so the error E can be changed by adjusting the weight. The gradient descent method is used to train the network, and the weights and thresholds are adjusted continuously until the error or the training times can meet the requirements.

Optimization of Neural Network Based on Genetic Algorithm
Genetic Algorithms (GAs) is a random search algorithm based on natural selection and natural genetic mechanisms. Based on the principle of "survival of the fittest and survival of the fittest", the algorithm is a global optimization algorithm based on genetic manipulation such as selection, crossover and mutation.
Genetic algorithm is used to estimate the fitness of each individual in the population by minimizing the error between the numerical simulation and the measurement frequency to realize the regional large-scale search.
Because the genetic algorithm has a fast convergence speed, high efficiency, global search and suitable for parallel processing advantages, for the optimization of BP neural network to overcome the convergence speed of BP neural network is slow, is easy to be trapped into local minima. Through the BP neural network weights and threshold encode chromosome formation, then to simulate the evolution process of nature, reproduction, crossover and mutation of a chromosome, chromosome evolution, resulting in an optimal solution to the problems on behalf of the chromosome, so as to obtain the optimization of the network weights and threshold. And then the optimized network weights and thresholds are taken as the initial values of the BP neural network weight coefficients and thresholds. The specific operation of the hybrid algorithm to optimize the network weights using a genetic algorithm is as follow [20]: (1) Generates chromosomes and initial populations: The initial weight of the neural network is i W , which it is equivalent to a chromosome { } The error squared sum E and the fitness function F are functions of the weights and thresholds of the network, which is Therefore, the objective function of genetic algorithm optimization is to find a matrix of weights 0 W and threshold matrix 0 θ , so that: (3) Through the selection, crossover, mutation genetic algorithm obtain the new gene values after the mutation. After several iterations, a set of weights and thresholds that minimize neural network errors are obtained, hybrid neural network genetic algorithm flow chart is shown in Figure 2.

Damage Identification Based on Wavelet Transform of Rotation Mode
Mathematically, the general function is smooth or without singularity. If the function is not continuous in a certain position or in a certain order, if the function is called, the singular features or irregularities contained in the basic information of the signal must be considered. The singularity of the general signal is divided into two kinds of situations: 1) at a certain moment, the amplitude of the signal is abrupt, which causes the discontinuity of the signal. This kind of mutation is called the discontinuity point between the first class; 2) the appearance of the signal is very smooth, there is no abrupt change in amplitude, but the first order differential of the signal is abrupt and the first order differential is discontinuous, which is the second kind of discontinuity [21].
The existence of damage causes the local stiffness of the structure to change, which leads to the corresponding change of the natural frequency and mode shape of the structure. Under the general conditions, the structure damage location  Vertical displacement: Rotation mode: Bending moment: Shear force: It can be known that only the vertical displacement and its first derivative (Rotation mode) are continuous, but The singularity of the two order derivative is different, which can be used to identify the damage of the structure by using the rotation model.
In the detection of structural singularity, the damage of the structure can cause the micro mutation of the signal, which is called the singularity point. For the mother wavelet function, the continuous wavelet transform is performed on the corner of the damaged beam. The modulus maxima can occur in the wavelet coefficients at the change of stiffness, which can determine the location of the damaged area.

The Establishment of BP Neural Network
The greater is the fitness function value, the smaller is the error between the network predictive value and the actual value. Table 1 lists the basic parameters of the genetic algorithm. Neural network method has a strong ability of nonlinear mapping and pattern recognition, as long as there is some kind of rule between the data, using neural network can describe the law with good accuracy. The neural network establishes the mapping relationship between input and output, and the hidden layer really performs this kind of relationship. Hidden layer neurons number selection on the convergence speed is very important, increase the implicit layer number can further reduce the training error and improve the precision of the simulation, but will also make the complex network, thereby increasing the training time of the network. The study showed that the three layers BP neural network can simulate the nonlinear relationship between any input and output. The complexity depends on the number of neurons in the hidden layer, and the number of neurons in each layer varies according to the specific problem. In this paper, the three layer BP network is used to input the parameters of the model, the maximum value of the wavelet coefficient modulus of the front n mode of the damaged beam, and the output parameters are the degree of damage. So the network structure is established for 4-8-1 (single point damage), as shown in Figure 3. The output layer neuron uses the pure linear

Numerical Example
In order to verify the validity and accuracy of the algorithm, the three-span continuous beam is used as an example to simulate the structural damage by reduc-

Single Damage Event
The assumption is that there is damage in the span of the continuous beam (i.e., element 50). Figure 5 and Figure 6 respectively show the continuous beam damage before and after the basic mode of vibration mode and the first order rotation mode.
The fundamental mode of the continuous beam is wavelet transform, as shown in Figure 7. It can be seen that there is no obvious singularity characteristic of the basic modal parameters after the damage. Due to the near the bearing is boundary condition effects, the presence of vertical constraint and the deflection approximate to zero and around the corner there is a marked difference.
Therefore, the use of basic mode shape wavelet transform can not accurately   scales are used as the input parameters of the neural network, the output parameter is the damage degree of the element, and the network structure is 4-8-1.
Then, the wavelet transforms maxima of the 7 kinds of damage cases (damage degree) of the element 50 in the first 4 scales are used as the network test data.
As shown in Table 2, the corresponding network test results and relative error can be obtained, and then the mean value of the network prediction error of the damage degree of the continuous beam (element 50) is calculated.

Multiple Damage Event
In order to simulate the multi loss working conditions, the damage is assumed to exist in the element 20 and element 50 in the continuous beam. Figure

Results Analysis and Discussion
The numerical calculation example shows that the structural damage identification method based on the rotation mode is not only more sensitive to damage, but also avoids the singularity caused by the effect of the boundary condition of

Conclusions
In this paper, combining continuous wavelet transform and neural network, an optimization method for structural damage identification based on rotation mode and wavelet neural network is proposed. This method can significantly improve the calculation efficiency of damage identification location and degree, so as to avoid or reduce detection errors to the greatest extent. The results also show some other research conclusions.
The combination of the proposed wavelet and neural network optimization method can provide a valuable tool in the damage identification of the health monitoring of a variety of railway vehicle structures and infrastructures, which is beneficial to the fault diagnosis and monitoring of the vehicle status in a complex load environment. Considering that the practicability of this neural network optimization tool is enhanced by the genetic algorithm optimization theory, the theory makes the analysis equally applicable to other structures. Therefore, the extension of this research is a very meaningful research topic. It not only combines the local characteristics of wavelet with the self-learning characteristics of neural network, but also can be applied to engineering practice.