Interval Analysis of the Eigenvalues of Closed-Loop Control Systems with Uncertain Parameters

: Uncertainty caused by a parameter measurement error or a model error causes difficulties for the implementation of the control method. Experts can divide the uncertain system into a definite part and an uncertain part and solve each part using various methods. Two uncertainty problems of the control system arise: problem A for the definite part—how does one find out the optimal number and position of actuators when the actuating force of an actuator is smaller than the control force? Problem B for the uncertain part—how does one evaluate the effect of uncertainty on the eigenvalues of a closed-loop control system? This paper utilizes an interval to express the uncertain parameters and converts the control system into a definite part and an uncertain part using interval theory. The interval state matrix is constructed by physical parameters of the system for the definite part of the control system. For Problem A, the paper finds out the singular value element sensitivity of the modal control matrix and reorders the optimal location of the actuators. Then, the paper calculates the state feedback gain matrix for a single actuator using the receptance method of pole assignment and optimizes the number and position of the actuators using the recursive design method. For Problem B, which concerns the robustness of closed-loop systems, the paper obtains the effects of uncertain parameters on the real and imaginary parts of the eigenvalues of a closed-loop system using the matrix perturbation theory and interval expansion theory. Finally, a numerical example illustrates the recursive design method to optimize the number and location of actuators and it also shows that the change rate of eigenvalues increases with the increase in uncertainty.


Introduction
As a result of the complexity of structural, manufacturing and installation errors, and the inaccuracy of measurements, structural parameters are inevitably uncertain. This uncertainty may lead to a deterioration in vibration control performance. Therefore, the uncertainty of structural parameters plays an important role in vibration control [1]. However, uncertainty is not equal to stochasticity, and the stochastic method is not the only method for studying uncertain problems. In recent years, non-stochastic methods have been considered to analyze uncertain problems. The convex model theory is one such approach and has been applied by Dief to linear programming and system theory [2]. The convex model theory has attracted the attention of control theory and structural analysis workers. For example, Lindberg [3] discussed the dynamic response and failure problems of structures with the convex model, and Ganzerli [4] discussed the optimization problem of structure with uncertain parameters. Since the 1960s, the interval method [5] proposed by Moore has become an important tool for the analysis of uncertainty problems. Since the interval method requires little structural information and only needs to know the upper and lower bounds of structural parameters, it has attracted the attention of many scholars. The interval finite element method [6] is proposed by combining interval analysis with the finite element method, which makes it easier for interval analysis to deal with the analysis of complex structures. Great attention has also been paid to the problem of structural vibration control in engineering practice. Experts mainly discuss robustness when studying uncertain system control. Mori and Kokame [7] discussed the necessary and sufficient conditions for the dynamic stability of an uncertain system. Rachid [8] discussed the robust control of uncertain parameter systems.
The essence of vibration control is the assignment of poles and zero. The active control method has been applied to the assignment of poles and zero since the 1960s. Wonham [9] concluded that arbitrary poles can be assigned by state feedback when the structure is controllable. Datta [10] showed the explicit solution of the feedback gain matrix for calculating part of the eigenvalue assignment. In this algorithm, arbitrarily selected poles can be assigned for any controllable structure while the remaining poles remain unchanged. Duan [11] established two methods to obtain simple parametric equations of feedback gain and closed-loop eigenvalues with the aid of singular value decomposition and the factorization method. Ram and Mottershead [12] devised the receptance method of pole placement, where the vibration of the structure can be controlled and suppressed using the Sherman-Morrison formula and by means of measurement.
The configuration position of the actuator/sensor affects the performance of the closed-loop control system. Having fewer actuators and sensors is better for the sake of economic benefit. Therefore, it is of great value to determine the optimal number of actuators and sensors and to configure them at the optimal position. Many in-depth studies have been undertaken in this field [13][14][15][16][17][18][19][20][21][22][23][24] and can be divided into two aspects: one is to determine the optimal allocation criterion, or objective function to be optimized; the other is to select an appropriate optimization method. The optimal configuration criteria include the following: the maximum configuration criterion based on the singular value of the controllable matrix; the maximum configuration criterion of controllability/observability defined by the Grimm matrix; the configuration standards based on the minimum energy consumption of the system; the minimum configuration criterion based on control/observation overflow; the configuration criteria based on failure and reliability; the configuration standards based on the maximum response value in the time domain or the maximum response value in the frequency domain, etc. The optimal allocation algorithms are as follows: the non-linear planning optimization method, the sequential method, the simulated annealing method, the genetic algorithm, and the particle swarm optimization algorithm.
The objectives of the paper are as follows: (1) to establish the finite element equations for the uncertain system and divide them into a definite part and an uncertain part where uncertain parameters are treated as an interval using the interval theory in Section 2; (2) to find out the optimal position for each actuator according to the singular value element sensitivity of the modal control matrix in Section 3; (3) taking into consideration that the actual force is smaller than the ideal control force, to optimize the location and position for each actuator using the recursive method in Section 4; (4) to analyze the effect of uncertain parameters on the eigenvalues of the closed-loop system in Section 5; (4) to carry out an example to illustrate the proposed methods in Section 6.

Finite Element Equations for Uncertain Systems
The following is the finite element equation of the intelligent structure with uncertain parameters [25][26][27][28][29][30]: Divide the mass matrix, damping matrix, and stiffness matrix into definite parts and uncertain parts [25] to obtain: where 0 M , 0 C , and 0 K are the definite parts of the mass matrix M , the damping matrix C , and stiffness matrix K , respectively; and  M ,  C , and  K are the uncertain parts of the aforementioned matrices, respectively.
Divide the displacement vector   t Ξ , control vector   t u , control matrix D , and mechanical into definite parts and uncertain parts in the same way [25], and we obtain: are the uncertain parts of the above.
All through the paper, assume that the uncertainty is small, and that all uncertain parameters are expressed by the first-order Taylor approximation method according to the natural interval expansion principle and interval operation. Substituting Equations (2)-(8) into Equation (1), we have: The uncertain system equation (Equation (1)) is expressed by both the definite part equation (Equation (9)) and the uncertain part equation (Equation (10)).

Interval Mass Matrix, Interval Damping Matrix, and Interval Stiffness Matrix
Assume that the structural parameter γ has some degrees of error or uncertainty and is expressed as [5]: where 1 n is the number of uncertain parameters. γ can be expressed by the following interval: Equation (10) can also be expressed as:

State Equation of the Uncertain System with/without Feedback Control
The state equation for the uncertain structure equation (Equation (1)) is: where The state equation (Equation (20)) of the uncertain structure in the closed-loop system can also written as: The corresponding eigenvalue problem for Equation (21) is: where ˆj S is the -th j order complex eigenvalue of the closed-loop system with an uncertain structure; ˆj W and ˆj V are the corresponding right and left eigenvectors; the superscript H is the conjugate transpose; and ˆj S is the conjugate of ˆj S .
The control system of the uncertain structure is as follows: where 0 A and 0 B are the definite parts of A and B , respectively; and  A and  B are the uncertain parts: The state vector and feedback gain matrix G are as follows: where   0 t X and 0 G are the definite parts of the state vector and the feedback gain matrix, respectively; and   t  X and  G are the uncertain parts.
By substituting Equations (24)- (27) into Equation (21) and ignoring the terms of the higher order, we obtain: We divided the control system into the definite parts, as in Equations (28) and (29), and the uncertain parts, as in Equations (30) and (31).
From Equations (28) and (29), the state equation of the definite part of the closed-loop system can be written as: The corresponding eigenvalue problem for the definite part of the closed-loop system is:   where 0 j S is the complex eigenvalues of the -th j order for definite closed-loop systems; 0 j W and 0 j V are the corresponding right and left eigenvectors; the superscript H is the conjugate transpose; and 0 j S is the conjugate of 0 j S .
Systems with uncertain parameters can be considered as perturbations of systems with certain parameters. The interval response where   0 t X is determined by Equation (32).
The interval state matrix and interval control matrix are [25]: The interval gain matrix I G is: in Equation (35) is determined by the following equation: In order to study the response of the system with uncertain parameters, we treat the control problem of the definite parameter system using Equation (32), and then study the influence of the uncertain parameters on the control using Equation (39).
In Section 3, the paper proposes a method to judge which position is the optimal position for a single input. In Section 4, a recursive design method is introduced to optimize the location and number of actuators when the actual control force is smaller than the ideal control force. In Section 5, the paper discusses the effect of uncertainty on eigenvalues of the closed-loop system using interval analysis.

The Measure of the Modal Controllability of Intelligent Structures
Equation (9) is the vibration control equation of the system with certain parameters. Although the original equation has n modes, the vibration of sub-systems with lower eigenvalues is much larger than that with higher eigenvalues, and we only need to control the vibration of the sub-systems with lower r eigenvalues and ignore the vibration of those with higher eigenvalues. Generally the (r + 1)-th eigenvalue is much larger than exciting frequency. Transfer Equation (9) into the modal coordinate system for the purpose of modal control and, according to the modal expansion theorem, we have: where is the mechanical forces in the modal coordinates.
A problem is that Equation (41) cannot reflect the direct relationship between control force   0 t u and mode controllability. To solve this problem, decompose modal control matrix * 0 r n R   D using the Singular Value Decomposition (SVD) method: where the columns of the orthogonal matrix 0 A new set of coordinates is introduced to measure the controllability of modes by 0 Σ : Substitute Equation (44) into Equation (41), and we obtain: where ) are the mass matrix, damping matrix, stiffness matrix, and mechanical force in the 0 0 p Φ space, respectively; . As

Singular Value Element Sensitivity of Modal Control Matrix
When the number of controlled modes is r and the number of actuators is n , from Equation (42), we have [17,31]: Remove the actuator on the -th i unit from the current structure, and Equation (46) becomes: It is necessary to perform a perturbation calculation for Equation (46) in order to obtain the optimal position of the actuator using the topology optimization method. According to perturbation theory, the eigenvalues and eigenvectors represented by Equation (48) can be expressed as the power series of the small parameter  : where 01 i j  is the element sensitivity contribution to the singular value of the -th j modal control matrix when the actuator on -th i element is removed.

Control Position Optimization Criteria
For intelligent structures, when the -th i actuators are removed, the singular value and singular vector of each order will change. When the actuator on the -th i unit is removed, we define ai  as follows [17]: which is a measure of actuator position optimization.  dimensionless. The criteria for intelligent structure actuator location optimization is: where Ω  is a given threshold. If the optimization criteria in Equation (56) is met, the actuator on the -th i element is preserved. Otherwise, the actuator on the -th i element is removed from the structure. Let 1 i i   and repeat the above steps to determine the optimal number and location of actuators. Obviously, the bigger Ω  is, the fewer actuators are placed.

Required Control Force Calculation Based on the Receptance Method
The state equation for the definite closed-loop system is Equation (32). In this section, the receptance method will be used to calculate the control force and feedback gain matrix for a single input [32].

Considering the condition of a single input,
where   0 t ξ is modal coordinate vector, Transfer Equation (58) into the Laplace domain, and we have: The Sherman-Morrison formula [12] is as follows:

H S B G H S H S H S G H S B
, Substitute Equation (62) into Equation (29), and we have: and When the structure is fully controllable for the single-input system, from Equations (63) and (64), we can calculate the required control force   0 t u and the state feedback gain matrix 0 G .

Optimal Number and Location of Actuators when the Actual Control Force Is Smaller than the Required Control Force
Generally, the actuators are the same and the range of execution force for each one is: where a1 is the maximum execution force. We will discuss the optimal number and location of the actuators, with special consider for the condition that the actual control force is smaller than the required force. u t  , a multi-input control should be introduced to make the structure fully controllable [33].
We will determine the minimum number and optimal location of actuators required using the recursive design method for a multi-input modal controller.
Adopt the state feedback control, and we have: Substitute Equation (67) into Equation (66), and we have: The complex eigenvalue is 0 Λ  , compare it to the eigenvalue that is to be configured. When , stop the recursive design and  is the minimum number of actuators.

Interval Analysis of the Robustness of Closed-Loop Systems
A closed-loop control system with uncertain parameters can be written as [25]: where A  is actually the state matrix with actuators. For any I  γ γ , the eigenvalue problem is: Apply the interval expansion method, replace the real variate I γ with the corresponding interval variate I γ , and replace the real operation with the corresponding interval operation, and Equations (71) and (72) become:   In other words, we obtain the following set: According to central representation of the interval, the interval matrix Substituting Equations (83)-(86) into Equations (81) and (82), the first-order perturbations of complex eigenvalues are: We can obtain the interval eigenvalue from Equation (83): The first-order perturbation of the eigenvalue (Equation (89)) is: Complex eigenvalues 1 j S   can be written as: The upper and lower bounds of the real and imaginary parts of complex eigenvalues are: The following condition can be used to estimate the stability of the control system: where j  is the real part of the eigenvalues for an uncertain closed-loop system. When designing the feedback control, if is large enough, we can guarantee the stability of the uncertain parameter system.

Numerical Example
Without loss of generality, we will control the first modal of a composite cantilever beam model, as shown in Figure 1. The length, width, and thickness of the main structure are 720 mm, 20 mm, and 10 mm, respectively. The mass density of the main structure is 7850 kg/m 3 and the Young's modulus is 210 GPa. Two thin layers of piezoelectric material (PVDF) are completely glued to the upper and lower surfaces to function as actuators and sensors. The parameters of the PVDF material are as follows: the thickness of each layer is 0.5 mm; the mass density is 1800 kg/m 3 ; the Young's modulus is 2 × 10 9 Pa; the dielectric constant (polarization along the depth of the cantilever beam) is 33  = 0.1062 × 10 -9 ; and the piezoelectric constant is e31 = 0.046. The example will show that the location of the actuator affects the performance of the vibration control. Assuming that the actuating force for the actuator is 50% of the ideal control force, which means the actual gain is 50% of the ideal gain, the example also shows how to find out the optimal location and number of actuators required to control the vibration using the recursive design method. The example also demonstrates the impact of uncertain parameters on the structure using interval analysis.  Finite element equations for main structure are established using Euler-Bernoulli beam theory, as the ratio of length to height is 72. The main structure is divided into nine uniform beam elements with 10 nodes along the direction of length. Each node of this finite element model has two degrees of freedom, one translation and one rotation. The rotational degree of freedom is associated with the cross-section rotation. The left node is fixed and the others are free to move. Furthermore, we take the damping matrix 0.0012  C K .
This means that when controlling the first-order mode, the optimal position of the control force is arranged on the first element 1 b ; the second optimal position is arranged on the second element 2 b ; the third optimal position is arranged on the third element 3 b ...; the worst position is on the ninth element 9 b .
The first two complex eigenvalues of the open-loop system are: 10 6.1531 101.2676i This means that the change in frequency of the control force in the optimal position 11   is 15,610 times that of the worst position 91   . In other words, if the pole assignments are the same, the energy consumed when placing the actuator in the worst position is much greater than when placing the actuator in the optimal position. This illustrates the conclusions in Section 3.1. Therefore, it is important to select the position of the control force in the structural vibration control in practice. 2. The optimal number and location of actuators.
Assume the actual gains provided by each actuator are which means the actual force is 50% of the ideal control force. The optimal number and location of the actuators will be calculated using the recursive method proposed in Section 4.
Compare 11  and 12  with 1  and 2  , and we have −9.2065 > −12.2600 and 111.0296 < 120, an additional actuator is required and the second step is as follows:  The second step: The second actuator   2 u t  is arranged in the optimal position 2 b , and we have 21 11.3349 117.3545i 22 11.3349 117.3545i Compare 21  and 22  with 1  and 2  , and we have −11.3349 > −12.2600 and 117.3545 < 120, an additional actuator is still required and the third step is as follows:  The third step: is arranged in the optimal position 3 b , and we have: 31 12.7140 121.2772i 32 12.7140 121.2772i Compare 31  and 32  with 1  and 2  , and we have −12.7140 < −0.1226 and 140.8476 > 140, and stop the calculation. Three actuators are needed to completely achieve the vibration control of the structure. In addition, they are arranged in positions 1 b , 2 b , and 3 b .
3. The effect of uncertain parameters on eigenvalues of the closed-loop control system. Now, we have determined the number and location of the required actuators. We can substitute them into the system and discuss the impact of uncertain parameters on the structure to illustrate Section 5.
Suppose that the Young's modulus , system with uncertain parameters using perturbation theory and interval analysis and is illustrated with a numerical example.
In conclusion, this paper solves the problem of the optimal number and position of actuators when the control force is greater than the range of the actuator's execution force. The structure can be fully controllable by increasing the number of actuators. Interval theory is used to discuss the effect of uncertain parameters on eigenvalues of closed-loop control systems.