Research Article  Open Access
Dynamic Responses of Simply Supported Girder Bridges to Moving Vehicular Loads Based on Mathematical Methods
Abstract
For dynamic responses of highway bridges to moving vehicles, most of studies focused on singlefactor analysis or multifactor analysis based on full factorial design. The defect of the former one is that it has no consideration of interaction effects, while that of the latter one is that it has large calculation. To avoid these defects, simplified theoretical derivations are presented at first; then some numerical simulations based on the proposed method of the orthogonal experimental design in batches have been done by our own program VBCVA. According to simplified theoretical derivations, three factors (κ, γ, and α) are proved as the most important factors to determine dynamic responses. Based on the modal synthesis method, the program VBCVA has been introduced in detail. Then on the basis of the orthogonal experimental design, both main effects and interaction effects are studied. The results show that, for different indices of dynamic responses, the influences of each factor are not the same. Additionally, the interaction effects have proved to be so small that they can be neglected. In the end, this method provides a good way to obtain more rational empirical formulas of the DLA and other dynamic responses, which may be adopted in the revision of codes for design and evaluation.
1. Introduction
Research on the dynamic analysis of the vehiclebridge coupled vibration system is an important issue in civil engineering [1]. As for this problem, three methods are used in general, including the theoretical derivation, dynamic loading test, and numerical simulation.
As moving vehicles on the bridge vary in both time and space, the problem becomes more complex. And it has been noted more than 100 years before. A lot of researchers tried to obtain the effects of moving load on various elements, components, and structures based on theoretical derivations. It had been well reviewed by Frýba [2]. In 1989, an analyticalnumerical method was presented that could be used to determine the dynamic behavior of beams, with different boundary conditions, carrying a moving mass [3].
The specification about impact factor or dynamic load allowance in most of design codes was obtained based on the dynamic loading test. The first thorough investigation of highway bridge dynamic loading was conducted from 1922 to 1928 by an ASCE committee (10 bridges) [4]. In the years 1958 to 1981, the Section Concrete Structure of the EMPA performed load tests on 356 bridges, and 226 static and dynamic tests on slab and beamtype highway bridges were of interest at last [5]. In Canada, three large scale series tests were completed. Firstly, a group of 52 bridges known to vibrate was selected for test in the years 1956 to 1957. Secondly, many tests were completed on continuous concrete bridges in the years 1969 to 1971. Thirdly, a total of 27 bridges were selected in 1980, which was the main basis of the design code OHBDC [6]. In China, the specification in current code was obtained from the test data of 7 simply supported girder bridges [7]. In Korea, from 1995 to 2007, a total of 256 bridges, of which span varied from 10 m to 160 m, were tested and analyzed to suggest new design criteria for an impact factor which was based on a natural frequency rather than the span length [8].
However, not all researchers can study the problem by dynamic fielding tests due to the budgets limitations. In recent years, with the development of computers, numerical simulation studies are widely used. The dynamic behavior of simplespan [9] and continuous [10] highway bridges under moving vehicles were studied based on a method of analysis which idealized the bridge as a single beam and represents the vehicle as a multiaxle sprung load. Wang et al. [11] studied the dynamic behavior of multigirder bridges, which was modeled as a grillage beam system, due to vehicles with 7 and 12 degrees of freedom (H2044 truck and HS2044 truck) moving across rough bridge decks. A general and efficient method was proposed for the resolution of the dynamic interaction problem between a bridge, discretized by a threedimensional finite element model, and a dynamic system of vehicles running at a prescribed speed [12]. To study the interaction problem of large complicated bridges with various types of running vehicles, a fully computerized approach for assembling equations of motion of any types of coupled vehiclebridge system [13]. Apart of this, many researchers focused on more complicated modeling of the vehiclebridge system [14–16], which are much more similar with the actual conditions.
According to the review of existing literature, the defects and shortages are listed as follows.(i)The vehiclebridge coupled model is oversimplified, especially in theoretical derivation. And it is largely different from the actual conditions.(ii)The selection of influence factors is subjective, and it lacks sufficient theoretical basis. As for one problem, the understanding from different engineers may be largely different, even if opposite.(iii)The interaction effects are seldom studied in the dynamic analysis of the bridge to moving vehicles. Most of studies focused on main effects.(iv)Abundant calculations are needed in the method of full factorial design, which has been used for multifactors problem before.(v)The impact factor in current code is the function of single parameter. It has been proved not rational and should be revised.
Therefore, the method of orthogonal experimental design in batches 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. To make a good understanding of the vehiclebridge system and to obtain the basis of the selection of some important factors, two simplified models are discussed at first. As a basic calculation tool, which will be more consistent with the actual condition, our own program VBCVA is described in detail. At last, using the proposed method, the influences of twelve common factors and some of their interaction effects on dynamic responses of the vehiclebridge coupled system are discussed.
For safety of bridges, the impact factor is usually adopted to account for the dynamic responses induced by moving vehicles. However, there are many different names [17], which are always confused in the research and engineering application. So we should restate it here. In this study, two terms are used, the dynamic load allowance (DLA, ) and the impact factor (IF, ). And the definitions of them can be given by in which the and denote, respectively, the maximum value of the dynamic response and the static response. It has to be noted that the position of the vehicle where the maximum dynamic response occurs is different from that related with the maximum static response in general.
In addition, for comfort analysis of pedestrians and passengers or drivers, the vibration accelerations of the bridge and the vehicle are adopted, because a large number of researchers have obtained the conclusion that the comfort of people is mainly determined by the acceleration [18–22].
2. Theoretical Derivation on Simplified Models
Theoretical derivation on simplified models does not quite agree with the actual situation, but it is one of the best ways to determine the key parameters and their influences on the dynamic responses. And, according to the derivation, the physical meanings may be more obvious, which will be good guidance for design and evaluation on dynamic performance of the highway bridge to moving vehicular loads. Consequently, in the beginning of this study, the dynamic responses of a simply supported girder bridge traversed by a single moving constant load and a moving sprung mass which are deduced and discussed, respectively.
2.1. Simply Supported Girder Bridges Traversed by a Single Moving Constant Load
Of the wide range of problems involving vibration of structures subjected to a moving load, the easiest one to tackle is that of dynamic responses in a simply supported girder bridge (or a simply supported beam), traversed by a constant force moving at uniform speed [23]. This classical case was first solved by Krylov [24], then by Timoshenko [25]. Other solutions worthy of mention are those by Inglis [26] and Koloušek and McLean [27].
If the weight of the vehicle is far smaller than that of the bridge, the inertia force of the vehicle can be ignored. Then the vehiclebridge coupled system can be simplified as a simply supported girder bridge traversed by a single moving constant load (Figure 1).
Based on the theory of the structural dynamics, the governing equation of the bridge can be given by in which the denotes the constant force, , , , , , , , , and v denote, respectively, the span length, mass per unit length, damping, moment of inertia, elasticity modulus, displacement, location, time, and the moving speed. And is the Diracdelta function.
The initial displacement and velocity are assumed as zero. As for simply supported girder bridge, the sine functions can be assumed as the mode shapes . Also, the dynamic displacement and the bending moment in the support section are zero. Based on the mode superposition method, the dynamic displacement of the bridge can be given by where is the damping ratio corresponding to the th mode shape and , , and are, respectively, the th natural frequency of the bridge without damping and with damping and the disturbance frequency which is dependent on the speed of the moving vehicle and the th mode shape of the bridge. And is the number of all considered mode shapes. Furthermore, an important nondimensional parameter is introduced. The relations of these parameters are listed as follows:
According to (3), if the damping of the bridge is ignored, the dynamic displacement can be simplified as It can be seen from (5) that there are two parts in the dynamic displacement of the simply supported girder bridge traversed by a single moving constant load. One part is the forced vibration related to the speed of the moving load, and the other part is the free vibration of the bridge itself.
If the speed of moving load is zero, but the force vibration in (5) is still retained, the static displacement of the bridge () is obtained as where and are, respectively, the locations of the bridge section and the load.
In accordance with theory of structural statics [28], the accurate static displacement can be given by Considering the symmetry, four sections of the bridge are selected, and the comparison of (6) and (7) is plotted in Figure 2.
Figure 2 shows that they are almost the same. And the assumption of sine functions being thought as mode shapes is verified rational enough in another aspect:
2.2. A Simply Supported Girder Bridge Traversed by a Moving Sprung Mass
Due to neglect of the mass of the vehicle in the section above, coupled vibration in the vehiclebridge system is not considered. To make a good comprehension of this effect, a simply supported girder bridge traversed by a moving sprung mass has been adopted in this section [29] (Figure 3).
In Figure 3, the moving vehicle has been simplified as a sprung mass, including the mass of vehicle body and the mass of the wheel . The vehicle body and the wheel are connected with a spring and a damper. Also, it is assumed that the wheel is contacted with the bridge all the time. In addition, the bridge is thought to be vibrated based on only the first mode shape. Likewise, there is no initial vibration of the bridge, and the damping of the bridge is ignored.
Respectively, based on the theory of structural dynamics, the governing equations of the vehicle body and the bridge can be given bywhere is the displacement of the contact point and , , and are, respectively, generalized mass, stiffness, and load of the bridge. And all of these parameters can be calculated by the following equations:
Substituting ((10c)~(10h)) into (9b) and combining (9a) and (9b), then
In general, the mass of the wheel is much less than that of the vehicle body, so it can be ignored. And when the damping of the vehicle is neglected, (11) can be simplified as
Obviously, the fundamental frequency and the maximum static displacement of the vehicle and the bridge can be given by
In accordance with the current codes, the dynamic load allowance (DLA, ) is defined as the ratio of the maximum dynamic response and the maximum static response. And the impact factor (IF, ) in the code of China (JTG D602004) is equal to DLA subtracted by one as follows: where is the ratio between the mass of the vehicle and that of the bridge. Besides, another two important nondimensional parameters are introduced here:
Then (12) can be written as
It can be seen from (17) that the most significant parameters that influenced the IF or the DLA of the bridge are , , , and .
Now, the influence of the natural frequency of the vehicle is discussed in detail. Equation (17) shows that the influence effect has two parts, including the term in the mass matrix and the term in the stiffness matrix. Obviously, when the natural frequency changes, due to the characteristics of the sine function, only the time corresponding to the maximum response has changed other than the maximum value itself. In addition, for fixed value of , , and , the maximum dynamic load allowance can be obtained. And then the corresponding accelerations are written aswhere and are only dependent on the values of , , and , other than the frequency . Then the first part of (17) is not related to the frequency, as the frequency has disappeared when the mass matrix is multiplied by the acceleration vector.
Totally, based on simplified model and derivation above, the most significant influence parameters may be , , and , defined in (16a) to (16c).
3. Program for VehicleBridge Coupled Vibration Analysis
Based on their own advantages of existing generalized commercial software ANSYS and MATLAB, the program VBCVA (VehicleBridge Coupled Vibration Analysis) has been developed by our own research group. It is used for vibration and dynamic analysis of the vehiclebridge coupled system. Upon the modal synthesis method, the bridge model is found by ANSYS, while the vehicle model and the roughness model are established by MATLAB.
3.1. Bridge Equations
According to the theory of structural dynamics [30], the governing equation of the bridge can be given by in which the denotes the load vector induced by the moving vehicles and [], [], and [] denote, respectively, the mass matrix, damping matrix, and stiffness matrix. Furthermore, is the displacement of the bridge, the first derivative of the displacement is the vibration velocity, and the second derivative of the displacement is the vibration acceleration. It is worth noting that all these symbols and this equation are described in Cartesian coordinate system.
As there are various types of highway bridges and they are much more complicated with the increasing technology, the accurate modeling of the bridge may be difficult to realize by our own written program. In addition, this will hinder the generalization and the development of that program. And the limitation can be obviously seen when the equations of various types of bridges are different. Therefore, the modal synthesis method has been adopted. Another advantage of this translation is the reduction of the degree number: where [Φ] is the mode shape matrix of the bridge and [] is the coordinate in the modal coordinate system, which denotes the contribution of every mode shapes.
Substituting (20) into (19), premultiplied by the transposed matrix of the mode shapes, the following equation is obtained: where [], [], [], and are mass matrix, damping matrix, stiffness matrix, and load vector in the modal coordinate system. They are given by
For convenience in application, the matrix of mode shape obtained from ANSYS is normalized. Then the following matrixes are used in the program VBCVA: where and are the damping ratio and the natural frequency of the th mode shape.
3.2. Vehicle Equations
Actually, the moving vehicle can be looked at as a vibration system with multidegrees of freedom. In this study, D’Alembert’s principle is used for deducing dynamic equations of spatial vehicle models. And the position induced by the static weight of the vehicle is selected as the reference position [31].
The schematic plot of the vehicle model can be seen in Figure 4. If the number of axles is , the degree of freedom (dof) is , including the vertical dof , pitching dof , swung dof of the vehicle body, and the vertical dof of every wheels. It can be seen from the Figure 3 that in which and are the coordinates of the centreofgravity of the vehicle and and denote, respectively, the distance from the wheel to the centre in longitudinal and transverse direction. Furthermore, is the distance from the wheel to the front axle in longitudinal direction, and is the wheel space.
(a) Elevation view
(b) Front view
There are some assumptions on the vehicle model. The wheel and the bridge will contact with each other all the time. Only vertical effects between the vehicle and the bridge are considered, while longitudinal and transverse effects are ignored. The vehicle body and all wheels are assumed as rigid bodies with corresponding mass, while the spring and the damper are linear [32].
Then the vehicle equation can be given by in which the denotes the load vector induced by the bridge and [], [], and [] denote, respectively, the mass matrix, damping matrix, and stiffness matrix of the vehicle. Furthermore, is the displacement of the vehicle, the first derivative of the displacement is the vibration velocity, and the second derivative of the displacement is the vibration acceleration. They are listed as follows. It has to be noted that the stiffness matrix [] is just similar with the damping matrix [], replacing the symbol with the symbol :
3.3. Coupled Equations
The interaction force between the bridge and the vehicle can be considered from two aspects. One is the force of the vehicle induced by the vibration of the bridge, and the other one is the force of the bridge caused by the vehicle, including the static weight and the dynamic force:where denotes the sum of two parts, the first part is the gravity of the th wheel itself, and the other part is the distribution of the weight of the vehicle body on the th wheel. And is the relative displacement of the contact point between the th wheel and the bridge. Also, the first derivative of the displacement is the vibration velocity. They can be obtained by in which the denotes the vertical displacement of the th wheel and and denote, respectively, the vertical displacement and the rotation of the bridge at the location of the contact point. And, is the vertical roughness, and is the speed of the moving vehicle.
Coupled equations of the vehiclebridge vibration system can be obtained by the combination of the bridge equation and the vehicle equation above. They are listed as follows:
where the denotes the number of selected mode shapes and , , and denote, respectively, the total number of dof of all vehicles, the number of dof of the vehicle body of all vehicles, and the number of dof of all wheels. Other symbols are the same as before.
3.4. Solution of Coupled Equations
When a vehicle goes across the bridge, the position of the contact point changes with time. Therefore, the coupled equations are timevarying system of differential equations. And it is difficult to obtain the closedform solutions. But they can be solved by some numerical methods, such as the central difference method, Newmark method, Wilson method, and so on [33–35]. Wilson method has been adopted in this study. It is a selfstabilization method, which is independent on the integration step and the shortest period [36, 37]. And it also can be looked at as the modified linear acceleration method.
It is assumed that the acceleration during the period [, ] is linearly varied. Firstly, the vibration of the system at the time of is calculated based on the linear acceleration method, and then the vibration of the system at the time of is deduced using the interpolation formula. And it is worth noting that this method has been verified unconditionally stable if the parameter is larger than 1.37.
Based on the assumption of linear acceleration, the acceleration during this period [, ] is given by According to the integration, the vibration velocity and the dynamic displacement can be shown as follows: When the time is equal to , the corresponding vibration velocity and the dynamic displacement areTo solve (32a) and (32b), the acceleration and the dynamic displacement at the time of can be listed as follows:
Of course, at the time of , the following vibration equation of the system should be satisfied: in which the load vector can be obtained by the method of linear extrapolation:
Substituting (33a), (33b), and (35) into (34) and rearranging, the following equation is obtained: where
According to solving (36), the dynamic displacement at the time of () can be known. Substituting into (33a), the acceleration at the time of can be obtained. Then substituting it into (31a) and using the designation , the acceleration at the time of is gained as
Similarly, using the designations and , the vibration velocity and the dynamic displacement at the time of are listed as follows:
3.5. Roughness Model
In both design and evaluation, pavement roughness is the primary factor affecting the dynamic performances of the bridge traversed by the designated vehicles [38]. And, in the dynamic responses analysis of highway bridges under moving vehicles, it is considered to be a cause of vehicle vibration [39]. Deng and Cai have proposed that, due to the road surface deterioration of existing bridges, the calculated impact factors form field measurements could be higher than the values specified in design codes that mainly target new bridge design [15]. So it is necessary to include the pavement roughness in our own program.
There are two methods to take the roughness into account, field measurement and numerical simulation. As for the former way, it is measured generally by one of the following two methods, that is, (1) by using a profilometer or (2) by calculating pavement roughness backwards from vibration data of the wellresearched dynamic properties of the vehicle [36]. Nowadays, based on much data from the field measurement, more researchers begin to admit the fact that the roughness is a realization of a random process that can be described by a power spectral density (PSD) function. And the pavement roughness model proposed by Hwang and Nowak [40] has been adopted in this study.
Typical PSD function can be approximated by an exponential function: in which the denotes roughness coefficient and , , , and denote, respectively, spatial frequency (m^{−1}), lower limit, upper limit, and spectral shape index. It is assumed that the pavement roughness can be modeled as a stationary Gaussian random process. Therefore, it can be generated by an inverse Fourier transform: where is PSD function and is random number uniformly distributed from 0 to 2π.
The process of generating pavement roughness is listed as follows:in which the denotes the position in the longitudinal direction of the bridge, is the roughness, and is the sampling number.
Considering common seen speeds of vehicles in highway bridges, the lower and upper limits of the spatial frequency are specified as 0.05 m^{−1} and 3.00 m^{−1} [7]. Furthermore, is usually equal to 512, and the spectral shape index is 2 in general. Also, the selected value of the roughness coefficient can be seen in Table 1.

As for different conditions, five typical pavement roughness examples are obtained using the program VBCVA. And the results are plotted in Figure 5.
(a) Very good ()
(b) Good ()
(c) Average ()
(d) Poor ()
(e) Very poor ()
In addition, the relation between the maximum amplitude of the roughness sample and the roughness coefficient is studied. It can be seen from Figure 6 that they are square relation, which is consistent with ((41a), (41b), (41c), (41d), and (41e)).
3.6. Flowchart of the Program
Due to various types of bridges, the modal synthesis method is adopted in the program VBCVA for general use. At first, the dynamic characteristics of the highway bridge, including natural frequencies and mode shapes, are obtained based on the finite element model (FEM) built by the commercial software ANSYS. Meanwhile, the data files for vehicles are prepared by another commercial software MATLAB. Then the coupled equations are calculated using the Wilson method [41]. It is developed following the flowchart in Figure 7.
To supplement the application of this program, it is necessary to emphasize that the data files preparation includes bridge data (frequencies, mode shapes, damping ratios, and coordinates of nodes), pavement roughness model, and vehicle data (vehicle type, speed, number, and initial location). Also, controlling parameters mainly means in the method of Wilson and the integration step. After all calculations, the postprocessing is done by MATLAB, including plots of dynamic responses of the designated degree of the bridge or vehicles.
The program can be used for calculating the cases of multilanes and multivehicles (both in longitudinal direction and in transverse direction). Also, the number of axles or vehicles is not limited. Furthermore, the validity and the rationality have been verified by some numerical and experimental results from the existing papers [42]. Therefore, this program has to be thought of as convenient and powerful enough for the analysis of the vehiclebridge coupled vibration problem.
4. Numerical Simulations Based on the Orthogonal Experimental Design
Concerning dynamic responses of highway bridges to moving vehicular loads, three methods are commonly adopted, including theoretical derivation, dynamic loading test, and numerical simulations. As we know, the first method can only be appropriate for some simplified models, while the last two methods can be applied to the actual situation in general. However, due to the economy and many other causes, most of studies focused on dynamic responses induced by much fewer factors. In addition, influence of interactions between these factors on dynamic responses is barely studied. Therefore, based on orthogonal experimental design, influences of many factors and their interactions on the dynamic responses of simply supported girder bridges to moving vehicular loads are discussed at length in this section.
4.1. Orthogonal Experimental Design
As for the multifactor experimental problem, the full factorial design was widely accepted in the early stage. There are two primary benefits of this method. Firstly, it reveals whether the effect of each factor depends on the levels of other factors in the experiment. And one factorial experiment can show “interaction effects” that a series of experiments each involving a single factor cannot. Secondly, it provides excellent precision for the regression model parameter estimates that summarize the combined effects of the factors [43].
However, when the variables are more, the size of the full factorial design is much extremely larger. For example, with three levels of every factor, increasing the number of factors from 2 to 4 increases the size of the full factorial design from 3^{2} to 3^{4}, 9 times larger. So fractional factorial designs have come out to largely reduce the work in the analysis of multifactor problem. In statistics, fractional factorial designs consist of a carefully chosen subset of the experimental runs of a full factorial design. That is chosen so as to exploit the sparsityofeffects principle to expose information about the most important features of the problem studied, while using a fraction of the effort of a full factorial design in terms of experimental runs.
Obviously, there are two principal contradictions in the fractional factorial design. One is the contradiction between the larger runs in the full factorial design and the expected smaller runs in actual operation. The other one is that between the smaller runs and the expected whole information, which is just the same as that obtained from the full factorial design [44]. To solve these two contradictions, the orthogonal experimental design is one of the best methods. According to a more rational arrangement, such as the orthogonal array, a minority of runs that are most typical can be selected, which can solve the first contradiction. Using some scientific methods of data processing, including the range analysis and the variance analysis, many reasonable conclusions can be obtained based on the minority of experiments, which can solve the second contradiction. In a word, the orthogonal experimental design has been verified as a good method.
The orthogonal array is the basis of the orthogonal experimental design. And it has been constructed based on these two mathematical courses, combinatorics and probability. Of course, it is not necessary for engineers to know how to construct it. We only need to transplant existing orthogonal arrays into engineering areas. The orthogonal array is usually denoted by , in which is the code of orthogonal array, is the number of experimental runs, is the number of levels, and is the number of factors. In the recent period, this type of construction thought has become widespread.
It has been proved that the “interaction effects” are always seen in the multifactor experiment. They are defined as the influence of combined factors on the test index. In fact, the “interaction effects” are reflections of mutual promotion or inhibition, and these effects are more or less existed in all physical phenomenon. In the statistical analysis of the results obtained from factorial experiments, the sparsityofeffects principle states that a system is usually dominated by main effects (single factor) and loworder interactions. Therefore, it is most likely that main effects and twofactor interactions are the most significant responses in an experiment [45].
A majority of studies have shown that the dynamic responses of vehiclebridge system are influenced by so many factors, including vehicle characteristics, pavement roughness, and bridge characteristics [46]. In this study, to efficiently reduce the experimental runs, the conventional orthogonal design is divided into two phases. In the first phase, main effects (single factor) are analyzed without any interaction effects. Based on the results from the first phase, the interaction effects of some of the most important factors are discussed in the second phase.
4.2. Samples Selection
As we know, there are many types of highway bridges, including the girder bridge, the arch bridge, the cablestayed bridge, and the suspension bridge. Of course, the simply supported girder bridge is the most common type. Based on the standard drawings issued by the Ministry of Transport of the People’s Republic of China in 2008, a 30 mspan bridge with smallbox section has been selected as the fundamental sample of bridges. The damping ratio is assumed as 0.05. The crosssection is in Figure 8 and its other parameters are listed in Table 2.
 
: length, : width, : moment of inertia, : modules of elasticity, and : density. 
In highway bridges, the types of moving vehicles are various, which is significantly different from the type of vehicles in railway bridges. And it makes the problem of the vehiclebridge coupled system become much more complex. Usually, there are three types of vehicles, including the small car, the bus or coach, and the loading truck. Apparently, due to their light weight, the first two types of vehicles may not be the main contributions on the dynamic responses. So the 3axle loading truck has to be selected in this study. Its weight is 30 tons. It is seen in Figure 9 and Table 3.

Pavement roughness is the excitation source of the vehiclebridge coupled vibration. For new bridges, the condition of the pavement is very good or good. So the maximum amplitude of the pavement roughness is assumed as 1cm () in this study.
Based on so many existing research results [2, 4, 5, 7–9], twelve factors are considered and discussed, including the pavement roughness, the span length, width, mass, stiffness, and damping of the bridge, the mass, upper stiffness (suspension system), upper damping, lower stiffness, lower damping, and speed of the vehicle. And they are successively denoted by capital letters from A to L. All of these factors are listed in Table 4.

Considering computational efficiency, three levels are given for every factor. The value of parameters in fundamental models described above is thought as the second level. Then the first level and the third level are taken by 20% lower and 20% higher, respectively. For instance, three levels of the span length are 24 m (, ), 30 m (), and 36 m (, ). Of all these factors, the factors of width and speed have to be emphasized as follows. The simply supported girder bridge is assembled by many prefabricated smallbox section girders. Three levels of the width are 9 m (3 girders), 12 m (4 girders), and 15 m (5 girders). And the speeds of the moving vehicle are assumed as 20 m/s, 25 m/s, and 30 m/s, which are common in actual operation.
4.3. Numerical Simulations without Interaction
Firstly, the interaction has been ignored in this section. As for twelve factors and three levels for every factors, the orthogonal array (3^{13}) has been adopted. If only main effects (single factor) are considered, there is no mixture involved in the arrangement of factors. So they can be arranged successively. It is seen in Table 5.
 
1: 1st level, 2: 2nd level, 3: 3rd level, and Non: no factor in this column. 
All of these runs are analyzed by our own program VBCVA. The fundamental frequency () of the bridge is the most important index of dynamics. The dynamic load allowance (DLA) is the common index to take account of the dynamic performance of the bridge traversed by moving vehicles in design and evaluation. For simplification, the comfort of pedestrians on the bridge and passengers in the vehicle is evaluated by the vibration accelerations of the bridge () and the vehicle () in general. Therefore, these four indices are obtained from a large number of results, and they are listed in Table 6.

In most of design codes of highway bridges, the dynamic load allowance (DLA) is defined as the function of the fundamental frequency of the bridge. Taking the current code (JTG D602004) [47] in China as an example, the values of DLA defined in the code are compared with that calculated by VBCVA (Figure 10).
Figure 10 shows that the differences between them are much more significant, especially for the frequency ranging from 2 Hz to 4 Hz. Also, the DLAs obtained in different runs are not the same, even if the frequencies are almost the same. So the DLA is dependent on a lot of factors, other than just the fundamental frequency. At another point of view, it has been proved again that dynamic responses of the vehiclebridge system should be thought as the multifactor problem.
In statistics, two methods are always used for data processing, range analysis and variance analysis. The first method has advantage of low cost, simple thought, and convenience to be popularized and applied. But, compared with variance analysis, the range analysis has two defects. Firstly, the error is hardly estimated. Secondly, the reliability cannot be determined. And it also cannot be used for regression analysis and design. Of course, the method of various analyses can avoid both defects. Then these two methods are described as follows, respectively.
The range of analysis is completed by two steps. At first, the ranges of every factor are calculated by [48] where is the mean value of the th factor with the th level, and is the range of the th factor. Then the trend between factors and the test index is plotted. The results can be seen in Table 7 and Figure 11.

(a) Trend plot between factors and the DLA of the bridge.
(b) Trend plot between factors and the acceleration of the bridge.
(c) Trend plot between factors and the acceleration of the vehicle body.
Variance analysis is more rigorous. There are four steps to realize the variance analysis [48, 49].(i)Step 1: as for every factors, calculate the sum of square of deviations (), the degree of freedom (dof, ), and the variance estimation ().(ii)Step 2: estimate the variance of the error ().(iii)Step 3: obtain the test static and compare the with its critical value for given significance level .(iv)Step 4: for simplicity, the variance analysis table is listed, including the process and the results.
All of them are calculated as follows: in 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 is resulted from all of the vacant columns in the orthogonal array. Also, the accuracy increases with the increasing dof of the error. Therefore, if the significance level of one factor is larger than 0.25, it can be included as the error. Of course, that may not be the same in different situations. All of the results are listed in Table 8.
(a) Dynamic load allowance (DLA) of the bridge  
 
(b) Vibration accelerations of the bridge () and the vehicle ()  

It can be concluded that the influence of factors on different indices are listed as follows.(i)For DLA of the bridge, (ii)For vibration acceleration of the bridge, (iii)For vibration acceleration of the vehicle,
Obviously, for different indices, the most important influence factors are different. The sensitive factors are the pavement roughness, the span length, and the mass of the vehicle for the indices of the DLA, the vibration accelerations of the bridge and vehicle, respectively. It has to be noted that the DLA is the relative value of the dynamic response and the static response, while the vibration acceleration is the absolute value. It is the result of influence differences between the DLA and the vibration acceleration of the bridge, even though they are correlated in some extent.
4.4. Numerical Simulations considering 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. As a result, the orthogonal design is introduced to solve this problem.
Based on the results from the above section, the first five most important factors influencing the DLA are selected. They are the pavement roughness (), the upper stiffness of the vehicle (), the span length of the bridge (), the lower damping of the vehicle (), and the mass of the vehicle (). Meanwhile, due to the calculation cost and the existing orthogonal array, the number of levels is determined as two for each factor. The second level in the above section is taken as the first level in this section, which is based on the fundamental models. And the second level is defined as 20% higher than the first level. According to the sparsityofeffects principle, only the lower order interaction (two factors) is considered.
To avoid the mixture, the most important factor has to be arranged at first. It can be seen in Table 9.
 
1: 1st level, 2: 2nd level (be equal to the 2nd and the 3rd level in above section, resp.), and Non: no factor in this column. 
The methods are the same as that in the above section. The results of range analysis and variance analysis are listed in Tables 10 and 11, respectively.
