Research Article  Open Access
Qingfei Gao, Zonglin Wang, Hongyu Jia, Chenguang Liu, Jun Li, Binqiang Guo, Junfei Zhong, "Dynamic Responses of Continuous Girder Bridges with Uniform CrossSection under Moving Vehicular Loads", Mathematical Problems in Engineering, vol. 2015, Article ID 951502, 29 pages, 2015. https://doi.org/10.1155/2015/951502
Dynamic Responses of Continuous Girder Bridges with Uniform CrossSection under Moving Vehicular Loads
Abstract
To address the drawback of traditional method of investigating dynamic responses of the continuous girder bridge with uniform crosssection under moving vehicular loads, the orthogonal experimental design method is proposed in this paper. Firstly, some empirical formulas of natural frequencies are obtained by theoretical derivation and numerical simulation. The effects of different parameters on dynamic responses of the vehiclebridge coupled vibration system are discussed using our own program. Finally, the orthogonal experimental design method is proposed for the dynamic responses analysis. The results show that the effects of factors on dynamic responses are dependent on both the selected position and the type of the responses. In addition, the interaction effects between different factors cannot be ignored. To efficiently reduce experimental runs, the conventional orthogonal design is divided into two phases. It has been proved that the proposed method of the orthogonal experimental design greatly reduces calculation cost, and it is efficient and rational enough to study multifactor problems. Furthermore, it provides a good way to obtain more rational empirical formulas of the DLA and other dynamic responses, which may be adopted in the codes of design and evaluation.
1. Introduction
The dynamic response of bridges due to moving vehicular loads has been a subject of interest to engineers more than 100 years. Owning to the uncertain, transient, and nonlinear characteristics of the vehiclebridge coupled vibration system, the phenomenon of dynamic response is one that requires fairly complex mathematical treatment to adequately characterize the motion and forces within the structure. However, for simplicity, it will instead be conceptual treatment using a common format, one that describes dynamic response as a fraction or multiple of the response that would be obtained if the same forces or loads were applied statically. This is, in fact, the same approach that most design specifications use to calculate the effects of dynamic live loading. When using such an approach, the amount of response above the static value is typically called the dynamic increment and is found by multiplying the static value by an “impact fraction.” Alternatively, the total response can be expressed as a multiple of the static value using an impact factor [1]. In recognition of the complex behavior associated with the dynamic response from vehicular loading, several authors have observed that the term “impact” is too limited and therefore not descriptive of the actual behavior [2–5]. Instead, the current trend is to replace the term “impact fraction” with “dynamic load allowance,” which represents the response from all types of vehicular dynamic effects, not solely impact. In this study, the term “dynamic load allowance” is adopted.
It has been proved that many different definitions of dynamic load allowance are adopted by various researchers and engineers [6, 7]. Based on the design theory of the bridge in Chinese current code, the dynamic load allowance in this paper is defined as follows:where is the location of the selected section, denotes the dynamic load allowance of the selected section of the bridge, and and , respectively, denote the maximum dynamic response and the maximum static response of the selected section. Furthermore, and mean the location of the vehicle when the dynamic response and the static response reach the maximum value. It should be noted that and are not the same in most of the time.
The dynamic response of a bridge to a crossing vehicle is a complex problem affected by the dynamic characteristics of both the bridge and the vehicle and by the bridge surface conditions [1]. Extensive research and development has been carried out to understand the vibration of bridges as a result of natural sources of vibration and to determine the dynamic allowance as a function of the fundamental frequency due to its uniqueness [8–11]. So it is important to find a reliable method to estimate the natural frequency of bridges. Komatsu et al. evaluated the natural vibration responses of straight and curved I or boxgirder bridges using Vlasov’s beam theory [12, 13]. Cantieni conducted an experimental study to establish a reliable relationship between the fundamental frequency and maximum span length of bridges [14]. Meanwhile, Billing presented tables of natural frequencies for straight, continuous symmetric multispan uniform beams of two through six spans, for various span ratios [15]. Additionally, in order to reduce the significant difference between the estimation of the natural frequency obtained from the codes and theoretical methods or numerical simulation, Gao et al. proposed a numerical improved method for straight bridges [16, 17]. Also, the empirical expressions are established by Mohseni et al. using regression analysis, to determine the fundamental frequencies of skew continuous multicell boxgirder bridges [18]. However, some limitations still exist for these methods to actually be used in practice.
Over the last 60 years, a significant amount of research has been conducted in the area of bridge dynamics, and this research has been both analytical and experimental in scope. Following the collapse of several railway bridges in Great Britain, the first study of vehiclebridge interaction was conducted in 1849 by Willis [19]. A major experimental investigation of bridge impact loads was initiated in 1958 by American Association of State Highway Officials (AASHO). Eighteen simply supported bridges traversed by the twoaxle trucks and the threeaxle tractortrailer combinations were investigated, each with a span of 15 m [20]. Wright and Green reported on a series of experimental investigations of vehiclebridge interaction made on 52 highway bridges in Ontario from 1956 to 1957 [21]. Csagoly et al. conducted a second series of tests on bridges (11 bridges) in Ontario in 1969 to 1971 [22]. Page [23] and Leonard [24] reported on a series of dynamic tests conducted on highway bridges (30 bridges) by the Transport and Road Research Laboratory (TRRL) in England. Shepard and Sidwell [25], Shepard and Aves [26], and Wood and Shepard [27] reported on fourteen bridge tests conducted in New Zealand. In 1980, a third series of dynamic testing of highway bridges in Ontario was performed. The tests were conducted on 27 bridges of various configurations of steel, concrete, and timber construction [5, 28]. Cantieni presented results from the dynamic load tests of 226 highway bridges in Switzerland from the mid1950s to the early 1980s [14, 29]. As for the numerical simulation, in 1970, Veletsos developed a numerical method for computing the dynamic response of highway bridges to moving vehicular loads [30]. Guant et al. investigated analytically the effects on bridge accelerations of several parameters, including the properties of the bridge and the vehicle as well as the initial conditions of the roadway [31]. Inbanathan and Wieland performed an analytical investigation of the dynamic response of a simply supported box girder bridge due to vehicle moving across the span [32]. Coussy et al. presented an analytical study of the influence of random surface irregularities on the dynamic response of bridges [33]. Recently, Bakht and Pinjarkar presented a literature review of bridge dynamics, with a particular focus on the dynamic testing of highway bridges [7]. And Paultre et al. presented a comprehensive review of analytical and experimental findings on bridge dynamics and the evaluation of the dynamic load allowance [19]. All these studies have proved that many of the parameters interact with one another, further complicating the issue, and consequently, many research studies have reported seemingly conflicting conclusions.
Therefore, the method of orthogonal experimental design in batches [34–38] has been proposed in this paper. It can be used for studying both the main effects and the interaction effects. Meanwhile, due to the processing in batches, it greatly reduces the calculation cost.
In this study, the natural frequencies of the simply supported girder bridge and the continuous girder bridge with uniform crosssection are firstly investigated based on the theoretical derivation and the numerical simulation. And then the effects of different parameters on the dynamic responses of the vehiclebridge coupled vibration system are discussed by our own program VBCVA. Finally, the orthogonal experimental design method is proposed for the dynamic responses analysis, including the case without interaction and the case considering interaction.
2. Natural Frequencies of Multispan Continuous Girder Bridges
In civil and transportation engineering, dynamic characteristics of structures always consist of natural frequencies, corresponding vibration modes, and damping ratios. The natural frequency is one of the most important dynamic characteristics, due to its key role in determining dynamic responses of the structure under dynamic loads, such as wind, earthquake, and moving vehicles. Also, the impact factor or dynamic load allowance, which is used for considering the dynamic increments generated by moving vehicular loads, is given as the function of natural frequencies in most codes and specifications for bridge design. In addition, the serviceability of bridges, especially for pedestrian bridges, is closely related to the natural frequency. Therefore, the natural frequency is firstly studied in this paper, and it is the basis of dynamic responses analysis.
2.1. Simply Supported Girder Bridges
The boundary of the simply supported girder bridge is assumed as pinned joint in both ends. However, according to inspections of existing bridges, we found that the ends could not completely rotate with the degeneration of rubber bearings. For convenience of engineers and researchers, three types of boundary are investigated, which are shown in Table 1.
 
Note: means the deflection is zero, means the rotation is zero, and means the moment is zero. 
As for the Euler beam with uniform crosssection, if the damping is ignored, the vibration governing equation can be simplified aswhere the denotes the force and , , , , , and denote, respectively, the mass per unit length, elasticity modulus, moment of inertia, displacement, location, and time.
Based on the method of separation of variables, the solution can be listed as follows:where is the th natural frequency (rad/s), is the th vibration mode, and is the constant parameter determined by both the crosssection and the boundary conditions. In addition, , , , and are constant parameters, which are determined by the boundary condition.
Submitting (4) into different types of boundaries in Table 1, then the results can be seen in Table 2.

It can be seen from Table 2 that the natural frequency is higher when there are more constraints and the corresponding stiffness is higher. The th natural frequency of the singlespan bridge with different boundary conditions can be summarized aswhere is used for denoting different boundary conditions and is the span length. When the boundary is hinged and hinged, the parameter . As for the other two conditions, and in turn.
For common bridges, the moment of inertia is positively correlated with the factor ( is the height of the beam), and the area is positively correlated with the hight . Therefore, the th natural frequency is negatively related to the span length . And it successfully explains the phenomenon that most empirical formulas regressed by fielding test for estimating natural frequencies are negatively correlated with the span length .
2.2. Continuous Girder Bridges with Uniform Spans
As for 2span, 3span, and 4span continuous girder bridges with uniform crosssection and uniform spans, the natural frequency and corresponding vibration modes are given by Gao et al. [16]. Based on a similar method, more general frequency equation can be listed as follows:where
For common bridges, the number of spans is not more than eight. Equation (6) is solved by the general software MATLAB, and the frequency equations and the natural frequencies are listed in Tables 3 and 4.


It can be seen from Table 4 that the first natural frequency of the multispan continuous girder bridge with uniform spans is completely the same as that of the corresponding simply supported girder bridge. Also, for bridges with different spans, the fundamental vibration modes are retained. However, the number of other vibration modes among fundamental vibration modes is more for the bridge with more spans.
In the Chinese code D602004 for bridge design, the impact factor is defined as the function of the natural frequency. For continuous girder bridges, the impact factor of the positive moment in midspan is different from that of the negative moment in piertop section, which is determined by the curvature of different modes. According to the analysis of vibration modes distribution, for positive moment in midspan section, the first mode is significant for all bridges. However, for negative moment in piertop section, the curvature of the second vertical vibration mode is larger for evenspan bridges, and the curvature of the third vertical vibration mode is larger for oddspan bridges (Figure 1).
So the second natural frequency in the evenspan bridge and the third natural frequency in the oddspan bridge are selected for calculating the impact factor of negative moment in piertop section. To make it clear, both of them are denoted as . Based on Table 4, the empirical formula of can be given bywhere is the modified factor, and it can be gained bywhere is the number of spans.
2.3. Continuous Girder Bridges with Nonuniform Spans
The multispan continuous girder bridge is usually seen in urban areas and other areas with special terrain. For urban bridges, when the road extends under the bridge, the span above the road will be larger than others. For bridges in special terrain areas, the piers cannot be constructed at the idealized location due to the bad terrain occasionally, so the spans should be nonuniformly arranged.
To our experience of bridge design, the length of side span is always smaller in the multispan continuous girder bridge with nonuniform spans. The span ratio is introduced here. And it has been assumed that only the length of side span is different from others (Figure 2):where is the length of side span and is the length of middle span.
The frequency equation of this type of bridge can be given by
When the number of spans is from 2 to 5, the frequency equations are listed in Table 5.

The empirical formulas for estimating natural frequencies of the bridge with nonuniform spans are listed aswhere and are modified factors, which can be seen in Figure 3.
(a) Modified factor
(b) Modified factor
It can be seen from Figure 3(a) that all of modified factors are not less than 1.0. In other words, the first natural frequency of the multispan continuous girder bridge is larger than that of the corresponding simply supported girder bridge. Because the side span length is smaller than the middle span length, the stiffness has been increased by the side span. In addition, the modified factor of the bridge with two spans is completely the same as that of the bridge with four spans. This is due to the free rotation of the middle support of the bridge with four spans, and the first vibration mode of fourspan bridge is equal of that of two twospan bridges.
It can be seen from the Figure 3(b) that the modified factor of 4span bridge is almost the same as that of 5span bridge. They are largely different from others and should be noted.
Finally, it has been emphasized that the empirical formulas and modified factors proposed in this paper can be used for estimating natural frequencies of bridges with uniform crosssection. Also, according to these formulas, the factors influencing natural frequencies are more clearly seen, which may give a better guideline for the design and evaluation of dynamic performance of bridges.
3. Dynamic Responses Analysis by Traditional Method
Research over the last 40 years has shown that the dynamic responses of bridges under moving vehicles are influenced by many parameters, such as the dynamic characteristics of the vehicle, the dynamic characteristics of the bridge, and variations in the surface conditions of the bridge and approach roadways [1]. However, most of the studies focused on the dynamic load allowance in midspan of the bridge. In this part, the dynamic load allowances (DLAs) in the side span and the middle span are discussed, respectively, to obtain the differences of DLAs in variance sections. Also, for comfort analysis, the accelerations of the bridge and the vehicle are investigated.
Based on the existing research results, thirteen parameters are considered here. According to the traditional method, every parameter is studied, respectively. In other words, if this parameter is changing, all of the other parameters are fixed. Therefore, the fundamental model of bridgevehicle coupled system is assumed at first, and then every parameter is varied from the basis of this model. They are listed in Tables 6, 7, and 8.
 
: length, : width, : moment of inertia, : modules of elasticity, : density, : damping ratio, and : maximum magnitude of the roughness. 

 
: the whole weight of the vehicle, : upper stiffness, : upper damping, : lower stiffness, and : lower stiffness. The values of all the parameters are obtained from the fundamental model (Table 7). 
According to the numerical simulation, the effect of different parameters on the dynamic response of the vehiclebridge system is investigated by our own program VBCVA (VehicleBridge Coupled Vibration Analysis). It has to be noted that the pavement roughness model proposed by Hwang and Nowak has been adopted in this program [39].
3.1. Effect of Roughness
For actual bridges, the pavement roughness is an inevitable factor. The effect of roughness on dynamic responses of the vehiclebridge system is shown in Figure 4.
(a) DLA of the bridge
(b) Acceleration of the bridge
(c) Acceleration of the vehicle
It can be seen from Figure 4 that, with the deterioration of bridge pavement conditions, the DLA of the bridge, the acceleration of the bridge, and the acceleration of the bridge grow in a linear fashion. Moreover, the dynamic response of the middle span is smaller than that of the side span. In addition, the growth rate in the middle span is a little smaller than that of the side span.
A great deal of research shows that the bridge pavement is one of the most vulnerable components due to its smaller stiffness [40, 41]. The results obtained from this study show that the roughness may significantly promote the vibration of the bridge traversed by moving vehicular loads. This will set up a vicious circle of increasing vibration, poorer roughness, and increased stress. Additionally, the larger vibration acceleration of the vehicle induced by the poorer roughness will make the passengers uncomfortable. Therefore, we should pay more attention to the bridge pavement condition in actual operation. Prompt repair of the damaged bridge pavement will reduce the unnecessary cost.
3.2. Effect of Span Length
As for the continuous girder bridges with uniform crosssection, the common span length ranges from 20 m to 40 m. It has to be emphasized that the crosssections of bridges with different span length are the same. The effect of span length on dynamic response of the vehiclebridge system is shown in Figure 5.
(a) DLA of the bridge
(b) Acceleration of the bridge
(c) Acceleration of the vehicle
It can be seen from Figure 5 that as the span length increases, the DLA of the middle span firstly decreases and then keeps steady, and the DLA of the side span is on the rise, which is largely different from the former one. But the change of the vibration acceleration in the middle span is the same as that of the vibration acceleration in the side span. Both of them slightly decrease and then increase. Additionally, the vibration acceleration of the vehicle rises at first and goes down later with the increasing of the span length.
For most dynamic responses in Figure 5, there is an extreme point with the increasing of the span length. It may be due to the resonance phenomenon between the bridge and the vehicle. The vibration natural frequency of the bridge with 25 m long span may be more close to that of the vehicle.
3.3. Effect of Span Ratio
From the aspect of mechanical characteristics due to the static load, the length of side span is not larger than that of middle span and not smaller than half of the length of the middle span. So the span ratio of the side span and the middle span ranges from 0.5 to 1.0 in this study. The effect of span ratio on dynamic responses of the vehiclebridge system is shown in Figure 6.
(a) DLA of the bridge
(b) Acceleration of the bridge
(c) Acceleration of the vehicle
It can be seen from Figure 6 that the change of the DLA is almost the same as that of the acceleration of the bridge. In the side span, they slightly decrease and then increase and then keep steady. In the middle span, they firstly go up and then go down. In addition, the vibration acceleration of the vehicle body keeps steady at first, and then it goes down. Moreover, it can be found that all of these dynamic responses reach their extreme value when the span ratio is 0.9. In China, as for the fabricated girder bridge, the span ratio is equal to 1.0 in common. And, as for the cast in place, the span ratio usually ranges from 0.6 to 0.8.
To determine the span ratio in design of the bridge, both the static and dynamic characteristics should be paid attention to.
3.4. Effect of Spans Number
For middle and smallspan continuous girder bridges, the number of spans is always ranging from 3 to 7, owning to the expansion or contraction according to the change of temperature. The effect of spans number on dynamic responses of the vehiclebridge system is shown in Figure 7.
(a) DLA of the bridge
(b) Acceleration of the bridge
(c) Acceleration of the vehicle
It can be seen from Figure 7 that the effect of spans number on the DLA of the bridge is not significant, especially for the DLA in the side span. However, the vibration acceleration of the bridge decreases with the increasing of the spans number. Moreover, the change of the vibration acceleration in the side span is much more similar to that in the middle span. Additionally, the vibration acceleration of the vehicle body is little influenced by the spans number of the bridge, as the fundamental natural frequency stays the same for bridges with different number of spans.
Based on the conclusion in this study, the dynamic performance is better for the bridge with more spans. However, the expansion or contraction induced by the change of the temperature may be larger when there are more spans. As a result, it has to be synthetically considered in the design.
3.5. Effect of Bridge Weight
The bridge weight is changed according to amending the density of the material. This method can make sure that the area and the stiffness of the crosssection of the bridge are invariable. The effect of bridge weight on dynamic responses of the vehiclebridge system is shown in Figure 8.
(a) DLA of the bridge
(b) Acceleration of the bridge
(c) Acceleration of the vehicle
It can be seen from Figure 8 that as the bridge weight increases, the DLA in the side span keeps steady, while the DLA in the middle span jumpily decreases. However, the change of the vibration acceleration in the side span and that in the middle span are almost the same. They slightly increase and then decrease. Additionally, the vibration acceleration of the vehicle body goes up and down, and the range is small.
In natural, it may be related to the dynamic characteristics of the bridge. The natural frequency is lower when the bridge weight is bigger. As for the vehiclebridge coupled vibration system, the resonance phenomenon will appear when its natural frequencies are so approximate.
3.6. Effect of Bridge Stiffness
The bridge stiffness is changed according to amending the elasticity modules of the material. This method can make sure that the area and the weight of the bridge are invariable. The effect of bridge stiffness on dynamic responses of the vehiclebridge system is shown in Figure 9.
(a) DLA of the bridge
(b) Acceleration of the bridge
(c) Acceleration of the vehicle
It can be seen from Figure 9 that, with the increasing of the bridge stiffness, the DLA in the side span continues to rise, while the DLA in the middle span slightly decreases and then increases. However, the change of the vibration acceleration in the side span and that in the middle span are almost the same. They fluctuate and the range is small. Additionally, the bridge stiffness has little influence on the vibration acceleration of the vehicle body.
Similarly, the natural frequency is higher when the bridge stiffness is bigger. Also, as for the vehiclebridge coupled vibration system, the resonance phenomenon will appear when its natural frequencies are so approximate.
3.7. Effect of Bridge Damping
As for the concrete structure, the damping ratio is assumed to be 0.05. In this study, it floats up and down at 40 percent. The effect of bridge damping on dynamic responses of the vehiclebridge system is shown in Figure 10.
(a) DLA of the bridge
(b) Acceleration of the bridge
(c) Acceleration of the vehicle
It can be seen from Figure 10 that the damping ratio has little influence on the DLA of the bridge. But the vibration acceleration of the bridge significantly decreases with the increasing of the damping ratio. Similarly, the vibration acceleration of the vehicle body continues to go down as the damping ratio increases.
Damping is one of the most important dynamic characteristics of the bridge. Up till now, it cannot be obtained by calculation. The only way to obtain the damping of the structure is by measuring in the fielding test. However, according to the results in this paper, the damping is not so important for calculating the DLA. But it is closely related to the vibration acceleration of the bridge and the vehicle, which may influence the inertia force of the bridge and the riding comfort. Therefore, the damping should still be noticed.
3.8. Effect of Vehicle Weight
To investigate the effect of the vehicle weight on dynamic responses of the bridge, the factor of the vehicle weight is introduced, and the fundamental weight is assumed to be 35 tons. In this study, the factor ranges from 0.6 to 1.4 at the step of 0.2. The effect of vehicle weight on dynamic responses of the vehiclebridge system is shown in Figure 11.
(a) DLA of the bridge
(b) Acceleration of the bridge
(c) Acceleration of the vehicle
It can be seen from Figure 11 that when the factor is smaller than 1.0, the DLA in the side span slightly decreases and then increases, while the DLA in the middle span firstly decreases and then keeps steady. When the factor is bigger than 1.0, the DLA in the side span decreases with the increasing of the vehicle weight, but the DLA in the middle span increases with the increasing of the vehicle weight. However, as the vehicle weight increases, the vibration acceleration in the side span goes down, while that in the middle span goes up. In addition, the vibration acceleration of the vehicle body decreases with the increasing of the vehicle weight.
Obviously, the effect of vehicle weight on dynamic responses in different positions is not the same. Limiting the vehicle weight can effectively reduce the static stress of the bridge, but its influence on the dynamic responses is not clearly enough. In particular, the change of the vibration acceleration in different positions of the bridge is completely the opposite.
3.9. Effect of Upper Stiffness
To investigate the effect of the upper stiffness on dynamic responses of the vehiclebridge system, the factor of the upper stiffness is introduced. In this study, the factor ranges from 0.6 to 1.4 at the step of 0.2. The effect of upper stiffness on dynamic responses of the vehiclebridge system is shown in Figure 12.
(a) DLA of the bridge
(b) Acceleration of the bridge
(c) Acceleration of the vehicle
It can be seen from Figure 12 that the DLA of the bridge almost goes down with the increasing of the upper stiffness of the vehicle. However, as the upper stiffness of the vehicle increases, the vibration acceleration of the bridge goes up. Moreover, the change of the acceleration in the side span is faster than that in the middle span. In other words, the effect of the upper stiffness on the dynamic response of the side span is much more sensitive. In addition, the vibration acceleration of the vehicle body continues to go up with the increasing of the upper stiffness of the vehicle.
As a result, the upper stiffness of the vehicle not only influences the dynamic responses of the bridge but also significantly affects the vibration acceleration of the vehicle body, which is closely related to the comfort of passengers and the safety of the goods. Therefore, to find the friendly stiffness of the vehicle is urgent and efficient in the area of riding comfort analysis.
3.10. Effect of Upper Damping
To investigate the effect of the upper damping on dynamic responses of the vehiclebridge system, the factor of the upper damping is introduced. In this study, the factor ranges from 0.6 to 1.4 at the step of 0.2. The effect of upper damping on dynamic responses of the vehiclebridge system is shown in Figure 13.
(a) DLA of the bridge
(b) Acceleration of the bridge
(c) Acceleration of the vehicle
It can be seen from Figure 13 that all of the dynamic responses of the vehiclebridge system decease with the increasing of the upper damping of the vehicle. The change of the DLA in the side span is faster than that in the middle span. However, the change of the vibration acceleration in the side span is almost the same as that in the middle span. It has to be emphasized that the influence of the upper damping on the vibration acceleration of the vehicle is the most significant. In other words, the shock absorber system of the vehicle is so important in its design.
3.11. Effect of Lower Stiffness
To investigate the effect of the lower stiffness on dynamic responses of the vehiclebridge system, the factor of the lower stiffness is introduced. In this study, the factor ranges from 0.6 to 1.4 at the step of 0.2. The effect of lower stiffness on dynamic responses of the vehiclebridge system is shown in Figure 14.
(a) DLA of the bridge
(b) Acceleration of the bridge
(c) Acceleration of the vehicle
It can be seen from Figure 14 that as the lower stiffness of the vehicle increases, the change of the DLA in the side span and that in the middle span are entirely the opposite. The former one slightly increases and then decreases, while the latter one slightly decreases and then increases. However, the vibration acceleration of the bridge goes up with the increasing of the lower stiffness. Moreover, the change of the vibration acceleration in the side span is faster than that in the middle span. Additionally, the vibration acceleration of the vehicle body goes up firstly and then keeps steady.
Comparing Figure 12 with Figure 14, it can be found that the effect of the upper stiffness is almost the same as the effect of the lower stiffness. Moreover, the former one is more significant than the latter one.
3.12. Effect of Lower Damping
To investigate the effect of the lower damping on dynamic responses of the vehiclebridge system, the factor of the lower damping is introduced. In this study, the factor ranges from 0.6 to 1.4 at the step of 0.2. The effect of lower damping on dynamic responses of the vehiclebridge system is shown in Figure 15.
(a) DLA of the bridge
(b) Acceleration of the bridge
(c) Acceleration of the vehicle
It can be seen from Figure 15 that as the lower damping increases, all of dynamic responses of the vehiclebridge system are falling. But the changing range of them is small. Also, the dynamic responses in the side span are significantly larger than those in the middle span, especially for the DLA.
Comparing Figure 13 with Figure 15, it can be found that the effect of the upper damping is almost the same as the effect of the lower damping. Moreover, the former one is more significant than the latter one.
3.13. Effect of Vehicle Speed
As for the loading truck, the speed ranges from 20 km/h to 120 km/h at the step of 20 km/h. The effect of vehicle speed on dynamic responses of the vehiclebridge system is shown in Figure 16.
(a) DLA of the bridge
(b) Acceleration of the bridge
(c) Acceleration of the vehicle
It can be seen from Figure 16 that as the speed increases, both the DLA in the side span and the DLA in the middle span increase firstly and then decrease. But the critical speeds of them are not the same, and the former one is bigger than the latter one. However, the vibration acceleration of the bridge goes up with the increasing of the vehicle speed. Similarly, the vibration acceleration of the vehicle body also increases.
Therefore, it has to be noted that limiting the speed directly is not a rational way to ensure the safety of the bridge, which is the common way during the maintenance of the old or damaged bridge.
Obviously, the effect of every parameter on the dynamic responses of the vehiclebridge system can be studied by the program VBCVA. When one parameter is studied, the other parameters should be fixed. But, according to this method, the results may not occasionally be the same as the actual condition, because the interaction effects among different parameters are not considered.
4. Numerical Simulations Based on the Orthogonal Experimental Design
There are two ways for the analysis of factors influencing the dynamic responses of the vehiclebridge coupled vibration system, full factorial designs and fractional factorial designs. Full factorial designs reveal whether the effect of each factor depends on the levels of other factors in the experiment. This is the primary reason for multifactor experiments. One factorial experiment can show “interaction effects” that a series of experiments each involving a single factor cannot [42]. However, fractional factorial designs permit investigation of the effects of many factors in fewer runs than a full factorial design.
Data from fractional factorial designs are interpreted based on the following two assumptions. The first one is the sparsity of important effects. Only a few of the many possible effects are prominent. Even if many effects are nonzero, we expect a few to stand out as much larger than the rest. The other one is the simplicity of important effects. Main effects and/or twofactor interactions are more likely to be important than highorder interactions. This is also known as the hierarchical ordering principle.
Therefore, considering the high efficiency, the orthogonal experimental design, one common type of the fractional factorial design, is adopted in this study, including the numerical simulation without interaction and the numerical simulation with interaction.
4.1. Orthogonal Experimental Design
According to the orthogonal experimental design, two objectives can be attained. The first objective is to select fewer typical combinations for the test. And the other objective is to obtain the correct conclusions by the scientific method based on the limited combinations. The orthogonal experimental design consists of the design of test plan and the process of test data.
As for the design of test plan, the main procedures are listed as follows [43]. Firstly, the objective and the test index are clearly proposed. Secondly, the factors and their levels are determined based on the comprehension in this research field other than the knowledge of statistics. Thirdly, the rational orthogonal table is selected, including the table considering interaction and the table without interaction. Fourthly, the table header is designed. Finally, the test plan is formed.
As for the process of the test data, two methods are adopted, including the range analysis and the variance analysis [44]. The range analysis is much simpler and more convenient. However, the test error cannot be estimated, and the reliability of the results cannot be determined. Also, it cannot be applied in the fields of regression analysis and regression design. But the variance analysis can remedy the defects above.
The method of range analysis is also called the visual analysis method and the (the abbreviation of the range) method. It consists of two procedures, calculation and judgment. The diagram can be seen in Figure 17.
In Figure 17 the is the sum of the test index with the th level in the th column and is the average of the . In addition, the is the difference between the maximum and the minimum value of the : The is represented for the changing amplitude of the test index with the changing of the levels of the th factor. Therefore, it can be concluded that the influence of the th factor on the test index is much more significant if the value of is bigger. Sometimes, for more clarity, the trend plot is given.
Variance analysis is more rigorous. There are four steps to realize the variance analysis [45, 46].
Step 1. As for every factor, calculate the sum of square of deviations (), the degree of freedom (dof, ), and the variance estimation ().
Step 2. Estimate the variance of the error ().
Step 3. Obtain the test static and compare with its critical value for given significance level .
Step 4. For simplicity, the variance analysis table is listed, including the process and the results. Considerin which is the index result of the th run and is the index result of the th factor with the th level. Also, and denote, respectively, the sum of square of deviations and the degree of freedom of the error. It has to be noted that the error has resulted from all of the vacant columns in the orthogonal array. Also, the accuracy increases with the increasing dof of the error [47]. Therefore, if the significance level of one factor is larger than 0.25, it can be included as the error.
In this study, the orthogonal table (3^{13}) is used for the numerical simulation without interaction, while the orthogonal table (2^{15}) is used for the numerical simulation with interaction.
4.2. Numerical Simulation without Interaction
There are thirteen factors possibly affecting dynamic responses of the vehiclebridge coupled vibration system. A large number of studies [8] have proved that the roughness is the most significant factor. Therefore, the roughness is not arranged in the orthogonal table. The other twelve factors are arranged in the orthogonal table, once in a column. And the residual column is thought as the test error.
The results obtained from the VBCVA are listed in Table 9. The range analysis of the data can be seen in Table 10 and Figure 18. The variance analysis of the data is shown in Table 11.


(a) Dynamic load allowance in the side span ()  
 
(b) Dynamic load allowance in the middle span ()  
 
(c) Vibration acceleration in the side span ()  
 
(d) Vibration acceleration in the middle span ()  
 
(e) Vibration acceleration of the vehicle body ()  

(a) Trend plot between factors and the DLA in side span
(b) Trend plot between factors and the DLA in middle span
(c) Trend plot between factors and the acceleration in side span
(d) Trend plot between factors and the acceleration in middle span
(e) Trend plot between factors and the acceleration of the vehicle body
It can be concluded that the influence of factors on different indices is listed as follows.(i)For DLA in the side span of the bridge, consider(ii)For DLA in the middle span of the bridge, consider(iii)For vibration acceleration in the side span of the bridge, consider(iv)For vibration acceleration in the middle span of the bridge, consider(v)For vibration acceleration of the vehicle body, consider
It has to be noted that the roughness is assumed as the significant factor. Obviously, among other factors, the most important factors affecting the different dynamic responses are not the same.
4.3. Numerical Simulation with Interaction
As mentioned earlier, the interaction almost exists in all of physical phenomena. When the interaction is so small, it can be ignored in application. However, as for research, we do not know whether the interaction can be ignored or not at first.
Based on the results from the above section, eight most important factors influencing the DLA are selected, and some interactions between them are investigated. They are the pavement roughness (A), the length of the main span (B), the ratio between the side span and the middle span (C), the number of spans (D), the weight of the bridge (E), the weight of the vehicle (H), the upper stiffness of the vehicle (I), and the upper damping of the vehicle (J). Meanwhile, due to the calculation cost and the existing orthogonal array, the number of levels is determined as two for each factor. To avoid the mixture, the most important factor is arranged at first. It can be seen in Table 12.

The results obtained from the VBCVA are listed in Table 13. The range analysis of the data can be seen in Table 14. The variance analysis of the data is shown in Table 15.


(a) Dynamic load allowance in the side span ()  
 
(b) Dynamic load allowance in the middle span ()  
 
(c) Vibration acceleration in the side span ()  
 
(d) Vibration acceleration in the middle span ()  
 
(e) Vibration acceleration of the vehicle body ()  

It can be concluded that the influence of factors on different indices is listed as follows.(i)For DLA in the side span of the bridge, consider