Buffeting performance of long-span suspension bridge based on measured wind data in a mountainous region

Long-span suspension bridge increases rapidly in size as a result of bridge construction in a mountainous region, in addition, more and more long-span suspension bridges are in process of preparation. The bridge stiffness decreases with the increase of bridge span length, and hence the buffeting performance of bridge is sensitive to external factors. In this paper, the Cuntan Yangze Bridge located in a mountainous region is taken as the background to study the effect of different power spectrums on the buffeting performance. A three-dimensional finite element model is set up on the ANSYS platform. The fitted power spectrum of extreme strong wind is recorded and taken as the sample to analyze the buffeting performance. The results are compared with the specified power spectrum in the time and frequency domains. Different from existing studies, buffeting performances with the fitted power spectrum are larger than those with the specified power spectrum on the whole. Two kinds of power spectrum are coincidental in the overall tendency in the frequency domain and are distinct in the low frequency region. Structure performance of long-span suspension bridge in the mountainous region should be the subject of specially paid attention.


Introduction
Long-span bridges have been built throughout the whole country of China with the rapid development of the economy and technology.At the design stage of long-span bridges, suspension bridges are always taken as the priority selection from all the types of bridges because of a large span capacity.The bridge stiffness decreases with the increase of bridge span length, and hence the long-span bridge is sensitive to wind-induced vibration.Wind-induced vibration is a key of long-span bridge considered by designers and bridge safety regulators.To successfully manage the bridge, buffeting suppression is the objective of aerodynamic control [1].

Suspension bridge and wind characteristic in a mountainous region
With the socioeconomic development of the Western China taken up mainly by large mountainous areas, there is a great need for good transport.Therefore, the long-span suspension bridge increases rapidly in size as a result of the bridge construction in the mountainous region, in addition, more and more long-span suspension bridges are in process of preparation.Affected by the topography, wind characteristics in the mountainous region are different from the standard wind characteristics.The influence of the topography on the wind field and the associated effects on the displacement response are remarkable [2,3].In order to get the buffeting performance closer to the realistic environment, the wind speed and related wind characteristic parameters of the bridge are derived from the wind measure system at the bridge site.

Buffeting performance of bridge
Buffeting performance of bridge is a forced vibration caused by wind turbulence.Fatigue damage of bridge components and discomfort for vehicles and pedestrians can be interfered with the buffeting performance.
Buffeting control work of a long-span bridge has always been a hot discussion researched by many scholars.Buffeting-induced mechanisms are the foundation of buffeting control, such as structural parameters [3], hanger oscillation of suspension bridge [4], non-stationary characteristics induced by typhoon [5].The buffeting suppression method includes active and passive buffeting control [6][7][8].Optimal buffeting control of multiple tuned mass dampers [9][10][11][12] and mechanically driven flaps [13][14][15] has been a hot area of research for recent years.
As to long-span bridge engineering, researchers have been working to effectively improve the calculating precision on the basis of classical theories and meet the engineering requirements.The improvements may be concluded in three aspects, the bridge with more details considered, and the outer environment more agree to the real condition and more reasonable analysis theories.Influences of buffeting performance with more bridge details were given including mid-tower, continuous span, catwalk, slotted deck and central buckle [16][17][18][19].Outer environment conditions more correspond to the reality and considered such astopography, extreme value of typhoon, turbulent spatial correlation coefficient, skew wind, design and measured power spectrum, non-stationary and stochastic excitation [2,5,[18][19][20][21][22][23][24][25].Improved analysis theories were applied in the buffeting performance including linear and non-linear regression analysis, evolutionary power spectral density, varying frequency-increment sweeping method, three-dimensional simulation [19,21,24,[26][27][28][29][30].
In conclusion, the buffeting performance with consideration of wind characteristics can reflect the response of bridge closer to the real condition.However, most of previous researches were concentrated in sea-crossing and river-crossing bridges in design or specified power spectrum of turbulent wind.The buffeting performance of long-span bridges based on the measured wind data is still insufficient.The long-span suspension bridge is about to appear in large numbers in the mountainous region, and hence its buffeting performance in extreme strong wind can provide references for other bridges in a similar region.
In this paper, the long-span suspension bridge and wind measure system were introduced.The real-time wind data recorded from wind measure system was recorded and analyzed.A three-dimensional finite element model was set up on the ANSYS platform.The extreme strong wind record in a wind database was taken as a sample to analyze the buffeting performance of long-span suspension bridge, and the results were compared with the specified power spectrum in the time and frequency domains.

Project overview
Cuntan Yangtze Bridge including main section and approach bridge locates in Chongqing, China.It's a streamlined box-girder suspension bridge with twin towers.Span arrangement is 250.0 m+880.0m+250.0 m.The main girder of the Cuntan Yangtze Bridge is a streamlined closed, flat box girder.The total width of the girder is 42.0 m, and the height at the centerline is 3.5 m.Rise-to-span ratio in this bridge is 1/8.8, and the distance of two main cables is 39.2 m.The bridge elevation is shown in Fig. 1(a), and the design sketch after construction is shown in Fig. 1(b).The north tower is 194.5 m, and the south tower is 199.5 m in height.The main cable is made up of 127 high-strength galvanized steel wires.

Wind measure system
The topography influence on the wind field and the associated effects on the wind-induced vibration are outstanding [31,32].The Cuntan Yangtze Bridge has two main towers, and the wind measure system is installed on the south main tower.The wind measure system is a ten-storied system with two kinds of wind speed measuring instruments.There is a ten-meter distance between two observation layers in vertical.In order to obtain the turbulent wind properties at the mountainous bridge site, two Yong81000 ultrasonic anemometers are used to collect wind data.Ten cup anemometers are installed in the wind measurement system to get the mean wind properties at the mountainous bridge site.The arrangement of wind measurement system is shown in Fig. 2.

Extreme wind record
The wind measurement system collected wind data from 1st May 2014 to 30th April 2016, the collected wind data was partitioned by 10 minutes, and the invalid data points were deleted.The wind data participated in analysis with an effective rate above 90 %.To ensure the integrity and consistency of the subsample, the position of deleted invalid data was filled by the interpolation method.Twenty-five times strong wind with a speed of above 10.0 m/s was recorded during the observation period.The largest recorded wind speed was on 2:25-2:35, 26th August 2014, and the wind data was given in Fig. 3.

Modal of bridge
A three-dimensional finite element model was set up on the ANSYS platform.Spatial beam4 element was used to simulate the main girder, and the beam44 element was used to simulate main towers.Link10 element with three degrees of freedom was used to simulate main cables and suspenders.Pavement and railings' stiffness contributions were neglected, and their lumped masses were equally distributed to the main girder with mass21 element.Combin14 element was selected to simulate damper.The numerical simulation of the Cuntan suspension bridge was constructed according to the data shown in 2.1 part of this paper.According to the design, material properties and elements' characteristics were added in the model and given in Table 1.Ernst equation of equivalent modulus of elasticity was used to calculate linearized stiffness of back cables' nonlinearity stiffness.According to the bridge design, the deck and main towers were coupled in three degrees of freedom, including the vertical displacement, transverse displacement, and rotation around longitudinal direction.Two main cables were fixed on the top of towers, and the bottoms of main cables were fixed at the bases.Two main towers were also fixed at the bases without considering the soil-pile structure interaction.The bridge was dispersed into 818 elements and 721 nodes.Fig. 4 gave the three-dimensional finite element model of the Cuntan Yangtze Bridge.

Specified power spectrum of turbulent wind
To bridge engineering, it is of great significance to investigate the wind properties acting on the long-span bridge.An accurate power spectrum mode of turbulent wind can give the exact response of long-span bridge.During the past research about power spectrum of turbulent wind, it is found that the Karman power spectrum in the along-wind direction and Panofsky power spectrum in the vertical-wind direction are normalized power spectrums [33].Two power spectrums are given in Eq. ( 1), Eq. ( 2) and Eq. ( 3): For the Karman power spectrum in the along-wind direction: (1) For the Panofsky power spectrum in the vertical-wind direction: where is the auto-power spectral density of along-wind turbulence, and is the auto-power spectral density of vertical turbulence. is the wind speed frequency; is the height of the wind speed. is the mean wind speed at 10 m height.is the Moning coordinate.and are the wind friction speeds.

Measured power spectrum of turbulent wind
The Hamming window and piecewise smoothing method were used to analyze the wind speed data in order to get the power spectrum of turbulent wind.The Hamming window was adopted to reduce the signal leakage in the frequency domain.The piecewise smoothing method was used to reduce the random error of spectrum estimates.
Measured power spectrums in the vertical-wind direction and along-wind direction were given in Fig. 5.It was found that the measured power spectrum lines did not match well with the Karman and Panofsky power spectrums.The relatively apparent deviations can be witnessed both in the along-wind and vertical-wind directions.Measured power spectrum line in the along-wind direction was a bit lower than the Karman power spectrum in a lower frequency and on the contrary in a higher frequency.Measured power spectrum in the vertical-wind direction was a bit higher than the Panofsky power spectrum in a lower frequency and on the contrary in a higher frequency.The relatively apparent deviations can be witnessed in the along-wind and vertical-wind directions.
It was witnessed that measured power spectrums had an apparent deviation with the specified ones.Therefore, fitted power spectrums were put forward to analyze the buffeting performance of suspension bridge.The Karman power spectrum and Panofsky power spectrum were used as the fundamental formulae.The fitted formulae can be expressed in Eq. ( 4) and Eq. ( 5):

Dynamic characteristic analysis
Analyses were performed using modal approach, solving the bridge dynamics equations of equilibrium in generalized coordinates through a state space transformation based on the first twenty modes.The vibration mode of the Cuntan Yangtze Bridge was shown in Fig. 6.The first twenty natural frequencies were listed in Table 2. From Fig. 6 and Table 2, it can be summarized that: The fundamental frequency is 0.1122 Hz, corresponding to the symmetric lateral vibration of girder, and the basic period of structure is short.
The second frequency is 0.1162 Hz, corresponding to anti-symmetric vertical vibration of the main girder.This vibration mode of second frequency conforms to the basic rule of the dynamic performance of the flexible structure.In the first twenty-order vibration mode, vibrations of main cable and girder are taken as the principal vibration modes, without appearing vibrations of main tower.It is probably due to that stiffness of main tower is considerably larger than stiffness of main cable and girder.The distribution range of first twenty-modal frequency is 0.1122 Hz-0.4702Hz, and distribution of overall frequency is relatively wide.

Measurement of aerostatic force coefficients
The aeroelastic effects on the bridge deck were measured by aerostatic force experimental coefficients.Fig. 7  Wind tunnel test requires that the section model is similar to the real bridge in geometric dimensions, as well as frequency and damping ratio.But actually, the section model of bridge cannot have similarity with a prototype model in all aspects.Deviation is allowable in the wind tunnel test.Allowable damping ratio deviation should be controlled in less than 10.0 %.From Table 3, it can be calculated that the deviation is 4.30 % on vertical bending damping ratio and 3.80 % on torsion damping ratio, and other parameters keep the same as the prototype modal.As a consequence, the results of experiment are effective.
The results of aerostatic force coefficients at different wind attack angles were shown in Fig. 8. Aerostatic force coefficients and their derivatives of , and were used to calculate the buffeting loads.Buffeting loads at per-unit span length are expressed in Eq. ( 6), Eq. ( 7) and Eq. ( 8): where , and represent lift, drag and moment from buffeting effects, respectively.( ), ( ) and ( ) are dimensionless lift, drag and moment coefficients at a specified wind attack angle, and wind attack angle is taken as 0° in this buffeting performance analysis.( ) and ( ) are wind velocity fluctuations in the horizontal and vertical directions, respectively., and are aerodynamic admittances.

Calculation of buffeting force
Wind loads acting on the bridge deck can be divided into two parts, including self-excited aerodynamic loads and turbulent loads.Buffeting loads can be expressed in Eq. ( 9), Eq. ( 10) and Eq. ( 11) [25,[33][34][35][36][37]: where , and represent lift, drag and moment from self-excited motions, respectively.Self-excited and buffeting forces are shown schematically in Fig. 9.The lift, drag and moment of the self-excited part of the force at the per-unit length of the deck can be expressed in Eq. ( 12), Eq. ( 13) and Eq. ( 14): where represents air density.represents the bridge deck width.represents the reduced frequency, and = ⁄ , represents the circular frequency of the bridge motion.* , * , * ( = 1, 2, 3, 4, 5, 6) represent flutter derivatives of the bridge deck measured through the wind tunnel tests in the Southwest Jiaotong University.The identified flutter derivatives were inputted as the coefficients in the aeroelastic stiffness and damping matrices of Matrix27.represent the mean wind velocity.ℎ, and represent vertical, horizontal and rotational displacements of the bridge deck, respectively.
The buffeting forces were applied to the finest element model by APDL programming.The buffeting forces acting on the main towers were neglected.

Buffeting performance in time domain and root mean square (RMS) responses
Based on the wind database, the mean wind speed value of the whole analysis period was taken as 19.3 m/s and the analysis period was 600 s.The buffeting responses of the Cuntan Yangze Bridge were calculated using the Kaimal power spectrum in the along-wind direction, the Panofsky power spectrum in the vertical-wind direction (which was named as specified power spectrum) and the fitted power spectrums.The buffeting responses of the bridge deck were critical in evaluating comfort during buffeting, determining the dynamic wind loads, and calculating the fatigue reliability of steel box girders.The girder of midspan was constrained at least by comparing with other sections of girder.Therefore, buffeting responses at the girder midspan were taken as a typical section of girder.Buffeting responses with different power spectrum were shown in Fig. 11-Fig.13.
Comparisons of RMS buffeting displacements at the midspan and 1/4 span of girder were given in Table 4. From the deviation in Table 4, the average deviation in three directions was 16.67 % at the midspan and 10.67 %, 10.22 % at 1/4 span.Therefore, it can be summarized that RMS buffeting displacements at the midspan were affected more apparently than 1/4 span in three directions.Deviations in vertical and transverse direction were 3.53 % and 3.20 %, 2.03 % and YANG YANG, YAO GANG, FUJIA WEI, WEIHE QIN 2.44 % at 1/4 span.Effects of power spectrum in the vertical and transverse directions were almost the same at 1/4 span.Deviations in the vertical and transverse directions were 11.16 % and 5.26 % at midspan.The deviation in the vertical direction was more than double in the transverse direction at midspan.Deviations in the torsional direction respectively were 33.57%, 25.25 % and 26.20 %.Effects of power spectrum on torsional direction were most obvious in three directions.In order to summarize the effect of power spectrum on the main girder, RMS responses at YANG YANG, YAO GANG, FUJIA WEI, WEIHE QIN critical parts were displayed in Fig. 14.
From RMS responses in Fig. 14, the following conclusions can be summarized: RMS in three directions showed a symmetrically tendency on the whole centering on the midspan.In transverse and torsional direction, the RMS responses increased moderately from endspan to midspan.In vertical direction, the RMS responses varied slowly from 1/4 span to midspan but increased at a quick speed from endspan to 1/4 span.The max RMS responses in transverse and torsional direction both appeared in the midspan, while max RMS responses in vertical direction appeared in at 1/4 span.Both of fitted and specified power spectrum RMS responses showed the same tendency in three directions, evidencing that numerical modeling accuracy of the suspension bridge was credible.RMS responses with fitted power spectrum were obvious larger than the specified ones at midspan in three directions, and discrepancy in torsional direction was the largest.

Buffeting performance in frequency domain
The buffeting performance in the time domain changes from time to time.For analysis convenience, the buffeting performances in the time domain was chosen for the auto-power spectrum analysis by fast Fourier transform (FFT).The size of FFT was 1024, and a nonlinear least-squares fitting technique was used to reduce the random error of spectral estimates.In order to find out the contribution of vibration mode to buffeting responses, vertical, horizontal and torsional direction responses were transformed in the frequency domain.Power spectrum density with specified power spectrum and fitted power spectrum of midspan were given in Fig. 15, including vertical, horizontal and torsional directions.
The following conclusions can be summarized from Fig. 15: Specified power spectrum and fitted power spectrum were coincidental in the tendency and location of peak points.For a frequency range higher than 0.40 Hz, the agreement was positive owing to that the deviation of along bridge distance (m) specified power spectrum fitted power spectrum wind power spectrum was exceedingly small.The specified power spectrum was consistent with fitted power spectrum better in the high frequency region, and there was a distinct discrepancy in the low frequency region.The energy of wind was concentrated within the range of 0.02 Hz-1.0 Hz, and the energy contribution of low frequency was more significant.The first peak point corresponded to the forth mode in the vertical direction, and that was the third symmetric vertical vibration of the main girder.The first peak point corresponded to the sixth mode in the transverse direction, and that was the second anti-symmetric lateral vibration of the main girder.The first peak point corresponded to the fourteenth mode in the torsional direction, and that was the first symmetric torsional bending mode.This implied that the contribution of the first mode was not always the largest.

Conclusions
A comprehensive research was carried out to study the buffeting performance of long-span suspension bridge in a mountainous region.The buffeting responses at different critical sections were displayed.The effects of wind power spectrum on buffeting responses were analyzed via displacement and power spectrum.Contribution of frequency on buffeting responses was revealed via the power spectrum.Conclusions can be summarized as follows: fitted power spectrum specified power spectrum 1) Fitted power spectrum in extremely strong wind has a significant difference from the specified power spectrum, and has great effects on the buffeting performance.As to this long-span suspension bridge, the basic frequency is high, and the overall mode frequency has a relative wide distribution.
2) The RMS responses with the fitted power spectrum overcome the responses with the specified power spectrum, and this phenomenon is outstanding in the torsional direction.Nonetheless, previous studies have shown the responses with the fitted power spectrum overcome the responses with the specified power spectrum in most of regions.Therefore, a particular attention should be paid to the buffeting performance of long-span bridges in the mountainous region.
3) As to the overall tendency and location of peak points, the buffeting performances in the frequency domain are coincidental.However, two kinds of power spectrum are consistent better in the high frequency region, and there is a distinct discrepancy in the low frequency region.For this bridge, the first mode contribution is not always the largest.

5 .
parameters and calculated by the anonymous function in MATLAB.Value of , , and are 30.21,5.32, 28.16, 6.56 respectively.a) Along-wind direction b) Vertical-wind direction Fig.Comparison of measured and specified power spectrum

6 .
a) First order vibration mode b) Second order vibration mode Fig. Vibration mode of Cuntan Yangtze bridge

7 .
showed the experimental aerostatic force in the wind tunnel test.The bridge section model was made of wood.Pedestrian guardrails, anti-collision guardrails and center separation band guardrails were manufactured in plastic plates by machine.The section model was 2.1 m in length, 0.7 m in width, and 0.0583 m in height with a scalar of 1/60 to the real bridge.It was fixed by eight springs on supports, providing a vibration system with two degrees of freedom which can simulate vertical and torsional vibration.The distance of springs was 108.0 cm.Two laser displacement sensors placed under the section model were used in testing the displacement of bridge.The distance between laser displacement sensors was 40.0 cm.a) 0° wind attack angle b) +12° wind attack angle Fig. Test of aerostatic force coefficients

Fig. 9 .
Fig. 9. Aerodynamic forces on bridge deck PERFORMANCE OF LONG-SPAN SUSPENSION BRIDGE BASED ON MEASURED WIND DATA IN A MOUNTAINOUS REGION.

11 . 12 . 13 .
a) Specified power spectrum b) Fitted power spectrums Fig. Vertical buffeting responses with different power spectrums a) Specified power spectrum b) Fitted power spectrums Fig. Horizontal buffeting responses with different power spectrums a) Specified power spectrum b) Fitted power spectrums Fig. Torsional buffeting responses with different power spectrums PERFORMANCE OF LONG-SPAN SUSPENSION BRIDGE BASED ON MEASURED WIND DATA IN A MOUNTAINOUS REGION.

Fig. 14 .
RMS responses along bridge distance

Table 1 .
Type of material

Table 2 .
Frequencies of Cuntan Yangtze bridge

Table 3 .
Design parameters of section mode