# Dynamics and Control of High-Rise Buildings under Multidirectional Wind Loads

Authors: Aly Mousaad Aly; Alberto Zasso; Ferruccio Resta
Journal: Smart Materials Research
 (2011)

---

## Abstract

This paper presents a procedure for the response prediction and reduction in high-rise buildings under multidirectional wind loads. The procedure is applied to a very slender tall building that is instructive. The structure is exposed to both cross-wind and along-wind loads obtained from pressure measurements on a rigid model (scaled 1 : 100) that was tested in a wind tunnel with two different configurations of the surroundings. In the theoretical formulation, dynamic equations of the structure are introduced by finite element and 3D lumped mass modeling. The lateral responses of the building in the two directions are controlled at the same time using tuned mass dampers (TMDs) and active tuned mass dampers (ATMDs) commanded by LQR and fuzzy logic controllers, while the effects of the uncontrolled torsional response of the structure are simultaneously considered. Besides their simplicity, fuzzy logic controllers showed similar trend as LQR controllers under multidirectional wind loads. Nevertheless, the procedure presented in this study can help decision makers, involved in the design process, to choose among innovative solutions like structural control, different damping techniques, modifying geometry, or even changing materials.


---

## Body

## 1. Introduction

Civil engineering structures are an integral part of our modern society. Traditionally, these structures are designed to resist static loads. However, they may be subjected to dynamic loads like earthquakes, winds, waves, and traffic. Such loads can cause severe and/or sustained vibratory motion, which can be detrimental to the structure and human occupants. Because of this, safer and more efficient designs are sought out to balance safety issues with the reality of limited resources. Wind-induced vibrations in buildings are of increasing importance, as the use of high-strength, lightweight materials, longer floor spans, and more flexible framing systems results in structures that are more prone to vibrations. In tall buildings, wind-induced vibrations may cause annoyance to the occupants (especially in the upper floors), impaired function of instruments, or structural damage.Traditionally, wind-induced response of tall buildings in the along-wind direction are evaluated using some codes and formulas [1–4]. However, these standards provide little guidance for the critical cross-wind and torsional responses. This is due to the fact that the cross-wind and torsional responses result mainly from the aerodynamic pressure fluctuations in the separated shear layers and the wake flow fields, which made it difficult to have an acceptable direct analytical relation to the oncoming velocity fluctuations [5, 6]. In addition, these methods may have some limitations, especially for accounting of surrounding tall buildings. Moreover, responses are restricted to some few modes, and the process of evaluating such response depends on much assumption. On the other hand, wind tunnel pressure measurements and finite element modeling (FEM) of the structures are an effective alternative for determining these responses. Wind tunnel tests have been industry wide accepted reliable tools for estimating wind loads on tall buildings. For tall buildings, there are two types of testing: (1) high-frequency base balance (HFBB) and (2) high-frequency pressure integration (HFPI). Inherent in the HFFB approach is the fact that only the global wind loads at the base of the test model are known. The test results from the HFBB measurements can be analyzed using frequency-domain or time-domain techniques to get the building responses. The frequency domain approach has been dominant over time domain approach for its lesser requirement of computational power though it involves more approximations compared to the time-domain approach. Nevertheless, with the current technology where computational power is no longer a problem, the time-domain method is becoming a popular analysis technique. The time-domain method allows determination of wind responses directly from the equation of motion using the measured time history, thereby avoiding all the simplifying assumptions used in the frequency domain technique [7]. However, even if the more accurate time-domain approach is used for the analysis of the response, the three-dimensional (3D) mode shapes found in complex tall buildings complicate the use of the HFBB test results for predicting the response [8, 9]. In general, mode shape correction factors for the HFBB technique are necessary for the assessment of wind-induced responses of a tall building. This is to account for the significant uncertainties in the prediction of generalized forces due to the nonideal mode shapes as well as presumed wind loading distributions [10, 11]. HFPI with the time-domain approach can be more accurate provided that enough coverage of pressure taps on the model’s outer surface is performed.HFPI technique is based on simultaneous pressure measurements at several locations on a building’s outer surface. Pressure data can be used for the design of the claddings as well as the estimation of the overall design loads. The HFPI technique cancels out any inertial effects that may be included in the overall loads measured by the base balance if the HFBB technique was used. Time histories of wind forces at several levels of tall building models can be obtained in the wind tunnel with a multichannel pressure scanning system. This enables the building responses to be computed directly in the time domain for buildings with simple or complex mode shapes. Finite element models (FEMs) can be used for describing the dynamic behavior of the structures. HFPI with FEM have the advantages of considering complex shapes of structures with nonuniform mass distribution and can easily account for any required number of mode shapes to be considered in the response analysis.Preliminary analysis of tall buildings in their preliminary design stages help the designer to make decision by modifying the design or adding passive, active, or semiactive control techniques. Structural control has recently been the subject of much discussion among structural designers. Structural control can potentially provide safer and more efficient structures.  The concept of employing structural control to minimize structural vibrations was proposed in the 1970s [12]. The purpose of structural control is to absorb and to reflect the energy introduced by dynamic loads. The reduction of structural vibrations occurs by adding a mechanical system that is installed in a structure. The control of structural vibrations can also be done by various means, such as modifying rigidities, masses, damping, or shape, and by providing passive or active counter forces [13, 14]. Passive, active and semiactive [15, 16] types of control strategies have been proposed.McNamara [17] studied the tuned mass damper (TMD) as an energy-absorbing system to reduce wind-induced structural response of buildings in the elastic range of behavior. Active control techniques are studied intensively for the control of the response of tall buildings under wind loads [18–22]. The most commonly used active control device for tall buildings is the active tuned mass damper (ATMD). TMDs and ATMDs are shown to be effectiveness in the response reduction of tall buildings under wind loads [23–28].The aim of this study is to present practical procedure for the response prediction and reduction in a very slender high-rise building under multidirectional wind loads. The procedure is schematically presented in Figure1.  Wind loads were obtained from an HFPI experiment conducted in a wind tunnel. The tower responses in the two lateral directions combined with the torsional responses (effect of higher modes on the responses is studied) are evaluated. Two important voids associated with procedures to aid in the design are considered: the first is on the distributions of the wind loads; the second is on the effects of the higher modes. Consideration of these two problems needs wind tunnel pressure measurements on the surface of the building and FEM. The building is modeled using the finite element techniques and a 3D lumped mass model. The uncontrolled responses obtained using the two techniques of modeling are compared. Active control of the structure using LQR and fuzzy logic controllers under wind that is attacking from different directions is proposed. In this study, the lateral responses of the building in the two directions are controlled at the same time, while the effects of the uncontrolled torsional responses of the structure are simultaneously considered.Figure 1
Schematic representation of the proposed procedure for response prediction and reduction in tall buildings under wind loads.


## 2. 3D Modeling of the Building

A 48-strory steel tower proposed in Aly et al. [23] is used in this research. The FEM of the tower, along with the coordinate system, is shown in Figure 2. The full-scale building has a height of about 209 m and a rectangular cross section of B/D≈3 (B: chord length and D: thickness). The aspect ratio in the y-direction is about 11, which makes it very sensitive to strong winds. Modal parameters of the FEM for the first six modes are given in Table 1.Table 1
Modal parameters of the FEM model.
Mode number*Generalized mass × 107 (kg·m2)Generalized stiffness × 109 (N·m)Frequency (Hz)Modal damping11.29530.01470.16940.010220.99370.01780.21320.011230.49450.02220.33700.015040.87240.11150.56890.023450.82730.21530.81200.032660.35440.16001.06950.0426*Modes 1 and 4 are lateral displacements inx-direction; modes 2 and 5 are lateral disp. in y-direction while modes 3 and 6 are torsion.Figure 2
Finite element model with the coordinate system and wind direction (B=57.6 m and D=19 m).


### 2.1. Equations of Motion in Modal Form

Equations of motion governing the behavior of the structure under wind loads are(1)MẌ+CẊ+KX=F(t),
where X=[xy]T is a 2n×1 vector and n is the number of nodes, while x and y are vectors of nodal displacements in x and y directions, respectively. F(t)=[Fx(t)Fy(t)]T, in which Fx(t) and Fy(t) are n×1 vectors of external forces acting in x and y directions, respectively. Using the first N modes obtained by FEM with the next transformation(2)X=ΦQ,
where Φ is 2n×N matrix of eigenvectors and Q is N×1 vector of generalized displacements; that is,(3)Φ=[ϕ1(x1)ϕ2(x1)⋯ϕN(x1)ϕ1(x2)ϕ2(x2)⋯ϕN(x2)⋮⋮⋮ϕ1(xn)ϕ2(xn)⋯ϕN(xn)ϕ1(y1)ϕ2(y1)⋯ϕN(y1)ϕ1(y2)ϕ2(y2)⋯ϕN(y2)⋮⋮⋮ϕ1(yn)ϕ2(yn)⋯ϕN(yn)],Q={q1q2⋮⋮qN}.Substituting by (2) into (1) and premultiplying by ΦT, one obtains
(4)ΦTMΦQ̈+ΦTCΦQ̇+ΦTKΦQ=ΦTF(t).By assuming the damping matrix,C, to be proportional damping, (4) results into six uncoupled equations(5)m11q̈1+c11q̇1+k11q1=∑i=1nΦ1(xi)Fx,i(t)+∑i=1nΦ1(yi)Fy,i(t)=GF1,m22q̈2+c22q̇2+k22q2=∑i=1nΦ2(xi)Fx,i(t)+∑i=1nΦ2(yi)Fy,i(t)=GF2,⋮mNNq̈N+cNNq̇N+kNNqN=∑i=1nΦN(xi)Fx,i(t)+∑i=1nΦN(yi)Fy,i(t)=GFN,
where mii, cii, kii, and GFi are generalized mass, generalized damping, generalized stiffness, and generalized force of the ith mode, respectively. Using the measurements obtained by the pressure transducers, pressure coefficients (matrix Cp) are evaluated at each tap location as a function of time. These values are used with the full-scale model to give the pressure distribution on the surface. The pressure values on the surface of the prototype can be calculated as (6)P(space,time)=12ρU2Cp(space,time),
where P(space,time) is a matrix containing the pressure values on the surface of the full scale model as a function of space (x, y, and z) and time; ρ is the air density which is assumed to be 1.25 kg/m3 (according to [2]), and U is the prototype mean wind speed. The wind load at any node of the outer surface is the integration of the pressure over the surface area in the vicinity of the node as(7)F(nodes,time)=∫P(space,time)dA.This means that once the time history of the pressures on the outer surfaces is calculated, the external forces acting on the nodes of the surface can be computed. The excitation forces acting on the internal nodes are of course equal to zero. Theqi(t) are then solved from each of (5). SIMULINK is used for the numerical solution of these equations [29].


### 2.2. 3D Lumped Mass Model

For control purposes, a 3D lumped mass model is derived from the original FEM. In this model, the total mass of the building was assumed to be lumped at the positions of the floors, and it was assumed for the floors to perform a general 3D movement (each floor has two translations in thex and y directions in addition to the torsional rotation). The building alone (without the control devices) is modeled dynamically using a total of 144 degree-of-freedom. In general, the equations of motion for an n-story building moving in both the two transverse directions and in torsion are written as(8)Msẍ+Csẋ+Ksx=-F+Λf,
where x=[XYΘ]T. The terms X=[x1x2⋯xn] and Y=[y1y2⋯yn] are row vectors of the displacements of the centre of mass of each floor in the x and y directions, respectively, and Θ=[θ1θ2⋯θn] is the vector of the rotations of each floor about the vertical axis (z-axis) while n is the number of stories. The mass matrix, Ms, and the stiffness matrix, Ks, have the following form:(9)Ms=[M000M000I],Ks=[Kx000Ky000Kθ],
where M=diag([m1m2⋯mn]) is a diagonal matrix of lumped masses, I=diag([I1I2⋯In]) in which Ii is the moment of inertia of the ith floor, Kx, Ky, and Kθ are the stiffness matrices in the transverse directions (x and y) and the torsional direction, respectively. The stiffness matrix of the spatial model (3D lumped mass model) is obtained by assuming the stiffness between floors as a combination of cantilever and shear rigidities. MATLAB codes (MATLAB 2008; [29]) were written and used to derive the best stiffness matrix that gives the closest mode shapes to those of the FEM and the same first six natural frequencies. In (8), the disturbance F=[FxFyT]T is a vector of excitation in which Fx  and Fy are two vectors of the horizontal loads acting in the x and y directions, respectively, and T is a vector of the external  torsional wind loads. Also, f is the vector of control forces, where its coefficient matrix Λ is the matrix determined by the location of control devices.Wind loading vectors (Fx, Fy, and T) lumped at the position of floors are obtained from wind tunnel pressure tests conducted at the wind tunnel of Politecnico di Milano [30] on a scaled 1 : 100 rigid model of the tower. Such large-scale allows for the advantage of testing the model at high Reynolds number with minimum blockage due to the huge dimensions of the test section. Pressure taps were distributed on the outer surface of the test model. To allow for sufficient pressure measurements (see [7]), 400 taps were mapped on the outer surface of the model. Pressure taps were distributed to cover the entire outer surface with more intense at the upper part of the test model (see Figure 5(a)).  Pressure data were acquired at a frequency of 62.5 Hz using 448 pressure taps. Figure 6 shows a typical spectrum of measured pressure data at model scale. Pressure data were integrated on the outer surface of the building (see Figure 5(b)) to obtain the corresponding time history of the two directional wind loads at each floor in addition to torsion. For the estimation of the wind loads at each floor, the tributary area for each floor was gridded into smaller areas and the time history of the wind loads at each area was found by using the Cp records of the closest pressure tap (see Figure 7). Codes were written in MATLAB to estimate the time histories of the wind forces acting at the center of each smaller area. After that, the floor forces in the two directions were obtained from the summation of the forces in each lateral direction. The torsion at each floor was the resultant of the summation of the force moments about the floor center. The surrounding buildings within a radius of 500 m from the centre of the tower were also scaled 1 : 100 to be presented on the turning test table according to the type of the test configuration used. The wind profile represents a typical urban terrain as shown in Figure 3 [31].  The reference mean wind speed (Uref=15 m/s) was measured at a height of 1 m. Prototype reference mean wind speed is assumed to be 30 m/s. The target for the wind profiles is the Eurocode 1 [2]. The turbulence intensities in the longitudinal, lateral, and vertical directions are referred to by Iu, Iv, and Iw, respectively. Two different test configurations with the same wind profile are considered in this paper (see Figure 4). In the first configuration, the building is subjected to the wind load without the existence of other tall buildings, and is referred to as Config. no. 1. In the second configuration, the rigid model of the building was tested with the existence of all the surrounding buildings. This configuration is referred to as Config. no. 2. Further details about the wind tunnel experiment are given in Aly [24].Mean wind speed profile, turbulence intensity profiles, and wind spectra (L is the integral scale).
(a)(b)(c)Two different configurations are used.
(a)(b)Pressures on the outer surface of a scaled 1 : 100 model were obtained from a wind tunnel test: (a) pressure tap distribution (elevation and side view), (b) mean surface pressure coefficient distribution (for 292.5 deg).
(a)(b)Figure 6
Typical spectrum of measured pressure data (at model scale).Wind load estimation from pressure data: the tributary area of floorN was divided into smaller areas; pressure forces acting on each smaller area, Ai,j, were calculated based on pressure data at the nearest pressure tap, m.
(a)(b)Due to the fact that the building’s mass is symmetrical, and the study is based on the assumption that the structure is responding in the linear region, the lateral and the torsional behavior of the building may be studied alone, then the response time histories may be combined simultaneously. In this study, the plane motion of the structure in thex-direction is controlled using both the TMDs and ATMDs. However, due to the fact that the control of the response in the x-direction will not affect the response of the building in the y-direction, another TMD and ATMD are designed to control the lateral in-plane response in the y-direction. Following that, the uncontrolled torsional response is added simultaneously to the two lateral responses to give the overall response in the two lateral directions.The state reduction approach derived by Davison [32] and summarized later in Wu et al. [33] is used in this study (also see Lu et al. [20]). In this approach, the 48 degree-of-freedom (DOF) in-plane system is reduced to 15 DOF, where the first 30 modes are retained [23]. Note that the condition for this approach was that the response in terms of displacements and accelerations of the 15 DOF and 48 DOF are very much the same (see Section 4). This model is referred to as the reduced order system (ROS). The addition of the TMD increases the DOF to 16. The system with the TMD is referred to as ROS-TMD. In a similar task, ROS-ATMD refers to the ROS utilizing ATMD.The state equation of the ROS that corresponds to the full order system (FOS) in (8) can be expressed as(10)ż=Az+Bf+Ew,
in which z=[X̅′;X̅̇′] is the 32-dimensional state vector, X̅ is a vector of the in-plane displacements of floors 3, 6, 9, 12, 16, 19, 22, 26, 29, 32, 35, 38, 41, 44, and 48 in addition to the displacement of the inertial mass of the damper.  A is a (32×32) system matrix, B is a 32 location vector, and E is a 32 excitation vector. In the system reduction above, the wind loads acting on each of the 15 floors above are computed from the wind loads F acting on each of the 48 floors by lumping wind forces on adjacent floors at the locations that correspond to the 15 DOF model.The controlled output vector,yc, and the measured output, ym, of the ROS described by (10) can be expressed as(11)yc=Ccz+Dcf+Fcwx,ym=Cmz+Dmf+Fmwx+ν,
where Cc, Dc, Fc, Cm, Dm, and Fm are matrices with appropriate dimensions and ν is the measurement noise vector. The model used for controller design was further reduced as follows:(12)żr=Arzr+Brf+Erwx,ycr=Ccrzr+Dcrf+Fcrwx,ymr=Cmrzr+Dmrf+Fmrwx+νr,
where zr is a 6-dimensional state vector of the reduced order system, ycr is a controlled output vector identical of yc defined by (11), ymr is the measured output vector; νr is the measurement noise, and Ccr, Dcr, Fcr, Cmr, Dmr, and Fmr  are appropriate matrices.


## 2.1. Equations of Motion in Modal Form

Equations of motion governing the behavior of the structure under wind loads are(1)MẌ+CẊ+KX=F(t),
where X=[xy]T is a 2n×1 vector and n is the number of nodes, while x and y are vectors of nodal displacements in x and y directions, respectively. F(t)=[Fx(t)Fy(t)]T, in which Fx(t) and Fy(t) are n×1 vectors of external forces acting in x and y directions, respectively. Using the first N modes obtained by FEM with the next transformation(2)X=ΦQ,
where Φ is 2n×N matrix of eigenvectors and Q is N×1 vector of generalized displacements; that is,(3)Φ=[ϕ1(x1)ϕ2(x1)⋯ϕN(x1)ϕ1(x2)ϕ2(x2)⋯ϕN(x2)⋮⋮⋮ϕ1(xn)ϕ2(xn)⋯ϕN(xn)ϕ1(y1)ϕ2(y1)⋯ϕN(y1)ϕ1(y2)ϕ2(y2)⋯ϕN(y2)⋮⋮⋮ϕ1(yn)ϕ2(yn)⋯ϕN(yn)],Q={q1q2⋮⋮qN}.Substituting by (2) into (1) and premultiplying by ΦT, one obtains
(4)ΦTMΦQ̈+ΦTCΦQ̇+ΦTKΦQ=ΦTF(t).By assuming the damping matrix,C, to be proportional damping, (4) results into six uncoupled equations(5)m11q̈1+c11q̇1+k11q1=∑i=1nΦ1(xi)Fx,i(t)+∑i=1nΦ1(yi)Fy,i(t)=GF1,m22q̈2+c22q̇2+k22q2=∑i=1nΦ2(xi)Fx,i(t)+∑i=1nΦ2(yi)Fy,i(t)=GF2,⋮mNNq̈N+cNNq̇N+kNNqN=∑i=1nΦN(xi)Fx,i(t)+∑i=1nΦN(yi)Fy,i(t)=GFN,
where mii, cii, kii, and GFi are generalized mass, generalized damping, generalized stiffness, and generalized force of the ith mode, respectively. Using the measurements obtained by the pressure transducers, pressure coefficients (matrix Cp) are evaluated at each tap location as a function of time. These values are used with the full-scale model to give the pressure distribution on the surface. The pressure values on the surface of the prototype can be calculated as (6)P(space,time)=12ρU2Cp(space,time),
where P(space,time) is a matrix containing the pressure values on the surface of the full scale model as a function of space (x, y, and z) and time; ρ is the air density which is assumed to be 1.25 kg/m3 (according to [2]), and U is the prototype mean wind speed. The wind load at any node of the outer surface is the integration of the pressure over the surface area in the vicinity of the node as(7)F(nodes,time)=∫P(space,time)dA.This means that once the time history of the pressures on the outer surfaces is calculated, the external forces acting on the nodes of the surface can be computed. The excitation forces acting on the internal nodes are of course equal to zero. Theqi(t) are then solved from each of (5). SIMULINK is used for the numerical solution of these equations [29].


## 2.2. 3D Lumped Mass Model

For control purposes, a 3D lumped mass model is derived from the original FEM. In this model, the total mass of the building was assumed to be lumped at the positions of the floors, and it was assumed for the floors to perform a general 3D movement (each floor has two translations in thex and y directions in addition to the torsional rotation). The building alone (without the control devices) is modeled dynamically using a total of 144 degree-of-freedom. In general, the equations of motion for an n-story building moving in both the two transverse directions and in torsion are written as(8)Msẍ+Csẋ+Ksx=-F+Λf,
where x=[XYΘ]T. The terms X=[x1x2⋯xn] and Y=[y1y2⋯yn] are row vectors of the displacements of the centre of mass of each floor in the x and y directions, respectively, and Θ=[θ1θ2⋯θn] is the vector of the rotations of each floor about the vertical axis (z-axis) while n is the number of stories. The mass matrix, Ms, and the stiffness matrix, Ks, have the following form:(9)Ms=[M000M000I],Ks=[Kx000Ky000Kθ],
where M=diag([m1m2⋯mn]) is a diagonal matrix of lumped masses, I=diag([I1I2⋯In]) in which Ii is the moment of inertia of the ith floor, Kx, Ky, and Kθ are the stiffness matrices in the transverse directions (x and y) and the torsional direction, respectively. The stiffness matrix of the spatial model (3D lumped mass model) is obtained by assuming the stiffness between floors as a combination of cantilever and shear rigidities. MATLAB codes (MATLAB 2008; [29]) were written and used to derive the best stiffness matrix that gives the closest mode shapes to those of the FEM and the same first six natural frequencies. In (8), the disturbance F=[FxFyT]T is a vector of excitation in which Fx  and Fy are two vectors of the horizontal loads acting in the x and y directions, respectively, and T is a vector of the external  torsional wind loads. Also, f is the vector of control forces, where its coefficient matrix Λ is the matrix determined by the location of control devices.Wind loading vectors (Fx, Fy, and T) lumped at the position of floors are obtained from wind tunnel pressure tests conducted at the wind tunnel of Politecnico di Milano [30] on a scaled 1 : 100 rigid model of the tower. Such large-scale allows for the advantage of testing the model at high Reynolds number with minimum blockage due to the huge dimensions of the test section. Pressure taps were distributed on the outer surface of the test model. To allow for sufficient pressure measurements (see [7]), 400 taps were mapped on the outer surface of the model. Pressure taps were distributed to cover the entire outer surface with more intense at the upper part of the test model (see Figure 5(a)).  Pressure data were acquired at a frequency of 62.5 Hz using 448 pressure taps. Figure 6 shows a typical spectrum of measured pressure data at model scale. Pressure data were integrated on the outer surface of the building (see Figure 5(b)) to obtain the corresponding time history of the two directional wind loads at each floor in addition to torsion. For the estimation of the wind loads at each floor, the tributary area for each floor was gridded into smaller areas and the time history of the wind loads at each area was found by using the Cp records of the closest pressure tap (see Figure 7). Codes were written in MATLAB to estimate the time histories of the wind forces acting at the center of each smaller area. After that, the floor forces in the two directions were obtained from the summation of the forces in each lateral direction. The torsion at each floor was the resultant of the summation of the force moments about the floor center. The surrounding buildings within a radius of 500 m from the centre of the tower were also scaled 1 : 100 to be presented on the turning test table according to the type of the test configuration used. The wind profile represents a typical urban terrain as shown in Figure 3 [31].  The reference mean wind speed (Uref=15 m/s) was measured at a height of 1 m. Prototype reference mean wind speed is assumed to be 30 m/s. The target for the wind profiles is the Eurocode 1 [2]. The turbulence intensities in the longitudinal, lateral, and vertical directions are referred to by Iu, Iv, and Iw, respectively. Two different test configurations with the same wind profile are considered in this paper (see Figure 4). In the first configuration, the building is subjected to the wind load without the existence of other tall buildings, and is referred to as Config. no. 1. In the second configuration, the rigid model of the building was tested with the existence of all the surrounding buildings. This configuration is referred to as Config. no. 2. Further details about the wind tunnel experiment are given in Aly [24].Mean wind speed profile, turbulence intensity profiles, and wind spectra (L is the integral scale).
(a)(b)(c)Two different configurations are used.
(a)(b)Pressures on the outer surface of a scaled 1 : 100 model were obtained from a wind tunnel test: (a) pressure tap distribution (elevation and side view), (b) mean surface pressure coefficient distribution (for 292.5 deg).
(a)(b)Figure 6
Typical spectrum of measured pressure data (at model scale).Wind load estimation from pressure data: the tributary area of floorN was divided into smaller areas; pressure forces acting on each smaller area, Ai,j, were calculated based on pressure data at the nearest pressure tap, m.
(a)(b)Due to the fact that the building’s mass is symmetrical, and the study is based on the assumption that the structure is responding in the linear region, the lateral and the torsional behavior of the building may be studied alone, then the response time histories may be combined simultaneously. In this study, the plane motion of the structure in thex-direction is controlled using both the TMDs and ATMDs. However, due to the fact that the control of the response in the x-direction will not affect the response of the building in the y-direction, another TMD and ATMD are designed to control the lateral in-plane response in the y-direction. Following that, the uncontrolled torsional response is added simultaneously to the two lateral responses to give the overall response in the two lateral directions.The state reduction approach derived by Davison [32] and summarized later in Wu et al. [33] is used in this study (also see Lu et al. [20]). In this approach, the 48 degree-of-freedom (DOF) in-plane system is reduced to 15 DOF, where the first 30 modes are retained [23]. Note that the condition for this approach was that the response in terms of displacements and accelerations of the 15 DOF and 48 DOF are very much the same (see Section 4). This model is referred to as the reduced order system (ROS). The addition of the TMD increases the DOF to 16. The system with the TMD is referred to as ROS-TMD. In a similar task, ROS-ATMD refers to the ROS utilizing ATMD.The state equation of the ROS that corresponds to the full order system (FOS) in (8) can be expressed as(10)ż=Az+Bf+Ew,
in which z=[X̅′;X̅̇′] is the 32-dimensional state vector, X̅ is a vector of the in-plane displacements of floors 3, 6, 9, 12, 16, 19, 22, 26, 29, 32, 35, 38, 41, 44, and 48 in addition to the displacement of the inertial mass of the damper.  A is a (32×32) system matrix, B is a 32 location vector, and E is a 32 excitation vector. In the system reduction above, the wind loads acting on each of the 15 floors above are computed from the wind loads F acting on each of the 48 floors by lumping wind forces on adjacent floors at the locations that correspond to the 15 DOF model.The controlled output vector,yc, and the measured output, ym, of the ROS described by (10) can be expressed as(11)yc=Ccz+Dcf+Fcwx,ym=Cmz+Dmf+Fmwx+ν,
where Cc, Dc, Fc, Cm, Dm, and Fm are matrices with appropriate dimensions and ν is the measurement noise vector. The model used for controller design was further reduced as follows:(12)żr=Arzr+Brf+Erwx,ycr=Ccrzr+Dcrf+Fcrwx,ymr=Cmrzr+Dmrf+Fmrwx+νr,
where zr is a 6-dimensional state vector of the reduced order system, ycr is a controlled output vector identical of yc defined by (11), ymr is the measured output vector; νr is the measurement noise, and Ccr, Dcr, Fcr, Cmr, Dmr, and Fmr  are appropriate matrices.


## 3. Controllers and Limitations

In this study, both TMDs and ATMDs are used for the reduction of the lateral responses of the building. However, in order to make the design of such control systems more realistic and applicable, the following restrictions and assumption are applied.(i)
The mass of the TMD in thex-direction is 100 ton, while the mass of the TMD in the y-direction is 150 ton. Such restrictions are applied to avoid excessive weight on the roof (the overall mass on the roof is about 0.625% of the overall building’s mass).(ii)
The TMDs are tuned to the first vibrational mode in each corresponding lateral direction. The damping factor is taken to be 20% of the critical. This amount of damping is selected higher than the optimal value for the sake of restricting the stroke of the ATMDs.(iii)
The maximum stroke of the actuators is restricted to 1.5 m.(iv)
The maximum control force of the actuator in they-direction is restricted to 100 kN, and that in the x-direction is restricted to 25 kN.(v)
The computational delay and the sampling rate of the digital controller are 0.001 s.(vi)
Three acceleration measurements are available for each lateral direction (at floor 30, roof, and mass of the TMD).Note that the tower required a TMD with heavier mass and ATMD with higher control force in one lateral direction than the other, which was basically attributed to geometry. A linear-quadratic regulator (LQR) design with output weighting is selected to give the desired control force using the MATLAB function (lqry.m). The state-feedback law f=-Gzr minimizes the cost function (13)J(f)=∫0∞(ycr′Qycr+f′Rf)dt,
where G is the feedback gain matrix, zr is a 6-dimensional state vector of the reduced order system, ycr is the measured output vector, the symbol (′) denotes transpose, and Q and R are weighting matrices. Parametric studies were performed with various weighting matrices Q, corresponding to various regulated output vectors. The results of these parametric studies indicated that an effective controller could be designed by selecting a vector of regulated responses to include the velocities of each floor.For comparison reasons, fuzzy logic controllers are used in this study to command the actuators of the ATMDs (see Nguyen et al. [34]). From a design point of view, fuzzy logic controllers do not require the complexity of a traditional control system. The measured accelerations can be used directly as input to the fuzzy controller. The main advantages of using a fuzzy control algorithm are summarized in Battaini et al. [35] and Samali et al. [36]. According to Samali et al. [36], uncertainties of input data are treated in a much easier way by fuzzy control theory than by classical control theory. Since fuzzy controllers are based on linguistic synthesis, they possess inherent robustness. Fuzzy controllers can be easily implemented in a fuzzy chip with immediate reaction time and autonomous power supply. Furthermore, the design of fuzzy controller does not require state reduction or concerning about observers. Only two acceleration measurements were used (floor 30 and roof).The input variables to the fuzzy controller were selected as accelerations of floors 30 and 48 and the output as the control force. The membership functions for the inputs were defined and selected as seven triangles with overlaps as shown in Figure8. For the output, they were defined and selected as nine triangles with overlaps as shown in Figure 9. The fuzzy variables used to define the fuzzy space are ZR (zero), PVS (positive very small), PS (positive small), PM (positive medium), PL (positive large), PVL (positive very large), NVS (negative very small), NS (negative small), NM (negative medium), NL (negative large), and NVL (negative very large). The rule-base for computing the desired current is presented in Table 2 [36].Table 2
Control rule base [36].
Acceleration of 48th floorAcceleration of 30th floorNLNMNSZRPSPMPLNLPVLPVLPLPVSZRZRZRNMPLPLPMPVSZRZRZRNSZRNVSPMPSPVSZRZRZRZRZRNVSZRPVSZRZRPSZRZRNVSNSNMPVSZRPMZRZRZRNVSNMNLNLPLZRZRZRNVSNLNVLNVLMembership functions for the input measured accelerations in thex-direction (Acc-x-30, Acc-x-48) and the y-direction (Acc-y-30, Acc-y-48).
(a)(b)(c)(d)Membership functions for the output control force in thex-direction (Force-x) and the y-direction (Force-y).
(a)(b)


## 4. Results and Discussion

Table3  gives the response of the top corner of the building in the y-direction for an incident angle of 0° under different consideration of mode shapes. It is shown that the displacement response of this building is dominated by the first lateral mode in the y-direction (modes 1 : 2 in the table). Nevertheless, this underestimates the displacement response by 3% to 4.4% and the acceleration response by about 12% to 17%. Note that the aspect ratio of this building in the y-direction is about 11. This means that for very slender buildings, the solo consideration of the first lateral mode may lead to significant error in the estimation of the response, especially for the acceleration response. Table 4 lists the response of the top corner of the tower in the x-direction for an incident angle of 90° under different consideration of mode shapes. It is shown that the displacement and acceleration response are dominated by the first lateral mode in the x-direction (modes 1 in the table). Note that the aspect ratio of this building in the x-direction is about 3.6. This means that for buildings with low aspect ratio, the solo consideration of the first lateral mode may be sufficient for the estimation of the response. Figure 10 shows the power spectra of the acceleration response of the top corner of the building in the two lateral directions. The figure shows that the third mode (torsion) contributes significantly to the acceleration in the y-direction. In general, results given by the figure, Tables 3 and 4 show that the responses of tall buildings under winds are dominated by the first few modes (for this specific building, the first two lateral modes and the first torsional mode can be sufficient).Table 3
Response of the top corner of the tower in they-direction for an incident angle of 0°.
ModeRMS-disp. (m)Max-disp. (m)RMS-accel. (m/s2)Max-accel. (m/s2)10.000 (−100%)0.001 (−99.8%)0.000 (100%)0.001 (−99.9%)1 : 20.129 (−4.4%)0.587 (−2.8%)0.199 (−17.1%)0.855 (−11.8%)1 : 30.136 (0.7%)0.613 (1.5%)0.238 (−0.8%)0.980 (1.1%)1 : 40.136 (0.7%)0.613 (1.5%)0.238 (−0.8%)0.980 (1.1%)1 : 50.135 (0%)0.606 (0.3%)0.239 (−0.4%)0.966 (−0.3%)1 : 60.135 (0%)0.604 (0%)0.240 (0%)0.969 (0%)Table 4
Response of the top corner of the tower in thex-direction for an incident angle of 90°.
ModeRMS-disp. (m)Max-disp. (m)RMS-accel. (m/s2)Max-accel. (m/s2)10.188 (1.1%)0.646 (−0.5%)0.203 (−0.5%)0.654 (−3.5%)1 :  20.188 (1.1%)0.646 (−0.5%)0.203 (−0.5%)0.654 (−3.5%)1 : 30.187 (0.5%)0.648 (−0.2%)0.204 (0%)0.653 (−3.7%)1 : 40.186 (0%)0.649 (0%)0.204 (0%)0.676 (−0.3%)1 : 50.186 (0%)0.649 (0%)0.204 (0%)0.676 (−0.3%)1 : 60.186 (0%)0.649 (0%)0.204 (0%)0.678 (0%)Figure 10
Power spectra of the acceleration response of the top corner of the building in the two lateral directions.Figure11 gives displacement and acceleration responses of a point at the top corner of the building for the FEM, the 3D full order system (3D-FOS), and the 3D reduced order system (3D-ROS). The figure shows that the response in terms of displacements and accelerations of the three types of modeling are very much the same (see Section 4). This means that FE modeling, 3D lumped mass modeling, and 3D reduced order modeling of tall buildings under wind loads can give an accurate assessment of the response provided that the first dominant modes are retained. The figure shows also that the cross-wind response is higher than the along-wind response. This reveals the importance of the procedure proposed in this study as many design codes and formula may provide accurate estimate of the along-wind response but less guidance for the estimation of the critical cross-wind and torsional response. The results show that the building is very much vulnerable to wind loads. This may be due to its low weight and its low dominant frequencies.displacement and acceleration responses of a point at the top corner for FEM, 3D full order system (3D-FOS), and 3D reduced order system (3D-ROS).
(a)(b)(c)(d)The building required a TMD with heavier mass and ATMD with higher control force in one lateral direction than the other. This may be attributed to geometry. Figures11–14 show the controlled and uncontrolled responses of the tower under wind loads for two test configurations at different incident angles. Two examples of control are considered, TMDs and ATMDs with LQR and fuzzy logic controllers. For each example, the controlled responses in the x and y directions are plotted with the uncontrolled responses. The controlled and uncontrolled responses of the tower are evaluated by simulations  [29]. Four evaluation criteria are used to examine the performance of the proposed controllers. Evaluation criteria include: rms displacements, maximum displacements, rms accelerations, and maximum accelerations of the top corner of the tower in the two lateral directions. The figures are superimposed by ellipses indicating the position of the most unfavourable responses (uncontrolled, with TMDs, with ATMDs (LQR), and with ATMDs (fuzzy)) over the two configurations in both x and y directions. The amount of reduction in the highest response achieved by TMDs and ATMDs over the worst uncontrolled response is indicated in the figures.RMS-displacements of the top corner of the tower.
(a)(b)(c)(d)Maximum displacements of the top corner of the tower.
(a)(b)(c)(d)RMS accelerations of the top corner of the tower.
(a)(b)(c)(d)Figures11 and 12 give controlled and uncontrolled rms displacements and max displacements of the top corner of the tower in both the x and y directions. It is shown that TMDs have a great effect on the reduction of the displacement response of the building. Reductions achieved by TMDs in the displacements responses range from 22–30% over the worst uncontrolled response. Generally, TMDs are able to give good reduction in the rms displacements in both the x and y directions for all wind incident angles. Reductions achieved by ATMDs in the displacement responses range from 29%–43% over the worst uncontrolled response. ATMDs with fuzzy logic controllers are able to enhance the reduction in the displacement responses over LQR (perhaps Q and R could be more optimal) most of the time (by about 1% to 5%). They also have a general similar trend over all of the wind attack angles.Figures13 and 14 give controlled and uncontrolled rms accelerations and maximum accelerations (Figure 15) of the top corner of the tower in both the x and y directions. It is shown that the TMDs have a significant effect on the reduction of the acceleration response of the building. Reductions achieved by TMDs in the acceleration responses range from 16%–30% over the worst uncontrolled response.Maximum accelerations of the top corner of the tower.
(a)(b)(c)(d)Generally, TMDs are able to give good reduction in the rms displacements in both thex and y directions for all the wind incident angles. However, the performance is limited in reducing the along-wind maximum acceleration of the tower in the y-direction under Config. no. 2, when the wind direction angle is 90°. This may be due to the interference effects of two high-rise buildings in the oncoming wind (see Figure 4).  Results also show that ATMDs are able to enhance the reduction in the responses. Reductions achieved by ATMDs in the acceleration responses range from 21%–43% over the worst uncontrolled response. ATMDs with fuzzy logic controllers are able to enhance the reduction in the acceleration responses over LQR, and in general, they have a similar trend over all of the wind incident angles.As a general comment on Figures11–14, one can see that the performance of the controllers is much better in the x-direction. In addition, the capability of the controllers to reduce the responses (especially maximum accelerations at angles 0° and 180°) in the y-direction is limited. This may be due to the effect of vortex shedding on the across-wind responses. Moreover, the structure is slender in x-direction (see Figure 2). The structure is also stiffer in the y-direction (see Table 1 for natural frequencies). However, the procedure presented in this study permits the response of tall buildings to be assessed and controlled in the preliminary design stages which can help decision makers, involved in the design process, to choose among innovative design solutions like structural control, considering several damping techniques, modifying geometry, or even changing materials (e.g., from steel to concrete).


## 5. Conclusions

The paper presents practical procedure for the response prediction and reduction in high-rise buildings under wind loads. To show the applicability of the procedure, aerodynamic loads acting on a quasirectangular high-rise building based on an experimental approach (surface pressure measurement) are used with a mathematical model of the structure for the response prediction and reduction. The building represents a case study of an engineered design of a very slender tower that is instructive. The contributions of this paper can be summarized as follows.(1)
The methodology based on HFPI and FEM proposed for the estimation of the response of high-rise buildings under wind loads has the advantage of combining lateral along-wind, lateral cross-wind, and torsional responses altogether. The technique allows for the consideration of any number of modes.(2)
Results show that the responses of tall buildings under winds are dominated by the first few modes. Consequently, FEM, 3D lumped mass modeling, and reduced order 3D modeling of tall buildings under wind loads give an accurate assessment of the response provided that the first dominant modes are retained.(3)
Results show that the response of tall buildings in the cross-wind direction (lateral response combined simultaneously with torsion) can be higher than the response in the along-wind direction. This reveals the importance of the procedure proposed in this study as many design codes and formula may provide accurate estimate of the along-wind response but less guidance for the estimation of the critical cross-wind and torsional response.(4)
The building represents an engineered steel design of a structure that is very much vulnerable to wind loads. This may be due to its low weight as well as high flexibility related to the low dominant frequencies and the high aspect ratio.(5)
The building demands TMD with heavier mass and ATMD with higher control force in one lateral direction than the other. This may be attributed to geometry.(6)
For the purpose of the use of active control, LQR and fuzzy logic controllers are shown to be effective in enhancing the response reduction over the TMD. ATMDs with fuzzy logic controllers show similar trend like LQR controllers under multidirectional wind loads. In addition, from a design point of view, fuzzy logic controllers do not require the complexity of traditional control systems.(7)
The procedure presented in this paper permits the response of tall buildings to be assessed and controlled in the preliminary design stages. This can help decision makers, involved in the design process, to choose among innovative design solutions like structural control, considering several damping techniques, modifying geometry, or even changing materials.



---

*Source: 549621-2011-07-14.xml*