Evaluation of Wind-Induced Response Bounds of High-Rise Buildings Based on a Nonrandom Interval Analysis Method

+e traditional wind-induced response analysis of high-rise buildings conventionally considers the wind load as a stationary stochastic process. +at is, for a certain wind direction angle, the reference wind speed (usually refers to the mean wind speed at the building height) is assumed to be a constant corresponding to a certain return period. Combined with the recorded data in wind tunnel test, the structural response can be computed using the random vibration theory. However, in the actual typhoon process, the average wind speed is usually time-variant. +is paper combines the interval process model and the nonrandom vibration analysis method with the wind tunnel test and proposes a method for estimating the response boundary of the high-rise buildings under nonstationary wind loads. With the given upper and lower bounds of time-variant wind excitation, this method can provide an effective calculation tool for estimating wind-induced vibration bounds for high-rise buildings under nonstationary wind load. +e Guangzhou East tower, which is 530m high and the highest supertall building in Guangzhou, China, was taken as an example to show the effectiveness of the method. +e obtained boundary response can help disaster prevention and control during the passage of typhoons.


Introduction
e conventional methods of wind-induced vibration of tall buildings are based on the basic assumption that wind loads are stationary random processes [1][2][3][4].at is, assuming the reference wind speed of a certain wind direction is a constant.Combined with wind tunnel test data, this assumption can be used to calculate structural response and provide the peak response for a structural design [5].An actual typhoon has obvious nonstationary characteristics.e average wind speed at the reference height of a building is time-varying rather than constant.
e conventional structural windinduced response analysis framework cannot effectively predict structural response at each moment when a typhoon passes.In actual operation, this nonrandom vibration analysis method is used to predict typhoon wind speed from a small sample and obtain the upper and lower boundaries of the typhoon.Combined with high-frequency force balance (HFFB) wind tunnel test data, the boundary values of structural response at each time instant can be predicted, and the boundary response of a building can be estimated to guide disaster prevention and control.
During a real typhoon, average wind speed at the reference height of a building is time-variant.Researchers have developed methods to estimate the response boundaries of a system subjected to nonstationary loads.Qiu and Elishakoff proposed a subinterval perturbation method for estimating the static displacement bound of structures with interval parameters [6].Chen et al. [7] proposed a method to solve the static displacement problem with interval parameter structure uncertainty.Qiu and Wang used parametric perturbation methods to evaluate the dynamic response of structures with uncertain but bounded parameters [8].Moens and Vandepitte proposed a method of interval sensitivity evaluation in the framework of general numerical interval analysis [9].Verhaeghe et al. [10] proposed the concept of interval domains, considering static analysis of a structure with interval uncertainty using the finite element method.With the development of interval convex model theory, Jiang et al. [11,12] proposed a nonrandom vibration analysis method for nonprobabilistic convex models.It is only necessary to know the boundary excitation information at an arbitrary time point rather than the exact probability distribution.
e solution system is a responsive interval under uncertain excitation. is method reduces the dependence on large sample sizes and considers the correlation between interval variables and provides a new nonlinear interval programming method, which can be used to deal with uncertain optimization problems when there is a dependency relationship among interval variables.Gao et al. [13] proposed a mixed probability and interval reliability analysis method to solve engineering problems.Considering that uncertain factors widely exist in practical engineering, Wang et al. [14][15][16] proposed a novel methodology of reliability-based multidisciplinary design optimization under hybrid interval and fuzzy uncertainties, and a dimension-wise method and its improvement for multidisciplinary interval uncertainty analysis, and the methods of hybrid time-variant reliability estimation for active control structures under aleatory and epistemic uncertainties.
ese methods have important value in solving the uncertainty problem in practical engineering.
e HFFB technique is one of the most widely used techniques to evaluate effects of wind-induced vibration of tall buildings [17][18][19][20][21][22]. e nonrandom vibration theory of the interval process model is combined with HFFB wind tunnel tests to predict the boundary response, including displacement, velocity, and acceleration response, of highrise buildings subjected to nonstationary wind loads.
rough a case analysis of Guangzhou East Tower (GZET), the feasibility of the method is examined.

Interval Process Model. For an uncertain process
} is called an interval process model, denoted as X I (t) [23].For the interval process model X I (t), X U (t), and X L (t) are the upper and lower bound functions, respectively.e middle point function of X I (t) and X C (t), is defined as e half width function X W (t) of the interval process model X I (t) is defined as e variation coefficient function η X (t) of the interval process model X I (t) is defined as For an interval process X I (t), the self-correlation function between any two interval variables X I (t i ) and X I (t j ) at different times t i and t j is defined as r XX (t i , t j ) where

􏽱
and are the values of the pth sample function at time t i and t j , respectively; and X C i and X C j are the values of the middle point function of X I (t) at time t i and t j , respectively.

Nonrandom Vibration Analysis Method.
Based on the interval process model, Jiang et al. [23] proposed a nonrandom vibration analysis method.In this method, the timevarying uncertain excitation is treated as an interval process model, and the corresponding structural response is the interval form of the upper and lower boundaries.In practical applications, this nonrandom vibration analysis method, combined with HFFB technology, can estimate the structural response boundaries of tall buildings and give nonstationary reference wind speeds.is method requires only conventional wind tunnel tests which are carried out in a stationary wind field.

Under Deterministic Excitation.
As shown in Figure 1, consider a simple single degree of freedom (SDOF) oscillator model, with mass m, spring stiffness k, damping coefficient c, and external force F(t).e motion equation can be given in the form as follows: or where z, _ z, and € z are the displacement, velocity, and acceleration response for the first mode in each direction of the system, respectively.e natural frequency ω n and damping ratio ζ are given by 2 Shock and Vibration where c c is the critical damping coefficient.In this paper, the response of Equations ( 5) or ( 6) is solved using Duhamel's integral.Herein, the piecewise linear interpolation function [24] is first used to approximately fit the sequence of F(t i ) (denoted as F i ), where the rate of change of F(t) within arbitrary period of time is supposed to be linear.In this way, Duhamel's integral can be then applied to compute the response at time t � t j , j � 1, 2, 3, . . ., n: where ΔF j � F j − F j−1 , Δt j � t j − t j−1 , and .e velocity _ z j at time t � t j can be calculated as follows: Equations ( 8) and ( 9) can be rewritten in the following form: where Rewriting Equation (10) in matrix form gives Figure 1: SDOF spring-mass-damper vibration system.

Shock and Vibration
By using either Equations ( 10) or ( 12) repeatedly, the following recurrence sequences can be obtained: Hence, there is: Equation ( 12) can be expressed as where H is a 2 × (j + 3) constant coefficient matrix, F is the vector of force sequence, and z 0 is the initial condition.

Under Uncertain Excitation.
In nonrandom vibration analysis of the interval process model, the external excitation is treated as an interval process, and the autocorrelation function is r(t i , t j ).F I (t) can be represented by an interval vector, 16)).erefore, Equation ( 15) can be rewritten by eliminating correlation between variables and transforming them into an independent interval processes F I (t): F I (t) is obtained through Cholesky decomposition, by eliminating the correlation between variables and transforming them into independent interval processes, F I (t), with a corresponding interval vector, F I .e following equation can be obtained [23]: e corresponding interval middle point vector is F C , for which the following equation can be obtained: e following equation can be obtained for the corresponding interval half width vector F W : where L is a (j + 1) × (j + 1) lower triangular matrix, [|L|] denotes the matrix whose entries are absolute values of corresponding ones in the matrix L. rough this treatment, Equation ( 16) can be rewritten as [23]: 4

Shock and Vibration
where H F and H Z are block matrices of H, with dimensions 2 × (j + 1) and 2 × 2, respectively.Denote According to Equations ( 18)-( 20), F I can be rewritten as Substituting Equations ( 22), (18), and (19) in Equation ( 20) yields is implies that the upper and lower bounds of the displacement response at time t � t j are Similarly, the upper and lower bounds of velocity response at time t � t j can be obtained as Equation ( 5) can be rewritten as Rewrite Equations ( 26) and ( 27) in matrix form: , where z I j , _ z I j and € z I j are the interval displacement response, interval velocity response, and interval acceleration response, respectively.Equation (28a) can be rewritten as Equation ( 26) indicates that the upper and lower bounds of acceleration response can be obtained as where F U (t) is the upper bound of excitation force and F L (t) is the lower bound of excitation force.

Data Process Method for HFFB Wind Tunnel Test
In the wind tunnel test, the mean wind speed profile is adjusted to match the target wind profile as where v z is the mean wind speed at height z and α is the terrain roughness exponent.A wind tunnel test is a kind of experiment that uses reduced-scale models to represent the geometrical characteristics of a real building.Based on the principle of similitude, the measured wind forces and moments at the base of the building model can be transformed to the actual wind force acting on the building prototype.To establish the time, geometry, and wind speed similarity relationship between the model and the prototype, time scale ratio λ T , geometrical scale ratio λ L , and wind speed scale ratio λ V are defined as Equations (31a)-(31c), respectively:

Shock and Vibration
where T m , L m , and V m are the time, geometrical size, and wind speed for the model, respectively.T full , L full , and V full are time, geometrical size, and wind speed for actual buildings, respectively.According to the principle of similitude, the above three satisfies Since T � 1/f, the scale ratio T of frequency can be obtained as where f m and f full are the frequencies for the model and actual building, respectively.
In conventional data analysis of a wind tunnel test, time or frequency scale is a constant for a certain wind direction; however, if a real wind field is nonstationary, the time scale is time-variant.In this study, the entire considered time is separated into a series of time segments defined as Δt (i) , which can be formulated as Conventionally, aerodynamic force and moment coefficients are employed to implement the model-toprototype transform.
e torque coefficients μ x , μ y , and μ z are defined as where m x , m y , and m z are the base bending moments measured in the direction x, y, and z, respectively, ρ is the air density, V h is the mean wind speed at the model reference height, B m is the model characteristic width, and H m is the total height of the model.rough the torque coefficients μ x , μ y , and μ z , the lower and upper bounds of the base moment of the real building can be formulated as where M Ux , M Uy , and M Uz are the upper bounds of the base bending moments of the actual structure in x, y, and z direction, respectively; M Lx , M Ly , and M Lz are the corresponding lower bounds; V U and V L are the upper and lower bounds of the time-variant wind speed; and B full and H full are the characteristic width and total height of the real buildings, respectively.A high-rise building can be examined using a cantilever structure response analysis which satisfies the motional equation as where m(z), c(z), and EI(z) are the mass, damping coefficient, and flexural rigidity per unit length at height z, respectively, and p(z, t) is the horizontal wind load at height z and time instant t. e modal mass m p , modal stiffness k p , and modal damping coefficient c p for the first mode can be computed as follows: en, the structural dynamic equation in modal coordinates can be formulated as Assuming the modal shape for the first two sway modes to be linear with height, modal force can be expressed as where M I (t) is the base bending moment, which can be measured by the HFFB wind tunnel test.e whole procedure of the proposed method for estimation of boundary response of high-rise buildings subjected to nonstationary wind load is listed below: where F U x , F U y , and F U z are the upper bounds of the modal force of the actual structure in the x, y, and z directions, respectively, and F L x , F L y , and F L z are the corresponding lower bounds.

Case Study
In this paper, the feasibility of the interval process model and nonrandom vibration analysis method [23] in practical engineering is explained based on dynamic excitation as a stationary interval process, dynamic excitation as a nonstationary interval process, and the engineering example of GZET.
4.1.F I (t) Stationary.As shown in Figure 1, consider a spring-mass-damping system where the mass m � 1, damping coefficient c � 0.2, and spring stiffness k � 1. Herein, F(t) is treated as a stationary interval process F I (t), F(t) ∈ [10,20], where the middle point function is e autocorrelation coefficient function is assumed to take the form of r(τ) � e −λ|τ| .Figures 2-4 show the boundary response of displacement, velocity, and acceleration for the case of λ � 0.2.We find that the response bounds of displacement can be roughly divided into two phases: the transient response stage and the steady-state response stage.Before 30 s, the response bounds exhibit obvious oscillation.is kind of oscillation tends to reduce with time until a steady state can be reached.Figures 5-7 show the stationary interval process autocorrelation coefficient function r(τ) � e −λ|τ| .Cases with different λ values are considered to observe the effect of λ on the obtained boundary response.As can be seen from the results, as λ increases, the interval of the system's steady response phase will increase, indicating that the self-correlation function of the time-varying force, F I (t), may have a significant impact on dynamic responses, such as displacement and velocity, when the incentive interval is not changed.

F I (t) Is Nonstationary.
is case also focuses on the damped SDOF system in Figure 1, and the values of relevant parameters are the same as those in Section 4.1.e difference is F I (t).It is a nonstationary interval process, as shown in Figure 8, where the middle point function is F c (t) � 16 + 2 sin t and half width function is F r (t) � 2 + sin t.  show the displacement response, velocity response, and acceleration response of the nonstationary interval process.We find that because F I (t) exhibits a variation of periodicity, the response also exhibits a variation of periodicity, and the periodicity changes are very similar.
Figures 12-14 show the stationary interval process autocorrelation coefficient function r(τ) � e −λ|τ| .Cases with different λ values are considered to observe the effect of λ on the obtained boundary response.As can be seen from the results, as λ increases, the interval of the system's steady response phase will increase, indicating that the self-correlation function of the time-varying force, F I (t), may have a significant impact on dynamic responses, such as displacement and velocity, when the incentive interval is not changed.

4.3.
e GZET Case Study.In this section, a supertall building, GZET, is used as an example to show the effectiveness of the proposed method.GZET, which is 530 m high with 112 stories, is in Guangzhou, China.Currently it is the tallest building in Guangzhou and the second tallest building in South China.
An HFFB wind tunnel test of GZET was carried out in the wind tunnel laboratory of South China University of Technology.e experimental section of the wind tunnel laboratory is 5 m wide, 3 m high, and 24 m long.e wind speed of the test section can be adjusted continuously, and the maximum wind speed can reach 30 m/s.A rigid model with a geometric length scale of 1 : 500 was made to represent GZET.e model was made of balsa and foam and was mounted on a metal slab which was rigidly connected to an HFFB measurement system.Models of surrounding buildings that might affect the wind's effects on the tested building were also made and mounted on the ground of the wind tunnel to simulate the surrounding conditions, as shown in Figure 15.A Pitot tube was installed at the height of the building model (1.06 m above ground) to measure the reference wind speed.e wind direction was defined as an angle, β, from the east along an anticlockwise direction, varying from 0 °to 360 °with an increment angle of 10 °.
Before the formal HFFB wind tunnel test was carried out, the wind profile in the wind tunnel had been adjusted to match the ESDU's wind profile with the exponential parameter α � 0.22, as shown in Figure 16.In the wind tunnel test, the sampling frequency is 400 Hz with a sampling length of 40960.e natural frequency of GZET is 0.128 and 0.132 Hz  for the first two sway modes and 0.360 Hz for the first torsional mode.According to the actual engineering situation, the correlation coefficient function of wind load is obtained directly from wind tunnel test to compute the response boudary of the structure.
In the wind-induced vibration response analysis, using the wind speed of a certain typhoon provided by the meteorological department, 3,600 samples were taken to measure the typhoon wind speed, as shown in Figure 17, and the upper and lower envelope diagrams thereof were plotted according to Figure 17     Shock and Vibration

Conclusions
Based on the HFFB technique and interval process model analysis theory, this paper proposes a new method for estimating the response boundary of high-rise building structures based on nonstationary wind loads.With a given nonstationary reference wind speed, this method can provide the upper and lower bounds for displacement, velocity, and acceleration response for different time instants.With the predicted wind information of a typhoon, it is easy to know the possible largest response for different time segments, which helps the owner of a building take necessary measures to prevent a typhoon disaster and reduce damage.is new method will help  12 Shock and Vibration researchers and engineers prevent and control typhoon disasters in high-rise buildings.

Data Availability
e wind tunnel test data (in * .txtfiles) to support the findings of this study were supplied by Guangzhou University with license and so cannot be made freely available.
Request for access to these data should be made to the corresponding author after publication of this article.

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

( 1 )
Get the upper and lower bounds of the reference wind speed conventionally from meteorological agencies.(2) Measure the base forces and torques through a wind tunnel test.(3) Using the base overturning moments, compute the base overturning moment coefficient according to Equations (35a)-(35c).(4) Using Equations (36a)-(36c), compute the upper and lower bounds of structural modal force by the following equations: 6 Shock and Vibration

10
Shock and Vibration acceleration responses of the structure, respectively.To better observe the changing trends of response at each time, the envelopes of structure displacement, velocity, and acceleration responses were, respectively, fitted.In Figures19-21, the blue curve is the structural boundary response obtained by the interval analysis method, and the red curve is the envelope of the boundary response method which is more intuitive for applications.e envelope diagram can reflect the response at each moment more intuitively and clearly.

10 Figure 13 :
Figure 13: Dynamic velocity response under different values of parameter λ.

Figure 14 :
Figure 14: Dynamic acceleration response under different values of parameter λ.

Figure 15 :
Figure 15: Building model of the GZET in the wind tunnel.