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 (
) [
32] is used to assess the CNN’s performance. It is calculated as follows:
where
n is the number of samples,
is the calculated value,
is the true value, and
is the average value of the true value.
The smaller the
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:
where:
n is the number of samples,
is the calculated value, and
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.
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.