Refined Vehicle-Bridge Interaction Analysis Using Incompatible Solid Finite Element for Evaluating Stresses and Impact Factors

+e vehicle-bridge interaction can induce bridge vibration and consequently fatigue, durability deterioration, local damage, and even collapse of bridge structure. In this paper, a solid vehicle-bridge interaction (VBI) analysis method is developed to provide refined analysis on the bridge responses including displacement and local stress under vehicle loads. +e incompatible solid finite element (FE) is introduced tomodel the bridge, where the element shear locking is alleviated by incompatible displacement modes without sacrificing the computational efficiency. Benchmark example shows the incompatible solid element has superior computational efficiency compared to the conventional solid element. By virtue of the mass-spring-damper vehicle model, the interaction between vehicle and bridge is simulated with point-to-point contact assumption and the coupled dynamic equations are solved via nonlinear iteration. A case study on a simply supported T-girder bridge is conducted to validate the developed solid VBI analysis method and then the dynamic impact factor (DIF) of the bridge is evaluated based on the computed stress results and compared to code values. Results show that the solid VBI analysis method yields more accurate time-history bridge responses including displacement and stress under moving vehicles than the grillage method despite higher computational cost. Particularly, it can simulate realistic stress distribution and concentration along any concerned sections as well as in local components, which can provide detail information on the bridge behavior under dynamic loads. On the other hand, the DIF based on the computed stress result generally agrees well with the code values except for heavy vehicles where the stress-based DIF is slightly higher than the value in Chinese code while lower than that of AASHTO, suggesting the value specified by Chinese code may underestimate the DIF of heavy vehicles in certain circumstances to which more attention should be paid.


Introduction
e in-service bridges are subjected to moving vehicles every day and the dynamic interactions between vehicles and bridges can, on one hand, induce the bridge vibration and consequently the fatigue, durability deterioration, local damage, and even collapse of the bridge structure and on the other hand, lead to the influence on the vehicle safety as well as passengers' comfort. As a result, the vehicle-bridge interaction (VBI) has been extensively studied in the past few decades. It is well accepted that VBI is influenced by many factors including the bridge structure (e.g., span length and natural frequency), road surface roughness, vehicle suspension system, vehicle weight and vehicle speed, etc. Due to the complexity of VBI problem and the high cost of field test including difficulties in data collection, numerical VBI analyses have been widely employed to investigate bridge behaviors under dynamic vehicle loads [1][2][3][4].
Most numerical VBI analyses are based on the finite element method (FEM), in which the bridges are modeled by different types of element such as beam element, plate/shell element and solid element. In the beam-based bridge models [3,4], the bridge is simplified as a beam axis or a grillage with cross-sectional properties assigned to each beam element.
ey have the superiorities of simplicity and high computational efficiency [5]. Nevertheless, as the beam element yields the nodal displacement and the internal force on the cross-section that the node correspondingly represents, stress distributions on the cross-section cannot be obtained straightforward.
us, these models are usually preferred particularly in long-span vehicle-bridge analysis to analyze the overall performance of bridge structure. In plate/shell bridge models [6,7], thin plate/shell theories neglecting the out-of-plane normal and shear strains are assumed; thus, the complete stress-state of the modeled structure cannot be pictured. Meanwhile, as plate/shell elements approximate the bridge deck by its mid-surface, other element types are needed for non-plate components such as the girders and piers, which induces additional constraints on the structure interface [8,9]. Recently, with the advances in computer technologies, sophisticated bridge models using three-dimensional (3D) solid finite elements have been developed and applied into VBI analyses [10][11][12]. e solid element has three translational degrees of freedom (DOFs) per node and can be used to model structures of any geometry. With 3D kinematic, constitutive and equilibrium equations, it is possible to completely describe the spatial behaviors of the structure and produce stress distributions. Besides, as solid element does not involve any rotational DOF like beam and plate/ shell elements, cumbersome updating of rotational DOFs and convergence difficulties in nonlinear analyses are circumvented [13]. Most of the solid bridge models in VBI analyses are built in commercial finite element (FE) software, e.g., ANSYS and ABAQUS. e modal superposition technique, approximating the displacement by linear combination of vibration modes, is used to calculate the bridge response which is then inputted into VBI analysis [12]. is can significantly improve the computational efficiency; however, the simulation accuracy may be limited since only those lower-order modes are mainly combined. Besides, it is also worth mentioning that the modal synthesis method cannot be directly extended to nonlinear analysis of bridges. On the other hand, Kwasniewski et al. [10] developed a solid VBI model in which a solid bridge FE model, a detailed truck FE model and their couple vibration analysis are all implemented in LS-DYNA. Nevertheless, when more types of vehicles are simultaneously travelling on the bridge, the number of DOFs in the model would increase dramatically, causing excessive computational cost.
In this paper, a refined VBI analysis method is developed. e incompatible eight-node hexahedral (H8) element free of shear locking is used to model the bridge so that fewer elements are required along the bridge thickness to save the computational cost. e vehicle is represented as a multiaxle mass-spring-damper system which can model different types of vehicles with different number of axles. e pointto-point contact is assumed at the bridge-vehicle interface and the VBI dynamic equations are solved through a nonlinear iteration. e proposed solid VBI analysis method is expected to accomplish the transient displacement and stress analysis of bridges under moving vehicles. e paper is organized as follows: firstly, a solid VBI model is established, including the solid bridge FE model, the vehicle model, the vehicle-bridge coupling and the solution method; then, a case study of a simply supported T-girder bridge is conducted to validate the VBI analysis method; finally, it is applied to investigate the dynamic impact factors (DIFs) of the bridge under random traffic flow.

Solid Finite Element Model.
As lower-order elements are normally more inclined in engineering analyses due to their low construction cost [14], the trilinear H8 element portrayed in Figure 1 is selected to model the bridge structure. e mapping between the global Cartesian coordinates (x, y, z) and the natural coordinates (ξ, η, ζ) is expressed as in which (x i , y i , z i ) and (ξ i , η i , ζ i ) are, respectively, the global coordinate and the natural coordinate of the i-th node and N i (ξ, η, ζ) � (1 + ξ i ξ) (1 + η i η) (1 + ζ i ζ)/8 is the three-dimensional trilinear interpolation function. e isoparametric interpolation for the displacement of H8 element is as usual in which L is the strain-displacement operator; B = LN is the strain-displacement matrix. For homogenous, isotropic, and linear elastic materials, the stress-strain relation fulfills the generalized Hooke's Law, namely, in which C denotes the constitutive matrix. Based on the principle of minimum potential energy, the element stiffness matrix of the conventional H8 element is written as in which Ω and Ω e denote the problem domain and the element domain, respectively. is conventional H8 element is denoted as CH8. It is a well-known concern when the first-order elements are applied to model the structure with dominating bending deformation, undesired locking phenomena would be induced and result in excessively rigid response of the structure [15]. Numerous measures such as the mesh refinement, higher-order elements, reduced integration, nonconforming FE formulation and hybrid FE formulation have been proposed to eliminate the locking deficiency [16]. Among them, the non-conforming FE formulation in which the incompatible displacement modes are assumed by introducing internal DOFs is well accepted as it can considerably improve the element accuracy neither causing element instability nor sacrificing computational efficiency. For H8 element, the incompatible displacement modes proposed by Wilson et al. [17] are as follows: 2 Advances in Civil Engineering in which λ i denotes the assumed internal DOFs in the element and N and λ are self-defined. With equations (2) and (6), the element displacement is enhanced as and recalling equation (3) the element strain becomes accordingly, in which B � LN is the strain-incompatible displacement matrix. Following the usual treatment in [17], the element stiffness matrix of the incompatible H8 element is expressed as in which It should be mentioned that the above element stiffness matrix only corresponds to the compatible displacement DOFs, d, whilst the incompatible displacement parameters, λ, have been statically condensed out at the element level.
is incompatible H8 element is denoted as ICH8.
To show the accuracy and efficiency of ICH8 over CH8, a simply supported concrete plate with dimension 10 m × 10 m × 0.5 m is considered. e material parameters of the concrete are E = 3.0 × 10 4 MPa, v � 0.2 and ρ = 2 × 10 3 kg/m 3 . e central deflection and the normal stress at top center of the plate subjected to self-weight are predicted using CH8 and ICH8, respectively. In order to check whether the element suffers from shear locking, a single-layer mesh is designed along the thickness direction whilst different mesh densities of N × N elements are adopted for in-plane directions. e predicted central deflections and normal stresses are normalized by the thin-plate solutions [18] and plotted in Figure 2. It can be seen that ICH8 performs far better than CH8 for both deflection and stress predictions. For example, with the mesh density of 8 × 8, ICH8 produces the normalized deflection and stress about 0.92 and 0.96 whilst those of the latter are only about 0.43 and 0.59, respectively. ICH8 is validated to possess immunity to shear locking and superior accuracy for both deflection and stress predictions with coarser mesh compared to CH8. erefore, the ICH8 element is selected to model the bridge, which can save considerable computational cost without sacrificing accuracy.

Vehicle-Bridge Coupling.
e vehicle is represented as a multi-axle mass-spring-damper system with many DOFs depending on the number of axles. Taking the three-axle truck model in Figure 3 as an example [19], Z vr , Φ vr , and θ vr are, respectively, the vertical displacement, roll (rotation about the longitudinal axis or x-axis), and pitch (rotation about the transverse axis or y-axis) of the vehicle body, whilst Z i vaL and Z i vaR denote the vertical displacement of the left and the right wheel of i-th axle (i = 1, 2, 3), respectively. e equation of motion for the vehicle is expressed as where M v , C v and K v are, respectively, the mass, damping and stiffness matrices of the vehicle; u v is the vector of the vehicle's DOF; and F vg and F vb are the self-weight of the vehicle and the wheel-bridge contact force acting on the vehicle, respectively; as usual, the single over-dot and double over-dots denote, respectively, the first-and second-order time derivatives of the quantity below them.
For the bridge FE model using ICH8 elements, the equation of motion is expressed as    Figure 3: ree-axle vehicle model. and V(t) is the vehicle velocity. It should be noted that the measured road roughness condition which will be introduced later is employed in the present study. In equation (13), Z j bL(R) is the vertical displacement of the bridge at the contact point C and determined through the interpolation relation of ICH8 as in which L n (ξ, η) � (1 + ξ n ξ)(1 + η n η)/4 and Z k bn are the interpolation function and vertical displacement of the n-th node at contact surface of element k, respectively; can be obtained through similar interpolation. e vertical contact force exerted by the left (right) wheel of the j-th axle of the vehicle is distributed to the bridge deck through the nodal forces at the contact surface of element k as in which F k bvn is the vertical contact force distributed to the n-th node at contact surface of element k. Similarly, the vertical contact force by all axles of all vehicles is obtained and assembled to get the global wheel-bridge contact force applied on the bridge, F bv .
To solve the coupled vehicle-bridge model indicated in equations (11)-(15), the bridge and the vehicle equations are established independently as in equation (11) and equation (12), and then solved separately through a nonlinear iterative procedure in each time step according to the displacement compatibility and force equilibrium conditions at the vehicle-bridge contact interface. e main procedures for solving the VBI model are as follows: (i) Input basic data of the VBI model: at the time step, t−1, yielding the vertical stimulus sources at the current time step, t, of each tire for all vehicles; see equation (13). (iii) e forces acting on each vehicle at the current time step t induced by the vertical stimulus are calculated. (iv) e vehicle subsystem in equation (11) is independently solved via the Newmark integration method to obtain the initial vehicle response (1) e wheel-bridge contact force at the current time step, t, can be determined according to the obtained vehicle responses. (2) e load vector of the bridge is formed and the bridge subsystem in equation (12) is independently solved by the Newmark integration method and the current deformed position and internal forces are then obtained and updated.
(vi) e deck movements obtained in the step (v) are adopted for step (ii). (vii) e steps from (ii) to (vi) are repeated until the geometric compatibility and force equilibrium conditions at the wheel-bridge contact points are satisfied and then the calculations continue for the next time step, t + 1.
It is noteworthy that a grillage VBI model has been developed by the coauthor following similar procedure except that the bridge is modeled by the grillage method with beam elements. e grillage model has been successfully applied to investigate the bridge behavior including displacement and internal force under dynamic vehicle load [3,20] and is selected here as a compared counterpart of the present solid VBI model.

Validation of Solid Vehicle-Bridge
Interaction Model e solid VBI model established in Section 2 is validated through a case study on a simply supported T-girder bridge. Figure 4: Transfer of variables at the tire-bridge deck contact interface.

Description of Prototype Bridge. A prefabricated simply
supported T-girder bridge with a span length of 30 m (PCT-30), which represents the majority of concrete bridges in China, is considered here. e mid-span cross-section of PCT-30 is depicted in Figure 5. e bridge consists of five identical prestressed T-girders with uniform spacing of 2.35 m. Besides the end diaphragms, the other three intermediate diaphragms are also used to connect the five girders.

Field
Test. Both static and dynamic tests were conducted on the bridge. Two dump trucks, i.e., Truck-A and Truck-B, were used for loading. e detailed information of two trucks is given in Figure 6(a), where the axle weights were obtained in advance from a weighing system. e static test included two loading scenarios: (1) symmetrical loading of two dump trucks; (2) eccentric loading of two dump trucks, as illustrated in Figure 6(b). e longitudinal truck position in both scenarios was determined to yield the maximum displacement at the mid-span. e dynamic test considered the pass of Truck-A along the bridge center line at a uniform speed of 30 km/h, as shown in Figure 6(c), where GVW denotes the gross vehicle weight.
e bridge was equipped with displacement, strain and acceleration transducers to collect the bridge responses at the selected points that were expected to be significant. Displacements were measured independently at the five points located on the mid-span bottom surfaces of five girders, as shown in Figure 6(b). Strains were also measured at the same points as those of displacements, and stresses were obtained from those strain readings using the onedimensional linear stress-strain relationship. Acceleration histories were collected for the bridge using six accelerometers symmetrically instrumented on the L/4, L/2, and 3 L/4 (L denotes the span length) as shown in Figure 6(c). Readings of acceleration histories at the free vibration stage when the Truck-A had left the bridge were used for modal analysis of the bridge.

Bridge FE Model.
Based on the configuration of the bridge, a numerical model is established and shown in Figure 7(a). e structural components including main girders, diaphragms and barriers are all modeled using the incompatible solid element, ICH8. e prestressed tendons are not independently modeled whereas the prestressing effect is taken into consideration by modifying the elastic modulus of the concrete girders in the subsequent model updating process [21]. e cast-in-place joint is treated as part of the flange plate of T-girder due to its near-perfect connection stiffness. As ICH8 element is restrained from shear locking, only two layers of elements are distributed along the thickness direction of the bridge deck. Totally 7236 solid elements and 32797 nodes are involved in the bridge FE model. To better demonstrate the validity as well as superiority of the developed solid VBI model, the grillage model of the bridge is also built using 280 two-node beam elements as shown in Figure 7(b). Design values of concrete material properties in both bridge models are provided in Table 1.

Model Updating.
e deviation of FE model to field measured result is mainly attributed to various simplifications and approximations on the material parameters, the physical dimensions of the structure and the boundary conditions in the FE model [22,23]. Here, the boundary condition and the concrete material parameters including elastic modulus and density in solid and grillage bridge models are corrected according to the field measured results. e three-axle vehicle model in Figure 3 is used to simulate the two dump trucks given in Figure 6(a) and the mechanical properties which have been provided in [19] are not given here due to space limitation.

Boundary Conditions.
It is indicated in [9] that the design boundary conditions of bridge are sometimes not fully realized due to the actual situation of bearing and expansion joint. In order to determine the correct boundary conditions for FE models, eigenvalue analyses are conducted for four configurations of boundary condition listed in Table 2, where Case 1 refers to the design boundary condition. Axes x, y, and z denote the longitudinal, transverse and vertical directions of the bridge, respectively. U i and R i are the translation and rotation along the i-axis, respectively. Since there are only three translational DOFs per node for solid element, Case 1 and Case 2 are considered for the solid FE model while all the four cases are considered for the grillage model. e computed first vertical bending mode, first torsional mode and second vertical bending mode under each case of boundary condition are compared with the corresponding measured frequencies in Table 3. It is found that, for both grillage and solid model, the design boundary condition shows the best agreement with the actual condition.

Material Parameters of Concrete.
Because of uncertainty, the elastic modulus and density of girders and diaphragms, i.e., E G , ρ G , E D and ρ D , are selected as updating variables, among which modifying elastic modulus of concrete can take into consideration the prestressing effect [21]. e four parameters take initial values as the corresponding design values given in Table 1 and then updated according to the bridge responses obtained from both static and dynamic test, with the objective function: }denotes the modal assurance criteria [6] between the jth measured mode φ t and the j-th computed mode φ c ; χ i , α j and β j are weight factors to define the importance level of each residual term and are all taken as unity; n d = 5 is the total number of measured positions for displacement. e genetic algorithm [6] is implemented to resolve a set of parameters that minimize the objective function in equation (16) within a prescribed allowable tolerance, i.e., |f k static − f k−1 static | ≤ 10 − 5 . For reasonable updating results, it is prescribed that the change percent for E G and E D is within −5% and +20% whilst the absolute change percent for ρ G and ρ D is within 2%. e measured and computed vertical displacements before and after updating for two static loading scenarios are   Advances in Civil Engineering 7 compared in Figure 8. It can be seen that the computed vertical displacements from the updated bridge FE models agree well with the measured results for both loading scenarios. Meanwhile, the natural frequencies and correlation of mode shapes are also compared in Table 4. e differences between FE and measured frequencies are reduced below 1.61% and 3.45% for solid and grillage models, respectively, and the correlations of mode shapes are also well improved, higher than 92% and 90% for solid and grillage models, respectively, which is a significant improvement for the subsequent dynamic analysis of vehicle-bridge interaction. Besides, both static and dynamic results of the updated solid    bridge model are closer to the measured results than that of grillage counterpart, which indicates the solid element can approximate the real structure with higher accuracy than the beam element. e material parameters before and after updating are compared in Table 5. To check the physical meaning of the updated parameters, it can be found from Figure 8 and Table 4 that the computed deflections and frequencies from both initial FE models are larger than the corresponding measured values, indicating the initial FE models are less rigid than the real bridge. From Table 5, an increase in the value of stiffnessrelated parameters (elastic modulus) of main girders is accordingly observed while the values of mass-related parameters (density) are scarcely changed, as expected.

VBI Analysis of the Bridge.
With the updated bridge FE models, dynamic analysis on PCT-30 subjected to a Truck-A travelling along the bridge center line at a uniform speed of 30 km/h is conducted using the developed solid VBI model and the grillage VBI model. To better fulfill the real condition, the road surface roughness of the bridge has been measured and is depicted in Figure 9(a). e power spectral density (PSD) based on the measured road surface roughness is generated and compared to the three PSD curves corresponding to ISO levels, i.e., "Very good," "Good," and "Average," as shown in Figure 9(b). It is found that the overall road surface condition of PCT-30 is between "Good" and "Average." e computed time-history vertical displacements and bending stresses at the five points located in the mid-span bottom of the five T-girders are compared to the measured results from the dynamic test in Section 3.2. e stress of the solid VBI model is computed by nodal displacement via strain-displacement and stress-strain relations given in Section 2.1. e measured stress is obtained from the strain reading in dynamic field test using the one-dimensional linear stress-strain relationship. It is worth mentioning that the stress result of the grillage VBI model is calculated via the stress formula of pure bending theory because the bending moment on the cross-section containing node is the direct output from the grillage model. e result of the middle girder G3 is shown in Figure 10 as an example. It can be seen that the computed vertical displacement and bending stress generally agree well with the measurements except for larger fluctuations of the computed results at some time instants which might be induced by the deviation between the measured and the real road surface condition, data processing on the measured results, approximations on the dynamic characteristics of vehicle model, etc.
Furthermore, the bending stress at the top line of midspan cross-section when the front wheels of the truck is travelling to L/4, i.e., t � 0.9 s, is computed using both VBI     models and depicted in Figure 11(a). It can be seen that the bending stress obtained from the solid VBI model increases smoothly from the side girder to the middle girder, which shows reasonable variation as in the real situation. By contrast the stress obtained from the grillage model remains uniform in each girder and sharp stress change which should be avoided in engineering practice is obviously observed at the interface between girders, for the reason that the grillage model can only produce the overall internal force at the cross-section of each girder rather than the stress distribution. e stress distributions in local components are also simulated by the solid VBI model. Figure 11(b) shows the bending stress distribution in web plates when the front wheels of the truck are just travelling to the mid-span (L/2), i.e., t � 1.8 s. It can be seen that overall the stress increases from both ends to the mid-span; the value seems to be reduced at the position where the diaphragm is located to enhance the stiffness of girders; stress concentration is also observed at these locations.

Dynamic Impact Factor Spectrum of Bridge under Random Traffic Flow
In most of the existing researches, the DIF is traditionally calculated from the bridge responses in displacement or bending moment and some researchers also argued that the DIFs calculated using different bridge responses are different and thus should be treated differently [24,25]. To the authors' best knowledge, the DIF based on the load effect of stress is rarely reported. With the validated solid VBI model, the stress-based DIF of the prototype bridge, PCT-30, is investigated. e vehicle load characteristics of different vehicle category and the traffic loading scenarios are first obtained based on the monitored traffic data. e static and dynamic bending stresses at the mid-span bottom concrete of the girder G3 under each loading scenario are then computed by the solid VBI model with the measured road surface condition shown in Figure 9. e DIFs are calculated based on the computed static and dynamic bending stresses and the DIF spectrum of each vehicle category is investigated and compared to code values.

Traffic Data.
e traffic data of G104 highway in Shaoxing, China, was monitored using the weigh-in-motion (WIM) equipment and the vehicle load characteristics including the classification of vehicle type were also investigated in previous research of the authors [26]. Here, the WIM data recorded in May 13, 2013, i.e., 5338 sets of vehicle data, is employed to analyze the DIF spectrum of PCT-30 under random traffic flow. Figure 12(a) shows the traffic volumes of the 24 time blocks of May 13, 2013, in which the 1 st time block corresponds to the time interval from 0:00 am to 1:00 am and similarly the 24 th time block corresponds to the time interval from 11:00 pm to 0:00 am. Two traffic volume peaks located at the 10 th and the 17 th time block with the values 342 and 405, respectively, are observed in Figure 12(a). A local minimum value obviously appears in the 12 th time block for the potential reason that the corresponding time interval is just the lunch time. According to the previous research [26], the vehicle types of highway transportation can be merged into the following five categories: Category-I: passenger cars, pickups, etc. Category-II: medium-scale buses and medium-scale 2axle trucks Category-III: Large-scale buses and large-scale 2-axle trucks Category -IV: single unit 3-axle and 4-axle trucks Category-V: single trailer trucks Figure 12(b) shows the percentages of the five categories of vehicles. It can be seen that the vehicle quantity of Category-I is the highest, occupying 56.2%, whereas the quantity of Category-II is the lowest, occupying only 7.6%.

Random Traffic Flow Based on WIM Data.
e traffic loading scenarios are extracted according to the recorded velocity and arriving time of vehicles as well as the span length of the prototype bridge, in which all vehicles simultaneously loading on the prototype bridge are merged into one individual loading scenario [26]. Totally 4926 loading scenarios are generated, including 4514 single-vehicle loading scenarios and 412 double-vehicle loading scenarios. e detailed information of each loading scenario including the vehicle weight, velocity, and lateral position is determined according to the recorded WIM data.

Dynamic Impact Factor Spectrum.
e simulations of 4926 traffic loading scenarios are conducted using the developed solid VBI model to investigate the static and dynamic stress responses at the mid-span bottom concrete of the girder G3 with the measured road surface condition shown in Figure 9. en, the DIF for each traffic loading scenario is calculated according to the computed stress results.

Dynamic Impact Factor Spectrum of the Five Vehicle
Categories.
e DIF samples with respect to each vehicle category are obtained and the fitting analysis results are depicted in Figure 13. e samples of Category-II and Category-III are merged in the analysis because the former is relatively small. As indicated in Figure 13, the DIF distributions of the five vehicle categories can be described by extreme-I distribution or normal distribution. e 95 percentile points of the DIF distribution, i.e., DIF 0.95 , are also provided in Figure 13. e DIF 0.95 of Category-I is the highest with the value of 0.448, whilst the DIF 0.95 of Category-V is the lowest and its value is 0.207.

Dynamic Impact Factor Spectrum of Heavy Vehicles.
In China, the gross vehicle weight of 20 t is frequently taken to define heavy vehicles. Here, the traffic loading scenarios with a vehicle heavier than 20 t are selected and the sample Advances in Civil Engineering analysis of corresponding DIF is depicted in Figure 14. It can be seen that the DIF of heavy vehicles can be described by extreme-I distribution and the resulting DIF 0.95 is 0.279.

Comparison with the Code Values.
e DIF 0.95 of the five vehicle categories and the heavy vehicles are compared to the values from AASHTO and Chinese codes  [ 27,28], as shown in Figure 15. Note that those code values are calculated based on the natural frequency or span length of PCT-30. It can be seen from Figure 15 that the DIF 0.95 decreases for the vehicles from Category-I to Category-V for the underlying reason that the velocity usually decreases whilst the gross weight usually increases for vehicles from Category-I to Category-V. Moreover, the DIF 0.95 of these vehicle categories generally conforms well to the code values except for that of Category-I. Besides, the DIF 0.95 of heavy vehicles (GVW > 20 t) is slightly higher than the value in Chinese code while lower than that of AASHTO, suggesting the value specified by Chinese code may underestimate the DIF of heavy vehicles in certain circumstances to which more attention is paid.

Conclusion
In this paper, a solid VBI model is developed to provide full transient analysis on bridge responses including displacement and stress under random traffic loads. e incompatible solid finite element possesses higher computational efficiency than the conventional counterpart and thus is selected to model the bridge. rough the case study on a prototype bridge, PCT-30, the solid bridge model has better approximation on the real bridge and more accurate timehistory bridge responses including displacement and stress to moving vehicles are achieved directly from the solid VBI model, compared to the commonly used grillage VBI model. Specifically, the solid VBI model can simulate realistic stress distribution and concentration along any concerned sections as well as in local components, which can provide detailed information on the bridge structure behavior subjected to dynamic loads. e DIF of PCT-30 based on the stress results obtained from the validated solid VBI model is also investigated. According to the comparison with code values, it is found that the DIF based on the computed stress results generally agrees well with the code values except for heavy vehicles (GVW > 20 t) where the stress-based DIF is slightly higher than the value in Chinese code while lower than that of AASHTO, suggesting the value specified by Chinese code may underestimate the DIF of heavy vehicles in certain circumstances to which more attention is paid.
It is further expected to introduce the reinforcement and pre-stressed steel models to the present solid VBI model so that the local performance of reinforcement and pre-stressed steel as well as their effects on the static and dynamic responses of the bridge structure can be better observed. e solid VBI model is expected to provide a more refined tool to investigate the performance of small-to-medium-span bridges under dynamic service loads.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.