MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 762385 10.1155/2013/762385 762385 Research Article Fault Tolerant Control for Civil Structures Based on LMI Approach http://orcid.org/0000-0002-7569-0945 Qu Chunxu Huo Linsheng Li Hongnan Zheng Bo-Chao Faculty of Infrastructure Engineering Dalian University of Technology Dalian 116024 China dlut.edu.cn 2013 20 11 2013 2013 18 07 2013 22 09 2013 16 10 2013 2013 Copyright © 2013 Chunxu Qu et al. 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.

The control system may lose the performance to suppress the structural vibration due to the faults in sensors or actuators. This paper designs the filter to perform the fault detection and isolation (FDI) and then reforms the control strategy to achieve the fault tolerant control (FTC). The dynamic equation of the structure with active mass damper (AMD) is first formulated. Then, an estimated system is built to transform the FDI filter design problem to the static gain optimization problem. The gain is designed to minimize the gap between the estimated system and the practical system, which can be calculated by linear matrix inequality (LMI) approach. The FDI filter is finally used to isolate the sensor faults and reform the FTC strategy. The efficiency of FDI and FTC is validated by the numerical simulation of a three-story structure with AMD system with the consideration of sensor faults. The results show that the proposed FDI filter can detect the sensor faults and FTC controller can effectively tolerate the faults and suppress the structural vibration.

1. Introduction

The structural control system is commonly classified into four types: passive, active, hybrid, and semiactive control . An active control system contains sensors, controllers, and actuators. Sensors can measure the structural vibration data that is sent to the controllers. The controllers compute the desired force according to the control strategy and sensor data and then command the control devices to generate control force for vibration reduction. Active control can achieve many control decisions. Therefore, active control is an efficient method to resist external excitation, such as earthquake and strong wind. Control strategy is the key of active control, which is widely studied in recent decades . However, there are some controller design problems in civil engineering due to the characteristics of the civil structures, such as large number of degrees of freedom and uncertainties of system parameters, in which many researchers are interested .

In practical engineering, sensors and actuators may be destroyed due to the lack of maintenance or extreme loading such as strong wind or intense earthquake. These faulty sensors and control devices may degrade the performance of vibration suppressing and even enlarge the structural response. Therefore, the robustness is important for the control system. Fault tolerant control (FTC) technology is an efficient method to promote the system reliability. FTC can tolerate faults caused by the components of control system, such as sensors and actuators, and then guarantee the system performance. FTC can be divided into two categories, passive FTC and active FTC. The difference between them exists if the faults are identified online. Passive FTC cannot detect faults online and cannot be recognized as traditional robust control . The faults in passive FTC are assumed as system model uncertainties. Accordingly, passive FTC has limited tolerance ability because its robustness is only for the estimated faults, and the faults cannot be detected online. On the other hand, active FTC can detect faults online, which is called the fault detection and isolation (FDI) technology, and suppress excessive structural response.

This paper designs a fault detection and isolation (FDI) filter and reforms a fault tolerant controller (FTC) through LMI approach considering sensor faults in civil building structures with active mass damper (AMD). An estimated system is built to facilitate the design of FDI and FTC. In order to make the difference between the estimated and the practical system smaller, an FDI filter is designed, which can minimize the H∞ norm of the transfer function from all the inputs to the residual of the estimated and practical state. According to the bounded real lemma, solving the FDI filter problem is transformed to solve a matrix inequality, which can be deformed to a linear matrix inequality (LMI) by variable substitution method . LMI problem is convex and can be solved easily. After the FDI filtering, the new data is used to isolate the sensor faults and transmitted to the previous designed controllers so that the control strategy is reformed. The process executes the function of fault tolerant control (FTC). Finally, the numerical example of a three-story building is provided to validate the performance of FDI and FTC.

2. Formulation of AMD Control System

Considering n degrees of freedom civil structural system which is subjected to a unidirectional earthquake with q¨g(t) acceleration, the dynamic equation is described as follows: (1)Mq¨(t)+Cq˙(t)+Kq(t)=Tgq¨g(t)+Tuu(t), where M, C, and Kn×n are, respectively, mass, damping, and stiffness matrix; q(t), q˙(t) and q¨(t)n×1 are the relative displacement, relative velocity, and relative acceleration vector, respectively; u(t)mu×1 is the control force vector generated by actuators and Tun×mu and Tgn×1 are the location matrices of the control force and earthquake excitation, respectively.

For a structure with an AMD on the top floor as shown in Figure 1, the control force is exerted only on the top floor. The control force location matrix Tu and earthquake excitation location matrix Tg are represented as (2)Tu=[0,0,0,n-11]1×nT,Tg=-M{1}n×1.

The model of AMD structure system.

The dynamic equation of the AMD system can be formulated as (3)ma(q¨a(t)+q¨g(t))+ca(q˙a(t)-q˙n(t))+ka(qa(t)-qn(t))=Fa(t), where ma, ca, and ka1×1 are, respectively, mass, damping, and stiffness of AMD system and qa(t), q˙a(t), and q¨a(t)1×1 are the relative displacement, relative velocity, and relative acceleration of AMD system. qn(t) and q˙n(t)1×1 are relative displacement and relative velocity of the nth floor, which is also the top floor and can be represented as (4)qn(t)=Tnq(t),q˙n(t)=Tnq˙(t),Tn=[0,0,0,n-11]1×n.Fa(t)1×1 is the force between structural system and AMD mass. The motion equation of the structural system with AMD can be derived as (5)Msq¨s(t)+Csq˙s(t)+Ksqs(t)=Tgsq¨g(t)+TusFa(t).

The matrices Ms, Cs, Ks, qs, Tgs, and Tu are denoted as (6)Ms=[Mma]=[m1mnma](n+1)×(n+1),qs=[qqa]=[q1qnqa](n+1)×1,Tgs=[Tg-ma]=-Ms{111}(n+1)×1,Tu=[-Tu1]={00-11}(n+1)×1,Ks=[K+TukaTn-Tuka-kaTnka]=[k1+k2-k200-k2-kn00-knkn+ka-ka00-kaka](n+1)×(n+1),Cs=[C+TucaTn-Tuca-caTnca]=[c1+c2-c200-c2-cn00-cncn+ca-ca00-caca](n+1)×(n+1), where mi, ci, and ki are the ith floor mass, damping, and stiffness. Equation (5) can be represented in state space as follows: (7)X˙s(t)=AsXs(t)+B1sq¨g(t)+B2sFa(t), where As2(n+1)×2(n+1),B1s2(n+1)×1, and B2s2(n+1)×1 are system matrix, excitation matrix, and control matrix, respectively.

One can see that(8)As=[(n+1)×(n+1)[I](n+1)×(n+1)-Ms-1Ks-Ms-1Cs],B1s=[(n+1)×1Ms-1Tgs],B2s=[(n+1)×1Ms-1Tus].Xs(t)=[qs(t);q˙s(t)]2(n+1)×1 is the state vector, which contains relative displacement and relative velocity. This relative state vector can be transformed to another state vector which consists of interstory drift and interstory velocity as,(9)X(t)=[q1q2-q1qn-qn-1qa-qn,q˙1q˙2-q˙1q˙a-q˙n]T2(n+1)×1.

Define a transformation matrix Γ2(n+1)×2(n+1) to obtain X(t) such that X(t)=ΓXs(t). Substituting Xs(t)=Γ-1X(t) into (7) and left-multiplying this equation with Γ, the state equation expressed by X(t) becomes (10)X˙(t)=AX(t)+B1q¨g(t)+B2Fa(t), where A=ΓAsΓ-12(n+1)×2(n+1), and B1=ΓB1s2(n+1)×1, B2=ΓB2s2(n+1)×1, and Γ is shown as (11)Γ=[Γs(n+1)×(n+1)(n+1)×(n+1)Γs],Γs=[100-1100000-11](n+1)×(n+1).

3. Fault Detection and Isolation and Fault Tolerant Control

A simple and efficient way to design a controller for this structural control system with AMD is to design the H∞ state feedback gain G1×2(n+1) , such that (12)Fa(t)=GX(t).

However, when there are some faults in sensors or actuators, the performance of vibration suppressing may be degraded and even enlarge the structural response. So, it is necessary to detect and isolate faults such that the controller can tolerate the faults. This paper only considers the faults in sensors and designs a fault detection and isolation (FDI) filter. The designed FDI filter is intended to estimate structural system state, which can be used to reappear the structural sensor data with fewer faults. The difference between the measured output data and the estimated output data is the faults. The estimated output data can be used to isolate the faults transmitted to the designed H∞ state feedback gain G. So, the controller strategy is reformed and can perform the fault tolerant control performance.

For the faults in sensors, the sensor data including faults can be represented as (13)yv(t)=CyX(t)+v(t), where v(t)q×1 is the faults part in sensor data yv(t)q×1 and CyX(t) is the real part of sensor data, which is the system real output. The faults can be recognized as the additional data into the real sensor data. Combined with (10), the following equation can be derived:(14)X˙(t)=AX(t)+B1q¨g(t)+B2Fa(t),yv(t)=CyX(t)+v(t).

The faults estimator is shown in Figure 2 and its state space formulation can be described as (15)X^˙(t)=AX^(t)+B2Fa(t)+r(t),y^(t)=CyX^(t),r(t)=K(y^(t)-yv(t)), where X^(t)2(n+1)×1 is the estimated state for the structural state X(t), y^(t)q×1 is the estimated measurement output for the real part of yv(t)q×1, that is, CyX(t), r(t)2(n+1)×1 is the adjustment vector for estimated state vector, and K2(n+1)×q is the FDI filter to make the estimated state vector closed with structural system state vector. In this structural system with the sensor faults, the state is estimated so that the sensor data real part CyX(t) can be estimated as y^(t)=CyX^(t). The faults in sensor data yv(t)q×1 can be calculated by yv(t)-y^(t) and then can be isolated by replacing yv(t) with y^(t). This is the fault detection and isolation (FDI). According to (14) and (15), the residue equation is obtained as follows: (16)X~˙(t)=A~X~(t)+B~w(t),Z~(t)=C~X~(t)+D~w(t).A~, B~, C~, D~, X~(t), Z~(t), and w(t) are defined as(17a)A~=A+KCy2(n+1)×2(n+1),B~=[B1K]2(n+1)×(q+1),C~=I2(n+1)×2(n+1),D~=2(n+1)×(q+1)2(n+1)×(q+1),(17b)X~(t)=(X(t)-X^(t))2(n+1)×1,Z~(t)=X~(t)2(n+1)×1,w(t)=[q¨g(t);v(t)](q+1)×1,where Z~(t)=X~(t)2(n+1)×1 is the output of the residue equation and X~(t) is the state of the residue equation.

FDI diagram of AMD structure system.

In order to estimate the system state X(t) precisely, this paper designs the FDI filter K in (15) to minimize the H∞ norm of the transfer function from input w(t) to output Z~(t). This design can be recognized as H∞ optimization problem. It can be described as designing a static controller K to stabilize the system (16) such that (18)J=Hzw(s)=supω,ω(t)20(Z~(t)2w(t)2)<γ, where γ is a given scalar. According to the bounded real lemma , the following two statements are equivalent:

H∞ norm of the system in (16) is less than γ, and A~ is stable;

There exists a symmetric positive definite matrix Φ2(n+1)×2(n+1) such that the following matrix inequality holds: (19)[A~TΦ+ΦA~ΦB~C~TB~TΦ-γID~TC~D~-γI]<0.

Substituting (17b) into (19),(20)[(A+KCy)TΦ+Φ(A+KCy)Φ[B1K]2(n+1)×(q+1)I2(n+1)×2(n+1)[B1K](q+1)×2(n+1)TΦ-γI(q+1)×(q+1)(q+1)×2(n+1)I2(n+1)×2(n+1)2(n+1)×(q+1)-γI2(n+1)×2(n+1)]<0.

Because both Φ and K are unknown variables and multiplied together, this problem cannot be solved directly. This paper employs variable substitution method  to linearize this constraint.

Define a new variable Ξ=ΦK2(n+1)×q and take it into (20)(21)[ATΦ+ΦA+CyTΞT+ΞCy[ΦB1]2(n+1)×1Ξ2(n+1)×qI2(n+1)×2(n+1)[B1TΦ]1×2(n+1)-γI1×11×q1×2(n+1)Ξq×2(n+1)q×1-γIq×qq×2(n+1)I2(n+1)×2(n+1)2(n+1)×12(n+1)×q-γI2(n+1)×2(n+1)]<0.

Equation (21) is a linear matrix inequality with matrix variables Φ and Ξ. The optimal values of Φ and Ξ can be computed through the solver feasp and mincx in MATLAB LMI toolbox. Then, K can be obtained from K=Φ-1Ξ. So, the FDI filter in (15) is designed.

Through FDI filter, the estimated state X^(t) and the estimated sensor data y^(t)=CyX^(t) can be very closed to the system state X(t) and the sensor data real part CyX(t), respectively. The estimated state X^(t) can be recognized as the input of H∞ state feedback gain G. This reformed control system can be described in Figure 3. The reformed control strategy is different from the one in (12) (22)F-a(t)=GX^(t).

FTC diagram of AMD structure system.

As a summary, when the sensor data has faults, the control system cannot guarantee a good performance. Through FDI filter, the sensor data real part and system state can be estimated. The estimated information can be used to isolate the sensor faults. Then, put the estimated state into the previous designed H∞ state feedback controller to reform the control strategy such that the control system can perform the fault tolerant control (FTC).

4. Numerical Example

Considering a three-story numerical example , an AMD control device is set on the top of this structure as shown in Figure 4. Define the structural mass, damping, and stiffness as follows: (23)M=×103kg,C=[12.4-5.16-5.1612.4-4.59-4.597.2]×103N/(m/s),K=[3.4-1.8-1.83.4-1.6-1.61.6]×106N/m.

Three-story structure model with AMD system.

AMD control system mass, damping, and stiffness are as follows: (24)ma=360kg,ka=18819N/m,ca=365.3910N/(m/s).

Suppose the sensor measurement outputs are the system states. H∞ state feedback control strategy is employed to control this structural system with AMD. The H∞ state feedback gain G in (12) can be obtained as follows:(25)G=[3.21553.12756.7705-0.0127-0.1344-0.0877-0.1335-0.0188]×106.

The El Centro NS earthquake acceleration record q¨g(t) with its peak scaled to 200gal (2m/s2) is used as the ground excitation. The response suppression comparisons are shown in Figures 5 and 6 for the bare structure without control, passive tuned mass damper (TMD) structure, and active mass damper (AMD) structure without sensor faults. The bare structure is the structure without actuators. The TMD structure has actuators but does not generate control force. The control strategy for the AMD structure is H∞ state feedback control.

Comparison of the 3rd floor relative displacement between no control and TMD passive control.

Comparison of the 3rd floor relative displacement between TMD passive control and AMD-H∞ state feedback control.

In Figure 6, G represents H∞ state feedback controller. Figures 5 and 6 illustrate the relative displacement response on the 3rd floor. The maximum value of the relative displacement for the TMD structure is 24.1% less than the one for the bare structure, while AMD structure is 67.0% off. Therefore, TMD and AMD can both suppress the structural vibration. Furthermore, AMD structure can reduce the vibration 56.5% more than TMD structure.

If there exist some faults in sensors, the faulty signal data v(t) is supposed to be located in the third floor displacement sensors. And its mean value and variance are zero and 2.5×10-5, respectively. The comparison of the 3rd interstory drift between the AMD structure and the TMD structure is shown in Figure 7.

Comparison of the 3rd story interstory drift between TMD and AMD with H∞ state feedback control strategy considering the sensor faults.

From Figure 7, AMD structural control system with H∞ state feedback control strategy loses its control performance when there exist faults in the sensors. Accordingly, it is necessary to detect and isolate the faults in the sensors. Through proposed FDI filter design method, the gain K is obtained as follows:(26)K=[-6.4162-2.4851-0.91600.4926-92.0527-118.6123-28.5338-23.9462-1.0675-2.3501-0.60250.1106-42.207593.6605-82.5947-27.40980.3623-0.1068-1.8490-0.409739.9383-55.9578105.6514-90.3240]T.

Taking K into (15) and putting the estimated state X^(t) into state feedback gain G, the control strategy is reformed to perform the fault tolerate control.

Figure 8 compares the interstory drift in the top story between the real interstory drift and the estimated drift. The related coefficients are shown in Table 1. From Table 1, the dash line, which is the real drift, matches the solid line, which is the estimated drift by FDI. From Table 1, the related coefficients (RCs) are almost close to 1, which means that FDI can estimate drift precisely.

Correlation coefficient between FDI estimated state and the real state.

Comparison between the states estimated by FDI and structure system real states
Interstory drift Interstory velocity
q 1 cm q 2 - q 1 cm q 3 - q 2 cm q a - q 3 cm q ˙ 1 cm q ˙ 2 - q ˙ 1 cm q ˙ 3 - q ˙ 2 cm q ˙ a - q 3 ˙ cm
RC 0.8422 0.8064 0.8584 0.9989 0.7741 0.7483 0.9035 0.9985

Comparison between the real and FDI estimated interstory drift in the top story.

Through FTC strategy, the interstory drift is suppressed, and the third story drift is shown in Figures 9 and 10. Table 2 shows some details about vibration suppression effects for other stories.

Suppressing effect comparisons among the bared structure, TMD structure, AMD structure without faults, and AMD H∞ structure with faults.

Systems q 1 cm q 2 - q 1 cm q 3 - q 2 cm q a - q 3 cm F a ( t ) kN
Bare system 4.04 2.81 1.68
TMD system 2.69 2.22 1.57 25.08
AMD system without faults using G 1.88 1.08 0.86 78.38 15.46
AMD system with faults using G 1.90 1.37 2.30 108.7 209.1
AMD system with faults using FTC 2.06 1.35 0.71 83.36 16.00

G represents H∞ state feedback control.

Comparison of the interstory drifts in the 3rd story between the TMD and AMD-FTC control.

Comparison of the 3rd floor interstory drifts between the AMD-H∞ state feedback control and AMD-FTC control.

Figure 9 illustrates the comparison of the interstory drift in the 3rd story between TMD structure and AMD structure considering sensor faults. In AMD structure, the control strategy is FTC. AMD structure can suppress the structural vibration more effectively than TMD structure does. While only using H∞ state feedback control strategy, the vibration is worsened as shown in Figure 7. Figure 10 can illustrate the difference more directly.

From Table 2, the maximum values of all the interstory drifts for AMD structure with H∞ state feedback control strategy are 29.37%, 38.29%, and −46.50% less than the TMD structural drifts, respectively. While using FTC strategy, the drifts are 23.42%, 39.19%, and 54.78% less than the TMD structural drifts, respectively. The actuator force Fa(t) saves a lot of energy using FTC strategy than when only using H∞ state feedback control strategy. Therefore, FTC can tolerate faults in sensors and suppress structural vibration efficiently.

5. Conclusions

In practical engineering, structural control system may have faults which can degrade control system performance. This paper presents fault tolerant control (FTC) for civil structures. Considering the faults in sensors, fault detection and isolation (FDI) filter is designed. Transmitting the estimated state to H∞ state feedback controller designed previously, the reformed control strategy is obtained, which is called fault tolerant control (FTC). This paper also presents the derivation of FDI, which transforms the inequality matrix to linear inequality matrix using variable substitution method. A three-story numerical example validates FDI and FTC performance. From this study, there are some conclusions.

When there is no fault in sensors, H∞ state feedback control in AMD structure has better suppressing performance than TMD structure.

When there are faults in sensors, H∞ state feedback control in AMD structure cannot guarantee the control performance and even worsen the vibration.

FDI designed by the proposed method can detect and isolate the faults. It can estimate the system state precisely. According to the estimated state, the measurement output from sensors can also be estimated.

FTC designed by the proposed method can tolerate the faults in sensors. So, it can strengthen the robustness of the control system.

Acknowledgment

The authors are grateful for the support of the National Natural Science Foundation of China (Grant nos. 51121005 and 51261120375).

Soong T. T. Active Structure Control Theory and Practice 1990 New York, NY, USA Longman Scientific & Technical Soong T. T. Dargush G. F. Passive Energy Dissipation Systems in Structural Engineering 1997 London, UK John Wiley & Sons Yao J. T. P. Concept of structural control Journal of Structural Division 1972 98 1567 1574 Song G. Qiao P. Sethi V. Prasad A. Active vibration control of a smart pultruded fiber-reinforced polymer I-beam Smart Materials and Structures 2004 13 4 819 827 2-s2.0-4043095599 10.1088/0964-1726/13/4/020 Song G. Schmidt S. P. Agrawal B. N. Experimental study of active vibration suppression of flexible structure using modular control patch 1 Proceedings of the IEEE Aerospace Conference March 1998 189 201 2-s2.0-0031629554 Song G. Kelly B. Agrawal B. N. Active position control of a shape memory alloy wire actuated composite beam Smart Materials and Structures 2000 9 5 711 716 2-s2.0-0034296026 10.1088/0964-1726/9/5/316 Sethi V. Franchek M. A. Song G. Active multimodal vibration suppression of a flexible structure with piezoceramic sensor and actuator by using loop shaping Journal of Vibration and Control 2011 17 13 1994 2006 2-s2.0-80155171959 10.1177/1077546310393440 ZBL1271.74349 Sethi V. Song G. Multimodal vibration control of a flexible structure using piezoceramic sensor and actuator Journal of Intelligent Material Systems and Structures 2008 19 5 573 582 2-s2.0-42449085433 10.1177/1045389X07077853 Gu H. Song G. Active vibration suppression of a flexible beam with piezoceramic patches using robust model reference control Smart Materials and Structures 2007 16 4 1453 1459 2-s2.0-34547456517 10.1088/0964-1726/16/4/060 Huo L. Li H. Multi-Dimensional Structure Damping Control 2008 Beijing, China Science Press Song G. Kelly B. Agrawal B. N. Lam P. C. Srivatsan T. S. Application of shape memory alloy wire actuator for precision position control of a composite beam Journal of Materials Engineering and Performance 2000 9 3 330 333 2-s2.0-0033699838 10.1361/105994900770346006 Wang Y. Time-delayed dynamic output feedback H controller design for civil structures: a decentralized approach through homotopic transformation Structural Control and Health Monitoring 2011 18 2 121 139 2-s2.0-79952566478 10.1002/stc.344 Bai Y. Grigoriadis K. M. H collocated control of structural systems: an analytical bound approach Journal of Guidance, Control, and Dynamics 2005 28 5 850 853 2-s2.0-26244444533 10.2514/1.12568 Gu H. Song G. Malki H. Chattering-free fuzzy adaptive robust sliding-mode vibration control of a smart flexible beam Smart Materials and Structures 2008 17 3 2-s2.0-45749137257 10.1088/0964-1726/17/3/035007 035007 Bai Y. Grigoriadis K. M. Damping parameter design optimization in structural systems using an explicit H norm bound Journal of Sound and Vibration 2009 319 3–5 795 806 2-s2.0-55349097087 10.1016/j.jsv.2008.06.046 Gu H. Song G. Active vibration suppression of a composite I-beam using fuzzy positive position control Smart Materials and Structures 2005 14 4 540 547 2-s2.0-23444449586 10.1088/0964-1726/14/4/012 Bai Y. Grigoriadis K. M. Song G. Lindner D. K. Active fault tolerant control of a flexible beam 6523 Modeling, Signal Processing, and Control for Smart Structures . 2007 Proceedings of SPIE Song G. Cai L. Wang Y. Longman R. W. A sliding-mode based smooth adaptive robust controller for friction compensation International Journal of Robust and Nonlinear Control 1998 8 8 725 739 2-s2.0-0032119060 10.1002/(SICI)1099-1239(19980715)8:8&lt;725::AID-RNC345&gt;3.0.CO;2-D ZBL0913.93038 Song G. Ma N. Robust control of a shape memory alloy wire actuated flap Smart Materials and Structures 2007 16 6 N51 N57 2-s2.0-36448981411 10.1088/0964-1726/16/6/N02 Song G. Mukherjee R. A comparative study of conventional nonsmooth time-invariant and smooth time-varying robust compensators IEEE Transactions on Control Systems Technology 1998 6 4 571 576 2-s2.0-0032124044 10.1109/87.701355 Li L. Song G. Ou J. Adaptive fuzzy sliding mode based active vibration control of a smart beam with mass uncertainty Structural Control and Health Monitoring 2011 18 1 40 52 2-s2.0-79951786603 10.1002/stc.356 Stavroulakis G. E. Marinova D. G. Zacharenakis E. C. Robust control in structural dynamics, taking into account structural uncertainties Proceedings of the 5th GRACM International Congress on Computational Mechanics 2005 Limassol, Cyprus Yu L. Robust Control—Linear Matrix Inequalities Approach 2002 Beijing, China Tsinghua University Press Bai Y. Grigoriadis K. M. Sethi V. Static output feedback H control of collocated structural systems 1–12 Proceedings of the American Control Conference June 2006 828 833 2-s2.0-34047199289 Mackenroth U. Robust Control Systems: Theory and Case Studies 2004 Berlin, Germany Springer Wang Y. Lynch J. P. Law K. H. Decentralized H controller design for large-scale civil structures Earthquake Engineering and Structural Dynamics 2009 38 3 377 401 2-s2.0-66149115256 10.1002/eqe.862