Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Composition, Structure and Properties of Geopolymer Concrete Dispersedly Reinforced with Sisal Fiber
Previous Article in Journal
Critical Factors Driving Construction Project Performance in Integrated 5D Building Information Modeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Load Identification on Prefabricated Girder Bridges Based on a CNN and Dynamic Strain Data

1
Shanghai Urban Construction and Maintenance Management Co., Ltd., Shanghai 200023, China
2
School of Water Conservancy and Transportation, Zhengzhou University, Zhengzhou 450001, China
*
Authors to whom correspondence should be addressed.
Buildings 2024, 14(9), 2809; https://doi.org/10.3390/buildings14092809
Submission received: 12 July 2024 / Revised: 31 August 2024 / Accepted: 5 September 2024 / Published: 6 September 2024
(This article belongs to the Section Building Structures)

Abstract

:
The vehicle load on a bridge is a critical and dynamic variable. It adversely affects bridges, especially when overloading occurs. Bridges are prone to fatigue damage or collapse. Therefore, identifying the size and type of dynamic vehicle loads on bridges is critical for theoretical studies and practical applications, such as health monitoring, daily maintenance, safety assessment, and traffic planning. The paper proposes a method for identifying the dynamic load parameters based on a convolutional neural network (CNN) and dynamic strain data. The model is implemented in MATLAB. An initial finite-element model of a three-span precast beam bridge is established in the software ABAQUS and modified by combining the modal and experimental data derived from a segmental girder bridge. The dynamic strain response of the bridge under a moving vehicle load is simulated under different working conditions. The results are used as the training data of the CNN to identify the vehicle’s position, speed, and load on the bridge. The high prediction accuracy indicates the proposed model’s suitability for identifying the dynamic load parameters.

1. Introduction

Bridges are critical components of the national road network. Concrete bridges are the most widely used bridge type. However, advances in science and technology have resulted in the application of new construction technologies, such as precast beam bridges. They have many advantages over other bridge types, such as fast construction speed, wide spans, lower cost, and environmental protection. Therefore, precast beam bridges will be extensively used in future bridge construction [1,2]. The number of prefabricated girder bridges in China has increased, and the threat of collapses cannot be ignored [3,4]. Therefore, it is critical to determine vehicle parameters efficiently and accurately, especially dynamic vehicle loads, to perform safety assessments of bridges and obtain traffic statistics [5,6].
Structural health monitoring of long-span bridges is typically conducted by installing a dynamic weigh-in-motion (WIM) system [7,8]. It uses sensors and data processing instruments to calculate vehicle weight, speed, and other information as the vehicle passes through the measurement site. Research on dynamic weighing started only recently in China. Due to the high installation and maintenance costs of the WIM equipment and monitoring system and complex factors affecting the working performance of bridges, these systems are typically only installed on service bridges and rarely on conventional beam bridges [9]. Using a WIM system to obtain vehicle parameter information has the advantages of low cost, convenient use, and simultaneous monitoring of the bridge’s utilization [10]. Zhou et al. [11] proposed an integral time-domain method to discretize the time-domain integral equation of the load response. They replaced the unit impulse response function with an integral form to obtain the load inversion matrix and transformed the load identification problem into a linear equation problem. Chan et al. [12] performed a Fourier transform of the equation of motion represented by modal coordinates. The relationship between the structural response and the moving force was established, and the least squares method was used to calculate the force spectrum in the frequency range. The time history of the moving force was obtained by performing an inverse Fourier transform. Qian et al. [13] used the influence line to establish a relationship between the moving load and the bending moment. They proposed a method for identifying the moving load using the bending moment influence line. Law et al. [14] simplified the bridge as a simply supported Euler beam model and expressed the time-varying interaction force between the vehicle and the bridge as a step function in a short time interval, expressing the vibration equation in modal coordinates. The moving load identification equation was derived using the modal superposition in the time domain. Au et al. [15] used a multi-stage genetic algorithm to select the optimum parameters by minimizing the error between the measured acceleration and the acceleration reconstructed from the identification parameters. The variable search domain was reduced step by step to identify the moving vehicle parameters and vehicle loads on the beam.
Although the theory of moving load identification on bridges is relatively mature [16,17,18], several problems remain in the practical application due to the complex and dynamic interactions between vehicles and bridges, such as low identification accuracy and noise pollution. In addition, the mathematical model of the identified structure should be established before identifying the load. It is difficult to determine the exact parameters for most structures, complicating the establishment of the theoretical model because the parameters affect the accuracy of moving load identification.
With the rapid development of artificial intelligence and computer technology, some scholars have used artificial neural networks to identify moving loads [19,20]. Kromanis et al. [21] proposed a vector regression prediction model to predict the thermal response of the bridge based on distributed temperature measurements. Six years of temperature and strain monitoring data of a concrete bridge were used as samples for the model. The results demonstrated the model’s effectiveness, but large amounts of historical monitoring data were required. Kim et al. [22] used an artificial neural network and data from a bridge WIM (BWIM) system to extract the parameters of heavy vehicles from the strain time-domain data. Numerous vehicle driving tests were conducted on simply supported prestressed concrete girder bridges and cable-stayed bridges to verify the model’s feasibility. However, the identification error of individual factors was large. Li et al. [23] proposed a phased identification method based on a backpropagation (BP) neural network for identifying moving loads on bridges. This method accurately identified the vehicle speed and dynamic loads in real time. However, it does not consider the interaction between moving load parameters, which is unrealistic. Chen et al. [24] analyzed the basic principle and implementation steps of a neural network to identify moving loads. They used an improved particle swarm optimization algorithm to obtain the deviation and weight of the neural network, but its effectiveness was not verified by experiments. Yang et al. [25] used an artificial neural network to identify the load, position, speed, and wheelbase of vehicles crossing a simply supported beam bridge based on the dynamic strain response. The feasibility of the neural network method was verified by a model test but not by a case study.
In summary, artificial neural networks have strong logical reasoning ability and are suitable for inverse problems such as structural load identification. They are less expensive than the installation and maintenance of bridge health monitoring systems. However, neural networks may have difficulty converging, and many training samples are required. Moreover, current neural network recognition methods are mainly aimed at simple and continuous girder bridges, whereas few studies have been conducted on prefabricated girder bridges. In this context, finite-element analysis (FEA), as a powerful numerical computation tool, can effectively complement the methods of artificial neural networks. By conducting finite-element modeling and analysis of bridge structures, FEA can provide a wealth of simulated data for the neural networks, therefore overcoming the issue of insufficient real training samples.
Therefore, we propose a method for identifying moving load parameters of vehicles on prefabricated girder bridges using a convolutional neural network (CNN) implemented in MATLAB. An initial finite-element model of a three-span segmental assembled bridge is established using the software ABAQUS 2021 and modified using modal test data obtained from a case study. The dynamic strain response of the bridge under a moving vehicle load is simulated under different working conditions, and the results are utilized as the training data for the CNN to identify the position, speed, and load of vehicles crossing the segmental girder bridge.

2. Moving Load Parameter Identification Based on CNN

A CNN can have various structures [26] but consists of several convolution, pooling, fully connected, input, and output layers. The convolution layer is the core of the CNN. The architecture is shown in Figure 1 [27].
The objective is to convolve the input data using the convolution kernel and introduce the activation function to obtain the output feature map. This operation generally increases the number of extracted features. The pooling layer is another critical layer. Its function is to downsample the input feature map, reduce the dimensions of the feature map and the number of parameters to prevent overfitting and ensure the stability of the network while reducing the computational complexity. The convolution and pooling layers alternate in the network and are used for feature extraction. The fully connected layer performs classification, i.e., the features are integrated to enable the classifier to provide an output.

2.1. Identification Method of Moving Load Parameters

It is necessary to extract the dynamic strain response of the main girder at different positions, vehicle speeds, and load combinations to identify the moving load of segmental girder bridges using a CNN. The strain signal, which is sensitive to the load identification parameters, is selected as the input vector, and the identification parameters of the moving vehicle on the bridge are output to obtain the pre-processed data set. A sample library for neural network identification is created based on feature changes. An appropriate network structure and training method are selected. Known samples are used as the training set to train the network, and test samples are used to test the recognition accuracy of the network. This process is used to adjust the network’s nonlinear mapping relationship. The method for the moving load identification of prefabricated bridges based on a CNN and dynamic strain is illustrated in Figure 2.

2.2. Structural Analysis and Network Design

The strain response data of the measuring point under specific working conditions is used as the network input, and the CNN is implemented in MATLAB 2021 software. The moving load parameter identification problem is a regression prediction problem with multiple outputs. After repeated network adjustment, a CNN with 11 layers is obtained [28].
The CNN structure is shown in Figure 3. The CNN is trained with the strain response data for different moving load parameters and used to establish a mapping relationship between the strain response and the moving vehicle parameters [29] to identify and output unknown moving vehicle parameters. The details of the different layers of the CNN are as follows:
(1) Input Layer. This layer receives a two-dimensional numerical matrix composed of the strain response data of segmented assembled beam bridges under the action of moving loads. The size of the matrix is A1 × A2 with one channel, where A1 is the number of strain measurement points under a single working condition, and A2 represents the number of collected data for a single strain gauge corresponding to the output indicators.
(2) Convolution layer 1. The size of the convolution kernel is B1 × B2, the number is 32, the convolution step size is 1, and the Padding method is valid. B1 and B2 correspond to the length and width of the convolution kernel layer. The convolution layer is used with a batch normalization layer and the ReLU layer. The former normalizes the output data of the previous layer in batches and passes it to the next layer to accelerate the convergence speed during network training. Common activation functions are Sigmoid, Tanh, and ReLU. The ReLU nonlinear activation function is more suitable for backpropagation and parameter updating of the network and prevents gradient disappearance. Therefore, the ReLU activation layer was used.
(3) Pooling layer. Maximum pooling is used to extract the background features. The size of the pooling area is C1 × C2, and the pooling step size is 1.
(4) Convolution layer 2. The size of the convolution kernel is D1 × D2, the number is 64, the convolution step size is 1, and the convolution method is valid. The convolution layer is also used with a batch normalization layer and a ReLU layer.
(5) Dropout layer. Due to the small training data set, the CNN has a deep structure with many parameters. A random deactivation layer is used to prevent overfitting. This layer randomly sets the input elements to zero according to a predetermined probability, and the remaining elements remain unchanged. This method improves the network’s generalization ability. The random zero probability is 20%.
(6) Fully connected layer. The fully connected layer reduces the dimensions of the learned features, integrates the local feature information extracted from the convolution layer or the pooling layer, and assigns sample labels. The number of neurons in the fully connected layer in the last stage is equal to the number of output target indicators.
(7) Output layer. Since this paper performs regression prediction, this layer is a regression layer that outputs the regression prediction results.

2.3. Hyperparameter Settings

In a deep network, hyperparameters refer to parameters that must be determined before network training. They control the model’s structure, function, and efficiency. The purpose is to obtain the optimal mapping relationship. While referring to existing design codes [30], after repeated network adjustment, the training batch size is 25, the maximum number of training epochs is 500, the initial learning rate is 0.001, the learning rate decline factor is 0.1, the learning rate decline cycle is 300, and the data are disorganized after each training round. We employ the commonly used learning algorithm Adam [31], an adaptive learning rate optimizer, which combines the characteristics of stochastic gradient descent (SGD) and root mean square propagation (RMSProp). Although the Adam optimizer is more complex than the SGD optimizers, it automatically adjusts the learning rate, has high computational efficiency, and is not prone to falling into a locally optimal solution. Thus, it is widely used in training neural networks.

2.4. Network Evaluation Index

The coefficient of determination ( R 2 ) [32] is used to assess the CNN’s performance. It is calculated as follows:
R 2 = 1 i = 1 n ( y ^ i y i ) 2 i = 1 n ( y ¯ i y i ) 2
where n is the number of samples, y ^ i is the calculated value, y i is the true value, and y ¯ i is the average value of the true value.
The smaller the R 2 value, the worse the performance of the CNN and vice versa. In addition, aiming at the recognition effect of the network on the moving vehicle load position, the gap between the recognition position and the real position is measured. In addition, the correct recognition rate of the network is calculated. It is the proportion of samples with a relative error of less than or equal to 10% [33]. The relative output error of the network is the ratio of the distance between the output vector of the network and the ideal output vector to the norm of the ideal output vector:
R e = i = 1 n | | y ^ i y i | | / 1 n i = 1 n | | y i | |
where: n is the number of samples, y ^ i is the calculated value, and y i is the true value.

3. Simulation of Precast Beam Bridge

3.1. Target Bridge

Determining the input and output parameters is crucial in load identification using neural networks. A 3 × 43 m prefabricated segmental bridge was used as a case study. The bridge has a three-span continuous beam, with a 43 m length for each span. The layout diagram is shown in Figure 4. The prestressed concrete single-box single-compartment segmental girder bridge is adopted in this joint, and the mid-span cross-section of the main beam is shown in Figure 5. The numerical simulation model of the bridge was established in ABAQUS, and the moving load conditions were applied using the DLOAD subroutine to calculate the response data of the bridge’s control section under different working conditions. Based on specific vehicle parameters, internal functions such as time, travel speed V, and vehicle load F are written and adjusted in Fortran, enabling the representation of various load effects of different vehicles [34,35]. The influence of the initial displacement and strain caused by prestressing was eliminated in the initial data processing.

3.2. Numerical Simulation

A finite-element model of the entire segmental assembly bridge was established using ABAQUS software, with the main girder and prestressing tendons modeled using C3D8R and T3D2 elements, respectively [36]. The mesh utilized a structured partitioning type; however, due to the complexity of the finite-element model, a direct structured mesh could not be implemented. Therefore, the partition tool was used to segment the irregular areas of the model before proceeding with the mesh division, resulting in a more optimized mesh. Ultimately, the entire bridge consisted of 301,820 nodes and 204,635 elements. The finite-element mesh model of the main girder section is shown in Figure 6, and the lane layout is depicted in Figure 7.
The cohesive finite-element method with finite interface thickness was used to model the adhesive joints [37]. A cohesive element consisting of epoxy resin with a thickness of 1 mm was inserted between the prefabricated sections. The material parameters of the components are listed in Table 1, Table 2 and Table 3.
Embedded elements were used to simulate the interaction between the prestressed tendons and the concrete to ensure convergence and high computational efficiency [38]. The bottom of the No.0 cast-in-place block in the middle span was attached to the pier using a full constraint. The remaining bearings were one-way movable bearings to limit the displacement in the x and y directions and the rotation in the x and z directions [34]. The three-dimensional model of the bridge is shown in Figure 8. The unit length in the longitudinal direction of the bridge is 0.2 m

3.3. Sensor Layout

The model was validated by monitoring the vehicle load on the target bridge. The maximum positive bending moment and the maximum vertical displacement occurred in the main girder span of the bridge. The maximum negative bending moment and shear stress were generated near the main pier fulcrum. Therefore, three sections in the middle of the main beam span, the adjacent joints of the 1/4 beam segment, and the 3/4 beam segment were selected as monitoring and control sections. One-third points at the bottom of the beam in each control section were selected to measure the strain. The measurement point configuration is shown in Figure 6. The dynamic strain acquisition was carried out using the DHDAS strain acquisition system. Strain gauge sensors were used to measure the strain in the control section. The specific equipment is illustrated in Figure 9. The parameters of the sensors are listed in Table 4.

3.4. Model Validation

To process the strain data from the field test in a targeted manner, it is first necessary to do a statistical analysis of the traffic volume during the field test. The total traffic volume in one month is selected for statistical analysis, and the statistical results of the monthly average daily traffic volume are shown in Figure 10.
The daily maximum hourly traffic volume basically appears in the morning and evening peaks, while in the evening, the traffic volume is smaller, but through the monitoring records the general heavy vehicles can be seen basically appearing in the evening. Therefore, the vehicle load and speed information at night were screened, and it ultimately found that the load of 300 KN vehicles accounted for 42% of the total number of vehicles, vehicle speed distribution peak more in 40–100 km/h, the night due to fewer vehicles, so the speed will be relatively large, so in the verification of the results of the analysis of the vehicle conditions selected the interval of the driving parameter values.
A modal comparison analysis of the segmental girder bridge model was conducted. The first-order frequency of the bridge was 2.615 Hz through simulation calculation. The first-order frequency of the bridge was 2.637 Hz through measurement, and the ratio of the measured to the simulated value of the first-order natural frequency was 1.01. It shows that the finite-element model established in this paper can simulate the real structure of the bridge well. Its damping ratio is 0.03.
Based on the aforementioned actual operational statistical analysis of the bridge, two sets of control conditions were established. In the second lane, a 300 kN heavy vehicle was allowed to travel at uniform speeds of 15 m/s and 20 m/s without any obstruction. After applying the moving load, the maximum strain values were extracted from the central measurement points at the bottom of the cross-sections along the bridge at 1/4 L, 1/2 L, and 3/4 L. The fitted curves comparing the measured strains near the bottom edge of the beam slab corresponding to the traffic lanes with the simulation results obtained from ABAQUS are shown in Figure 11 and Table 5.
The measured and calculated values of the dynamic strain were in good agreement. A slight difference was observed when the vehicle reached the control section of the main beam under each working condition, indicating that the vehicle did not travel at the design average speed. However, the difference in the peak value of the dynamic strain between the measured and calculated results was only about 1 s, demonstrating a good agreement.
The relative error of the maximum strain for the 1/4 L, 1/2 L, and 3/4 L sections was within 5%, verifying the accuracy of the finite-element model. The maximum strain was smaller than the calculated value, indicating that the mechanical performance of the segmental assembled beam bridge was consistent with the theoretical results.

4. CNN Training Test

4.1. Analysis of Network Input Parameters

A sufficient number of strain data were obtained from the ABAQUS numerical simulation. Nonlinear mapping was performed from the input to the output parameters using the CNN. The significant eigenvalues of the bridge responses were selected as input variables to improve the accuracy and efficiency of network learning. Figure 12 and Figure 13 show the strain time history curves for the bottom center of the main beam’s control section and the four points at the bottom of the mid-span cross-section (the configuration of the measuring points is shown in Figure 6) as a vehicle with a standard axle load of 200 kN traveled on the third lane of the target bridge.
As shown in Figure 12, the vehicle load significantly affected the dynamic strain response of the segmental girder bridge. The regularity trend of the dynamic strain time history curve at 1/8 and 7/8 sections of the middle span is not as strong as that of other sections. The dynamic strain time history curves were similar in the 1/4, 3/4, and 1/2 sections. Therefore, only the strain time history curves in these sections were used to identify the vehicle parameters.
The vehicle position substantially influenced the strain response. As shown in Figure 13, the No. 4 measuring point was the closest to the lane’s center. The strain value of this measuring point was the highest in this section when the vehicle traveled in lane 3. Therefore, the strain time history curve of the measuring points in the same section must be considered in the load identification.

4.2. Dynamic Vehicle Parameter Identification

4.2.1. Dynamic Vehicle Load Category Identification

The strain time history data obtained from the 9 measuring points at the bottom of the 1/4, 1/2, and 3/4 sections of the mid-span of the bridge were used as input to the CNN. The sampling frequency was the ratio of the maximum design speed to the length of the two longitudinal cells; thus, 9 × 356 data points were extracted under one working condition. The travel time differed for different speeds. Therefore, the same number of data points were used to compare the results for different working conditions. The strain was 0 when there was no load on the bridge.
The vehicle load, driving speed, and lane were considered and simplified. According to the “Code for Design of Urban Bridges” [39], the plan dimensions of a standard heavy vehicle are shown in Figure 14. The factors and the design of the load conditions for a single vehicle are listed in Table 6, Table 7 and Table 8. A total of 240 samples were used (240 × 9 × 356); 70% of the samples were used as training data, and 30% were used as test data. The number of levels increased or decreased according to the test results and training conditions. The strain caused by the self-weight of the bridge and the prestress was ignored. The training data to obtain the optimal parameters of the neural network are listed in Table 9.

4.2.2. Identification of Moving Vehicle Position

The vehicle positions must be known to assess the strain responses at the measuring points. Therefore, it is necessary to change the features of the input layer sample dataset in Section 4.2.1. Since the distance between the front and rear axles of the same vehicle type is constant, the position of the front axle relative to the main span was determined, and the input data were filtered accordingly. The number of samples after reclassification was 30,240. One working condition was selected as the test data, and the other working conditions for the same vehicle type were used as the training data, resulting in 6 × 10 groups of working conditions. The number of training sets for each working condition was 378, the number of test sets was 126, and the ratio was 75%:25%. The training set was used to train the CNN, and the test set was used to assess the model’s performance to determine the vehicle’s position on the bridge.
After repeated debugging, the hidden layer of the model was used to fit the complex relationship between the feature matrix and the position information after the input signal was pooled by the convolution kernel. Its structure was similar to that of the neural network. The number of hidden layer neurons was 16, and the training batch size was changed to 30 in the hyperparameter setting of the deep learning model. The remaining settings were the same as described in Section 2.3 to obtain the optimal parameters of the CNN (Table 10). Increasing the number of convolutional or fully connected layers did not significantly increase the recognition accuracy, whereas reducing the number of layers decreased it.

5. Results and Discussion

5.1. Model Performance

The model results for the training and test sets are shown in Figure 15, Figure 16 and Figure 17. To facilitate subsequent analysis, we define the output variables—namely the lane in which the vehicle operates, the vehicle load type, and the vehicle speed—as index 1 to 3, respectively. The following is observed:
(1) The predicted lane position, vehicle type, and vehicle speed obtained from the CNN model were very close to the actual values. The fitting degree of the CNN model was high, indicating that it had high training performance.
(2) The results for the test set demonstrated the high generalization ability of the neural network. The difference between the predicted and actual values for the indices was very small, indicating that the neural network could accurately identify and predict the moving load parameters. In addition, the output parameters were suitable for characterizing the traffic on the bridge.
The R 2 values are listed in Table 11.
The following is observed in Table 11:
(1)
The R 2 values for the training set exceeded 0.9, indicating high model performance. Those for the test set differed but were greater than 0.9. The minimum R 2 value was 0.9068 for index 1, and the maximum was 0.9564 for index 2. The model had a high fit for the three moving load parameter indices, with the best performance for index 2.
(2)
The ranking of the indices based on the R 2 value in the test set was index 2 > index 3 > index 1. The CNN exhibited different fitting degrees for various types of moving loads, showing a better performance for vehicle load identification and a worse performance for lane identification.

5.2. Recognition Accuracy

Figure 18 shows the model’s recognition accuracy for identifying the front axle position of a class I vehicle passing through the middle span of the bridge at a speed of 10 m/s. Errors occurred when the vehicle entered and left the mid-span. When the vehicle was traveling on the bridge, the identified position was very close to the actual position.
Table 12 lists the correct recognition rates of the CNN for the vehicle’s position under different working conditions. The rates exceeded 80% [40], indicating a high accuracy in identifying the vehicle’s position. The lowest recognition rate was 80.27% for a two-axle vehicle traveling at 13 m/s, and the highest rate was 95.70% for a five-axle vehicle moving at 19 m/s.
The recognition rates for the vehicle speed and load type were similar, and good recognition results can be obtained for both small and large wheelbases. This is attributed to the permutation and method to select samples, resulting in uniform samples and high similarity in the vehicle conditions.
As the vehicle speed increased, the recognition rate of the vehicle’s position increased, reaching the maximum at 19 m/s, followed by a decrease in the recognition rate. After the vehicle reached the mid-span at a given speed, the recognition rate increased with the number of axles and the axle load of the vehicle. The reason for this is that the strain increases with an increase in the spacing between the front and rear axles [41]. The position recognition accuracy increased with the number of axles and the axle load and with the vehicle speed in a certain speed range. Subsequently, the recognition rate decreased slightly, but it met the engineering needs.

5.3. Effect of Noise on CNN Network Load Identification

During the data acquisition process, environmental factors and instrument noise will introduce certain levels of noise to the strain signals. Since these noise influences are challenging to eliminate, it is crucial to consider the impact of noise on the test signals. In this section, we first simulate actual measurement samples by adding 1%, 5%, and 10% random noise to simulated data [42]. Next, we train a CNN network following the training scheme in Section 4.2. Finally, 60 sets of representative strain signals were randomly selected under three different working conditions, namely the lane affected by vehicles, the type of vehicle load, and the vehicle running speed, to evaluate the recognition accuracy of the mobile vehicle parameters detected by the CNN network. The recognition prediction results for each category are shown in Figure 19 and Table 13.
The above results indicate that noise has a certain impact on the recognition of load parameters. As the level of noise increases, the recognition accuracy of the network gradually decreases. Among them, the addition of noise has the greatest impact on the identification of vehicle load categories, with a decrease in recognition accuracy of 20.19%. However, even when 10% of noise is added, the recognition accuracy of the network is still above 75%, indicating that the network remains effective under these conditions.

6. Conclusions

This paper establishes a single-vehicle moving load identification model for precast beam bridges based on convolutional neural networks (CNN) and dynamic strain data, performing regression predictions under various working conditions. The model demonstrated the highest performance in predicting vehicle load types, achieving an R2 value of 0.9564, while it showed the lowest performance in lane identification, with an R2 value of 0.9068. The recognition accuracy of vehicle position improved with the number of axles and axle loads. Within the speed range of 10 m/s to 19 m/s, the recognition accuracy increased with speed, reaching a maximum of 95.70% for a five-axle vehicle traveling at 19 m/s. Subsequently, the model’s accuracy slightly decreased but remained above 80%, meeting the engineering accuracy requirements.
This study only considered the effects of different vehicle speeds, weights, and lane positions of a single vehicle on the dynamic response of bridges. Future research should also take into account the impact of different vehicle fleets on dynamic performance. In summary, the segmental assembly beam bridge single-vehicle moving load parameter identification model based on CNN and dynamic strain exhibits good predictive performance, indicating that it is feasible to identify single-vehicle moving load parameters by measuring the strain response of segmental assembly beam bridges under vehicle traffic. This provides a technical reference for the operation, maintenance, and performance evaluation of segmental beam bridges.

Author Contributions

Investigation, S.Y.; Resources, X.Z.; Writing—original draft, W.W.; Writing—review & editing, P.G.; Visualization, L.L.; Project administration, Y.C.; Funding acquisition, L.Z. All authors have read and agreed to the published version of the manuscript.

Funding

The authors are grateful for the financial support from the Key Research Projects of Higher Education Institutions in Henan Province (24A560021 (Pan Guo)) and the Henan Postdoctoral Foundation (202102015 (Pan Guo)).

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

Authors Lun Zhao, Xuetao Zhang and Liang Li were employed by the company Shanghai Urban Construction and Maintenance Management Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Ahmed, G.H.; Aziz, O.Q. Shear strength of joints in precast posttensioned segmental bridges during 1959–2019, review and analysis. Structures 2019, 20, 527–542. [Google Scholar] [CrossRef]
  2. Shi, X.; Liu, Z.; Hu, K.; Zhou, Z. Full-scale Test of Bearing Capacity of a Complete External Prestressed Segmental Precast Continuous Girder Bridge. China J. Highw. Transp. 2018, 31, 163–173. [Google Scholar]
  3. Xu, F.Y.; Zhang, M.J.; Wang, L.; Zhang, J.R. Recent Highway Bridge Collapses in China: Review and Discussion. J. Perform. Constr. Facil. 2016, 30, 04016030. [Google Scholar] [CrossRef]
  4. Cook, W.; Barr, P.J. Observations and Trends among Collapsed Bridges in New York State. J. Perform. Constr. Facil. 2017, 31, 04017011. [Google Scholar] [CrossRef]
  5. Yu, L.; Chan, T.H.T. Recent research on identification of moving loads on bridges. J. Sound Vib. 2007, 305, 3–21. [Google Scholar] [CrossRef]
  6. Law, S.S.; Bu, J.Q.; Zhu, X.Q.; Chan, S.L. Moving load identification on a simply supported orthotropic plate. Int. J. Mech. Sci. 2007, 49, 1262–1275. [Google Scholar] [CrossRef]
  7. Richardson, J.; Jones, S. On the use of bridge weigh-in-motion for overweight truck enforcement. Int. J. Heavy Veh. Syst. 2014, 21, 83–104. [Google Scholar] [CrossRef]
  8. Haider, S.W.; Masud, M.M.; Selezneva, O.; Wolf, D.J. Assessment of Factors Affecting Measurement Accuracy for High-Quality Weigh-in-Motion Sites in the Long-Term Pavement Performance Database. Transp. Res. Rec. 2020, 2674, 269–284. [Google Scholar] [CrossRef]
  9. Li, S.L.; Liang, Z.Z.; Guo, P. A FBG pull-wire vertical displacement sensor for health monitoring of medium-small span bridges. Measurement 2023, 211, 112613. [Google Scholar] [CrossRef]
  10. Li, J.; Law, S.S. Damage identification of a target substructure with moving load excitation. Mech. Syst. Signal Process. 2012, 30, 78–90. [Google Scholar] [CrossRef]
  11. Zhou, H.C.; Li, H.N.; Yang, D.H.; Yi, T.H. Moving force identification of simply supported bridges through the integral time domain method. J. Sound Vib. 2022, 534, 117046. [Google Scholar] [CrossRef]
  12. Chan TH, T.; Yu, L.; Law, S.S. Comparative studies on moving force identification from bridge strains in laboratory. J. Sound Vib. 2000, 235, 87–104. [Google Scholar] [CrossRef]
  13. Qian, C.Z.; Chen, C.P.; Xiao, Y.G. Identification Method for Moving Loads over Continuous Beam Based on Bending Moment Influence Lines. In Proceedings of the 3rd International Conference on Civil, Architectural and Hydraulic Engineering (ICCAHE), Hangzhou, China, 30–31 July 2014. [Google Scholar]
  14. Law, S.S.; Chan TH, T.; Zeng, Q.H. Moving force identification—A frequency and time domains analysis. J. Dyn. Syst. Meas. Control-Trans. Asme 1999, 121, 394–401. [Google Scholar] [CrossRef]
  15. Jiang, R.J.; Au FT, K.; Cheung, Y.K. Identification of vehicles moving on continuous bridges with rough surface. J. Sound Vib. 2004, 274, 1045–1063. [Google Scholar] [CrossRef]
  16. Deng, L.; Cai, C.S. Identification of Dynamic Vehicular Axle Loads: Demonstration by a Field Study. J. Vib. Control 2011, 17, 183–195. [Google Scholar] [CrossRef]
  17. Deng, L.; Cai, C.S. Identification of Dynamic Vehicular Axle Loads: Theory and Simulations. J. Vib. Control 2010, 16, 2167–2194. [Google Scholar] [CrossRef]
  18. Zhu, X.; Cao, M.S.; Ostachowicz, W.; Xu, W. Damage Identification in Bridges by Processing Dynamic Responses to Moving Loads: Features and Evaluation. Sensors 2019, 19, 463. [Google Scholar] [CrossRef]
  19. Zhang, Y.; Hearn, G.E.; Sen, P. Neural network approaches to a class of ship control problems. I. Theoretical design. Elev. Ship Control Syst. Symp. 1997, 111, 115–133. [Google Scholar]
  20. Cui, C.; Sang, N.; Gao, C.X.; Zou, L. Vehicle Re-identification by Fusing Multiple Deep Neural Networks. In Proceedings of the 7th International Conference on Image Processing Theory, Tools and Applications (IPTA), Montreal, QC, Canada, 28 November–1 December 2017; IEEE: Piscataway, NJ, USA, 2017. [Google Scholar]
  21. Kromanis, R.; Kripakaran, P. Support vector regression for anomaly detection from measurement histories. Adv. Eng. Inform. 2013, 27, 486–495. [Google Scholar] [CrossRef]
  22. Kim, S.; Lee, J.; Park, M.S.; Jo, B.W. Vehicle Signal Analysis Using Artificial Neural Networks for a Bridge Weigh-in-Motion System. Sensors 2009, 9, 7943–7956. [Google Scholar] [CrossRef]
  23. Zhongxian, L.I.; Feng, C.; Bo, W. A BP Neural Network-Based Stage Identification Method for Moving Loads on Bridges. Eng. Mech. 2008, 25, 85–92. [Google Scholar]
  24. Chen, X.-H.; Sun, L.-L.; Wang, D. Identification of moving load on the bridge based on improved bacteria foraging optimization algorithm. J. Qingdao Technol. Univ. 2013, 34, 25–28, 66. [Google Scholar]
  25. Yang, H.; Yan, W.M.; He, H.X. Parameters Identification of Moving Load Using ANN and Dynamic Strain. Shock Vib. 2016, 2016, 8249851. [Google Scholar] [CrossRef]
  26. Fan, Z.; Lin, H.B.; Li, C.; Su, J.; Bruno, S.; Loprencipe, G. Use of Parallel ResNet for High-Performance Pavement Crack Detection and Measurement. Sustainability 2022, 14, 1825. [Google Scholar] [CrossRef]
  27. Jadeja, V.; Rao AL, N.; Srivastava, A.; Singh, S.; Chaturvedi, P.; Bhardwaj, G. Convolutional Neural Networks: A Comprehensive Review of Architectures and Application. In Proceedings of the 2023 6th International Conference on Contemporary Computing and Informatics (IC3I), Gautam Buddha Nagar, India, 14–16 September 2023; pp. 460–467. [Google Scholar]
  28. Zhou, Y.; Zhao, Y.; Hao, G.; Fang, L. Based on the strain signal time-frequency analysis with CNN the vehicle load identification method of the network. J. Hunan Univ. Nat. Sci. 2022, 49, 21–32. [Google Scholar]
  29. Xu, B.; Wu, Z.S.; Chen, G.D.; Yokoyania, K. Direct identification of structural parameters from dynamic responses with neural networks. Eng. Appl. Artif. Intell. 2004, 17, 931–943. [Google Scholar] [CrossRef]
  30. Gu, Y.; Lu, W.; Qin, L.; Li, M.; Shao, Z. Short-term prediction of lane-level traffic speeds: A fusion deep learning model. Transp. Res. Part C 2019, 106, 1–16. [Google Scholar] [CrossRef]
  31. Chen, T.Y.; Sun, Y.J.; Yin, W.T. Solving Stochastic Compositional Optimization is Nearly as Easy as Solving Stochastic Optimization. Ieee Trans. Signal Process. 2021, 69, 4937–4948. [Google Scholar] [CrossRef]
  32. Yin, W.D. Research on Community Risk Prediction Model and Management Based on Deep Learning. Wirel. Commun. Mob. Comput. 2022, 2022, 2373284. [Google Scholar] [CrossRef]
  33. Yang, H.; Li, F.; Yan, W.; He, H. Moving Load Identification and Experimental Verification of Beam Bridge Based on Dynamic Strain and ANN. J. Vib. Meas. Diagn. 2018, 38, 305–313. [Google Scholar]
  34. Han, Y.; Zhou, Z.; Zhang, Y.; Liu, X.; Wang, J.; Zhou, W. Dynamic Deflection Simulation of Bridge in Driving Process Based on ABAQUS Software. Technol. Earthq. Disaster Prev. 2023, 18, 118–126. [Google Scholar]
  35. Dong, Z.; Liu, M.; Zheng, H.; Gong, X. Dynamic Mechanical Analysis of Asphalt Pavement Based on Cross-isotropic Properties. China J. Highw. Transp. 2012, 25, 18–23. [Google Scholar]
  36. Zivaljic, N.; Nikolic, Z.; Smoljanovic, H.; Munjiza, A. Numerical simulation of reinforced concrete structures under impact loading. Mater. Und Werkst. 2019, 50, 599–610. [Google Scholar] [CrossRef]
  37. Sheng, Z.; Wang, Y.; Zhang, Z. Local cohesive finite element method with finite interface thickness. Chin. J. Comput. Mech. 2021, 38, 802–809. [Google Scholar]
  38. Li, F.; Zhang, W.; Ni, J. Study on seismic performance of reinforced concrete frame structure with embeded steel plates. Build. Struct. 2020, 50, 106–113. [Google Scholar]
  39. CJJ 11-2011; Code for Design of Urban Bridges. Shanghai Municipal Engineering Design & Research Institute: Shanghai, China, 2011; 109P.
  40. Wu, D.; Zhao, R. Research on the Method of Load Identification for Concrete Bridge Based on Neural Network. China Railw. Sci. 2002, 23, 25–28. [Google Scholar]
  41. Deng, L.; Luo, X.; Ling, T.; He, W. Bridge Weigh-in-motion Algorithm Considering Multi-vehicle Based on Convolutional Neural Network. J. Hunan Univ. Nat. Sci. 2022, 49, 33–41. [Google Scholar]
  42. Zhang, C. Research on Mobile Load Identification Method Based on BP Neural Network. Zhejiang University: Hangzhou, China, 2007. [Google Scholar]
Figure 1. Typical convolutional neural network structure diagram.
Figure 1. Typical convolutional neural network structure diagram.
Buildings 14 02809 g001
Figure 2. Flow diagram for the identification of moving load parameters.
Figure 2. Flow diagram for the identification of moving load parameters.
Buildings 14 02809 g002
Figure 3. CNN recognizes the network structure diagram and related information.
Figure 3. CNN recognizes the network structure diagram and related information.
Buildings 14 02809 g003
Figure 4. General layout (Unit: m).
Figure 4. General layout (Unit: m).
Buildings 14 02809 g004
Figure 5. Layout drawing of a cross-section of the main girder mid-span(Unit: mm).
Figure 5. Layout drawing of a cross-section of the main girder mid-span(Unit: mm).
Buildings 14 02809 g005
Figure 6. Finite-Element Mesh Model of Main Girder Segment.
Figure 6. Finite-Element Mesh Model of Main Girder Segment.
Buildings 14 02809 g006
Figure 7. Four-lane driving band diagram.
Figure 7. Four-lane driving band diagram.
Buildings 14 02809 g007
Figure 8. Finite-element model of segment-assembled bridge.
Figure 8. Finite-element model of segment-assembled bridge.
Buildings 14 02809 g008
Figure 9. DHDAS strain acquisition system equipment.
Figure 9. DHDAS strain acquisition system equipment.
Buildings 14 02809 g009
Figure 10. Monthly average daily hourly traffic volume.
Figure 10. Monthly average daily hourly traffic volume.
Buildings 14 02809 g010
Figure 11. The strain time history curve of the key section of the main beam at different speeds. (a) The strain time history curve of the key section of the main beam at a speed of 15 m/s. (b) The strain time history curve of the key section of the main beam at a speed of 20 m/s.
Figure 11. The strain time history curve of the key section of the main beam at different speeds. (a) The strain time history curve of the key section of the main beam at a speed of 15 m/s. (b) The strain time history curve of the key section of the main beam at a speed of 20 m/s.
Buildings 14 02809 g011
Figure 12. The strain response time history curve of each section of the main beam.
Figure 12. The strain response time history curve of each section of the main beam.
Buildings 14 02809 g012
Figure 13. The strain time history curve of each measuring point of the mid-span section of the main beam.
Figure 13. The strain time history curve of each measuring point of the mid-span section of the main beam.
Buildings 14 02809 g013
Figure 14. Standard dimensions of heavy vehicles (unit: m). (a). Plan dimensions of a dual-axle heavy vehicle. (b). Plan dimensions of a triple-axle heavy vehicle. (c). Plan dimensions of a quad-axle heavy vehicle. (d). Plan dimensions of a five-axle heavy vehicle.
Figure 14. Standard dimensions of heavy vehicles (unit: m). (a). Plan dimensions of a dual-axle heavy vehicle. (b). Plan dimensions of a triple-axle heavy vehicle. (c). Plan dimensions of a quad-axle heavy vehicle. (d). Plan dimensions of a five-axle heavy vehicle.
Buildings 14 02809 g014
Figure 15. Index 1 Recognition results of the training set and test set.
Figure 15. Index 1 Recognition results of the training set and test set.
Buildings 14 02809 g015
Figure 16. Index 2 Recognition results of the training set and test set.
Figure 16. Index 2 Recognition results of the training set and test set.
Buildings 14 02809 g016
Figure 17. Index 3 Recognition results of the training set and test set.
Figure 17. Index 3 Recognition results of the training set and test set.
Buildings 14 02809 g017
Figure 18. Vehicle position recognition results.
Figure 18. Vehicle position recognition results.
Buildings 14 02809 g018
Figure 19. Identification results of load parameters under different noise influences.
Figure 19. Identification results of load parameters under different noise influences.
Buildings 14 02809 g019aBuildings 14 02809 g019b
Table 1. Concrete material characteristics.
Table 1. Concrete material characteristics.
Material TypeModulus of Elasticity (MPa)Poisson’s RatioWeight Capacity (kN·m−3)Compressive Strength (MPa)Tensile Strength (MPa)
C603.6 × 1040.226.527.52.04
Table 2. Material properties of steel strand.
Table 2. Material properties of steel strand.
Material TypeModulus of Elasticity (MPa)Compressive Strength (MPa)Tension Control Stress (MPa)Coefficient of Thermal Expansion (°C−1)Weight Capacity (kN·m−3)
Prestressed bar1.95 × 1041860.01395.0 (1305.0)27.578.5
Table 3. Performance index of segmental assembly epoxy resin adhesive.
Table 3. Performance index of segmental assembly epoxy resin adhesive.
Colloidal propertiesTensile strength (MPa)40.0
Tensile modulus of elasticity (MPa)3200.0
Flexural Strength (MPa)50.0
Compressive strength (MPa)80.0
Elongation (%)1.5
bonding abilitySteel-Steel T Impact Stripping Capacity (mm)20.0
Standard value of steel-steel tensile shear strength (MPa)15.0
Positive tensile bond strength of steel as base material (MPa)15.0
Tensile bond strength to concrete (MPa)2.5
Table 4. Sensor parameter.
Table 4. Sensor parameter.
Model Number120-50AA
resistance120 ± 0.5 Ω
Sensitivity coefficient2.11 ± 1%
Strain limit20,000 μm/m
Operating temperature−20–80 °C
Table 5. The maximum strain error is calculated by field measurement and ABAQUS.
Table 5. The maximum strain error is calculated by field measurement and ABAQUS.
Machine Speed (m/s)Location of Cross-SectionMeasured ValueCalculated ValueRelative Error %
Time of Arrival (s) Maximum   Strain   Value   ( μ ε )Time of Arrival (s) Maximum   Strain   Value   ( μ ε )
151/4 L3.607.883.568.102.8
1/2 L4.2013.964.2914.272.2
3/4 L4.957.964.998.122.0
201/4 L2.558.142.678.372.8
1/2 L3.2014.903.2315.141.6
3/4 L3.758.233.748.422.3
Note: Relative error = |Maximum measured dynamic strain − Maximum simulated dynamic strain|/Maximum measured dynamic strain.
Table 6. Vehicle loading type.
Table 6. Vehicle loading type.
Vehicle TypeStandard Load (B1) (kN)Overload 50% (B2) (kN)Overload 100% (B3) (kN)
Two-axle vehicle (A1)200 (I)//
Tri-axle vehicle (A2)300 (II)450 (III)600 (IV)
Four-axle vehicle (A3)400 (V)600 (VI)800 (VII)
Five-axle vehicle (A4)500 (VIII)750 (IX)1000 (X)
Table 7. Single-lane traffic lanes.
Table 7. Single-lane traffic lanes.
Lane
Traffic lanesFirst laneSecond lane Third laneFourth lane
Lane type numberingC1C2C3C4
Table 8. Vehicle speed (m/s).
Table 8. Vehicle speed (m/s).
Vehicle Speed(m/s)
10 13 16 19 22 25
Speed type numberingD1D2D3D4D5D6
Table 9. Vehicle load category identification network structure parameter setting table.
Table 9. Vehicle load category identification network structure parameter setting table.
Layer Name Parameter Setting
Input layer9 × 356 × 1
Convolution layer 1(2 × 2 × 16) [2 × 2]
Convolution layer 2(2 × 2 × 32) [1 × 1]
Fully connected layer3
Output layer1 × 1 × 3
Note: In the convolutional layer, “16” in “(2 × 2 × 16) [2 × 2]” represents the number of convolutional nuclei; “(2 × 2)” represents the size of the convolution kernel; “[2 × 2]” represents the size of the core used for the maximum pooling operation. The parameter settings for the fully connected layer represent the number of neurons.
Table 10. Vehicle Position identification network structure parameter setting table.
Table 10. Vehicle Position identification network structure parameter setting table.
Layer NameParameter Setting
Input layer9 × 1 × 1
Convolution layer 1(3 × 1 × 16) [2 × 1]
Convolution layer 2(3 × 1 × 32) [1 × 1]
Fully connected layer 116
Fully connected layer 21
Output layer1 × 1 × 1
Note: As shown in Table 9.
Table 11. The prediction effect of each output index of the convolutional neural network ( R 2 value).
Table 11. The prediction effect of each output index of the convolutional neural network ( R 2 value).
Data SetIndex 1Index 2Index 3
Training sets0.94060.96320.9845
Testing set0.90680.95640.9433
Table 12. Qualified rate of vehicle load position identification under different vehicle types and speeds.
Table 12. Qualified rate of vehicle load position identification under different vehicle types and speeds.
Qualified Rate/%Speed (m/s)
101316192225
Vehicle loading typeI84.8380.2790.6491.7388.1883.65
II86.0985.9391.2391.2392.1186.76
III85.2189.4489.4789.4789.4786.76
IV85.2388.5790.3592.1190.3585.89
V89.9189.7988.3390.0192.590.62
VI89.2489.2994.1788.3385.8385.62
VII90.7489.7993.3394.1791.6788.96
VIII90.3191.8193.9493.3793.7989.09
IX91.8692.5894.5095.7095.0291.26
X90.3192.5894.1693.1794.5690.42
Table 13. The impact of different noise levels on identification effectiveness.
Table 13. The impact of different noise levels on identification effectiveness.
Pass Rate1% Noise Level5%1% Noise Level10%1% Noise Level
Vehicle lane94.16%87.55%81.71%
Vehicle loading type96.37%89.36%76.18%
Vehicle speed96.95%84.71%78.52%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, L.; Wu, W.; Zhang, X.; Li, L.; Guo, P.; Yang, S.; Cai, Y. Dynamic Load Identification on Prefabricated Girder Bridges Based on a CNN and Dynamic Strain Data. Buildings 2024, 14, 2809. https://doi.org/10.3390/buildings14092809

AMA Style

Zhao L, Wu W, Zhang X, Li L, Guo P, Yang S, Cai Y. Dynamic Load Identification on Prefabricated Girder Bridges Based on a CNN and Dynamic Strain Data. Buildings. 2024; 14(9):2809. https://doi.org/10.3390/buildings14092809

Chicago/Turabian Style

Zhao, Lun, Wenqi Wu, Xuetao Zhang, Liang Li, Pan Guo, Shaolin Yang, and Yingchun Cai. 2024. "Dynamic Load Identification on Prefabricated Girder Bridges Based on a CNN and Dynamic Strain Data" Buildings 14, no. 9: 2809. https://doi.org/10.3390/buildings14092809

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop