A Study on the Impact Coefficient of Traffic Flows on a Large Cable-Stayed Bridge in a Windy Environment

By developing a wind-vehicle-bridge coupling dynamic analysis system, the impact coefficient of traffic flows on a cable-stayed bridge and its probabilistic distribution are studied in the paper. Some important influence factors of the impact coefficient, such as the vehicle-bridge coupling effect, traveling velocity, and traffic density, are analyzed.The influence of the coupling effect between a bridge deck’s buffeting and vehicles’ vibration on the impact coefficient of traffic flows on a cable-stayed bridge is paid due attention to. The results show that vehicle-bridge coupling effect and surface roughness exert important influence on the impact effect of a cable-stayed bridge. Considering that impact coefficient of traffic flows on a bridge bears the characteristics of randomness, the paper, on the basis of a large number of samples, interprets the influence of traveling velocity and traffic density in view of its probabilistic distribution. In view of results derived from the coupling dynamic analysis, it is suggested that the impact coefficient be 25% for all spans of a long cable-stayed bridge like the Sutong Yangtze Bridge.


Introduction
The impact coefficient is an important problem in bridge design.Currently, the code provisions of impact coefficients in bridge design specifications could be classified into two categories: one gives the correlation between impact coefficients and bridges' lengths, while the other provides the correlation between impact coefficients and bridges' basic frequencies.Owing to the limitation of theoretical studies and calculating technology, many early studies on impact coefficients rely on field measurements.And on the basis of these measurements, code provisions were offered in design specifications.However, the measurements mainly concentrated on short or medium-length bridges.As we know, there are many factors such as driving velocity, bridge dynamic characteristics, and surface roughness, which will influence the impact effect of a bridge.A field measurement entails substantial human and financial resources.To study the impact coefficient problem of bridge components and so forth, bridge motion equations under the action of moving loads were developed into a vehicle-bridge coupling vibration analysis system.
The development of vehicle-bridge coupling vibration analysis technology has made it possible to study impact coefficients through numerical simulation.Such simulation could bring various influential factors into consideration and save both labor and financial resources.Previous studies [1][2][3][4] analyzed impact coefficients by using a vehicle-bridge coupling analysis system and considering the influence of surface roughness on impact coefficients.
But for a long cable-stayed bridge, some special features that do not exist in short-span bridges appear, such as simultaneous presence of various trucks and significant sensitivity to the wind.In recent years, there have been a number of studies on coupled wind-vehicle-bridge systems.For example, the studies on the coupled vibration between road vehicles and cable-stayed bridges in a cross wind environment could be found in the literature [5,6].In addition, vehicle-bridge coupling analysis under the action of a highway stochastic traffic flow was paid more attention to.Au et al. [7] studied the bridge motion under a traffic flow, in which axle loads whose statistical characteristics accord with the real situations were applied to a bridge model, but vehicle vibration and coupling 2 Shock and Vibration effect were neglected, and space intervals among vehicles were simplified into definite values.
A study on the impact coefficient of a large cablestayed bridge is conducted by a newly developed windvehicle-bridge coupling analysis system.Given that there are various types of vehicles in an actual traffic flow, the mathematic models of 4-axle, 5-axle, and 6-axle vehicles are deduced and incorporated into current wind-vehiclebridge coupling analysis system.Corresponding calculating programs are written in Fortran 90.With the Sutong Yangtze Bridge as an example, the impact effect of the bridge's three spans of different lengths is studied based on the developed wind-vehicle-bridge system.The differences between impact coefficients of traffic flows on the bridge and those of single vehicles on it are discussed.The coupling effect between vibrations of vehicles and the buffeting of the bridge is studied and impact coefficients of traffic loads on a cable-stayed bridge in a strong windy environment are paid due attention to.Impact coefficients from the coupling dynamical analysis in this paper are further compared with the results derived from bridge design specifications.

Vehicle-Bridge Coupling Vibration Analysis
System with a Traffic Flow and Wind Loads 2.1.Dynamical Model of a Bridge.The analytical model of a long-span cable-stayed bridge is established through the finite element method (FEM) using different kinds of elements such as beam elements and link elements in this study.
According to structural dynamics and finite element theory, the movement equation of a bridge model could be expressed as where M b , C b , K bT denote the mass matrix, the damping matrix, and the tangent stiffness matrix of a bridge, respectively.In this study, the material is homogeneous, isotropic, and linearly elastic.The Rayleigh damping matrix is adopted.F b () denotes equivalent node loads caused by external loads.
With regard to the wind-vehicle-bridge coupling system, the bridge subsystem's equivalent loads include where F bg , F vb , F stb , F bub , F seb denote self-weight, vehicle load, static wind load, buffeting force, and self-excited force, respectively.
A study on the impact coefficient of a cable-stayed bridge is conducted by taking the Sutong Yangtze Highway Bridge, a large cable-stayed bridge, as an example.The total length of the bridge is 8,206 meters (26,923 ft).Its main channel is a cable-stayed bridge with a span of 1,088 meters (3,570 ft) and had been with the longest cable-stayed bridge in the world until 2012.Its two side spans are 300 meters (980 ft) long each, and there are another four small cable spans (see Figure 1(a)).The bridge received the 2010 Outstanding Civil Engineering Achievement Award (OCEA) from the American Society of Civil Engineers.Two towers of the bridge are 300.4meters (986 ft) high and are the second tallest in the world until now.The tower is an inverted Y-shaped reinforced concrete structure with one connecting girder between the tower legs.From each tower, 34 stay cables radiate downwards in each of the four planes to support the bridge deck.The interval between every two cables is 16.2 m at the height of the bridge deck.The outer cable is up to 577 meters (1,885 ft) and is the longest in the world.The plane view size of foundation of pile group is 113.75 × 48.1 square meters and is the largest globally.The bridge deck is a steel box girder with transverse and longitudinal internal diaphragms and has fairing noses on both sides.The total width of the bridge deck is 41 meters including the fairing noses (see Figure 1(b)).Six standard highway lanes in two directions are designed, the traveling velocity of navigation channels is designed at 100 km/h, and that of nonnavigation channels is designed at 120 km/h.In Figure 1, the points A, B, and C denote the mid-spans of the 2nd, 3rd, and 4th spans.
A three-dimensional finite element model is established for the main channel of the Sutong Yangtze Bridge.Compared with the bridge length, the bridge width of 41 m is still very narrow.The existence of various stiffening ribs inside the steel box can ensure that the applicability of the plane crosssection assumption is applicable and a three-dimensional beam element model can embody the main vibration modes of the whole bridge.On the basis of the plane cross-section assumption, the main girder model adopted in this paper is the backbone model.That is to say, the bending, torsion stiffness, and the mass of the girder are all concentrated over the middle line consisting of three-dimensional beam elements.Rigid cross-girders with no mass are used to realize force transmission between the cables and the main girder.The two bridge towers are also discretized by using threedimensional beam elements.The sag effect of inclined cables exerts nonignorable influence on the cable stiffness, which is considered by equivalent elastic modulus in this paper.The simplification will not influence the dynamic characteristics of the whole bridge and thus will not influence the response of the bridge deck.Fixed constraints are adopted at the bottom of piers and towers.Vertical and lateral coupling constraints are used at the joints of the piers and the main beam.Since the bridge is a longitudinal floating structural system, only vertical coupling constraints are adopted at the joints of the towers and the main beam.The influence of the bridge deck pavement on dynamic characteristics is included by using lumped masses.The total finite element model contains 866 nodes, 1122 elements, 16 supporting conditions, 28 coupling constraints, and 636 sections.Some modes and corresponding natural frequencies of the bridge are given in Table 1.The 3rd-4th modes are the first two vertical bending modes that mainly present the main span's bending while the first two vertical modes embodying side spans' bending are the 47th-48th modes.Figure 2 shows the FE model of the bridge and some important vertical bending modes.In the dynamical analysis, the damping ratio of the bridge is assigned at 0.5%.the rigidity and damping of the upper suspension system of the vehicle in vertical and lateral directions, respectively;   ( = 1, . . ., 4) denotes the horizontal distance of the th axle from the center of gravity of the vehicle body;  1 denotes half of the transverse distance between tires; ℎ V denotes the vertical distance of the center of gravity of the vehicle body from the ground.The dynamics equation of the vehicle models could be deduced by the Lagrange principle, which is expressed in the following form:

Dynamical Models of
where u k is vector of vehicle displacement; M k , C k , K k are mass matrix, damping matrix, and stiffness matrix of the vehicle, respectively; F stk , F buk are static wind loads and  3, which are available in the literature [5].

Calculation Method of Vehicle-Bridge Coupling System.
It is assumed that the surface roughness of a bridge deck is not that rough to make a vehicle jump from the deck's surface and that there is no lateral sliding between the vehicle and the bridge deck; thus all tires of the vehicle always remain in contact with the bridge deck.The movement of a bridge deck produces additional roughness, which makes the vehicle vibration couple with the bridge deck's movement.A separate iteration method is used to solve the vehicle-bridge coupling system in this paper.Dynamical equations for each vehicle and the bridge are solved individually, and the consistent geometry relationship between the vehicles and the bridge is satisfied by repeatedly solving the two subsystems.The specific procedures are as follows.
(1) Combine surface roughness and the movements of the bridge deck at the last time step and calculate the forces acting on vehicles running on the bridge deck.
(2) Solve individually the vehicle subsystem using Newmark integration method to obtain the initial vehicle response ( ü  V , u  V ,   V ) at current time step.(3) Calculate the static wind forces, buffeting forces, and self-excited forces acting on the bridge deck according to the movement state of the bridge deck ( ) at the last time step, and solve individually the bridge subsystem and obtain the responses of the bridge ( ü   , u   ,    ) at the current time step.(4) The movements of the bridge deck obtained in step (3) are applied to step (1).Then steps (1) to (4) are repeated until the conditions of geometric compatibility at contact points are well satisfied and then the calculations for next time step continue.
The aforementioned analysis process is programmed using Fortran 90, in which two self-defined structural data types are used to describe dynamical models and motion state of vehicles, respectively.All the data of the bridge model, vehicle models, traffic flow, surface roughness, turbulent wind, and so forth are saved in text files and read by programs.
In order to verify the precision of the separate iteration method, which is used for the dynamical analysis of the vehicle-bridge coupling system in this paper, the responses of a simple-supported beam subjected to a moving mass and a moving sprung mass are calculated, respectively.The simple beam is 30 m long.Its modulus of elasticity and geometrical moment of inertia are 3.45 + 10 and 10 m 4 , respectively.The moving mass weighs 30 tons and stiffness of spring is 399 × 4 KN/m.The moving velocity of the mass or the sprung mass is assigned the value of 10 m/s.The damping of the beam and surface roughness on the beam are not included in this example.Biggs [9] gave the motion equations for the vibration of the coupling system of a simple beam subjected to a moving sprung mass on the basis of modal analysis method.
The matrix forms of the motion equations are as follows: where   () denotes the th generalized degree of freedom,   denotes the th frequency of the simple beam, and sin(V/) is its corresponding modal function;  denotes the number of the used modal functions;  V and  V denote sprung mass and unsprung mass, respectively;  V denotes stiffness of spring; V denotes moving velocity.When  V and  V are set to zero, the responses of the simple beam subjected to a moving mass could be obtained.When  V is set to zero, the responses subjected to a moving sprung mass in this example would be obtained.The vertical displacement and moment are calculated as follows: Figures 4(a)-4(d) give the responses of the simple beam subjected to a moving constant load, a moving mass, and a moving sprung mass, respectively.The figures show that the displacement and moment responses of the simple beam at mid-span obtained by the calculation programs in this paper are in accordance with those obtained by the method in the literature [9].The first 10 modals are included in the calculation.

Wind Forces on the Sutong Yangtze Bridge.
The bridge deck is a blunt body section, and the aerodynamic force it experiences is complicated.The three major parts of the wind force acting on a bridge deck, namely, the steady-state forces due to mean wind, buffeting forces due to wind turbulence, and wind-excited forces due to the aeroelastic interaction between bridge motion and fluent field, are included in the dynamical analysis in this paper.Generally speaking, the aerodynamic forces of cables are not comparative with those experienced by the bridge deck, and the cables' vibration exerts little influence on the bridge deck.Thus, like analyses of a bridge deck's vibration in most relevant literature, the aerodynamic forces of cables are ignored in this paper.The Sutong Yangtze Bridge is so famous that its aerodynamic force coefficients and flutter derivatives of the bridge deck are available in some published literatures.As we know, Reynolds number effect will influence aerodynamic force coefficients evidently.Here, the measured aerodynamic force coefficients of the Sutong Yangtze Bridge on the condition of high Reynolds number (1.0 × 10 6 ∼1.38 × 10 6 ) in the literature [10] are adopted.The flutter derivatives of the bridge come from the literature [11].

Simulation of Surface Roughness on the Bridge Deck.
As we know, surface roughness is generally regarded as a Gaussian stationary stochastic process.Simulation technology of surface roughness has been developed by some previous studies [12,13].The power spectral density (PSD) functions, which were developed by Dodds and Robson [12] to describe the road surface roughness, are as follows: where  is roughness coefficient (m  roughness exponents.There are three different classes of roads, and each class has either two, three, or four road conditions, namely, very good, good, average, and poor.The values of roughness coefficient  and exponents  1 and  2 are dependant on the class of road roughness, which were given in the literature [12].Here, the class of surface roughness is supposed to be "good" or "average" of "Principal roads." The roughness coefficient  equals 16 × 10 −6 or 64 × 10 −6 m 3 /cycle, respectively.The exponents  1 and  2 are set to 2 to simplify the problem. The minimum value of frequency for vehicle vibration is 0.5 Hz while the maximum is 30 Hz; the lower limit of traveling velocity is supposed to be 0.5 m/s and the upper limit is 35 m/s; the lower limit of space frequency of roughness is 0.5/35 = 0.014 cycle/m and the upper limit is 30/0.5 = 60 cycle/m.Longitudinal intervals are assigned at 0.1667 m.Considering that surface roughness on a road or a bridge deck varies along the width of the bridge, multipath roughness is simulated.Each path roughness represents a roughness at a point along the bridge width.The simulated roughness and their PSD (power spectral density) are shown in Figures 5 and  6.The simulated spectra are well accordant with their target spectra.A simulated roughness is periodic and the roughness only in one period is given in Figure 5. Since the bridge deck of the Sutong Yangtze Bridge is 2088 meters long, which is longer than the length of surface roughness in one period, the simulated roughness needs to be expanded according to its periodicity.

Impact Coefficient of a Large
Cable-Stayed Bridge The contacting forces of the two tires of the same axle are somewhat different because the multipath surface roughness on the bridge deck is adopted.
The ratio of the maximum dynamic response to the maximum static one is the dynamic magnification coefficient (1 + ) and  is the impact coefficient.Under the action of single vehicles, vertical displacements of the 1000-meterlong span are much larger than those of the 100-meter-long span.However, the former's impact coefficients are smaller than those of the latter.The surface roughness and vehicles' vibration make some contribution to impact coefficients of single vehicles on the 100 m span.Impact coefficients of different types of vehicles on the bridge are different.In view of the example, impact coefficients of 2-axle, 4-axle, and 6axle vehicles are relatively larger than others.The discussion of the influence of vehicle types and its regularity needs more bridge examples and should be paid more attention to when impact coefficients of ordinary beam bridges, which are more sensitive to the action of single vehicles, are studied.In general, with regard to the large cable-stayed bridge, impact coefficients of single vehicles are smaller than 5%.It could be deemed as safe enough to set the impact coefficient of a single vehicle on a cable-stayed bridge at 5%.But large cable-stayed bridges bear some features which short-and medium-span bridges do not.They are burdened with heavier traffic flows which consist of multiple vehicles.They have clustered frequencies and are sensitive to wind loads.The authors are more interested in whether the impact effect of a traffic flow on a cable-stayed bridge is evidently different from that of single vehicles on it.In a strong windy environment, vehicles' vibration excited by surface roughness couples with the bridge deck's buffeting caused by wind loads, which possibly exerts important influence on impact coefficients of traffic loads on a large cable-stayed bridge.

Impact Coefficients of Traffic Loads on the Bridge under the Action of Traffic
Flows.The information of a traffic flow is often obtained from on-the-spot traffic investigations, which consume lots of human and financial resources.Modern highways are all equipped with tolling stations and the detailed traffic flow data are recorded in the database there.In fact, these data provide a convenient access to the traffic flow information, which is adopted in the paper.Based on the data during a day from the tolling station of the Sutong Yangtze Bridge, traffic volumes within 24 hours are given in Figure 8, in which north (south) net data denote the number of vehicles from the north (south).The figure shows that traffic volume from north is evidently larger than that from south in this day.There are two peaks in hourly traffic volume curves.One occurs from 9:00 to 10:00 a.m. and the other from 3:00 to 5:00 p.m.The former is the rush hour in a day.According to the pass time of vehicles, the traffic flow data during the rush hour are shown in Figure 9, in which vertical coordinates 2-6 denote 2-6 axle vehicles.Of course, an actual traffic flow is complex, and the vehicles' traveling velocity and intervals Coupling dynamic response (I = 2.4%)

Static response
Vertical moment between them vary all the time.Figure 8 can only represent the traffic flow situation on one day and it is very possible that the traffic presents different characteristics on other days.However, the paper focuses on impact effect of traffic flows and wind loads on a cable-stayed bridge.It is unnecessary to pursue the traffic flow information of absolute precision.In the paper, the vehicles' traveling velocity is deemed as a constant.The influence of a traffic flow's intensity on impact coefficients will be discussed using the data during a day, which will be verified later.
This paragraph, by setting the traveling velocity at 100 km/h, calculates the dynamic and static responses of the bridge during the rush hour.The grade of surface roughness on the bridge deck is set as "good." The calculated results are partly shown in Figure 10.The figure shows that the impact effect of a 100-meter-long span under the action of a traffic flow is evidently larger than that under the action of a single vehicle and that the coupling effect exerts an evident influence on impact effect.The impact coefficients of vertical displacement and vertical moment responses are as large as 13% according to a part of the coupling dynamic responses given in Figure 10.However, the impact coefficients of the 1000-meter-long span are still smaller than 5% under the action of a traffic flow.
To understand the reason why the impact effect caused by a traffic flow is apparently larger than that caused by a single vehicle (a 4-axle vehicle), the power spectrum densities of responses are analyzed.Figure 11(a) shows that there are two evident peaks in the 100-meter-long span's response spectrum caused by either the traffic flow or the single vehicle.One is between 0.07 and 0.2 Hz; the other is around 1.0 Hz.Combining with Table 1, one can know that the former corresponds to low frequency vertical bending modes while the latter corresponds to high frequency ones.However, the power spectrum density of the response produced by a traffic flow is remarkably different from that produced by a single vehicle.The maximum value of the first peak band is about 10 times that of the second one in the 100 m long span's response spectrum produced by a single vehicle, while it is about 100 times in the spectrum caused by a traffic flow.It indicates that a traffic flow could excite more low frequency vertical bending modes, which accounts for the reason why a traffic flow can cause larger impact coefficients.Figure 11(b) shows that the 1000-meter-long span's response spectrum caused by the traffic flow presents a linearly declining trend and that there is no peak even around the basic frequency of the bridge (0.19 Hz) in the response spectrum, which indicates that there are little vibration ingredients in the 1000 m span's response.

Probabilistic Characteristics of Impact Coefficients of Traffic
Flows on the Bridge.In the aforementioned discussions, the impact effect of several traffic flow samples on the bridge is analyzed, which could illustrate the influence of vehiclebridge coupling effect and surface roughness on impact effect.However, owing to the randomness of impact coefficients, the limited number of samples is not convincing enough to account for the influence of some parameters like traveling velocity and so forth on impact coefficients.Here, the traffic flow during a whole day is used and the traffic flow in every 100 sec is taken as a sample.A total of 36×24 sample points of impact coefficients are obtained by coupling dynamic theory and then its probabilistic distribution is studied.
The traveling velocity is assigned at 100 km/h and the grade of roughness is set as "good." The impact coefficients of the displacement at mid-span of the 100-meter-long span are given in Figure 12, in which the probabilistic distributions of impact coefficients in different time segments are provided.Figures 12(a is the heaviest over the bridge.However, it is also proved reasonable to suppose that the probabilistic distribution of impact coefficients of a traffic flow on a bridge is independent on its intensity.So when enough data about traffic flows were not provided at a tolling station, one could randomly investigate the traffic flow data within 2-3 hours.Figure 13 presents the distributions of vertical displacement impact coefficients of the bridge under the conditions of different surface roughness.Figure 14 presents the probabilistic distributions of vertical moment impact coefficients at point A. The results further make it clear that vehiclebridge coupling effect and the grade of surface roughness exert remarkable influence on impact coefficients of traffic loads on the 100 m span of the cable-stayed bridge.When the grade of surface roughness is "good, " impact coefficients of traffic loads are generally no more than 20%.When the grade is "average, " the maximum impact coefficient of traffic loads on the 100 m span of the bridge is even up to 30%.However, the impact coefficients of traffic loads on the 1000-meterlong span are smaller than 5% in a windless environment.The reason is that the frequencies of the 1000 m span are much lower than those of vehicles' vibration.The 1000 m span presents noticeably different characteristics in a strong windy environment, which will be analyzed in detail later in this paper.
Figure 15 presents the distributions of impact coefficients of traffic flows on the 100 m span at different traveling velocities.Although the impact coefficient bears the feature of randomness, its probabilistic distribution illustrates that traveling velocity exerts regular influence on the impact effect.With traveling velocity rising, both the mean values and standard deviations of impact coefficients increase.

Impact Coefficients of Traffic Loads on the Bridge in a
Strong Windy Environment.To calculate the time history responses of the bridge under the action wind loads, the stochastic wind velocity field needs to be simulated.The spectrum representation method developed by Deodatis [14] is adopted here.By referring to the literature [15], the lateral wind spectrum adopted is in Kainmal's form; the vertical spectrum adopted is in the form presented by Lumley and Panofsky.The coherence function adopted is in Davenport's form in which the correlation coefficient is set to 7.0.The dynamic responses of the bridge under the action of wind loads at a mean wind velocity of 49.8 m/s, which is the design wind velocity for the Sutong Yangtze Bridge, are analyzed.The results are partly shown in Figure 16.The figure shows that both the vertical and lateral responses at point C under the action of wind loads fluctuate around their mean values, respectively.However, their mean values are not zero and the nonzero mean values are caused by static wind loads.The vertical static response caused by wind loads is smaller than the lateral one.Yet the amplitude of the vertical dynamic response caused by turbulent wind is much larger than that of the lateral one.The reason is that the lift item of the bridge deck's aerodynamic coefficients is more sensitive to wind attack angle than the drag and torque items and the derivative of the lift coefficient with respect to the wind attack angle is larger than those of the other two.The bridge deck's selfexcited force, which is caused by the mechanism of fluid-solid coupling, makes the vertical buffeting decrease evidently.It indicates that fluid-solid coupling restraints the buffeting of the bridge deck at the design wind velocity.
By using one sample of the traffic data in 10 min during the above-mentioned rush hour, the responses produced by wind loads and caused by the traffic flow are compared.Figure 17 shows that when the mean wind velocity is 20 m/s, the main span's buffeting responses caused by turbulent wind are comparative with that caused by the traffic flow.When the mean wind velocity increases to 30 m/s or 49.8 m/s, the main span's buffeting responses caused by turbulent wind become much larger, even larger than that caused by the traffic flow.In these conditions, the wind loads become one of the main loads.Figure 17(c) offers the contacting force responses between a 4-axle vehicle and the bridge deck in the strong windy environment.Figure 17(d) shows roughness and additional roughness, which is caused by the bridge deck's motion, of a contacting point of one tire.All contacting responses of the tires are positive, which presents the correctness of assumption that all tires of vehicles always keep in touch with the bridge deck.Moreover, it could be observed that the frequencies of the 1000-meter-long span's buffeting are much higher than those of its response under the action of a traffic flow.So there exists a coupling phenomenon between the bridge deck buffeting and vehicle vibrations.The wind-vehicle-bridge coupling vibration is further analyzed.To analyze the problem clearly, the buffeting response point C caused by wind loads is subtracted from the dynamic response produced by the joint action of wind loads and the traffic flow.The subtraction is also the response produced by the traffic flow, in which the coupling between the bridge buffeting and vehicle oscillation is taken into consideration.Figure 18 compares the bridge deck's response caused by the traffic flow in a windy environment and that in a windless environment.When the mean wind velocity is 20 m/s, there is no obvious difference between the results from the wind-vehicle-bridge coupling analysis and the vehicle-bridge coupling one.However, when the mean wind velocity rises to 30 m/s, the maximum impact coefficient obtained from the wind-vehicle-bridge analysis is 13%, much larger than that from the vehicle-bridge one.When the mean wind velocity is 49.8 m/s, the maximum (2) Impact coefficients bear the feature of randomness, so the probabilistic distribution of impact coefficients, which are based on large numbers of samples, could interpret the influence of various kinds of factors more clearly.It is presented that standard deviations of impact coefficients of heavy traffic flows are a little larger than those of sparse traffic flows while their mean values are almost the same.So it is reasonable to suppose that the probabilistic distribution of traffic flows' impact coefficients on a bridge is independent on the traffic intensity.With traveling velocity rising, both mean values and standard deviations of traffic flows' impact coefficients on a bridge increase.Vehicle-bridge coupling effect and the grade of surface roughness exert important influence on impact coefficients of traffic loads on a cable-stayed bridge.When the grade of surface roughness is "good, " the impact coefficients are no more than 20%.When the grade is "average, " the maximum impact coefficient on the 100 m span is even up to 30%.However, the impact coefficients of traffic loads on the 1000-meter-long span are smaller than 5% in a windless environment.
(3) The cable-stayed bridge's static vertical response produced by wind loads is smaller than the lateral one.But the amplitude of the bridge's vertical buffeting response caused by turbulent wind is much larger than that of the lateral one.When the mean wind velocity is more than 30 m/s, the 1000-meter-long span's buffeting response is much larger than its response caused by a traffic flow.And there exists a coupling phenomenon between the vibration of vehicles and the buffeting of the bridge deck, which makes that impact coefficients of traffic loads on the 1000 m span of the cable-stayed bridge become as large as 25% in a strong windy environment.Grade of surface roughness exerts little influence on impact coefficients of traffic loads on the 1000 m span and thus the vehicle-bridge coupling effect becomes a decisive factor for impact effect of traffic loads on a cablestayed bridge in a strong windy environment.
(4) Impact coefficients calculated by the coupling dynamic analysis system are compared with the results derived from bridge design specifications.According to results obtained from coupling dynamic analysis in the paper, a preliminary conclusion is drawn that the impact coefficient should be 25% for a large cable-stayed bridge.Impact effect of a large cable-stayed involves many factors.Some important factors and relevant regularities are discussed and some conclusions are drawn in the paper.However, only one cable-stayed bridge is taken as an example.More examples should be studied before rational suggestions about impact coefficients of traffic loads on a cable-stayed bridge are provided to bridge designers.The problems analyzed in the paper should also be paid due attention to when the impact effect of other types of bridges is studied.

Figure 2 :
Figure 2: FE model of the bridge and its vertical bending modes: (a) FE model; (b) the 3rd mode; (c) the 4th mode; (d) the 47th mode; (e) the 48th mode.

Figure 4 :
Figure 4: Dynamic response of a simply supported beam at mid-span: (a) vertical displacement under the action of a moving constant load; (b) vertical displacement under the action of a moving mass; (c) vertical displacement under the action of a moving sprung mass; (d) vertical moment under the action of a moving sprung mass.

3. 1 .
Vehicle-Bridge Coupling Vibration Analysis.Suppose a vehicle runs on the inner lane and its traveling velocity is 100 km/h.The grade of surface roughness on the bridge deck is supposed to be "average." Time histories of vehiclebridge coupling dynamic responses of the Sutong Yangtze Bridge under the action of single vehicles are given in Figures7(a)-7(j), in which the static responses are obtained by neglecting all dynamic effect and surface roughness and dynamic responses are obtained by regarding all axles as moving constant loads.Point A represents the mid-span of a 100-meter-long span and point C represents the midspan of the 1088-meter-long span, which will be referred to as the 1000 m span in the following context.Every axle is supported by two tires on the bridge deck and the contacting force of every tire fluctuates around half of its axle load.

2 Figure 8 :
Figure 8: Traffic flow in every hour during a day.

Figure 10 :Figure 11 :
Figure 10: Time histories of responses of the bridge deck under the action of a traffic flow: (a) vertical displacement at point A; (b) vertical displacement at point C; (c) vertical moment at point A; (d) vertical moment at point C.

Figure 14 :
Figure 14: Probabilistic distribution of vertical moment impact coefficients at point A of the bridge under the conditions of different surface roughness: (a) good; (b) average.

Figure 15 :
Figure 15: Probabilistic distributions of vertical displacement impact coefficients at point A of the bridge at different traveling velocities: (a) 80 km/h; (b) 120 km/h; (c) 80 km/h; (d) 120 km/h.
Figure 19: Probabilistic distributions of vertical displacement impact coefficients at point C of the bridge in a windy environment: at the mean wind velocity of (a) 30 m/s; (b) 49.8 m/s; (c) 30 m/s; (d) 49.8 m/s.

Figure 20 :
Figure 20: Probabilistic distributions of vertical moment impact coefficients at point C of the bridge in a windy environment: at the mean wind velocity of (a) 30 m/s; (b) 49.8 m/s.
[8]tiaxle Road Vehicles.According to the literature[8], nationwide surveys have been conducted in China and data of almost 70,000 vehicles have been

Table 2 :
Typical style for vehicle loads.

F bk is tire-deck
interactions which are induced by deck roughness and deck motion.The dynamical parameters of vehicle models, such as stiffness and damping, are given in Table