# Discrete-Time Multioverlapping Controller Design for Structural Vibration Control of Tall Buildings under Seismic Excitation

Authors: Francisco Palacios-Quiñonero; Josep Rubió-Massegú; Josep M. Rossell; Hamid Reza Karimi
Journal: Mathematical Problems in Engineering
 (2012)

---

## Abstract

In this paper, a computationally effective strategy to obtain multioverlapping controllers via the Inclusion Principle is applied to design discrete-time state-feedback multioverlapping LQR controllers for seismic protection of tall buildings. To compute the corresponding control actions, the proposed semidecentralized controllers only require state information from neighboring stories. This particular configuration of information exchange allows introducing a dramatic reduction in the transmission range required for a wireless implementation of the communication system. To investigate the behavior of the proposed semidecentralized multioverlapping controllers, a proper simulation model has been designed. This model includes semiactive actuation devices with limited force capacity, control sampling times consistent with the communication latency, time-delayed state information, and communication failures. The performance of the proposed multioverlapping controllers has been assessed through numerical simulations of the seismic response of a 20-story building with positive results.


---

## Body

## 1.  Introduction

Over the last decades, problems of ever increasing complexity have been considered in the field of Structural Vibration Control (SVC). Current SVC systems for seismic protection of tall buildings can involve a large number of sensors and actuation devices and a wide and sophisticated communication network [1–3]. Semidecentralized control strategies, which can operate using only state information from neighboring stories, are especially relevant for wireless implementations of the communication system. Semidecentralized state-feedback LQR controllers were proposed by Wang and Lynch in [4], and the study was extended to state-feedback H∞ controllers by Wang et al. in [5]. The numerical and experimental results obtained in these works clearly indicate that the proposed semidecentralized control strategies are specially suitable for SVC of tall buildings with wireless communications. It should be highlighted, however, that important computational difficulties can arise when applying these control design strategies to large buildings. The LQR controller design presented in [4] uses a variant of the heuristic iterative procedure proposed by Lunze in [6], and the H∞ controller design in [5] is based on a Linear Matrix Inequality formulation (LMI). For large-dimensional problems, a great computational effort is required by the iterative procedure used in the LQR design. Analogously, solving large-dimensional convex optimization problems with LMI constraints is also a costly computational task.In this context, the design of semidecentralized controllers using multioverlapping decompositions based on theInclusion Principle (IP) is a very interesting option [7–12]. Broadly speaking, the IP allows decomposing the original large-dimensional problem into a set of low-dimensional decoupled problems. This decomposition takes advantage of the particular structure of the original system and can help to significantly reduce the computational effort. Examples of successful applications of the IP to SVC can be found in [13–15]. Recently, an effective computational strategy to design semidecentralized multioverlapping controllers based on a sequential application of the IP was presented by Palacios-Quiñonero et al. in [16]. In that work, semidecentralized multioverlapping LQR controllers are designed for seismic protection of a four-story building with positive results. However, it has to be noted that all these applications of the IP to SVC have been conducted using small buildings, continuous-time models, and assuming highly idealized conditions, such as active force actuators with unrestricted force capacity and communication systems with no failures nor delays.The main contribution of the present paper is to present a large-scale application of the IP to the design of semidecentralized controllers for SVC, paying special attention to some aspects of practical relevance. More specifically, the computational strategy proposed in [16] is applied to design discrete-time state-feedback multioverlapping LQR controllers to mitigate the seismic response of a 20-story building. Moreover, to gain a meaningful insight into the behavior of the proposed multioverlapping controllers, the models used in the numerical simulations include some factors of practical relevance such as control sampling rates, realistic implementation of the control actions, time-delayed state information, and communication latency and failures. One of the main difficulties encountered when applying the IP to discrete-time controller design is that the natural structure of the continuous-time model is lost in the discretization process. To overcome this difficulty, the discretization process has been carried out on the expanded decoupled subsystems. The results obtained in the numerical simulations confirm the excellent characteristics of the proposed semidecentralized multioverlapping controllers for SVC of large buildings.The organization of the paper is as follows: In Section2, a detailed derivation of the continuous-time state-space model for an n-story building is presented. Section 3 begins with a summary discussion on the design of discrete-time state-feedback centralized LQR controllers. Next, the main ideas involved in the design of discrete-time state-feedback multioverlapping LQR controllers are briefly presented. In Section 4, three mathematical models used to conduct the numerical simulations of the building seismic response are presented: (i) Basic building model, which consists in a discrete-time approximation of the continuous-time state-space model with a small sampling time and no control action. (ii) Centralized control model, which implements the discrete-time centralized LQR controller with perfect state knowledge, small sampling time, and ideal semiactive force actuators with limited force capacity. (iii) Multioverlapping control model. This case implements a discrete-time multioverlapping controller considering semiactive actuation devices with limited force capacity, a control sampling time consistent with the communication latency, time-delayed state information, and communication failures. Finally, in Section 5, the control design methodology presented in Section 3 and the simulation models introduced in Section 4 are applied to a particular 20-story building to assess the performance of the proposed multioverlapping controllers.


## 2. Continuous-Time Building Model

Let us consider then-story building schematically displayed in Figure 1, which is modeled as a lumped-mass planar system with displacements in the direction of the ground motion. The building motion can be described by the second-order differential equation:
(2.1)Mq¨(t)+Cq˙(t)+Kq(t)=Tuu(t)+Tωω(t),
where
(2.2)q(t)=[q1(t),…,qn(t)]T
is the vector of story displacements with respect to the ground, and qi(t) represents the displacement of the ith story. M, C, K are the mass, damping, and stiffness matrices, respectively. The vector of control actions is
(2.3)u(t)=[u1(t),…,un(t)]T,
where ui(t) represents the control force exerted by the actuation device ai (see Figure 2), and Tu is the control location matrix. The seismic ground acceleration is ω(t), and Tω denotes the disturbance input matrix. The mass and stiffness matrices have the following structures:
(2.4)M=[m1⋯⋯mn],(2.5)K=[k1+k2-k2-k2k2+k3-k3⋯⋯⋯⋯⋯⋯-kn-1kn-1+kn-kn-knkn],
where mi and ki represent, respectively, the mass and stiffness of the ith story. When the values of the story damping coefficients ci, 1≤i≤n, are known, a damping matrix C with the same structure as K can be obtained by replacing ki by ci in (2.5). Alternatively, a tridiagonal damping matrix in the form
(2.6)C=α0M+α1K
can be computed following the Rayleigh dampingapproach by setting the damping ratio values for two selected natural frequencies [17]. For the actuation system schematically depicted in Figure 2, the control location matrix has dimensions n×n and the following structure:
(2.7)Tu={[Tu]i,i=1,for1≤i≤n,[Tu]i,i+1=-1,for1≤i<n,[Tu]i,j=0,otherwise,
where [Tu]i,j denotes the element in the ith row and jth column of Tu. Finally, the disturbance input matrix is
(2.8)Tw=-M[1]n×1,
where [1]n×1 is a column vector of dimension n with all its entries equal to 1.Figure 1
Building lumped-mass model.Figure 2
Actuation scheme.Now, we take the state vector(2.9)xI(t)=[q(t)q˙(t)]
and derive a first-order state-space model
(2.10)SI:x˙I(t)=AIxI(t)+BIu(t)+EIω(t).
The state, control, and disturbance input matrices are, respectively,
(2.11)AI=[[0]n×nIn-M-1K-M-1C],BI=[[0]n×nM-1Tu],EI=[[0]n×1-[1]n×1],
where [0]r×s represents a zero matrix of dimensions r×s, and In is the identity matrix of dimension n. Next, we consider a new state vector
(2.12)x(t)=[x1(t),…,x2n(t)]T,
which groups together the interstory drifts and interstory velocities in increasing order
(2.13)x(t)={x1(t)=q1(t),x2(t)=q˙1(t),x2i-1(t)=qi(t)-qi-1(t),for1<i≤n,x2i(t)=q˙i(t)-q˙i-1(t),for1<i≤n.
Using the change of basis
(2.14)x(t)=PxI(t),
defined by the 2n×2n matrix
(2.15)P={p1,1=1,p2,n+1=1,p2i-1,i-1=-1,p2i-1,i=1,for1<i≤n,p2i,n+i-1=-1,p2i,n+i=1,for1<i≤n,pi,j=0,otherwise,
we obtain the new state space model:
(2.16)S:x˙(t)=Ax(t)+Bu(t)+Eω(t),
where
(2.17)A=PAIP-1,B=PBI,E=PEI.
In this work, we restrict our attention to the interstory drifts as output variables. The output vector can then be obtained as
(2.18)y(t)=[y1(t),…,yn(t)]T=Cyx(t),
where Cy is a matrix of dimensions n×2n with the following structure:
(2.19)Cy={[Cy]i,2i-1=1,for1≤i≤n,[Cy]i,j=0,otherwise.


## 3. Control Design

To design a discrete-time centralized state-feedback LQR controller for then-story building model presented in the previous section, we begin by considering the continuous-time system:
(3.1)Sc:x˙(t)=Ax(t)+Bu(t),
obtained from (2.16) by removing the disturbance term Eω(t). The discrete-time system corresponding to the zero-hold approximation of (3.1) with sampling time τ is
(3.2){Sc}τ:x[k+1]=Aτx[k]+Bτu[k],
where
(3.3)Aτ=eAτ,Bτ=∫0τ‍eAtBdt.
Next, we consider the discrete-time state-feedback controller:
(3.4)u[k]=-Gτx[k]
and the quadratic index:
(3.5)J(x,u)=∑k=0k=∞‍x[k]TQx[k]+u[k]TRu[k],
where Q is a symmetric positive semidefinite matrix, and R is a symmetric positive definite matrix. The control gain matrix Gτ that minimizes (3.5) under constraints (3.2) and (3.4) can be computed as
(3.6)Gτ=(R+BτTPBτ)-1BτTPAτ,
where P is the solution of the discrete-time Riccati equation:
(3.7)AτTPAτ-P+Q-AτTPBτ(R+BτTPBτ)-1BτTPAτ=0.To design a multioverlapping controller that is able to compute the control actionsui[k] using only state information corresponding to neighboring stories, we consider the n-story building decomposed into a sequence of n-1 two-story overlapped subsystems
(3.8)S(i)=[si,si+1],1≤i≤n-1,
where si represents the ith story. This overlapping decomposition is schematically depicted in Figure 3. Following the sequential multioverlapping decomposition strategy proposed in [16], the initial continuous-time system (3.1) can be conveniently expanded to form a new continuous-time system:
(3.9)S~:x~˙(t)=A~x~(t)+B~u~(t),
where the state matrix A~ and the control input matrix B~ are block diagonal. The expanded system S~ can then be decomposed into a sequence of decoupled continuous-time subsystems:
(3.10)S~(i):x~˙(i)(t)=A~(i)x~(i)(t)+B~(i)u~(i)(t),1≤i≤n-1.
For the continuous-time subsystems S~(i), we compute discrete-time zero-hold approximations with sampling time τ^:
(3.11){S~(i)}τ^:x~(i)[k+1]=A~τ^(i)x~(i)[k]+B~τ^(i)u~(i)[k],1≤i≤n-1,
where
(3.12)A~τ^(i)=eA~(i)τ^,B~τ^(i)=∫0τ^‍eA~(i)tB~(i)dt
and consider the local quadratic indexes:
(3.13)J~(i)(x~(i),u~(i))=∑k=0k=∞‍{x~(i)[k]}TQ(i)x~(i)[k]+{u~(i)[k]}TR(i)u~(i)[k],1≤i≤n-1,
to compute local discrete-time LQR controllers
(3.14)u~(i)[k]=-G~τ^(i)x~(i)[k],1≤i≤n-1,
which minimize the indexes (3.13) under constraints (3.11) and (3.14). Finally, the sequence of expanded local control matrices G~τ^(i) is contracted back to a control gain matrix G^τ^ in order to define a discrete-time multioverlapping controller:
(3.15)u^[k]=-G^τ^x[k]
for the original discrete-time system (3.2).Figure 3
Decomposition in two-story overlapping subsystems.Remark 3.1.
The expansion-contraction procedure associated to the design of multioverlapping controllers for large buildings is only outlined in this section. For clarity and simplicity, a detailed account of this procedure has not been included in the paper. However, a complete presentation of this background material together with some practical applications to SVC of small buildings can be found in [15, 16].Remark 3.2.
The expanded block-diagonal system (3.9) can only be computed when the matrices of the initial state-space system have a suitable zero-nonzero block structure. For the building model (2.1), the initial state-space system (3.1) has a proper structure. However, this structure is lost in the discretization process and the expansion-decoupling process can no longer be applied to the discrete-time state-space system (3.2). To overcome this difficulty, the expansion-decoupling process is first completed for the continuous-time system and, after that, the discretization process is carried out on the continuous-time expanded decoupled subsystems (3.10) to obtain the discrete-time expanded decoupled subsystems (3.11).Remark 3.3.
It should be noted that the control gain matrixGτ=[(gτ)i,j] given in (3.6) is a full matrix of size n×2n, and the full state is required to compute the control action for the actuation device ai:
(3.16)ui[k]=∑j=1j=2n‍(gτ)i,jxj[k],1≤i≤n.
In contrast, the multioverlapping control matrix G^τ^=[(g^τ^)i,j] has a block-tridiagonal structure and only requires a reduced number of 4–6 states to compute the control action for each actuation device. More specifically, we have
(3.17)u^[k]={u^1[k]=∑j=1j=4‍(g^τ^)1,jxj[k],u^i[k]=∑j=2i-3j=2i+2‍(g^τ^)i,jxj[k],for1<i<n,u^n[k]=∑j=2n-3j=2n‍(g^τ^)n,jxj[k].Remark 3.4.
For clarity and simplicity, the controllers presented in this section have been computed following an LQR approach. However, it has to be highlighted that other control strategies are also possible. For example, an application of the IP to the design of semidecentralized static output-feedback controllers for SVC can be found in [14].


## 4. Simulation Models

One of the main objectives of the present work is to gain a meaningful insight into the behavior of semidecentralized multioverlapping controllers through numerical simulations. To this end, the simulation models have to include some relevant factors such as sampling rates, realistic implementation of the control actions, time-delayed state information, and communication latency and failures. Trying to achieve a proper balance between simplicity and accuracy, we have considered a simulation framework formed by three different models: (i)Basic building model, which consists in a discrete-time approximation of the continuous-time state-space model (3.1) with a small basic sampling time τ and no control action. (ii) Centralized control model, which implements the discrete-time LQR controller given in (3.4) with perfect state knowledge and the basic sampling time τ. The actuation devices, however, are assumed to be ideal semiactive force actuators with limited force capacity. (iii) Multioverlapping control model. This case implements the discrete-time multioverlapping controller given in (3.15) considering semiactive actuation devices with limited force capacity, a control sampling time τ^>τ consistent with the communication latency, time-delayed state information, and communication failures. In all the cases, the interstory drifts are taken as output variables.The basic building model is:(4.1)ℳτ:{x[k+1]=Aτx[k]+Eτωτ[k],y[k]=Cyx[k],
where Aτ is the discrete-time state matrix in (3.3); Eτ is the discrete-time disturbance input matrix, which can be computed as
(4.2)Eτ=∫0τ‍eAtEdt,
with A and E representing, respectively, the state and disturbance input continuous-time matrices; and ωτ[k]=ω(kτ) is the sampled disturbance. The vector of interstory drifts is computed with the output matrix Cy given in (2.19). A good approximation of the uncontrolled seismic response of the building can be obtained using the basic building model ℳτ with a small sampling time τ.The centralized control model is:(4.3)ℳ-τ:{u[k]=-Gτx[k],x[k+1]=Aτx[k]+Bτσ(u[k])+Eτωτ[k],y[k]=Cyx[k],
where the sampling time τ, the matrices Aτ, Eτ, and Cy, and the sampled disturbance ωτ[k] are the same as those used in (4.1); and Bτ is the discrete-time input-control matrix defined in (3.3). The vector of control actions
(4.4)u[k]=[u1[k],…,un[k]]T,
is computed using the discrete-time centralized control gain matrix Gτ given in (3.6), and the vector of control forces is
(4.5)σ(u[k])=[σ(u1[k]),…,σ(un[k])]T.
In this section, the actuation devices ai are modeled as ideal semiactive force actuators with maximum actuation force [fmax]i. For a given control action ui[k], the actual control force exerted by the actuation device ai is
(4.6)σ(ui[k])={ui[k]ifui[k]·vi[k]<0,and|ui[k]|≤[fmax]isgn(ui[k])·[fmax]iifui[k]·vi[k]<0,and|ui[k]|>[fmax]i0ifui[k]·vi[k]≥0,
where vi[k] is the corresponding interstory velocity, and sgn(x)=x/|x| is the signum function. In the centralized control model, we assume an ideal communication system, which can provide a perfect knowledge of the full state vector x[k]. This model is used as a reference in the performance assessment of the multioverlapping controller.In the multioverlapping control model, we consider the actuation-communication system schematically depicted in Figure4, consisting in an actuation device ai, a sensor unit s^i, a local control unit c^i, and a wireless communication unit w^i. The actuation device ai produces the semiactive implementation of the control actions defined in (4.6). The sensor unit s^i is an ideal sensor that provides an exact measurement of the local state
(4.7)x^i[k]=[yi[k]vi[k]],
where yi[k] and vi[k] denote the local interstory drift and velocity, respectively. To model the operation of the local control unit, we introduce the controller sampling time
(4.8)τ^=r^τ,r^>1,
where τ is the basic sampling time used in (4.1) and (4.3), and r^ can be understood as the maximum number of sampling steps that can be spent by c^i to collect state information from neighboring stories through the communication unit w^i. We also consider the update-control times
(4.9)k^j=jr^,j≥1,
and define the interval of state-information gathering
(4.10)Δg[k^j]=[k^j-1+1,…,k^j],
and the interval of control-action holding
(4.11)Δh[k^j]=[k^j,…,k^j+1-1],
which are schematically represented in Figure 5. The operation of the local control unit c^i has been modeled in accordance with the following set of basic principles: (P. 1)
The local control actionu^i[k] is updated at the sampling times k=k^j, j≥1.(P. 2)
The local controller unitc^i has direct access to the sensing unit s^i; consequently the local state x^i[k^j] is assumed to be always available.(P. 3)
The state information of neighboring stories obtained through the wireless communication unitw^i has the form x^i'[k-δi′[k^j]], where i'=i±1 and the time delay satisfies 0≤δi'[k^j]<r^.(P. 4)
For a given time intervalI with length τ^, the events Ei,i′(I,τ^)=[c^i obtains the neighboring state x^i′ in the time interval I of length τ^] are independent random events with common probability
(4.12)Prob[Ei,i′(I,τ^)]=pτ^.(P. 5)
IfI, I' are non-overlapping time intervals with respective lengths τ^ and τ^′, then Ei,i'(I,τ^) and Ei,i'(I',τ^′) are independent random events.(P. 6)
Through the time intervalΔg[k^j], the local control unit c^i tries to collect the state information required to compute the control action. If this state information is successfully acquired, the flag variable ϕi[k^j] is set to 1 and the new control action u^i[k^j] is computed; otherwise, the flag variable ϕi[k^j] and the control action u^i[k^j] are both set to 0.(P. 7)
The control action computed at the sampling timek=k^j is held through the time interval Δh[k^j]. According to the previous principles, the vector of control actions
(4.13)u^[k]=[u^1[k],…,u^n[k]]T
can be computed by setting the initial value
(4.14)u^[0]=[0]n×1,
and, for k>0, using the expression
(4.15)u^[k]={-Fϕ[k](DG^τ^x[k]+LG^τ^x[k-δ[k]])ifmod(k,r^)=0,u^[k-1],otherwise,
where mod(k,r^) represents the integer remainder after division, Fϕ[k] is the diagonal matrix
(4.16)Fϕ[k]=[ϕ1[k]⋱ϕn[k]],DG^τ^ is the block-diagonal matrix
(4.17)DG^τ^=[[G^τ^]1,1⋱[G^τ^]n,n]
formed by the diagonal blocks of G^τ^(4.18)[G^τ^]i,i=[{g^τ^}i,2i-1,{g^τ^}i,2i],1≤i≤n,
and the matrix
(4.19)LG^τ^=G^τ^-DG^τ^
contains the out-of-diagonal blocks of the block-tridiagonal multioverlapping control matrix G^τ^. The notation x[k-δ[k]] represents the delayed state
(4.20)x[k-δ[k]]=[x^1[k-δ1[k]]⋮x^n[k-δn[k]]],
where δi[k] is the delay in the local state x^i.Figure 4
Actuation-communication system for the multioverlapping control model.Figure 5
Time intervals for state information gathering and control action holding.The multioverlapping control model can now be obtained by completing (4.14) and (4.15) with the state and output equations
(4.21)x[k+1]=Aτx[k]+Bτσ(u^[k])+Eτωτ[k],y[k]=Cyx[k].
It should be noted that expressions (4.16) and (4.20) are only evaluated at the update-control times k^j=jr^. Moreover, according to (P. 1), (P. 4) and (P. 6), the information-state flag variables ϕi[k^j] are independent random variables with Bernoulli distributions ℬ(p). In particular ϕ1[k^j] and ϕn[k^j] have a Bernoulli distribution ℬ(pτ^), and ϕi[k^j] has distribution ℬ(pτ^2) for 1<i<n. Finally, it should also be noted that the probability of gathering the neighboring state information x^i' by the controller unit c^i in a time interval I' of length 2τ^ is
(4.22)Prob[Ei,i′(I′,2τ^)]=2pτ^-pτ^2.
This formula can be easily obtained by writing the time interval I' as union of two non-overlapping intervals I'=I1′∪I2′ of length τ^. According to (P. 5), Ei,i'(I1′,τ^) and Ei,i'(I2′,τ^) are independent random events, and the probability of failing to acquire the state x^i' in the whole interval I' is (1-pτ^)2. Analogously, it can be shown that the corresponding probability for a time interval I′′ of length (1/2)τ^ is
(4.23)Prob[Ei,i′(I′′,12τ^)]=1-1-pτ^.


## 5. Numerical Simulations

In this section, the behavior of discrete-time multioverlapping LQR controllers is investigated through numerical simulations of the seismic response of a 20-story building. The parameter values for this particular building are collected in Table1 and are similar to those used in [5]. The damping matrix has been computed as a Rayleigh damping matrix by setting a 5% of damping ratio for the 1st and 18th natural frequencies. The actuation system ai implemented between the (i-1)th and ith stories (see Figure 2) is assumed to be formed by a number of identical actuation devices that work coordinately as a single device. The force saturation level of a single actuation device has been taken as 1.2×106 N. The total number of actuation devices and the maximum actuation force for the actuation systems ai, 1≤i≤n, is also presented in Table 1.Table 1
Particular parameter values for the 20-story building.
Story
1–5
6–11
12–14
15–17
18-19
20
Mass (×106 Kg)
1.10
1.10
1.10
1.10
1.10
1.10
Stiffness (×106 N/m)
8.62
5.54
4.54
2.91
2.56
1.72
Number of actuation devices
4
2
2
1
1
1
Max. actuation force (×106 N)
4.8
2.4
2.4
1.2
1.2
1.2
Natural damping
5%For this 20-story building, three different discrete-time LQR controllers are designed: a centralized controllerGτ, which has been obtained using the basic sampling time τ=10-3 s; and two multioverlapping controllers G^τ^ and G^τ'^, computed with sampling times τ^=40τ, and τ^′=20τ, respectively. The particular values of the weighting matrices in (3.5) used to design the centralized controller are Q=I40 and R=10-17.5×I20. For the multi-overlapping controllers, the weighting matrices in (3.13) used to compute the local expanded LQR controllers G~τ^(i) and G~τ^′(i), 1≤i<20, have been taken as Q(i)=I4 and R(i)=10-17.5×I2, for 1≤i<20.In the numerical simulations, the maximum absolute interstory drifts have been computed for different control configurations. The basic building model given in (4.1) with sampling time τ=10-3 s has been used to compute the uncontrolled seismic response. The controlled response corresponding to the centralized controller Gτ has been obtained with the centralized control model presented in (4.3). Finally, the multioverlapping control model defined in (4.14), (4.15), and (4.21) has been used to compute the controlled response for the multioverlapping controllers G^τ^ and G^τ^′. In all the cases, the full scale 1995 Kobe North-South seismic record has been taken as ground acceleration (see Figure 6). This seismic record, obtained at the Kobe Japanese Meteorological Agency station during the Hyogoken-Nanbu earthquake of January 17, 1995, is a near-field record that presents large acceleration peaks which are extremely destructive to tall structures [18, 19]. In the multioverlapping control model, the reference value of pτ^=0.95 has been set for the probability of obtaining state information from neighboring stories in a time interval of length τ^=40ms. According to (4.23), for the multioverlapping controller G^τ^′ with control sampling time τ^′=20ms, the probability of successfully gathering state information from neighboring stories can be taken as pτ^′=0.78.Figure 6
Full scale Kobe 1995 North-South seismic record.In Figure7(a), the red line with asterisks (Overlap. fail in the legend) presents the maximum absolute interstory drifts corresponding to the multioverlapping controller G^τ^ with controller sampling time τ^=40ms, probability of successful communication pτ^=0.95, and state delay δ=39ms. The blue line with circles (Overlap. in the legend), displays the values obtained with no communication failures, that is, with pτ^=1. The interstory drifts peak values corresponding to the multioverlapping controller G^τ^′ with controller sampling time τ^′=20ms, and state delay δ=19ms are presented in Figure 7(b). Here, the red line with asterisks corresponds to the probability pτ^′=0.78, and the blue line with circles presents again the results for pτ^′=1. In both cases, the graphics corresponding to the uncontrolled response (black line with triangles), and the controlled response for the centralized controller Gτ, with controller sampling time τ=1ms, with no communication failures nor delays (black line with squares) have been included as reference. In Figure 8, the red line with asterisks (Over. delay in the legend) displays the maximum absolute interstory drifts corresponding to the multioverlapping controller G^τ^ with τ^=40ms, pτ^=0.95, and state delay δ=39ms, while the green line with circles (Over. no delay in the legend) presents the response obtained with null state delay.Maximum absolute interstory drifts for the 1995 Kobe North-South seismic record. Simulations with maximum state delay.
(a)
Controller sampling time 40 ms
(b)
Controller sampling time 20 msFigure 8
Maximum absolute interstory drifts for maximum and minimum state delays (controller sampling time 40 ms).The graphics in Figure7(a) show the excellent performance of the proposed multi-overlapping controller for controller sampling times τ^ compatible with moderate values of communication latency and also compatible with moderate rates of communication failures. Moreover, the graphics in Figure 7(b) clearly illustrate the trade-off between the controller sampling time τ^ and the probability of successful communication pτ^. Certainly, taking a smaller sampling time τ^ allows a more accurate implementation of the control actions; however, this also implies a reduction of the probability pτ^ which, in the end, may result in an overall loss of performance.Finally, the graphics in Figure8 show the moderate influence of the state delay δ in the multioverlapping controller performance for reasonable values of the controller sampling time τ^.Remark 5.1.
It is worth to be mentioned that the behavior of the ideal discrete-time centralized controllerGτ is very similar to the behavior exhibited by an ideal continuous-time centralized LQR controller. As mentioned in Remark 3.3, full state information is required by centralized controllers and this fact makes them unsuitable for SVC of large buildings with wireless communication systems. A detailed discussion of this point can be found in [5].Remark 5.2.
The proposed semidecentralized controllers can operate using only state information from neighboring stories. This fact makes it possible for them to successfully collect the required state information in a relatively small time interval. As a side effect, state delays are also small and have no significant impact on the controller performance.Remark 5.3.
Force saturation is an important issue in SVC. For large seismic excitations, the required control actions frequently exceed the force capacity of the actuation devices. Consequently, force actuation constraints should be considered when studying the controllers behavior. All the numerical simulations of the controlled responses presented in this paper have been conducted using the force saturation values displayed in Table1.


## 6. Conclusions and Future Directions

In this paper, a computationally effective strategy has been used to design discrete-time state-feedback multioverlapping LQR controllers for seismic protection of tall buildings. This strategy, based on a sequential application of the Inclusion Principle, produces a block-tridiagonal control gain matrix that allows computing the corresponding control actions using only state information from neighboring stories. Due to this particular information exchange configuration of the multioverlapping controllers, the transmission range and the control sampling frequency in wireless implementations of the communication system can be dramatically improved. To investigate the behavior of the proposed semidecentralized multioverlapping controllers, a proper simulation model has been designed, which allows including semiactive actuation devices with limited force capacity, control sampling times consistent with the communication latency, time-delayed state information, and communication failures. To assess the performance of the proposed multioverlapping controllers, numerical simulations of the seismic response for a 20-story building model have been conducted with positive results.For clarity and simplicity, the controllers presented in this paper have been designed following an LQR approach. In future works, further research effort should be addressed at exploring the effectiveness of the proposed control design strategy in more complex scenarios, which can involve issues of practical interest such as structural information constraints [20], actuator saturation [21], actuation and sensor failures [22], and limited frequency domain [23]. Other natural extensions of the present work should include a deeper treatment of some important practical aspects related to the communication system such as missing measurements [24–28], stochastic uncertainties [29], and stochastic nonlinearities [30–33].



---

*Source: 636878-2012-10-22.xml*