Abstract
The design problem of a discrete controller with time delay and acceleration feedback for a single-link flexible manipulator system is addressed in this paper. The dynamical model of a single-link flexible manipulator system is presented by the adoption of the finite element method and Lagrange’s equation. Based on the random-walk process and the discrete reduction method, an augmented discretized delay-free state derivate space equation containing the random noise is established. An acceleration-based Kalman filtering method is developed in order to estimate the system state and external excitation necessary for the controller design. In light of the estimated augmented states, a hybrid controller that combines a feedback control algorithm and a feedforward control algorithm is designed according to optimal control theory and Moore–Penrose theory. Numerical simulation results show that the proposed controller can damp out the vibration response of the flexible manipulator system effectively upon external excitations. Moreover, it is further revealed that the control performance of the presented method can be improved by adding the time delay appropriately.
1. Introduction
Over the past few decades, manipulators constructed with flexible arms have attracted much more attention because flexible manipulators have the advantages of lighter weight, lower energy consumption, and wider operating range as compared to the traditional rigid manipulators. However, these advantages are generally obtained at the expense of the stiffness of the links. Due to the complexity of link deformation and the rigid-flexible coupling effect [1], establishing an accurate dynamic model and suppressing the unwanted vibration have become a key research challenge in practical design. Comprehensive literature reviews on the dynamic analysis and control schemes of manipulators with flexible links or flexible joints are summarized in a number of survey papers [2–5]. Nevertheless, vibration control is still a key issue in flexible manipulators, especially when the controlled system suffers from a complicated environment.
Various control schemes can be applied to discrete flexible manipulators after the implementation of a series of truncating techniques, such as the finite element method (FEM) [6] or the assumed mode method (AMM) [7]. Based on the FEM and Newmark’s solution, Karagulle et al. [8] developed a Matlab code for a two-link manipulator with flexible members and simulated end point vibrations of the system. They designed trapezoidal-velocity-profile inputs to minimize the residual vibrations. Meng et al. [9] presented a stable control strategy to suppress the vibration of a single-link flexible manipulator built by using AMM. In their work, the performance of the system was improved by optimizing the design parameters of the controller using an online fuzzy genetic optimization algorithm. When the flexible manipulators conduct a space operation, the robust controllers are required to cope up with the unavoidable existence of parameter uncertainties and unknown bounded external interferences. Sayahkarajy et al. [10] proposed a composite controller to control the joint angles and vibrations of a two-link manipulator with elastic links, which uses a command preshaping filter as a compensator and uses a multiple-input and multiple-output H∞ loop-shaping controller as the feedback controller. Yang et al. [11] developed the closed-form dynamics of a two-link flexible manipulator under unknown disturbances by using the AMM and Lagrangian methods. A hybrid controller was then designed based on the singular perturbation theory in which an adaptive sliding control was designed for the slow subsystem. Steve et al. [12] proposed an adaptive trajectory control of space flexible joint manipulators under parametric uncertainties and modeling errors. Friction, inertial coupling, and nonlinear characteristics are considered in their dynamic model. Luo et al. [13] designed an active vibration control of flexible manipulators by combining input shaping technique (IS) and an autodisturbance rejection controller (ADRC).
The Kalman filter is a powerful tool to estimate the states of a state-space plant model having the presence of modeling uncertainties or the measurement noise. A linear quadratic Gaussian (LQG) controller associated with the LQR and Kalman–Bucy filter (KBF) is adopted to control the vibration of flexible manipulators by Castri and Messina [14]. They demonstrated that the control spillover and the observation spiller could be reduced by an appropriate design of the regulator and the KBF. Kumar [15] used a positive position feedback (PPF) controller at piezoelectric (PZT) actuators to minimize the end effector vibrations of the flexible link. The negative effect of hardware filters was compensated by a modified Kalman filter. Moreover, in previous works, the accelerometer is generally used to achieve a high vibration feedback control efficiency for the flexible manipulators with good accuracy. The measured acceleration signals can be directly adopted for feedback control to avoid direct current component and trend terms induced by integral actions [16, 17]. As a result, it is desirable to design an acceleration-based Kalman filter in a straightforward manner for the state observation and control implementation.
In practice, low-frequency vibration is one of the main factors that affect the motion accuracy of space flexible manipulators. From the viewpoint of signal processing, it is necessary to introduce a lower filter in the control structure to eliminate the influence of measurement noise and high frequency signals. As a result, the inevitable time delay induced by the lower filter must be taken into account in the modeling and control of flexible manipulators. Chen et al. [18] was the first to introduce time delay caused by the friction into flexible joint robotic manipulators. Xu et al. [19–21] conducted a series of studies on the stability and vibration of a flexible joint with transmission delay. Li and Yin [22] investigated the adaptive position-force control for multiple flexible joint coordinated robot, where the time delays were caused by gear clearance. Although the time delays will bring undesirable control performance or even lead to instability to the system, many researchers have utilized it as a control parameter to improve the efficiency of the controller. A type of time-delayed optimal controller was designed by Qiu and Wu [23] to compensate for the unknown time delay in a rotating flexible manipulator. A computed torque-time delay approach was presented in Reference [24] to damp out the end point vibration of a single-link flexible manipulator. According to the measured deflection point, the delay time was chosen to get a fast decay of the end point vibration. Saeed et al. [25, 26] proposed a time-delayed position-velocity feedback controller to suppress the vibration of the nonlinear system. They made a clarification of the conditions that the time-delayed controller works as a damper or an exciter. An et al. [27, 28] designed a time-delayed acceleration feedback controller for a smart cantilever beam. The proposed controller not only changed the mass properties of the controlled system but also improved the damping coefficient of the system.
In the present study, a time-delayed acceleration feedback controller considering the effects of the external disturbance is applied to a single-link flexible manipulator. The dynamic model of the control system considering input delay and random excitation was first established by the use of the FEM and Lagrange’s equation. The corresponding continuous system is transformed into its discrete delay-free form through a peculiar integral transformation. An augmentation discrete-state derivate space model including system states and external force is constructed, where the unknown external disturbance can be seen to be a random-walk process. An acceleration-based feedback Kalman filtering equation is then developed to estimate the system state and external force, which are necessary for the controller design. In light of the estimated augmented states, a hybrid controller that combines the feedback control algorithm and the feedforward control algorithm is designed by the use of optimal control theory and Moore–Penrose theory. Numerical simulations are performed to demonstrate the ability of the proposed controller in reducing the vibration response.
2. Dynamic Modeling
2.1. Equation of Motion
A schematic diagram of a single-link flexible manipulator system is shown in Figure 1. The entire system is driven by an electric motor, and a pair of PZT actuators is embedded in the root segment of the flexible link. The mounting base is considered to be a part of the flexible link in the present study. An accelerometer is attached to the distal part of the flexible link, which can be regarded as a payload with mass and inertia.

The link can be maneuvered in the horizontal plane by a revolute joint with one degree of freedom. The rigid motion and the flexible deformations of the link occur in the same plane. The FEM is applied to discretize the flexible link, which satisfies the complicated cross-sectional and boundary conditions. Euler–Bernoulli beam theory is applied to the modeling of flexible links by assuming that the link is long and slender. The coordinate systems are used to develop the dynamic model, as shown in Figure 1. Here, OXY is the inertial reference frame, and its origin is located at the driving hub. O1X1Y1 is the local coordinate system attached to the rigid hub, where O1X1 is the tangent to the undeformed link. After dividing the flexible link into ne elements, the transverse deflection at position x in the ith (i ≤ ne) element with respect to its local coordinate system can be expressed in terms of nodal displacement multiplied by shape functions aswhere is the shape function of the ith element, denotes the deflection and the flexural rotation at nodes i and i + 1, respectively. The displacement of point x with respect to the inertial frame can be expressed aswhere T is the rotational transform matrix from the local frame to the inertial frame, and is the displacement of point x at the ith element in the local frame.
Then, the kinetic energy of the flexible link can be written aswhere li is the length of the ith element, and are the density and cross-sectional area of the ith element, respectively.
In addition, the total system kinetic energy includes the kinetic energy of the rotors,and the kinetic energy of the accelerometer,where Jm is the rotary inertial of the hub; mp and Jp are the mass and moment of inertia of the payload, respectively; is the angular position of the link; and is the nodal rotation of the neth element in the flexible link.
The potential energy is generated only by the elastic deformation of the flexible link in the following form:where Ei and Ii are Young’s modulus and inertial moment of the ith element in the flexible link, respectively.
The work done by the nonconservative external torque is given by
2.2. Dynamics of PZT Actuators
Two PZT actuators are embedded symmetrically in the mounting base, which is considered as a part of the flexible link. The cross sections are shown in Figure 2. By assuming that p1 and p2 are the nodal numbers of the flexible link on which the actuators act, the bending moments produced by the PZT actuators can be expressed aswhere u is an externally driven voltage; ; ; and are the deflections at nodes p1 and p2, respectively; Ep, Ap, and lp are the Young modulus, cross-sectional area, and length of the piezoelectric actuators, respectively; hp is the distance from the controlling point to the neutral axis of the flexible link; and is the constant describing the relationship between the control voltage and the deformation of the actuator.

2.3. Dynamics of the Flexible Manipulator
Substituting the expressions for kinetic and potential energy into Lagrange’s equations and performing the necessary operations by ignoring the nonlinear terms [29], one obtains the following equations:wherewhere is the equivalent moment of inertia of the flexible link; is the Boolean operation matrix of the ith element; is the control voltage input matrix; and is the control delay.
It should be noted that equation (9) describes the rigid motion of the flexible link, while equation (10) represents the vibration of the flexible link subject to control force. Substituting equations (9) into (10) and considering the structural damping, the vibration equation of the flexible link subject to motor torque can be expressed aswhere , , and ; c is the structural damping matrix, which has the form of the assumed Rayleigh damping; and the corresponding proportional damping coefficient can be determined by a modal test, as done in our previous study [30].
Then, the state-space equation of the flexible manipulator reads
where , , , .
The acceleration at the tip of link can be expressed aswhere is the output acceleration matrix and n = ne.
2.4. Discrete Delay-Free System
Both the continuous state equation (12) and the measurement (14) are discretized by using a time interval ∆T. Furthermore, the measurement noise is added to the discretized observation equations by taking into consideration the uncertainties and disturbances in the practical environment. The discrete time forms of system can be written aswhere t = k∆T, , , and are the kth step of the discrete vector , , , and , respectively; the measurement noise vector is assumed to be white noise with zero mean, and covariance is ; is the relationship between the time delay and the sampling interval, where h is a positive integer; and
A peculiar transformation, as done in common practice [31], is introduced to transform the time-delay system into one which contains no time delay apparently. Let
Substituting (16) into (15), the discrete delay-free system can be derived aswhere .
The random-walk process has been widely used in estimation and identification of unknown time-varying system parameters [32]. The unknown input of the system can be represented to be a random-walk process as follows:where is a zero-mean white Gaussian random process with covariance .
Substituting (19) into (18) yields the following extended discrete delay-free system:where
3. Acceleration-Based Kalman Filtering Equation
In the previous section, the external force is regarded as an auxiliary state into the delay-free state system by using random-walk processes. Therefore, the Kalman filtering method can be used directly to estimate the system state and external force at the same time.
According to Reference [33], the extended discrete delay-free system is a linear deterministic-stochastic one, which can be decomposed into the following linear deterministic system:and a linear stochastic system,
with and . It is easy to obtain the solution of in the linear deterministic system by using a transition equation. Hence, it is sufficient to derive the optimal estimate of in the stochastic state space (20), so that
becomes the optimal estimate of the state vector in the original system (18).
Assume that a priori estimate of the state vector at time k +1 satisfies
Then, the a priori estimate of the measurement can be written as
The residual measurement can be obtained according to (23) and (26) as follows:
Then, the posteriori state estimate can be corrected using the following expression:where is the Kalman gain matrix.
The error between the actual state and its prediction and the error between the actual state and the updated state estimate are defined at time k + 1, respectively, as
Substituting equation (28) into equation (29) and combining equations (23) and (25), the following state estimation error is obtained:where equation (30) is a recursive equation for the state estimation error. Choosing the covariance of as the optimal objective function and expanding the terms in the equation above, we can express the updated estimate covariance as follows:
When the derivative of estimate covariance with respect to the Kalman gain matrix equals zero, the trace is minimized. Therefore, the optimal Kalman gain is derived as follows:
In summary, the presented acceleration-based Kalman filtering algorithm can be determined by using the following equations:
4. Control Design
In relation to the random-walk process, the system state and input force can be estimated simultaneously by using the Kalman filtering method. Hence, a hybrid controller using the estimated system state and input disturbance force is designed as follows [34]:where is the feedback gain with respect to the state estimate and is used for suppressing the transient vibration, whereas is the feedforward gain with respect to the external force estimate and is used for damping the stationary vibration.
Substituting equation (32) into equation (18), we have
The regular gain can be obtained by using the optimal control method when is controllable. The optimal performance index is defined as follows:where and are all symmetric weighting matrices. Then, the feedback gain is given bywhere P is the solution of the Riccati equation,
In order to attenuate the effect of the input force, is designed to satisfy the convergence condition , according to the second term on the right of (34). There exists the only Moore–Penrose inverse since matrix has the same size of matrix Buz, which is a row full-rank matrix. According to the theory of Moore–Penrose pseudoinverse, can be obtained from the following equation:
In addition, is a least squares solution of minimum norm of so that
Finally, the control law with an input delay can be written as
The schematic diagram of the time-delayed acceleration feedback control based on the Kalman filter is shown in Figure 3.

5. Numerical Simulations
5.1. Parameters
The applicability of the proposed control strategy has been examined by simulation on a rotating flexible manipulator with variable cross sections, which is simplified from an experimental model with reasonable accuracy [30], as shown in Figure 4. The parameters of the PZT actuator and accelerometer are given in Table 1. According to the parameters in Table 2, the flexible link has been finely meshed to ensure the reliability of the finite element division. In relation to the cross-sectional characteristics, the follow-up unit and the mounting base of the actuators can be divided into five elements, as shown in Figure 5.


In order to achieve higher effective control with a lower computational burden, an iterated improved reduced (IIRS) method is applied [30]. The original model is divided into 20 elements with the dimension of full-order model being 40. By using the IIRS theory, PZT acting nodes and tip nodes are chosen as master coordinates, while the dimension of the reduced model is finalized to be 4. The first three eigenvalues and impulse response of the reduced system are compared with that of the original one, as illustrated in Table 3 and Figures 6 and 7, respectively. It can be seen from them that the reduced model has a good agreement with the original model. In addition, the corresponding computational time is decreased from 33 s to 1 s.


5.2. Effect of Time Delay on Controller
To investigate the effect of the time delay on the proposed controller, the controlled system with zero initial conditions is analyzed upon harmonic excitation containing random white noise. The time delay λ is chosen to range from 0 to 6∆T, as listed in Table 4, with the state weighting matrix Q = Q0 × I6 × 6, where Q0 = 108, and the control weighting R = 1. The root mean square (RMS) ratio between the acceleration response of the controlled system ac and the counterpart of the uncontrolled system ao is shown in Figure 8. With the increase in λ, the RMS ratio is decreased from 10% to 4%, and it is not nearly altered when λ is bigger than 3∆T. Hence, a good control performance can be achieved when λ is set to be 3∆T. It can be found from Figure 9 that the vibration acceleration is further attenuated when the time delay is actively introduced in the proposed controller. These results show that the active utilization of time delay will become a favorable factor in vibration suppression, and the control performance will be improved as the time delay is chosen within a reasonable range [26].


5.3. Effect of Weighting Matrix on Controller
To study the influence of the weighting matrices on the controller, the state weighting Q0 is first selected to range from 100 to 1012 with the control weighting R and time delay λ being fixed. The RMS ratio of the acceleration is reduced to 4% with the increase of Q0 to 108, and it is not nearly altered when Q0 is bigger than 108, as shown in Figure 10. Therefore, a good control performance can be obtained when Q0 is set to be 108. Subsequently, by changing the control weighting R in the range from 0 to 30 with Q0 = 108 and λ = 3∆T, the influence of control weighting on the control performance is evaluated, as illustrated in Figure 11. It can be found that the RMS ratio will exhibit larger values with the increase in the control weighting R, and the minimum one is achieved when R = 1, subject to the physical limitation. Finally, a good combination of the state weighting matrix and the control weighting one, which leads to a low RMS ratio, can be found in a 3D colormap surface, as shown in Figure 12.



5.4. Simulation Cases
In order to demonstrate the effectiveness of the proposed method, the controller presented in Section 4 was implemented and verified on a single-link flexible manipulator containing a PZT actuator (as shown in Figure 4) for vibration reduction at its tip node upon different excitations, as listed in Table 5.
A multifrequency torque excitation is first taken into account, which contains the first two resonance frequencies that dominate the vibration response of the structure. By using the proposed acceleration-based Kalman filtering method, the torque excitation along with the system state has been estimated accurately, as indicated in Figures 13 and 14, respectively. The control behavior of the proposed hybrid controller is then assessed by simulation, and the vibration response of the flexible manipulator can be effectively suppressed, as validated in Figure 15. By comparing the result of the hybrid control with that of the control without compensating for an external torque, it can be observed that the former in control performance is clearly superior to the latter, where the RMS ratio of acceleration can be decayed from 10% to 3.8%, as shown in the subplot of Figure 15.



Subsequently, the torque excitation with a random white noise is considered to verify the robustness of the proposed controller. As mentioned in Section 3, the state and external force are estimated through a Kalman filter being designed in the state derivate space, and as a result, the robustness of the control system is guaranteed. It is easy to observe that the estimated state and external torque with the noise are in reasonable agreement with the actual counterparts, as shown in Figures 16 and 17, respectively. In turn, the controller can keep the vibration response at tip node converging to a low level, as validated in Figure 18, and the stability of the control system is enhanced under uncertain random disturbance.



In final, a quintic polynomial trajectory, which is proved to be smoother for tracking purpose [35], is chosen as the external motion excitation of the flexible link in the simulation, as depicted in Figure 19. The vibration response of the flexible manipulator subject to this excitation can be further reduced by the use of the proposed controller, as shown in Figure 20.


6. Conclusion
In this paper, a time-delayed acceleration feedback controller for reducing the vibration of a single-link flexible manipulator is presented. An augmented discretized delay-free state derivate space system is established by employing a discrete transformation and random-walk process. A hybrid controller that uses a combination of the feedback signal and feedforward one is designed to cancel out external excitation and damp out vibration at the same time. Compared with the conventional LQG method, a discrete acceleration-based Kalman filtering method is developed to estimate the system states and external excitation with good accuracy for the improved control behavior, regardless of the effect of white noise. Moreover, it can be further confirmed from simulations that the introduction of time delay will become a favorable factor in vibration suppression, and the control performance will be improved as the time delay is chosen within a reasonable range. The simulation results demonstrated the effectiveness of the proposed control scheme. This study provides a basis for validating its practical feasibility, which is an open research topic.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.