New Active Control Method Based on Using Multiactuators and Sensors Considering Uncertainty of Parameters

New approach is presented for controlling the structural vibrations. The proposed active control method is based on structural dynamics theories in whichmultiactuators and sensors are utilized. Each actuator force is modeled as an equivalent viscous damper so that several lower vibration modes are damped critically. This subject is achieved by simple mathematical formulation. The proposed method does not depend on the type of dynamic load and it could be applied to control structures with multidegrees of freedom. For numerical verification of proposed method, several criterions such as maximum displacement, maximum kinetic energy, maximum drift, and time history of controlled force and displacement are evaluated in two, five, and seven-story shear buildings, subjected to the harmonic load, impact force, and the Elcentro base excitation. This study shows that the proposed method has suitable efficiency for reducing structural vibrations. Moreover, the uncertainty effect of different parameters is investigated here.


Introduction
Smart structures are systems that can teach and protect themselves against the external excitation such as wind and earthquake.Analyzing and designing of smart structures is based on set of sciences including materials science, applied mechanics, electronics, biomechanics, and structural dynamics.In this procedure, maintaining the structural performance against the external hazards is a very important issue called control system.Many studies have been performed in the field of structural control.These methods can be categorized into three groups, that is, passive, semiactive, and active procedures [1].Due to the simplicity, low cost of assembly and no need to the external power, the passive control systems are numerous.However, the constant control feature makes these systems fail during the earthquakes.In other words, these systems are designed to work only for a certain excitation and limited frequency bound.
The passive control system tries to remove the kinetic energy from the structure.Because of the mentioned constraints in passive algorithms, active control is highly regarded systems to cope with the earthquake.These techniques have suitable efficiency in different excitation, so that they could exactly sense and adopt the structural vibrations.To achieve this goal, each active control method is constructed based on algorithm, which verifies its efficiency and accuracy.The application of such systems began in 1989.In these systems an external power source is required so that this applied force affects the structural equilibrium equation.This applied force may lead to instable vibrations if the active control algorithm is not suitable.
Hence, the complexity, the calculations volume, the instability risk, and the uncertainty factor are some difficulties that arise from active control systems.It should be noted that good performance of active methods depends on some parameters such as the reliable algorithm and the suitable positions for both sensors and actuators [2][3][4][5].
There are several active control mechanism proposed by different researchers that deal with such subjects.In this way, Bayard and his coworkers present the D-optimal design principle which is chosen by the maximum determinant of Fisher information matrix as the criteria function [6].This paper 2 Advances in Civil Engineering simplifies the selected modes into a unitary form by a simple method, so that the suitable position of the piezoelectric elements is achieved.Moreover Kamada and his coworkers modeled a four-story building, which was controlled by piezoelectric actuators, utilizing different algorithms [7].Such study shows that floor accelerations could be reduced up to seventy percent.The subject of determining the suitable position of the piezoelectric actuators was investigated by some other researchers.For example, Han and Lee present a controllable Grammian matrix for the smart composite plates, in which the maximum eigenvalue is used as performance function [8].In this study, the genetic algorithm was utilized to find the effective locations of piezoelectric sensors and actuators.Moreover, Sadri and his coworkers presented some criteria for determining the optimal position of piezoelectric actuators using the controllability of modes [9].In the other research, Gao and his coworkers considered a vibration suppression problem so that the total radiated acoustic power or acoustic potential energy is minimized [10].They used genetic algorithm with immune diversity to evaluate the suitable positions of actuators.Simultaneously, Zhang and his coworkers invested a performance function based on maximizing the dissipated energy that arises from the control action [11].According to this study, a float-encoded genetic algorithm was presented which is capable of solving this optimization problem.Cao and his coworkers used the element sensitivities of singular values to identify the suitable positions of actuators, based on running an optimization process [12].By using topology optimization Kögl and Silva presented an approach to design the piezoelectric plates and shell actuators [13].In this method, the optimization problem consists of distributing the piezoelectric actuators in such a way as to achieve a maximum output displacement in a given direction at a given point of the structure.Moreover, an experimental study on piezoelectric actuators was performed by Sethi and Song [14].They controlled the vibrations of three-story frame, using a piezoelectric patch sensor and actuator.Also, they implemented pole placement modal mechanism to control all vibrations modes.Furthermore than using piezoelectric, there are other kinds of active control algorithms.For example, Song and his coworkers presented the active control mechanism for space truss, using a Lead zirconate Titanate stack actuator [15].The common active control algorithms have been listed in Table 1, followed by the main idea used in each method [16].Finally, the semiactive procedures are achieved based on modifying the passive control systems in combination with active mechanisms [17].
It should be noted that the common active control algorithms use the mathematical concepts such as the genetic algorithms, the Fuzzy logical approaches, the optimization techniques, and other mathematical theories.In these methods, the fundamental principles of structural dynamic which introduce the dynamic behavior of structures were ignored.For this reason, common active control schemes are consistent with structural behavior.The proposed method tries to solve this defect so that the new algorithm is achieved based on a well-known structural dynamics theory, that is, critical damping concept.Based on this theory and also using multiple actuators, a new method is presented here for active control of structures.For this purpose, some fundamental theories of structural dynamics are utilized so that the actuators are modeled as additional viscous dampers in dynamic equilibrium equation.This procedure leads to the actuators forces.Moreover proper positions of both actuators and sensors are determined by an innovate technique.Efficiency of the proposed control method is also evaluated by solving some numerical examples.

The Proposed Active Control Concept Based on Multiactuators and Sensors
Dynamic equilibrium equation of structure can be implemented with various methods such as the Hamilton principle [28]: where M, C, and S are mass, damping, and stiffness matrices of structure, respectively.Furthermore P and D are external force and the nodal displacement vectors, respectively.Also, super dots (⋅) denote differential with respect to time.In active control case, the dynamic equilibrium equation is incorporated into the following relationship: where F a is the equivalent actuator force vector which is generated from the active control mechanism.From the structural dynamics point of view, the structural vibrations damp in the lowest possible time if the structure behaves in the critical damping condition.This is the main concept used here for designing new active control procedure.To explain this idea, consider a multidegree of freedom structure i.e.q.
In modal dynamic analysis, this structure has  vibrations mode.For exposing the structure to behave in fully critical damping conditions, all vibration modes should be in critical damping conditions.From this point of view,  actuators should be attached to structure which causes all elements of the actuator force vector, F a , to be nonzero.It is clear that high number of actuators increases the cost of control process which is not suitable.To prevent this difficulty, few numbers of actuators, that is, ( less than ) is used.Also, the effect of lower dynamic modes is more than higher modes.Therefore  actuators forces are calculated so that  lower modes oscillate in critical damping conditions.Based on the above discussion, the critical damping theory is utilized to determine the equivalent actuators forces.For this purpose, (2) is transformed to the modal space as follows: where   ,   , and   are mass, damping, and stiffness of ith modal coordinates, that is,   , respectively.Also,  is the number of degrees of freedom and   is the ith mode shape vector of free vibration of the structure.If the mode's rank in (3) increases, its effect on dynamic response decreases.Therefore, the primary mode has the highest effect on the dynamic response compared with other vibrations modes.

Active control method Main concept Linear optimal control
Minimize the performance index (Yang 1975 [18]) Pole assignment technique Minimize the performance index (Abdel-Rohman and Leipholz 1978 [19]) Independent modal space control Minimize the modal control force (Meirovitch and Oz 1980 [20]) Instantaneous optimal control Control by minimizing the energy function at each instant of time (Yang et al. 1987 [21]) Bounded state control Keeping the range of responses allowed (Reinhorn et al. 1987 [22]) Nonlinear control Minimizing higher-order function (Wu et al. 1995 [23]) Generalized feedback Control force will be a function of a displacement, velocity, and acceleration (Yang et al. 1991 [24]) Sliding mode control (SMC) Creating the sliding surface (Yang et al. 1994 [25]) Time delay compensation Enter the time delay between the response and performance of control (Abdel-Rohman 1985 [26]) Database and rule base Active control using neural network and fuzzy logic (Tani et al. 1998 [27]) Genetic algorithm Active control using genetic algorithm (Akutagawa et al. 2004 [1]) This principle is utilized to obtain the equivalent actuators forces.Since there are  actuators, attached to the structure, the modal equations could be written as follows: where   is the th element of the th modal shape vector and    is the equivalent actuator force attached to th degree of freedom.The actuators act as additional viscous dampers.This is a model used in the mathematical formulation of the proposed active control method.This model leads to an effective actuator's force which has suitable compatibility with structural behavior.It should be emphasized that modeling the actuator as an equivalent viscous damper does not mean that the actuator's force should be applied to the structure by viscous dampers.In other words, any device which produces such forces is suitable for using in the practical cases.Therefore, the proposed algorithm only presents/calculates the suitable value of the actuator's force in each time of dynamic analysis.Then, this force could be generated by any device of the power source such as piezoelectric.In the following, ( 4) is transformed to the following relationship: where  *  is the th equivalent coordinate damping which is formulated based on both natural structural damping and the effect of the actuators forces: If the equivalent damping coordinates are equal to the critical damping, the structural oscillations damp in the lowest time: where   is the th natural frequency of the structure.It is clear that using  actuators lead to  unknown actuators forces which should be determined at each second.For this purpose, (7) presents a system of equations which could be solved at each analysis time.In the case of existing one actuator substituting ( 6) into ( 7) leads to the following result: Equation ( 8) is completely consistent with the results presented by Alamatian and Rezaeepazhand [29] in which the active control process is formulated based on using only one actuator.Therefore, the proposed method is much more general than existing methods, so that the structural vibrations could be controlled by several actuators.In the case of structural control with two actuators, (7) leads to the below set of equations: where  and  are two degrees of freedom which actuators are attached to them.By solving system (9), the two actuators forcesare obtained: Similar approach could be utilized to formulate the proposed active control method with three actuators.In this case, the Advances in Civil Engineering first three damping coordinates are equal to their corresponding critical values, which leads to the following system of equations: where , , and  are the three degrees of freedom; actuators are attached to them.Based on proposed method the actuators forces can be updated at any second of analysis just by solving a system of equations.There are some unknown parameters in previous systems of equations, that is, ( 9) and ( 11), such as coordinate's velocities ( Ż  ).Based on the structural dynamics theory, modal coordinate velocities depend on the nodal velocities.The coefficient of proportionality is elements of inverse modal shapematrix.Consider where  inv  is th element of the inverse modal shape matrix.To determine the modal velocities, it is necessary to determine both number of sensors and their locations.By increasing number of sensors, the accuracy of modal velocities increases.For example, if there are  sensors attached to the structure, (12) can be written as follows: Here, , , . ..,  are degrees of freedoms, sensors attached to them.It is clear that the proposed method presents the actuator's forces by solving a set of simultaneous equations in each time step.The dimension of this set of equation is equal to the number of actuators.Since number of actuators, attached to the structure, is finite (for practical cases), the dimension of the obtained set of simultaneous equations is quite small so that it solves in a small amount of time, compared with other calculations.Therefore, the required time, spent for calculating the actuator's forces is negligible.In other words, the proposed algorithm could run in real-time fast enough so that it can be utilized for active control of the realistic structures.

The Proper Actuator and Sensor Locations
In this section, the suitable degrees of freedoms actuators and sensors could be attached to them are evaluated.Here fundamental structural dynamics theories are utilized.For example the first mode usually has the highest portion in dynamic response.Therefore, degrees of freedom with high effect in the first mode are suitable for attaching actuators.
In other words, the proper degrees of freedom for actuator locations are those which their corresponding values in the first mode shape are higher than others.For determining the sensors locations, ( 13) is considered.The main criterion utilized for judgment about sensors locations is that reliable values of Ż  are obtained.These quantities correspond to both nodal velocities and elements of  −1 .The nodal velocities are sensed from the dynamic structural response and will be unknown.The only available parameters are nodal velocity coordinates and elements of inverse of modal shape matrix.
The proper sensor locations are determined based on those quantities so that sensors are attached to the degrees of freedom with higher corresponding values in the first row of the inverse matrix of mode shape.For example, in the case of two sensors, they should be installed in degrees of freedom associated with the two highest elements in the first row of the inverse matrix of mode shape.Since there are two sensors in the smart structure, the first and second modal coordinate's velocity can be calculated as follows: where  1 and  2 are th element in the first and second row of the inverse modal shape matrix.Also, Ḋ  and and Ḋ are the velocity of th and th degree of freedom which the sensor are attached to them.

The Proposed Active Control Algorithm
To verify the proposed new active control method, some numerical dynamic analyses are performed.For this purpose, the analysis time is divided into limited number of time steps.In each time step, numerical time integration scheme will be utilized to achieve structural responses.Here, Newmark method with linear acceleration is used.The main steps of proposed active control process are as follows.(f) Compute the displacements vector of the current time steps using Newmark method with linear acceleration [30].
(g) Calculate the structural velocity vector based on the Newmark method [30].
(h) Calculate the acceleration vector of th time step by solving the following linear system:  (k) Print the results and end.

Numerical Study
To verify the validity of the proposed active control method, some numerical examples are presented.For this purpose, the suggested active control process is combined with numerical dynamic analysis methods, that is, Newmark integration using linear acceleration scheme.

Two Degree of Freedom System
. Figure 1 shows a linear two-DOF system, subjected to impact load that is,  2 defined as follows: where Δ is the time step.The procedures, which leads to the optimum sensor and actuators locations, are summarized in Table 2.By applying the proposed method, the second degree of freedom is the suitable location for attaching both the actuator and sensor, that is, S2-A2, due to its highest corresponding values of  inv 12 = 0.8167 and  21 = 0.8594.After determining the sensor and actuators locations, the proposed active control method is applied to this structure.Figures 2 and 3 show the displacement time responses of the first and second degrees of freedom for different control process, respectively.It is clear that all control methods are stable.On the other hand, the proposed method reduces the vibrations amplitude of both degrees of freedom in a short time.Therefore, the suggested process has suitable efficiency in vibration control.If the actuator and sensor are attached to the first degree of freedom (case S1-A1), the efficiency of the control process is reduced.Also, the control case S1-A2 has lower efficiency than the two other cases.In fact S2-A1 is incapable of controlling the system.This subject clearly proves that the proposed algorithm for determining the actuator and sensors locations works property.

Five Story Shear Building.
Figure 4 shows a five-story shear building, modeled by lumped mass and lateral stiffness (five horizontal degrees of freedom).This structure is analyzed in two different cases.In the first analysis, a harmonic load that is, () = 50 sin(10 ), is applied to the fifth-story building.The second analysis is performed when the structure is excited by the Elcentro ground acceleration record.
In both analyses, the damping ratio of 5% for the first mode is assumed for constructing the Rayleigh damping matrix with two factors [31].To control the vibrations of the structure, the location of the sensors and actuators are determined based on the proposed algorithm.Table 3 presents the details which lead to the optimum locations of sensors and actuators.Using the results presented in Table 3, the various control algorithms with regard to the number of actuators and sensors have been inserted in Table 4.
Three factors are considered in numerical evaluation of the proposed method: the maximum displacements of stories, the maximum kinetic energy of the structure, and the maximum drift.Table 5 shows the maximum structural displacement for the harmonic load.Structure with one actuator (degree of freedom of 5) and one sensor (degree of freedom of 4) A 5,4 -S 4 Structure with 2 actuators and one sensor A 5,4,3 -S 4 Structure with 3 actuators and 1 sensor A 5,4,3 -S 4,3 Structure with 3 actuators and 2 sensor A 5,4,3 -S 4,3,5 Structure with 3 actuators and 3 sensor A 5,4,3 -S 4,3,5,2 Structure with 3 actuators and 4 sensor A 5,4,3 -S 4,3,5,2,1 Structure with 3 actuators and 5 sensor According to the result of Table 5, in the case of using one actuator and one sensor, the displacement of upper floors can be reduced to about 65%.If the structural vibrations are controlled by using two actuators and one sensor the displacement of upper floors can be diminished about 65%.Moreover, the mentioned reduction rate will reach to 75% in the case of three actuators and two sensors.If three sensors and actuators are used, the maximum displacement in fourth and fifth floor will be reduced to about 80 and 75 percent, respectively.Another factor, utilized for the numerical evaluation of the proposed method is the maximum kinetic energy of structural dynamic analysis.Table 6 shows the maximum kinetic energy of fifth-story shear building subjected to the harmonic load.Accordingly, the proposed control approach has suitable performance in reducing the kinetic energy of the system so that by increasing number of sensors and actuators the kinetic energy is reduced considerably.Also, the variation of the fifth-story drift is different in various control cases, inserted in Table 7.It is clear that the performance of the proposed method increases by using more numbers of actuators.
To assess the proposed algorithm against the broadband earthquake excitations, this structure is analyzed for the seismic load, that is, the Elcentro earthquake accelerograms.The maximum displacements of the fourth and fifth floors, caused by the El Centro Earthquake, have been inserted in Table 8. Results of Table 8 show that if the structural vibration is controlled by one actuator and one sensor, the displacements of upper floors are reduced to about 42%.Moreover, the mentioned reduction will reach to about 70% if three actuators and two sensors are used.
Figure 5 shows the time history displacement of the fifthstory for different control processes.Moreover, the variation of the fifth floor acceleration for both controlled and controloff system has been plotted in Figure 6.Based on Figures 5 and 6 the proposed algorithm, presented for determining the locations of the sensor and actuator and also calculating the actuator's force, has a suitable applicability in numerical results.
Moreover, the time history for both of the earthquake and the actuator forces (for the case of one actuator and one sensor) have been plotted in Figure 7.It is clear that the maximum actuators force is less than fifty percent of the floors weight.Moreover, Figure 7 shows that there is logical balance between the actuator's force and the earthquake load.

Seven-Story Shear Building.
Here, a seven-story shear building is analyzed so that the lumped mass and lateral stiffness of each story are 0.2591 kg and 100 N/cm, respectively.This structure is excited by the Elcentro acceleration.To control the vibration of the structure, the sensors and   actuators locations are determined on proposed algorithm.The details of such procedure have been inserted in Table 9.Using the various control algorithms with regard the number of actuators and sensors, the maximum displacements of the sixth and seventh have been in Table 10.Results that the proposed control approach is suitable for reducing the maximum displacements of the system.Moreover, time history displacement the structure has been plotted in Figure 8 for different control cases.Figure 8 shows that efficiency of the proposed active control scheme increases if more number of sensor and actuator is used.
All of the above dynamic analyses are completely stable so that any instability the control process of proposed algorithm is not considered.It is worth that  the proposed active control mechanism does not impose any additional condition for stability of the analysis, because it is formulated based on the well-known structural dynamic theories.From this point of view, the prepared algorithm has suitable compatibility with dynamic behavior of structures.
This subject reduces the instability potential that exists in any active control method.

The Effect of Uncertainty in the Proposed Method
There are many variables in the analysis and design of structures such as loads, capacity of elements, and material properties.Uncertainty in each of these variables has a significant effect on structure safety.Stories mass in this section are considered as random variables.To evaluate the uncertainty of stories mass in the proposed method, for each story, we considered a In the range of +15% and −15%, a set of numbers is generated by MATLAB as stories masses.In the stories masses in each of the cases considered purely (case1 until case 30) and each story may be reduced or increased to the base 11 shows the results of the proposed method using a sensor and an actuator, in various states of random masses in the range 15 percent compared to the base masses.According to the received results, changes in masses in the range of 15% caused displacement increment up to 14.6% and displacement reduction up to 10% for fourth floor.Also, in this condition, the displacement of fifth floor, increased up to 18% and decreased up to 14%.

Applicability of the Proposed Method
In the case of applicability of the proposed method in realistic structures, using of the piezoelectric stack could be useful for a lot of reasons.In this way the method which has been introduced by Kamada et al. [7] is proposed.In this method piezoelectric stack actuators are placed at the bottom of firststory columns in both sides.If the induced voltage is applied in reverse phase to these actuators, the concentrated bending moments   are produced at the bottom of columns.By applying the basic structural calculations, the equivalent shear force at each story level is determined as follows: where ,   , and   are the columns height, actuators produced force per unite voltage, and equivalent shear force produced by actuators in first floor level, respectively.Other parameters have been shown in Figure 9.

Conclusions
In this paper, a new method was developed for active control of structure.This process was based on the theory of structural dynamics.In the proposed model, the structural fluctuations are controlled by multiple actuators and sensors and the actuator has been modeled as an additional viscous damper.In this way, the actuators forces are calculated in such a way that a greater number of low structural damping of the various modes of vibration is critical.Actuators forces are calculated in such a way that more first equivalent modal coordinates damping is critical.To consider the proposed method, the changes in maximum displacement, the maximum kinetic energy, and total maximum relative displacement of stories were considered in various control cases in a five-story structure.Based on the numerical results, in the case of using an actuator and one sensor, the displacement of upper floors can be reduced to about 65% compared to the case of lack of control system.Also, the displacement in upper stories reduced in the case of using two actuators and one sensor.Moreover, in the case of three actuators and two sensors at the proper locations the maximum upper stories displacement reduced up to 75%.In term of maximum kinetic energy, using one actuator and one sensor reduced this item up to 90%.Beside, using three actuators and sensors caused the total of stories diminish up to 70%.In addition, the uncertainty effect the mass of stories on the floor's displacement was (a) Set  = 0 and select the time step of dynamic analysis.(b) Construct the stiffness, mass, and damping matrices.(c) Determine the actuators and sensors locations using modal shape matrix and its inverse.(d) Calculate the modal velocities.(e) Solve set of equation to determine the actuators forces.

Figure 2 :
Figure 2: The displacement-time response to the first degree of freedom of 2D system.

Figure 3 :
Figure 3: The displacement-time response to the second degree of freedom of 2D system.

Figure 5 :
Figure 5: The time history displacement for fifth-story of shear building excited by the Elcentro earthquake.

Figure 6 :
Figure 6: history acceleration for fifth-story of shear building excited by the Elcentro earthquake.

Figure 7 : 4 Figure 8 :
Figure 7: Comparison between the actuator and the earthquake forces in the case of using one actuator and one sensor.

Figure 9 :
Figure 9: The piezoelectric actuator's mechanism (a) front view and (b) side view.

Table 1 :
Common Active Control algorithms.

Table 2 :
Suitable sensor and actuator locations.

Table 3 :
The actuator and sensor locations for the five-story shear building.

Table 4 :
Different control cases for the five-story shear building.

Table 5 :
The maximum structural displacements for the harmonic load.

Table 6 :
The maximum kinetic energy for the five-story shear building subjected to the harmonic load.

Table 7 :
The story drift of the fifth floor for the five-story shear building subjected to the harmonic load.

Table 8 :
The maximum displacement of the fourth and fifth floors in the five-story shear building excited by the Elcentro earthquake.

Table 9 :
The actuator and sensor locations for seven-story shear building.

Table 10 :
The maximum displacements of the seven-story shear building excited by the Elcentro earthquake.

Table 11 :
The effect of uncertainty on fourth and fifth floor displacement.