Passive Control System for Mitigation of Longitudinal Buffeting Responses of a Six-Tower Cable-Stayed Bridge

This paper presents an investigation of mitigation of longitudinal buffeting responses of the Jiashao Bridge, the longest multispan cable-stayed bridge in the world. A time-domain procedure for analyzing buffeting responses of the bridge is implemented in ANSYS with the aeroelastic effect included. The characteristics of longitudinal buffeting responses of the six-tower cable-stayed bridge are studied in some detail, focusing on the effects of insufficient longitudinal stiffness of central towers and partially longitudinal constraints between the bridge deck and part of bridge towers. The effectiveness of viscous fluid dampers on the mitigation of longitudinal buffeting responses of the bridge is further investigated and amultiobjective optimization designmethod that uses a nondominating sort genetic algorithm II (NSGA-II) is used to optimize parameters of the viscous fluid dampers. The results of the parametric investigations show that, by appropriate use of viscous fluid dampers, the top displacements of central towers and base forces of bridge towers longitudinally restricted with the bridge deck can be reduced significantly, with hampering the significant gain achieved in the base forces of bridge towers longitudinally unrestricted with the bridge deck. And the optimized parameters for the viscous fluid dampers can be determined from Pareto-optimal fronts using the NSGA-II that can satisfy the desired performance requirements.


Introduction
For long-span cable-stayed bridges, the multispan cablestayed bridges with three or more towers have been a recent design trend [1,2].Typical examples of this bridge type are the Millau Viaduct Bridge in France, the Maracaibo Bridge in Venezuela, the Rion-Antirion Bridge in Greece, the Mezcala Bridge in Mexico, the Dongting Lake Bridge in China, and the Ting Kau Bridge in Hong Kong [1][2][3][4].The longest multispan cable-stayed bridge in the world is Jiashao Bridge in China, which is a six-tower cable-stayed bridge with the total length of 2680 m.
Compared with a conventional three-span cable-stayed bridge with two towers, there are two major problems in the design of multispan cable-stayed bridges.One is the insufficient longitudinal stiffness that arose from the central tower(s).In the conventional cable-stayed bridges with two towers, each of the towers is connected through outermost stay cables to the fixed anchorage or anchor pier, which can provide effective support to the towers [2].However, in a multispan cable-stayed bridge, the beneficial effect of the fixed anchorage or anchor pier diminishes for the central tower(s), which results in the insufficient longitudinal stiffness of central towers.Hence, the structural responses of a multispan cable-stayed bridge such as deflections of main girders are increased significantly under unbalanced live loadings compared with a conventional cable-stayed bridge.Another important problem in the design of multispan cablestayed bridges is the large temperature deformation in the long bridge deck.The commonly used structural measure to reduce the temperature effects is the application of partially longitudinal constraints between the bridge deck and part of bridge towers such as the Millau Viaduct Bridge in France and the Jiashao Bridge in China.The application of the partially longitudinal constraint system can improve the static performance of the multitower cable-stayed bridge under temperature action.
To date, an amount of research work had been done on investigating the static performance of the multispan cablestayed bridge under living loadings and temperature action [1,3,4].However, studies of the wind-resistant performance of the multispan cable-stayed bridge considering the insufficient longitudinal stiffness of central towers and partially longitudinal constraint system are meager.It is especially desirable to understand the effects of the insufficient longitudinal stiffness of central towers partially longitudinal constraint system on the characteristics of longitudinal buffeting responses of the bridge.And it is also necessary to provide a mitigation scheme for the longitudinal buffeting responses of the multispan cable-stayed bridge.Based on the aforementioned motivation, this paper investigates the performance of viscous fluid dampers for mitigation of longitudinal buffeting response of the Jiashao Bridge, the longest multispan cablestayed bridge in the world.The specific objectives of this study are to (i) investigate the characteristics of longitudinal buffeting responses of the multispan cable-stayed bridge, focusing on the effects of insufficient longitudinal stiffness of central towers and partially longitudinal constraints between the bridge deck and part of bridge towers; (ii) investigate the effects using viscous fluid dampers for mitigation of longitudinal buffeting response of the multispan cable-stayed bridge; and (iii) arrive at the values of optimal design parameters of the viscous fluid dampers for mitigation of longitudinal buffeting responses of the bridge using a multiobjective optimization design method.

Bridge Description
The subject of this study is Jiashao Bridge shown in Figure 1(a), which is a six-tower cable-stayed bridge that crosses the Hangzhou Bay, along the highway between Jiaxing and Shaoxing in China.The total length of the bridge is 2680 m with the span arrangement of 70 m + 200 m + 5 × 428 m + 200 m + 70 m, which is the longest multispan cablestayed bridge in the world.Figure 1(b) shows the schematic elevation view of the Jiashao Bridge.In order to overcome the problem of large temperature deformation in the long bridge deck, the partially longitudinal constraint system is applied in the design of Jiashao Bridge.As shown in Figure 1(b), the longitudinal constraints are applied to restrict the bridge deck from moving in the longitudinal direction at bridge towers number 2 and number 5, respectively.And there are no longitudinal constraints between the bridge deck and other bridge towers.The role of such partially longitudinal constraint system is to reduce the temperature-induced longitudinal deformation in the bridge deck.Figure 1(c) shows the typical deck cross section of the bridge.

Dynamic Finite Element Model
A three-dimensional finite element model of the Jiashao Bridge has been developed by use of the commercial software package ANSYS [5].The finite element model involves 1402 nodes and 1872 elements, as shown in Figure 2. In this model, a double-girder model is used to simulate the bridge deck system when conducting dynamic analysis.The two steel box girders are modeled as Timoshenko's beam elements with 6 degrees of freedom (DOFs) at each node, which account for transverse shear deformation, biaxial bending, and axial strain.And the cross-beams connecting two steel box girders are modeled as rigid links at 30 m intervals as shown in Figure 3(a).The bridge towers and piers are also modeled as Timoshenko's beam elements.A 2-node truss element is used to simulate the stay cables, which accounts for only tension and no compression based on the real condition.Considering the geometric stiffness of stay cables under dead loading, the Ernst equivalent elastic modulus for stay cables is adopted [6].Rigid links at 15 m intervals are used to connect the cables in four planes to the girders as shown in Figure 3(b).Figure 3(c) further shows the case of rigid links simultaneously connecting two steel box girders and four stay cables.
The detailed finite element modeling for the constraints between the girders and towers (piers) are described as follows.
(i) Modeling of constraints between the girders and towers: As shown in Figure 4(a), constraints are applied to restrict the deck from moving in the longitudinal, lateral, and vertical directions at bridge towers number 2 and number 5. Thus, the constraint equations of DOFs for the nodes N1∼N3 shown in Figure 4(a) can be defined as It should be noted that, for other towers number 1, number 3, number 4, and number 6, constraints are applied to restrict the deck from moving in the lateral and vertical directions.Thus, the constraint equations of DOFs for the nodes N1∼N3 can be defined as (ii) Modeling of constraints between the girders and piers in each side span: As shown in Figure 4(b), constraints are applied to restrict the motion of bridge deck at all bridge piers including auxiliary pier and end pier to allow only longitudinal displacement  and rotations about the  and  axes.Thus, the constraint equations of DOFs for the nodes N1∼N6 shown in Figure 4(b) can be defined as The modal analysis of the Jiashao Bridge is essential to determine the buffeting responses in succeeding dynamic analysis subjected to wind loadings.It should be noted that  the modal analysis of bridge is followed starting from the static equilibrium state due to dead load and cable pretension, which is obtained by geometrically nonlinear analysis of the bridge.The LANCZOS eigenvalue solver is then adopted for modal analysis.Main vibrations modes of finite element model are listed in Table 1 and part of them are shown in Figure 5.

Longitudinal Buffeting
Responses of the Bridge where the design wind speed 39.3 m/s at the height of 10 m has a return period of 100 years.As for the design wind spectra, the autospectra adopted for longitudinal and vertical components are Kaimal spectrum and Panofsky spectrum, respectively.Kaimal spectrum: Panofsky spectrum: where   and   are spectral-density function of longitudinal velocity fluctuation and vertical velocity fluctuation, respectively;  is fluctuation frequency;  = /() is the Monin coordinate;  * is the friction velocity.
The following empirical exponential coherence function is adopted for the velocity cross-spectrum [7]   () where   and   are coordinates in the and -directions, respectively;   and   are 8.3 and 6.4 for longitudinal turbulence components and 3.8 and 4.8 for vertical turbulence components.It should be noted that the coherence between the longitudinal and vertical components is disregarded in the present study.

Wind-Field Simulation at Bridge Deck and Towers.
The wind velocity at each simulation point of bridge deck and towers is the sum of the mean wind velocity and turbulent fluctuation around the mean wind velocity.In the present study, the fluctuating wind velocity is taken as realization of a multivariate Gaussian vector process having zero mean and known fixed spectral densities at all simulation points [8].Thus, the longitudinal and vertical fluctuating velocities in the overall wind field are modeled as 191-variate vector process, u(), and a 95-variate vector process, w(), respectively.By using the improved Deodatis simulation method [8,9], the longitudinal and vertical fluctuating wind-velocity samples at 95 points on the bridge deck and 96 points on the bridge towers are simulated.The parameters used in the simulation are as follows: the cut-off frequency   = 5rad/s, the number of fast-Fourier transform  = 1024, and time interval Δ = 0.2 s.Figures 6 and 7 illustrate the autospectra and correlation functions of the longitudinal fluctuating velocities for two points, Point 1 and Point 15, at the bridge tower number 3 with a distance of 150 m.As can be seen from the figure, the simulated autospectra and correlation functions are in well agreement with the targets.In this section, the buffeting analysis method is introduced briefly.Consider a deck section in a turbulent flow and in motion at an instant time .Based on the quasi-steady theory, the aerodynamic forces per unit span in the global axes are composed of three wind loading components, namely, the steady aerodynamic forces, the buffeting forces due to the fluctuating wind components  and , and the selfexcited forces induced by interaction between bridge and wind motions [10][11][12].The steady aerodynamic forces are expressed as

Longitudinal Buffeting
where  st  ,  st  , and  st  are three components in the horizontal, vertical, and rotational directions;  is the air density and  = 1.225 kg/m 3 in this study; (= 55.6 m) is the deck width;  is the mean wind velocity;   ,   , and   are the drag force, lift force, and pitching moment coefficients.The aerodynamic force coefficients are measured through the wind-tunnel tests.In this study   ,   , and   are taken as 0.0943, −0.1627, and 0.0088, respectively, at 0 ∘ wind attack angle.The self-excited forces are expressed as where  se  ,  se  , and  se  are three components in the horizontal, vertical, and rotational directions;  = / is the reduced circular frequency of motion and  is the circular frequency;   ,   , and   are the horizontal, vertical, and torsional displacements, respectively;  *  ,  *  , and  *  ( = 1, 2, . . ., 6) are the dimensionless flutter derivatives, which are characteristic cross-sectional properties expressed in terms of the reduced wind speed Ũ = /() and  is the natural frequency.The flutter derivatives at 0 ∘ wind attack angle under different reduced wind speed Ũ are obtained through the wind-tunnel tests and listed in the literature [13].
The buffeting forces are expressed as where    ,    , and    are the derivates of   ,   , and   with respect to wind attack angle, respectively; () and () are the longitudinal and vertical turbulence wind components, respectively;   ,   ,   ,   ,   , and   are the six complex aerodynamic admittance functions which are functions of reduced frequency.In this study, the    ,    , and    are taken as 0.1934, 5.5119, and 0.8337, respectively, at 0 ∘ wind attack angle.It should be noted that the effect of the aerodynamic admittance is disregarded because it was not available from wind-tunnel tests.Thus, the values of six aerodynamic admittance functions are all taken as unity conservatively.
Based on (5)∼(7c), investigations on time-domain buffeting analysis have been done by several researchers in recent years [11,12,14].In the present study, the buffeting analysis is carried out in the software package ANSYS with the following steps [12,15]: (1) A 3D finite element model of the bridge is established using the software package ANSYS.(2) The selfexcited aerodynamic forces acting on unit span of the bridge deck according to (6a)-(6c) can be modeled as elemental aeroelastic stiffness and aeroelastic damping matrices by the user-defined Matrix27 element in ANSYS [11,12,14].The detailed procedure to calculate the self-excited aerodynamic forces using user-defined Matrix27 element in ANSYS is demonstrated in detail by Wang et al. [12].(3) The steady aerodynamic forces and buffeting forces are calculated by ( 5), (7a), (7b), and (7c), respectively.(4) The self-excited aerodynamic forces, buffeting forces, and steady aerodynamic forces are applied as external loadings to the structural model to analyze the buffeting responses of the bridge in time domain [11,12].

Longitudinal Buffeting Responses of the Bridge.
The buffeting response analysis of the Jiaoshao Bridge was carried out at the design wind velocity of 49.72 m/s at the bridge deck for the return period of 100 years.The longitudinal buffeting response quantities selected are listed in Table 3.It should be noted that considering the dynamic coupling between vertical bending of the bridge deck and longitudinal bending of bridge tower as shown in Table 1 the vertical   displacement of the bridge deck is also selected as longitudinal buffeting response quantity.Figure 8 shows the timehistories of vertical buffeting displacement in the middle of the girder span between bridge towers number 3 and number 4. The root-mean-square (RMS) buffeting responses are further calculated to evaluate the longitudinal windresistance performance of the multispan cable-stayed bridge.
Figure 9 shows the RMS longitudinal buffeting responses of bridge towers.Table 4 shows RMS vertical displacements of bridge deck.
From Figure 9 and Table 3 the following can be seen.(i) The RMS base shears and RMS base moments of the bridge towers number 2 and number 5 in the partially longitudinal constraint system are approximately 3.6 times as those in the other towers.Thus, part of bridge towers longitudinally restricted with the bridge deck results in the concentration effects of base shear and moment.(ii) The RMS top displacements of central towers (towers number 2, number 3, number 4, and number 5) are approximately 1.7 times as those of side towers (tower number 1 and number 6).This is because the beneficial effect of the anchor piers diminishes for the central towers and hence the longitudinal stiffness of the central towers is smaller than that of the side towers [2].(iii) The top displacements of bridge towers are highly correlated with vertical displacements of bridge deck.Figure 10 further shows the power spectrum of longitudinal buffeting displacement on the top of bridge tower number 3 and power spectrum of vertical buffeting displacement in the middle of the bridge deck.It can be seen that the first and second peaks in both two power spectra correspond to the 1st and 9th vibration mode, which are all coupling modes of vertical bending of the bridge deck and longitudinal bending of bridge tower.

Parametric Study on Viscous Fluid Dampers.
Passive dampers have attracted the particular attention of structural control engineers due to their reliability and ease of implementation.Previous investigations have demonstrated the effectiveness of passive dampers in reducing the seismic responses of long-span cable-stayed bridges [17][18][19].And Domaneschi and Martinelli [20,21] provide a passive control method for mitigation of transverse buffeting responses of a long-span suspension bridge.Therefore, it would be interesting to investigate the application of passive dampers such as viscous fluid dampers in reducing the longitudinal buffeting responses of the multitower cable-stayed bridge.In this section, a parametric study is performed to investigate the effects of viscous fluid dampers on the mitigation of longitudinal buffeting responses of the multispan cablestayed bridge.In a viscous fluid damper, the damping force  is described as [22]  =  |V|  sgn (V) , where  is the damping coefficient of the damper;  is the velocity exponent of the damper; v is the relative velocity across the damper.From (8), the damping coefficient and the velocity exponent are the two parameters that govern the selection of a viscous fluid damper.In the present study, the velocity exponent  is varied from 0.1 to 1.0 with the interval of 0.1.And the damping coefficient  is varied from 2000 to 12000 kN⋅(s/m)  with the interval of 1000 kN⋅(s/m)  .Considering the configuration of the multispan cable-stayed bridge with the partially longitudinal constraint system, the viscous fluid dampers are placed at the bridge towers longitudinally unrestricted with the bridge deck, that is, towers number 1, number 3, number 4, and number 6 in Jiashao Bridge as shown in Figure 11.The buffeting responses of the bridge using the viscous fluid dampers are calculated.For instance, Figure 12 shows the comparisons of time-histories of vertical buffeting displacement in the middle of the girder span between bridge towers number 3 and number 4 with and without viscous fluid dampers with regard to the parameters of the dampers  = 0.3 and  = 5000 kN⋅(s/m) 0.3 .A parametric study is performed to investigate the effects of variations in damping coefficient  and velocity exponent  of viscous fluid dampers on the buffeting responses of the bridge.To evaluate the effectiveness of passive viscous dampers, evaluation criteria  1 - 4 are selected, where  1 = RMS base shear of towers,  2 = RMS base moment of towers,  3 = RMS top displacement of towers, and  4 = RMS vertical displacement of bridge deck.To compare the response quantities of the bridge with and without passive viscous dampers, criteria reduction factor   is further defined as where   and  0 are th evaluation criteria of the bridge with passive viscous dampers and the uncontrolled bridge, respectively.
Figure 13 shows the plots of the criteria reduction factors  1 of the bridge with passive viscous dampers.It can be seen from Figures 13(a) and 13(c) that, for all the velocity exponents , the criteria reduction factors  1 for the RMS base shears of bridge towers number 1 and number 3 decrease with increasing damping coefficient  of the viscous fluid dampers.And for small velocity exponents such as  = 0.1 and 0.2, the criteria reduction factors  1 of the towers number 1 and number 3 decrease significantly, which indicates that the installation of viscous fluid dampers results in the significant increases in the base shears of bridge towers longitudinally unrestricted with the bridge deck.It can be seen from Figure 13(b) that, (1) for large velocity exponents such as  = 0.6∼1.0, the criteria reduction factor  1 increases with an increasing in damping coefficient ; (2) for medium velocity exponents such as  = 0.2∼0.5, the criteria reduction factor  1 increases with an increase in damping coefficient  up to a certain level, giving maximum values for the criteria reduction factor  1 ; (3) for small velocity exponents such as  = 0.1, the criteria reduction factor  1 decreases with an increase in damping coefficient .On the whole, the maximum criteria reduction factor  1 is approximately the same for different damping exponent .Moreover, the smaller the damping exponent , the smaller the damping coefficient  corresponding to the maximum criteria reduction factor  1 .The maximum criteria reduction factor  1 for the RMS base shear of tower number 2 is 57.01%, which indicates that the installation of viscous fluid dampers results in a reasonable reduction in base shears of bridge towers longitudinally restricted with the bridge deck.Figure 14 shows the plots of the criteria reduction factors  2 of the bridge with passive viscous dampers.It can be seen from Figure 14(a) that, for small velocity exponents such as  = 0.1 and 0.2, the criteria reduction factor  2 of the tower number 1 decreases significantly, which will result in the significant increase in the base moment of side tower number 1.However, for other velocity exponents the criteria reduction factor  2 of the tower number 1 increases slightly.
It can be seen from Figures 14(b) and 14(c) that the effects of variations in parameters of the viscous fluid dampers on the criteria reduction factors  2 of the towers number 2 and number 3 are similar to the criteria reduction factor  1 of tower number 2. The maximum criteria reduction factors  2 for the RMS base moment of towers number 2 and number 3 are 57.88% and 28.63%, respectively.
Figures 15 and 16 show the plots of the criteria reduction factors  3 and  4 of the bridge with passive viscous dampers.It can be observed from Figures 15 and 16 that the effects of variations in parameters of the viscous fluid dampers on the criteria reduction factors  3 and  4 are similar to the criteria reduction factor  1 of tower number 2. The maximum criteria reduction factors  3 for the RMS top displacement of towers number 1, number 2, and number 3 are 7.71%, 18.20%, and 23.90%, respectively.And the maximum criteria reduction factors  4 for the RMS vertical displacement in the middle of three girder spans between towers number 1 and number 2, number 2 and number 3, and number 3 and number 4 are 13.43%, 21.74%, and 24.13%, respectively.
Table 5 summarizes the extreme values and the average values for the criteria reduction factors  1 ∼  4 .The analytical investigations for the mitigation of longitudinal buffeting responses of the multispan cable-stayed bridge using viscous fluid dampers indicate the following: (1) The large base shears and moments of bridge towers longitudinally restricted with the bridge deck are significantly reduced.However, inappropriate parameters  and  of the viscous fluid dampers will result in significant increases in the base shears and moments of bridge towers longitudinally unrestricted with the bridge deck.(2) The RMS top displacements of central towers and vertical displacements of the bridge deck can be effectively reduced so that the large longitudinal buffeting responses of the bridge due to insufficient longitudinal stiffness of the central towers can be improved using viscous fluid dampers.

Multiobjective Optimization Problem.
The optimization design of vibration control with passive dampers is necessary for practical application [23].As shown in Figures 13∼16, individual optimal solutions for each criteria reduction factor are obviously different for the mitigation of longitudinal buffeting responses of the multispan cable-stayed bridge using viscous fluid dampers.For example, as shown in Figure 13 it is impossible to meet minimum base shears of bridge towers longitudinally restricted with the bridge deck such as tower number 2 and minimum base shears of bridge towers longitudinally unrestricted with the bridge deck such as towers number 1 and number 3. Thus, none of the optimal solutions allows simultaneous optimizations for the base shears of all of bridge towers simultaneously, which indicates that multiobjective optimization is necessary to obtain acceptable solutions.Multiobjective optimization provides a family of equally valid solutions (Pareto-optimal solutions) rather than a single optimal solution.As shown in Figure 17, a trade-off exists between the objectives such that an improvement in one objective cannot be achieved without detriment to another [16].The nondominating sort genetic algorithm II (NSGA-II), which is well-adopted in multiobjective optimization, is used in the present study.Deb et al. [24] have proposed an elitist, fast, and parameterfree multiobjective optimization genetic algorithm (NSGA-II).Specific details of the NSGA-II can be referenced in Deb et al. [24].

NSGA-II-Optimized Parameters for Viscous Fluid
Damper.In this section, the optimization design method based on the NSGA-II is proposed to determine the optimized parameters  and  of the viscous fluid dampers for the mitigation of longitudinal buffeting responses of the multispan cable-stayed bridge.As discussed in Section 5.2.1, the optimization design objective is to yield maximum reductions in the top displacements of central towers and base forces of bridge towers longitudinally restricted with the bridge deck, with hampering the significant gain achieved in the base forces of bridge towers longitudinally unrestricted with the bridge deck.Meanwhile, the favorable controlling effect on the vertical displacements of bridge deck is also worthy of attention.Therefore, a summary of the objective functions used in this study is given in Table 6.When the optimization is started, NSGA-II generates an initial population of chromosomes, each of which represents a set of Mathematical Problems in Engineering parameters  and  of the viscous fluid dampers.In this study, the initial population size is taken to contain 25 individuals.Next, NSGA-II creates a new generation of chromosomes based on multiobjective ranking using several simple GA operators such as selection, crossover, and mutation.As the number of generations increases, the control performance of the elite (i.e., nondominated) individuals is improved because of elitism [16].In this study, an upper limit on the number of generations is specified to be 50.The final number of nondominated individuals obtained after optimization is 19.
Pareto-optimal fronts between different objectives are shown in Figure 18.
As shown in Figures 18(a)∼18(d), the objective  1 is equivalent to the objectives  2 ,  3 ,  4 , and  5 , respectively.Thus, after multiobjective optimization the maximum reduction in the top displacement of central tower number 2 can also results in the maximum reductions in the vertical displacements of bridge deck, base forces of tower number 2, and base moment of tower number 1.Such optimization results are largely due to the coupling effects of vertical bending of the bridge deck and longitudinal bending of bridge tower.However, as shown in Figures 18(b), 18(d), and 18(e), the base shear of tower number 1 is in competition with the top displacement of tower number 2 and base forces of tower number 2. Here, one set effectively controls the base forces of bridge towers longitudinally restricted with the bridge deck and the top displacements of central towers but sacrifices control of base shears of bridge towers longitudinally unrestricted with the bridge deck.Thus, there is no individual that can appropriately reduce base shears of all of bridge towers.And the objectives  1 ,  2 ,  3 ,  4 , and  5 can be improved at the cost of the degraded objective  6 .In summary, for mitigation of longitudinal buffeting response of the multispan cable-stayed bridge, a design engineer can choose a set of proper parameters  and  of the viscous fluid dampers from Figure 18 that can satisfy the different performance requirements.If an engineer wants to minimize top displacements of central towers and base forces of bridge towers longitudinally restricted with the bridge deck regardless of base shear of tower number 1, the most appropriate parameters  and  of the viscous fluid dampers can be selected as the set A presented in Figure 18.If an engineer wants to control the base shear of tower number 1 with no increase or medium increase, the sets C or B presented in Figure 18 can be selected.Table 7 shows the control effects using three sets of A, B, and C.

Conclusions
The Jiaoshao Bridge in China is a six-tower cable-stayed bridge, which is the longest multispan cable-stayed bridge in the world.In this paper, the longitudinal buffeting responses of the bridge were obtained using a time-domain procedure  for analyzing buffeting responses implemented in ANSYS and the performance of viscous fluid dampers for the mitigation of longitudinal buffeting responses is further investigated.
From the analytical results of the present study the following conclusions are drawn: (1) Due to the insufficient longitudinal stiffness of the central towers, the longitudinal top displacements of central towers are significantly larger than those of side towers in the multispan cable-stayed bridge.The installations of viscous fluid dampers are beneficial for the reductions in the longitudinal top displacements of central towers.Due to the coupling modes of vertical bending of the bridge deck and longitudinal bending of bridge tower, the vertical displacements of the bridge deck can also be effectively reduced using viscous fluid dampers.
(2) In the multispan cable-stayed bridge, part of bridge towers longitudinally restricted with the bridge deck results in the concentration effect of longitudinal base shears and moments.The base shears and moments of bridge towers longitudinally restricted with the bridge deck are significantly larger than those of bridge towers longitudinally unrestricted with the bridge deck.The large base shears and moments of bridge towers longitudinally restricted with the bridge deck are significantly reduced using viscous fluid dampers.However, inappropriate parameters of the viscous fluid dampers will result in significant increases in the base shears and moments of bridge towers longitudinally unrestricted with the bridge deck.

Mathematical Problems in Engineering
(3) The optimization design objective is to yield maximum reductions in the top displacements of central towers and base forces of bridge towers longitudinally restricted with the bridge deck, with hampering the significant gain achieved in the base forces of bridge towers longitudinally unrestricted with the bridge deck.To this end, a multiobjective optimization design method that uses a nondominating sort genetic algorithm II (NSGA-II) is used to optimize parameters of the viscous fluid dampers.Optimization results reveal that a design engineer can choose a set of proper parameters of the viscous fluid dampers from Pareto-optimal fronts that can satisfy the different performance requirements.
Constraints between the girders and side towers number 2 and number 5 Auxiliary pier and end pier in each side span

Figure 4 :
Figure 4: Detailed finite element modeling for the constraints between the girders and towers (piers).

Figure 5 :
Figure 5: Mode shapes of the bridge.(a) 1st symmetric vertical bending of bridge deck + symmetric longitudinal bending of bridge tower, (b) 1st antisymmetric vertical bending of bridge deck + antisymmetric longitudinal bending of bridge tower, (c) 1st symmetric lateral bending of bridge tower, (d) 1st antisymmetric lateral bending of bridge tower, (e) 1st symmetric lateral bending of bridge deck + symmetric lateral bending of bridge tower, and (f) 1st symmetric torsion of bridge deck.

Responses of the Bridge 4 . 2 . 1 .
Buffeting Analysis Method.The simulated fluctuating wind speed time-histories are used as inputs to obtain the buffeting responses of the Jiashao Bridge in the time domain.

Figure 9 :Figure 10 :
Figure 9: RMS longitudinal buffeting responses of bridge towers.(a) RMS base shear, (b) RMS base moment, and (c) RMS top displacements of bridge towers.

Figure 11 :
Figure 11: Installation of viscous fluid dampers in the bridge.

Figure 12 :
Figure 12: Comparisons of time-history of vertical buffeting displacement with and without dampers.

Figure 13 :
Figure 13: Criteria reduction factor  1 for the RMS base shears of towers.(a)  1 of tower number 1, (b)  1 of tower number 2, and (c)  1 of tower number 3.

Figure 14 :
Figure 14: Criteria reduction factor  2 for the RMS base moments of towers.(a)  2 of tower number 1, (b)  2 of tower number 2, and (c)  2 of tower number 3.

Figure 15 :
Figure 15: Criteria reduction factor  3 for the RMS top displacements of towers.(a)  3 of tower number 1, (b)  3 of tower number 2, and (c)  3 of tower number 3.

Figure 16 :
Figure 16: Criteria reduction factor  4 for the RMS vertical displacements of bridge deck.(a)  4 in the middle of girder span between towers number 1 and number 2, (b)  4 in the middle of girder span between towers number 2 and number 3, and (c)  4 in the middle of girder span between towers number 3 and number 4.

Table 1 :
Vibration modes of finite element model of Jiashao Bridge.

Table 2 :
Simplified simulation of the wind field.

Table 4 :
RMS vertical displacements of the bridge deck (unit: m).

Table 5 :
Statistical results for the criteria reduction factors  1 ∼ 4 .

Table 6 :
Six objective functions. 4 for the RMS vertical displacement in the middle of girder span between towers number 2 and number 3  3  2 for the RMS base moment of tower number 2  4  2 for the RMS base moment of tower number 1  5  1 for the RMS base shear of tower number 2  6  1 for the RMS base shear of tower number 1