LQG Control of Along-Wind Response of a Tall Building with an ATMD

1 The Department of Architectural Engineering, Chonbuk National University, 567 Baekje-daero, deokjin-gu, Jeonju-si, Jeollabuk-do 561-756, Republic of Korea 2 KOCEDWind Tunnel Center, Chonbuk National University, 567 Baekje-daero, deokjin-gu, Jeonju-si, Jeollabuk-do 561-756, Republic of Korea 3The Department of Architectural Engineering, Chonbuk National University and Long-Span Steel Frame System Research Center, 567 Baekje-daero, deokjin-gu, Jeonju-si, Jeollabuk-do 561-756, Republic of Korea


Introduction
The majority of today's tall buildings are made from lighter construction materials that are both stronger and more flexible.Although structural safety is greatly improved as a result of using these materials, excessive wind-induced vibration can occur, resulting in occupant discomfort [1].It is necessary to reduce the wind-induced displacement and acceleration at the top floor level of tall buildings in order to increase occupant comfort and structural safety [2][3][4].Although much research has been undertaken to estimate and mitigate the wind-induced structural vibration of tall buildings, the complicated mechanism, interacting with the fluctuating atmospheric flow and building sides, has not been developed successfully [5][6][7].The fluctuating along-wind load, which has almost the same property of approaching turbulent flow velocity, could be predicted theoretically, so the along-wind response of a tall building could be estimated by the use of the gust factor method [8,9].The across-wind vibration response of a tall building is much more sensitive and larger than that of along-wind response.And the across-wind response is mainly due to the vortex-shedding reattachment and wake flow which occur on the side face of leeward side and cannot be formulated theoretically [6].
Studies on mitigating such an excessive wind-induced vibration have been conducted over several decades [2][3][4].In recent years, researchers have utilized the modern optimal control theory and devices to obtain a required control force and reduced vibration response [3,4,10,11].In 1972, Yao introduced modern control theory into the vibration control of civil engineering structures.Due to the fact that modern tall buildings are subjected to fluctuating wind load, oscillations at the fundamental natural frequency of the building occur.Modern control theory and auxiliary devices have been applied to control wind-induced excitations of tall buildings.One of the common control devices is a tuned mass damper (TMD) system.The concept of TMD and its relation to reducing structural responses comes from Den Hartog [12].It consists of an auxiliary mass with properly tuned spring and damping devices, which increase damping in the main structure to reduce the response of the main structure.Many studies have been conducted on the behavior and effectiveness of TMD and a number of TMDs have been installed in tall buildings for the control of wind-induced vibration [3,10].The Center Point Tower in Sydney is one of the first examples of TMD use in a building.When wind load is modeled as a stationary Gaussian white noise, the TMD parameter for reducing wind-induced excitations in a building could be derived by McNamara [2].At that time, 400-ton TMD had been installed for the Citicorp Center, a 274 m tall office building in New York.Another TMD has been designed in the John Hancock Tower, Boston.Both of these TMD examples have been installed to reduce wind-induced vibrations [3,10].It was accepted as common knowledge that the performance of TMD could be increased by incorporating a feedback controller through the use of an active control force in the design of TMD.This is called an ATMD [3,11,13,14].An ATMD design for mitigating the wind-induced vibrations endured by tall buildings with an LQR controller that uses a deterministic harmonic wind load was presented by Chang and Soong [10].Theirs is the first active control study on reducing the wind-induced vibration of tall buildings.Following this, many studies on optimal control algorithms, including LQR, LQG H2, and Hinf, have been developed to obtain the optimal control force for minimizing the wind-induced vibration of a tall building [11,[14][15][16][17][18][19][20][21][22][23].However, it is still believed that ATMD could be superior to TMD in reducing wind-induced vibrations endured by tall buildings [11].In terms of the modern optimal control theory, fluctuating along-wind load acting on a tall building can be considered a stationary Gaussian white noise process that has a constant power spectral density function.Therefore, a linear control system that has a system noise and measurement noise (a Gaussian white noise) is considered to be a stochastic linear control system.An LQG controller is used to investigate the effectiveness of ATMD in reducing the along-wind response of a tall building [14,18,19].Fluctuating along-wind load, acting on a tall building, was simulated numerically in the time domain using the along-wind load spectra proposed by Solari [24].The simulation procedure used in this study was taken from Schueller and Shinozuka [26] and Deodatis [25].Using this simulated fluctuating along-wind load, the estimated along-wind responses of a tall building were undertaken with and without an ATMD, while the LQG controller revealed the effectiveness of the ATMD and the LQG controller in reducing the along-wind response of a tall building.

Equations of Motion
A tall building, which has ATMD installed at the top floor level with an active control force (such as an actuator), is shown in Figure 1.The building is modeled as an equivalent single degree of freedom system with a generalized mass  1 , damping  1 , and stiffness  1 , which corresponds to the first mode modal mass, damping, and stiffness of the building.  ,   , and   are the corresponding quantities of mass, damping, and stiffness constant of TMD, while () represents the along-wind force.() is the active control force of the LQG controller.The linear dynamic equations of the motion of the system can be written as [10,27,28] Here, () =   () −  1 () is the displacement of   relative to that of  1 .This equation can be written in terms of the state-space formulation, where denotes the state vector of the system, is a system dynamic matrix, Mathematical Problems in Engineering 3 is a location vector of (), and is a location vector of () [29].

Optimum Parameters of ATMD
While the basic concept of TMD on reducing the main body of structural response has been well established, the optimal parameters of TMD could be different for different structures and external loading conditions [30].Warburton investigated the effectiveness of TMD on reducing the random responses of the main structure under a random load, which has Gaussian white noise spectra of constant value.Warburton derived the optimum parameters for the tuning frequency ( opt ) and the damping ratio ( opt ) of the passive TMD, under stationary random excitations, which have white noise spectra.The optimum TMD parameters of tuning frequency and damping ratio were determined for a mass ratio  [30] as Here,  = mass ratio of m d /m 1.
The optimal tuning frequency ratio is greater for random excitation than for harmonic excitation [31].However, the optimal damping ratio is smaller for random excitation than for harmonic excitation.The optimum parameters of ATMD, meanwhile, are the same as those of passive TMD.

Numerical Simulation of Fluctuating Along-Wind Load
Fluctuating along-wind load (treated as a random process of stationary Gaussian white noise) can be simulated numerically in the time domain using along-wind load power spectral density data.This is particularly useful for some response estimations, which are more or less narrow banded random processes, such as the along-wind response of a tall building.The numerical simulation procedure presented in this work is taken from Scheuller and Shinozuka [26] and Deodatis [25] as Here,   ( 1 ) = the value of the spectral density of the along-wind load, corresponding to the first modal resonant frequency.The along-wind load power spectral density, used in (6), is by Solari [24].It is as follows:

Δ𝜔 = (𝜔
where where where  = /(),   ,   = lateral and vertical exponential decay coefficients,   = cross-correlation coefficient of pressure acting on the windward and leeward face,  V (ℎ) = integral length scale of turbulence at height ℎ,  = air density, B = width of building, H = height of building, h = reference height of building,   = drag coefficient,   ,   = absolute values of mean pressure coefficients on windward and leeward face,  = mean wind velocity,  V = standard deviation of longitudinal turbulence, n = frequency, and   () = power spectrum of first fluctuating modal force.

Linear Quadratic Gaussian Controller
A tall building that is subjected to a fluctuating along-wind load and has a stationary white noise spectrum can be considered as the linear dynamic system, which has a system and measurement noise.If the system and measurement noises are zero-mean Gaussian white noise that has a constant covariance intensity, they are stationary white noise spectra matrices.If the external random load (a fluctuating alongwind load, for instance) acting on a tall building is considered a system noise that has a constant power spectral density of Gaussian white noise, we can formulate a plant model dynamic with an LQG controller in modern optimal control theory.This is as follows [27]: () =  () +  () + V () .
Here, we let  = 0 for simplicity.() and () are the state and output vectors, while () and V() are the system and measurement noises, respectively.These are assumed to be uncorrelated zero-mean Gaussian white noises that have covariance intensity matrices of  and , respectively.That is,
It is known that the optimal controller () in ( 12) is obtained when all states () of the system and the output () are a combination of all states.However, in practice, all states () are not available and system and output measurements are driven by stochastic disturbance (called noises), which have a constant power spectral density of Gaussian white noise.In these situations, we need a state estimator or observer to estimate all states of the system.The design of the state observer could then be performed using a Kalman filter, an optimal state observer for a stochastic dynamic system [27].
If we take () to be the state estimate and () = () − X() to denote the estimation error, then the state can be estimated using a Kalman filter.This is described as Ẋ () =  X () +  () +  ( () −  X ()) , (15) where the observer gain matrix () is given by where and where Γ satisfies the filter algebraic Riccati equation (FARE) In the LQG controller problem, the optimal controller () in (12), which minimizes the cost functional  of (19), is subjected to constrained (15).This is determined, separately, as the deterministic LQR controller problem.
As with the deterministic LQR problem, the optimal control law that minimizes the same quadratic cost functional has a tradeoff between the state cost and control cost [28].However, as with a stochastic dynamic system, such as an optimal control problem of wind-induced vibration of a tall building, the cost functional for the deterministic LQR problem cannot be employed due to the stochastic nature of the state space formulation of the stochastic differential equation, as shown in (12).That is, the ensemble average over all possible realizations of the excitation is considered, so the cost functional  is given by Here,  and  are positive semidefinite and positive definite weighting matrices.The term ()  (), in (19), is a measure of control accuracy and the term ()  () is a measure of control effort.Minimizing , while keeping the system response and the control effort close to zero, requires an appropriate choice of the weighting matrices  and  [15,27].If it is desirable that the system response be small, then large values for the elements of  should be chosen.The selected matrix should be diagonal and the diagonal element should have a large value for any respective state variable to be small [15].If the control energy is required to be small, then large values of the elements of  should be chosen.The unique state-feedback optimal controller (), which minimizes the cost functional  of (19), is determined as follows: Here,  = , where  is the unique symmetric, positive semidefinite solution to the algebraic Riccati equation (ARE), given by It is noted that we can determine the LQR controller feedback gain matrix  and the observer gain Kalman filter matrix , independently.This is the so-called separation principle.In the LQR control problem, the optimal controller is revealed by tuning some weighting matrix  and .If the value of  is relatively large compared to , the state vector () is small, relative to the control () [28].That is, more control force is applied to the main structure.This means that larger values of  result in the poles of the closed-loop system matrix ( −  * ) far left in the s-plane, so the state decays faster to zero state [28].

Numerical Example
This numerical example is from "Numerical Examples" in Solari [24].The tall building's height  = 180 m, width  = 60 m, depth  = 30 m, first modal natural frequency 1 = 0.27 Hz, critical damping ratio = 0.015, and so forth.ℎ = 120 m, (ℎ) = 40.96m/s,  V (ℎ) = 5.39 m/s,  V (ℎ) = 582.48m,   = 16,   = 10,   = 0.8,  ℓ = 0.5,   = 0.5, and so forth.Additional data for the along-wind load and properties of the building were seen in [24].The optimum parameters of ATMD are considered to be of the same value as the passive TMD.The optimum parameters of TMD have a mass ratio of  = 0.01, a tuning frequency of  opt = 1.0, and damping ratio of  opt = 0.05.The numerically simulated along-wind load, along with response without ATMD, is shown in Figures 2 and 3. RMS displacement response, without ATMD, is shown in Figure 3 and is 0.021, a good approximation to that of Solari's closed form response of 0.027 m.It was known that the relative displacement response of ATMD to that of the main structure was five to ten times larger than that of the main structure, demonstrating that Using those values of  and  and an assumed value of  = 1.0 − 008, the estimated reduced response with ATMD, using LQG controller, is shown in Figure 4.The RMS value of the estimated reduced response is 0.0190 m, which shows around a 10% reduction of the response without ATMD.The estimated relative displacement response of ATMD to that of the main structure is shown in Figure 5, with RMS value of 0.2217 m.This is roughly 12 times larger than that of the main structure.The active control force of the LQG controller is shown in Figure 6.If the value of  = 1.0 − 004 is used as an increased measurement noise, then the estimated reduced response with ATMD, using the LQG controller, is shown in Figure 7.The RMS value of the estimated reduced response is 0.0072 m, which is about a 65% reduction of the response without ATMD.The estimated relative displacement response of ATMD to that of the main structure is shown in Figure 8, with an RMS value of 0.0815 m that is roughly 11 times larger than that of the main structure.The active control force of the LQG controller is shown in Figure 9.If a value of  = 1.0 is used, the estimated reduced response of ATMD, using the LQG controller, is shown in Figure 10.The RMS value of the estimated reduced response is 0.00055026 m, which shows around a 97% reduction of the response without ATMD.The estimated relative response of ATMD to that of the main structure is shown in Figure 11, where the RMS value is 0.0041 m.This is around 7 times greater than that of the main structure.The active control force of the LQG controller is shown in Figure 12.As shown in the above results, the reduced along-wind responses of a tall building with ATMD using LQG controller are obtained as LQG control design parameters are varied.Therefore, LQG control for reducing along-wind response of a tall building with ATMD is effective for reducing wind-induced vibrations endured by tall buildings.

Conclusions
The optimal control technique of LQG for obtaining the reduced along-wind responses of a tall building with ATMD has been investigated.The fluctuating along-wind load, acting on a tall building, was simulated numerically and using this simulated along-wind load, we calculated the along-wind responses of a tall building without ATMD and with ATMD, using the LQG controller.The estimated along-wind response without ATMD is good when considered to the closed form response.The reduced across-wind responses, controlled with ATMD using the LQG controller, were estimated as varying the LQG control design parameters.Therefore, ATMD system using LQG controller is effective and useful for the design of mitigating wind-induced vibrations endured by tall buildings.