A Quick Assessment and Optimization Method for a Flutter Aerodynamic Measure of a Typical Flat Box Girder

Flutter is one of the most serious wind-induced vibration phenomena for long-span bridges andmay cause the collapse of a bridge (e.g., the Old Tacoma Bridge, 1940). )e selection and optimization of flutter aerodynamic measures are difficult in wind tunnel tests. It usually takes a long time and consumes more experimental materials. )is paper presents a quick assessment and design optimization method for the flutter stability of a typical flat box girder of the long-span bridges. Numerical analysis could provide a reference for wind tunnel tests and improve the efficiency of the test process. Based on the modal energy exchange in the flutter microvibration process, the global energy input and local energy input are analyzed to investigate the vibration suppression mechanism of a flat steel box girder with an upper central stabilizer. Based on the comparison between the experimental and numerical data, a quick assessment method for the optimization work is proposed. It is practical to predict the effects of flutter suppression measures by numerical analysis. )us, a wind tunnel test procedure for flutter aerodynamic measures is proposed which could save time and experimental materials.


Introduction
As the spans of bridges have become larger, the bridge deck section of long-span bridges suffers aeroelastic phenomena of flutter and buffeting, as well as VIV or galloping, and the aerodynamic stability requirements of long-span bridges have become greater [1]. Flutter is one of the most serious wind-induced vibration phenomena for bridges, which causes the entire structure to fail. Researchers have carried out many studies to investigate the mechanism of flutter such as identification of flutter derivatives and the flutter aerodynamic measures [2,3].
Based on the coupling effect of aerodynamic forces, Yang et al. [4][5][6] established the relationship between torsional, vertical, and lateral vibration parameters and flutter derivatives (e.g., using the two-dimensional three-degree-of-freedom coupled flutter analysis method) and studied the flutter mechanism of two typical bridge sections: the closed box girder and two-isolated-girder sections. e flutter of these two sections is caused by negative aerodynamic damping, and the negative aerodynamic damping of the streamlined section is the result of the coupling effect of torsion and vertical motion. With the rapid development of CFD (computational fluid dynamics) methods [7], numerical simulation technology has been applied in the field of bridge wind engineering [8][9][10]. Vairo [11] used the RNG turbulence model to simulate the self-excited force acting on the typical main beam section. e identified flutter derivatives are basically consistent with the wind tunnel test results. Shirai and Ueda [12] used a nonlinear turbulence model k-ε to identify the flutter derivatives of the slotted section and the section with the central stabilizer, and the flutter mechanism was explained by its calculated flow field morphology and unsteady wind pressure distribution. Sarwar et al. [13] used the LES turbulence model to investigate the aerodynamic characteristics of a box girder bridge section. e flow around the streamline box girder section was analyzed, and the effect of section details on the aerodynamic characteristics of the bridge section was evaluated. Sun et al. [14] explored the turbulence model for the numerical calculation of fluid-solid coupling and calculated 18 flutter derivatives of the rectangular section with aspect ratios of B/D � 4 using the standard k-ω model; in this way, the influence of inflow turbulence intensity on the flutter derivatives was studied. Bai et al. [15] proposed an improved CFD method based on block-iterative coupling, performing 2D and 3D CFD simulations on two generic bridge deck sections to determine their aerodynamic force coefficients and flutter derivatives.Šarkić et al. [16] used the SST k-ω turbulence model to identify the flutter derivatives of a streamlined section, and their results are in good agreement with the experimental data.
In order to investigate the mechanism of the occurrence of flutter and the mechanism of vibration suppression by aerodynamic measures, many scholars have focused on energies and pressures. Scanlan [17] proposed a multimode flutter analysis technique from the viewpoint of the energy of the system. Based on the motion law of the vortex in discrete vortex calculation, Larsen [18] proposed a simplified analysis model that describes the motion of the vortex of the bridge section during a torsional motion cycle and estimates the total work done via the aerodynamic force generated by the vortex on the bridge section; the model then estimates the critical wind speed of the flutter. Chen et al. [19] studied the mechanism of multimodal coupled flutter from the perspective of system energy. Liu et al. [20] combined structure and airflow as a system to establish a new calculation method for the flutter analysis of suspension bridges from the perspective of energy balance. Liu et al. [21] used CFD to simulate the flutter process of the box girder and gave a definition for the aerodynamic energy calculation formula of the airflow input to the vibration model and the aerodynamic energy ratio of the surface area of the model. e flutter stability criterion was given from the viewpoint of energy balance. Yang et al. [22] used a method of flutter energy analysis to investigate the components of energy input and energy dissipation, and their results showed that components associated with flutter derivatives A * 1 and A * 2 produce the maximum positive and negative energy inputs, respectively. e maximum energy dissipation has a strong correlation with mechanical damping force. Tang et al. [23] studied the effects of different horizontal and vertical aerodynamic countermeasures, and flutter performance was judged by calculating the input energy by aerodynamic forces.
Although scholars have performed many wind tunnel tests and numerical simulations in the study of flutter mechanisms and vibration suppression measures, the selection and optimization of flutter aerodynamic measures are difficult to obtain in wind tunnel tests. In this study, the wind tunnel test method is used to study the energy input of the vibration suppression mechanism of the upper central stabilizer. e characteristics of global energy input and local energy input were investigated, which together demonstrate the mechanism from a macroscopic and microcosmic perspective. en, the numerical simulation method is used to accomplish bending-torsion coupling flutter for the box girder, and a numerical analysis of the effect of the upper central stabilizer on flutter stability was conducted, so the numerical simulation method can be applied to the subject, and the obtained result is determined to be applicable.

Wind Tunnel Test Setup and Data
Processing Method  Figure 1 shows the schematic diagram of the models. Meanwhile, 102 pressure taps were arranged at the middle section of the models. e schematic diagram of the pressure measuring point of the model is shown in Figure 2. e skeleton of the model was welded by aluminum alloy tubes. e surface of the model was made from aerospace wood and electrical boards.
e edges of the model were equipped with end plates made of acrylic. A picture of the sectional model for the construction stage is shown in Figure 3.
For the simulation of aeroelastic parameters in the wind tunnel test, a first-order symmetrical vertical bending mode and a first-order symmetric torsional mode were simulated for the girder of a long-span suspension bridge. Table 1 shows the model's design parameters.

Wind Tunnel Test Procedure.
e wind test was carried out in the CA-01 wind tunnel at Chang'an University. e test section was 15.0 m long, and the cross section was 3.0 m wide and 2.5 m high; the adjustable wind speed ranged from 0.2 to 53 m/s, and the turbulence of the uniform field was less than 0.3%. e wind tunnel test included synchronous vibration measurement, force measurement, and pressure measurement. In the free vibration tests (2-dof only-vertical and torsional), two laser displacement meters were installed at each end of the model to measure the displacement. During aeroelastic vibration, the electronic pressure scanner valves made by PSI were used to obtain the pressure, and the sampling rate is 312.5 Hz, as shown in Figure 4. In the aerodynamic force tests, the strain-gauge balance was used when the sectional model was fixed. rough the vibration test, the effects of different upper central stabilizers on the bending-torsion coupling flutter were measured. rough the wind force test, the aerodynamic force coefficients of the girder were obtained. Moreover, the wind pressure distribution on the surface of the model was measured synchronously to study the influence of different aerodynamic measures on the surface pressure distribution.

Data Processing Method.
e mean aerodynamic force coefficients are nondimensional parameters describing the magnitude of the static wind loads acting on the structure. rough a force test of the sectional model, this paper measured the drag force coefficient C D , the lift force coefficient C L , and the torque coefficient C M of each static model with different aerodynamic measures, which are defined as follows: where F D , F L , and M are the drag, lift, and torque of the sectional girder model, respectively, ρ is the air density, U is the mean longitudinal wind speed, and B and L are the height and width of the sectional girder model, respectively. e wind pressures on the models are expressed in the form of a nondimensional pressure coefficient, defined as where C P is the wind pressure coefficient at tap i and time t on the outer surface of the model and P 0 is the static reference pressure.  4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

Relationship between Amplitude and Wind Speed.
Based on the relationship between vibration amplitude and wind speed, it can be judged whether flutter occurred, and the critical wind speed of the flutter can be determined.

Shock and Vibration
During the critical state of the flutter, the vertical displacement and torsional displacement of the vibration system, during a short-term vibration period, are regarded as the simple harmonic vibrations of different phases at the same frequency, which are expressed as follows: where h 0 and a 0 are the vertical bending amplitude and the torsional amplitude, ω is the vibration circular frequency when the flutter occurs, and φ 1 and φ 2 are the initial phases of the vertical displacement and the torsional displacement.
(2) Time History of Global Aerodynamic Force. In this study, the time history of global self-excited aerodynamic force is obtained by full-section pressure integration. Using MAT-LAB programming to perform the integration process for pressures at several moments in several vibration cycles, the aerodynamic time history F L (t) and M(t) can be obtained to calculate the aerodynamic work.
(3) Global Aerodynamic Force Doing Work. For the twodegree-of-freedom bending-torsion coupling vibration system, if the calculation time t 0 is selected as a reference point, when the section vibrates until time t 1 , the cumulative work of the self-excited aerodynamic force for the vertical and torsional degrees of freedom is where W H se (t 1 ) and W T se (t 1 ) are the work done by the selfexcited aerodynamic force on the vertical and torsional degrees of freedom, respectively, and F L (t) and M(t) are the lift force and torque time history, respectively. us, the total work done by the self-excited aerodynamic force in the vibration system is (4) Mechanical Energy Increment. In the fluid-solid coupling system, the vibration system is not conservative, and the mechanical energy of the system is not a constant. e work done by the aerodynamic force is partially dissipated by structural damping, and the other portion is converted into the mechanical energy of the system. If one vibration period is considered, the mechanical energy increment of the system can be expressed as where W se and W D are the work done by the self-excited aerodynamic force and structural damping, respectively.

(5) Criterion for the Effectiveness of Flutter Aerodynamic
Measures Based on the Energy Method. From the perspective of energy balance, while the vibration system is in a complete vibration cycle, the mechanical energy increment ΔE M is greater than 0. If the sum of the aerodynamic work and the structural damping dissipative work is greater than 0, the flutter has a tendency to be excited. When the mechanical energy increment ΔE M is less than 0, the sum of the aerodynamic work and the structural damping dissipative work is less than 0, and the flutter has a tendency to be suppressed. us, when the mechanical energy increment ΔE M is equal to 0, the flutter is at a critical state. However, under normal circumstances, due to the small damping of the bridge structure, the structural damping dissipation work is also small. us, this phenomenon can be largely ignored. e positive or negative values of the aerodynamic cumulative work in a vibration cycle will directly determine whether the aerodynamic measures are effective. Figure 11 shows a comparison of the global aerodynamic cumulative work of the section with an upper central stabilizer in a microvibration cycle, including a 0°wind attack angle and a wind speed of 14.1 m/s. e cumulative work is the sum of the work, including each time step. Since the figure shows the cumulative work, the positive or negative slopes of the curve in the figure represent the positive or negative work during this period, respectively. If the slope is positive, positive work is done during this time. e positive or negative of the accumulated work value represents the input energy or the dissipated energy. erefore, the cumulative total work at the last moment of the vibration period in this figure will be used as the criterion for the effectiveness of the aerodynamic measures based on the energy method. It can be seen from the figure that the cumulative total work of the original section at the last moment is the largest among these three girder sections. In comparison, the cumulative total work of the section with a 0.3 H upper central stabilizer is relatively small, and the cumulative total work of the section with a 0.4 H upper central stabilizer decreases until it is negative.

Analysis of Calculation Results.
us, it can be judged that the 0.40 H upper central stabilizer is effective. When the cumulative total work of the section with a 0.30 H central stabilizer shows a small positive value at the last moment and no flutter occurs during the test, it is recommended to perform a plurality of cycles to calculate the work and determine the effectiveness of the aerodynamic measures. (2) e Local Aerodynamic Force Doing Work. For the twodegree-of-freedom bending-torsion coupling vibration system, if the calculation time t 0 is selected as the reference point, when the section vibrates until time t 1 , the cumulative work of the wind pressure for the local area x of the section surface is

Calculation Method
where P(x, t) is the static pressure at the section's surface, cos P → , v → h is the cosine of the angle between the directions that are perpendicular to the pressure tap and the direction of the vertical velocity, and cos P → , v → a is the cosine of the angle between the directions that are perpendicular to the pressure tap and the direction of the torsional linear velocity.

Analysis of the Calculation Results.
In order to compare the effects of different upper central stabilizers on the local energy input characteristics of the original section, the energy inputs acting on the upper and lower surface areas for one vibration period for different girder sections are calculated and are shown in Figures 13 and 14, respectively. It can be found that, for the upper surface of the model, the cumulative values of the original section are mostly negative, except in the downstream fairing region. e aerodynamic force dissipates energy during vibrations. In the upstream region, the cumulative work is all positive, and the self-excited aerodynamic force in this region continuously inputs energy into the vibration, which causes the original section to flutter. Once the 0.30 H upper central stabilizer is arranged on the girder section, the cumulative work in the upstream region reduces to nearly 0, and only the downstream region shows positive work. For the girder section with a 0.40 H central stabilizer, the cumulative work in the upstream region reduces to a negative value. e aerodynamic force in this region continuously dissipates energy during the vibration. At the same time, compared with the 0.3 H central stabilizer, the positive work occurring in the downstream region reduced as well, which indicates that the effect of the 0.4 H central stabilizer is better than that of the 0.3 H central stabilizer. e aerodynamic force in the upstream region of the upper surface of the model played a role in dissipating energy during vibration, making it the "critical region" for vibrations.

CFD Modeling.
In order to improve the efficiency of the wind tunnel test, the effectiveness of the upper central stabilizer is studied through numerical simulation software CFD. By establishing a CFD numerical model, simulating the bending-torsion coupling vibration, verifying the accuracy of the numerical model, and calculating the energy input, the effectiveness of the measure was quantitatively judged. e length of the full flow field was about 30 times the bridge width, the width of the full flow field was more than 20 times the beam height, and the centroid of the twodimensional bridge section was located at the flow field of about 1/3 the entrance. e number of mesh grids is 53,243. e SST k-ω turbulence model is adopted. In order to better simulate the flow field at the wall's surface and the wake region of the girder section, the grid was properly encrypted at the wall and wake. e time step is 0.0032, the spatial discretization format is second-order precision, and the SIMPLE algorithm is used. Grid meshing is shown in Figure 15. e bending-torsion coupling vibration was realized by the UDF program driving the moving grid to realize the forced coupling vibration.
To verify the accuracy of the CFD model, the comparison of mean aerodynamic force coefficients between the wind tunnel test and numerical simulation results is listed in Table 2. is comparison shows that the numerical simulation results are close to the wind tunnel test, which explains that the CFD model is reasonable.

Comparison of the Characteristics of Pressure Distribution.
e wind tunnel test and numerical simulation results of the pressure distribution on the upper and lower surfaces of different sections are compared in Figures 16 and 17, respectively. is comparison demonstrates that the numerical simulation results are consistent with the trend of the wind tunnel test results. e mean value of the pressure coefficient in the numerical simulation results is slightly higher than that of the wind tunnel test results.

A Comparison of the Local Energy Input in a Vibration
Period. e local energy input of the original section and the section with different upper central stabilizers at the upper and lower surface areas in a vibration cycle is given in Figures 19 and 20  Moreover, for the lower model surface of these three cases, the sum of all local work is also negative. erefore, these results are consistent with those in Section 5.3.1, and the central stabilizers with heights of 0.15 H, 0.30 H, and 0.40 H are effective.

Concluding Remarks
e selection and optimization of flutter aerodynamic measures are difficult to accomplish in wind tunnel tests; based on the angle of the energy exchanged in the microvibration process, both experimental and numerical methods are studied. e global and the local energy input characteristics from a macroscopic perspective were also defined. Furthermore, the mechanisms of the effectiveness of different upper central stabilizers for coupled bending and torsional flutter were studied. Moreover, according to the energy input calculation method proposed in this study, a numerical judgment of the effectiveness of the aerodynamic measures of a box girder section was discussed. e purpose of this study was to optimize aerodynamic measures through numerical analysis when the flutter of a bridge occurs in a wind tunnel test, which could provide a reference for wind tunnel tests and improve the efficiency of the test process. e main conclusions are as follows: (1) By calculating the global energy input characteristics, the cumulative total work of the model in a complete vibration cycle can be obtained. When the cumulative total work is negative, the global aerodynamic force dissipates energy in the vibration, and the aerodynamic measures can be judged to be effective. e local energy characteristics could explain the global energy input characteristics from a microscopic perspective, thereby providing a basis for optimizing aerodynamic measures. stabilizer may produce cumulative work in the upstream region of the upper surface of the girder. e aerodynamic force in this region plays a role in dissipating energy in the vibration system. us, this region is a "critical region" for flutter. (3) e numerical judgment method for the effectiveness of the upper central stabilizer has a certain reference value. When the upper central stabilizer is tested in a wind tunnel test, the recommended measures can first be obtained by numerical study, which may save time and facilitate the process.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.