Model shape correction method for highfrequency force balance technique
Ruohong Zhao^{1} , An Xu^{2} , Weixing Sun^{3} , Xiaohua Lan^{4}
^{1, 2, 3, 4}Engineering Technology Research and Development Center for Structural Wind Resistance and Health Monitoring in Guangdong Province, Guangzhou University, Guangzhou, China
^{2}Corresponding author
Journal of Vibroengineering, Vol. 19, Issue 3, 2017, p. 16651679.
https://doi.org/10.21595/jve.2017.17938
Received 3 November 2016; received in revised form 14 March 2017; accepted 17 March 2017; published 15 May 2017
JVE Conferences
This study proposed the mode shape correction formulae for the high frequency force balance technique. The presented formulae respectively give the mode shape for translational modes and torsional modes, meaning that there is no need to distinguish the traditionally concerned along wind case and crosswind case anymore. This treatment makes the mode shape correction much simpler especially when wind blows to a building obliquely. Comparison with some experimental results of other researchers shows the reliability of the proposed formulae. The mode coupling effect is also analyzed and discussed through an actual supertall building, the Guangzhou East Tower. The results show that the mode shape coupling effect is significant and cannot be ignored.
Keywords: tall building, wind effect, highfrequency force balance.
1. Introduction
The highfrequency force balance technique (HFFB) was firstly introduced by Tschanz and Davenport (1983) for the measurements of wind effects on tall buildings [1], and now it has become one of the standard windtunnel methods by which overall wind loads and responses such as accelerations, displacements and velocities are determined for tall buildings at the design stage [24]. By assuming the fundamental mode shape of a supertall building being linear with height, the HFFB technique allows users easily obtaining the modal force for the first three modes from the measured horizontal force, overturning moments and torsion at the base of building model. After that, the windinduced structural response can be computed. For its usability and comparatively lower test costs, the HFFB technique has now become one the most widelyused techniques for investigating the wind effects on tall buildings.
Researchers have noticed that the linear modal shape assumption in the classic HFFB technique has now been facing a challenge: the structural systems of modern supertall buildings are becoming more and more complex due to their innovative architectural and structural design; therefore, the mode shape of these modern supertall buildings might not be approximately linear with the height anymore. In other words, the linear mode shape assumption of the classic HFFB technique may not be suitable for these supertall buildings, and consequently, the classic HFFBbased data analysis method needs the socalled mode shape correction. Moreover, the mode coupling effect is not taken into consideration in the classic HFFB technique, however, if the geometrical axes and modal vibration axes have an obvious drift angle, just as the example illustrated later, the mode coupling effect would be significant. Obviously, the classic HFFB technique doesn’t consider the mode coupling effect. Therefore, there is a need to improve the HFFB technique to satisfy the requirement of wind effect analysis for modern supertall buildings.
There were many efforts to find out an effective solution to the aforementioned problems. Holmes et al. (1987) investigated some mode shape correction methods and found that these methods are conservative [5]. He derived the mode shape correction factor for two extreme cases: low correlation case by assuming no correlation between the fluctuating sectional forces at any pair of different height levels on the building, and conversely, the full correlation case. Holmes also proposed an intermediate correction factor in terms of mode shape exponent between the low and high correlation limits. Boggs and Peterka (1989) introduced a load correction factor, which is essentially related to the base bending moment. However, because the base moment factor is not sensitive to the mode shape, this mode shape correction factor is not suitable for estimating wind load effects other than the base bending moment. Besides, Boggs and Peterka only considered the full correlation limit, but assumed that the fluctuating forces varied with the height as a power law [6]. Vickery et al. (1985) made some direct measurements of the mode shape correction, using time histories of sectional forces based on pressure measurements [7]. Xu and Kwok (1993) modified the low and high limits of Holmes et al. (1987) to account for other spectral density variations with the height [8], in that study, the mode correction factors were given in different formulae for longitudinal and lateral cases. Zhou et al. (2002) proposed a series of mode correction factors for different kinds of wind load effects, including generalized wind loads, base shear response and base bending moment responses, displacement response, and so on. By separating the total structural response into background and resonant parts, these mode factor computation equations are quite complicated and cannot consider the model coupling effect. In summary, some of the existing mode correction methods used different mode correlation factors for the cases of longitudinal wind direction and lateral wind direction. While in some cases, wind may blow to the building with an oblique angle, and it’s hard to determine such a case belongs to the longitudinal or lateral wind cases. On the other hand, to take the mode coupling effect into consideration, the correlation between the background response and resonant response should not be neglected.
In this study, a uniform mode shape correction factors are proposed for all the translational mode cases, meaning that there is no need to differ longitudinal or lateral wind cases anymore, and consequently, the proposed method is easier to be used than the exiting methods. The proposed mode shape correction method can not only consider the actual mode shapes of buildings rather than adopting the linear mode shape assumption, but also provides an approach to take the mode coupling effect into consideration. A supertall building, the Guangzhou East Tower, is taken as an example to show the validity of the proposed formulae.
2. Methodology
2.1. Mode shape correction factor for translational modes
In general, the first sway mode in the horizontal direction can be formulated as:
where $z$ is the height of current place, $h$ is the total height of the building and $\beta $ is the exponent depending on the actual mode shape. The power spectral density of the generalized force can be expressed as:
where $n$ denotes frequency, ${C}_{0}\left({z}_{1},{z}_{2},n\right)$ is the cospectrum of fluctuating wind forces which can be expressed as:
where ${S}_{f}\left(n\right)$ is the fluctuating wind force spectrum at the standard height, $R\left({z}_{1}{z}_{2},n\right)$ is the coherence function in terms of the distance between ${z}_{1}$ and ${z}_{2}$, and $W\left(z\right)$ is the adjusting factor at the height of $z$ as:
where $\alpha $ is the power law for the terrain, and $\nu $ is a constant for different wind case. According to the suggestion of Boggs and Peterka (1989), the value of $\nu $ is taken as 1 for the longitudinal wind case and 2 for lateral or torsional wind cases [6]. The mode shape correction factor is defined as the ratio of the generalized force spectrum for the real mode shape to that for the linear mode shape as:
For a low correction level, the correlation of wind loads falls off rapidly with the increasing distance $\left{z}_{1}{z}_{2}\right$, therefore, the crosscorrelation function $R\left({z}_{1}{z}_{2},n\right)$ can be assumed to take the form as:
For along wind responses:
For crosswind or torsional responses:
For a high correlation level, a full correlation is assumed for the crosscorrelation over the distance between ${z}_{1}$and ${z}_{2}$:
For along wind response:
For crosswind response:
Fig. 1. Mode shape correction factor for high and low correlation assumptions
a) Along wind
b) Crosswind
Fig. 1 shows the relation between the mode shape correction factor and the mode shape exponent for different terrains. The mode shape correction method described above has at least two problems. Firstly, the mode shape correction formulae given by Eqs. (711) are intended for two extreme cases: full correlation or nocorrelation. Actually, the correlation of the wind forces on different levels of tall buildings usually decreases over the distance in a power law manner [9]. Neither the full correlation assumption nor the nocorrelation assumption is consistent with the real situation. On the other hand, when investigating the wind effects on an actual tall building, the full wind angle cases should be considered. In wind tunnel tests, the data for 24 or 36 uniformly distributed wind direction angles are commonly measured and then the windinduced vibrations are analyzed. This means that besides the along wind and crosswind angles, other cases that wind blows obliquely to the building need to be analyzed as well. It’s difficult to discern if such cases belong to the along wind cases or crosswind cases. Consequently, Eqs. (711) cannot be applied directly to these cases.
In this study, a uniform correction formula for translational vibration analysis of any wind angles is proposed by referring the form of Eq. (11) as:
The coefficients ${a}_{1}$, ${a}_{2}$, ${a}_{3}$ and ${a}_{4}$ must satisfy Eq. (1) for any value of $\alpha $, ${\eta}^{2}$ should be equal to 1 when $\beta $ equals to 1, and Eq. (2) the correction factor given by Eq. (12) falls intermediately in the correction factor curves given by Eq. (11). Obviously, the former can be automatically satisfied, and the second condition requires repeatedly adjusting the four coefficients to make the computed correction factor falling between the curves of high and low correlation cases for any values of $\alpha $ and $\beta $. In this regard, the four coefficients can only be obtained empirically, any coefficients that satisfy the aforementioned two conditions are acceptable. The final form of Eq. (12) we suggest for practical use is:
Fig. 2 shows the comparison of the proposed correction factors with those given by Eqs. (711). It can be seen that for both the along wind and crosswind cases, the computed correction factors are located between the curves determined by high correlation and low correlation assumptions. Holmes (1987) proposed an intermediate function without considering the category of terrain as:
Fig. 2. Proposed mode shape correction factors for different categories of terrains
a)$\alpha =$ 0.12
b)$\alpha =$ 0.15
c)$\alpha =$ 0.22
d)$\alpha =$ 0.30
By assuming the wind flow field being uniform over height, a comparison is made between the correction factors proposed by Holmes and this study, as shown in Fig. 3. It can be seen that in the practical range of $\text{0.5}\le \alpha \le \text{2}$, the correction factor proposed in this study shows a good agreement with the results of Holmes. Moreover, the proposed method in this study can take the power law of terrain into consideration.
Fig. 3. Comparison of correction factor proposed in this study and by Holmes in 1987
Based on the mode shape correction factor for the generalized forces, the correction factor of the generalized acceleration response can be easily deduced by assuming the mass of each floor being constant, as shown in Eq. (15):
In which:
where ${m}^{*}$ and ${m}_{1}^{*}$ are the modal masses for a real mode shape and assumed linear mode shape, respectively. In a similar way, the mode shape correction factor for the resonant overturning base moment response can be formulated as:
If $\alpha $ and $\beta $ are taken as 0.22 and 1.5, respectively, the value of ${\eta}_{M}$ is 0.94, indicating that the computed overturning base moment response is not sensitive to the mode shape correction. In other words, the mode shape correction can be ignored when computing the windinduced overturning base moment response.
2.2. Mode shape correction factor for torsional modes
Windinduced torsional forces acted on tall buildings are caused by unsymmetrical wind pressure distributed on the building surface. In a similar way as before, the correction factor for the generalized torque can be deduced as
where ${S}_{T}\left(n\right)$ and ${S}_{{T}_{1}}\left(n\right)$ are the generalized torque for the real mode shape and assumed mode shape of $\beta =\text{0}$, respectively. For a low correlation level, the correction factor for the generalized torque can be deduced as:
For a high correlation level, the corresponding correction factor is:
Fig. 4 gives the torsional mode shape correction factors for the high and low correlation cases.
Fig. 4. Torsional mode shape correction factors for high and low correlation cases
Xu et al. (1993) suggested that the mode shape correction factors for the low correlation case can be used for practical applications [8], but Fig. 4 clearly shows that such a treatment would overestimate the mode correction factor. As before, here an intermediate mode correction factor is proposed for the generalized force in the torsional direction as Eq. (20):
and the corresponding plots are shown in Fig. 5.
It can be seen from Fig. 5 that Eq. (20) is capable of giving appropriate corrections of torsional modes for different terrains. If $\alpha $ is taken as zero, referring to a uniform flow field case, the mode shape correction factor given by Eq. (20) is exactly the same as suggested of Holmes. Li et al. [10] proposed a formula for estimating the mode shape correction factor of torsional modes as:
Table 1. Comparison of torsional mode shape correction factors
$\alpha =0$

Study of Li

This study

Experimental results
of Tallin and Ellingwood

$\beta $

${\eta}_{t}^{2}=\frac{15\alpha \beta +6}{15\alpha +11\beta +6}$

${\eta}_{t}^{2}=\frac{2\alpha +1}{2\alpha +2\beta +1}$


0.0

1.0

1.0

–

0.75

0.37

0.40

–

1.0

0.29

0.33

0.32

1.25

0.24

0.29

–

1.5

0.2

0.25

0.26

Tallin and Ellingwood [11] in 1985 carried out an experimental study on the value of mode shape correction factors. By assuming that $\alpha $ equals to zero, the torsional mode shape correction factors estimated by this study and the one from the study of Li are compared to the experimental results given by Tallin and Ellingwood, as shown in Table 1. It can be seen that the torsional mode shape correction factor given by this study is closer to the experimental results, which demonstrates the applicability of Eq. (20).
Fig. 5. Mode shape correction factors for torsional generalized force
a)$\alpha =\text{0.12}$
b)$\alpha =\text{0.15}$
c)$\alpha =\text{0.22}$
d)$\alpha =\text{0.30}$
In a similar way like the procedure of obtaining the correction factor of the acceleration response for the sway modes, the correction factor for the acceleration response in torsional direction can be easily deduced as:
where:
in which ${m}^{*}\text{,}$${m}_{1}^{*}$ are the modal masses for the real mode shapes and ideal linear modes, respectively.
The mode shape correction factors for resonant base overturning moment response can subsequently be deduced based on the concept of inertia force as:
Substituting a pair of normally used value of $\alpha $ and $\beta $, for instance 0.22 and 0.15 respectively, into Eq. (23), the mode shape correction factor is obtained as equal to about 0.83, meaning that the effect of mode shape on the base torsional moment response cannot be ignored.
2.3. Mode coupling effects
As it is known to researchers and engineers, the classic HFFB theory ignores the mode correlation effect. This treatment is applicable for buildings and structures with a simple structural system and low mode coupling effects. However, as the development of building techniques, more and more modern supertall buildings with the complex structure system have coupling modes. Ignoring the mode coupling effect would lead to unreliable computational results of windinduced vibration. In this section, the mode coupling effects on structural response are discussed. Here only the first mode in each direction (two sway directions along the orthogonal horizontal axes and one torsional direction along the vertical axis) are considered.
The power spectral density of the modal force can be formulated in Eq. (24) by neglecting the item of cospectrum.
where ${\eta}_{1x}$, ${\eta}_{1y}^{2}$ and ${\eta}_{1\theta}^{2}$ are the mode shape correction factors for the modal force ${F}_{1}$ in the direction of $x$, $y$ and $\theta $, and so on. The method described in Eq. (24) is called method I in this study. In fact, the modal forces can be written in a more accurate form as:
$+2{\eta}_{1x}{\eta}_{1y}\left(\frac{1}{{h}^{2}}\right){S}_{{M}_{y}{M}_{x}}\left(n\right)+2{\eta}_{1x}{\eta}_{1\theta}\left(\frac{1}{{h}^{2}}\right){S}_{{M}_{y}M\theta}\left(n\right)+2{\eta}_{1y}{\eta}_{1\theta}\left(\frac{1}{{h}^{2}}\right){S}_{{M}_{x}M\theta}\left(n\right),$
$+2{\eta}_{2x}{\eta}_{2y}\left(\frac{1}{{h}^{2}}\right){S}_{{M}_{y}{M}_{x}}\left(n\right)+2{\eta}_{2x}{\eta}_{2\theta}\left(\frac{1}{{h}^{2}}\right){S}_{{M}_{y}M\theta}\left(n\right)+2{\eta}_{2y}{\eta}_{2\theta}\left(\frac{1}{{h}^{2}}\right){S}_{{M}_{x}M\theta}\left(n\right),$
$+2{\eta}_{3x}{\eta}_{3y}\left(\frac{1}{{h}^{2}}\right){S}_{{M}_{y}{M}_{x}}\left(n\right)+2{\eta}_{3x}{\eta}_{3\theta}\left(\frac{1}{{h}^{2}}\right){S}_{{M}_{y}M\theta}\left(n\right)+2{\eta}_{3y}{\eta}_{3\theta}\left(\frac{1}{{h}^{2}}\right){S}_{{M}_{x}M\theta}\left(n\right).$
Unlike Eq. (24), Eq. (25), called method II, took the cospectrum into consideration. Obviously, Eq. (25) took the mode shape correction and the mode coupling effect into consideration, and theoretically it is the most accurate form to calculate the spectrum of modal forces. In the next section, the Guangzhou East Tower will be taken as an illustrative example to discuss the computational results using different methods.
3. Case study
3.1. Introduction to Guangzhou east tower
In this section, an iconic supertall building, the Guangzhou East Tower (GZET), is taken as an example to show the effectiveness of the proposed mode shape correction factors and discuss the mode coupling effect. The GZET has the hight of 530 m above ground and is now the tallest building in Guangzhou, China. Guangzhou, located in the southeast coastal region of China, is one of the regions most seriously affected by strong typhoons. Hence, the GZET could endure wind effects induced by typhoons and can be regarded as ideal examples for investigating typhoon effects on supertall buildings. The architectural design of GZET adopted a socalled ‘setback’ scheme in the upper part of the building to reduce its response to strong wind. Its elevation and typical plan are illustrated in Fig. 6.
Fig. 6. Elevation and typical planes of GZET (unit: m)
A rigid model with a geometric length scale of 1:500 was made to represent the GZET. The model was made of balsa and foam, and was mounted on a metal slab which was rigidly connected to a highfrequency force balance (HFFB). The models of surrounding buildings that might affect the wind effects on the tested building were also made and mounted on the ground of the wind tunnel to simulate the surrounding conditions. Fig. 7 shows the models mounted in the wind tunnel. A Pitot tube was installed at height of GZET model (1.06 m above the wind tunnel ground) to measure the reference wind speed. The wind direction was defined as an angle $\beta $ from the east along an anticlockwise direction varied from 0° to 360° with the increment of 10°. In the wind tunnel test, the forces at the bases of the building model were measured by an ATItype HFFB. The ATItype HFFB is capable of measuring six force components including the shear forces along two horizontal directions, axial force along the vertical direction, bending moment along two horizontal directions and torque along the vertical direction. The building model of GZET was installed on the ATI force balances by ensuring that the coordinate system of the building matched the coordinate system of the balance. The sampling frequency was set as 400 Hz with the sampling length of 40960 for each wind direction.
Fig. 7. Simulated mean wind speed and turbulence intensity profile in the wind tunnel
a) Mean wind speed
b) Turbulence intensity
Fig. 8. Models in wind tunnel
3.2. Results and discussions
3.2.1. Power spectrum of modal forces
The properties of the modal force on a structure can significantly affect its response. Here the spectrum of modal forces computed by the classic HFFB method, the method considering the mode shape correction factors and the method considering the mode coupling effects using method I, is presented here in Fig. 9. It can be seen that, if the mode coupling effect is ignored, the computed power spectrum of modal forces changed a little bit when the mode shape correction factor was considered. Actually, the computed mode shape correction factor is about 0.7 for this case, meaning that it’s necessary to make the mode shape correction before computing the structural response. On the other hand, the mode shape coupling effect is significant for the wind direction of 0°, indicating that the effect of modal shape coupling effect on the structural response should not be ignored.
Fig. 10 gives the spectrum of modal force computed by the method I and II. The result computed by the classic HFFB method without considering the mode shape correction and mode coupling effect is also given as a reference. It can be seen that spectrum computation results are close to each other in the wind direction of 0°, whereas in the wind direction of 90°, the results obtained by method I and II are quite different from the spectrum computed by the classic HFFB method. This means that it’s necessary to carry out the mode shape correction and consider the mode coupling effect when estimating the windinduced response. The difference of the results given by method I and II cannot be clearly identified through Fig. 10, and it will be discussed in the next section.
Fig. 9. Comparison of computed power spectrum of modal force using different methods
a) 0°
b) 90°
Fig. 10. Comparison of results given by method I and II
a) 0°
b) 90°
3.2.2. Displacement response
Displacement is usually the basic variable in the dynamic equation. In this section, the structural displacement response is computed to investigate the influence of the cospectrum item in Eq. (25). Fig. 11 shows the displacement at the top of GZET under all wind directions. It can be seen that the displacements given by method I and II are different. In most of wind direction cases, the results given by method II, which considered the cospectrum item, are greater than that given by method II. For instance, under the wind direction of 100°, the displacement given by method II is about 13.4 % greater than that given by method I. For the case of GZET, the mode coupling effect is significant and considering the cospectrum item might tend to provide more safety to the structure.
The significant mode coupling effect for this case can be attributed to the drift angle between the geometry axes and the modal vibration direction, as shown in Fig. 12. Just as aforementioned, the axes of ATI balance matched the geometrical axes at the wind tunnel test, therefore, the wind tunnel test and the subsequent windinduced response analysis are based on the geometrical coordinate system. For common supertall buildings, the modal vibration directions for the first two modes are usually consistent with two horizontal geometrical axes. However, there is an obvious drift angle of 32.6° between the geometrical axes and the modal vibration directions for the first two modes. In other words, the modal vibrations for the first two modes project components on both the geometrical axes of $x$ and $y$. Therefore, the mode coupling effect is significant. Since the unsymmetrical structural system became more and more popular in modern supertall building construction, the mode coupling effect of these buildings might be more obvious than those built in the past, and consequently, the mode coupling effect should be taken into consideration when estimating their windinduced vibrations.
Fig. 11. Displacement response computed by method I and II
a) Direction $x$
b) Direction $y$
Fig. 12. drift angle between geometrical axes and modal vibration direction
a) First order
b) Second order
4. Conclusions
In this study, the mode shape correction method for the HFFB technique and the mode coupling effect on windinduced response were analyzed and discussed through a practical engineering example, the GZET. Major conclusions include.
Unlike the traditional methods, the uniform formulae proposed by this study directly compute the mode shape factor for translational modes and torsional modes, and there is no need to differentiate the along wind case and crosswind case anymore.
The computing method of mode shape correction factors proposed in this study has been compared with formulae by other researchers and some experimental results, it can be seen that the proposed formulae is capable of giving reliable results.
The case study shows that the mode shape correction is necessary for the response computation. For the reason that there is a drift angle between modal vibration directions and geometrical axes, the mode coupling effect is significant, and therefore, it’s necessary to make the mode shape correction and take the mode coupling effect into consideration.
In summary, this paper provides a possible way to extend the classic HFFB technique to the supertall buildings with nonlinear mode shapes. It’s helpful to analyze wind effects on supertall buildings.
Acknowledgements
The research described in this paper was financially supported by Grants from the National Science Foundation of China (51478130, 51208127), the Project of Science and Technology of Guangdong, China (2016B050501004), the Program for Changjiang Scholars and Innovative Research Team in University (IRT13057) and the Project of Science and Technology of Guangzhou, China (201707010285). The authors gratefully acknowledge such a financial support.
References
 Tschanz T., Davenport A. G. The base balance technique for the determination of dynamic wind loads. Journal of Wind Engineering and Aerodynamics, Vol. 13, 1983, p. 429439. [Publisher]
 Xie Z. N., Gu M. Simplified formulas for evaluation of windinduced interference effects among three tall buildings. Journal of Wind Engineering and Aerodynamics, Vol. 95, 2007, p. 3152. [Publisher]
 Xu A., Xie Z. N., Gu M., Wu J. R. A new method for dynamic parameters identification of a modelbalance system in highfrequency forcebalance wind tunnel tests. Journal of Vibroengineering, Vol. 5, 2015, p. 26092623. [Search CrossRef]
 Xu A., Xie Z. N., Gu M., Wu J. R. Amplitude dependency of damping of tall structures by random decrement technique. Wind and Structures, An International Journal, Vol. 21, 2015, p. 159182. [Search CrossRef]
 Holmes J. D. Mode shape corrections for dynamic response to wind. Engineering Structures, Vol. 9, 1987, p. 210212. [Publisher]
 Boggs D. W., Peterka J. A. Aerodynamic model tests of tall buildings. Journal of Engineering Mechanics, Vol. 115, 1989, p. 618635. [Publisher]
 Vickery P. J., Steckley A., Isyumov N., Vickery B. J. The effect of mode shape on the windinduced response of tall buildings. 3rd. U.S. National Conference on Wind Engineering, Lubbock, Texas, 1985. [Search CrossRef]
 Xu Y. L., Kwok K. C. S. Mode shape corrections for wind tunnel tests of tall buildings. Engineering Structures, Vol. 15, 1993, p. 387392. [Publisher]
 Xu A., Xie Z. N., Ni Z. H. Investigation of equivalent wind load on tall buildings using instantaneous pressure measurement method. China Civil Engineering Journal, Vol. 37, 2004, p. 1116. [Search CrossRef]
 Li Y. Y. Research on the Influence of Mode Shapes on WindInduced Response of Tall Buildings. Hunan University, 2014. [Search CrossRef]
 Tallin A., Ellingwood B. Analysis of torsional moments on tall buildings. Journal of Wind Engineering and Industrial Aerodynamics, Vol. 18, 1985, p. 191195. [Publisher]