Field Measurement of Wind Speeds and Wind-Induced Responses atop the Shanghai World Financial Center under Normal Climate Conditions

Field measurement data on wind velocities and wind-induced acceleration responses at the top of the 492m high Shanghai World Financial Center (SWFC) under normal climate conditions are studied. Characteristics of the mean wind speeds and turbulence intensities, gust factors, power spectral densities, and turbulence integral scales of the fluctuating wind speed are analyzed in different observation time intervals. Power spectral densities of wind-induced acceleration are also investigated. The basic natural frequencies and structural damping ratios of the building are identified based on Hilbert-Huang transform method and random decrement method. The field measurement results of wind-induced responses of the SWFC are finally compared with those from the corresponding high-frequency force balance wind tunnel test study.


Introduction
Field measurement of wind velocities and wind-induced responses is an essential aspect of wind engineering.In 1960, a field measurement study on mean wind speed at different heights and time intervals was conducted by Durst [1], who also proposed an empirical formula for different terrain categories.A long-term observation of wind speeds and wind directions with ultrasonic anemometers in Tokyo and Yokohama was performed by Choi [2] and Kato et al. [3].Several databases on local wind characteristics, for example, the offshore observation databases in Canada and the United Kingdom [4,5] and the Froya database in Norway [6,7], were established successively with field measurement data.The characteristics of the mean and fluctuating wind velocity of typhoon Muifa were analyzed by An et al. [8] through field measurement at the top of the Shanghai World Financial Center.The field measurement data of the windinduced responses at the top of a 254 m high bridge tower was researched by Larose et al. [9] and compared with the corresponding wind tunnel test study results in frequency domain.The dynamic characteristics of a 200 m high building were identified from its field measured wind-induced vibration data by Li et al. [10], and then the correlation between structural damping and acceleration responses was discussed.Field measurements of the strong wind characteristics and wind-induced responses at the 384 m high Shenzhen Diwang Mansion and the 420.5 m high Shanghai Jin Mao Tower were conducted by Xu and Zhan [11] and Li et al. [12], respectively, and relevant empirical formulas were fitted.Field measured acceleration responses under strong winds at the 391 m high Guangzhou Citic Tower and the 450 m high Guangzhou West Tower were shown by Fu et al. [13,14] and compared with the corresponding wind tunnel test.Bai et al. [15] studied the wind velocity, acceleration responses, and base bending strains on parts of a transmission tower.Yi et al. [16] and Li et al. [17] applied GPS technology to field measurement of high-rise structures and studied the displacement responses.
The previous studies on the field measured wind speeds and wind-induced responses were mainly focused on typhoon climate conditions, while those under normal climate conditions were barely found.Based on the wind speed and wind-induced acceleration responses measured by the structural health monitoring system (SHMS) installed on the SWFC, the mean and fluctuating wind speeds at a level near 500 m under normal climate conditions are thus analyzed in the present study.The maximum wind speeds and the turbulence characteristics for different observation time intervals are discussed.Frequency domain and timefrequency domain analysis of the field measured acceleration responses at the topmost part of the SWFC are also conducted.The dynamics of the building are identified with random decrement method (RDM).The wind-induced acceleration responses from the field measurement are compared with those from the corresponding high-frequency force balance (HFFB) wind tunnel test study.Some conclusions are drawn based on the results.

Field Measurement Profile
Figure 1 presents the structural health monitoring system installed in the SWFC, which mainly consists of a wind velocity monitoring system and a structural vibration monitoring system.The wind velocity monitoring system is installed on the northeast and southwest end of the top of the building at a level of 495 m high from the ground (Figure 2).As shown in Figure 3, each side has a Windmaster Pro triaxial ultrasonic anemometer (Gill Co. Ltd, UK).The distance between the two sites of instruments is 72 m.The -, -, and -axes of the anemometer face the north, west, and upward direction, respectively (Figure 4).The wind angle increases along the counterclockwise direction; thus, the measured wind velocity is 0 ∘ for the south wind and 90 ∘ for the east wind.Since they are exposed outdoors, the anemometers might be influenced by the climate like rainfall and temperature.The equipment used here has a permissible rainfall less than 300 mm/hr, and the operating temperature is between −40 and 70 ∘ .Thus, the data collection is reliable under normal climate conditions.The measuring range of the anemometer is from 0.01 m/s to 65 m/s.The maximum dynamic response frequency is 40 Hz and the actual sampling frequency is 10 Hz.All the data are collected in real time and stored by using the CR1000 data collection system (Campbell Co. Ltd, USA).In order to prevent the anemometers from being struck by lightning,  the anemometers are installed on the beam and fixed on the lightning rod with a flange on each side of the topmost part of the building.Meanwhile, the anemometers are several meters higher than the top of the building to reduce the effects of the building itself, and CFD method was used to identify its effective angles.
Computational fluid dynamics simulation of the wind field around the target building shows that the mean and fluctuating wind speeds approaching the building are influenced by the building itself.However, the influences are negligible when the included angle between the approaching wind and the -axis of the building is smaller than 15 ∘ (Figure 5), and those angles are defined as the effective wind angle.Effective wind angle ranges at the northeast and southwest ends are 120 ∘ to 150 ∘ and 300 ∘ to 330 ∘ , respectively.The following analyzed samples are those whose 10 min mean wind speeds are relatively large (over 8 m/s) and wind angles are in the effective wind angle ranges.The acceleration data are measured by the acceleration sensors installed on the 98th floor of the building (475 m).The acceleration sensors are SLJ-100 balance accelerometer produced by China Earthquake Administration.The basic frequency of the accelerometer is 100 Hz, the measurement range is ±2 g, and the frequency band is 0-100 Hz.The two horizontal normal axes of the sensors are parallel to the body axes of the building (Figure 5).Vibration signals are simultaneously measured from the three orthogonal directions by the sensors.Two different groups of windinduced responses data with a length of 10 min are selected for further analysis.In those two data groups, the mean wind speed and the corresponding mean wind angle are 16.33 m/s, 139 ∘ and 9.48 m/s, 123 ∘ , respectively.Apparently, the mean wind angles fall in the effective wind angle range at the northeast end.

Characteristics of Wind Speed
The characteristics of the measured mean and fluctuating wind speeds are analyzed in this section.The wind velocity samples are shown at first.Then, the mean wind speeds under different average time intervals are analyzed.The turbulence intensity, gust factors, and turbulence scales of the fluctuating wind speeds are also discussed.

Mean Wind Speed.
Three orthogonal components of wind speed,   (),   (), and   (), are measured simultaneously with the ultrasonic anemometers from the three directions of the -, -, and -axes, respectively.Ignoring the vertical component,   (), the mean wind speeds, , and directions, , are calculated through the vector splitting method as follows: The longitudinal and lateral fluctuating components, () and V(), are calculated with the following formula, respectively:  () =   () cos  +   () sin  − , V () = −  () sin  +   () cos . ( The analyzed samples are collected from the anemometers at the southwest and northeast ends of the SWFC rooftop (Table 1).Figure 6 shows the maximum values of mean wind speeds in 1 hr with the durations of 3 s, 10 min, and 1 hr, respectively, and the mean wind angles with a duration of 1 hr.The mean wind angle has no apparent change within the analytical time intervals.The maximum values of mean wind speeds within 1 hr decrease with the increase of the durations, but the maximum 10 min wind speeds within 1 hr are close to the mean wind speed of 1 hr.
Figures 7 and 8 show the relationship between the maximum 10 min mean wind speeds within 1 hr and the 1 hr mean wind speed, as well as that between the maximum 3 s mean wind speeds within 10 min and the 10 min mean wind speeds, respectively.The figures also show a linear fitting in which the length of the error bar in Figure 7 represents the standard deviation, and the fitting demonstrates a good result.As shown in Figure 7, the ratio of the maximum 10 min mean wind speeds within 1 hr and the 1 hr mean speed is approximately 1.04.In most cases, the standard deviation of the maximum 10 min mean wind speeds within 1 hr does not exceed 1 m/s.In Figure 8, the correlation between the   maximum 3 s mean wind speeds within 10 min (instant maximum wind speed) and the 10 min mean wind speeds can be expressed as  = 0.85 + 4.5.

Turbulence Intensity and Gust
Factor.Turbulence intensity is defined as the ratio of the standard deviation and the mean values of the approaching wind speeds.The analyzed time interval is 10 min (all average time intervals mentioned subsequently are 10 min, unless specified otherwise).The turbulence intensity is calculated as follows: where   is the standard deviation of the fluctuating wind speeds;  and V represent the longitudinal and lateral components of the fluctuating wind speeds, respectively.Gust factor is defined as the ratio of the maximum wind speed within the gust duration,   , and the 10 min mean wind speeds , which can be expressed as where max((  )) and max(V(  )) represent the maximum mean wind speeds of the longitudinal and lateral fluctuating wind speeds, respectively, within the gust duration time   () during the observation time interval.Figures 9(a) and 9(b) illustrate the relationships between turbulence intensities and the 10 min mean wind speeds as well as between the gust factors and the 10 min mean wind speeds, respectively.The longitudinal and lateral turbulence intensities both decrease as the mean wind speeds increase (Figure 9(a)).The longitudinal turbulence intensity is in [0.01, 0.30] and the mean value is 0.085.The lateral turbulence intensity is in [0.01, 0.25] and the mean value is 0.075.The ratio of those two values is   :  V = 1 : 0.88.The longitudinal turbulence intensity recommended by AIJ [18] for the same terrain category (V type) and the same height (495 m) is 0.11 and that based on the Chinese load code [19] is 0.12.As shown in Figure 9(b), the longitudinal and lateral gust factors both decrease as the 10 min mean wind speeds increase.Regulations of turbulence intensity and gust factors were also found by Fu et al. [13] through field measurement of a 391 m high-rise building in south China under typhoon condition.Then, the gust factors remain almost the same after the mean wind speed reaches approximately 20 m/s.When the gust duration time   is equal to 3 s, the longitudinal gust factor   is in [1.02, 1.59] and the mean value is 1.15; the lateral gust factor  V is in [0.01, 0.62] and the mean value is 0.17,   :  V = 1 : 0.15.
Figure 10 shows the correlation between the longitudinal gust factor and the turbulence intensity as well as between the lateral gust factor and the turbulence intensity for a gust duration time of   = 3 s.The longitudinal and lateral gust factors both increase as the turbulence intensity increases.A linear fitting is shown for the relationship between the turbulence intensity and the gust factor.Then, based on a previous study [20], a nonlinear fitting is used in accordance with the following formula: where  represents the mean time interval.Choi [21] defined  as 1 hr and suggested that  1 = 0.62 and  2 = 1.27.As shown in Figure 10, the linear fitting relations between the longitudinal gust factors and turbulence intensity as well as between the lateral gust factors, and turbulence intensity are   = 1 + 1.6   and  V = 2.4  V , respectively.The results of the linear and nonlinear fitting are almost similar, which implies that the gust factor and turbulence intensity have a linear relationship.By comparing the measurement results, the empirical formula proposed by Choi [21], as well as the linear and nonlinear fitting, the measured results of the gust factors and the results from the said empirical formula are relatively close when the turbulence intensity is smaller than Mathematical Problems in Engineering   of measured fluctuating wind speeds are in accordance with the von Karman spectra, particularly if the mean wind speeds are over 20 m/s.In the power spectral densities of the fluctuating wind speed, both longitudinal and lateral components at high frequency severely deviate from the von Karman spectra and exhibit an uptrend characteristic.Similar phenomenon that power spectral densities match well with the von Karman spectrum at low wind speeds and are uptrend at high wind speeds was also found in the research of Fu et al. [14] at a 450 m high-rise building in south China.The deviation and uptrend characteristic may be caused by the nonstationary original data.Since fast Fourier transform (FFT) is effective only for stationary data, the power spectral density determined through FFT method cannot reflect the overall characteristics of the original data.

Turbulence Integral Scale.
The turbulence integral scale can be expressed as follows: where   (0) is the value that corresponds to the power spectral density of the , V, and  components of approaching wind speed for a frequency of 0.
Figure 15 shows the relationship between turbulence integral scales and the 10 min mean wind speeds for longitudinal and lateral components of approaching wind.Longitudinal and lateral turbulence integral scales have an increasing trend as the mean wind speeds increase.This trend was also found by Li et al. [10]

Acceleration Response of the Building in Normal Climate Wind and Identification of Dynamic Parameters
The time histories and the power spectral densities of the acceleration responses for strong wind in normal climate are discussed at first in this section.Then, time-frequency analysis is conducted with the Hilbert-Huang transform (HHT) method.The dynamic characteristics of the building, that is, basic mode frequency and structural damping ratio, are identified by using the random decrement method.The dynamic characteristics identified from various methods are compared with one another and with the finite element analysis results.Finally, the acceleration responses from the field measurement are analyzed and compared with those from the corresponding wind tunnel test.

Time History of Acceleration Response. Figures 16 and 17
show the time histories of acceleration responses with particular wind speeds.The lateral wind-induced accelerations ( direction) are significantly larger than the longitudinal ones ( direction).As shown in Figure 16, both lateral and longitudinal acceleration responses are relatively slight during the first 5 min, but the lateral responses increase sharply in the last 5 min.According to wind velocity records, the increase may be caused by the difference of the wind directions between the first and last 5 min.During the first 5 min, the included angle between the approaching wind and the -axis of the building is 7 ∘ .However, it is 0 ∘ in the last 5 min when the vibration in the  direction (lateral direction) is stronger.
In another record of wind-induced acceleration responses as shown in Figure 17, there is no apparent difference between      are perceptible in the power spectral densities, while they are negligible because of the significance of contributions from the basic modes.

Time-Frequency Analysis of Acceleration
Response.Timefrequency analysis of the acceleration responses is conducted on several periods with the HHT method [22].Based on the HHT analysis, the measured record of acceleration responses corresponding to the fastest 10 min mean wind speed (16.33 m/s) is picked out for the following analysis.Empirical mode decomposition (Figures 22 and 23) is conducted at first.The intrinsic mode functions (IMF) of the first six orders have narrow ranges and high-frequency components but have narrow ranges and low-frequency components in the eleventh to thirteenth orders.Figure 24 shows that the original signals match well with the superposition of IMF in the eighth to tenth orders.Figure 25 indicates that in the low-frequency domain, the power spectral densities of the superimposed IMF in the eighth to tenth orders almost match those of the original data.The matched power spectral densities can effectively reflect the frequency domain characteristics in the first five orders of the building.When the frequency is over 2 Hz, the difference is more apparent because the high-frequency component is filtrated by the IMFs of the eighth to tenth orders.
Figure 26 shows the time-frequency series of acceleration responses.The Hilbert-Huang marginal spectra of acceleration signal in the field measurement are presented in Figure 27.The basic mode frequency in the  direction is 0.152 Hz and 0.468 Hz for the second mode.

Identification of Dynamic Structural Characteristics in
Time Domain.The random decrement method (RDM) [23], a time averaged method, is adopted to identify the natural frequencies and structural damping ratios of the basic mode in  and  directions of the building.
Figures 28 and 29 show the identified results of the basic natural frequencies and damping ratios in  and  directions versus the intercepting threshold values.The basic natural frequencies in  and  directions are 0.155 Hz and 0.153 Hz, respectively, which change hardly with the increase of the intercepting threshold values.When the intercepting thresholds reach a certain value, the damping ratios are almost constant.The values in  and  directions fluctuate at approximately 2.4% and 1.6%, respectively.The difference of those two values may be attributed to the dissimilar foundation in  and  directions, which may be caused by the construction of the Shanghai Center Building located in the southwest side of the target building.Therefore, the soilstructure interaction changes the wind-induced responses of SWFC [24] and thus affects the identification of its damping ratios.
Finite element model was established in a former seismic review of SWFC.The profile and vibration mode in  and  directions of that model are shown in Figure 30.The translational vibration modes of the first two orders according to finite element analysis are in the  and  directions, and the corresponding basic mode frequencies are 0.146 Hz and 0.153 Hz [25].Table 2 lists the results of the finite element analysis and this field measurement.The natural frequencies of the first two orders identified by the power spectral density, HHT marginal spectra, and RDM match the natural frequencies from the finite element analysis.

Comparison between Field Measured Wind-Induced Responses and Those from HFFB Wind Tunnel Test Study
The Boundary Layer Wind Tunnel Laboratory at the University of Western Ontario, Canada, once conducted a wind tunnel test on SWFC by using the HFFB technology [26].The laboratory defined two test cases as the present one and the (Figures 31 and 32).When the mean wind speeds are over 14.5 m/s, the field measured data in the  direction (longitudinal direction) properly match the tested ones, while in the  direction (lateral direction) the field measured data are slightly larger than the tested ones.The differences between those two directions may be related to three issues: the difference of the wind directions, lateral aerodynamic damping, and the interference effects from surrounding buildings.Firstly, the wind angles corresponding to the measured data are variable within the interval of 107 ∘ to 170 ∘ , but the wind angle of the wind tunnel test is stable at 135 ∘ (Table 3).The wind-induced responses in the  direction (lateral direction) are sensitive to the wind direction changes; however, the longitudinal direction is a different case.Secondly, the calculation results of the wind tunnel are obtained by applying the aerodynamic force tested from the rigid model on the structural finite element model without considering effects of the aerodynamic damping, whereas the lateral aerodynamic damping of this real super-tall building may be shown as a negative value [27].Last, numerous differences also exist between the current surroundings and the wind tunnel test case.

Conclusions
Detailed analysis of the wind velocities at 500 m upper air in Shanghai and wind-induced responses of the SWFC under normal climate conditions is conducted by collecting data through the SHMS at the top of this building.The field measured results of the wind-induced responses are compared with the corresponding HFFB wind tunnel test results.The following conclusions are drawn.
(1) Through linear fitting, the ratio of maximum 10 min mean wind speeds within 1 hr and the 1 hr mean speed is approximately 1 : 1.04.The relationship between the maximum 3 s mean wind speeds within 10 min and the 10 min mean wind speeds is expressed as  = 0.85 + 4.5.
(2) When mean speed is small, both the turbulence intensity and the gust factor decrease with the increase of the mean wind speed.When the mean wind speed is over 20 m/s, the turbulence intensity and the gust factor are not sensitive to the changes of the mean (3) The power spectral densities of the longitudinal and lateral fluctuating wind speeds match well with the von Karman spectra when the mean wind speed is small.However, when the mean wind speed is over 20 m/s, the power spectral densities are larger than the von Karman spectra.
(4) The longitudinal and lateral turbulence integral scales have an increasing trend with the increase of mean wind speed.The ratio of the longitudinal and lateral turbulence integral scale is approximately 1 : 0.8.
(5) The natural frequencies of the first two modes identified by using the power spectral density, HHT marginal spectra, and RDM match properly with the results from the finite element method.The damping ratio in  and  directions is about 2.4% and 1.6%, respectively.
(6) The measured longitudinal and lateral responses match properly with the results from the wind tunnel test when the mean wind speeds are small.When the mean wind speed is over 20 m/s, the longitudinal responses of measurement and test match each other well, while the measured lateral responses are slightly larger than the tested ones.This phenomenon may be related to the differences of the approaching wind directions, lateral aerodynamic damping, and the building surroundings.

Figure 2 :
Figure 2: Sketch of location of the anemometers.

Figure 4 :Figure 5 :
Figure 4: Definition of and -axes of the anemometer.

Figure 6 :
Figure 6: Maximum mean wind speeds in 1 hr for different durations and 1 hr mean wind directions.

Figure 8 :
Figure 8: Maximum 3 s mean wind speeds versus 10 min mean wind speeds.

Figure 9 :
Figure 9: Relationships among turbulence intensities, gust factors, and 10 min mean wind speeds.(a) Turbulence intensities versus mean wind speeds and (b) gust factors versus mean wind speeds.
in the typhoon measurement on the top of Jin Mao building, which is close to SWFC.The mean values of the longitudinal and lateral turbulence integral scales are    ≈ 200 m and   V ≈ 160 m, respectively.Thus,    :   V = 1 : 0.8.

Figure 16 :Figure 17 :
Figure 16: Time histories of acceleration responses with a 9.48 m/s wind speed.

FrequencyFigure 18 :
Figure 18: PSD of acceleration responses in  direction with 16.33 m/s wind speed.

Figure 19 :Figure 20 :
Figure 19: PSD of acceleration responses in  direction with 16.33 m/s wind speed.

Figure 21 :
Figure 21: PSD of acceleration responses in  direction with 9.48 m/s wind speed.

Figure 29 :Figure 30 :
Figure 29: Natural frequencies and damping ratios versus amplitude in  direction.

Table 1 :
Record of measured samples.

Table 3 :
10 min mean wind speeds and mean wind angles for comparison during each period.