Flutter and Buffeting Control of Long-span Suspension Bridge by Passive Flaps: Experiment and Numerical Simulation
Flutter and Buffeting Control of Long-span Suspension Bridge by Passive Flaps: Experiment and Numerical Simulation
International Journal of Aeronautical and Space Sciences. 2013. Mar, 14(1): 46-57
Copyright ©2013, The Korean Society for Aeronautical Space Science
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : January 08, 2013
  • Accepted : March 15, 2013
  • Published : March 30, 2013
Export by style
Cited by
About the Authors
Duc-Huynh, Phan
Department of Civil Engineering and Applied Mechanics, University of Education and Technology, Hochiminh city, Vietnam
Ngoc-Trung, Nguyen
Department of Mechanical & Biomedical Engineering, Kangwon National University, Chuncheon-si, Gangwon-do, South Korea

Flutter stability and buffeting response have been the topics of most concern in the design state of long-span suspension bridges. Among approaches towards the aerodynamic stability, the aerodynamic-based control method which uses control surfaces to generate forces counteracting the unstable excitations has shown to be promising. This study focused on the mechanically controlled system using flaps; two flaps were attached on both sides of a bridge deck and were driven by the motions of the bridge deck. When the flaps moved, the overall cross section of the bridge deck containing these flaps was continuously changing. As a consequence, the aerodynamic forces also changed. The efficiency of the control was studied through the numerical simulation and experimental investigations. The values of quasi-steady forces, together with the experimental aerodynamic force coefficients, were proposed in the simulation. The results showed that the passive flap control can, with appropriate motion of the flaps, solve the aerodynamic instability. The efficiency of the flap control on the full span of a simple suspension bridge was also carried out. The mode-by-mode technique was applied for the investigation. The results revealed that the efficiency of the flap control relates to the mode number, the installed location of the flap, and the flap length.
1. Introduction
Flutter has been realized to be one of the most dangerous aeroelastics of structures in civil engineering, especially after the disaster of the Tacoma Narrow Bridge. Flutter relates to both wind loads and structural response, from the perspective of control structure interactions. The flutter vibration mode is caused by self-excitation force from wind loads. When flutter occurs, the oscillatory motions regarding all degrees of freedom in the structure couple to create a single frequency called the flutter frequency . Also, the displacement is infinite by linear theory. Flutter could occur in both smooth and turbulent flows. The response of the bridge in turbulent flow in a wide range of wind speeds is called buffeting response . This response is normally used to determine the size of the structural members and requires careful evaluation at the design stage.
The question regarding the effective flutter and buffeting control is of growing interest as the span length of suspension bridges increases. This has resulted in a variety of techniques to reduce or suppress the aeroelastic response. The most common control methods can roughly be divided into the structure-based method and the aerodynamics-based method. The former approach exploits stiffened girders and/ or eccentric mass to control the flutter phenomenon passively whereas the latter approach focuses on generating stabilizing aerodynamic forces from the flow to solve the flutter problem.
A high torsion rigidity truss-stiffening girder was widely used in the design of a suspension bridge. For instance, the Akashi Kaikyo Bridge has a vertical stabilizer in the center span located along the centerline of the truss-type stiffening girder to improve aerodynamic stability. However, truss sections usually exhibit reasonably high wind forces (drag loading) which must be resisted by the bridge structure leading to the increase in the construction cost. The development of box-girders with a nearly streamlined section could also achieve the wind stability of a long span bridge. A twin-deck configuration with a center gap has been proven to be an effective means of improving the flutter characteristics (Fok et al., 2004a, 2004b; Kwok et al., 2012). However, the adequate stability cannot be assured beyond a certain limit in regard to length of spans (Brancaleoni, 1992; Ostenfeld and Larsen, 1992).
Research on aerodynamic control by using winglets and flaps (or being called control surfaces) was proposed and developed (Kobayashi et al., 1992, 1996, 1998, 2001 and 2005). The study was continued with experiments on passive flap control in turbulent and smooth flow (Phan and Kobayashi, 2011). However, the numerical simulation was not presented. Furthermore, the effects of flap control on the full span of the suspension bridge were not carried out. An extensive theoretical study on the active control of bridge flutter using a model similar to that proposed by Kobayashi was presented (Wilde et al., 1998; Preidikman and Mook, 1998; Nissen et al., 2004). The active model was then translated to a passive one (Wilde et al., 1999). Omenzetter et al. (2000, 2002) also proposed the passive flap control with springs and supplementary cables.
In order to study the effects of wind on structures, the time domain analysis may be the most effective approach. This approach can be done if forces acting on structures are expressed in the time domain. The famous self-excited forces suggested by Scanlan and Tomko (1971) contain frequencydependent coefficients. These forces are functions of flutter derivatives and reduced frequency. These flutter derivatives are found by experimental testing. To use these forces in the time domain, they have to be transformed into time dependency forces. This can be achieved by using indicial functions (Costa and Borri, 2006; Scanlan et al., 1974) or rational functions (Chen et al., 2000; Chen and Kareem, 2001; Cao and Sarkar, 2012). The simpler approach is to use quasi-steady theory. It is valid in both the time and frequency domain. These forces are the functions of the aerodynamic force coefficients: lift, drag, and moment coefficients.
In this study, the efficiency of mechanically driven flaps for controlling the aerodynamic stability of a bridge deck was investigated with both numerical simulation and experiment approaches whereby, a flap was mechanically controlled by the pitching motion of the deck and the aerodynamic damping or suppressing aerodynamic force was then produced. As shown in Figure 1 , the flaps were installed on both sides of the bridge deck and allowed the pitching motion. Providing the appropriate pitching motion of, the aerodynamic suppressing forces produced by them overcome the aerodynamic excitation, thus the flutter and buffeting were suppressed. Sectional model tests were also conducted to investigate the effect of the mechanically driving flap on the aerodynamic control. To explore the effects of flap control on both two-dimensional and three-dimensional cases, the quasi-steady theory has been selected with the related force coefficients obtained from the experiment in the wind tunnel. In addition, the numerical study on flutter and buffeting passive control for a super-long span suspension bridge (3,000 m central span) has been taken for the analysis. A three-dimensional model in conjunction with mode-by-mode analysis has been considered in an attempt to predict the flutter speed and buffeting response for full span with and without flaps. The method of analysis in essence is based on the quasisteady forces method as used in the two-dimensional analysis.
2. Two-dimensional analysis by quasi-steady forces
The bridge deck section in this paper was a shallow box girder and was aerodynamically treated as a flat plate. The controlling flaps at both edges of the bridge deck were also assumed to be flat plates. Thus, the numerical model of the bridge deck with the controlling flaps is described as shown in Figure 1 .
- 2.1 Control method
Figure 2 presents the bridge deck system having flaps which were mechanically driven by the rotational motion of
Lager Image
A bridge deck with flaps installed at leading and trailing edges
the deck. The flap was fixed with a hinge at the longitudinal edge of the bridge deck. Two hanger cables from a main cable were fixed to an anchor beam at its ends. A cross beam of the deck system was fixed to the anchor beam with a hinge. The gears Gs and Gf were fastened to the anchor beam and flap, respectively. Both gears were connected by a driving belt. When a bridge deck caused a torsion motion in its natural oscillation mode, hanger ropes nearly kept their vertical figure. The relative rotation angle was transmitted to the rotation of the flap through the driving belt with a given amplifying factor G = Rs/Rf where Gs and Gf are the radiuses of the gear Gs and Gf , respectively. Thus, the flap was driven just following the pitching motion of the bridge deck mechanically.
The same mechanism of the flap driving system was installed at the leeward side of the bridge deck. When the flaps were driven by the above mentioned manner, if we set the head-up motion of the bridge deck as α , the turning down of the leading edge flap was γ = and the turning up the trailing edge flap was β =- The same flap motions were seen when the wind direction was reversed.
- 2.2 Numerical analysis
The governing equations of motion with respect to the static equilibrium position of a two-dimensional bridge deck excited by aerodynamic forces are given by
Lager Image
In which m is the mass per unit length; I is the moment of inertial per unit length; ξh , ξα are the damping ratios of heaving motion h and pitching motion α , respectively; ωh , ωα are the natural circular frequencies of heaving and pitching motion, respectively; L is the total aerodynamic lift force (positive when upward); M is the total aerodynamic
Lager Image
Bridge deck with the mechanically driven flaps.
moment (positive when clockwise). L and M are defined in Eqs. (2).
By using the classical thin airfoil theory (Anderson 1984), the static aerodynamic lift L and the aerodynamic moment M are computed by:
Lager Image
In which ρ is the air density, B = 2 b is the bridge deck width (see Figure 1 ), U is the wind speed, W is the vertical wind, α is the angle of attack of the bridge deck, β is the angle of the trailing edge flap and γ is the angle of the leading edge flap.
The vertical wind is defined as follows:
Lager Image
Lager Image
is the mean vertical wind speed, which is assumed to be zero; w ( t ) is the vertical wind speed. The W was taken from a gust simulation of the wind tunnel. The coefficients C' are the slopes of the aerodynamic force coefficients. They are explained as follows:
Lager Image
where θ 0 =cos -1 c and c indicates the location of a hinge of the flap as shown in Figure 1 .
The flaps are passively driven by the pitching motion:
Lager Image
With the known values of all parameters, Eqs. (1) were solved by using the Runge-Kutta method to obtain the critical speed and the buffeting responses.
- 2.3 Numerical results based on the experimental data
The model dimensions and the dynamic properties of the investigated system are shown in Table 1 . In still air, the damping ratios and the natural frequencies were measured when the flaps were fixed.
The vertical wind speed W was actively simulated in the Eiffel type wind tunnel based on the Karman’s spectrum with the turbulence intensity of 5% (Kobayashi, et al. 1992, 1994). The results shown in Figure 3 were at the reduced wind speed U /( ωαb )=4.3, where ωα is the circular frequency of pitching motion (see Table 1 ).
The sensitivity numerical study in regard to the critical speed was performed with the combination of the coefficients X = G and Y =- G where G varied from -5.5 to 5.5. The instability phenomenon of the bridge deck, which could be flutter, divergence, or large buffeting response, was defined when the pitching response exceeded 0.17 radians in root mean square (RMS). The wind speed at this phenomenon is called the critical speed Ucr . The equations of motion in the turbulent flow were solved by the Runge
Model dimensions and dynamic properties of the investigated system
Lager Image
Model dimensions and dynamic properties of the investigated system
The slopes of lift and moment coefficients
Lager Image
The slopes of lift and moment coefficients
Lager Image
Time serials of the vertical wind velocity w at U/(ωαb)=4.3
Lager Image
Effect of G on the flutter control
Lager Image
Effect of G on buffeting control at U/(ωαb) = 3.08
Kutta method. The slopes of the aerodynamic force coefficients for the numerical simulation were obtained from experiments done in smooth flow. These values are listed in Table 2 .
Figure 4 , in which the vertical axis shows the critical speed Ucr compared with that of no control case U cr,0 , shows the variation of the critical wind speed due to the gain of the control surfaces, G . Increase in positive gain resulted in the improvement of the critical wind speed. The maximum ratio Ucr / U cr,0 was approximately 1.86 when G varied from -5.5 to 5.5. In the case of a further negative decrease of the gain G , the critical wind speed became smaller than that of the bridge deck without the control.
To check the effect of the amplifying factor G on buffeting, the percentage of reduction effect for the buffeting response ( PRE ) is introduced:
Lager Image
where RMScc and RMS c0 are the root mean squares of the root mean squares buffeting response for control cases and non-control cases, respectively; the subscript x indicates the heaving h or pitching α . Figure 5 shows the effect of the amplifying factor G on the buffeting at U /( ωαb )=3.08. When G increased, PREx increased. PREx reduced and received a negative value when G < 0.
3. Wind tunnel test
- 3.1 Wind tunnel model
The two-dimensional wind tunnel tests were conducted to investigate the effects of the mechanical flaps on the controlling of the wind-induced oscillations. The springmounting of the model and flaps are shown in Figure 6 . Details of the cross sectional shapes are also included.
Lager Image
The wind tunnel model
The model was fixed with the hinge at the point H α of the supporting frame. The frame was supported by the vertical spring kh which allowed only vertical motion of the frame including the model. The model was elastically supported by the pitching springs kα from the frame allowing only a pitching oscillation about the hinge H α . The sway motion and rotation of the frame was inhibited by the connecting rods AC and BD. Thus, the model has two degrees of freedom for oscillation, heaving and pitching.
The flaps were installed at both sides of the model. The leading-edge flap had a driving arm FH γ to its end of the rotation axis. During a heaving motion of the modelmounted system, this driving arm did not move. The vertical tie bar EF connected a given point of the driving arm and the frame. The pitching motion of the flap was given by the pitching motion of the bridge deck through the driving arm. The trailing-edge flap also had the same mechanism. Thus, the control manner described in the above section ( Figure 2 ) was realized in this model.
The magnitudes of the flaps depend on the location of the tie bars EF and JK as follows:
Lager Image
If l 1 = l 2 and G =( cb + l 1 )/ l 1 , the controlling angles are β =- , γ = .
- 3.2 Wind-induced responses
The static aerodynamic force coefficient by the flap was measured in a smooth wind when the angle of attack regarding the bridge deck was 0 degrees. The slopes of the coefficients (derivatives) of the lift and moment among these angles were estimated and compared with the theoretical values on the flat plate. The results are shown in Table 2 . A very small lift slope was obtained from the trailing edge flap. It is assumed to be caused by the following reasons. The trailing flap was in a region of the separated flow different from the theoretical flow. In addition, the small gap between the deck and the flap also reduced the aerodynamic force by the flap resulting in the differences between the theoretical and measured values as provided in Table 2 . The damping ratios and the frequencies of the model with the controlling flaps were measured in still air. The results are listed in Table 3 . These values changed when factor G changed.
The control effect of the flap control with several values of G was investigated in a turbulent flow. The results are shown in Figure 7 . In this figure, the numerical results are also included. The numerical simulation was made with parameters listed in Table 3 figure, ωα is the natural frequency of ωα pitching motion in the case of without control (see Table 1 ); b is the half-width of the bridge deck. In the experiment case, the flutter motion of the model without control took place at U /( ωαb )=3.69. This speed could be improved to Ucr /( ωαb )=6.46 by the flap control with G = 5.0. The flap motion by G = -1.4 was considered to give a reversed
Damping and frequency in smooth flow test
Lager Image
Damping and frequency in smooth flow test
Lager Image
Reduced wind speed at G on flutter control
Lager Image
Time trace, without control G = 0 and with control G = -1.4 at U/(ωαb)=2.77
effect for the aerodynamic response. The experimental result by G = -1.4 demonstrated that flutter appeared at a slightly lower speed compared with the case without control. The trends for the numerical and simulation results are similar.
An example for the time history of the response in regard to the model with and without flap control in a turbulent flow is shown in Figures 8 - 10 . The time experiment responses
Lager Image
Time trace, without control G = 0 and with control G = -1.4 at U/(ωαb)=3.38
Lager Image
Time trace, without control G = 0 and with control G = 5 at U/(ωαb) = 3.38
of heaving and pitching for the reverse effect are shown in Figure 8 . For the case without control, G = 0; and for the case with control G = -1.4. The reduced wind speed was Ucr / αb ) = 2.77. The amplitudes increased after controlling. The reversed effect happened. ωα was the value at G = -1.4. At the flutter speed for the case of G = -1.4, as shown in Figure 9 , the responses increased exponentially. If the value of G = 5,nse amplitudes decreased after controlling as shown in Figure 10 . The pitching motion completely reduced, but the heaving motion was not affected. If the model was controlled, the flaps were driven with the amplification factor of G = 5.0, the pitching response was suppressed, however almost no effect on the heaving response was observed.
4. The three-dimensional mode-by-mode analysis
- 4.1 The aeroelastic full model of suspension bridge
To apply the flap control to a full span bridge, the following problems were to be investigated. First, the relation between the install flap length and the critical wind speed was examined. Second, the influence of the vibration mode number on the control result was checked. Finally, the buffeting response of the suspension bridge was also tested.
In this numerical study, a simple span suspension bridge with span length l = 3000 m was employed, which is shown in Figure 11 . Details of the bridge model and the flaps for
Lager Image
The suspension bridge model
Structural properties of the full bridge model and the flaps
Lager Image
Structural properties of the full bridge model and the flaps
this simulation are listed in Table 4 . The damping ratios of bending and torsion modes were assumed to be 0.01 for the mode-by-mode study. This assumption was made for all modes.
The vibration mode shapes and the natural frequencies of bending and torsion motion were obtained by using the finite element method (Abdel-Ghaffar, 1976). The results are shown in Table 5 for the first twenty modes. The lateral motion was neglected in this study.
- 4.2 Simulation of turbulent wind field
The simulation of the vertical turbulent wind field acting on the bridge was carried out. The spectra representation method was used (Paola 1998, Deodatis 1996, Cao et al . 2000, Ding et al . 2006) with the Von-Karman’s spectrum formula. The Davenport’s coherence function (Davenport 1968) between point j and m was obtained by the following formula:
Lager Image
Lager Image
Sample of vertical turbulent wind at U = 10 m/s
In which ω is the frequency; λ is the dimensionless decay factor; Δ jm =Δ| j - m | is the horizontal distance from point j to m . △ is the distance interval between successive points. These parameters are listed in Table 6 . Figure 12 presents the first 1000 sec samples of the vertical turbulent wind field at point P 1 (left support), P 2, P50 (midpoint), P 51, P 100 and P 101 (right support). The mean wind speed in this case was 10 m/s.
- 4.3 Numerical analysis
The bridge deck motions were assumed by the combination of one bending mode and one torsion mode as follows:
Lager Image
where ϕb ( x ), ϕt ( x ) are the mode shapes of bending and torsion, respectively; gb ( t ), gt ( t ) are the generalized coordinates of bending and torsion, respectively; x is the span-wise coordinate. In this study, the first three modes of
The first twenty natural frequencies of the bridge
Lager Image
The first twenty natural frequencies of the bridge
Characteristic values of the wind field simulation
Lager Image
Characteristic values of the wind field simulation
bending and torsion were examined. The lateral mode was not considered. The response calculation was conducted based on the modal analysis with the second order equations of motion composed of one of the bending modes and one of the torsion modes.
Referring to Eqs. (1), the equations of motion in modal analysis become
Lager Image
where I 0b and I 0t are the bending and torsion generalized inertias, respectively
Lager Image
The buffeting forces and the self-excited forces are the sum of the forces on the bridge deck and the flaps. U ( x , t ) is assumed to be constant with time and space. The vertical turbulent wind w ( x , t ) has a span-wise correlation. The flap with length △ lF is installed at xj , j = 1, 2... N , where N is the divided number of the span.
By replacing α in Eqs. (5) by Eqs. (9), with X = G and Y = -G , the flap angles at position x and at time t are driven by the following equations:
Lager Image
In Eqs. (12), the coefficients G = 5 is used for the control case. The three methods of the flap arrangement as shown in Figure 13 are employed here. In method 1, the flap is installed partially at the center of the span. The installed portion corresponds to the hoop of the first natural vibration mode shape. The locations of the flap of method 2 and 3 correspond to the hoops of the natural vibration modes for the second and the third, respectively.
Bridge deck motions are assumed by the combination of one bending mode and one torsion mode, bending and torsion of bridge deck are defined by the Eqs. (9) and (10). If the flap angles are derived by Eqs. (12), the equations of motion for the full span of the bridge are:
Lager Image
Lager Image
Flap arrangement.
Critical speeds – without flap
Lager Image
Critical speeds – without flap
Lager Image
where ψF is a forcing function taking the value of 1 where flap is installed and is 0 otherwise.
- 4.4 Numerical results
At first, the critical speed of the bridge deck without controlling was checked. The result is shown in Table 7 . In which, modes( i-j ) means the combination of the i th bending mode and the jth torsion mode. The term ω α1 is the torsion natural frequency of the first mode. For the mode combinations (1-1), (2-2) and (3-3), the critical phenomenon was flutter. In addition, the divergence motion took place in the mode combinations (1-2), (2-1), (2-3) and (3-2). In the design of a bridge, only the lowest critical speed, in this case the mode combinations (1-1), is important. However, in this study, all the mode combinations shown in Table 7 were investigated.
The critical speeds were checked for the arrangement methods with a variation of the flap length ratio, and modes( i-j ). The flap length ratio FLR is defined by the ratio of the flap length on one side and the span length. The results are depicted in Figure 14 . In case of modes(1-1) without control, the flutter speed was 1.26. After controlling, the flaps installed for about 50% of the span length could increase the flutter speeds of the modes(1-1) up to 2.28, approximately 1.8 times compared to without control. From the results of this analysis and the plots in Figure 14 , it was determined that the flaps could improve the critical speed of the bridge. They controlled the flutter speed well.
The arrangement method 1 could increase the flutter speeds of the modes(2-2) up to 4.19, approximately 2.34 times compared to without control. It was also effective for divergence. Modes(1-2), (2-1), (2-3), (3-2) were the divergence cases. Modes(1-2) and (3-2) were of the same effect, because of the same pitching mode 2. In the real case, the arrangement method 2 is more realistic since the length of the hung cable was long enough to produce the relation motion between the flap and bridge, as discussed above. This method also increased the flutter speed and the divergence speed. The same effect could be seen in the arrangement method 3.
5. Conclusion and Discussion
A mechanically turning flap to control the aerodynamic stability was studied with the quasi-steady force approach. We have demonstrated that the quasi-steady force is easy to be applied in the time domain. In the two-dimensional study, the experimental force coefficients and vertical wind gust were used to both simulate and find the best control. The amplification factor G was investigated in the viewpoint of the large flutter speed and the small buffeting response. The following remarks were addressed:
  • - The quasi-steady force was reasonably accepted to
Lager Image
Flutter speed improvement by FLR with different modes
  • investigate the efficiency of the amplification factorG. The control manner for the leading-edge flap and the training-edge flap such asγ=Gαandβ=-Gαcan be easily set up in a real case. The positive values ofGcan be applied to suppress the buffeting response and improve the flutter speed or to improve the aerodynamic stability of the bridge deck. The negative values ofGhad the reverse effect. The numerical results using the experimental frequencies and damping ratios were compared to the experiment cases to check the efficiency of the flap control and to validate the quasi-steady force method.
  • - In a turbulent flow, control byG= 5.0 also improved the flutter speed approximately 2 times and suppressed the divergence phenomenon up to its wind speed. Buffeting in the pitching motion was effectively suppressed but the heaving motion was not.
The quasi-steady force method was also applied for full span bridge to check the efficiency of the flap control. Effects of the flap location, the flap length, and the mode shapes on the control were investigated. The mode-by-mode technique of coupled motions was considered. The flutter speed and the buffeting response for with and without the control were computed. The parameters used in the numerical analysis were the same with those in the two-dimensional cases. The following results were found by the mode-by-mode technique:
  • - The flaps installed for about 50% of the span length could increase the flutter speed of the modes(1-1) up to approximately 1.8 times.
  • - For divergence, the passive flap control was proven to be effective.
Anderson J.D 1984 Fundamentals of Aerodynamics McGraw-Hill New York
Brancaleoni F. 1992 “The construction phase and itsaerodynamic issues, Aerodynamics of Large Bridges (ed.Larsen)” Balkema 147 - 158
Brown W. C. 1996 “Development of the deck for the 3300mspan Messina” 15th IABSE Congr. Rep., IABSE Zurich 1019 - 1030
Brown W.C. 1999 “Long span bridge project - a personalview” Springer Proc. of the international seminar on Long-Span Bridges and Aerodynamics 3 - 19
Cao B. , Sarkar P. P. 2012 “Identification of RationalFunctions using two-degree-of-freedom model by forcedvibration method” Engineering. Structures 43 21 - 30    DOI : 10.1016/j.engstruct.2012.05.003
Chen X. , Matsumoto M. , Kareem A. 2000 “Aerodynamic coupling effects on flutter and buffeting ofbridges” J. of Eng. Mech. 126 (1) 7 - 26    DOI : 10.1061/(ASCE)0733-9399(2000)126:1(7)
Chen X. , Kareem A. 2001 “Nonlinear response analysisof long-span bridges under turbulent winds” J. Wind Eng. Ind. Aerodyn. 89 (14) 1335 - 1350    DOI : 10.1016/S0167-6105(01)00147-7
Costa C. , Borri C. 2006 “Application of indicialfunctions in bridge deck aeroelasticity” J. Wind Eng. Ind. Aerodyn. 94 (11) 859 - 881    DOI : 10.1016/j.jweia.2006.06.007
Dung N. , Miyata T. , Yamada H. 1996 “Application ofrobust control to the flutter of long – span bridges” J. Struct. Eng. 42A 847 - 853
Fok C.H. , Kwok K.C.S. , Qin X.R. , Hitchcock P. A. 2004 “Sectional pressure tests of a twin-deck bridge: part1: experimental techniques and effects of angle of windincidence” Proc. of 11th AWES Workshop Darwin, Australia
Fok C.H. , Kwok K.C.S. , Qin X.R. , Hitchcock P.A. 2004 “Sectional pressure tests of a twin-deck bridge: Part 2:effects of gap-width on a twin-deck configuration”. Proc. of 11th AWES, Workshop Darwin, Australia
Ge Y. , Zou X. , Yang Y. 2009 “Aerodynamic stabilizationof central stabilizers for box girder suspension bridges” Wind and Structures 12 (4) 285 - 295
Gua M. , Chang C.C. , Wua W. , Xiang H.F. 1998 “Increase of critical flutter wind speed of long-span bridgesusing tuned mass dampers” J. Wind Eng. Ind. Aerodyn. 73 111 - 123    DOI : 10.1016/S0167-6105(97)00282-1
Gua M. , Chena S.R. , Chang C.C. 2001 “Parametricstudy on multiple tuned mass dampers for buffeting controlof Yangpu Bridge” J. Wind Eng. Ind. Aerodyn. 89 987 - 1000    DOI : 10.1016/S0167-6105(01)00094-0
Gua M. , Chena S.R. , Chang C.C. 2002 “Control ofwind-induced vibrations of long-span bridges by semiactivelever-type TMD” J. Wind Eng. Ind. Aerodyn. 90 111 - 126    DOI : 10.1016/S0167-6105(01)00165-9
Kobayashi H. , Nagaoka H. 1992 “Active control of flutter of a suspension bridge” J. Wind Eng. Ind. Aerodyn. 41 (1 - 3) 143 - 151    DOI : 10.1016/0167-6105(92)90402-V
Kobayashi H. , Hatanaka A. 1992 “Active generationof wind gust in a two-dimensional wind tunnel” J. Wind Eng. Ind. Aerodyn. 42 (1-3) 959 - 970    DOI : 10.1016/0167-6105(92)90102-G
Kobayashi H. , Hatanaka A. , Ueda T. 1994 “Activesimulation of time histories of strong wind gust in a windtunnel” J. Wind Eng. Ind. Aerodyn. 53 315 - 330    DOI : 10.1016/0167-6105(94)90089-2
Kobayashi H. , Nitta Y. 1996 “Active flutter control ofsuspension bridge by control surfaces” Third International Conference on Motion and Vibration Control Chiba 1 - 6
Kobayashi H. , Ogawa R. , Taniguchi S. 1998 “Activeflutter control of a bridge deck by ailerons” Proc, 2nd World Conf. on Structural Control Kyoto
Kobayashi H. , Mitani K. , Ogawa 2001 “Activebuffeting control by flaps” The Fifth Asia-Pacific Conference on Wind Engineering
Kobayashi H. , Phan D.-H. 2005 “Bridge deck fluttercontrol by control surfaces” Proc, 6th Asia-Pacific Conf. on Wind Engineering Seoul, Korea
Körlin R. , Starossek U. 2007 “Wind tunnel test ofan active mass damper for bridge decks” J. Wind Eng. Ind. Aerodyn. 95 267 - 277    DOI : 10.1016/j.jweia.2006.06.015
Kwon S. D. , Jung S. , Chang S.P. 2000 “A new passiveaerodynamic control method for bridge flutter” J. Wind Eng. Ind. Aerodyn. 86 187 - 202    DOI : 10.1016/S0167-6105(00)00010-6
Kwon S. D. , Park K.-S. 2004 “Suppression ofbridge flutter using tuned mass dampers based on robustperformance design” J. Wind Eng. Ind. Aerodyn. 92 919 - 934    DOI : 10.1016/j.jweia.2004.05.006
Kwok K.C.S. , Qin X.R. , Fok C.H. , Hitchcock P.A. 2012 “Wind-induced pressures around a sectional twindeckbridge model: Effects of gap-width on the aerodynamicforces and vortex shedding mechanisms” J. Wind Eng. Ind. Aerodyn. 110 50 - 61    DOI : 10.1016/j.jweia.2012.07.010
Lin Y. Y. , Cheng C.-M. , Lee C.-H. 2000 “A tuned massdamper for suppressing the coupled flexural and torsionalbuffeting response of long-span bridges” Eng. Struct. 22 1195 - 1204    DOI : 10.1016/S0141-0296(99)00049-8
Miyata T. , Yamada H. , Dung N. , Kazama K. 1994 “Onactive control and structural response control of the coupledflutter problem for long span bridges” 1st World Conf. on Structural Control Los Angeles, California, USA
Nissen H. D. , Sørensen P. H. , Jannerup O. 2004 “Activeaerodynamic stabilisation of long suspension bridges” J. Wind Eng. Ind. Aerodyn. 92 829 - 847    DOI : 10.1016/j.jweia.2004.03.012
Okada T. , Honke K. , Sugii K. , Shimada S. , Kobayashi H. 1998 “Suppression of coupled flutter of a bridgedeck by tuned pendulum damper” Proc, 3rd World Conf. on Structural Control Kyoto
Omenzetter P. , Wilde K. , Fujino Y. 2000 “Suppressionof wind-induced instabilities of a long span bridge by apassive deck-flaps control system. Part I: Formulation” J. Wind Eng. Ind. Aerodyn. 87 (1) 61 - 79    DOI : 10.1016/S0167-6105(00)00016-7
Omenzetter P. , Wilde K. , Fujino Y. 1 “Suppressionof wind-induced instabilities of a long span bridge bya passive deck-flaps control system. Part II: Numericalsimulations” J. Wind Eng. Ind. Aerodyn. 87 (1) 81 - 91    DOI : 10.1016/S0167-6105(00)00017-9
Omenzetter P. , Wilde K. , Fujino Y. 2002 “Studyof passive deck-flaps flutter control system on full bridgemodel. I: theory” J. Engrg. Mech. 128 (3) 264 - 279    DOI : 10.1061/(ASCE)0733-9399(2002)128:3(264)
Omenzetter P. , Wilde K. , Fujino Y. 2002 “Studyof passive deck-flaps flutter control system on full bridgemodel. II: results” J. Engrg. Mech. 128 (3) 280 - 286    DOI : 10.1061/(ASCE)0733-9399(2002)128:3(280)
Ostenfeld K. H. , Larsen A. 1992 “Bridge engineeringand aerodynamics, Aerodynamics of Large Bridges, LarsenA. (ed.)” Balkema Rotterdam
Peidikman S. , Mook D.T. 1998 “On the developmentof a passive-damping system for wind-excited oscillation oflong-span bridges” J. Wind Eng. Ind. Aerodyn 77 (78) 443 - 456    DOI : 10.1016/S0167-6105(98)00163-9
Phan D.-H. , Kobayashi H. 20011 “An experimentalstudy of flutter and buffeting control of suspension bridge bymechanically driven flaps” Wind and Structures 14 (2) 152 - 163
Sato H. , Kusuhara S. , Ogi K. , Matsufuji H. 2000 “Aerodynamic characteristics of super long-span bridgeswith slotted box girder” J. Wind Eng. Ind. Aerodyn 88 (2-3) 297 - 306    DOI : 10.1016/S0167-6105(00)00055-6
Sato H. , Hirahara N. , Fumoto K. , Hirano S. , Kusuhara S. 2002 “Full aeroelastic model test of a super longspanbridge with slotted box girder” J. Wind Eng. Ind. Aerodyn. 90 (12-15) 2023 - 2032    DOI : 10.1016/S0167-6105(02)00318-5
Scanlan R. H. , Béliveau J-G. , Budlong K. S. 1974 “Indicial aerodynamic functions for bridge decks”. J. The Eng.Mech. Division 100 (4) 657 - 672
Scanlan R. H. , Tomko J. J. 1971 “Indicial aerodynamicfunctions for bridge decks”. J. The Eng. Mech. Division 97 (6) 1717 - 1737
Shubov M. A. 2004 “Mathematical modeling and analysisof flutter in long-span suspension bridges and in bloodvessel walls” J. Aero. Eng 17 (2) 70 - 82
Songpol P. 1998 “Analytical study on flutter suppression byeccentric mass method on 3D full suspension bridge model” Proc, 3rd World Conf. on Structural Control Kyoto
Ubertini F. 2010 “Prevention of suspension bridge flutterusing multiple tuned mass dampers”. Wind and Structures 13 (13) 235 - 256
Wilde K. , Fujino Y. , Kawakami T. 1999 “Analyticaland experimental study on passive aerodynamic control offlutter of a bridge deck” J. Wind Eng. Ind. Aerodyn. 80 105 - 119    DOI : 10.1016/S0167-6105(98)00196-2
Wilde K. , Fujino Y. 1998 “Aerodynamic control ofbridge deck flutter by active surfaces” J. Eng. Mech. 124 (7) 718 - 727
Wilde K. , Fujino Y , Prabis V. 1996 “Effects of eccentricmass on flutter of long span bridge” Proc, 2nd Int. Workshop on Structural Control Hong Kong
Zhang C. W. , Li J. L. , Li H. , Ou J. P. 2011 “PreliminaryNumerical Study on TRID System for Flutter VibrationControl of Bridge Structure” Procedia Engineering 14 2796 - 2806    DOI : 10.1016/j.proeng.2011.07.352