Abstract
The bridge-vehicle interaction (BVI) system vibration is caused by the vehicles passing through the bridge. The road roughness has a great impact on the system vibration. In this regard, poor road roughness is known to affect the comfort of the vehicle crossing the bridge and aggravate the fatigue damage of the bridge. Road roughness is usually regarded as a random process in numerical calculation. To fully consider the influence of road roughness randomness on the response of the BVI system, a random BVI model was established. Thereafter, the random process of road roughness was expressed by Karhunen–Loeve expansion (KLE), after which the moment method was used to calculate the maximum probability value of the BVI system response. The proposed method has higher accuracy and efficiency than the Monte Carlo simulation (MCS) calculation method. Subsequently, the influences of vehicle speed, roughness grade, and bridge span on the impact factor (IMF) were analyzed. The results show that the road roughness grade has a greater impact on the bridge IMF than the bridge span and vehicle speed.
1. Introduction
The bridge-vehicle interaction (BVI) vibration caused by vehicles passing through the bridge will have a direct impact on the working state and service life of the bridge. It is an important design measure to evaluate the structural design parameters such as the stability and safety of vehicles running on the bridge [1, 2]. The main causes of the BVI vibration are vehicle dynamic characteristics, road roughness, vehicle braking, vehicle speed, dynamic characteristics of the bridge’s structure, bridgehead overlap, bridge deck local potholes, vehicle driving position, etc.
Road roughness is an important factor affecting the BVI vibration. The influence of road roughness, which has strong randomness, cannot be ignored when calculating the BVI vibration. Road roughness is the main excitation source of the BVI vibration. Generally, road roughness is considered a stationary Gaussian random process with zero-mean value wherein roughness can be characterized by power spectrum density (PSD). Oliva et al. [3] proposed an algorithm for generating road roughness samples with correlation based on Fourier transform. On the basis of this method, Zou et al. [4] studied the influence of the sample correlation coefficient of roughness on the vibration response of the BVI system, pointing out that there are obvious differences in the calculated responses of the BVI system under different road roughness.
At present, there are many achievements in the study of BVI vibration. Zou et al. [5] proposed a method for predicting the bounds of the BVI responses with uncertain bridge and vehicle parameters, whereas Wang et al. [6] estimated the road roughness based on the BVI system. Similarly, Zhong et al. [7] analyzed the BVI dynamic response considering the foundation settlement, while Liu et al. [8] analyzed the safety of vehicles running on the bridge during earthquakes. Xu et al. [9, 10] analyzed the vehicle-track coupled system based on the multi-finite-element method. Xin et al. [11] analyzed the resonance of bridge considering the randomness of parameters, and Xin et al. [12] examined the uncertainty and sensitivity of parameters of the vehicle-bridge system. Meanwhile, Xiang et al. [13] analyzed the vehicle-bridge coupled dynamic system after considering the creep effect.
Generally, the response of the BVI system in consideration of road roughness can be calculated from the aspect of the random dynamic theory. In this context, Wu and Law [14] proposed a new method for dynamic analysis of the BVI considering uncertainty, in which the Karhunen–Loeve expansion (KLE) was applied to deal with the random field of road roughness. In addition, the statistical moment of system response was solved by the spectral method; the results showed that the method has high accuracy and efficiency. Then, Wu and Law [15] used a similar method to calculate the dynamic response of the BVI that considered the stochastic field of bridge structure parameters [15], vehicle axle load identification [16], and random system response calculation [17]. On the other hand, Chen et al. [18, 19] and Li et al. [20] proposed an efficient method to calculate the static and random response of the structural system. Xu et al. [21–25] carried out extensive research on the vehicle-bridge coupling vibration considering track irregularity. In the studies conducted by Liu et al. [26–28], KLE was also used to represent the stochastic process and combined with the point estimate method (PEM) [29, 30]. Thereafter, a new stochastic FEM approach and the dynamics response of vehicle and bridge were analyzed.
In order to study the dynamic response of vehicles crossing the bridge that fully considers the randomness of the road roughness, the moment method is applied to calculate the BVI system, in which KLE is used to mathematically express the random process of road roughness, whereas the PEM is used to calculate the statistical moment of the response. The system response can be evaluated quickly and easily by combining KLE and PEM methods (called KLE-PEM).
2. Dynamic Equation of the Bridge-Vehicle System
When a vehicle crosses a bridge, the two interact. Therefore, in the dynamic response analysis of the vehicle passing the bridge, the two are usually regarded as a system for analysis, which is referred to as the BVI system [31].
As shown in Figure 1, a car moves at a speed on a simply supported beam. The model of the vehicle is regarded as a mass-spring-dashpot system, which includes one car body, two wheel-axles, and two wheels. The car body has two degrees of freedom (DOFs), namely, vertical and nod displacement. Each wheel-axle has only one DOF. The wheel is supposedly in close contact with the road, which means that the wheel has no independent DOFs. Therefore, the whole vehicle system has four DOFs.

The FEM is used to model the bridge, in which the bridge damping matrix adopts Rayleigh damping. According to the energy law and considering the damping force of the system, the total potential energy of the BVI system can be obtained, which is expressed as follows:where Ui is the strain energy of the BVI system; Vm refers to the negative work done by the inertia force of the system and ; Vc denotes the negative work done by damping force and ; VF is the negative work done by Coulomb friction of the system and ; Vp signifies the negative work done by external force of the system and ; and Vg is the system gravitational potential energy and .
The dynamic equation of the system can be obtained quickly by combining the first-order variation of the total potential energy and the principle of elastic potential energy invariance [32]. The dynamic equation of the BVI can be written as follows:where M, C, and K represent mass matrix, damping matrix, and stiffness matrix, respectively; , and X denote acceleration vector, velocity vector, and displacement vector, respectively; the subscript vv or denotes vehicle; bb or b signifies bridge; and vb or bv is the coupling part of bridge and vehicle. F denotes the force vector, whereas Fv signifies the force vector of the vehicle caused by the displacement and the first derivative of road roughness, which can be written aswith and , where and denote the roughness of front and rear axle positions at time t, and and denote the first derivative of the roughness of front and rear axle positions at time t. FB denotes the force vector of the bridge caused by road roughness and axle load of the vehicle, which can be expressed aswithwhere denotes the shape function of the bridge at the position of front or rear axle, which can be written aswhere is the distance between the axle and the left node of the element and le signifies the length of the element.
Equation (2) will be changed by the movement of the vehicle. Therefore, the system is a time-varying system, whereas the Newmark-β integral method can be used to calculate the BVI system response.
3. Dynamic Statistics Calculation Approach
The Karhunen–Loeve expansion (KLE) is applied as an expression method to simulate the random process of road roughness. The basic expression can be written as follows:where denotes the mean value of road roughness process; and refer to the ith eigenvalue and eigenfunction of the covariance kernel, respectively; and denotes the ith random variable, which is a set of uncorrelated random variables and can be expressed as follows:
Because the road roughness is assumed to be a zero-mean Gaussian random process in this paper, will be a set of uncorrelated random normal variables. Equation (7) can be transformed as
After obtaining enough samples of road roughness, we can use the method introduced in [26] to solve and numerically.
According to the theory of statistical moment calculation, the first two central moments of the random system can be calculated by the following formula:where MEAN and VAR denote mean value and variance value, respectively.
Because the integral in (10) is continuous, the result of MEAN and VAR can be obtained by the Gauss integral. By utilizing the dimension reduction technique, the problem of calculating the statistical moment of a system with multiple random variables can be transformed into the problem of calculating the statistical moment of a composite system with a single random variable [27]. Thereafter, the calculation of (10) can be converted to the calculation of the following formula:withwhere r denotes the number of quadrature points, whereas xGH,l and represent the abscissa and weight of Gaussian–Hermite (G-H) quadrature, respectively. The detailed values of G-H are listed in Table 1. A more detailed derivation process can be found in [27, 34].
By combining the KLE and PEM (called KLE-PEM), the response of the BVI system can be easily and quickly obtained.
4. Simulation of Roughness
Road roughness spectrum recommended in Chinese code GB7031-2005 [35] was used as the roughness power spectral density (PSD) function to simulate the roughness samples, which can be written as follows:where n0 is the spatial reference frequency and the value is 0.1 m−1; denotes the road roughness PSD at the spatial frequency of n0, and its value is related to the road roughness class; is the frequency index, which determines the frequency structure of road roughness spectrum, wherein generally = 2; and n represents a spatial frequency in the effective frequency band of spatial frequency and signifies the number of cycles of waves contained in each meter of length, with its bandwidth being (n1, n2). Meanwhile, n1 and n2 are the upper and lower limits of the effective frequency band, respectively. The bandwidth should ensure that the vehicle vibration kinetic energy caused by road roughness includes the main natural frequency of vehicle vibration when the vehicle is driving at the average speed. In the code GB7031-2005, according to the road PSD, the road surface is divided into eight classes in accordance with the roughness. In this paper, five roughness grades were taken for calculation, as shown in Table 2.
Using the trigonometric series method (TSM) [16], the PSD function can be transformed into spatial samples of road roughness. Notably, 20,000 road roughness samples corresponding to the PSD of class A, class B, class C, class D, and class E were simulated with the length of 120 m (which can meet the total length of vehicles before entering, during passing, and after leaving the bridge). The roughness samples of different road roughness classes are shown in Figure 2. It can be seen that the roughness range of class A was about −3 to 3 mm. Similarly, the roughness range of class B was about -7 to 7 mm, class C was about −10 to 10 mm, class D was about −20 to 20 mm, and class E was about −50 to 50 mm.

(a)

(b)

(c)

(d)

(e)
The eigenvalues and eigenfunctions of the road roughness samples under different road classes were calculated to obtain the KLE mathematical expression of the random process of road roughness. The first 160 items of the characteristic value of the roughness under each road class are illustrated in Figure 3. Through calculation, the first 115 items of KLE can meet the accuracy requirement of 98% under different road roughness conditions. Therefore, in the subsequent calculations, 115 was used as the number of truncation terms for each road roughness class; that is, the BVI system has 115 independent random variables.

(a)

(b)

(c)

(d)

(e)
5. Numerical Simulation
5.1. Verification of the BVI Model
In order to verify the accuracy of the BVI model, the results of the classical numerical case in [36] were used to verify the model. In this case, a moving mass-spring system passing through a simply supported beam was calculated. The parameters of the BVI model were set to be consistent with the parameters of the case, subsequent to which the results of the vertical displacement time history response of the midspan are compared with the results of the case, as shown in Figure 4. It can be seen that the results of the current BVI model were essentially consistent with the results of the classic case, thereby verifying the accuracy of the BVI model.

5.2. Comparison with MCS
To verify the applicability of KLE-PEM in the calculation of BVI model, an example was taken in which the bridge was a one-span simply supported beam. Its parameters are shown in Table 3, and the first four mode shapes are shown in Figure 5. The vehicle was a two-axle vehicle, and its parameters are displayed in Table 4. The speed of the vehicle was 20 m/s.

(a)

(b)

(c)

(d)
MCS can be used to calculate the response of a stochastic system as the exact solution to evaluate the accuracy of the new method in a stochastic system. Here, MCS and KLE-PEM were used to calculate the BVI system with class C road. The number of times BVI programs were called by MCS was 20,000. Since the number of KLE truncation terms of a random process of road irregularity was 115 times, the number of the BVI programs called by KLE-PEM with three estimation points was 115 × 2 + 1 = 231. The comparison of the mean and variance of the bridge midspan displacement response obtained by calculation is shown in Figure 6. As evidenced by the figure, the results obtained by KLE-PEM coincide with the findings obtained by MCS in terms of mean or variance. In addition, the mean value of response of bridge displacement increased gradually with the passing of vehicles and then decreased gradually after the vehicles left the bridge, finally floating along zero. The variance of bridge displacement response increased when the vehicle moved forward and decreased when the vehicle moved forward. The time when the maximum value appears was essentially consistent with the time when the maximum absolute value of the mean displacement appeared.

(a)

(b)
The statistical moment results of the vertical acceleration time history response of the car body obtained by the two methods are demonstrated in Figure 7. The statistical moment obtained by KLE-PEM was essentially consistent with that obtained by the MCS method. In addition, for the mean value of response of the vehicle body, the mean value was zero before the vehicle entered the bridge. This is because the road roughness was assumed as a zero-mean random process.

(a)

(b)
5.3. Bridge Dynamic Response
The dynamic responses of the vehicle passing through the bridge at different vehicle speeds under different road grades were calculated, and the response’s maximum probability value was determined by the triple standard deviation criterion [27]. The maximum probability value of vertical displacement response in the middle span of the bridge was obtained, as shown in Figure 8. It can be seen that the bridge displacement response did not change significantly with the vehicle speed and that the road roughness class had a greater impact on the bridge displacement response. The bridge displacement response increased with the increase of road roughness grade. For class A road, the bridge displacement response was less than 6 mm, while the bridge displacement response was 10 mm for class D.

The maximum probability value of vertical acceleration response in the middle of the bridge span is shown in Figure 9. Unlike the displacement response, the acceleration response in the middle of the bridge evidently changed with the change of the vehicle speed, whereas the acceleration response first increased and then decreased with the change of the vehicle speed. For bridges with different roughness grades, the maximum acceleration response appeared when the vehicle speed is 20 m/s. In addition, the acceleration response of the bridge increased with the increase of road roughness grade.

5.4. Vehicle Dynamic Response
The maximum value of the vertical acceleration of the car body calculated by KLE-PEM under different roughness classes without vehicle speed is shown in Figure 10. It can be seen that in the case of different roughness levels, the acceleration of the car changes with the speed. More specifically, it first shows a trend of a decline followed by an increase. When the speed was 20 m/s, the maximum acceleration of the car body was found to be the smallest. Similarly, the roughness grade of a road surface was found to have a great impact on the acceleration of the car body.

5.5. Impact Factor of the Bridge
Dynamic impact factor (IMF) is generally defined as the ratio of dynamic load effect and static load effect of the bridge under vehicle load; i.e.,where Sd is the maximum dynamic response of the bridge under moving vehicle load, SS denotes the maximum static response of the bridge under corresponding vehicle load, and μ signifies the IMF.
IMF is not only an important index to characterize the impact effect of moving vehicle load on the bridge, but also a parameter affected by multiple factors. In traditional research studies, the dynamic IMF has been calculated by collecting the bridge response from the real bridge test, whereas the empirical calculation formula of IMF was obtained by regression analysis of the collected samples. However, on the one hand, the test method is expensive and difficult to implement; on the other hand, different bridges have different stress conditions, and the measured results of a limited number of bridges may not be widely representative. In contrast, IMF analysis based on numerical simulation has the characteristics of low cost and flexible simulation of different bridges and different working conditions, which is why it is eliciting widespread attention.
In Section 5.3, it can be inferred that road roughness has a great effect on the dynamic response of the bridge, thereby implying that road roughness influences the IMF of the bridge as well. To systematically discuss the IMF of the bridge, the IMFs with different vehicle speeds, different bridge spans, and different road grades were calculated from the random aspect. The range of vehicle speed was from 10 to 35 m/s; the bridge spans were 20 m, 22.5 m, 25 m, 27.5 m, and 30 m, whereas their first natural frequencies were 8.79, 6.94, 5.62, 4.65, and 3.90 Hz, respectively.
The calculated IMF under different cases is shown in Figure 11. For different bridge spans and road roughness grades, the IMF of the bridge generally increased with the increase of vehicle speed. For different bridge spans, the change of IMF was not particularly obvious, while the influence of road roughness on IMF was strong. For class A and class B roughness, the IMF was smaller than 1, while the IMF value was larger than 5 for class E. With the growth of the bridge operation time, it was found that the pavement roughness gradually increases, which is why the IMF of the bridge will also gradually increase over time, thus jeopardizing the bridge’s structural safety. Therefore, when maintaining the bridge, measurement and maintenance of the pavement roughness should be given due attention.

(a)

(b)

(c)

(d)

(e)
6. Conclusions
To comprehensively analyze the influence of road roughness on the response of the BVI system, a random BVI model considering the randomness of road roughness was established. The vehicle model was simulated by a mass-spring-damping system, and the bridge was simulated by the FEM theory. The time-varying system equations of the two were obtained by the energy variational principle. Karhunen–Loeve expansion was used to express the road roughness, and the response’s statistical moment was calculated by the PEM. Subsequently, the response of the BVI system with different speeds, different spans, and different roughness was analyzed. The results are as follows:(1)The KLE-PEM method can obtain the random response of the BVI quickly and precisely, and the computational efficiency is two orders of magnitude higher than that of MCS method.(2)The sensitivity of bridge displacement and acceleration response to the road roughness is higher than the vehicle speed. In addition, the coarser the road surface, the more obvious the influence of vehicle speed change on the bridge response.(3)Vehicle speed has a great influence on vehicle acceleration, but it does not increase linearly with the increase of vehicle speed; on the other hand, road roughness has a great influence on vehicle acceleration.(4)Compared with different speeds, different bridge spans, and different pavement roughness, road roughness has a strong influence on the IMF of the bridge, which should elicit sufficient attention when carrying out bridge engineering maintenance.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by a grant from the National Natural Science Foundation of China (Grant no. 12072375).