Enhancing the Classical Closed-Loop Algorithm in terms of Power Consumption

An approach is suggested to reduce the peak and average control forces of actively controlled structures. In this method, responses of an actively controlled building should be much smaller than the responses of the same building controlled by the similar passive control mechanism.This approach leads to a time varying gain matrix, which is not restricted by external excitation but it is related to the selection of a scalar function. Extensive numerical analyses by using various scalar functions show that the proposed strategy effectively can reduce the need of the required control force consumptions.


Introduction
Significant efforts have been directed toward the possibility of employing protective systems in civil engineering structures including passive, semiactive, active, and hybrid control systems [1][2][3][4][5].Among these, active control systems are the most advanced technique to achieve this issue.It makes a control mechanism act in a manner that automatically supplies the forces to control the responses of a structure depending upon the severity of the external forces and the instant state of the structure during vibrations.The direction and magnitude of the controlling forces in an active system are estimated by a competent algorithm.General concept of the modern control methods to reduce structural vibrations due to external excitations was proposed in the 1960s [6].Thereafter, a vast amount of efforts has been devoted to develop a variety of control algorithms based on different control design criteria, like the classical optimal control method [7], the pole assignment method [8], the predictive control method [9], and the instantaneous optimal control method [10,11].Although the linear quadratic regulator algorithm, LQR [6], is known as the one of the most famous and historic methods with widespread applications in different fields (because of its simple procedure and ease of implementation on actual large-scale systems), only the classical closed-loop methods control is feasible for structural control applications.This means that there are some other methods which perform better than the closed-loop control algorithm; however, since they need foreknowledge of entire external disturbance time history, they are not applicable.On the other hand, since the Riccati equation is obtained by ignoring the earthquake excitation term, classical closed-loop control algorithms are approximately optimal methods and do not entirely satisfy the optimality conditions [12].To overcome the shortcomings of the classical optimal control algorithms, several methods have been proposed until now such as the instantaneous optimal control method [10,11], the generalized optimal active control (GOAC) [13], and the instantaneous optimal Wilson- method [11].These algorithms, in spite of their suitable performances, are sensitive to change of time increment.Also, they could not significantly improve the performance of the controlled structure with respect to the classical optimal control methods.Basu and Nagarajaiah proposed a wavelet-based adaptive linear quadratic regulator formulation for optimal control problems [14].They also presented a method for the control of time varying systems 2 Advances in Civil Engineering based on wavelet transformation [15].Bagheri and Amini proposed a procedure based on pattern search method and wavelet to acquire a more efficient control scheme than the LQR [16].
In this paper, an approach called new linear quadratic regulator (NLQR) is proposed, in which a relationship is established between the active and passive structural responses, via a proper transformation matrix.Then, this relationship is linked to a scalar function (considering an additional assumption) which causes control force computations to be quite simple.In order to evaluate effectiveness and performance of the proposed method, it is applied to an 8-story shear type building structure to reduce its seismic vibration responses during the El-Centro earthquake excitation.A few proper scalar functions are derived using the El-Centro earthquake ground acceleration and then the performance of this building is evaluated during some other earthquake excitations.

Classical Optimal Closed-Loop Algorithm (CCLQR)
Consider a building equipped with an active control system excited by a strong ground motion.The governing dynamic equation of motion may be written as the following matrix form: where  is the -dimensional displacement vector and dots state the derivative of  with respect to time as the velocity and acceleration vectors; , , and  are the  ×  mass, damping, and stiffness matrices of the structure, respectively,  is the  × 1 influence vector of the ground acceleration on the building masses,  is the  ×  location matrix of the control forces affecting the structure, and () is the  × 1 control force vector applied by the  actuators.With some manipulations, the equation of motion may be rewritten in terms of the state-space variables  as follows: in which  0 is the initial time instant, () is the vector of state variables, and  depicts the system matrix, respectively.Vector  and matrix  are defined as follows: In addition, matrix  and vector  are given as In the classical linear optimal control, a performance index () is defined in order to minimize control forces to achieve the best structural performances as where  is a 2 × 2 positive semidefinite weighting matrix related to the structural response,  is an × positive definite weighting matrix related to the active control force, and   indicates the terminal time that should be longer than the earthquake duration.
In classical closed-loop quadratic regulator (CCLQR), optimal control force is achieved by minimizing a performance index () with respect to the vectors () and () as follows: where () is the Riccati matrix, obtained by solving the following nonlinear matrix equation: There are two assumptions to solve (7); first, the external disturbance () is set to zero (it is assumed as a stochastic white noise process); second, the Riccati matrix is a constant one over the time.Although the second assumption is almost acceptable, the first assumption is not fulfilled in all situations.Anyway, (7) can be rewritten as the following one: By selecting appropriate  and  weighting matrices, this equation can simply be solved and, according to (6), the constant gain matrix is determined as

Proposed Linear Quadratic Regulator (NLQR)
As we know, by removing control forces of active mass drivers (AMDs), an active control system may passively perform as tuned mass dampers (TMDs).Hence, from control mechanisms point of view, first we assume that the building is controlled by a virtual passive controller which includes a few single degree of freedom systems as TMDs tuned to the building.Then, by employing the similar active control system, in which TMDs were replaced by AMDs, control forces are determined in the way that the responses of the active controlled building are a portion of the similar responses in the passively controlled building.The new proposed method or the new linear quadratic regulator (NLQR) is formed based on this idea.In the first case, we assume that the dynamic equation of motion of the passively controlled building in terms of the state-space variables, shown by , is Now, suppose that in practice the building is actively controlled by AMDs with the same location and characteristics of the passive control system.AMDs apply active control forces to the building such that the building responses are much reduced with respect to the structural responses of the first case, as where () is an admissible transformation matrix.It should be mentioned that the passive control system is a virtual system; therefore, there is no need to know the responses of the passively controlled building.In other words, the transformation matrix is not tuned for a specific external excitation, but it is considered based on the building properties in a way where the building performs better against external disturbances.By using first derivative of ( 11) with respect to time, we have By combining ( 2) and ( 10) in ( 12) and using the transform relation of ( 11) and ( 6), after some manipulations, ( 13) is obtained: For further convenience, assume that the transform matrix of ( 11) is chosen as an identity matrix premultiplied by a scalar time dependent function as follows: Then, the external excitation term in ( 13) is directly related to the passive state-space variable through a simple form: This equation represents the relation between external excitation and passively controlled responses of the structure.By inserting (15) in optimization process of (7) and using (11), as well as ignoring the first derivative of  matrix, the following equation is obtained: This equation is a modified form of (8), in which matrices  and  are changed based on the selected scalar function, ().Now, by proper adjustment of the scalar function, (), and solving (16), time dependent gain matrix is achieved as

Limitations of the Scalar Function 𝛼(𝑡)
With a glance at (16), we find out that the proposed method modifies definition of matrices  and  in (8) to the new forms as follows: Equation ( 18) dictates that the scalar function, (), has to satisfy the following two important criteria: (1) we know that, based on the Lyapunov stability criteria in the continuoustime state-space formulation, an active controlled building will be stable if the eigenvalues of its system matrix have negative real parts.On the other hand, according to the shift method, by using  + (), we forced our system to have eigenvalues with real parts less than −().Hence, if the coefficient of identity matrix, (), added to  in ( 18) is positive, stability of the whole system will be ensured during control time; (2)  matrix is a positive definite matrix.Necessarily, therefore, the coefficient of () in (18) should be positive over the time to satisfy this criterion.

Numerical Example
In order to compare performance of the proposed method (NLQR) in comparison with the classical closed-loop optimal algorithm (CCLQR), an eight-story shear type building is considered.All the properties of the floors are identical; the floor mass of each story is equal to 345.5 tons, the elastic stiffness of each story is 3.404 × 10 5 kN/m, and the internal damping capacity of each story is 2937 tons-sec/m.The active control system includes an active mass damper/driver (AMD) mechanism, which is installed at the roof.The characteristics of AMD are as follows: its mass is 29.63 tons, its tuned frequency is 98% of the predominant vibration frequency of the building without control, and its damping is 25 tonssec/m.

Ground Accelerations.
Performance of the active controlled building is considered during four different ground motions, which include the El-Centro earthquake, the Hyogoken-Nanbu (Kobe) earthquake, the Landers earthquake, and the Parkfield earthquake.These acceleration time histories are shown in Figures 1 to 4. The active control system is designed such that the controlled building presents the best performance during the El-Centro earthquake excitation.Then, behavior of the controlled building during other earthquake excitations is evaluated and their performances are compared together.

Selection of a Proper Scalar Function 𝛼(𝑡).
Based on Section 4, three different functions are presented here as suitable options for the scalar function.The first scalar function is as follows: This function is arbitrarily selected to achieve desirable performances against the El-Centro ground motion.Appropriate value for variable  is equal to 2.
As a second case, another scalar function is proposed as follows: Appropriate value for variable  is set to 0.001.Parameter  0 is an auxiliary variable to set the necessary conditions presented in Section 3, which should be much longer than the length of the control time duration.As a rule,  0 may be almost two or three times greater than the control time.Different than the above functions, the last scalar function is proposed as an acceleration time history dependent function.This function uses time history of the recorded earthquake until that time to find the instantaneous Riccati matrix in each time instant.The proposed function is as follows: where  () = {∫ After extensive analysis, the appropriate values for variables  and  in (21) are set to be 0.1 and 1.5, respectively.This function satisfies the necessary conditions of Section 4 but as it was previously mentioned, the rates of change of the two variables () and () are very rapid.Hence, the Riccati matrix independent of the previous and the next time instant should be found.

Evaluation Criteria.
To evaluate the efficiency of control algorithms, researchers have employed various indices such as the maximum displacement, velocity, or acceleration of the stories, the drift ratios of the adjacent floors, and the maximum base shear of the structures.In this paper, different criteria are developed in order to painstakingly compare the efficiency of different algorithms.Table 1 summarized these categories of criteria, in which indices  1 through  3 represent the criteria for the maximum displacement, velocity, and acceleration responses of the top story, which are normalized to their corresponding uncontrolled values, that is, the structure without any active or passive control systems.
The performance index  4 represents the normalized maximum base shear of the controlled building with respect to the uncontrolled one.Indices  5 through  8 show the root mean square (RMS) of the maximum story responses, such as the displacement, velocity, acceleration, and the story drifts, with respect to the corresponding response quantities in the uncontrolled case.Finally, indices  9 and  10 represent the maximum and average amounts of the required control forces with respect to the classical closedloop optimal control algorithm.Meanwhile, an additional parameter, called efficiency index (E.I.), is defined as the average of the indices  1 through  8 .All of these indices help to show an overall insight on the performances of the various control systems.

Selection of Weighting Matrices.
In many researches [17], variety of  weighting matrix has been suggested.Proper arrangement of the weighting matrix elements and increasing or decreasing its values may significantly affect the stability of the controlled structure and the performance of the control system to achieve the objectives.In order to investigate the performance of the proposed method using various scalar functions, it is assumed that the arrangements of the  weighting matrix are the same for the proposed method and the classical optimal algorithm: where  and  are the matrices, with dimensionless numerical values corresponding to the stiffness and mass matrices of the controlled building, omitting the stiffness and mass values of the active mass damper/driver.Weighting matrix related to the control force, , is assigned to be equal to 1 for all algorithms.Notice that matrix  is a scalar quantity, because only one AMD is installed at the roof level.

Stability Diagram of the Actively Controlled Building.
Since the proposed method results in a time varying gain matrix, stability of the whole controlled building may be changed over the time.In order to inspect this issue, stability diagrams of the building equipped with active and passive systems are depicted together in Figures 5 to 7. Stability diagram of the building, using scalar function  1 (), is plotted in Figure 5.At the beginning of ground motion the structure is more stable while as time goes away poles of actively controlled building move toward passively controlled building poles.On the other hand, in Figure 6, the stability diagram for the case corresponding to scalar function  2 () shows a reverse treatment in comparison with scalar function  1 ().Finally, in Figure 7, the stability diagram of the controlled building using  3 () is depicted.In this figure, a manner similar to case one can be seen.

Response of the Controlled Building during the El-Centro
Earthquake Excitation.In order to examine performance of the proposed control method, linear time history analysis is carried out for all mentioned cases.Maximum acceleration, velocity, and displacement responses of all floors including  passive and active control systems using CCLQR and NLQR with different scenarios during the El-Centro earthquake excitation are compared with the responses of the uncontrolled building as shown in Figure 8 in which the proposed method using various scenarios can significantly decrease all responses of the floors.The best responses belong to the third scalar function  3 (), which is almost the same as the CCLQR algorithm.
Numerical results of the responses of the controlled building are briefly tabulated in Table 2.The maximum and the root mean square of the top story responses, maximum and average required control force, maximum values of the base shear, and the root mean square of the floor drifts for different control systems are presented.Comparing the results shows that the best performance belongs to CCLQR and the results of the new proposed method using the third scalar function  3 () are also very close to it.The only difference is in consuming control forces, which may be observed in columns 8 and 9 of Table 2. Decreasing in maximum control force is only about 1.2% but decreasing in the average (or total) required control force is about 28%.This is the major capability of the proposed method.It is noted   that selecting proper scalar function plays an important role in decreasing the need of the total power supply of the system.

Active Control Forces against the El-Centro
Earthquake Excitation.Active control force time histories of the proposed method using three various scalar functions are compared to that of CCLQR as shown in Figures 9 and 10.In Figure 9, it is observed that the control system using the first scalar function  1 () totally determines active control forces lower than CCLQR.Figure 10 shows that at the start time of control duration the system using the second scalar function  2 () determines active control forces lower than CCLQR, but after a while by changing its poles, control forces are increased.Therefore, it is expected that the performance of the control system using  2 () is better than the previous one but not as good as CCLQR.
Figure 11 shows that the control system using the third scalar function requires the same control forces as CCLQR when the acceleration ground motion is large.Additionally, as the severity of the earthquake is decreased, the control force demand decreases too.
Briefly regarding the above discussion, we can say NLQR with  3 () may present the best performance of the control  system since instant control force follows the severity of the ground excitations.

Response of the Controlled Building during the Other
Earthquake Excitations.Performances of various control systems, which are designed to achieve the high efficiency during the El-Centro earthquake, are examined under a few other earthquakes excitations including the Hyogoken-Nanbu earthquake, the Landers earthquake, and the Parkfield earthquake.After computing, performance indices of the controlled building for each earthquake are listed in Tables 3 to 5. The indices of controlled building using NLQR- 1 () show the lowest efficiency in comparison with the other cases.High values of its separate indices or its averaged index, E.I., confirm this statement.
Generally, the control system using  2 () presents the middle performance between the three scalar functions, except the Landers earthquake, in which the performance is better than NLQR- 3 ().By noticing the earthquake records in Figures 1 to 4, it is seen that the strong part of the Landers record is distributed over the whole duration of the signal, while for the other earthquakes it is concentrated at a short length beyond the beginning.Since NLQR- 2 () increases

Conclusion
In this paper, an approach, called NLQR, is suggested for reducing the peak and average required control forces in actively controlled structures employing optimal control methods.In order to evaluate the effectiveness of the proposed method, the responses of an 8-story shear type

Figure 5 :
Figure 5: Stability diagram of the controlled building using  1 () scalar function, the El-Centro earthquake.

Figure 6 :
Figure 6: Stability diagram of the controlled building using  2 () scalar function, the El-Centro earthquake.

Figure 7 :
Figure 7: Stability diagram of the controlled building using  3 () scalar function, the El-Centro earthquake.

Figure 11 :
Figure 11: Control force time history of NLQR using  3 () scalar function, during the El-Centro earthquake excitation.

Table 2 :
Structural peak response of the controlled building during the El-Centro earthquake excitation.

Table 4 :
Performance indices for building, subjected to Landers ground excitation.
building structure subjected to different ground motions are evaluated.Different types of scalar functions are also utilized for this purpose.The extensive analyses show that NLQR can decrease all responses of the building similar to CCLQR

Table 6 :
Averaged performance indices of the controlled building, subjected to the El-Centro, Hyogoken-Nanbu, Landers, and Parkfield earthquakes.