LQR Control of Wind Excited Benchmark Building Using Variable Stiffness Tuned Mass Damper

LQR control ofwind inducedmotion of a benchmark building is considered.Thebuilding is fittedwith a semiactive variable stiffness tuned mass damper adapted from the literature. The nominal stiffness of the device corresponds to the fundamental frequency of the building and is included in the system matrix. This results in a linear time-invariant system, for which the desired control force is computed using LQR control. The control force thus computed is then realized by varying the device stiffness around its nominal value by using a simple control law. A nonlinear static analysis is performed in order to establish the range of linearity, in terms of the device (configuration) angle, for which the control law is valid. Results are obtained for the cases of zero and nonzero structural stiffness variation. The performance criteria evaluated show that the present method provides displacement control that is comparable with that of two existing controllers. The acceleration control, while not as good as that obtained with the existing active controller, is comparable or better than that obtained with the existing semiactive controller. By using substantially less power as well as control force, the present control yields comparable displacement control and reasonable acceleration control.


Introduction
Active control devices, such as the Active Tuned Mass Damper (ATMD), require substantial input power and could also destabilize the system if the controller is improperly designed.On the other hand, passive control devices are less effective in the presence of stochastic disturbances and/or structural property variations.Semiactive control devices do not possess these disadvantages and thus appear to be sound alternatives to active and passive devices [1][2][3][4].Such devices provide control forces by varying their mechanical properties, based on feedback.The variable stiffness damper is a semiactive device with good potential for controlling wind/earthquake generated response.Kobori et al. [5] and Nasu et al. [6] considered an active variable stiffness (AVS) system, comprising an on-off type two ended hydraulic damper, to make the structure nonresonant during an earthquake.Nemir et al. [7] considered a variable stiffness bracing and obtained rapid dissipation by way of energy redistribution to higher modes.Such AVS systems, while effective, cause abrupt switching of stiffness.Yang et al. [8] proposed a sliding mode controller for an AVS system.A resetting control algorithm, involving the release of potential energy of the device followed by a quick resetting of the device to its full-stiffness state, was considered by Yang et al. [9].Yang et al. [2] proposed a 76-storey building in Melbourne as a benchmark structure for evaluating algorithms for wind induced response control.Results using LQG control with an ATMD were obtained in their study.
Nagarajaiah [10] developed a semiactive variable stiffness (SAVS) device and studied its performance using a scaled model.The SAVS-TMD has been shown to be effective for structures that are subjected to force/base excitation [11].Varadarajan and Nagarajaiah [12,13] studied the wind response control of the benchmark building [2].They used Empirical Mode Decomposition-Hilbert Transform method [12] and Short Time Fourier Transform (STFT) method [13] in order to track the dominant response frequency.The SAVS-TMD was then tuned to this frequency.Wu and Yang [14] studied the performance of Linear Quadratic Gaussian control (LQG),  ∞ control, and continuous sliding mode control applied to an active mass driver, for acceleration  [13].
reduction of the wind excited Nanjing tower.Using a variable stiffness TMD, Collins et al. [15] considered bangbang control combined with semiactive control.However, they did not consider the actuator dynamics.Semiactive controller designs using other devices are also available.For magnetorheological devices, Yuen et al. [16] used reliability based robust linear control with a clipped control law, and Karamodin and Kazemi [17] used LQG control and a semiactive neural controller with acceleration/velocity feedback.Sohn et al. [18] studied the semiactive control of a suspension system, by estimating the road profile using the extended least squares method and then applying LQG control.Gaul et al. [19] studied the control of a truss with semiactive friction joints.They used two methods, that is, one with a controller for each joint and another with a single clipped-optimal controller.
In the present study, the SAVS-TMD of [12,13] is deployed in order to control the wind excited benchmark building [2] by using a Linear Quadratic Regulator (LQR) controller.The nominal stiffness of the device, corresponding to the fundamental frequency of the structure, is included in the system matrix.This results in a linear time invariant system, for which algorithms suitable for real-time control applications can be employed.One such algorithm is LQR control wherein gains are computed offline, thus making it suitable for real-time control.The desired control force is computed using LQR control and then realized by changing the stiffness of the device within limits specified around the nominal stiffness.This is done using a simple control law which changes the configuration of the device-within specified limits-by means of an electromechanical actuator.The dynamics of the actuator are excluded from this study.A nonlinear static analysis is performed in order to obtain the operational range of the device configuration angle.This ensures a linear force-displacement behavior for the device, and hence validity of the control law.In order to assess robustness of control, the controller thus designed is implemented on the structure having ±15% stiffness variation [2].The goals of the paper are (i) implementation of STLC, that is, the SAVS-TMD device with LQR control, with a simple control law that is valid within the operational range of the device, (ii) comparing the performance of STLC with that of the controller in [2], which is based on ATMD with LQG control (ALC), and the controller in [13] which is based on SAVS-TMD with STFT control (STSC).This is done for all cases of structural stiffness variation, in order to assess control robustness.In contrast to STLC where gains are computed offline (i.e., computed only once), STFT involves a time-varying system with online computations for real-time frequency tracking during its control law implementation, so as to tune the device to the tracked frequency.This involves intensive online computations which increase the controlloop time and thus renders STFT less suitable for real-time control.On the other hand, ALC, being an active method, requires substantially more power and control force than STLC.Thus, the present study provides a new power-efficient controller design for the benchmark problem, that is, one which is suitable for real-time control and which can be readily extended to output-based feedback control in order to further decrease the loop time.

Semiactive Variable Stiffness TMD
The SAVS-TMD is fitted at the top of the 76-storey benchmark building [2], as shown in Figure 1(a).The device comprises a rhombus of four springs, each having stiffness   and unstretched length   .The springs are pin connected at sliding joints (Figures 1(b) and 2).The masses of sliders and springs, and the effect of friction, are neglected.Joint-3 and joint-4 slide along a horizontal guide-rail fixed on the floor.Joint-2 slides along the -directed groove which is present at the bottom of the TMD mass.The TMD mass moves along the -directed rails that are fixed on the floor.By using a controlled actuator, joint-1 can be made to move along a -directed guide that is fixed on the floor.This causes the stiffness of the device to vary, due to variation in the device configuration angle .The coordinates of joints 1, 3, and 4 are denoted  1 ,  3 , and  4 , respectively.The  and  coordinates of joint-2 are denoted  drel and  2 , respectively, with  drel being the TMD displacement measured relative to the top storey (Figure 2).The joint coordinates are measured as per directions shown in Figure 2. The TMD displacement results in force   at joint-2, measured positive rightward.Equilibrium of joints 2, 3, and 4 yields  23 sin  2 +  24 sin  3 = 0;  23 cos  2 +  13 cos  1 = 0;  24 cos  3 +  14 cos  4 = 0; where are the spring forces.Here  1 , . . .,  4 are expressed in terms of joint coordinates.For example, cos  2 = ( 3 +  drel )/( √ ( 3 +  drel ) 2 +  2 2 ).For  1 held fixed and  drel = 0, the equilibrium configuration angles are equal and denoted as ; that is,  1 =  2 =  3 =  4 =  = sin −1 ( 1 /  ).Thus, , or equivalently  1 , represents the device configuration.Choosing device configuration  1 (which is equivalent to choosing ),    drel , by solving (1) using MATLAB fsolve.Thus, the forcedisplacement behavior of the SAVS-TMD, that is,   versus  drel , is obtained as shown in Figure 3. Here,   = 2m,   = 650 kN/m, and various device configurations , such that 5 ∘ ≤  ≤ 50 ∘ , are considered.The behavior is almost linear up to  = 25 ∘ , beyond which the linearity holds within a range of  drel .This range reduces as the device configuration  (or  1 ) increases.The nonlinearity becomes significant as the device configuration angle increases.The nonlinear behavior begins, for example, at a TMD displacement of around 0.55 m for  = 40 ∘ and at a TMD displacement of around 0.45 m for  = 45 ∘ .Thus, for large TMD-displacement, the device behaves nonlinearly when the device configuration angle is also large.Assuming that  drel ≪   , that is, the TMD displacement is small relative to the unstretched spring length, the model of [12,13] can be used.This model implies a linear forcedisplacement relation for a given configuration , with device stiffness given as Here,   () is the stiffness provided by the device for a given configuration .Based on the force-displacement behavior shown in Figure 3, the maximum configuration angle is restricted to  max = 41 ∘ (i.e., open position) in order to maintain device linearity.This corresponds to around 10% nonlinearity at  drel = 0.68 (Figure 3).Further, the minimum configuration angle is restricted to  min = 7 ∘ (i.e., closed position) due to mechanical constraints during flattening of the rhombus.The nominal configuration angle of the device is set as   = 28.2∘ in order that, in this configuration, the device is tuned to the fundamental frequency of the structure.Depending on the control force required, the device stiffness is increased/decreased by varying the device configuration, , about   , such that  min ≤  ≤  max .With the linear forcedisplacement relation of (3) used in the controller design, the RMS value of TMD displacement is obtained as less than 26 cm, as seen from the controlled responses in Table 1.When both the control-force required and the TMD-displacement are large and have opposite signs, the device angle  required is large and the force-displacement relation is nonlinear (Figure 3).However, instances of this happening are few, as is evident from the relatively small RMS values of  drel (given in Table 1) when compared to the linearity limit of 0.55 m for  = 40 ∘ .Hence, the model given by ( 3) is henceforth used in the control law.The stiffness of the device can also be written as Here,   is the nominal stiffness that corresponds to the TMD being tuned to the fundamental frequency of the structure.Further, Δ  () is the additional stiffness required to attain the desired control force.The additional stiffness is obtained by varying the device configuration, , such that

Reduced Order Model
The equation of motion of the wind excited building, with SAVS-TMD at the top storey, is written as Here, X = [ 1 ,  2 . . .,  −1 ,   ]  is the displacement vector measured relative to the ground, where   ,  = 1, . . ., ( − 1), is the lateral displacement of the th storey and   is the TMD displacement, () is the control input (i.e., force) that is realized via the additional stiffness provided by the SAVS device, M = diag[M  ,   ] is the system mass matrix, C = diag[C  , 0] + diag[0, C  ] is the system damping matrix, and is the system stiffness matrix, where subscript  refers to the structure and subscript  refers to the damper, ,   , and   , are the mass, damping coefficient, and nominal stiffness, respectively, of the TMD, D = [0, 0, . . ., 1, −1]  is the control force placement vector, and W() is the wind excitation vector, with its last element being zero.Since the system stiffness matrix contains only the nominal stiffness of the TMD, the system is timeinvariant and the control force required is given as The state space representation of ( 5) is where The output vector for assessing control effectiveness is given by where C y , D y , and F are as defined in Section 5.In order to reduce the computation time, for real-time control applications, the model is reduced by using the method of Davison [20] as done in [2,14].This entails choosing a reduced set of 2 states that are representative of the system response and then expressing this reduced set of states in terms of the first 2 eigenmodes that dominate the response.Let Γ represent the matrix of eigenvectors of A, with the eigenvectors arranged in decreasing order of eigenvalue dominance (i.e., the first eigenvector corresponds to the eigenvalue lying closest to the origin and the last eigenvector corresponds to the eigenvalue lying farthest from the origin).Rearrange is the reduced order 2-dimension state vector and X  is the reduced order dimension displacement vector.Here, Z  comprises the 2( − ) states, of the full order system, that are excluded from the reduced order system.The A, B, C y , E, Γ matrices are now rearranged according to Z = [Z   Z   ]  .Thus, the approximation of the states in terms of the first 2 eigenmodes yields [14 Z  (10) and it also yields the reduced order system

Controller Design
Linear quadratic regulator (LQR) control [21] is considered for controller design.The reduced order model, that is, (11), is considered for the plant dynamics.Thus, the reduced states are assumed to be available for feedback.Using LQR control, the control force   is obtained such that the performance index is minimized.The performance index represents the total energy of the system (i.e., the structure and the SAVS-TMD).
Here,  is the positive scalar weighting of the control effort, and Q is the (2 × 2) positive semidefinite state weighting matrix; that is, Here,  1 ,  2 , and , are chosen such that effective control is achieved without exceeding the limits prescribed on the TMD displacement,  drel , and the device configuration angle .Equations ( 10) and ( 14) yield the (2 × 2) positive semidefinite state weighting matrix for the reduced order system as Here, Q has been rearranged according to Z = [Z   Z   ]  .Minimizing  subject to the constraint represented by (11) considered without wind excitation, one obtains the desired optimal control force as [21] Here, P is the solution of the algebraic Riccati equation given as Considering ( 4), (6), and the configuration limits of the SAVS device, one obtains the control law.This yields the position of joint-1, that is required to realize the desired control force   , as Equations ( 3) and (4) yield the control force thus realized; that is, An alternative procedure to determine the device configuration  1 that involves the solution of the nonlinear static equations ( 1) is as follows.For a known desired control force,   , and TMD displacement,  drel , the total force required from the TMD is obtained as   =    drel +   .Using   and  drel in (1), the device actuation,  1 , as well as  2 ,  3 , and  4 , can be obtained, subject to the minimum values permissible for  1 and  4 .In this manner, the linear forcedisplacement relation, (3), and the resulting control law, (19), are not used.However, as discussed in Section 2 (on the basis of RMS values of the TMD displacement), the linear forcedisplacement relation suffices for the present application.Hence, this alternative procedure is not adopted.

Results
The 76-storey benchmark building is modeled using 76 translational degrees of freedom, as considered in [2].The TMD, having mass   = 500, is placed at the top storey.The mass matrix, M  , stiffness matrix, K  , and damping matrix C  for the structure, as well as the across-wind data, are considered from [22].A ±15% variation in structural stiffness is also considered in order to assess the effectiveness of the controller [2].The damping ratio for the TMD is considered as 7% [12].As one of the aims of this study is to compare results from STLC with those using ALC [2] and STSC [13], the reduced order model of [2] is considered.This is a 24-degreeof-freedom model, with the reduced state, X  , comprising the displacements at storeys 3, 6, 10, 13, 16, 20, 23, 26, 30, 33, 36, 40, 43, 46, 50, 53, 56, 60, 63, 66, 70, 73, and 76, and the TMD displacement, all measured relative to ground.The wind force vector is obtained by lumping wind forces at the reduced DOFs, with the wind location matrix E  modified appropriately [2].Thus, the third part of ( 12) is not considered when obtaining the wind force vector.Figure 4 shows the resulting across-wind load that acts on storeys 50 and 73 for the reduced order model.
The nominal stiffness of the SAVS device,   , is tuned to the fundamental frequency of the structure.The fundamental frequency is 0.16 Hz, resulting in   = 505.12kN/m.Using (3), the stiffness of the spring used in the device is chosen as   = 650 kN.This ensures that the nominal stiffness of the device is the average of the stiffness values of the device at its operational limits  min = 7 ∘ and  max = 41 ∘ .Thus, the limits on the additional stiffness that can be provided by the device are Δ  max = 135.223and Δ  min = −134.89kN/m.
The parameters for the controller design are chosen as  1 = 900,  2 = 1, and  = 2 × 10 −2 .The reduced order plant, that is, (11), is integrated using MATLAB ode-45.The initial conditions for the state and the initial control input are considered as zero (i.e., Z  = 0 and  = 0).Thus, Z  is obtained at the end of each time step.Subsequently, the desired control force,   , is obtained from (17), the position of joint-1,  1 , is obtained from the control law, that is, (19), and the control force, , applied at the start of the next time step is obtained from (20). Figure 5 shows the block diagram for the control loop simulation.
Equation ( 20) Equation ( 19) Equation ( 11) Performance criteria denoted by  1 , . . .,  12 are evaluated [2].These are defined in terms of controlled responses that are suitably normalized wherever indicated.The uncontrolled structure with zero variation in stiffness is considered when obtaining the normalizing quantity which, unless noted otherwise, corresponds to the response being normalized.The displacement and acceleration of storeys 1, 30, 50, 55, 60, 65, 70, 75, and 76 and the displacement and velocity of the TMD measured relative to storey 76 are considered.Hence, the matrices C y , D y , and F, appearing in the output equation, that is, (13)  1 denotes the maximum RMS acceleration, normalized with the RMS acceleration of storey 75, with storey 76 and the TMD being excluded from the maximum. 2 denotes the average of the normalized RMS accelerations, with the average being taken over storeys 50 to 75. 3 denotes the normalized RMS displacement of storey 76.  4 denotes the average of the normalized RMS displacements, with the average being taken over storeys 50 to 76.  5 denotes the RMS of  drel (i.e., the TMD displacement relative to storey 76) normalized with the RMS displacement of storey 76;  6 denotes the RMS of  ẋ drel , that is, the average input power.Note that in [13] the "average power" is defined as the RMS of ẋ drel normalized with the RMS velocity of storey 76, that is, a ratio of velocities.This is denoted here as Ĵ6 and is used for comparison with the results from STSC.This essentially represents the RMS of the TMD velocity.Performance criteria  7 to  12 and Ĵ12 are defined analogous to  1 to  6 and Ĵ6 , respectively, by replacing the RMS values with corresponding peak values.The constraints stipulated in the benchmark problem are   ≤ 100 kN for the RMS control force,   drel ≤ 30 cm for the RMS TMD-displacement, max || ≤ 300 kN for the peak control force, and max | drel | ≤ 95 cm for the peak TMD-displacement.
The present STLC results are compared with the ALC results and the passive TMD results obtained by Yang et al. [2,22] for 20% damping ratio and with the STSC results obtained  by Nagarajaiah and Varadarajan [13] for 7% damping ratio.Figure 6 compares the time histories of displacement and acceleration of storey75.The case of 0% stiffness variation is considered.The present STLC yields displacements that are comparable with those from ALC and STSC.The acceleration control obtained from STLC is marginally better than that from STSC, but it is somewhat inferior when compared to ALC results.The passive TMD is the least effective of the three controllers.
Peak responses from the four controllers, for the case of 0% stiffness variation, are compared in Figure 7.The passive TMD reduces peak displacements by 17-21% as compared to the uncontrolled structure.The STLC, STSC, and ALC methods yield a further reduction of 7-10% in peak displacements, as compared to the passive-TMD controlled structure.Peak displacements from STLC, STSC, and ALC are comparable, with the present STLC being marginally lower than ALC and STSC.The present STLC yields a reduction of 21-32% in peak accelerations when compared to passive TMD control.It yields marginally lower peak accelerations compared to STSC.However, STLC mostly yields an increase of up to 10% in peak accelerations when compared to ALC.Thus, for 0% stiffness variation, STLC provides storeywise peak responses that are comparable with results of STSC and ALC, except in respect of accelerations for which STLC is inferior as compared to ALC.Note that ALC is generally better than STLC or STSC in acceleration control, since it is an active method whereby the desired control force can always be attained (up to a 300 kN limit) by the actuator and is independent of the ATMD displacement.In contrast, the force attained by the SAVS-TMD at a particular instant is limited by the range of device stiffness (i.e., 505 ± 135 kN/m) and the TMD displacement at that instant.RMS responses from the four controllers, for the case of 0% stiffness variation, are compared in Figure 8.The passive TMD yields up to 32% reduction in displacements and up to 42% reduction in accelerations, as compared to the uncontrolled structure.The present STLC attenuates displacements by 15% and accelerations by 21-26%, as compared to passive TMD control.It yields a marginal 3% reduction in displacements and up to 7% reduction in acceleration, as compared to STSC.It is comparable in displacements but yields accelerations that are mostly higher, by up to 9%, when compared to ALC.Thus, for 0% stiffness variation, STLC provides comparable-to-moderately-better storeywise RMS responses as compared to STSC and ALC, except for RMS accelerations for which it is inferior as compared to ALC.
Peak and RMS responses from the four controllers, for the case of +15% stiffness variation, are compared in Figures 9 and 10, respectively.When compared to passive-TMD control, the present STLC yields an increase of up to 4% in peak displacements, a reduction of 8-14% in peak accelerations, a reduction of around 5% in RMS displacements, and a reduction of 4-10% in RMS accelerations.When compared to STSC, STLC yields an increase of up to 5% in peak displacements and 9% in peak accelerations.When compared to ALC, STLC yields an increase of up to 6% in peak displacements, 8-25% in peak accelerations, up to 3% in RMS displacements, and 6-22% in RMS accelerations.Thus, for +15% stiffness variation, STLC provides mostly inferior storeywise peak and RMS responses as compared to STSC and ALC.
Peak and RMS responses from the four controllers, for the case of −15% stiffness variation, are compared in Figures 11 and 12, respectively.When compared to passive-TMD control, the present STLC yields a reduction of around 19% in peak displacements, 14-25% in peak accelerations, 13% in RMS displacements, and 18-24% in RMS accelerations.When compared to STSC, STLC yields a reduction of around 4% in peak displacements and up to 8% in peak accelerations.When compared to ALC, STLC yields an increase of around 5% in peak displacements, but a marginal reduction in RMS displacements.The peak and RMS accelerations show a mixed trend, that is, the performance of STLC ranges from a reduction of 5% to an increase of 12% in peak accelerations, and a reduction of 12% to an increase of 20% in RMS accelerations.Thus, for −15% stiffness variation, STLC provides storeywise displacement responses comparable to results from STSC and ALC and moderately better storeywise peak accelerations as compared to STSC and a mixed trend in storeywise peak and RMS accelerations vis-a-vis ALC.Performance criteria  1 - 12 obtained using the present STLC are compared with those from ALC (shown in parentheses) and STSC (shown in square brackets) in Table 1.STLC yields RMS and peak displacements, that is,  3 ,  4 ,  9 ,  10 , that are comparable with results from ALC and STSC.The average RMS acceleration,  2 , obtained from STLC is comparable with results from ALC for the 0% and −15% stiffness variation cases, 10% higher than the ALC results for the 15% stiffness variation case, comparable with results from STSC for the 0% and 15% stiffness variation cases, and 9% lower than the STSC results for the −15% stiffness variation case.The maximum RMS acceleration,  1 , obtained from STLC is 18%, 22%, and 16% higher than the ALC results for the 0%, +15%, and −15% stiffness variation cases, respectively, comparable with results from STSC for the 15% stiffness variation case, and 7% and 10% lower than the STSC results for the 0% and −15% stiffness variations, respectively.The average peak acceleration,  8 , obtained from STLC is comparable with results from ALC for the 0% and −15% stiffness variation cases, 12% higher than the ALC results for the 15% stiffness variation case, comparable with results from STSC for all three stiffness variation cases.The maximum peak acceleration,  7 , obtained from STLC is 16%, 24%, and 7% higher than the ALC results for the 0%, +15%, and −15% stiffness variation cases, respectively, comparable with results from STSC for the 0% and 15% stiffness variation cases and 8% lower than the STSC results for the −15% stiffness variation case.
Thus, STLC provides displacement control that is comparable to that provided by ALC/STSC.However, the acceleration control obtained from STLC is not as good as that obtained from ALC; that is, the average-peak and average-RMS accelerations for the +15% stiffness variation case and the max-peak and max-RMS accelerations for all three stiffness variation cases show lesser attenuation when using STLC as compared to ALC.The reason for this, as discussed previously, is that ALC is an active control method that yields superior acceleration control.The acceleration control obtained from STLC is better than that obtained from STSC for the −15% stiffness variation case, and otherwise comparable with the STSC results.
The RMS relative displacement of the SAVS-TMD, that is,   drel , obtained from STLC is 6% higher than the ALC result for the +15% stiffness variation case, 8% lower than the ALC result for the −15% stiffness variation case, comparable with results from ALC for the 0% stiffness variation case, and 9% higher than the STSC results for the 15% stiffness variation case, comparable with results from STSC for the 0% and −15% stiffness variation cases.Note that, for ALC, the actuator stroke is considered instead of the SAVS-TMD displacement.The RMS displacement of storey 76 is 10.14 cm for the 0% stiffness variation case.Since  5 is obtained by merely normalizing   drel with this value, it is omitted from Table 1.The peak relative displacement of the SAVS-TMD, that is, max | drel |, is the highest when using STLC, that is, higher by as much as 18% and 22% vis-a-vis ALC and STSC, respectively.The peak displacement of storey 76 is 32.30cm for the 0% stiffness variation case.Since  11 is obtained by merely normalizing max | drel | with this value, it is omitted from Table 1.
The average power expended,  6 , and the peak power expended,  12 , when using STLC are at least 49% lower visa-vis ALC, since the latter is a fully active control method.The RMS control force,   , is at least 32% lower, and the peak control force, max |()|, is at least 21% lower, when using STLC vis-a-vis ALC.Thus, the semiactive STLC requires substantially less power and control force in order to achieve comparable displacement control and acceptable acceleration control vis-a-vis ALC.Note that the RMS and the peak values of both the SAVS-TMD displacement and the control force provided by the SAVS-TMD are within limits prescribed by the benchmark problem.Comparisons between STLC and STSC-for power and control force-are excluded, since the values of  6 ,  12 ,   , and max |()| are not available for STSC.The index Ĵ6 obtained from STLC varies between 6% lower and 7% higher than the ALC results, and it is up to 6% higher than the STSC results.The index Ĵ12 obtained from STLC is up to 11% higher than the ALC results and up to 8% higher than the STSC results.However, as noted previously, Ĵ6 and Ĵ12 represent the RMS and peak values, respectively, of the TMD velocity; that is, they do not actually represent the power input.
Figure 13 shows the time history of the applied and the desired control forces for the 0% stiffness variation case.The applied force depends on the relative displacement and the stiffness of the SAVS-TMD.Since the device stiffness available is limited, due to the constraints imposed on the configuration angle, the semiactive device is not always able to produce the desired force.The RMS value of the difference between applied and desired forces is 38.6 kN. Figure 14 shows the time history of the device position, that is,  1 (Figure 2).The horizontal portions represent the limits on the device configuration (19).The device position, and hence its stiffness, varies smoothly except when there is a change in  the direction of its relative displacement  drel (19).When the device position changes abruptly, the applied control force is very small due to  drel being small at these instants.However, the applied control force does not change abruptly at these instants.

Conclusion
LQR control is employed for a wind excited benchmark building that is fitted with a semiactive variable stiffness TMD device.The nominal stiffness of the device is tuned to the fundamental frequency of the structure and included in the system stiffness matrix.The additional, time-varying, component of the device stiffness is obtained via LQR control and a suitable control law and utilized to apply the control force.A nonlinear static analysis is done to establish the operational limits on the device configuration angle, so as to ensure near-linear behavior of the device and thus yield a simple control law.
Comparison of the present STLC with ALC and STSC permits the following conclusions.
(i) The performance criteria show that STLC generally provides displacement control that is comparable with that of ALC and STSC.However, the acceleration control from STLC is not as good when compared with ALC results.The acceleration control from STLC is better than that from STSC for the −15% stiffness variation case, and otherwise comparable with STSC results.
(ii) For the 0% stiffness variation case, STLC provides comparable-to-moderately-better peak and RMS storeywise responses, except that, for accelerations, it is inferior vis-a-vis ALC (since the latter is Shock and Vibration a fully active method).For the +15% stiffness variation case, STLC generally provides inferior peak and RMS storeywise responses.For the −15% stiffness variation case, STLC provides storeywise displacement responses comparable with ALC and STSC results, moderately better peak storeywise accelerations vis-a-vis STSC, and a mixed trend in peak and RMS storeywise accelerations vis-a-vis ALC.
(iii) The STLC requires substantially less power and control force to achieve comparable displacement control and acceptable acceleration control vis-a-vis ALC.
Future work would involve LQG and Static Output Feedback controller designs that reduce the number of sensors required.Further, knowing the SAVS-TMD displacement (i.e.,  drel ) and the desired control force (i.e.,   , obtained by linear control methods), one can obtain the device actuation required (i.e.,  1 ) via the nonlinear static equation (1) or their dynamic counterparts (in case the masses of the sliders are substantial).This would eliminate the limiting of the device configuration angle and would thus reduce the error between desired and applied control forces.

Figure 4 :
Figure 4: Time history of wind force on the 50th and 73rd storeys [22].

Figure 13 :
Figure 13: Time history of desired and applied control force (0% stiffness variation).