# A Damage Prognosis Method of Girder Structures Based on Wavelet Neural Networks

Authors: Rumian Zhong; Zhouhong Zong; Jie Niu; Sujing Yuan
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

Based on the basic theory of wavelet neural networks and finite element model updating method, a basic framework of damage prognosis method is proposed in this paper. Firstly, a damaged I-steel beam model testing is used to verify the feasibility and effectiveness of the proposed damage prognosis method. The results show that the predicted results of the damage prognosis method and the measured results are very well consistent, and the maximum error is less than 5%. Furthermore, Xinyihe Bridge in the Beijing-Shanghai Highway is selected as the engineering background, and the damage prognosis is conducted based on the data from the structural health monitoring system. The results show that the traffic volume will increase and seasonal differences will decrease in the next year and a half. The displacement has a slight increase and seasonal characters in the critical section of mid span, but the strain will increase distinctly. The analysis results indicate that the proposed method can be applied to the damage prognosis of girder bridge structures and has the potential for the bridge health monitoring and safety prognosis.


---

## Body

## 1. Introduction

Currently, Structural Health Monitoring (SHM) systems are widely used in bridge monitoring and maintenance management. However, they focus on the data accumulation, safety assessment, and so forth. Then the SHM systems can not satisfy the demands of proprietors, because they are generally limited to the damage detection (DD) level, in which the damage in a structural or mechanical system is herein defined as intentional or unintentional changes to its material and/or geometric properties, including variations of its boundary conditions, which adversely affect its performance [1]. According to Farrar et al. and Inman et al. [2–4], damage prognosis (DP) is the theme of future structural health monitoring and can be defined as the estimate of a system’s remaining useful life. In other words, damage prognosis attempts to forecast the performance of a system by assessing its current state of damage through NDE measurements, by estimating the future loading environments for that system, and by predicting through simulation and past experience its remaining useful life. The damage prognosis problem includes three main critical areas: (i) sensing and processing hardware, (ii) data interrogation, and (iii) modeling and simulation. The essential differences between DP and DD can be distinguished as follows: DD can only find out the existing damage; on the other hand, DP can predict damage which will occur and its consequences. However, DP is more practical than DD, and its successful implementation will allow the owners to take the necessary protections and maintenances in time to avoid the catastrophic failures, rather than waiting until the occurrence of structural damage or failure [5].The structural uncertainty and prognosis reliability are the basic topics of DP. Current prognosis methods for predicting structural failures are classified as conventional reliability models, model-based prognosis models, and data-driven prognosis model [3, 6, 7]. The conventional reliability models are used in complex engineering systems, such as rotating machinery [8] and aircraft structures [9, 10]. A physical-based prognostics model begins by constructing a finite element model, and then it is verified or validated using structural monitoring information, and these information along with outputs from the physics-based models will be used to assess the current state of the structure (e.g., existence, location, and type and extent of damage), which has been proved to predict effectively [11–15]. For most structural applications, physics-based models may not be the most practical solution since the fault type in question is often not unique from component to component. However, physics-based models may be the most suitable approach for cost-justified applications in which accuracy outweighs most other factors and physics models remain consistent across structural systems, such as bridge structures. They also generally require less data than data-driven models.Data-driven approaches attempt to derive models directly from routinely collecting condition-monitored (CM) data instead of building models based on comprehensive system physics and human expertise. They are built based on historical records and produce prediction outputs directly in terms of CM data [16], including the DP methods based on wavelet artificial neural network [17], Bayesian framework [18–20], Kalman estimator [21, 22], fuzzy theory [23], and probability analysis [24, 25]. The DP methods based on data-driven can detect the structural damage with limited input and output measurement signals; however, the integrity of the data cannot be guaranteed. The wavelet artificial neural network method, which combines denoising ability of wavelet analysis and forecasting ability of neural network, is a fast and rapidly convergent iterative optimization algorithm, and it is widely used in various fields [26–29].In this paper, the theoretical basis of wavelet neural network algorithm and its implementation process are described. Then, the wavelet neural network method and model updating method are combined to conduct the DP of structure. Subsequently, the framework of DP is proposed and applied to the DP of a damaged I-steel beam and six-span continuous beam bridge, which means the initial realization of the second step of DP-predicting the future structural loads and structural properties. It also laid a solid foundation to the realization of the third step-predicting the remaining life of the structure.


## 2. The Wavelet Neural Networks

### 2.1. The Wavelet Neural Networks Model

Many problems are nonlinear in the field of civil engineering. The relationships between variables are so complicated that a large number of practical engineering problems are difficult to have exact solutions in math. However, the wavelet neural networks model is relying on an artificial intelligence algorithm, which uses wavelet unit to replace neuron. The model is composed of many interconnected wavelet units, which are used to simulate the formal structure of the human mind. Then the nonlinear dynamical system based on wavelet network can be obtained. Due to the learning ability, the weight parameters of the network can be updated continuously according to the setting rules, and then the wavelet neural network has the abilities of fitting function and dealing with information [30, 31].A wavelet neural networks model is shown in Figure1; it has n nodes in the input layer, l nodes in hidden layer, and only one in the output layer; the network topology is n-l-1, in which the relationship between the inputs and output can be expressed as
(1)y=∑k=1lνkψk=∑k=1lνkψ(∑i=1nwkixi-bkak),
where xi is the input value of the wavelet neural network; wki is the weight function between the kth neuron in hidden layer and the ith unit in input layer; νk is the weight parameter between the kth neuron in hidden layer and the unit in output layer; ψ(·) is the excitation function of wavelet function; ak and bk are the scaling and translation parameters; y is the output value of wavelet neural network.Figure 1
The wavelet neural network model.


### 2.2. A Learning Algorithm of Wavelet Neural Network

The learning algorithm is the core of wavelet neural network, and its purpose is to establish mapping relationship between inputs and outputs by function fitting. A widely applied iterative algorithm called Stochastic Gradient Algorithm (e.g., [26, 32]) is introduced in this paper.


#### 2.2.1. Principle of the Stochastic Gradient Algorithm

A learning sample is given as (x1s,x2s,…,xis,…,xns,ds), where the parameter i is the vector components serial number of input value, the parameter s is the learning sample serial number, and ds is the expected value. Setting ys as the actual output corresponding to the inputs (x1s,x2s,…,xns), the error function E of the wavelet neural network can be defined as follows:
(2)E=12∑s=1S(ys-ds)2=12∑s=1S(∑k=1lνkψ(∑i=1nwkixis-bkak)-ds)2.
The principle of wavelet neural network is to get the minimum of E by updating the parameters of wavelet neural network, and then the output ys which is obtained by fitting can be approximately employed as mapping function to contact inputs and outputs.


#### 2.2.2. The Learning Algorithm

As shown in (1) and (2), the network parameters can be solved based on the Stochastic Gradient Algorithm, and they can be calculated as follows:
(3)wki=wki+αΔwkii=1,2,…,n;k=1,2,…,l,νk=νk+βΔνkk=1,2,…,l,ak=ak+γΔakk=1,2,…,l,bk=bk+ηΔbkk=1,2,…,l,
where α,β,γ,η are learning rates. For simplicity,
(4)u=∑i=1nwkixis-bkak.
Then,
(5)Δwki=∂E∂wki=(y-d)(∑k=1lνkψ′(u))xisakΔνk=∂E∂νk=(y-d)(∑k=1lψ(u))Δak=∂E∂ak=(y-d)(∑k=1lνkψ′(u))(-uak)Δbk=∂E∂bk=(y-d)(∑k=1lνkψ′(u))-1ak,
where ψ′(u) is the first derivative of ψ(u), and then the calculated error, relative error, and the new network parameters can be obtained analytically by substituting (3) into (5). If the relative error can meet the requirements of calculation, the output value will be the optimal solution. The calculation process can be summarized as in Figure 2.Figure 2
The flow chart of wavelet neural network. Note:θ0 and θi are the network parameters matrices and m is determined by the owner’s requests.


#### 2.2.3. Wavelet Networks and Their Parameterizations

(1) The Excitation Functions of Wavelet Neural Network. There are various wavelet functions, such as Harr wavelet, Mever wavelet, Mever Daubechies wavelet, Morlet wavelet, and spline wavelet. Due to the flexibility of choices, there is no uniform standard. Morlet wavelet is a Gaussian wave, and it is widely used and allows adjustment of the relative resolution in time and scale [33, 34]. Morlet wavelet basis function (e.g., [35]) is chosen as the excitation function of wavelet neural network in this paper.(2) The Number of Wavelet Units. At present, there is no theoretical guidance about choosing the number of wavelet units in hidden layer [25]. However, there are lots of empirical formulas, and they can be expressed as(6a)l=m×n(6b)l=m+n+t(6c)l=0.43mn+0.12n2+2.54m+0.77n+0.35+0.51(6d)l=mn+(1/2)n(n2+n)-1m+n,where n is the number of input nodes, m is the number of output nodes, and t is a constant between 1 and 20.Theoretically, there should be an optimal number of wavelet units in hidden layer. However, these empirical formulas are always gained from experiments in other fields, and they can only get a range or approximation instead of an exact figure. Therefore, the experience formulas can only be taken as a reference to determine the number of wavelet units in hidden layer in the field of civil engineering.(3) Parameters Initialization of Wavelet Neural Network. As an example, l is the number of nodes in hidden layer, n is the number of nodes in input layer, wki is the weight parameters between the kth neuron of hidden layer and the ith neuron of input layer, and νk is weight parameters between the kth neuron of hidden layer and the unit of output layer. Initial parameters of wki,νk,ak,bk can be obtained by the experiential analysis of signals; reasonable parameters can reduce the number of iterations; however, the empirical formulas have not been verified widely. Compared with the circulation signals in mechanical engineering, the measured signals from real civil engineering structures are more discrete and stochastic, so initial parameters cannot be defined as those from mechanical engineering [26, 32, 35]. The empirical formulas have been tried to select the initial parameters, but the results are not satisfactory. It is not easy to describe a parameter which contains discreteness and randomness by a deterministic formula. A standard normal distribution random matrix is herein used as the initial parameters; the results were conformed to be accurate and effective.


## 2.1. The Wavelet Neural Networks Model

Many problems are nonlinear in the field of civil engineering. The relationships between variables are so complicated that a large number of practical engineering problems are difficult to have exact solutions in math. However, the wavelet neural networks model is relying on an artificial intelligence algorithm, which uses wavelet unit to replace neuron. The model is composed of many interconnected wavelet units, which are used to simulate the formal structure of the human mind. Then the nonlinear dynamical system based on wavelet network can be obtained. Due to the learning ability, the weight parameters of the network can be updated continuously according to the setting rules, and then the wavelet neural network has the abilities of fitting function and dealing with information [30, 31].A wavelet neural networks model is shown in Figure1; it has n nodes in the input layer, l nodes in hidden layer, and only one in the output layer; the network topology is n-l-1, in which the relationship between the inputs and output can be expressed as
(1)y=∑k=1lνkψk=∑k=1lνkψ(∑i=1nwkixi-bkak),
where xi is the input value of the wavelet neural network; wki is the weight function between the kth neuron in hidden layer and the ith unit in input layer; νk is the weight parameter between the kth neuron in hidden layer and the unit in output layer; ψ(·) is the excitation function of wavelet function; ak and bk are the scaling and translation parameters; y is the output value of wavelet neural network.Figure 1
The wavelet neural network model.


## 2.2. A Learning Algorithm of Wavelet Neural Network

The learning algorithm is the core of wavelet neural network, and its purpose is to establish mapping relationship between inputs and outputs by function fitting. A widely applied iterative algorithm called Stochastic Gradient Algorithm (e.g., [26, 32]) is introduced in this paper.


### 2.2.1. Principle of the Stochastic Gradient Algorithm

A learning sample is given as (x1s,x2s,…,xis,…,xns,ds), where the parameter i is the vector components serial number of input value, the parameter s is the learning sample serial number, and ds is the expected value. Setting ys as the actual output corresponding to the inputs (x1s,x2s,…,xns), the error function E of the wavelet neural network can be defined as follows:
(2)E=12∑s=1S(ys-ds)2=12∑s=1S(∑k=1lνkψ(∑i=1nwkixis-bkak)-ds)2.
The principle of wavelet neural network is to get the minimum of E by updating the parameters of wavelet neural network, and then the output ys which is obtained by fitting can be approximately employed as mapping function to contact inputs and outputs.


### 2.2.2. The Learning Algorithm

As shown in (1) and (2), the network parameters can be solved based on the Stochastic Gradient Algorithm, and they can be calculated as follows:
(3)wki=wki+αΔwkii=1,2,…,n;k=1,2,…,l,νk=νk+βΔνkk=1,2,…,l,ak=ak+γΔakk=1,2,…,l,bk=bk+ηΔbkk=1,2,…,l,
where α,β,γ,η are learning rates. For simplicity,
(4)u=∑i=1nwkixis-bkak.
Then,
(5)Δwki=∂E∂wki=(y-d)(∑k=1lνkψ′(u))xisakΔνk=∂E∂νk=(y-d)(∑k=1lψ(u))Δak=∂E∂ak=(y-d)(∑k=1lνkψ′(u))(-uak)Δbk=∂E∂bk=(y-d)(∑k=1lνkψ′(u))-1ak,
where ψ′(u) is the first derivative of ψ(u), and then the calculated error, relative error, and the new network parameters can be obtained analytically by substituting (3) into (5). If the relative error can meet the requirements of calculation, the output value will be the optimal solution. The calculation process can be summarized as in Figure 2.Figure 2
The flow chart of wavelet neural network. Note:θ0 and θi are the network parameters matrices and m is determined by the owner’s requests.


### 2.2.3. Wavelet Networks and Their Parameterizations

(1) The Excitation Functions of Wavelet Neural Network. There are various wavelet functions, such as Harr wavelet, Mever wavelet, Mever Daubechies wavelet, Morlet wavelet, and spline wavelet. Due to the flexibility of choices, there is no uniform standard. Morlet wavelet is a Gaussian wave, and it is widely used and allows adjustment of the relative resolution in time and scale [33, 34]. Morlet wavelet basis function (e.g., [35]) is chosen as the excitation function of wavelet neural network in this paper.(2) The Number of Wavelet Units. At present, there is no theoretical guidance about choosing the number of wavelet units in hidden layer [25]. However, there are lots of empirical formulas, and they can be expressed as(6a)l=m×n(6b)l=m+n+t(6c)l=0.43mn+0.12n2+2.54m+0.77n+0.35+0.51(6d)l=mn+(1/2)n(n2+n)-1m+n,where n is the number of input nodes, m is the number of output nodes, and t is a constant between 1 and 20.Theoretically, there should be an optimal number of wavelet units in hidden layer. However, these empirical formulas are always gained from experiments in other fields, and they can only get a range or approximation instead of an exact figure. Therefore, the experience formulas can only be taken as a reference to determine the number of wavelet units in hidden layer in the field of civil engineering.(3) Parameters Initialization of Wavelet Neural Network. As an example, l is the number of nodes in hidden layer, n is the number of nodes in input layer, wki is the weight parameters between the kth neuron of hidden layer and the ith neuron of input layer, and νk is weight parameters between the kth neuron of hidden layer and the unit of output layer. Initial parameters of wki,νk,ak,bk can be obtained by the experiential analysis of signals; reasonable parameters can reduce the number of iterations; however, the empirical formulas have not been verified widely. Compared with the circulation signals in mechanical engineering, the measured signals from real civil engineering structures are more discrete and stochastic, so initial parameters cannot be defined as those from mechanical engineering [26, 32, 35]. The empirical formulas have been tried to select the initial parameters, but the results are not satisfactory. It is not easy to describe a parameter which contains discreteness and randomness by a deterministic formula. A standard normal distribution random matrix is herein used as the initial parameters; the results were conformed to be accurate and effective.


## 2.2.1. Principle of the Stochastic Gradient Algorithm

A learning sample is given as (x1s,x2s,…,xis,…,xns,ds), where the parameter i is the vector components serial number of input value, the parameter s is the learning sample serial number, and ds is the expected value. Setting ys as the actual output corresponding to the inputs (x1s,x2s,…,xns), the error function E of the wavelet neural network can be defined as follows:
(2)E=12∑s=1S(ys-ds)2=12∑s=1S(∑k=1lνkψ(∑i=1nwkixis-bkak)-ds)2.
The principle of wavelet neural network is to get the minimum of E by updating the parameters of wavelet neural network, and then the output ys which is obtained by fitting can be approximately employed as mapping function to contact inputs and outputs.


## 2.2.2. The Learning Algorithm

As shown in (1) and (2), the network parameters can be solved based on the Stochastic Gradient Algorithm, and they can be calculated as follows:
(3)wki=wki+αΔwkii=1,2,…,n;k=1,2,…,l,νk=νk+βΔνkk=1,2,…,l,ak=ak+γΔakk=1,2,…,l,bk=bk+ηΔbkk=1,2,…,l,
where α,β,γ,η are learning rates. For simplicity,
(4)u=∑i=1nwkixis-bkak.
Then,
(5)Δwki=∂E∂wki=(y-d)(∑k=1lνkψ′(u))xisakΔνk=∂E∂νk=(y-d)(∑k=1lψ(u))Δak=∂E∂ak=(y-d)(∑k=1lνkψ′(u))(-uak)Δbk=∂E∂bk=(y-d)(∑k=1lνkψ′(u))-1ak,
where ψ′(u) is the first derivative of ψ(u), and then the calculated error, relative error, and the new network parameters can be obtained analytically by substituting (3) into (5). If the relative error can meet the requirements of calculation, the output value will be the optimal solution. The calculation process can be summarized as in Figure 2.Figure 2
The flow chart of wavelet neural network. Note:θ0 and θi are the network parameters matrices and m is determined by the owner’s requests.


## 2.2.3. Wavelet Networks and Their Parameterizations

(1) The Excitation Functions of Wavelet Neural Network. There are various wavelet functions, such as Harr wavelet, Mever wavelet, Mever Daubechies wavelet, Morlet wavelet, and spline wavelet. Due to the flexibility of choices, there is no uniform standard. Morlet wavelet is a Gaussian wave, and it is widely used and allows adjustment of the relative resolution in time and scale [33, 34]. Morlet wavelet basis function (e.g., [35]) is chosen as the excitation function of wavelet neural network in this paper.(2) The Number of Wavelet Units. At present, there is no theoretical guidance about choosing the number of wavelet units in hidden layer [25]. However, there are lots of empirical formulas, and they can be expressed as(6a)l=m×n(6b)l=m+n+t(6c)l=0.43mn+0.12n2+2.54m+0.77n+0.35+0.51(6d)l=mn+(1/2)n(n2+n)-1m+n,where n is the number of input nodes, m is the number of output nodes, and t is a constant between 1 and 20.Theoretically, there should be an optimal number of wavelet units in hidden layer. However, these empirical formulas are always gained from experiments in other fields, and they can only get a range or approximation instead of an exact figure. Therefore, the experience formulas can only be taken as a reference to determine the number of wavelet units in hidden layer in the field of civil engineering.(3) Parameters Initialization of Wavelet Neural Network. As an example, l is the number of nodes in hidden layer, n is the number of nodes in input layer, wki is the weight parameters between the kth neuron of hidden layer and the ith neuron of input layer, and νk is weight parameters between the kth neuron of hidden layer and the unit of output layer. Initial parameters of wki,νk,ak,bk can be obtained by the experiential analysis of signals; reasonable parameters can reduce the number of iterations; however, the empirical formulas have not been verified widely. Compared with the circulation signals in mechanical engineering, the measured signals from real civil engineering structures are more discrete and stochastic, so initial parameters cannot be defined as those from mechanical engineering [26, 32, 35]. The empirical formulas have been tried to select the initial parameters, but the results are not satisfactory. It is not easy to describe a parameter which contains discreteness and randomness by a deterministic formula. A standard normal distribution random matrix is herein used as the initial parameters; the results were conformed to be accurate and effective.


## 3. The Framework of Damage Prognosis Based on Wavelet Neural Network and Model Updating

Damage prognosis combines model-based method with date-driven method. Namely, the structural responses are calculated based on the updated finite element model and used as sample values which are needed by data-driven method [2–4]. Thus, the advantages of the two methods can be combined effectively in order to get a suitable DP model. It can be applied to evaluate the structure safety and to put forward the best maintenance plan. The basic process is shown in Figure 3. The process begins by collecting as much initial system information as possible. This information is used to develop model-based numerical model of the structure as well as to define the system that will be used for state awareness assessments and whatever sensors needed to monitor operational and environment conditions. As the data become available from the SHM system, they will be used to validate and update the FE model, and the updated and validated model can be used to develop the DP model based on the wavelet neural network method, which is the key part of the framework; the other part is damage identification, where the data of SHM system will also be used to assess the current state of the structure (existence, location, type, and extent of damage). Data from the operational and environmental sensors will be used to predict the future system loading. The output of the future loading model, state awareness model, and the DP model will be used to estimate the future state of the structure. As is indicated, the solution process will be iterative, relying on experience gained from past predictions to improve future predictions.Figure 3
A general damage prognosis solution procedure based on wavelet neural network and model updating.


## 4. Damage Prognosis of Steel Beam Model

At present, there are rarely experiments about DP in the field of civil engineering. In this paper, through a damaged I-steel beam model testing, the damage evolution of the dynamic property is explored. Finally, the prediction results based on the wavelet neural networks method are compared with the ambient vibration testing results in order to verify the validity of the proposed DP method.


### 4.1. The Model Testing

This testing includes four I-steel beams, each one with a length of 3 meters. The density of the steel material is 7800 kg/m3, and the elasticity modulus is 2.1 × 105 MPa. In the I-section, area is 14.33 cm2, inertia moment is 223 cm4, and the damage section inertia moment is 10.5 cm4. Among the four beams in this testing, one beam is undamaged, marked as beam-0; one has a notch in the L/2 of the beam, marked as beam-1; one has two notches in the L/2 and the L/4 of the beam, marked as beam-2; and the last one has 3 notches in the middle, a quarter, and 3 quarters of the beam, marked as beam-3. The notches of the beams got the same length of 0.1 m and a height of 0.05 m (half of the beam height) (Figure 4). The ambient vibration testing of beam-2 is shown in Figure 5; data was sampled at 1000 Hz, and each setup for all tests was recorded for duration of 15 minutes. The vertical frequencies can be obtained by testing, as shown in Table 1.Table 1
Comparison of updated frequencies and measured frequencies, MAC values (unit, Hz).
Model number
Measured frequencies ①
Updated frequencies ②
Relative error(②−①)/① (%)
MAC values (%)
beam-0
First order
33.62
32.78
−2.5
93.1
Second order
135.03
130.78
−3.1
91.4
Third order
293.03
293.02
0
95.5
beam-1
First order
20.71
22.02
6.3
95.1
Second order
131.02
130.04
−0.7
93.4
Third order
231.18
232.69
−0.7
90.9
beam-2
First order
18.90
19.46
2.9
96.7
Second order
85.94
90.67
5.5
94.3
Third order
202.54
203.07
0.3
90.8
beam-3
First order
17.02
17.29
1.6
94.7
Second order
66.69
67.14
0.7
94.9
Third order
148.33
149.68
0.9
90.2Figure 4
The notch of steel beam.Figure 5
The environmental vibration testing.


### 4.2. The Finite Element Model Updating of Steel Beam

Using ANSYS finite element (FE) analysis software, these FE beam models are established by beam-3 elements, and it is composed of 60 elements and 61 nodes, as shown in Figure6. The first three order frequencies are selected as the response features; I1, I2, and I3 are the inertia moment of damage elements in the middle, a quarter, and 3 quarters of the beam. Then FE model updating is conducted based on third-order response surface method [36], and the updated frequencies, measured frequencies, and MAC values are shown in Table 1.Figure 6
Numerical model (unit, m).


### 4.3. Determine the Network Structure and Initial Parameters

Wavelet neural network training is conducted according to the method of this paper, and Morlet wavelet is selected as the excitation function; the network training is effective when the number of wavelet units in hidden layer is 3 (as (6b), l=3+1+1=3, and t=1), the number of training is set as 10000 times, error m=0.001, and the training rate α=β=0.01;γ=η=0.001.According to the standard normal distribution, the initial parameters are selected randomly as follows:(7)W={-0.8242.288-0.8521.8740.3750.072-4.427-1.9394.941};V={1.125-1.3523.647};A={1.807-0.0420.404};B={0.572-0.1701.985}.
Note: “W,” “V,” “A,” and “B” are the parameters of wki, νk, ak, and bk.


### 4.4. Wavelet Neural Network Training and Damage Prognosis Model

The future fundamental frequency of a damaged I-steel beam was selected as the output value in this paper. Firstly, the experimental samples were selected by D-optimal design method, and the damage degree was simulated by cutting section stiffness. Then fundamental frequencies of beams under different damage degrees were calculated through the updated FE model. In order to validate the network, three beams were employed as the test samples, as shown in Tables2 and 3.Table 2
Training samples of the experimental design.
N
Inputs of section stiffness (×10−6 m4)
Output of the fundamental frequency (Hz)
Damage section-1
Damage section-2
Damage section-3
1
0.93
0.92
2.2
30.67
2
0.05
0.05
1.58
14.07
3
1.21
1.4
0.45
29.91
⋮
⋮
⋮
⋮
⋮
28
0.05
2.23
0.72
20.18
29
0.05
0.73
2.23
19.91
30
0.05
0.05
0.05
12.56Table 3
The test samples.
N
Inputs of section stiffness (×10−6 m4)
Outputs of fundamental frequency (Hz)
1
2.23
0.105
2.23
22.018
2
0.105
0.105
2.23
19.455
3
0.105
0.105
0.105
17.288
Note: “N” is the number of the samples.The training samples and test samples are put into the wavelet neural networks solution procedure based on MATLAB 7.0 software. Then adaptive training was realized, and suitable weight parameters can be obtained. Furthermore, the DP function model can be established.


### 4.5. The Result of Damage Prognosis

As the DP function model has been established, the fundamental frequencies of these three steel beams can be predicted relying on the function model. The DP results based on wavelet neural network are shown in Table4. It is shown that the predicted results based on wavelet neural network and the measured results from model testing are very well consistent, and the error is less than 5%.Table 4
The results of damage prognosis.
N
The measured valuesf0
Theoretical valuesf1
Outputvalues f2
The error rate with theoretical values|(f2-f1)/f1|
The error rate with measured values|(f2-f0)/f1|
1
20.713
22.018
21.72
1.35%
4.86%
2
18.903
19.455
19.59
0.70%
3.63%
3
17.017
17.288
16.76
3.02%
−1.51%


## 4.1. The Model Testing

This testing includes four I-steel beams, each one with a length of 3 meters. The density of the steel material is 7800 kg/m3, and the elasticity modulus is 2.1 × 105 MPa. In the I-section, area is 14.33 cm2, inertia moment is 223 cm4, and the damage section inertia moment is 10.5 cm4. Among the four beams in this testing, one beam is undamaged, marked as beam-0; one has a notch in the L/2 of the beam, marked as beam-1; one has two notches in the L/2 and the L/4 of the beam, marked as beam-2; and the last one has 3 notches in the middle, a quarter, and 3 quarters of the beam, marked as beam-3. The notches of the beams got the same length of 0.1 m and a height of 0.05 m (half of the beam height) (Figure 4). The ambient vibration testing of beam-2 is shown in Figure 5; data was sampled at 1000 Hz, and each setup for all tests was recorded for duration of 15 minutes. The vertical frequencies can be obtained by testing, as shown in Table 1.Table 1
Comparison of updated frequencies and measured frequencies, MAC values (unit, Hz).
Model number
Measured frequencies ①
Updated frequencies ②
Relative error(②−①)/① (%)
MAC values (%)
beam-0
First order
33.62
32.78
−2.5
93.1
Second order
135.03
130.78
−3.1
91.4
Third order
293.03
293.02
0
95.5
beam-1
First order
20.71
22.02
6.3
95.1
Second order
131.02
130.04
−0.7
93.4
Third order
231.18
232.69
−0.7
90.9
beam-2
First order
18.90
19.46
2.9
96.7
Second order
85.94
90.67
5.5
94.3
Third order
202.54
203.07
0.3
90.8
beam-3
First order
17.02
17.29
1.6
94.7
Second order
66.69
67.14
0.7
94.9
Third order
148.33
149.68
0.9
90.2Figure 4
The notch of steel beam.Figure 5
The environmental vibration testing.


## 4.2. The Finite Element Model Updating of Steel Beam

Using ANSYS finite element (FE) analysis software, these FE beam models are established by beam-3 elements, and it is composed of 60 elements and 61 nodes, as shown in Figure6. The first three order frequencies are selected as the response features; I1, I2, and I3 are the inertia moment of damage elements in the middle, a quarter, and 3 quarters of the beam. Then FE model updating is conducted based on third-order response surface method [36], and the updated frequencies, measured frequencies, and MAC values are shown in Table 1.Figure 6
Numerical model (unit, m).


## 4.3. Determine the Network Structure and Initial Parameters

Wavelet neural network training is conducted according to the method of this paper, and Morlet wavelet is selected as the excitation function; the network training is effective when the number of wavelet units in hidden layer is 3 (as (6b), l=3+1+1=3, and t=1), the number of training is set as 10000 times, error m=0.001, and the training rate α=β=0.01;γ=η=0.001.According to the standard normal distribution, the initial parameters are selected randomly as follows:(7)W={-0.8242.288-0.8521.8740.3750.072-4.427-1.9394.941};V={1.125-1.3523.647};A={1.807-0.0420.404};B={0.572-0.1701.985}.
Note: “W,” “V,” “A,” and “B” are the parameters of wki, νk, ak, and bk.


## 4.4. Wavelet Neural Network Training and Damage Prognosis Model

The future fundamental frequency of a damaged I-steel beam was selected as the output value in this paper. Firstly, the experimental samples were selected by D-optimal design method, and the damage degree was simulated by cutting section stiffness. Then fundamental frequencies of beams under different damage degrees were calculated through the updated FE model. In order to validate the network, three beams were employed as the test samples, as shown in Tables2 and 3.Table 2
Training samples of the experimental design.
N
Inputs of section stiffness (×10−6 m4)
Output of the fundamental frequency (Hz)
Damage section-1
Damage section-2
Damage section-3
1
0.93
0.92
2.2
30.67
2
0.05
0.05
1.58
14.07
3
1.21
1.4
0.45
29.91
⋮
⋮
⋮
⋮
⋮
28
0.05
2.23
0.72
20.18
29
0.05
0.73
2.23
19.91
30
0.05
0.05
0.05
12.56Table 3
The test samples.
N
Inputs of section stiffness (×10−6 m4)
Outputs of fundamental frequency (Hz)
1
2.23
0.105
2.23
22.018
2
0.105
0.105
2.23
19.455
3
0.105
0.105
0.105
17.288
Note: “N” is the number of the samples.The training samples and test samples are put into the wavelet neural networks solution procedure based on MATLAB 7.0 software. Then adaptive training was realized, and suitable weight parameters can be obtained. Furthermore, the DP function model can be established.


## 4.5. The Result of Damage Prognosis

As the DP function model has been established, the fundamental frequencies of these three steel beams can be predicted relying on the function model. The DP results based on wavelet neural network are shown in Table4. It is shown that the predicted results based on wavelet neural network and the measured results from model testing are very well consistent, and the error is less than 5%.Table 4
The results of damage prognosis.
N
The measured valuesf0
Theoretical valuesf1
Outputvalues f2
The error rate with theoretical values|(f2-f1)/f1|
The error rate with measured values|(f2-f0)/f1|
1
20.713
22.018
21.72
1.35%
4.86%
2
18.903
19.455
19.59
0.70%
3.63%
3
17.017
17.288
16.76
3.02%
−1.51%


## 5. Damage Prognosis Analysis of Xinyihe Bridge

### 5.1. Engineering Background

Xinyihe Bridge (Figure7) of the Beijing-Shanghai Highway is located in Shuyang, Jiangsu province. The upper structure is a six-span continuous beam bridge, and design loads are car-super level 20 and trailer-120. It was built in 2001 and has served for 12 years.Xinyihe Bridge.
(a)
Xinyihe Bridge panorama photo
(b)
The weigh-in-motion (WIM)


### 5.2. Load Prediction

The traffic data used were monthly average number of cars during June 2012 to June 2013 (Figure8); the DP model based on wavelet neural network method could provide short-term predictive traffic information (one month). According to the historical traffic data, it is obvious that the traffic intensity in a given month (n) is related to the last 3 months’ (n-1, n-2, and n-3) traffic intensity (Qn-1, Qn-2, and Qn-3) and the given month (n), so the 4 parameters (Qn-1, Qn-2, Qn-3, and n) could be the input values, and the traffic intensity Qn will be the output value. Based on the wavelet neural network method and the traffic data, the traffic intensity prediction could be realized. As shown in Figures 8 and 9, within the next year and a half, the traffic intensity will increase and seasonal differences will decrease; the prediction model can be employed for the bridge safety evaluation of Beijing-Shanghai Highway. As the WIM continues to collect the new data, the samples will be refreshed constantly, and the load prediction model could be updated continuously.Figure 8
Historical traffic intensity (2012.06–2013.06).Figure 9
The future traffic intensity (2013.07–2014.12).


### 5.3. The Model Updated of Xinyihe Bridge

According to the design papers, the FE model was built. Based on the twice ambient vibration testing in 2012 and 2013, the FE model was updated according to the third-order response surface method [36]. The updated frequencies are shown in Table 5, and the updated parameters are shown in Table 6.Table 5
Parameters change before and after updating based on the response surface method.
Updated parameters
E
0
E
1
E
2
E
3
K
1
K
2
K
3
(×104 MPa)
(×104 MPa)
(×104 MPa)
(×104 MPa)
(×106 N/m)
(×106 N/m)
(×106 N/m)
2012
3.6
3.5
3.2
2.5
0.6
3.1
3.1
2013
3.3
2.8
3.2
2.2
0.9
4.7
1.7
Rate (%)
−8.8
−19.8
−1.9
−12.4
50.0
48.4
−43.3
Notice:E1 is the elasticity modulus of webs in the 2nd, 3rd, and 5th spans; E2 is the elasticity modulus of web in the 4th span; E3 is the elasticity modulus of bottom slab in the 2nd, 5th spans; E0 is the elasticity modulus of the other parts; K1 is transverse spring stiffness at the support and expansion joints; K2 is the longitudinal spring stiffness at the support; K3 is the longitudinal spring stiffness at the expansion joints.Table 6
Comparison of updated frequencies and measured frequencies (unit, Hz).
Vibration mode
Frequency (Hz)
Measuredvalue-1 (SSI)
Updatedvalue-1
Error-1 (%)
Measuredvalue-2 (SSI)
Updatedvalue-2
Error-2 (%)
Vertical
First order
2.89
2.88
0.43
2.87
2.77
3.48
Second order
3.03
3.081
1.84
3.16
2.89
8.66
Third order
3.79
3.61
0.51
3.67
3.34
8.99
Fourth order
4.26
4.32
1.27
4.25
4.02
5.37
Transverse
First order
0.83
0.73
0.16
0.92
0.87
4.90
Second order
1.43
1.47
2.84
1.45
1.48
1.45
Longitudinal
First order
1.79
1.74
2.45
1.30
1.28
1.61
Notice: measured value-1, updated value-1, and error-1 were the updated results in February 2012; measured value-2, updated value-2, and error-2 were the updated results in June 2013.As mentioned above, the FE model could be updated continuously based on the results of ambient vibration testing in the future. If there are enough frequencies samples, the wavelet neural network will be used for the model parameters prediction, namely, the prediction of the FE model. Due to that there is no real-time dynamics monitoring, we could not have enough samples, and then the model parameters change rate between the two updated FE models will be used as the actual structure performance degradation rate in this paper.


### 5.4. The Damage Prognosis Model

The maximum values of displacement and strain were separately selected as outputs of the prognosis model in this paper, and the parameters (E0, E1, E2, and E3) and load (F) were selected as the inputs. Similarly to the previously damaged I-steel beams, the parameters samples were selected by D-optimal design method, and the training samples could be obtained based on the finite element analysis model, as shown in Tables 7 and 8.Table 7
Training samples.
N
E
0 (×104 MPa)
E
1 (×104 MPa)
E
2 (×104 MPa)
E
3 (×104 MPa)
Vehicle loadF
The maximum value of displacement
The maximum value of strain
1
2.5
1.5
2
1.5
2.0
35.700
443.51
2
2.5
3.3
2
3
2.0
28.316
347.75
⋮
⋮
⋮
⋮
⋮
⋮
⋮
30
3.5
1.5
2.6
1.5
1.1
16.164
198.27
31
2.5
3.3
3.5
1.5
2.0
28.969
352.76
Notice: the value of vehicle load (F) is the n times of highway-I load.Table 8
Verification of the results of damage prognosis model.
Measure value
Measured values
Outputs of network
Error (%)
The maximum value of displacement (mm)
15.37
16.77
9.11
The maximum value of strain (με)
534
475.2
11.01The training samples are substituted into the MATLAB model based on wavelet neural network; the DP model can be established relying on the self-training of the prognosis model, and the goal is the error parameterE<m.


### 5.5. Validation and Application of Prognosis Results

The accuracy of the prognosis model is verified by comparing the prediction results with the static test results in June 2013, as shown in Figure10. The maximum values of displacement and strain under highway-I load can be obtained based on the static test, and this damage degree is similar to the state of the second ambient vibration test. The results were compared with the prognosis results based on wavelet neural network, as shown in Table 8. It is clear that the prediction results of DP and the measured results of static test are very well consistent, and the errors of the maximum values of displacement and strain are less than 12% (seen in Table 8).Figure 10
The static and dynamic test of Xinyihe Bridge.As shown in Table8, the DP model can be used in the bridge health monitoring system, and then based on the load prognosis model and the updated model, the maximum values of displacement and strain under working load conditions can be predicted from July 2013 to December 2014, as shown in Figures 11 and 12. The results show that the traffic volume will increase and seasonal differences will decrease in the next year and a half; the displacement has a slight increase and seasonal characters in the critical section of mid span, but the strain increases distinctly, and timely remedial actions need to be taken.Figure 11
The maximum value of displacement (2013.07–2014.12).Figure 12
The maximum value of strain (2013.07–2014.12).


## 5.1. Engineering Background

Xinyihe Bridge (Figure7) of the Beijing-Shanghai Highway is located in Shuyang, Jiangsu province. The upper structure is a six-span continuous beam bridge, and design loads are car-super level 20 and trailer-120. It was built in 2001 and has served for 12 years.Xinyihe Bridge.
(a)
Xinyihe Bridge panorama photo
(b)
The weigh-in-motion (WIM)


## 5.2. Load Prediction

The traffic data used were monthly average number of cars during June 2012 to June 2013 (Figure8); the DP model based on wavelet neural network method could provide short-term predictive traffic information (one month). According to the historical traffic data, it is obvious that the traffic intensity in a given month (n) is related to the last 3 months’ (n-1, n-2, and n-3) traffic intensity (Qn-1, Qn-2, and Qn-3) and the given month (n), so the 4 parameters (Qn-1, Qn-2, Qn-3, and n) could be the input values, and the traffic intensity Qn will be the output value. Based on the wavelet neural network method and the traffic data, the traffic intensity prediction could be realized. As shown in Figures 8 and 9, within the next year and a half, the traffic intensity will increase and seasonal differences will decrease; the prediction model can be employed for the bridge safety evaluation of Beijing-Shanghai Highway. As the WIM continues to collect the new data, the samples will be refreshed constantly, and the load prediction model could be updated continuously.Figure 8
Historical traffic intensity (2012.06–2013.06).Figure 9
The future traffic intensity (2013.07–2014.12).


## 5.3. The Model Updated of Xinyihe Bridge

According to the design papers, the FE model was built. Based on the twice ambient vibration testing in 2012 and 2013, the FE model was updated according to the third-order response surface method [36]. The updated frequencies are shown in Table 5, and the updated parameters are shown in Table 6.Table 5
Parameters change before and after updating based on the response surface method.
Updated parameters
E
0
E
1
E
2
E
3
K
1
K
2
K
3
(×104 MPa)
(×104 MPa)
(×104 MPa)
(×104 MPa)
(×106 N/m)
(×106 N/m)
(×106 N/m)
2012
3.6
3.5
3.2
2.5
0.6
3.1
3.1
2013
3.3
2.8
3.2
2.2
0.9
4.7
1.7
Rate (%)
−8.8
−19.8
−1.9
−12.4
50.0
48.4
−43.3
Notice:E1 is the elasticity modulus of webs in the 2nd, 3rd, and 5th spans; E2 is the elasticity modulus of web in the 4th span; E3 is the elasticity modulus of bottom slab in the 2nd, 5th spans; E0 is the elasticity modulus of the other parts; K1 is transverse spring stiffness at the support and expansion joints; K2 is the longitudinal spring stiffness at the support; K3 is the longitudinal spring stiffness at the expansion joints.Table 6
Comparison of updated frequencies and measured frequencies (unit, Hz).
Vibration mode
Frequency (Hz)
Measuredvalue-1 (SSI)
Updatedvalue-1
Error-1 (%)
Measuredvalue-2 (SSI)
Updatedvalue-2
Error-2 (%)
Vertical
First order
2.89
2.88
0.43
2.87
2.77
3.48
Second order
3.03
3.081
1.84
3.16
2.89
8.66
Third order
3.79
3.61
0.51
3.67
3.34
8.99
Fourth order
4.26
4.32
1.27
4.25
4.02
5.37
Transverse
First order
0.83
0.73
0.16
0.92
0.87
4.90
Second order
1.43
1.47
2.84
1.45
1.48
1.45
Longitudinal
First order
1.79
1.74
2.45
1.30
1.28
1.61
Notice: measured value-1, updated value-1, and error-1 were the updated results in February 2012; measured value-2, updated value-2, and error-2 were the updated results in June 2013.As mentioned above, the FE model could be updated continuously based on the results of ambient vibration testing in the future. If there are enough frequencies samples, the wavelet neural network will be used for the model parameters prediction, namely, the prediction of the FE model. Due to that there is no real-time dynamics monitoring, we could not have enough samples, and then the model parameters change rate between the two updated FE models will be used as the actual structure performance degradation rate in this paper.


## 5.4. The Damage Prognosis Model

The maximum values of displacement and strain were separately selected as outputs of the prognosis model in this paper, and the parameters (E0, E1, E2, and E3) and load (F) were selected as the inputs. Similarly to the previously damaged I-steel beams, the parameters samples were selected by D-optimal design method, and the training samples could be obtained based on the finite element analysis model, as shown in Tables 7 and 8.Table 7
Training samples.
N
E
0 (×104 MPa)
E
1 (×104 MPa)
E
2 (×104 MPa)
E
3 (×104 MPa)
Vehicle loadF
The maximum value of displacement
The maximum value of strain
1
2.5
1.5
2
1.5
2.0
35.700
443.51
2
2.5
3.3
2
3
2.0
28.316
347.75
⋮
⋮
⋮
⋮
⋮
⋮
⋮
30
3.5
1.5
2.6
1.5
1.1
16.164
198.27
31
2.5
3.3
3.5
1.5
2.0
28.969
352.76
Notice: the value of vehicle load (F) is the n times of highway-I load.Table 8
Verification of the results of damage prognosis model.
Measure value
Measured values
Outputs of network
Error (%)
The maximum value of displacement (mm)
15.37
16.77
9.11
The maximum value of strain (με)
534
475.2
11.01The training samples are substituted into the MATLAB model based on wavelet neural network; the DP model can be established relying on the self-training of the prognosis model, and the goal is the error parameterE<m.


## 5.5. Validation and Application of Prognosis Results

The accuracy of the prognosis model is verified by comparing the prediction results with the static test results in June 2013, as shown in Figure10. The maximum values of displacement and strain under highway-I load can be obtained based on the static test, and this damage degree is similar to the state of the second ambient vibration test. The results were compared with the prognosis results based on wavelet neural network, as shown in Table 8. It is clear that the prediction results of DP and the measured results of static test are very well consistent, and the errors of the maximum values of displacement and strain are less than 12% (seen in Table 8).Figure 10
The static and dynamic test of Xinyihe Bridge.As shown in Table8, the DP model can be used in the bridge health monitoring system, and then based on the load prognosis model and the updated model, the maximum values of displacement and strain under working load conditions can be predicted from July 2013 to December 2014, as shown in Figures 11 and 12. The results show that the traffic volume will increase and seasonal differences will decrease in the next year and a half; the displacement has a slight increase and seasonal characters in the critical section of mid span, but the strain increases distinctly, and timely remedial actions need to be taken.Figure 11
The maximum value of displacement (2013.07–2014.12).Figure 12
The maximum value of strain (2013.07–2014.12).


## 6. Conclusions

(1)
A damage prognosis framework for bridge structure is proposed combining the wavelet neural network method with the finite element model updating. The effectiveness of the proposed method is verified through a damaged I-section steel beam testing, and the maximum error is less than 5%.
(2)
The prognosis method can be used to predict future traffic load volume based on the traffic load monitoring. The traffic intensity predictions of Xinyihe Bridge show that the traffic volume will increase and seasonal differences will decrease in the next year and a half.
(3)
The prognosis prediction results agree very well with those from static load testing of Xinyihe Bridge. The errors between the predicted and measured results were within the range of 12%. The displacement has slight seasonal characters in the mid span, and the strains increase distinctly. The prognosis model can be used for the bridge safety prognosis and future maintenance.
(4)
This paper only solves the problem of the predicting of the future structural loads and structural properties. The remaining life prediction of real bridge structure is the next attention of our research. The proposed damage prognosis method can be incorporated into the structure health monitoring system, for the purpose of the online safety prognosis, and cost-efficient condition based maintenance of bridge structures.



---

*Source: 130274-2014-04-06.xml*