3.2. On-Site Measurement Campaign before Damage
- a.
Description of experiment
In the framework of this study, a comprehensive survey and measurement campaign was carried out on span 08 of Chuong Duong bridge. The campaign includes geometric surveys and vibration measurements of the entire span 08. After finishing the geometric survey, conduct vibration measurement under random stimuli (wind, current, surrounding loads, vehicles crossing the bridge, among others). Eight highly sensitive sensors, ranging in sensitivity from 1054 to 1083mV/m/s2 were used. With a sampling frequency of 1651 Hz, each setup’s average acquisition time was 30 minutes. Due to a shortage of sensors and in order to determine the global vibration modes of the bridge. To achieve the vibration modes of the entire span, a measuring grid of 34 points covering all truss nodes was established. However, due to the limited number of sensors, the measuring grid is divided into 8 setups (
Figure 4). The reference point at node 103 is selected to link the data from the sensors. Other moving points are located at truss nodes on the bridge. The measurement procedure is controlled by a laptop, which also gathers and saves dynamic responses.
In space, the Cartesian coordinate system is used to determine the direction of each measuring point. The x-axis is in the longitudinal direction, the y-axis is in the bridge’s transverse direction, and the z-axis is vertical with a positive upward direction. Two sensors with x and y axes were mounted at each bearing. The sensors were positioned on the y-axis, z-axis, or both y and z-axis at other points.
- b.
Data processing and feature extraction
The MACEC toolbox [
46] was used to analyze all of the measurement data. First, the data needs to be pre-processed. A measuring grid is created on the MACEC system. The measuring points are assigned and numbered, corresponding to the actual measuring points. Assign input parameters such as sensor label, sensitivity, data and measure to each measuring point. The measured signal data is often skewed and does not coincide with the balance axis; the remove-offset function removes these components from the measurement data. The obtained dynamic signal was taken from the time domain and represented in the frequency domain using the Fast Fourier Transform (FFT) (
Figure 5).
Based on input data and data obtained from pre-processing (noise removal, data classification into corresponding nodes), a model with complete data for measurement points is formed. System identification is accomplished using the covariance-based stochastic subspace identification (SSI-COV) technique. Based on knowledge from numerous similar constructs, the following criteria were selected to concretize and characterize the modality: frequency stabilization (1%), damping ratio stabilization (5%), and mode shape stabilization (1%). The stabilization diagrams (
Figure 6) are constructed. The choice of stable poles (in red dash-line) were based on their obvious appearance in 8 setups.
After data processing, seven identified modes shape is identified from the campaigns.
Figure 7 displays mode shapes from 8 setups:
From the different setups, 7 modes could be identified. Their natural frequencies and damping ratios can be found in the
Table 1.
The standard deviation of the natural frequency is calculated to assess the defined modes’ effectiveness. Because the values of std.f are low, each setting’s system recognition quality is high. Modal phase alignment (MPC) measures the mode shape’s departure from actual values, MPC=1 corresponds to the pure real mode. Every MPC value is higher than 0.998. A structure with light and/or proportional damping physics modes is likely to be realistic, so the elevated MPC result typically indicates a mode that has been precisely defined.
3.3. FEM Creation and Updating
A FE model of Chuongduong bridge is built while taking into account the bridge’s structure (
Figure 8). Top chords, bottom chords, cantilevers, gate frames, and verticals, all of which were modeled using three-dimensional beam elements, make up the principal structural components. Other components such as wind bracing, stiffening frame, and longitudinal linkage are modeled with truss elements. The global X-axis of the bridge is in its longitudinal direction, the Z-axis is vertical, and the Y-axis is transverse (to the direction of the river flow). The built model consists of 619 elements, which includes 461 beam elements and 158 truss elements. Six degrees of freedom (DOF) are available for each element node, and these DOFs correspond to translational and rotational displacements in the X, Y, and Z axes.
The input parameters of the material (Young’s modulus, specific gravity) as well as of the section (area, moment of inertia) are referenced from the as-built records. Specifically: Young’s modulus of steel (beams, truss rods, cantilever) E
s = 200 Mpa; the density of steel
ρs
= 7850kg/m
3. For non-structural elements such as bridge deck, balustrades, lighting systems, and plumbing are included in the model as additional mass. Typical cross-sections of some truss members are shown in
Table 2.
Table 2.
Cross-sectional of truss members.
Table 2.
Cross-sectional of truss members.
No |
Truss members |
Area (mm2) |
Moment of inertia Iy (mm4) |
Moment of inertia Iz (mm4) |
1 |
Bridge gate frame |
4.27104
|
2.9109
|
1.15109
|
2 |
Top lateral bracing |
4.75104
|
3.31109
|
1.75109
|
3 |
Bottom lateral bracing |
4.75104
|
3.31109
|
1.75109
|
4 |
Struts |
1.83104
|
1.03109
|
5.29107
|
5 |
Diagonal chords |
4.17104
|
2.82109
|
1.04109
|
6 |
Vertical chords |
1.83104
|
1.03109
|
5.29107
|
7 |
Top chords |
1.83104
|
1.03109
|
5.29107
|
8 |
Bottom chords |
1.83104
|
1.03109
|
5.29107
|
Boundary Conditions: The Dirichlet boundary conditions of the numerical model are created to accurately reflect the boundary conditions of the actual structure. The span of 8 bridges in Chuong Duong includes 2 types of bearings (
Figure 8). Based on the survey results, the displacement constraints of the model are made corresponding to the actual displacement capacity of the bearing.
Utilizing the Block-Lancios method, The FE model’s dynamic analysis is carried out.
Figure 10 and
Table 3 displays some mode forms’ mode shapes and natural frequencies.
The Modal Assurance Criterion (MAC), a statistical indicator, is particularly sensitive to significant discrepancies in the mode shapes and is comparatively insensitive to smaller differences [
47]. This results in a reliable statistical indication and consistency between the mode shapes. This study uses the MAC value to evaluate and update the finite element model. The formula 26 is used to calculate the MAC value:
Where n is number of modes shape considered, MAC is modal assurance criterion, ϕk, 𝜑̃k is modes shape FEM and experimental, ωk ; ω;̃k is frequencies FEM and experimental; T represents transposed matrix.
The MAC values are determined through formula (26) between the FE model results and the actual measurements. The first 3 MAC values greater than 0.85 show good agreement between each pair of mode shapes. However, other MAC values do not reach this minimum value. The correlation between the calculated and measured mode shape vectors are not guaranteed. The frequency values are also significantly different. This is a common situation for initial FE models, most of which have not been able to extract modes with high accuracy. Meanwhile, depending on the calculation requirements of the structure, some structures need high accuracy to structure health monitoring, diagnose damage, and make major predictions accuracy for maintenance. There are many uncertain parameters such as material properties, stiffness parameters. For this reason, it is recommended to perform a model update procedure to reduce errors.
The particle swarm optimization (PSO) algorithm was established and developed on the ideas of swarm intelligence to find solutions for optimization problems in a particular search space [
48,
49]. To understand the PSO algorithm better, observe a simple example of a flock of birds foraging. The foraging space is now the entire three-dimensional space. At the beginning of the search, the whole flock flies in a certain direction, it can be very random. However, after a while of searching, some individuals in the herd began to find a place to contain food. Depending on the amount of food just searched, the individual sends a signal to other individuals searching in the vicinity. This signal propagates throughout the population. Based on the information received, each individual will adjust its flight direction and speed in the direction of where there is the most food. Such communication is often viewed as a phenotype of herd intelligence. This mechanism helps the whole flock of birds to find out where there is the most food in the extremely large search space.
In swarm optimization, each particle searches a space by itself, remembering the best value, and informing other individuals. Other instances will receive the information and decide to continue the search or report its location so that other instances continue to act. So that, values in search space will be done quickly and accurately. There are two parameters that are particularly important, the location of an instance and the search velocity. These two parameters are expressed through the formulas for updating the position and updating the velocity of the instance:
Where xi is the position of instance i at different times (t and t+1); vi is the speed of individual i; w is the parameter of inertial weight; C1 and C2 represent the population’s cognitive coefficient; r1 and r2 are random numbers in the range 0,1]; pi(t) is the best position of each individual; Gbest is the best location of the entire population. Each individual is characterized by its velocity vector and its position in space.
To evaluate the similarity between the FE model and practical structure, an objective function is built based on the natural frequency and the mode shape of the structure:
In the case of Chuong Duong bridge, determining material parameters requires many experiments. At the same time, the masses of non-structural parts are also difficult to determine accurately. Some uncertain parameters are selected to update the numerical model. Based on experience
The updated uncertainty parameters and FE model after using the PSO are presented in
Table 5.
The results after updating show that updating the model lowers the discrepancy between the calculated and measured natural frequencies, and the MAC value reaches a good level.
The MAC values of the FE model before and after the update are shown in
Figure 9. After being updated, the accuracy of the model has increased a lot. These MAC values show good agreement between the FE and the actual models.
Figure 10 shows that before updating the model, the calculation results and natural frequency measurement results have high errors, after updating the model, the error of the natural frequency between calculation and measurement decreases significantly, the model is highly accurate and reliable.
The calculation results show that, after about 30 iterations performed by PSO, the model parameters begin to converge and give good results, minimizing the process of testing parameters in the modeling. After updating, the numerical model is accurate and almost the same as the actual object. This model is employed for creating data to train the ANN.
3.4. Generate Data and Train the ANN Model
a. Single damages
Three layers make up the ANN’s architecture in this work: an input layer, a hidden layer, and an output layer. The first seven modes’ frequencies are used as input data with various damage situations, and the output data comprises damage locations and levels (
Figure 11).
The updated model’s modal analysis creates input and output data for the ANN. The elements’ stiffness is decreased to create scenarios of structural damage. With a 1% interval, the elements’ stiffness decreases from 0% to 50%. Damages are only truly dangerous and sensitive enough when they occur on major structural components. Therefore, only the main truss rods are considered in this case. The equation determines the quantity of input data needed to train the network:
Where: ne – total of elements is considered, ns - the number of damage scenarios occurring for a single element. Consider 84 main elements of the truss. The input data will consist of 4200 samples
The ANN is configured to the settings for training the network after the data generation. The results of the network training process are significantly influenced by the number of neurons in the hidden layer. If the ANN has too few hidden layers, the ANN is too simple and difficult to deal with the problems to be solved. In contrast, ANN has too many hidden layers, the network is too complicated. Computer resources consume a lot, easily leading to overfitting. The loop is used to choose the ideal number of hidden layers. Within the range of 1 to 50, the number of hidden layers that will be selected by the loop. Additionally, the impact of noise is evaluated in all situations with a level of 2% for natural frequencies. The Levenberg-Marquardt backpropagation algorithm is used by the ANN to train the network. For damage identification, data split in the training procedure by 70%-15%-15% is used. There are a maximum of 1000 epochs. In the case of epochs greater than 1000 but still not reaching the best value, it is necessary to implement network optimization solutions.
After training, the following figures demonstrate how the ANN model performed:
All training cases with regression values more than 0.99 are displayed in
Figure 12. The training, evaluation, and test datasets are located along the target line (45degree line). This demonstrates that the real value and predicted value are almost similar. The regression values (R) in linear regression models always range from 0 to 1. The estimated and desired outcomes are the same if R is near to the upper bound (1).
Figure 13 shows the histogram of the calculated and intended output errors. There is extremely little variance between the target and the output.
Figure 14 shows the training performance in the datasets. The best validation performance value of 3.0179 is achieved at epoch 190. Epoch number does not exceed 1000, network optimization is not necessary. The graph also shows that the values in the data sets are quite convergent, no overfitting phenomenon occurs. From the performance graph of R-values, MSE (Tolerance), error histogram, the model has clearly been successfully trained, which can be used to apply to the contruction
b. 2 damaged elements
Damages are simultaneously generated by two random elements. For each element, a damage level is assigned from 0% to 50% with a 1% interval. In this study, the failures at the elements are assumed to be the same. The amount of data is calculated according to the formula:
The total amount of data generated is: 174300 samples
Figure 15 shows that
– values of the network using ANN is 0.986. The training, validation, and test datasets follow the regression line.
According to
Figure 16 and Figure 19, there are small deviations from the zero error line between computed and desired values. Results obtained show that a good agreement between predicted and actual outputs.
3.5. The Service of the Trained ANN Model in Actual
a. Single damage
Normally, to detect and locate the damage of the bridge structure, a comprehensive survey investigation will be carried out. After the survey, experts will make a computational model and evaluate. This work requires a lot of human resources and costs a lot of money. Only consider measuring the vibration of the entire span structure, the implementation of the setups measurements as described above also requires a lot of work. Instead, using an artificial neural network offers a huge advantage.
Using a trained network applied to the actual span structure of Chuong Duong bridge. During the service, there was an accident on Chuong Duong bridge that affected the bridge structure (
Figure 18).
A simple experiment was conducted since only the structure’s natural frequency needs to be ascertained. A single vibration measurement point in 3 directions is located on Chuong Duong bridge (
Figure 18). With only 1 measuring point, the natural frequency of Chuong Duong bridge has been identified. Although it is not possible to identify the mode shapes of the bridge, with the trained ANN it is possible to locate and quantify the damage.
The results of the nature frequencies of the first 7 mode shapes determined through the experiment are [1.229; 2.4167; 4.1954; 4.2623; 4.5925; 6.2147; 7.7865]. This result is used as input to feed the trained network. Results after putting data into the network return [25; 6]. This means that 25% damage is detected at element 6 (truss number 6). It is easy to see that the result of locating damage at element 6 is correct. According to the report of the UCT company who performed the calculation of damage assessment, damage level of truss rod is 30%. Results using ANN are relatively accurate. This is a case of damage that is easy to identify and recognize in practice. However, for cases where the damage location is difficult to detect (for example, the damaged location is located in too high, difficult to reach locations, in the middle of the river), the ANN will show its advantages.
b. 2 damaged elements
In this case, a hypothetical failure is generated on the updated finite element model. Two elements 10 and 15 are assumed to be 40% damaged. Put the data into the trained ANN model, the network detects and returns the result 39% and occurs at element 10 and 15. Although not 100% accurate, the trained network detected the location and was relatively close to the simulation results.