Abstract

This paper presents an investigation of the wind-induced buffeting responses of the Jiashao Bridge, the longest multispan cable-stayed bridge in the world. A three-dimensional finite element model for the Jiashao Bridge is established using the commercial software package ANSYS and a 3D fluctuating wind field is simulated for both bridge deck and towers. A time-domain procedure for analyzing buffeting responses of the bridge is implemented in ANSYS with the aeroelastic effect included. The characteristics of buffeting responses of the six-tower cable-stayed bridge are studied in some detail, focusing on the effects including the difference in the longitudinal stiffness between the side towers and central towers, partially longitudinal constraints between the bridge deck and part of bridge towers, self-excited aerodynamic forces, and the rigid hinge installed in the middle of the bridge deck. The analytical results can provide valuable references for wind-resistant design of multispan cable-stayed bridges in the future.

1. Introduction

For long-span cable-stayed bridges, the multispan cable-stayed 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 [14]. 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 2680m. Compared with a conventional three-span cable-stayed bridge with two towers, large temperature-induced deformation in the long bridge deck is one of the major problems in the design of multispan cable-stayed bridges. The commonly used structural measures to reduce the temperature effects include two aspects: (i) 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 and (ii) rigid hinge to release the temperature-induced deformation in the bridge deck such as the Jiashao Bridge in China. The aforementioned structural measures can improve the static performance of the multispan cable-stayed bridge under temperature action. However, studies of the wind-resistant performance of the multispan cable-stayed bridge using the partially longitudinal constraints and rigid hinge are meager. In addition, it is also necessary to investigate some unique features of the aerodynamic behavior of the multispan cable-stayed bridge considering the insufficient longitudinal stiffness of the central towers and self-excited force effect.

To date, an amount of research work had been done on investigating the buffeting response analysis of long-span bridges. The buffeting of long-span bridges is a type of vibration motion induced by wind turbulence. The buffeting responses are related not only to the wind field property and structural dynamic properties but also to the geometrical shape of bridge sections and the interaction between wind and structural motions [5, 6]. With the increase of the bridge length, buffeting analysis for slender bridge structures under strong wind actions becomes increasingly important. Buffeting analysis can be carried out in either frequency domain or time domain. Many early investigations were focused on the frequency domain analysis methods, in which statistical quantities in terms of spectral density and variance of structural response are obtained [7]. Compared with the frequency domain analysis methods, time-domain assessment of the buffeting response is convenient, since structural and aerodynamic nonlinearities may be taken more easily into account in the time domain. Hence, investigations of time-domain buffeting analysis have been done by several researchers in recent years [5, 813]. These investigations conducted into the application of time-domain buffeting analysis of long-span bridges provide valuable references for the wind-induced buffeting responses of the multispan cable-stayed bridges in the present study.

This paper investigates the wind-induced buffeting responses of the Jiashao Bridge, the longest multispan cable-stayed bridge in the world. The specific objectives of this study are to (i) perform the time-domain buffeting responses of the multispan cable-stayed bridge implemented in ANSYS with the aeroelastic effect included; (ii) investigate the distribution characteristics of buffeting responses of the bridge deck and bridge towers; and (iii) investigate the effects including the difference in the longitudinal stiffness between the side towers and central towers, partially longitudinal constraints between the bridge deck and part of bridge towers, self-excited aerodynamic forces, and the rigid hinge installed in the middle of the bridge deck on the buffeting responses of the multispan cable-stayed bridge.

2. 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 2680m with the span arrangement of 70m + 200m + 5 × 428m + 200m + 70m, which is the longest multispan cable-stayed 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-induced longitudinal deformation in the bridge deck, two important structural measures are applied in the design of Jiashao Bridge as shown in Figure 1(b): (i) The “drawer-type” rigid hinge shown in Figure 2 is installed in the midspan of the bridge deck. The rigid hinge can automatically adapt to changes in longitudinal deformation due to ambient temperature variations and meet the requirements of loading capacity under normal traffic conditions. (ii) 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 have been 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.

3. 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 (ANSYS 2009). The finite element model involves 1402 nodes and 1872 elements, as shown in Figure 3. 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 crossbeams connecting two steel box girders are modeled as rigid links at 30m intervals as shown in Figure 4(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 [14]. Rigid links at 15m intervals are used to connect the cables in four planes to the girders as shown in Figure 4(b). Figure 4(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) is described as follows:(i)Modeling of constraints between the girders and towers: as shown in Figure 5(a), constraints are applied to restrict the deck from moving in the , , and directions at bridge towers number 2 and number 5. Thus, the constraint equations of DOFs for the nodes shown in Figure 5(a) can be defined as , , and . 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 and directions. Thus, the constraint equations of DOFs for the nodes can be defined as and .(ii)Modeling of constraints between the girders and piers in each side span: as shown in Figure 5(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 shown in Figure 5(b) can be defined as , , , , , and .

Considering that the structural dynamic properties including modal frequencies and mode shapes form the basis of buffeting response analysis of the bridge, it is especially desirable to understand the dynamic characteristics of the bridge. Hence, the modal analysis of the Jiashao Bridge is conducted with the developed finite element model. The static equilibrium state of the bridge, which is the initial configuration for modal analysis, is achieved by geometrically nonlinear analysis of the bridge under dead loadings. The LANCZOS eigenvalue solver is adopted for modal analysis. It should be noted that, before implementing main FE modal analysis, the inaccuracy in FE modeling due to model meshing error should be eliminated considering mesh dependency. The mesh size of the FE model is determined by the convergence test of the mesh size and the analysis results reveal that the modeling meshing error caused by the discretization of the steel box girder can be ignored when the steel box girder is meshed according to the locations of the node of stay cables. And the number of FE elements of bridge towers is also determined using the convergence test. Figure 6 shows the relative errors of model frequencies of the first symmetric bending of bridge tower in direction and directions with regard to various numbers of FE elements of bridge towers, where the model frequencies with the 10 elements are assumed to be reference values. It can be seen that the modeling meshing error is reduced significantly with the increase of the number of FE elements. When the number of FE elements of bridge towers is larger than 8, the relative error is smaller than 0.05%. Thus, 8 elements of bridge towers are used in the FE model.

Main vibrations modes of finite element model are listed in Table 1 and part of them are shown in Figure 7. In order to verify the representing capabilities of the bridge dynamic behavior using the developed finite element model, Table 1 further shows the comparison of the calculated frequencies of the finite element model and experimentally identified frequencies using the ambient vibration test. It is clearly shown that a good agreement between calculated and identified frequencies of the Jiashao Bridge is obtained. Thus, the developed finite element model will serve as a basis for nonlinear buffeting analysis of the bridge.

The modal analysis results reveal the following dynamic properties of the bridge: (i) The floating vibration mode in direction does not occur due to the partially rigid system with longitudinal constraints between the bridge deck and two bridge towers number 2 and number 5, respectively. (ii) The natural frequency of the first symmetric bending of bridge tower in direction is found to be 0.2894Hz which is the third global vibration mode. Hence, the bridge towers are flexible with small lateral bending stiffness. (iii) The first bending mode of the bridge deck in direction with a frequency of 0.7087Hz occurs at the 21st global mode, which is higher than the corresponding value 0.2274Hz of the first bending mode in direction. Thus, the in-plane stiffness of bridge deck is weaker than the out-plane stiffness. (iv) The first torsional mode of the bridge deck with a frequency of 1.1361Hz occurs at the 43rd global mode. Thus, the higher torsional stiffness of the bridge deck is favorable for wind-resistant performance.

4. Nonlinear Buffeting Response of the Bridge

4.1. Simulation of the Wind Field

To perform the nonlinear buffeting analysis of the bridge in the time domain, eight one-dimensional independent multivariable stochastic wind fields are simulated based on the structural configuration as listed in Table 2. One of them is along the bridge deck, one is across the bridge deck, and the other six are across the six bridge towers, respectively. The bridge deck has 95 simulation points, spaced at 30m from left to right. And each bridge tower is divided into equal intervals of 10m from top to bottom, which results in 16 simulation points.

The longitudinal and vertical nonstationary fluctuating wind speeds in the overall wind field are modeled as 191-variate evolutionary vector process, , and a 95-variate evolutionary vector process, , respectively. The autospectra adopted for longitudinal and vertical components are Kaimal spectrum and Panofsky spectrum, respectively:where and are spectral-density function of longitudinal velocity fluctuation and vertical velocity fluctuation, respectively; is fluctuation frequency; is the Monin coordinate, in which is mean velocity at height and ; and is the friction velocity.

The following empirical exponential coherence function is adopted for the velocity cross-spectrum:where and are coordinates in the and directions, respectively, and 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.

The longitudinal and vertical fluctuating wind-velocity samples at 95 points on the bridge deck and 96 points on the bridge towers are simulated by using a modified Deodatis method, with a cut-off frequency of , number of fast-Fourier transform , and time interval . Table 3 summarizes the bridge’s input condition caused by wind when the wind-induced buffeting is analyzed. Figures 8 and 9 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 150m. It should be noted that wind velocities at different points on the bridge seldom reach the maximum value simultaneously. Thus, cross-correlation functions between two different points are important to describe the characteristics of spatial wind field. The correlation value decreases as the distance between the two points increases. As can be seen from the figure, the simulated autospectra and correlation functions are in well agreement with the targets.

4.2. Wind-Induced Buffeting Response 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. The buffeting analysis is carried out in the software package ANSYS with the following steps [12, 15]: A 3D finite element model of the bridge is established using the software package ANSYS. In order to account for self-excited aerodynamic forces, element aeroelastic stiffness and aeroelastic damping are both incorporated in buffeting analysis through the user-defined Matrix27 element in ANSYS [5, 10, 12]. The steady aerodynamic forces and buffeting forces are calculated by the measured aerodynamic coefficients and the Devaport formula, respectively [16]. 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 [5, 12]. It should be noted that the effect of the aerodynamic admittance is disregarded because it was not available from wind-tunnel tests. And the aerodynamic force coefficients and flutter derivatives are measured through the wind-tunnel tests.

4.2.2. Buffeting Responses of the Bridge Deck

The buffeting responses of the Jiashao Bridge were calculated at the design wind velocity of 49.72m/s at the bridge deck for the return period of 100 years. Figure 10 shows the time histories of vertical, horizontal, and torsional buffeting displacements of the 33rd simulation point at the bridge deck, which is located in the middle of the girder span between bridge towers number 2 and number 3. The root-mean-square (RMS) buffeting displacements of all simulation points at the bridge deck are further calculated as shown in Figure 11. It can be seen from the figure that the vertical and torsional buffeting displacements are approximately symmetrically distributed along the bridge deck. However, there is no symmetry in the distribution of the horizontal buffeting displacements. Moreover, the horizontal buffeting displacements at the girder midspan are significantly larger than those of other girder spans due to the installation of the rigid hinge in the middle of the midspan of the bridge deck, which will be discussed in the subsequent Section 5.2.

In order to investigate the contribution of the responses at different frequencies to the overall buffeting response, the power spectra of vertical, horizontal, and torsional buffeting displacements are calculated and Figure 12 shows the displacement spectra at the girder midspan. The following can be seen:(i)In the spectrum of vertical buffeting displacement, the first and second peaks correspond to the first and second symmetric bending modes of the bridge deck in direction, respectively. Thus, the vertical buffeting displacement at the girder midspan mainly depends on the symmetric bending mode of the bridge deck in direction.(ii)In the spectrum of horizontal buffeting displacement, the first and second peaks correspond to the 2nd symmetric bending of bridge tower in direction and the 1st symmetric bending of bridge deck in direction, respectively. Thus, the horizontal buffeting displacement at the girder midspan mainly depends on the coupling effect of the symmetric bending modes of bridge tower and bridge deck in direction.(iii)In the spectrum of torsional buffeting displacement, the first peak corresponds to the couple mode of the antisymmetric bending of bridge deck in direction and antisymmetric bending of bridge tower in direction, which indicates that the torsional buffeting displacement at the girder midspan mainly depends on the coupling effect of antisymmetric bending mode of bridge tower and bridge deck in direction.

4.2.3. Buffeting Responses of the Bridge Towers

The RMS buffeting responses of all simulation points at each bridge tower are calculated in the longitudinal and transverse directions, respectively. Figures 13(a) and 13(b) show the longitudinal RMS buffeting displacements and bending moments along the bridge towers, respectively. It can be seen that the longitudinal buffeting responses of bridge towers are approximately symmetrically distributed. And the longitudinal buffeting responses of the side towers (towers number 1 and number 6) are significantly smaller than those of other towers. This is because side towers are connected through outermost stay cables to the anchor piers, which can provide effective support to the side towers [2]. However, 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. In addition, the base bending moments of the towers number 2 and number 5 are significantly larger than those of towers number 3 and number 4 due to the partially longitudinal constraints at bridge towers number 2 and number 5. Thus, part of bridge towers longitudinally restricted with the bridge deck results in the concentration effect of longitudinal base moments. Figure 13(c) further shows the power spectrum of longitudinal buffeting displacement on the top of bridge tower number 3. It can be seen that the first, second, and third peaks correspond to the 1st, 9th, and 10th modes, which are all coupling modes of bending of the bridge deck in direction and bending of bridge tower in direction.

Figures 14(a) and 14(b) show the transverse RMS buffeting displacements and bending moments along the bridge towers, respectively. There is no symmetry in the distribution of the transverse buffeting responses of bridge towers. The transverse buffeting responses of towers number 2 and number 4 are maximum and the transverse buffeting response of tower number 6 is minimum. And there is no concentration effect of transverse base moments due to partially longitudinal constraints. Therefore, the effect of partially longitudinal constraints should be taken into account especially for analyzing the buffeting base moments of the bridge towers longitudinally restricted with the bridge deck. Figure 14(c) further shows the power spectrum of transverse buffeting displacement on the top of bridge tower number 3. It can be seen that the first peak corresponds to the 5th mode, namely, the 2nd symmetric bending of bridge tower in direction.

5. Discussions

5.1. Influence of Self-Excited Forces on the Buffeting Response

The self-excited aerodynamic forces are induced by interaction between bridge and wind motions. In this section, the influence of self-excited forces on the buffeting responses of Jiashao Bridge is discussed. Figure 15 shows the RMS of displacement responses at the bridge deck with and without considering self-excited forces. It can be shown in Figure 15(a) that the vertical buffeting displacements obtained with considering self-excited forces are significantly smaller than those obtained without considering self-excited forces. The RMS vertical displacement at the girder midspan obtained with considering self-excited forces is 2.57 times as that obtained without considering self-excited forces. However, Figures 15(b) and 15(c) show that the horizontal and torsional buffeting displacements obtained with considering self-excited forces are a little smaller than those obtained without considering self-excited forces.

Generally, the self-excited forces acting on the bridge towers can be disregarded. However, it is necessary to investigate the effect of self-excited force acting on the bridge deck on the buffeting responses of bridge towers. Figure 16 shows the RMS of displacements and bending moments of bridge tower number 3 with and without considering self-excited forces. It can be seen that the self-excited forces acting on the bridge deck have little effect on the transverse buffeting displacements and bending moments of the bridge tower. However, the longitudinal buffeting displacements and bending moments of the bridge tower are significantly reduced due to the self-excited force effect. The reason is that the mitigating effect on the vertical buffeting displacements of the bridge deck can also reduce the longitudinal buffeting responses of the bridge towers due to the dynamic coupling between the bending of bridge towers in direction and bending of the bridge deck in direction. Therefore, the self-excited force effect should be taken into account especially for analyzing the vertical buffeting responses of the bridge deck and longitudinal buffeting responses of the bridge towers.

5.2. Influence of Rigid Hinge on the Buffeting Response

The rigid hinge installed in the middle of the bridge deck is very important to release the temperature-induced deformation in the bridge deck. In this section, the effect of rigid hinge on the buffeting responses of the bridge is investigated. The contrast analytical model without rigid hinge is developed and the modal analysis result of the contrast model is shown in Table 4. It can be seen that the modal frequency of the first symmetric bending mode of bridge deck in direction of the contrast model without rigid hinge is 0.8610Hz, while the modal frequency of the original model with rigid hinge is 0.7087Hz. The relative variation is about 21.50%. And the rigid hinge has a limited effect on the other global vibration modes.

The buffeting response analysis of the contrast model without rigid hinge is carried out. The comparisons between the buffeting responses obtained with and without rigid hinge indicate that the installation of rigid hinge has little effect on the vertical and torsional buffeting responses of the bridge deck and buffeting responses of bridge towers. However, the installation of rigid hinge has a significant effect on the horizontal buffeting displacements of the bridge deck as shown in Figure 17(a). It can be seen that the horizontal buffeting displacement in the girder midspan obtained with rigid hinge is approximately 2 times as that obtained without the rigid hinge. Figure 17(b) shows the spectra of horizontal buffeting displacements obtained with and without rigid hinge. It can be shown that the frequency of the second peak corresponding to the 1st symmetric bending of bridge deck in direction increases from 0.7087Hz with rigid hinge to 0.8610Hz without rigid hinge. Moreover, the spectrum amplitude of the second peak obtained without rigid hinge is smaller than that obtained with rigid hinge. Therefore, the effect of the installation of the rigid hinge in the middle of the bridge deck should be taken into account for analyzing the horizontal buffeting responses of the bridge deck.

6. Conclusions

The Jiashao 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 nonlinear buffeting responses of the bridge were obtained using a time-domain procedure for analyzing buffeting responses implemented in ANSYS, focusing on the characteristics of buffeting responses of the bridge. From the analytical results of the present study the following conclusions are drawn:(1)The vertical and torsional buffeting displacements along the bridge deck are approximately symmetrically distributed. However, there is no symmetry in the distribution of the horizontal buffeting displacements. The power spectrum analysis results indicate that, for the girder midspan, the vertical buffeting displacement is dominated by the symmetric bending mode of the bridge deck in direction; the horizontal buffeting displacement is dominated by the coupling effect of the symmetric bending modes of bridge tower and bridge deck in direction; the torsional buffeting displacement is dominated by the coupling effect of antisymmetric bending mode of bridge tower and bridge deck in direction.(2)The longitudinal buffeting responses of bridge towers are approximately symmetrically distributed. And there is no symmetry in the distribution of the transverse buffeting responses of bridge towers. The power spectrum analysis results indicate that the longitudinal buffeting displacement is dominated by the coupling effect of bending of the bridge deck in direction and bending of bridge tower in direction; the transverse buffeting displacement is dominated by the bending of bridge tower in direction.(3)Due to the small longitudinal stiffness of the central towers, the longitudinal buffeting responses of the central towers are significantly larger than those of side towers in the multispan cable-stayed bridge. And part of bridge towers longitudinally restricted with the bridge deck results in the concentration effect of longitudinal base moments. The base moments of bridge towers longitudinally restricted with the bridge deck are significantly larger than those of bridge towers longitudinally unrestricted with the bridge deck.(4)The self-excited force effect has a significant mitigating effect on reducing the vertical buffeting displacements of the bridge deck and longitudinal buffeting responses of the bridge towers due to the dynamic coupling between the bending of bridge towers in direction and bending of the bridge deck in direction. And the self-excited force effect has a little effect on the horizontal and torsional buffeting displacements of the bridge deck and transverse buffeting responses of the bridge towers.(5)The installation of rigid hinge in the middle of the bridge deck has little effect on the vertical and torsional buffeting responses of the bridge deck and buffeting responses of bridge towers. However, the horizontal buffeting displacement in the girder midspan obtained with rigid hinge is approximately 2 times as that obtained without the rigid hinge. The reason is that the installation of rigid hinge reasonably reduces the modal frequency of the first symmetric bending mode of bridge deck in direction and increases the corresponding spectrum amplitude, which deserves special attention in practice.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors gratefully acknowledge the National Basic Research Program of China (973 Program) (no. 2015CB060000), the National Science and Technology Support Program of China (no. 2014BAG07B01), the National Natural Science Foundation of China (no. 51438002 and no. 51278106), the Program of “Six Major Talent Summit” Foundation (no. 1105000268), and Key Program of Ministry of Transport (no. 2013319223120).