Dynamic Responses of Simply Supported Girder Bridges to Moving Vehicular Loads Based on Mathematical Methods

For dynamic responses of highway bridges to moving vehicles, most of studies focused on single-factor 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 somenumerical simulations based on the proposedmethod 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 themodal synthesis method, the programVBCVAhas 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 goodway to obtainmore rational empirical formulas of the DLA and other dynamic responses, which may be adopted in the revision of codes for design and evaluation.


Introduction
Research on the dynamic analysis of the vehicle-bridge 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 beam-type 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 simple-span [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 (H20-44 truck and HS20-44 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 three-dimensional 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 vehicle-bridge system [14][15][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 vehicle-bridge 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 vehicle-bridge 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  = 1 +  =  ,max  ,max (1) in which the  ,max and  ,max 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][19][20][21][22].

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.

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 vehicle-bridge 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 Dirac-delta 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:

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 vehicle-bridge 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  V 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 by where   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 D60-2004) 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  V , , , 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 1/ 2 V in the mass matrix and the term  V  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 as where  and  are only dependent on the values of , , and , other than the frequency  V .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).

Program for Vehicle-Bridge Coupled Vibration Analysis
Based on their own advantages of existing generalized commercial software ANSYS and MATLAB, the program VBCVA (Vehicle-Bridge Coupled Vibration Analysis) has been developed by our own research group.It is used for vibration and dynamic analysis of the vehicle-bridge 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.

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.

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 2 + 3, including the vertical dof  1 , pitching dof  2 , swung dof  3 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 centre-ofgravity 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.
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 : ] . (26d)

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.
where the  denotes the number of selected mode shapes and V, , 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.

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 time-varying system of differential equations.And it is difficult to obtain the closed-form solutions.But they can be solved by some numerical methods, such as the central difference method, Newmark method, Wilson- method, and so on [33][34][35].
Wilson- method has been adopted in this study.It is a self-stabilization 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 are To 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  = 1 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 well-researched 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:  () =  −   <  <   (39) 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: =   + ( − 0.5) Δ  = 1, 2, . . ., , (41b) 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.
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)).

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 vehicle-bridge coupled vibration problem.

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.

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 sparsity-of-effects 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 sparsity-of-effects principle states that a system is usually dominated by main effects (single factor) and low-order interactions.Therefore, it is most likely that main effects and two-factor interactions are the most significant responses in an experiment [45].
A majority of studies have shown that the dynamic responses of vehicle-bridge 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.

Samples Selection.
As we know, there are many types of highway bridges, including the girder bridge, the arch bridge, the cable-stayed 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 m-span bridge with small-box section has been selected as the fundamental sample of bridges.The damping ratio is assumed as 0.05.The cross-section is in Figure 8 and its other parameters are listed in Table 2.
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 vehicle-bridge 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 3-axle 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 vehicle-bridge 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 ( = 0.2345) 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 ( 1 , 30 m × 0.8), 30 m ( 2 ), and 36 m ( 3 , 30 m × 1.2).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 small-box 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.

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  27 (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.
All of these runs are analyzed by our own program VBCVA.The fundamental frequency ( 1 ) 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 ( V ) 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 D60-2004) [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  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]   = max ( 1 ,  2 , . . .,   ) − min ( 1 ,  2 , . . .,   ) , (42) 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.
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 (  ).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.
It can be concluded that the influence of factors on different indices are listed as follows.(iii) For vibration acceleration of the vehicle,  >  >  >  >  >  >  >  >  >  >  > .(46) 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.

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 sparsity-of-effects 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.
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.
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, apart from the lower damping of the vehicle (), other four main factors are still the most important influences on the test indices, including the DLA of the bridge, the vibration accelerations of the bridge and the vehicle.And the interaction effects are less.Maybe they can be neglected in engineering application.In addition, there are some differences for the importance order of factors between the results from this section and the above section.However, general trend has been little changed.Of course, the later one is more accurate and reliable, due to considering the interaction other than just main effects.
Similarly, for different testing indices, the influence of each factor is not the same.

Results and Discussions.
According to numerical simulations, the sensitivity of each factor for different indices is not the same.Therefore, we should pay more attention to various factors for different purpose in the design of bridges and vehicles.For example, if the bridge is mainly for pedestrians, the acceleration of the bridge is related with their comfort.Then the span length and the pavement roughness should be noted.If the riding comfort is the primary consideration, maybe we should focus on the pavement roughness and the mass of the moving vehicle.However, as for safety of the bridge, the important factors include the pavement roughness, span length of the bridge, mass, and upper stiffness of the vehicle.In short, it has to be noted that the pavement roughness is the essential factor.
Comparing the acceleration of the bridge with that of the vehicle, the former one is more related with the span length, while the latter one is more related with the mass of the vehicle.Actually, the span length of the bridge is significantly related with its fundamental frequency.And the mass of the vehicle will influence natural frequencies of the vehicle.Therefore, in nature, the determined factors are natural frequencies of the bridge and moving vehicles.Also, it is consistent with the conclusion obtained from the simplified theoretical derivation in Section 2. Based on this consistency, the proposed method, orthogonal experimental design, has been verified.And the results are reliable.
For vehicle design, stiffness is more important than the damping of the suspension system, especially the upper stiffness.According to the structural dynamics, the damping of the vehicle has less influence on its natural frequencies than the stiffness.Meanwhile, it has been proved that the frequency is the main factor.
Using the method of orthogonal experimental design, the interaction effects are so small that they can be neglected in engineering application.
Furthermore, the plot of relation between DLA and the fundamental frequency shows that the DLA is dependent on many factors, and more synthetical and accurate formulas should be proposed in the code.In addition, for comfort of pedestrians and passengers, the empirical formulas of accelerations of the bridge and the vehicle could be added based on the proposed method in this study.

Conclusions
To get a good understanding of the dynamic responses of the simply supported girder bridge to moving vehicular loads or vehicles, theoretical derivations are presented based on simplified models of the vehicle-bridge coupled vibration system; then some numerical simulations based on the proposed method of the orthogonal experimental design in batches have been done by our own program VBCVA.
Firstly, according to simplified theoretical derivations, three factors are proved as the most important factors to determine dynamic responses.They are the mass ratio of the vehicle and the bridge (), the frequency ratio of the vehicle and the bridge (), and the ratio of disturbance frequency related with the speed and the frequency of the bridge ().It has to be noted that the roughness has not been considered in simplified models.
Secondly, based on the modal synthesis method, the thought of our own program VBCVA has been given in detail.Combining  cases of multilanes and multivehicles (both in longitudinal direction and transverse direction).Also, the number of axles or vehicles is not limited.Therefore, this program has to be thought of as convenient and powerful enough for the analysis of vehicle-bridge coupled vibration problems.Thirdly, on the basis of the orthogonal experimental design, both the main effects and the interaction effects can be studied.For different indices of dynamic responses, the influences of each factor are not the same.However, the pavement roughness is assumed as the most important factor.In addition, dynamic responses of the bridge and the vehicle are much more dependent on their own fundamental frequency, respectively.Furthermore, the interaction effects have proved to be so small that they can be neglected in engineering application.
In the end, it has to be emphasized that the proposed method of the orthogonal experimental design in batches greatly reduces calculation cost.And it is efficient and rational enough to study multifactor problems.
The proposed method is used for not only the analysis of influence factors, but also the analysis of regression.And it can be applied in all types of bridges, other than just the simply supported bridge.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.

Figure 1 :
Figure 1: A single moving constant load.

Figure 4 :
Figure 4: Schematic plot of the vehicle model.
Maximum amplitude h (cm)

Figure 8 :
Figure 8: The cross-section of the bridge (cm).

Figure 10 :
Figure 10: Relation between the DLA and the fundamental frequency.

Figure 11 :
Figure 11: Trend plot between factors and the test index.

Table 2 :
Parameters of the bridge.: width, : moment of inertia, : modules of elasticity, and : density.
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 vehicle-bridge system should be thought as the multifactor problem.

Table 6 :
Results of numerical simulations without interaction.

Table 7 :
Range analysis of numerical simulations without interaction.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.

Table 8 :
Various analyses of numerical simulations without interaction.Vibration accelerations of the bridge (  ) and the vehicle ( V )

Table 9 :
Orthogonal array  16 (215) and the arrangement of factors.

Table 10 :
Range analysis of numerical simulations considering interaction.
both advantages of the commercial software ANSYS and MATLAB, the bridge model is found by ANSYS, while the vehicle model and the roughness model are established by MATLAB.It can be used for calculating the (c) Continued.