JAM Journal of Applied Mathematics 1687-0042 1110-757X Hindawi Publishing Corporation 10.1155/2015/507306 507306 Research Article Synchronized Control for Five-Story Building under Earthquake Loads Mesbahi Javad 1 Malek Alaeddin 2 Su Ray K. L. 1 Department of Mathematics Payame Noor University Tehran 19395-4697 Iran pnu.ac.ir 2 Department of Applied Mathematics Faculty of Mathematical Sciences Tarbiat Modares University P.O. Box 14115-134, Tehran 1411713116 Iran modares.ac.ir 2015 752015 2015 22 02 2015 14 04 2015 15 04 2015 752015 2015 Copyright © 2015 Javad Mesbahi and Alaeddin Malek. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Synchronized control is implemented for a five-story building under earthquake loads and its capabilities are investigated for protection of building under earthquake. In this regard, we applied H2/LQG control algorithm in form of synchronized control for structural vibration reduction. Simulation results of modeling indicated that not only the provided control is able to reduce the responses of vibrations for the structure, but also it is even capable of supplying the objectives of synchronized control at the same time. Numerical results for uncontrolled, traditional H2/LQG control and synchronized control coupled with H2/LQG algorithm are presented. It is shown that for El Centro and Bam earthquakes the synchronized control is more efficient to reduce damage to the given structures.

1. Introduction

In the recent years, applying structure controlling technologies for reducing the responses of structure, that is, speed, displacement, acceleration, and force under load of earthquakes or severe winds attracted a lot of attentions. Structure control is classified in different forms including active, passive, semiactive and hybrid controls . In passive control strategy, the vibration force is absorbed by damper. In active control systems, large actuators are used to perform direct control force but, in semiactive systems, the control force is applied indirectly to the structure. The semiactive control system gained more attention, because it can control the performance of the structure satisfactorily, and at the same time it needs less energy to achieve the control objectives. The semiactive actuators are including actuator variable stiffness (AVS) and dampers containing Electrorheological Fluids (ER) or Magnetorheological Fluids (MR) [1, 2].

Active and semiactive systems both have actuators, sensors, and controllers in the structure; therefore, they have more intelligent performance than the passive systems for change in the structural and environmental conditions. In this system, the sensor will collect data from structure during the dynamic loading and send it to the controllers. Then, after process by the controller’s algorithm the control force will be determined and will be sent to actuators and the vibrations of structure are controlled.

The objective of controlling algorithms is to determine the control force in optimized way and to deduct the vibration responses of buildings. For example, control algorithms like LQR, pole assignment, sliding mode control, bang-bang control, clipped control, LQG, independent model space control (IMSC), Fuzzy Logic Control (FLC), Genetic Fuzzy Logic Control (GFLC), and Adaptive Neuro-Fuzzy Inference System (ANFIS) can be named. Each one of these systems has specific particulars that are considered to be used based on condition of structures [1, 2].

On the other hand, the reduction of the vibration response should be in such a manner that the internal forces of the structure will reach their minimum amount. In this regard, the control algorithm shall be set in a way that it can minimize relative displacement between freedom levels synchronically. On the other hand, they shall be synchronized. For example, displacement and drift in a multistory structure shall be reduced at the same time.

Synchronization is a very important issue in most of the controlling systems and it has different definitions [3, 4]. Motion synchronization in multiagent systems with different coordinates and same objective is of great importance . In such applications of control, the performance of system is more dependent on motion synchronization than the correctness of motion in each independent part. The need for synchronized motion will feel when the parts with different coordinates are supposed to move at the same time and the relationship between their relative displacements is important . In most of the experiments while the start and stop times should be synchronized, also the relationship between location and speed between coordinates is very important. The need for improvement in such systems will lead to performing wide researches in this area. The initial researches are in field of motion and indicated that the general performance may be improved through adjustment of error between independent agents .

In most of the nonvibrational control models, some agents have the tracking objective but, in vibrational control models, the objective is regulation of outputs or reduction of structure vibration. In fact regulation is a type of tracking but the amount of objective or expected objective is zero . Therefore, based on the explanation, the regulated outputs of the model should be defined in such a manner that they can satisfy the objectives of vibration reduction and, at the same time, they can synchronize the targets.

In this study, a five-story model building was examined to compare control and a traditional active control algorithm H2/LQG against a mathematical synchronized algorithm coupled with H2/LQG. This novel form of synchronized control algorithm combines H2/LQG and mathematical synchronized parameters in a novel algorithm. By using this synchronized control algorithm, displacement and relative displacement of floors, as synchronized parameters under earthquake loading is minimized. With the knowledge of the authors there is no previous research in which a structural seismic control is regarded as motion synchronized control problem or in similar perspective.

2. Synchronized Control Theory

Consider a multiagent system in which n agents are participating in one process. Therefore we need a synchronizer for all of the agents. The target of synchronized control is synchronizing all of the agents in such a manner that the agents can maintain a kinematic relationship, that is, the objective of synchronizer .

Setting multiagents for maintain kinematic relationship may be performed in form of guiding and locating the agents along the boundary (or curves) in one specific path.

Definition 1.

Let S(ξ,t) be a function with boundary Γ(ξ,t) in a compact set ΩR2, where ξΩ stands for a state vector and t the time.

Definition 2.

Define xit and xid(t) as general state and desired state variables in the ith story.

Definition 3.

Let xid(t) be an arbitrary point on Γ, where Γ(xid,t)=0. Denote ei(t) as an error of the state variable in the ith story by(1)eit=xit-xidt,i=1,2,3,4,5.

Definition 4.

Let an earthquake happen in the period of t[t0,tmax], where t0 is an initial time and tmax is the final time of earthquake.

Definition 5.

The goal of synchronized control is achieved when for each agent the state variables xi(t) converge to the desired values xid(t), as ttmax; that is, e0 and Γ(xid(t),tmax)=0.

The formulization of synchronizing is different for various issues and for each issue it can be defined in a certain form .

3. <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M29"><mml:mrow><mml:msub><mml:mrow><mml:mi>H</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>/LQG Synchronizing Control Algorithm with Feedback of Acceleration

Equation of motion of a structure under dynamic loads is in the form of a second-order differential equation as follows :(2)Mq¨t+Cq˙t+Kqt=B0ut+Lwtin which M, C, and K are matrices of mass, damping, and stiffness of the system, respectively. q(t) is displacement vector, u(t) is control force vector, and w(t) is external disturbance vector (resulting from earthquake or wind load). B0 and L are, respectively, force control place and external stimulation matrices. And the state-space model of the equation of motion depends on the choice of coordinates.

In order to gain nodal model of (2), assuming that the mass matrix is not unit, the equation can be rewritten as follows :(3)q¨t+M-1Cq˙t+M-1Kqt=M-1B0ut+M-1Lwt.If we define displacement and velocity as variables of state space, xt=q(t)q˙(t), the state equation will be as follows:(4)x˙t=Axt+B2ut+B1wt,where A=0I-M-1K-M-1C is state matrix, B2=0M-1B0 is actuator stimulus place matrix, and B1=0M-1L is disturbance stimulation place matrix.

Consider the following controlling system. It uses feedback acceleration of a system as output :(5)x˙t=Axt+B1wt+B2ut,zt=C1xt+D11wt+D12ut,yt=C2xt+D21wt+D22ut+v,where z consists of the system evaluation parameters and y is an output of ideal sensors without any noise. D11, D12, and C1 are corresponding parameters of z(t) as a regulated output for minimizing the value of the cost function.

They are defined as acceleration, velocity, and displacement (or relative displacement) coefficients in matrix or vector forms. D21, D22, and C2 are corresponding matrices and a vector for sensors outputs. They are defined based on which sensors are selected:(6)J=limτ1τE0τzTtQzt+uTtRutdt,where E is the expected value and Q and R are positive semidefinite matrices to regulate the evaluation outputs and weighting the control force, respectively.

In the following instead of J, two Riccati equations are used; without loss of generality one may assume that the active control force u(t) is proportional to the estimated state variable x^t. That is,(7)ut=-Kcx^t,where Kc is defined in a way to minimize (6) by proper u(t). For this respect an algebraic Riccati equation which is corresponded with (6) should be solved; x^ is estimated state variable and Kc=B2P and B2 is defined in (5) and P is the solution of the following Riccati equation to gain controller gain matrix :(8)ATP+PA-PB2R-1B2TP+Q=0.This equation is solved, through introduction of known parameters in MATLAB Control Toolbox; by different Q matrices, a set of different controllers is gained. In order to gain suitable controller for system, it is needed to give proper weights to Q matrix. In this regard, we can weigh important evaluation parameters for the system and place them at one diagonal matrix. This matrix is a weighted matrix.

In the next step standard Kalman filter is used to estimate the system state variables :(9)x^˙t=A-KKFC2x^t+KKFyt+B2-KKFD22ut,where KKF=P~C2T+B1QKFD21TRKF+D21QKFD21T-1 is the estimator gain matrix and P~ is gained from the below Riccati equation to get the observer gain matrix:(10)AP~+P~AT-P~C2T+B1QKFD21T·RKF+D21QKFD21T-1C2P~+D21QKFB1T=-B1QKFB1T,where QKF is magnitude of stimulation spectrum density entering system Sdd(ω) and RKF is magnitude of noise spectral density of measurement. It is supposed that the stimulates entering system are d, sensor noises are v, and they are not correlated. Thus H2/LQG controller will be designed through gaining K and KKF through using MATLAB Control Toolbox. ut is gained using (7) and (9) as follows:(11)x^˙t=Acx^t+Bcyt,ut=Ccx^t,where(12)Ac=A-KKFC2-B2K+KKFD22K,Bc=KKF,Cc=-Kc,Dc=0,where Ac, Bc, Cc, and Dc are the controller dynamic. Now, the evaluation parameters of control algorithm shall be defined in such a manner that, beside the target of output regulation, the synchronization is gained as well. Considering the fact that in structure control the objective is to control structure vibrations, therefore the structure’s displacements in different degree of freedoms are considered as state variable error and it is introduced as part of the regulated outputs. The relative displacements between these degrees of freedoms are defined as objective of synchronized control. The objective of synchronized control will lead to synchronizing in displacements and reduction of internal forces.

4. The Study Model

In this part, a five-story structure model is considered that is similar to Kajima Shizuoka Building . This building has five active hydraulic actuators that are placed between each two adjacent floors (see Figure 1). First story height is 4.2 m and height of each of the remaining stories is 3.6 m. The motion equation of this building can be presented by (2), that is, a second-order differential equation.

A five-story structure model similar to Kajima Shizuoka Building [8, 9].

In this case, qtR5 is displacement vector, u(t)R5 is control force vector, and w(t)R is disturbance vector (caused by earthquakes or wind load). B0 and L are control force place and external disturbance place-matrices, respectively, M, B0, L, C, and K are defined as follows:(13)M=103×215.200000209.200000207.000000204.800000266.1,B0=1-100001-100001-100001-100001,L=-M11111,C=103×650.4-231.1000-231.1548.9-202.5000-202.5498.6-182.0000-182.0466.7-171.8000-171.8318.5,K=106×260-113000-113212-99000-99188-89000-89173-84000-8484.The natural damping is considered as %5 and the natural frequencies of the structure are as follows:(14)ωi=42.5423,36.4257,28.2253,17.7417,6.3343rad/sec,where the mass in kilograms (kg), damping coefficient in Ns/m, and roughness factor are N/m. In the modeling, only acceleration of first, third, and fifth stories is used as the partial-state feedback to compute control force.

In order to implement the H2/LQG control algorithm in the form of synchronized control, in this model, the error of state variables is determined and introduced them in regulated outputs. In this regard, we consider the displacement of each floor as state variable for error determination. The expected amounts for each component of the error vector are zero.

Therefore, in the optimistic case, the error of state variables in each moment is equal to the amount of displacement in each floor and thus the objective of control is to tend these values to zero as ttmax (Figure 2). Consider(15)e=x1,x2,x3,x4,x5.Let us to assume that ε5×1 is a vector that each component stays on the final desired curve as ttmax; that is, Γ(xid,tmax)=0. By this approach, the relevant displacement between the floors, ε is considered as synchronizing error (see Figure 2). Relevant displacement between the floors (drift) is determined as follows:(16)ε=x1,x2-x1,x3-x2,x4-x3,x5-x4.Reduction of relevant displacement between the floors will lead to reducing internal forces of structure. Note that the coupled position error is defined by(17)Ec=I+αTe,where T satisfies(18)ε=Teand without loss of generality one may assume(19)T=1-100001-100001-100001-100001,where α is an identical matrix. This means that (I+αT) is an invertible matrix and as e converges to zero Ec converges to zero too and vice versa. With respect to the definitions, as ttmax coupled position error Ec, e and ε can converge to zero simultaneously .

Displacement and relevant displacement of the floors.

Now, the synchronizing H2/LQG controller is designed by considering e and Ec as the evaluated outputs z from (5). With this one can simultaneously gain the output setting objectives together with synchronized goals. These regulated outputs defined in synchronization H2/LQG controllers are different from H2/LQG controllers. That is, in the traditional algorithm there are no synchronization parameters Ec and ε. These two algorithms are the same in using first, third, and fifth stories displacement, drifts, and acceleration to compute the regulated outputs. They are different in using different related parameters.

The abovementioned structure is modeled in MATLAB in Simulink environment (see Figure 3). In this model, at first the record of earthquake w(t) is inserted into structure model and then the vibration responses of structure are transferred to controller through acceleration sensors y(t). Finally, proper control force u(t) is applied to the structure and the excessive vibrations are controlled.

5. Numerical Results and Discussions

In order to test the designed control algorithm of the structure model in Figure 1, it is analyzed under El Centro Earthquake 1940 and Bam Iran Earthquake 2003 in two forms, with and without controller. The results of analysis are represented for output regulation and synchronizing objectives in structure’s time history response graphs. Figure 4 shows control force level under El Centro Earthquake for all stories. Figures 5, 6, and 7 are demonstrating displacement output, acceleration, and drift ratio (ratio of story’s drift to the height of the corresponded story) of 5th floor of abovementioned structure in two ways of controlled and uncontrolled under El Centro Earthquake. The reduction values are %92, %94, and %88.7 for the displacement, interstory drift, and acceleration of fifth floor, respectively.

Time history graphs of control force for all stories under El Centro Earthquake (1940), when synchronized control coupled with H2/LQG is used.

Time history graphs of displacement for the fifth floor for uncontrolled and synchronized control coupled with H2/LQG algorithm in El Centro Earthquake (1940).

Time history graphs of acceleration for the fifth floor for uncontrolled and synchronized control coupled with H2/LQG algorithm in El Centro Earthquake (1940).

Time history graphs of drift ratio of fifth floor for uncontrolled and synchronized control coupled with H2/LQG algorithm in El Centro Earthquake (1940).

Figure 8 shows control force level under Bam Earthquake for all stories. Figures 9, 10, and 11 are demonstrating displacement output, acceleration, and drift ratio of 5th floor of abovementioned structure in two ways of controlled and uncontrolled under Bam Earthquake. The reduction values are %73, %73.5, and %93.64 for the displacement, interstory drift, and acceleration of fifth floor, respectively.

Time history graphs of control force for all stories under Bam Earthquake when synchronized control coupled with H2/LQG is used.

Time history graphs of displacement of fifth floor for uncontrolled and synchronized coupled with H2/LQG algorithm in Bam Earthquake.

Time history graphs of acceleration of fifth floor for uncontrolled and synchronized control coupled with H2/LQG algorithm in Bam Earthquake.

Time history graphs of drift ratio of fifth floor for uncontrolled and synchronized control coupled with H2/LQG algorithm in Bam Earthquake.

Considering Figures 57 and Figures 911, it is clear that the proposed control algorithm was able to prevent excessive vibration of the structure and this reduces the internal forces in the structure. It also reduces the damages to structures during the period of earthquake.

On the other hand, the objectives of synchronized control are supplied in addition to regulation of structure output; that is, with convergence of E, e will vanish as ttmax. It is demonstrated in Figures 57 and Figures 911 that the uncontrolled and the controlled solutions for El Centro Earthquake and for Bam Earthquake vanish as ttmax=60 seconds.

Reduction of displacement in floors and convergence toward zero occurred for all of the floors as error of state variable in the same manner that happens for 5th floor. Figures 12 and 13 display the displacements of all five floors together when both under controlled and uncontrolled forms are considered. On the other hand, after converging, the relative displacement ε for all floors with the controller vanishes. Figures 14 and 15 demonstrated the drift ratio of all of the five floors in controlled and uncontrolled ways.

Displacements of all five floors for uncontrolled and synchronized control coupled with H2/LQG algorithm in El Centro Earthquake (1940).

Displacement of all five floors for uncontrolled and synchronized control coupled with H2/LQG algorithm in Bam Earthquake.

Drift ratio of all five floors for uncontrolled and synchronized control coupled with H2/LQG algorithm in El Centro Earthquake (1940).

Drift ratio of all five floors for uncontrolled and synchronized control coupled with H2/LQG algorithm in Bam Earthquake.

For comparison of synchronizing H2/LQG controller and traditional H2/LQG controller, the maximum absolute interstory drift and the maximum control force diagrams are supplied in Figure 16 for the El Centro (1940) earthquake and Figure 17 for the Bam earthquake. (Note that Figure 16 is not scaled, and for the scaled El Centro ground motion to 1 the response would be 3.417 times smaller.)

Maximum absolute interstory drifts and control efforts in El Centro Earthquake (1940).

Maximum absolute interstory drifts and control efforts in Bam Earthquake.

The synchronized H2/LQG controller demonstrates its superiority compared to traditional H2/LQG controller in controlling interstory drifts and also in control force.

6. Conclusion

In this study, we proposed a mathematical synchronized control algorithm coupled with H2/LQG for five-story buildings. The simulation results indicate that, defining the state variable error and synchronized error in regulated outputs, the objectives of synchronized control algorithm can be gained. In addition, after setting these two errors, displacement, drifts, and acceleration values of each floor will decrease significantly and, therefore, the amount of forces on each floor and internal forces of structure will decrease. In general it can be concluded that if conditions for synchronization and regulation are defined properly, the structure can be protected against damage more efficiently.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Housner G. W. Bergman L. A. Caughey T. K. Chassiakos A. G. Claus R. O. Masri S. F. Skelton R. E. Soong T. T. Spencer B. F. Yao J. T. P. Structural control: past, present, and future Journal of Engineering Mechanics 1997 123 9 897 971 10.1061/(asce)0733-9399(1997)123:9(897) 2-s2.0-0030724284 Fallah A. Y. Taghikhany T. Time-delayed decentralized H2/LQG controller for cable-stayed bridge under seismic loading Structural Control and Health Monitoring 2013 20 3 354 372 10.1002/stc.498 2-s2.0-84872932215 Zhao D. Zhu Q. Li S. Gao F. Synchronized control of mechanical systems: a tutorial Applied Methods and Techniques for Mechatronic Systems 2014 452 Berlin, Germany Springer 1 25 Lecture Notes in Control and Information Sciences 10.1007/978-3-642-36385-6_1 MR3205023 Sun D. Synchronization and Control of Multiagent Systems 2010 41 New York, NY, USA CRC Press Zhao D. Zhu Q. Position synchronised control of multiple robotic manipulators based on integral sliding mode International Journal of Systems Science 2014 45 3 556 570 10.1080/00207721.2012.724106 MR3172833 2-s2.0-84884887383 Sun D. Shao X. Y. Feng G. A model-free cross-coupled control for position synchronization of multi-axis motions: theory and experiments IEEE Transactions on Control Systems Technology 2007 15 2 306 314 10.1109/tcst.2006.883201 2-s2.0-33947609593 Sun D. Wang C. Shang W. Feng G. A synchronization approach to trajectory tracking of multiple mobile robots while maintaining time-varying formations IEEE Transactions on Robotics 2009 25 5 1074 1086 10.1109/TRO.2009.2027384 2-s2.0-70350217781 Kurata N. Kobori T. Takahashi M. Niwa N. Midorikawa H. Actual seismic response controlled building with semi-active damper system Earthquake Engineering and Structural Dynamics 1999 28 11 1427 1447 2-s2.0-0033230847 Rubió-Massegú J. Palacios-Quiñonero F. Rossell J. M. Decentralized static output-feedback H controller design for buildings under seismic excitation Earthquake Engineering & Structural Dynamics 2012 41 7 1199 1205 10.1002/eqe.1167 2-s2.0-84859954404 Gawronski W. Advanced Structural Dynamics and Active Control of Structures 2004 New York, NY, USA Springer Dyke S. J. Spencer B. F. Jr. Quast P. Kaspari D. C. Jr. Sain M. K. Implementation of an active mass driver using acceleration feedback control Microcomputers in Civil Engineering 1996 11 5 305 323 10.1111/j.1467-8667.1996.tb00445.x 2-s2.0-0030242672