Active Vibration Control of PID Based on Receptance Method

Active vibration control approaches have been widely applied on improving reliability of robotic systems. For linear vibratory systems, the vibration features can be altered by modifying poles and zeros. To realize the arbitrary assignment of the closedloop system poles and zeros of a linear vibratory system, in this paper, an active PID input feedback vibration control method is proposed based on the receptance method. The establishment and verification of the proposed method are demonstrated. The assignable poles during feedback control are calculated and attached with importance to expand the application of the integral control. Numerical simulations are conducted to verify the validity of the proposed method in terms of the assignment of closed-loop poles, zeros, and both. The results indicate that the proposed method can be used to realize the active vibration control of closed-loop system and obtain the desired damping ratio, modal frequency, and dynamic response.


Introduction
Nowadays, the demand for lower vibration and higher reliability in robotic systems is increasing [1]. Extensive studies have been conducted to optimize the controlling of various nonlinear systems [2][3][4][5][6]. A common vibration control approach is to modify system dynamics using controllers [7]. PID-based approaches have been widely applied on motion control of robotic systems [8,9]. Generally, the vibration characteristics of a linear vibratory system are determined by its zeros and poles, and the assignment of eigenvalues in a system is of great importance in improving system reliability.
Traditionally, the solutions of pole assignment are solved by using the information of mass, damping, and stiffness matrices acquired from the finite element (FE) models [10]. However, FE models are generally limited by the incapability in obtaining precise damping models of real structures. Ram and Mottershead [11] developed an active vibration control method for linear system based on receptance method to accurately assign the system poles and zeros to specified values. One considerable advantage of this method is that it is entirely based on data from modal testing rather than the exact mass, damping, and stiffness matrices. Mottershead et al. assigned part of the system poles to the predetermined values while keeping nonimportant poles unchanged by making the nonimportant poles uncontrollable or unobservable [12]. Zhang et al. presented the theory of applying acceleration feedback and position feedback to the nondamping system to assign part of poles [13]. Ouyang proposed a method of using acceleration feedback and velocity feedback in active vibration control [14]. Recently, neural network-(NN-) based pole-zero assignment methods have gained a lot of attention [15][16][17][18]. However, NN-based approaches may fail to yield desired results on practice due to the black-box nature of NN and the requirement of NN parameter optimization skills.
Previous studies indicate that various combinations of acceleration, velocity, and position make it flexible and applicable to select the type of feedback when implementing active vibration control. However, a vibration system will inevitably become a singular system when the acceleration feedback is applied, which leads to the unavailability of the receptance method. In this paper, integral control is introduced to overcome this limitation. The contributions of this paper are summarized as follows: (1) A PID input feedback active vibration control method is proposed based on the receptance method. (2) The numerical solutions for assignment of closed-loop poles and zeros by applying the proposed method are demonstrated.
The rest of the paper is organized as follows: The assignments of poles and zeros in PID control systems are demonstrated in Section 2 and Section 3, respectively. Then, numerical simulations are conducted to verify the validity of the proposed method in Section 4. Finally, Section 5 concludes the paper.

Pole Assignment of PID Control Systems
2.1. Formulation of Pole Assignment. The equation of a n degree freedom linear system can be described in secondorder form as follows according to [19]: where M, C, K ∈ R n×n are mass, damping, and stiffness matrices, respectively; xðtÞ ∈ R n×1 is the displacement vector. M = M T , C = C T , K = K T , and for any nonzero vector ν ∈ R n×1 , there is ν T Mν ≥ 0, ν T Cν ≥ 0, and ν T Kν ≥ 0.
Integrate PID feedback control to the system and the closed-loop system can be written as follows: where b ∈ R n×1 is the force distribution vector and uðtÞ is the input control vector: where g 1 , g 2 , g 3 ∈ R n×1 are input control gain vectors for displacement, differential, and velocity, respectively. Express Equation (2) in the Laplace space: Then, Equation (5) shows that the highest order of the system stiffness matrix is a positive definite matrix, so the closedloop system must be a nonsingular system. Therefore, the system receptance matrix exists.
The receptance matrix can be obtained by applying the Sherman-Morrison formula for Equation (5): where HðsÞ = ½Ms 2 + Cs + K −1 is the receptance matrix of the open-loop system, which can be obtained in practice by mea-suring the receptance HðiωÞ [20]. The characteristic polynomial pðsÞ can be derived from Equation (6) at the same time: Thus, the roots of Equation (8) are the pole of the closedloop system: Given HðsÞ, b, and Equation (9): The problem of assigning the poles of a closed-loop system to setpoints fμ 1 , μ 2 , ⋯, μ N g can be transformed to the solution of input control gain vectors g 1 , g 2 , and g 3 .
is closed under conjugation, N is the number of closed-loop system poles that can be assigned.

Minimum Number of Assignable Poles N.
Given a positive definite system with closed-loop conjugate poles fμ i g N i=1 , suppose that the adjoint matrix of HðsÞ is AðsÞ and the determinant is DðsÞ. pðsÞ = 0 can be expressed as follows: Since for any μ i ∈ fμ i g N i=1 , Dðμ i Þ ≠ 0, Equation (10) can be written as follows: The highest degree of s in qðsÞ is ð2n − 2Þ + 2 = 2n and the highest degree of s in DðsÞ is 2n. Since the number of closedloop poles N is equal to the degree of Equation (10) and the highest degree of s in Equation (10) is 2n + 1, the number of characteristic polynomial root is N = 2n + 1. It should be noted that due to the introduction of integral control, the closed-loop system will produce an extra nonzero pole.
In fact, the introduction of integral control into closedloop system will lead to a nonregular and generalized system. The closed-loop system is stable if (1) it is regular and (2) all of its assigned poles have negative real components [21]. The regularity of the closed-loop system in this method is guaranteed by the positive definiteness of the mass matrix. The values of the assigned poles can be arranged to stabilize the closed-loop system given that their real components are negative and Journal of Sensors 2.3. Calculate Vectors g 1 , g 2 , g 3 and Minimum Sensor Layout Method. The solution of g 1 , g 2 , g 3 can be realized by the following equation: where , and G ∈ R N×3n , g ∈ R 3n×1 , and γ ∈ R N×1 . The pole assignment problem raised by Equation (2) can be solved by applying Equation (13) and Equation (14). One solution for Equation (13) is as follows: where G + is the Moore-Penrose generalized inverse of G.
When the assigned closed-loop poles are self-conjugate, g 1 , g 2 , and g 3 can be calculated according to Equations (13)- (15). Considering the limitation of the additional mass in the actual control system, the number of sensors used to measure the motion parameters of the system should be as small as possible. Therefore, the value of some items in g 1 , g 2 , g 3 can be assigned to zero by setting the values of the corresponding columns of the G to be zero. In practice, the solution means that fewer sensors are involved. Due to the integration of PID input feedback control and N = 2n + 1, the number that can be allocated to zero in G has at most 3 n − N = n − 1 columns. Therefore, the number of sensors can be reduced accordingly in practice to minimize the effect of the additional mass on the system.

Validity Analysis.
When the value of the preset poles fμ i g N i=1 and the optional k ð0 ≤ k ≤ n − 1Þ column vector to be replaced by the 0 in the G are determined, a specific PID input feedback solution is derived, identified as W k . From Equations (13)-(15), g 1 , g 2 , g 3 can be calculated in this specific PID input feedback solution. Below we will show how to verify whether W k is valid or not.
The closed-loop receptance matrixĤðsÞ can be derived from Equation (5) as follows: where ΩðsÞ is as follows: The problem of verifying whether a particular W k is valid is equivalent to verifying whether ΩðsÞ is a nonsingular matrix for any s ∉ fμ i g N i=1 under the determined case. The verification process of W k is as follows: Firstly, the real vectors g 1 , g 2 , g 3 are calculated according to Equations (13)- (15). Then, substitute g 1 , g 2 , g 3 into Equation (17) and obtain the determinant wðsÞ of ΩðsÞ. For any s ∉ fμ i g N i=1 , if wðsÞ is zero, the closedloop system poles cannot be assigned to the designated values, meaning that W k is invalid. Otherwise, the closed-loop system poles can be assigned to the specified desired values, indicating that W k is valid.

Zero Assignment of PID Control Systems
3.1. Formulation of Zero Assignment. The characteristic polynomial for the zeros of receptanceĤ i j can be realized when the denominator matrix (i, j) of Equation (6) is zero: where e i , e j are the unit vectors obtained from the corresponding columns of the identity matrix. Therefore, the problem of assigning zeros of the closedloop system receptanceĤ i j to the setpoints fξ 1 , ξ 1 , ⋯, ξ r g can be realized as follows: Given HðsÞ, b, i, j, and Calculate the input control gain vectors g 1 , g 2 , g 3 . Where ξ k ∉ fξ i g r i=1 , r ≤ 2n − 1, r is the number of poles in a closedloop system that can be assigned.
Equation (19) can be written in the following form: where H i j ðsÞ is the open-loop receptance value at the open-loop system coordinate of (i, j). Transform Equation

3
Journal of Sensors (20) to We have Therefore, for preset closed-loop system zeros fξ 1 , ξ 1 , ⋯, ξ r g, we have where It can be seen from Equation (22) that the inverse of the matrix T can be obtained when the preset zeros fξ 1 , ξ 1 , ⋯, ξ r g are given, and then the real vectors g 1 , g 2 , g 3 can also be obtained. In particular, when i and j are overlapped, the distribution of zeros will affect the closed-loop pole assignment of the system, which is significant [22].

Pole and Zero Assignment.
We can assign the zeros and poles in the closed-loop system at the same time according to Equation (13) and Equation (22). In this case, the total number of zeros and poles follows r + N ≤ 2n + 1. Similarly, by assigning k ð0 ≤ k ≤ n − 1Þ columns in Equation (13) and Equation (22) to zero, we can get the optimized allocation scheme.
Obtain g 1 , g 2 , g 3 according to Equations (13) The above PID input feedback solution is valid according to Section 2.4. In order to visualize the effect of this method, the frequency response curve, zero-pole distribution, and Nyquist plot of the H 3,3 of original system and the modified system are given, shown in Figures 1-5. Figures 2 and 3 show that the poles of the modified system are assigned to the specified positions, and an additional pole is added to the modi-  Journal of Sensors reassignment, system features such characteristic frequency, damping, and stability have been changed, and the overall performances have been improved.

Minimum Sensor Configuration for Pole Assignment.
Based on the pole assignment in Section 4.1 and the minimum sensor layout in Section 2.3, a minimum sensor configuration scheme for pole assignment simulation is developed and verified according to Section 2.4. Table 1 shows all solutions that can assign the poles of the closed-loop system to preset values μ 1;2 = −1 ± 0:5i, μ 3;4 = −1 ± i, and μ 5;6 = −2 ± i, while keeping the closed-loop system stable at the same time for the case in Section 4.1. The PID feedback gains for the 2 solutions are given in Table 2 (only the control gains of the first method in each W k are given). The system performance may destabilize when the number of assigned poles is less than N (i.e., when the lacked poles are assigned, the number of assigned poles in Table 3 is 2n). This highlights the importance of determining the number of assignments. The specific results are demonstrated in Tables 1-3.       Considering the threedegree-of-freedom system in Section 4.1, we use zero assignment method in Section 3 to modify the closed-loop system zeros and the closed-loop response of the system. Set b = 1 1 1 ½ T and the closed-loop system H 3,3 ðsÞ zeros as ξ 1;2 = −1 ± 0:5i, ξ 3;4 = −2 ± 0:5i. Obtain g 1 , g 2 , g 3 according to Equation (22) Figure 4, Figure 6, and Figure 8 indicate that, after the assignment of zeros, the system stability increases.

Pole and Zero Assignment
Simulation. Based on the threedegree-of-freedom system in Section 4.1, set b = 1 1 1 ½ T , the closed-loop system H 3,3 ðsÞ zeros as ξ 1;2 = −1 ± 0:5i, and the closed-loop system poles as μ 1;2 = −1 ± 0:5i. Obtain g 1 , g 2 , g 3 according to Equations (13)- (15) and Equation (22): The closed-loop poles calculated at this time are −0:3563 ± 1:6213i, −0:9778 ± 0:9769i, −1 ± 0:5i, and −3. The poles are distributed in the left half plane of the imaginary axis and the closed-loop system is stable. The frequency response curve, the zero-pole distribution diagram, and the Nyquist diagram of the modified system H 3,3 are illustrated in Figure 9, Figure 10, and Figure 11, respectively. Figure 2 and Figure 10 indicate that the system zeros and poles are accurately assigned to the desired values. After the assignment, an increased system stability is observed according to Figure 4, Figure 9, and Figure 11. The validity of the proposed method is verified by the simulation results.

Conclusions
In this paper, a PID-active vibration control method is proposed based on the receptance method. The proposed method can be utilized to assign the poles and zeros of closed-loop systems in order to improve the performances. Due to the integration of integral control, the proposed method guarantees the positive definiteness of the system and avoids the appearance of a nonsingular system. The proposed method also extends the application of the dynamic flexibility method. The procedures of assigning poles and zeros to a system are demonstrated. Numerical simulations for pole and zero assignments are conducted. Simulation Table 3: Results under lacked poles.

Preset poles
Control gain Actual poles    Journal of Sensors results indicate that the proposed method is generally effective. The proposed method in this paper can be utilized to improve performance and reliability on robotic systems such as flexible robots, robot arms, and intelligent vehicles. In future work, the proposed method will be expanded to handle vibration control tasks for nonlinear systems with multi inputs and outputs.

Data Availability
No data was used to support this study.

Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.    Pole-zero map Figure 10: The pole-zero map for the modified system H 3,3 .