Research Article  Open Access
Workspace Tracking Control of TwoFlexibleLink Manipulator Using Distributed Control Strategy
Abstract
A distributed nonlinear control strategy for twoflexiblelink manipulators is presented to track a desired trajectory in the robotâ€™s workspace. The inverse dynamics problem is solved by transforming the desired trajectory from the workspace to the joint space using an intermediate space, called virtual space, and then using the quasistatic approach. To solve the nonminimum phase problem, an output redefinition technique is used. This output consists of the motorâ€™s angle augmented with a weighted value of the linkâ€™s extremity. The distributed control strategy consists in controlling the last link by assuming that the first link is stable and follows its desired trajectories. The control law is developed to stabilize the error dynamics and to guarantee bounded internal dynamics such that the new output is as close as possible to the tip. The weighted parameter defining the noncollocated output is then selected. The same procedure is applied to control and stabilize the first link. The asymptotical stability is proved using Lyapunov theory. This algorithm is applied to a twoflexiblelink manipulator in the horizontal plane, and simulations showed a good tracking of the desired trajectory in the workspace.
1. Introduction
Many control strategies for manipulators have been the focus of several studies in recent years. The robot manipulators consist of a sequence of links and joints in various combinations. In industrial applications, most of the existing manipulators use rigid links and joints and are known as rigid manipulators. Rigid manipulators are generally slow, extremely rigid, and massive, and the useful load is very low compared to their weight. To improve the performance of the robot manipulators, their links must be lighter, and therefore they become more flexible. Flexible manipulators present more advantages when compared to rigid manipulators: they are faster and less massive and consume less energy. Some flexible manipulators are used in different areas, for example, the aerospace applications [1] and medical applications [2].
For flexible manipulators, the problem of workspace tracking trajectory is less covered so far than that of joint space tracking. There are few solutions to the workspace tracking problem, particularly for manipulators with many flexible links. The workspace tracking trajectory is very important since most of the tasks are defined in the operational space, such as painting, welding, and assembly. The flexible link manipulators are a nonminimum phase system when controlling the position of the end effector [3]. Unlike rigid manipulators, the inverse kinematics of flexible manipulators is not sufficient to transform the desired trajectory from workspace to joint space because they are linked by kinematics and dynamics relationships. To solve this problem, many studies focused on the output redefinition technique [4â€“7]. This technique consists in selecting a new output as close as possible to the tip such that internal dynamics become bounded. In [5], this approach was used for a class of manipulators where the last link is flexible. An intermediate space between the joint space and the workspace called virtual space is used to transform the desired trajectory from the workspace to the joint space. The virtual space is linked with the workspace by a simple kinematics relation as in rigid manipulators. In [4], the output redefinition technique was used for one flexible link manipulator. This output consists of the motorâ€™s angle augmented with a weighted value of the angle of links extremity.
The nonlinear dynamics of flexible link manipulators combined with their underactuated nature (the deflection variables are not actuated) present a challenging control problem. Multiflexiblelink manipulators can be controlled as one MIMO system so a single controller is used for all joints and links or as a set of interconnected subsystems so each pair of joint and link is controlled by its own controller. For the first case, many control schemes were used. Several studies used linearization around a nominal configuration of flexible manipulators model [8, 9]. Some nonlinear effects such as the variation of inertia matrix around an operating point were taken into account in the control design methodology. In [10], a control method based on the inputoutput linearization was developed to track a desired trajectory in the workspace for a class of flexible link manipulators when the last link is flexible. To select an output near the tip that guarantees the stability of the zero dynamics, authors used the output redefinition technique. A robust adaptive controller was developed in [11] for a class of flexible link manipulators where the last link is flexible. This controller is based on feedback linearization and uses the virtual joint space that is kinematically related to workspace. For trajectory tracking control problems, many other techniques were covered in the literature such as singular perturbation technique [12], sliding mode control [13], inversionbased nonlinear control [14], and so forth. All the abovementioned control methods use a single controller for all joints and links as one MIMO system. Unfortunately, due to the complexity of the control structures, the realtime implementation in industrial applications is not easy [15]. A solution to this problem can be followed by considering the dynamics of the robot manipulators as interconnected subsystems (joint and link). Many control schemes consider this configuration. In [16], a decentralized control method for flexible link manipulators was used. The authors used a simple proportional derivative (PD) controller for the joints and a linear quadratic regulator (LQR) with output feedback for each link. In [17], a PD controller was used for each joint and the measurement of linear velocity of the tip position was used for controlling each link. A distributed control strategy was introduced in [18, 19] for rigid link manipulators. The distributed control strategy consists in controlling one joint at a time starting with the last joint and going backward until the first one. Lyapunov theory was used to prove the global stability of the error dynamics, and this controller was successfully applied on a 7 DOF manipulator. This control strategy was modified in [20] to take into account linksâ€™ flexibility. A good tracking performance in the workspace of a twoflexiblelink manipulator was obtained. However, the control law, based on the feedback linearization approach, ensures only local stability. In this paper, a nonlinear distributed control strategy is presented for a twoflexiblelink manipulator that ensures global stability. This strategy consists of controlling the second joint and link by assuming that the remaining joint and link are stable and follow their desired trajectories. Then, going backward, the first joint and link are controlled by following the same strategy. For the inversion dynamics problem, the virtual space and the quasistatic approach are used. The output redefinition technique is used to obtain the nearest point to the tip that ensures bounded internal dynamics. The Lyapunov approach is used to analyze stability of the tracking errors.
The paper is organized as follows. Section 2 presents the modeling of the twoflexiblelink manipulator and presents its main properties that will be used in the control law design. Section 3 presents the distributed control strategy. The stability analysis is given in Section 4. The control method is applied on a twoflexiblelink manipulator, and the simulations are given in Section 5. Finally, a conclusion is given in Section 6.
2. Modeling
2.1. System Description
Figure 1 shows a twoflexiblelink manipulator. This system moves in the horizontal plane and consists of two motors that generate two torques, two flexible links with mass , length , linear density , and rigidity , and a payload that has a mass . The first link is attached to the first motor, and the second link is clamped to the rotor of the second motor. The flexible links are supposed uniform and are modeled as EulerBernoulli beams, and the deformations are assumed to be small.
Using Lagrange equations, the dynamical model of an DOF flexible manipulator is given by [21] where is the inertia and mass matrix, is the Coriolis and centrifugal forces vector, is the friction matrix, and is the rigidity matrix. represents the vector of the generalized coordinates and is the vector of the applied torques. For the rigid coordinates and flexible links, the deformation of the th flexible link is given by the following equation: where is the th generalized flexible coordinate, is its th shape function, and is the number of the retained flexible modes of the th flexible link. Note that the total number of the flexible modes is and the number of the rigid modes is . The dynamical model (1) can be written as follows: where and are mass and inertia matrices for the rigid and flexible part, respectively. is a coupled element. is the stifness diagonal matrix, and is the damping diagonal matrix of the flexible part. The subscripts and denote the rigid and flexible modes.
2.2. Proprieties and Problem Formulation
The dynamical model of the flexible link manipulators has the following properties that will be used in the control law development:(P1)â€‰â€‰and are symmetric positive definite matrices [21].(P2)From (P1), we can deduce that the diagonal elements of are positive [22]: (P3)There exists a matrix such that (P4)The inertiamass matrix and the Coriolis matrix satisfy the following skewsymmetric property: (P5)The propriety (P4) is preserved for the diagonal elements of and . Then, we can write
The new noncollocated output of the th link defined by the motorâ€™s angle augmented with a weighted value of the endpoint angular position is given as follows: where , with being the length of the th flexible link, and is the shape function of link and mode . In this paper, we consider only the first flexible mode of each link (i.e., ). Thus, we have .
In the distributed control strategy, when controlling the th link, we assume that the other links follow their desired trajectories. For the th link, a new generalized coordinate is defined such that the th coordinate is the controlled one and the other links follow their desired trajectories such that
where () are the rigid and flexible modes associated with the th motor and link, respectively. The objective of this work is to track a desired trajectory defined in the workspace of a twoflexiblelink manipulator using distributed control. Three steps are followed to achieve this objective.
Step 1. Transform the desired trajectory from the workspace to the joint space using inverse kinematics and quasistatic approach.
Step 2. Develop the control law for the second and first links to stabilize the error dynamics and to guarantee bounded internal dynamics such that the output is as close as possible to the tip.
Step 3. Study the global stability.
3. Distributed Control Strategy
3.1. Inverse Dynamics
To achieve the objective of workspace tracking trajectories, we need to transform the desired trajectories from the workspace to the joint space. The flexible manipulator is a nonminimum phase system when the end effector is used as the output. In this case, kinematic and dynamic relationships link the workspace and the joint space. To overcome this problem, an intermediate space called virtual space can be used. Then the desired workspace trajectory is transformed to the virtual space using an inverse kinematics relation as in rigid manipulators. To transform the desired trajectories from the virtual space to the joint pace, the quasistatic approach can be used to solve a nonlinear equation for the flexible part.
Using inverse kinematics as in rigid manipulators, the generalized coordinates in virtual space can be easily found. The deformation is assumed to be small. According to Figure 2, we can write In the virtual space, the generalized coordinate is given as follows: where and is the shape function of link .
Using the Jacobian matrix as in rigid manipulators, the velocity and acceleration in the virtual space can be deduced. Then, for a twoDOF manipulator, the inverse kinematics is given by the following equation: where and are the workspace desired positions; ; and is the Jacobian matrix. ; ; ; ; is length of link 1 and is length of link 2. The rigid and flexible coordinates are derived via quasistatic approach and using the generalized coordinates in the virtual space.
The dynamical model (3) can be written as a function of the desired coordinates as follows: The generalized flexible and rigid coordinates in the joint space are deduced using the quasistatic approach. In the first step, the nonlinear equation of the flexible part (14) is solved to find the generalized flexible coordinates. The quasistatic approach neglects the desired velocity and acceleration of the flexible coordinates (). Then, the equation of the flexible part (14) can be written as follows: where Using (16), the generalized flexible coordinates are found by solving (15). Then, the generalized rigid coordinates are given by (16).
3.2. Control Strategy
For twoflexiblelink manipulator, the distributed control strategy consists of controlling and stabilizing the last joint and flexible link by assuming that the first joint and flexible link are stable and follow their desired trajectories. Then, we move backward and apply the same procedure to the first joint and link. For each step, a new noncollocated output and a control law are developed. The weighting parameter characterizing the noncollocated output is calculated such that the tracking error is asymptotically stable. Thus, the system becomes a minimum phase with the selected new weighted outputs.
In this paper, the twoflexiblelink model given in [21, 23] is modified by considering only the first flexible mode of each link [20]. Thus, we have â€‰ â€‰andâ€‰â€‰.
Using (P3), the dynamical model of the twoflexiblelink manipulator can be written as where
In rigid and flexible part decomposition, the dynamical model (17) can be written as where
To develop a control law, the dynamical model (17) is written as two interconnected subsystems. Each subsystem has a pair of joint and link. There exists a nonsingular matrix of transformation such as where is the original generalized coordinate, is the transformed one, and the matrix of transformation is given by Using the previous transformation, the dynamical model (17) can be written as follows: where and .
Equation (23) is equivalent to the following expression: where ; ; ; ; ; and .
The modified dynamical model (24) can be written as follows: where
Equation (25) can be written in the following form that will be used for control law: where ; ; ; ; ; ; ; .
The first and the second subsystems can be characterized by two mass and inertia symmetric positive definite matrices:
The new generalized coordinate, used for the second joint and link while the first joint and link are assumed stable, is given as follows: Note that the coordinates of the first subsystem are the desired ones, and the coordinates of the second subsystem are the controlled ones. Using (27), the equation of motion of the second joint and link can be written as The velocity is the time derivative of and the acceleration is the time derivative of .
According to (8), the new noncollocated output is given as follows: where when considering just one mode. is the weighted parameter characterizing the second noncollocated output.
Using (31), the equation of motion of the second subsystem is given as follows: where and .
The internal dynamics of the second link is deduced from (33) as follows: where ; ; ; ; ; ; ; , and is chosen such as .
Inserting (34) in (32), we get where ; ; ; ; ; ; ; and .
To control the second link, we propose the following control law: where is used for the rigid part, is for the flexible part, and is the setting term: where ;â€‰ â€‰ and ; .
Now going backward to the first subsystem and assuming that the second subsystem is stable, the new generalized coordinate associated to the second subsystem becomes Using (27), the equation of motion of the first subsystem (i.e., first link and joint) is given as follows: The new noncollocated output, defined as the angle of the first motor augmented by a weighted value of the linkâ€™s extremity angle, is given as follows: where for one mode.
Using (39) and (40), the equation of motion of the first joint and link is given as follows: Let us define and .
To find the internal dynamics of the first link, is given from (42) as follows: where ; ; ; ; ; ; ; , and is selected such as . Using the same strategy, the control law of the first subsystem is given as follows: where where ;â€‰ â€‰ and ; . , ; , ; ; . The distributed control strategy can be represented by Figure 3
4. Stability Analysis
This section presents the stability analysis of the tracking errors of the rigid and the flexible parts. First, we study the global stability of the rigid part by inserting the two control laws in the initial dynamical model. Second, the stability of the flexible part is studied to select the weighted parameters and guarantee bounded internal dynamics such that the new output is as close as possible to the tip.
4.1. Stability of Rigid Part
In compact form, the first (36) and second (44) control laws can be written as follows: where are given in .
To study the global stability, we insert the control law given in (46) in the dynamical model given in (1).
Proposition 1. The equation of motion of the twoflexiblelink manipulator (17) is equivalent to the following model: where ; is positive definite matrix; and ; and satisfy the following skewsymmetric property:
Proof. See Appendix A.
can be written as where is the diagonal elements and â€‰â€‰is the nondiagonal elements of . Then, the dynamical model (50) can be written as follows: Using (52), the control law (46) is equivalent to Now, the global stability is studied by inserting the control law (54) in the dynamical model (53). Then, the error dynamics are given by the following equation:
Proposition 2. When using the control law (54) and the dynamical model (53),â€‰â€‰the error dynamics given by (55) are equivalent to
Proof. See Appendix B.
To prove the asymptotical stability of the error dynamics, we propose the following Lyapunov function: The time derivative of is Using Proposition 1, we can conclude that , â€‰â€‰andâ€‰â€‰ are positive definite diagonal matrices. Then, using LaSalle theorem [24], the error dynamics are globally asymptotically stable.
4.2. Stability of the Flexible Part
In this section, asymptotical stability of the flexible part is studied using the internal dynamics. From the dynamical model given in (1) and (3), the equation of motion of the flexible part of the second link is given by Using the new generalized coordinate (29) and the new noncollocated output (31), the equation (60) is equivalent to The internal dynamics of the second link are given by subtractingâ€‰â€‰â€‰â€‰from (61) as follows: where ; ; ; ; ; , and is chosen such as .
The quasistatic approach, using the inverse dynamics, neglects . Then, the tracking errors of the rigid and flexible part of the second link are given by the following expressions: In the above section, the tracking error of the rigid part was proved to be asymptotically stable; then, as . Then, the tracking error of the flexible part can be written as In state space form, the tracking error of the flexible part takes the same form as the one obtained in [20] and is given by the following expression: Let ; .
The critical value of is selected such that is Hurwitz and the Nyquist diagram of is in the right half plan [20], where . and are given in (62).
The same procedure is now applied, proceeding backward, to the first jointlink subsystem. Using (1) and (3), the equation of motion of the flexible part of the first link is given by Using the new generalized coordinate (38) and the new noncollocated output (40), the equation of motion becomes Subtracting from (67), the internal dynamics are given as follows: where ; ; ; ; ; , and is chosen such as .
In the state space form, the tracking error of the flexible part of the first link can be written as follows: Let .
As already shown in the stability analysis of the rigid part, ,â€‰â€‰andâ€‰â€‰ as . The error dynamics of the flexible part (69) take the same form as the one obtained in [20].
The critical value of is selected when the eigenvalues of change their signs from negative to positive; for example, the internal dynamics become instable [20], where .
5. Simulation Results
The twoflexiblelink manipulator shown in Figure 2 is used to test the control strategy performance. Using the trial and error method, the controllersâ€™ gains are chosen as , , , and .
In Sections 3 and 4, the theoretical development of the control law and stability analysis are given in the general case, that is, for any higher order. In the simulation case, a finite order of Taylor series is fixed. In this section, for simplicity, the Taylor series is limited to the first order. Then, .
Table 1 presents the parameters of the twoflexiblelink manipulator.

Simulation results are given Figures 4â€“14.
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(a)
(b)
(a)
(b)
(a)
(b)
(c)
(d)
(a)
(b)
(a)
(b)
For twoflexiblelink manipulators, the desired workspace trajectory is chosen as lozenge form. The results of the inverse kinematics problem are given in Figures 4â€“7. Workspace trajectory is shown in Figures 4 and 5, virtual space trajectories are given in Figure 6, and joint space trajectories are shown in Figure 7. According to simulation results, the workspace desired trajectory was well transformed to the joint space using the virtual space and the quasistatic approach. For the second subsystem, the Nyquist diagram and the evolution of eigenvalues of are given in Figure 8. For the first subsystem, the eigenvalues evolution of â€‰â€‰is given in Figure 9. Using the system parameters given in Table 1, the critical value of the second and first links are and , respectively. The simulation results were obtained with and . Note that these values satisfy the conditions given in (30) and (38).
A good tracking performance of the new noncollocated output was obtained, as shown in Figure 10. This tracking is confirmed by Figure 11, which shows the tracking error. In the workspace, a good tracking performance is shown in Figure 12, and the tracking error, shown in Figure 13, does not exceed 2â€‰mm in position and position. The good tracking in the workspace is confirmed in Figure 14 that shows the tracking of the desired lozenge. Thus, we can conclude that the control strategy was effective to ensure the tracking of the noncollocated output and to reduce vibrations at the extremity; this explains the satisfactory results obtained using this approach.
6. Conclusion and Future Work
This paper presents a nonlinear distributed control for twoflexiblelink manipulator. For the inverse dynamics, a virtual space, linked with the workspace by a simple kinematics relation as in rigid manipulators, and a quasistatic approach were used. Using this transformation procedure, a workspace desired trajectory (lozenge) has been successfully transformed to the joint space. The distributed control strategy presented in this paper uses the output redefinition technique and consists of stabilizing the flexible manipulators starting with the last joint and flexible link and going backward until the first joint and link. Lyapunov theory was used to prove the asymptotical stability. An adaptive version of this control strategy will be investigated in future work.
Appendices
A. Proof of Proposition 1
Using (3) and (5), the dynamical model can be written as follows: Equation (A.1) can be written as Using the outputs redefinition (31) and (40), we can write where ; ; and . Then, the dynamical model becomes where ; ; ; .
Deducing from (A.5), we get the following internal dynamics: Inserting (A.6) in (A.4), we obtain where ; ; ; and .
The models using the collocated output (3) and the new output (A.4)(A.5) are linked by a nonsingular positive transformation matrix as follows: where and .
Using the transformation , the mass matrix becomes