Next Article in Journal
Chopstick Robot Driven by X-shaped Soft Actuator
Next Article in Special Issue
Modeling of Attractive Force by Magnetic Wheel Used for Mobile Robot
Previous Article in Journal
Development of Haptic Stylus for Manipulating Virtual Objects in Mobile Devices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

Key Laboratory of CNC Equipment Reliability (Jilin University), Ministry of Education, School of Mechanical and Aerospace Engineering, Jilin University, Changchun 130022, China
*
Author to whom correspondence should be addressed.
Actuators 2020, 9(2), 31; https://doi.org/10.3390/act9020031
Submission received: 14 March 2020 / Revised: 10 April 2020 / Accepted: 15 April 2020 / Published: 21 April 2020
(This article belongs to the Special Issue Actuators for System Identification, Vibration Analysis, and Control)

Abstract

:
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.

1. 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.

2. Fundamental Equations for the Interval Analysis of Intelligent Structures

2.1. 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]:
M Ξ ¨ ( t ) + C Ξ ˙ ( t ) + K Ξ ( t ) = D u ( t ) + F ( t ) ,
where M R n × n and K R n × n are the mass and stiffness matrices of the intelligent structure, respectively; C R n × n is the proportional damping matrix, C = α M + β K , where α and β are the mass and stiffness proportional damping constants, respectively (more information can be found in [29]); D R n × n is the control matrix determined by the position and number of actuators; u ( t ) R n × 1 is the control vector; Ξ ( t ) R n × 1 is the displacement vector; and F ( t ) R n × 1 is the mechanical force.
Divide the mass matrix, damping matrix, and stiffness matrix into definite parts and uncertain parts [25] to obtain:
M = M 0 + δ M ,
C = C 0 + δ C ,
K = K 0 + δ K ,
where M 0 , C 0 , and K 0 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 u ( t ) , control matrix D , and mechanical force F ( t ) into definite parts and uncertain parts in the same way [25], and we obtain:
Ξ ( t ) = Ξ 0 ( t ) + δ Ξ ( t ) ,
u ( t ) = u 0 ( t ) + δ u ( t ) ,
D = D 0 + δ D ,
F ( t ) = F 0 ( t ) + δ F ( t ) ,
where Ξ 0 ( t ) , u 0 ( t ) , D 0 , and F 0 ( t ) are the definite parts of the displacement vector Ξ ( t ) , control vector u ( t ) , control matrix D , and mechanical force F ( t ) , respectively; δ C , δ Ξ ( t ) , δ u ( t ) , δ D , and δ F ( t ) 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:
M 0 Ξ ¨ 0 ( t ) + C 0 Ξ ˙ 0 ( t ) + K 0 Ξ 0 ( t ) = D 0 u 0 ( t ) + F 0 ( t ) ,
δ M Ξ ¨ 0 ( t ) + M 0 δ Ξ ¨ ( t ) + δ C Ξ ˙ 0 ( t ) + C 0 δ Ξ ˙ ( t ) + δ K Ξ 0 ( t ) + K 0 δ Ξ ( t ) = δ D u 0 ( t ) + D 0 δ u ( t ) + δ F ( t ) .
The uncertain system equation (Equation (1)) is expressed by both the definite part equation (Equation (9)) and the uncertain part equation (Equation (10)).

2.2. 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]:
γ = [ γ 1 , γ 2 , , γ n 1 ] T ,
where n 1 is the number of uncertain parameters.
γ can be expressed by the following interval:
γ = γ c + δ γ γ I = [ γ , γ ] = [ γ c Δ γ , γ c + Δ γ ] = [ γ c , γ c ] + [ Δ γ , Δ γ ] = γ c + Δ γ I ,
where δ γ is the uncertain part of γ ; γ I = [ γ , γ ] is the interval parameter; γ = γ c Δ γ is the left end point of interval γ I ; γ = γ c + Δ γ is the right end point of interval γ I ; γ c = ( γ + γ ) / 2 is the midpoint of interval γ I ; Δ γ = ( γ γ ) / 2 is the radius of interval γ I .
For any γ γ I , when the interval parameter γ has a small disturbance around the median value γ c , the first-order approximations of M ( γ ) , C ( γ ) , and K ( γ ) are obtained by expansion using the Taylor approximation near γ c [25].
We also have [5]:
M ( γ I ) = M ( γ c ) + Δ M ( γ I ) ,
C ( γ I ) = C ( γ c ) + Δ C ( γ I ) ,
K ( γ I ) = K ( γ c ) + Δ K ( γ I ) ,
where M ( γ I ) , C ( γ I ) , and K ( γ I ) are the interval mass matrix, the interval damping matrix, and the interval stiffness matrix, respectively; M ( γ c ) , C ( γ c ) , and K ( γ c ) are the median mass matrix, the median damping matrix, and the median stiffness matrix of the structure, respectively. Δ M ( γ I ) , Δ C ( γ I ) , and Δ K ( γ I ) are the first-order increment matrix of the interval mass matrix, interval damping matrix, and interval stiffness matrix, respectively:
Δ M ( γ I ) = l = 1 n 1 M ( γ c ) γ l I ( γ l I γ l c ) ,
Δ C ( γ I ) = l = 1 n 1 C ( γ c ) γ l I ( γ l I γ l c ) ,
Δ K ( γ I ) = l = 1 n 1 K ( γ c ) γ l I ( γ l I γ l c ) ,
Equation (10) can also be expressed as:
Δ M I Ξ ¨ 0 ( t ) + M 0 Δ Ξ ¨ I ( t ) + Δ C I Ξ ˙ 0 ( t ) + C 0 Δ Ξ ˙ I ( t ) + Δ K I Ξ 0 ( t ) + K 0 Δ Ξ I ( t ) = Δ D I u 0 ( t ) + D 0 Δ u I ( t ) + Δ F I ( t ) ,

2.3. State Equation of the Uncertain System with/without Feedback Control

The state equation for the uncertain structure equation (Equation (1)) is:
X ˙ ( t ) = A X ( t ) + B u ( t ) + η ( t ) ,
where X ( t ) = [ Ξ T ( t ) , Ξ ˙ T ( t ) ] T is the state vector; u ( t ) = G T X ( t ) is the input vector of the state feedback control; and G R p × 2 n is the state feedback gain matrix. The state matrix A R 2 n × 2 n , and the distribution matrix of the actuators B R 2 n × p and η ( t ) are:
A = [ 0 I M 1 K M 1 C ] ,   B = [ 0 D ] ,   η ( t ) = [ 0 F ( t ) ] .
The state equation (Equation (20)) of the uncertain structure in the closed-loop system can also written as:
X ˙ ( t ) = ( A + B G T ) X ( t ) + η ( t ) .
The corresponding eigenvalue problem for Equation (21) is:
( A + B G T ) W ^ j = S ^ j W ^ j ,
( A + B G T ) H V ^ j = S ^ ¯ j V ^ j ,
where S ^ j is the j -th order complex eigenvalue of the closed-loop system with an uncertain structure; W ^ j and V ^ j are the corresponding right and left eigenvectors; the superscript H is the conjugate transpose; and S ^ ¯ j is the conjugate of S ^ j .
The control system of the uncertain structure is as follows:
A = A 0 + δ A ,
B = B 0 + δ B ,
where A 0 and B 0 are the definite parts of A and B , respectively; and δ A and δ B are the uncertain parts:
A 0 = [ 0 I M 0 1 K 0 M 0 1 C 0 ] , δ A = [ 0 I ( M 0 1 δ K + δ M 1 K 0 ) ( M 0 1 δ C + δ M 1 C 0 ) ] , B 0 = [   0 D 0 ] , δ B = [   0 δ D ] .
The state vector X ( t ) and feedback gain matrix G are as follows:
X ( t ) = X 0 ( t ) + δ X ( t ) ,
G = G 0 + δ G ,
where X 0 ( t ) and G 0 are the definite parts of the state vector and the feedback gain matrix, respectively; and δ X ( t ) and δ G are the uncertain parts.
By substituting Equations (24)–(27) into Equation (21) and ignoring the terms of the higher order, we obtain:
X ˙ 0 ( t ) = A 0 X 0 ( t ) + B 0 u 0 ( t ) + η 0 ( t ) ,
u 0 ( t ) = G 0 T X 0 ( t ) ,
δ X ˙ ( t ) = A 0 δ X ( t ) + δ A X 0 ( t ) + B 0 δ u ( t ) + δ B u 0 ( t ) ,
δ u ( t ) = G 0 T δ X ( t ) + δ G T X 0 ( t ) .
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:
X ˙ 0 ( t ) = ( A 0 + B 0 G 0 T ) X 0 ( t ) + η 0 ( t ) .
The corresponding eigenvalue problem for the definite part of the closed-loop system is:
( A 0 + B 0 G 0 T ) W ^ 0 j = S ^ 0 j W ^ 0 j ,
( A 0 + B 0 G 0 T ) H V ^ 0 j = S ^ ¯ 0 j V ^ 0 j ,
where S ^ 0 j is the complex eigenvalues of the j - th order for definite closed-loop systems; W ^ 0 j and V ^ 0 j are the corresponding right and left eigenvectors; the superscript H is the conjugate transpose; and S ^ ¯ 0 j is the conjugate of S ^ 0 j .
Systems with uncertain parameters can be considered as perturbations of systems with certain parameters. The interval response X I ( t ) is:
X I ( t ) = X 0 ( t ) + Δ X I ( t ) ,
where X 0 ( t ) is determined by Equation (32).
The interval state matrix and interval control matrix are [25]:
A I = A 0 + Δ A I ,
B I = B 0 + Δ B I ,
where A I = A ( γ I ) and B I = B ( γ I ) are the interval state matrix and interval control matrix; Δ A I = Δ A ( γ I ) is the first-order increment matrix of the interval state matrix; and Δ B I = Δ B ( γ I ) is the first-order increment matrix of the interval control matrix:
Δ A I = [ 0 I M 0 1 Δ K I ( Δ M I ) 1 K 0 M 0 1 Δ C I ( Δ M I ) 1 C 0 ] , Δ B I = [ 0 Δ D I ] .
The interval gain matrix G I is:
G I = G 0 + Δ G I .
Δ X I ( t ) in Equation (35) is determined by the following equation:
Δ X ˙ I ( t ) = A 0 Δ X I ( t ) + Δ A I X 0 ( t ) + B 0 Δ u I ( t ) + Δ B I u 0 ( t ) .
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.

3. Optimization of Actuator Position in Vibration Control of Intelligent Structure

3.1. 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:
Ξ 0 ( t ) = p 0 q 0 ( t ) ,
where q 0 ( t ) R r × 1 is the modal coordinate vector; p 0 R n × r is the modal matrix; and r is the number of modes intercepted.
Substitute Equation (40) into Equation (9) and multiply both sides of the equation by p 0 T to the left, and we have:
M 0 * q ¨ 0 ( t ) + C 0 * q ˙ 0 ( t ) + K 0 * q 0 ( t ) = D 0 * u ( t ) + F 0 * ( t ) ,
where M 0 * = p 0 T M 0 p 0 is the modal mass matrix, M 0 * R r × r ; C 0 * = p 0 T C 0 p 0 is the modal damping matrix, C 0 * R r × r ; K 0 * = p 0 T K 0 p 0 is the modal stiffness matrix, K 0 * R r × r ; D 0 * = p 0 T D 0 R r × n is the modal control matrix; and F 0 * ( t ) = p 0 T F 0 ( t ) is the mechanical forces in the modal coordinates.
A problem is that Equation (41) cannot reflect the direct relationship between control force u 0 ( t ) and mode controllability. To solve this problem, decompose modal control matrix D 0 * R r × n using the Singular Value Decomposition (SVD) method:
D 0 * = Φ 0 Σ 0 Ψ 0 T ,
where the columns of the orthogonal matrix Φ 0 R r × r ( Φ 0 T Φ 0 = I r ) are the left singular vector of D 0 * ; the columns of the orthogonal matrix Ψ 0 R n × n ( Ψ 0 T Ψ 0 = I n ) are the right singular vector of D 0 * ; Σ 0 r = d i a g ( σ 01 , σ 02 , , σ 0 r ) is the diagonal matrix of non-zero elements in Σ 0 ; σ 01 σ 02 σ 0 r > 0 ; σ 0 j ( j = 1 , 2 , , r ) is the singular value of D 0 * ; and I r is the identity matrix of the order r .
Multiply both sides of Equation (42) by Φ 0 T and consider D 0 * = p 0 T D 0 , and we obtain:
Φ 0 T p 0 T B 0 = Σ 0 V 0 T .
A new set of coordinates q ˜ 0 ( t ) is introduced to measure the controllability of modes by Σ 0 :
Ξ 0 ( t ) = p 0 Φ 0 q ˜ 0 ( t ) .
Substitute Equation (44) into Equation (41), and we obtain:
M ˜ 0 * q ˜ ¨ 0 ( t ) + C ˜ 0 * q ˜ ˙ 0 ( t ) + K ˜ 0 * q ˜ 0 ( t ) = Σ 0 u ˜ 0 ( t ) + F ˜ 0 * ( t ) ,
where M ˜ 0 * = Φ 0 T p 0 T M 0 p 0 Φ 0 ( M ˜ 0 * R r × r ), C ˜ 0 * = Φ 0 T p 0 T C 0 p 0 Φ 0 ( C ˜ 0 * R r × r ), K ˜ 0 * = Φ 0 T p 0 T K 0 p 0 Φ 0 ( K ˜ 0 * R r × r ), and F ˜ 0 * ( t ) = Φ 0 T p 0 T F 0 ( t ) ( F ˜ 0 * ( t ) R r × 1 ) are the mass matrix, damping matrix, stiffness matrix, and mechanical force in the p 0 Φ 0 space, respectively; u ˜ 0 ( t ) = Ψ 0 T u 0 ( t ) . As Ψ 0 T Ψ 0 = I n , u 0 ( t ) and u ˜ 0 ( t ) are equivalent in a modular sense. The controllability of q ˜ 0 ( t ) can be measured by σ 0 j in Σ 0 : the larger σ 0 j is, the stronger the control effect on q ˜ 0 ( t ) and the less energy is consumed.

3.2. Singular Value Element Sensitivity of Modal Control Matrix D 0 *

When the number of controlled modes is r and the number of actuators is n , from Equation (42), we have [17,31]:
D 0 * D 0 * T Φ 0 = Φ 0 Σ 0 Σ 0 T .
Remove the actuator on the i - th unit from the current structure, and Equation (46) becomes:
D ˜ 0 * i D ˜ 0 * i T Φ ˜ 0 i = Φ ˜ 0 i Σ ˜ 0 i Σ ˜ 0 i T ,
or
D ˜ 0 * i D ˜ 0 * i T Φ ˜ 0 j i = Φ ˜ 0 j i ( σ ˜ 0 j i ) 2 .
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 ε :
D ˜ 0 * i = D 0 * + ε D 01 * i ,
σ ˜ 0 j i = σ 0 j + ε σ 01 j i + ε 2 σ 02 j i + ,
Φ ˜ 0 j i = Φ 0 j + ε Φ 01 j i + ε 2 Φ 02 j i + ,
where D ˜ 0 * i is the modal control matrix when moving the actuator on the i - th element; the j - th column Φ ˜ 0 j i of orthogonal matrix Φ ˜ 0 i is the left singular vector of the j - th order; when removing the actuator on the i - th element, σ 0 j and Φ 0 j are the j - th order singular value and j - th order left singular vector of the original system, respectively; σ 01 j i and σ 02 j i are respectively the first-order and second-order perturbations of singular value σ 0 j when removing the actuator on the i - th element; and Φ 01 j i and Φ 02 j i are, respectively, the first-order and second-order perturbations of the j - th order left singular vector Φ 0 j when removing actuator on the i - th element.
To solve Equation (48), the first-order perturbation momentum of eigenvalues is:
σ 01 j i = Φ 0 j T ( B 0 B 01 i T + B 01 i T B 0 ) Φ 0 j ,
where σ 01 j i is the element sensitivity contribution to the singular value of the j - th modal control matrix when the actuator on i - th element is removed.

3.3. Control Position Optimization Criteria

For intelligent structures, when the i - th actuators are removed, the singular value and singular vector of each order will change. When the actuator on the i - th unit is removed, we define Ω a i as follows [17]:
Ω a i = j = 1 r σ 01 j i ,   Ω a = [ Ω a 1 , Ω a 2 , , Ω a n ] ,
which is a measure of actuator position optimization. Ω a i is the sum of the singular value sensitivity contributions to the modal control matrix. The dimensionless optimization metric Ω ˜ a i of the i - th actuator position and sorted in the form of matrix Ω ˜ a is as follows:
Ω ˜ a i = Ω a i / max ( Ω a 1 , Ω a 2 , , Ω a n ) ,   ( i = 1 , 2 , , n ) ,
Ω ˜ a = [ Ω ˜ a 1 , Ω ˜ a 2 , , Ω ˜ a n ] ,
where max ( Ω a 1 , Ω a 2 , , Ω a n ) is the largest absolute value of element in Ω a ; and Ω ˜ a is also dimensionless. The criteria for intelligent structure actuator location optimization is:
Ω ˜ a Ω ˜ ,
where Ω ˜ is a given threshold. If the optimization criteria in Equation (56) is met, the actuator on the i - th element is preserved. Otherwise, the actuator on the i - th element is removed from the structure. Let i = i + 1 and repeat the above steps to determine the optimal number and location of actuators. Obviously, the bigger Ω ˜ is, the fewer actuators are placed.

4. The Recursive Design Method of Multi-Input Modal Controller

4.1. 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, B 0 is the column vector and G 0 T is the row vector. With modal transformation X 0 ( t ) = W 0 ξ 0 ( t ) , we have the following modal governing equation:
ξ ˙ 0 ( t ) = ( Λ 0 + B 0 * G 0 * T ) ξ 0 ( t ) + η 0 * ( t ) ,
where ξ 0 ( t ) is modal coordinate vector, ξ 0 ( t ) R 2 r × 1 ; Λ 0 = V 0 H A 0 W 0 is eigenvalue matrix of the state matrix Λ 0 R 2 r × 2 r ; W 0 and V 0 are the right and left mode matrices of the state matrix A 0 , respectively; and V 0 H U 0 = I r ; B 0 * = V 0 H B 0 is modal control matrix in state space; and G 0 * T = G 0 T W 0 is modal gain matrix in state space.
Equation (57) can be written as:
ξ ˙ 0 ( t ) ( Λ 0 + B 0 * G 0 * T ) ξ 0 ( t ) = η 0 * ( t ) .
Transfer Equation (58) into the Laplace domain, and we have:
( S Λ 0 B 0 * G 0 * T ) ξ 0 ( S ) = η 0 * ( S ) .
The Sherman–Morrison formula [12] is as follows:
H ^ ( S ) = H ( S ) H ( S ) B 0 * G 0 * T H ( S ) 1 + G 0 * T H ( S ) B 0 * ,  
where H ^ ( S ) = ( S Λ 0 B 0 * G 0 * T ) 1 is the receptance matrix of the closed-loop system. H ( S ) = ( S Λ 0 ) 1 is the receptance matrix of the open-loop system. The characteristic polynomials of closed-loop systems are 1 + G 0 * T H ( S ) B 0 * . If we want to assign the poles to Λ 0 * = [ λ 01 * , λ 02 * , , λ 02 r * ] , let the denominator of Equation (60) be zero:
1 + G 0 * T H ( Λ 0 * ) B 0 * = 0 .
The feedback gain matrix G 0 * can be obtained by solving Equation (61). Multiply X 0 ( t ) = W 0 ξ 0 ( t ) by V 0 H , and we have:
ξ 0 ( t ) = V 0 H X 0 ( t ) .
Substitute Equation (62) into Equation (29), and we have:
u 0 ( t ) = G 0 * T V 0 H X 0 ( t ) = G 0 T X 0 ( t ) ,
and
G 0 T = G 0 * T V 0 H .
When the structure is fully controllable for the single-input system, from Equations (63) and (64), we can calculate the required control force u 0 ( t ) and the state feedback gain matrix G 0 .

4.2. 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:
u ˜ 01 ( t ) [ a 1   a 1 ] ,  
where a 1 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.
  • If | u 01 ( t ) | < a 1 , one actuator can provide the feedback gain required by the controlled structure in excess;
  • If | u 01 ( t ) | = a 1 , one actuator just provides the feedback gain required by the controlled structure;
  • If | u 01 ( t ) | > a 1 , the single-input control cannot make the structure fully controllable and multiple actuators are needed to provide sufficient feedback gain to make the structure fully controllable.
For the case of | u 01 ( t ) | a 1 , only a single input can satisfy the requirement. However, for the case of | u 01 ( t ) | > a 1 , 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. Assume b 1 , b 2 , …, b Θ are the optional, sub-optional, and the τ - th optional positions for τ actuators. Construct the state equation of the structure as:
X ˙ 0 ( t ) = A 0 X 0 ( t ) + τ = 1 Θ B 0 τ u ˜ 0 τ ( t ) + η 0 ( t ) .
Adopt the state feedback control, and we have:
u ˜ 0 τ ( t ) = G ˜ 0 τ T X 0 ( t ) .
Substitute Equation (67) into Equation (66), and we have:
X ˙ 0 ( t ) = ( A 0 + τ = 1 Θ B 0 τ G ˜ 0 τ T ) X 0 ( t ) + η 0 ( t ) ,
where u ˜ 0 τ ( t ) is the actual input vector in the τ - th step; B 0 τ is distribution matrix of the actuator in the τ - th step; G ˜ 0 τ is the actual gain provided by the actuator in the τ - th step. Then, we can solve the following state equation:
X ˙ 0 ( t ) = ( A 0 + τ = 1 Θ B 0 τ G ˜ 0 τ T ) X 0 ( t ) .
The complex eigenvalue is Λ ˜ 0 τ , compare it to the eigenvalue that is to be configured. When Λ ˜ 0 τ = Λ 0 * , stop the recursive design and τ is the minimum number of actuators.

5. Interval Analysis of the Robustness of Closed-Loop Systems

A closed-loop control system with uncertain parameters can be written as [25]:
X ˙ ( t ) = A ˜ X ( t ) + η ( t ) ,
where A ˜ is actually the state matrix with actuators.
For any γ γ I , the eigenvalue problem is:
A ˜ ( γ ) W ˜ j = S ˜ j W ˜ j ,
A ˜ H ( γ ) V ˜ j = S ˜ ¯ j V ˜ j ,
where A ˜ ( γ ) is the closed-loop system state matrix; S ˜ j is the j - th order eigenvalue of A ˜ ( γ ) ; W ˜ j is the j - th order eigenvector of A ˜ ( γ ) ; A ˜ H ( γ ) is the conjugate transpose of A ˜ ( γ ) ; S ˜ ¯ j is the conjugate of S ˜ j ; and V ˜ j is the j - th order eigenvector of A ˜ H ( γ ) .
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:
A ˜ ( γ I ) W ˜ j I = S ˜ j I W ˜ j I ,  
A ˜ H ( γ I ) V ˜ j I = S ˜ ¯ j I V ˜ j I ,  
where A ˜ ( γ I ) is the interval state matrix for a closed-loop system; S ˜ j I is the j - th order interval eigenvalue of A ˜ ( γ I ) ; W ˜ j I is the j - th order interval eigenvector of A ˜ ( γ I ) ; A ˜ H ( γ I ) is the conjugate transposition of matrix A ˜ ( γ I ) ; S ˜ ¯ j I is conjugate of S ˜ j I ; V ˜ j I is the j - th order interval eigenvector of A ˜ H ( γ I ) . Equations (73) and (74) are interval eigenvalue problems. The interval eigenvalue S ˜ j I gives the interval state matrix A ˜ ( γ I ) ; therefore, we cannot only obtain the minimum interval, but also obtain S ˜ j thus satisfying the equation A ˜ ( γ ) W ˜ j = S ˜ j W ˜ j . In other words, we obtain the following set:
Γ = { S ˜ j : A ˜ ( γ ) W ˜ j = S ˜ j W ˜ j , A ˜ ( γ ) A ˜ ( γ I ) } ,  
Γ = { S ˜ ¯ j : A ˜ H ( γ ) V ˜ j = S ˜ ¯ j V ˜ j , A ˜ H ( γ ) A ˜ H ( γ I ) } .
According to central representation of the interval, the interval matrix A ˜ ( γ I ) and A ˜ H ( γ I ) are:
A ˜ ( γ I ) = A ˜ ( γ c ) + Δ A ˜ ( γ I ) ,
A ˜ H ( γ I ) = A ˜ H ( γ c ) + Δ A ˜ H ( γ I ) ,
where A ˜ ( γ c ) = A ˜ 0 , Δ A ˜ ( γ I ) = Δ A ˜ . The interval eigenvalue problem of the system is:
[ A ˜ ( γ c ) + Δ A ˜ ( γ I ) ] W ˜ j I = S ˜ j I W ˜ j I ,  
[ A ˜ H ( γ c ) + Δ A ˜ H ( γ I ) ] V ˜ j I = S ˜ ¯ j I V ˜ j I .
For any γ γ I , there are always δ A ˜ ( γ ) Δ A ˜ ( γ I ) and δ A ˜ H ( γ ) Δ A ˜ H ( γ I ) , and the corresponding characteristic equations are:
[ A ˜ ( γ c ) + δ A ˜ ( γ ) ] W ˜ j = S ˜ j W ˜ j ,  
[ A ˜ H ( γ c ) + δ A ˜ H ( γ ) ] V ˜ j = S ˜ ¯ j V ˜ j .
According to complex mode perturbation theory, complex eigenvalues and complex eigenvectors of the structure are:
S ˜ j = S ˜ 0 j + δ S ˜ 1 j + δ S ˜ 2 j + ,  
S ˜ ¯ j = S ˜ ¯ 0 j + δ S ˜ ¯ 1 j + δ S ˜ ¯ 2 j + ,  
W ˜ j = W ˜ 0 j + δ W ˜ 1 j + δ W ˜ 2 j + ,  
V ˜ j = V ˜ 0 j + δ V ˜ 1 j + δ V ˜ 2 j + .
Substituting Equations (83)–(86) into Equations (81) and (82), the first-order perturbations of complex eigenvalues are:
δ S ˜ 1 j = V ˜ 0 j H Δ A ˜ ( γ ) W ˜ 0 j ,  
δ S ˜ ¯ 1 j = W ˜ 0 j H Δ A ˜ H ( γ ) V ˜ 0 j .
We can obtain the interval eigenvalue from Equation (83):
S ˜ j I = S ˜ 0 j + Δ S ˜ 1 j I .
The first-order perturbation of the eigenvalue (Equation (89)) is:
Δ S ˜ 1 j I = ( V ˜ 0 j I ) H Δ A ˜ ( γ I ) W ˜ 0 j I = [ Δ S ˜ 1 j , Δ S ˜ 1 j ] .
Complex eigenvalues Δ S ˜ 1 j can be written as:
Δ S ˜ 1 j = | Re ( Δ S ˜ 1 j ) | + i | Im ( Δ S ˜ 1 j ) | .
The upper and lower bounds of the real and imaginary parts of complex eigenvalues are:
Re ( S ˜ j ) = Re ( S ˜ 0 j ) | Re ( Δ S ˜ 1 j ) | ,   Re ( S ˜ j ) = Re ( S ˜ 0 j ) + | Re ( Δ S ˜ 1 j ) | ,
Im ( S ˜ j ) = Im ( S ˜ 0 j ) | Im ( Δ S ˜ 1 j ) | ,   Im ( S ˜ j ) = Im ( S ˜ 0 j ) + | Im ( Δ S ˜ 1 j ) | .
The following condition can be used to estimate the stability of the control system:
Re ( S ˜ 0 j ) | Re ( Δ S ˜ 1 j ) | Θ j Re ( S ˜ 0 j ) + | Re ( Δ S ˜ 1 j ) | ,  
where Θ j is the real part of the eigenvalues for an uncertain closed-loop system. When designing the feedback control, if Re ( S ˜ 0 j ) ( k = 1 , 2 , , r ) is large enough, we can guarantee the stability of the uncertain parameter system.

6. 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/m3 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/m3; the Young’s modulus is 2 × 109 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 C = 0.0012 K .
1.
The calculation of the optimal position of the control force that controls the first-order mode under the condition of realizable control.
From Equation (55) in Section 3.2, we see that the dimensionless sensitivity of the first-order modes of elements 1~9 is:
Ω ˜ a = [ 1.0000 ,   0.6970 ,   0.4517 ,   0.2652 ,   0.1356 ,   0.0565 ,   0.0167 ,   0.0025 ,   6 × 10 5 ] .
This means that when controlling the first-order mode, the optimal position of the control force is arranged on the first element b 1 ; the second optimal position is arranged on the second element b 2 ; the third optimal position is arranged on the third element b 3 …; the worst position is on the ninth element b 9 .
The first two complex eigenvalues of the open-loop system are:
λ 10 = 6.1531 101.2676 i ,
λ 20 = 6.1531 + 101.2676 i .
Suppose the first two complex eigenvalues of the open-loop system, λ 10 and λ 20 , are configured on the poles μ 1 and μ 2 of the closed-loop system:
μ 1 = 12.2600 120.0000 i ,
μ 2 = 12.2600 + 120.0000 i .
Adopt velocity and displacement feedback control, and the control force u 1 ( t ) should be applied at the optimal position b 1 , the velocity feedback gain at g c 1 = 166,860 , and the displacement feedback gain at g k 1 = 67,951 . If g k 1 is the displacement feedback gain for each actuator, when the control force acts at element b i , ω ˜ i 1 is the first frequency of the close-loop system, and it can be sorted as ω ˜ 1 :
ω ˜ 1 = [ 120.0000 ,   114.6481 ,   110.1239   106.5570 ,   104.0065 ,   102.4169 101.6097 ,   101.3195 ,   101.2688 ] .
Let Δ ω i 1 = ω ˜ i 1 ω 0 , where ω 0 = 101.2676 rad/s is the first frequency of the open-loop system:
Δ ω 1 = [ 18.7324 ,   13.3805 ,   8.8563 5.2894 ,   2.7389 ,   1.1493 0.3421 ,   0.0519 ,   0.0012 ] .
Consider Δ ω 11 = 18.7324 rad/s, Δ ω 91 = 0.0012 rad/s, we have:
Δ ω 11 Δ ω 91 = 18.7324 0.0012 = 15,610 .
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 g ˜ c 1 max = 0.5 g c 1 and g ˜ k 1 max = 0.5 g k 1 , 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.
  • The first step:
When the actuator u ˜ 1 ( t ) is arranged in the optimal position b 1 , we have
μ 11 = 9.2065 111.0296 i ,
μ 12 = 9.2065 + 111.0296 i .
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 u ˜ 2 ( t ) is arranged in the optimal position b 2 , and we have
μ 21 = 11.3349 117.3545 i ,
μ 22 = 11.3349 + 117.3545 i .
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:
A third actuator u ˜ 3 ( t ) is arranged in the optimal position b 3 , and we have:
μ 31 = 12.7140 121.2772 i ,
μ 32 = 12.7140 + 121.2772 i .
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 b 1 , b 2 , and b 3 .
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 E I = [ E c Δ E , E c + Δ E ] , mass density ρ I = [ ρ c Δ ρ , ρ c + Δ ρ ] , and height h I = [ h c Δ h , h c + Δ h ] are taken as interval parameters of the closed-loop system. Then, the upper and lower bounds of complex eigenvalues can be obtained according to interval stiffness, interval mass, and interval damping. Δ Re = | Re ( Δ S ˜ ) | / | Re ( S ˜ c ) | % and Δ Im = | Im ( Δ S ˜ ) | / | Im ( S ˜ c ) | % are the rates of change of the real and imaginary parts of eigenvalues, respectively. When the Young’s modulus of the main structure is Δ E = ( 2 % , 4 % , 6 % ) E c , the upper and lower bounds of complex eigenvalues of the first-order uncertain closed-loop system and the change rates of the real and imaginary parts of eigenvalues are respectively:
μ 311 = [ 12.8370 , 12.5909 ] + [ 122.1199 i , 120.4286 i ] ,   Δ Re 1 = 0.9678 % ,   Δ Im 1 = 0.6973 % ,
μ 312 = [ 12.9601 , 12.4679 ] + [ 122.9568 i , 119.5741 i ] ,   Δ Re 2 = 1.9357 % ,   Δ Im 2 = 1.3948 % ,
μ 313 = [ 13.0832 , 12.3448 ] + [ 123.7880 i , 118.7133 i ] ,   Δ Re 3 = 2.9039 % ,   Δ Im 3 = 2.0926 % .
When the mass density of the main structure is Δ ρ = ( 2 % , 4 % , 6 % ) ρ c , the upper and lower bounds of complex eigenvalues of the first-order uncertain closed-loop system and the change rates of the real and imaginary parts of the eigenvalues are:
μ 311 = [ 12.9734 , 12.4647 ] + [ 122.5085 i , 120.0823 i ] ,   Δ Re 1 = 1.9998 % ,   Δ Im 1 = 1.0001 % ,
μ 312 = [ 13.2437 , 12.2250 ] + [ 123.7780 i , 118.9221 i ] ,   Δ Re 2 = 3.9998 % ,   Δ Im 2 = 2.0008 % ,
μ 313 = [ 13.5255 , 11.9943 ] + [ 125.0879 i , 117.7948 i ] ,   Δ Re 3 = 6.0000 % ,   Δ Im 3 = 3.0027 % .
When the height of the main structure Δ h = ( 2 % , 4 % , 6 % ) h c , the upper and lower bounds of complex eigenvalues of the first-order uncertain closed-loop system and the change rates of the real and imaginary parts of the eigenvalues are:
μ 311 = [ 13.0813 , 12.3516 ] + [ 123.3007 i , 119.2539 i ] ,   Δ Re 1 = 2.8691 % ,   Δ Im 1 = 1.6684 % ,
μ 312 = [ 13.4535 , 11.9942 ] + [ 125.3242 i , 117.2306 i ] ,   Δ Re 2 = 5.7345 % ,   Δ Im 2 = 3.3368 % ,
μ 313 = [ 13.8307 , 11.6417 ] + [ 127.3478 i , 115.2075 i ] ,   Δ Re 3 = 8.5936 % ,   Δ Im 3 = 5.0052 % .
These results show that when uncertain parameters increase, the change rate of eigenvalues increases, that is, the error of eigenvalues increases with the increase in uncertain parameters. For example, when the Young’s modulus of the structure is Δ E = 2 % E c , the change rates of the real and imaginary parts of the first eigenvalue are 0.9678% and 0.6973%, respectively. When the Young’s modulus of the main structure is Δ E = 6 % E c , the change rates of the real part and imaginary part of the first eigenvalue are 2.9039% and 2.0926%, respectively. When the height of the main structure is Δ h = 2 % h c , the change rates of the real part and imaginary part of the first eigenvalue are 2.8691% and 1.6684%, respectively. When the height of the main structure is Δ h = 6 % h c , the change rates of the real part and imaginary part of the first eigenvalue are 8.5936% and 5.0052%, respectively. In addition, on the whole, the uncertainty of the height of the structure has a greater influence on the eigenvalues than the Young’s modulus or mass density.

7. Conclusions

  • This paper discusses the measurement of modal controllability of intelligent structures, deduces a control matrix singular value sensitivity formula, and determines the optimal position of piezoelectric actuators for an intelligent structure. In a numerical example, when the sizes of the actuators are exactly the same and the displacement feedback gains are exactly the same, it is very important to select the location of actuators in a structural vibration control;
  • When a single actuator is applied, the feedback gain matrix is calculated using the receptance method, and the ideal control force of actuator is calculated. When the control force required for complete controllability of the structure is much higher than the actual force of the actuator, a single-input control cannot fully control the structure, and a multi-input control is needed for the purpose of vibration control. The recursive design method of a modal controller is introduced to study the optimal number and location of the actuators;
  • The change rate of eigenvalues of a closed-loop control system increases with the increase in the uncertain parameters. Uncertain parameters are expressed as intervals. A method is proposed to solve the upper and lower bounds of the eigenvalues of a closed-loop control 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.

Author Contributions

All authors contributed equally and significantly in writing this article. Writing-original draft, J.-Z.Z., G.-F.Y., R.-Y.L., Y.-C.Z., K.-Y.G. and M.W.; Writing-review and editing, J.-Z.Z., G.-F.Y., R.-Y.L., Y.-C.Z., K.-Y.G. and M.W. All authors have read and agreed to the published version of the manuscript.

Funding

Jilin Provincial science and technology development program (20170520098JH) and Natural Science Foundation of Jilin Province (20160101280JC).

Acknowledgments

The authors thank CRRC Changchun Railway Vehicles Co., LTD. for supporting this study.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Lin, Y.K. Probabilistic Theory of Structural Dynamics; McGraw-Hill: New York, NY, USA, 1967. [Google Scholar]
  2. Dief, A.S. Sensitivity Analysis in Linear Systems; Springer: Berlin, Germany, 1986. [Google Scholar]
  3. Lindberg, H.E. Dynamic response and buckling failure measures for structures with bounded and random imperfections. Trans. ASME J. Appl. Mech. 1991, 58, 1092–1094. [Google Scholar] [CrossRef]
  4. Ganzerli, S.; Pantelides, C.P. Load and resistance convex models for optimum design. Struct. Optim. 1999, 17, 259–268. [Google Scholar] [CrossRef]
  5. Moore, R.E. Methods and Applications of Interval Analysis; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1979. [Google Scholar]
  6. Dessombz, O.; Thouverez, F.; Laine, J.P.; Jezequal, L. Analysis of mechanical system using interval computations applied to finite element methods. J. Sound Vibr. 2001, 239, 949–968. [Google Scholar] [CrossRef] [Green Version]
  7. Mori, T.; Kokame, H. Convergence property of interval matrices and interval polynomials. Int. J. Contr. 1987, 45, 481–484. [Google Scholar] [CrossRef]
  8. Rachid, A. Robustness of discrete systems under structural uncertainties. Int. J. Contr. 1989, 50, 1563–1566. [Google Scholar] [CrossRef]
  9. Wonham, W.M. On pole assignment in multi-input controllable linear system. IEEE Trans. Autom. Control 1967, 12, 660–665. [Google Scholar] [CrossRef] [Green Version]
  10. Datta, B.N.; Elhay, S.; Ram, Y.M. Orthogonality and partial pole assignment for the symmetric difinite quadratic pencil. Linear Algebra Its Appl. 1997, 257, 29–48. [Google Scholar] [CrossRef]
  11. Duan, G.R. Parametric eigenstructure assigment in second-order descriptor linear systems. IEEE Trans. Autom. Control 2004, 49, 1789–1795. [Google Scholar] [CrossRef]
  12. Mottershead, J.E.; Ram, Y.M. Receptance method in active vibration control. AIAA J. 2007, 3, 562–567. [Google Scholar]
  13. Ram, Y.M.; Mottershead, J.E. Multiple-input active vibration control by partial pole placement using the method of receptances. Mech. Syst. Signal Process. 2013, 40, 727–735. [Google Scholar] [CrossRef]
  14. Mokrani, B.; Palazzo, F.; Mottershead, J.E.; Fichera, S. Multiple-input multiple-output experimental aeroelastic control using a receptance-based method. AIAA J. 2019, 7, 3066–3077. [Google Scholar] [CrossRef]
  15. Tsai, S.H.; Ouyang, H.; Chang, J.Y. Identification of torsional receptances. Mech. Syst. Signal Process. 2019, 126, 116–136. [Google Scholar] [CrossRef]
  16. Gao, X. Research on Theory and Algorithm of Local Vibration Control of Structure. Ph.D. Thesis, Jilin University, Changchun, China, 2009. [Google Scholar]
  17. Chen, S.; Cao, Z. A new method for determining locations of the piezoelectric Sensor/Actuator for vibration control of intelligent structures. J. Intell. Mater. Syst. Struct. 2000, 11, 108–115. [Google Scholar] [CrossRef]
  18. Lim, T.W. Actuator/sensor placement for modal parameter identification of flexible structures. Modal Anal. Int. J. Anal. Exp. Modal Anal. 1993, 1, 1–13. [Google Scholar]
  19. Schulz, G.; Heimbold, G. Dislocated actuator/sensor positioning and feedback design for flexible Structures. J. Guid. Control. Dyn. 1982, 6, 361–367. [Google Scholar] [CrossRef]
  20. Rao, S.S.; Pan, T. Optimal placement of actuators in actively controlled structures using genetic algorithm. AIAA J. 1991, 5, 942–943. [Google Scholar] [CrossRef]
  21. Levine, W.S.; Athans, M. On the determination of the constant output feedback gains for linear multivariable systems. Trans. Autom. Control 1970, 15, 44–48. [Google Scholar] [CrossRef]
  22. Kang, Y.K.; Park, H.C.; Hwang, W.; Han, K.S. Optimum placement of piezoelectric sensor/actuator for vibration control of laminated beams. AIAA J. 1996, 9, 1921–1926. [Google Scholar] [CrossRef]
  23. Choe, K.; Baruh, H. Actuator placement in structural control. J. Guid. Control Dyn. 1991, 1, 40–46. [Google Scholar] [CrossRef]
  24. Baruh, H.; Merovitch, L. On the placement of actuators in the control of distributed-parameter systems. AIAA J. 1981, 611–618. [Google Scholar] [CrossRef]
  25. Chen, S.H. Matrix Perturbation Theory in Structural Dynamics; Science Press: Beijing, China, 2007. [Google Scholar]
  26. Zienkiewicz, O.C.; Taylor, R.L. The Finite Element Method; Butterwort-Heinemann: Oxford, UK, 2000. [Google Scholar]
  27. Craig, R.R. Structural Dynamics: An Introduction to Computer Methods; Wiley: New York, NY, USA, 1981. [Google Scholar]
  28. Rong, B.; Rui, X.; Lu, K.; Tao, L.; Wang, G.; Ni, X. Transfer matrix method for dynamics modeling and independent modal space vibration control design of linear hybrid multibody system. Mech. Syst. Signal Process. 2018, 104, 589–606. [Google Scholar] [CrossRef]
  29. Neto, M.A.; Ambrósio, J.A.C.; Roseiro, L.M.; Amaro, A.; Vasques, C.M.A. Active vibration control of spatial flexible multibody systems. Multibody Syst. Dyn. 2013, 30, 13–35. [Google Scholar] [CrossRef]
  30. Tzou, H.S.; Tseng, C.I. Distributed piezoelectric sensor/actuator design for dynamic measurement/control of distributed parameter system: A piezoelectric finite element approach. J. Sound Vib. 1990, 138, 17–34. [Google Scholar] [CrossRef]
  31. Preumont, A. Mechatronics Dynamics of Electromechnical and Piezoelectric Systems; Springer: Dordrecht, The Netherlands, 2006. [Google Scholar]
  32. Mokrani, B.; Batou, A. The minimum norm multi-input multi-output receptance method for partial pole placement. J. MSSP 2019, 129, 437–448. [Google Scholar] [CrossRef]
  33. Yao, G.; Gao, X. Controllability of regular systems and defective systems with repeated eigenvalues. J. Vib. Control. 2010, 17, 1574–1581. [Google Scholar]
Figure 1. Composite cantilever beam structure.
Figure 1. Composite cantilever beam structure.
Actuators 09 00031 g001

Share and Cite

MDPI and ACS Style

Zhao, J.-Z.; Yao, G.-F.; Liu, R.-Y.; Zhu, Y.-C.; Gao, K.-Y.; Wang, M. Interval Analysis of the Eigenvalues of Closed-Loop Control Systems with Uncertain Parameters. Actuators 2020, 9, 31. https://doi.org/10.3390/act9020031

AMA Style

Zhao J-Z, Yao G-F, Liu R-Y, Zhu Y-C, Gao K-Y, Wang M. Interval Analysis of the Eigenvalues of Closed-Loop Control Systems with Uncertain Parameters. Actuators. 2020; 9(2):31. https://doi.org/10.3390/act9020031

Chicago/Turabian Style

Zhao, Jing-Zhou, Guo-Feng Yao, Rui-Yao Liu, Yuan-Cheng Zhu, Kui-Yang Gao, and Min Wang. 2020. "Interval Analysis of the Eigenvalues of Closed-Loop Control Systems with Uncertain Parameters" Actuators 9, no. 2: 31. https://doi.org/10.3390/act9020031

APA Style

Zhao, J. -Z., Yao, G. -F., Liu, R. -Y., Zhu, Y. -C., Gao, K. -Y., & Wang, M. (2020). Interval Analysis of the Eigenvalues of Closed-Loop Control Systems with Uncertain Parameters. Actuators, 9(2), 31. https://doi.org/10.3390/act9020031

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop