*Article* **Investigation of Cyclicity of Kinematic Resolution Methods for Serial and Parallel Planar Manipulators**

**Maurizio Ruggiu 1,\*,† and Andreas Müller 2,†**


**Abstract:** Kinematic redundancy of manipulators is a well-understood topic, and various methods were developed for the redundancy resolution in order to solve the inverse kinematics problem, at least for serial manipulators. An important question, with high practical relevance, is whether the inverse kinematics solution is cyclic, i.e., whether the redundancy solution leads to a closed path in joint space as a solution of a closed path in task space. This paper investigates the cyclicity property of two widely used redundancy resolution methods, namely the projected gradient method (PGM) and the augmented Jacobian method (AJM), by means of examples. Both methods determine solutions that minimize an objective function, and from an application point of view, the sensitivity of the methods on the initial configuration is crucial. Numerical results are reported for redundant serial robotic arms and for redundant parallel kinematic manipulators. While the AJM is known to be cyclic, it turns out that also the PGM exhibits cyclicity. However, only the PGM converges to the local optimum of the objective function when starting from an initial configuration of the cyclic trajectory.

**Keywords:** kinematic redundancy; cyclicity; augmented Jacobian method; projected gradient method

### **1. Introduction**

Kinematically redundant serial kinematic manipulators (SKM) are gaining importance for dedicated solution in flexible automation due their advantageous properties. Only a few kinematically redundant parallel kinematic manipulators (PKM) have been proposed, and they have not yet seen a similar application in industry. Kinematic redundancy offers various advantages, and it can in particular be exploited to improve dexterity, to avoid singularities, to circumvent obstacles, and eventually to increase the workspace. These advantages are accompanied with an increased complexity of the control and with increased costs, however. The complexity is due to the non-uniqueness of the solution to the geometric inverse kinematics problem (IKP), which consists of finding a vector **<sup>q</sup>** <sup>∈</sup> <sup>V</sup>*<sup>n</sup>* of generalized coordinates (joint coordinates in SKM) for a given vector **<sup>p</sup>** <sup>∈</sup> <sup>R</sup>*<sup>m</sup>* of end-effector (EE) coordinates satisfying the forward kinematics relation

$$\mathbf{p} = f(\mathbf{q}) \tag{1}$$

where *f* is the forward kinematics mapping. The corresponding velocity IKP, in a given pose **q**, is to find a vector of generalized coordinates **q**˙ for given EE-velocity **p**˙ satisfying the velocity forward kinematics

$$
\dot{\mathbf{p}} = \mathbf{J}(\mathbf{q})\dot{\mathbf{q}}\tag{2}
$$

where **<sup>J</sup>**(**q**) <sup>∈</sup> <sup>R</sup>*m*×*<sup>n</sup>* is the forward kinematics Jacobian.

A manipulator is kinematically redundant if *r* = *n* − *m* > 0, where *r* is called the degree of redundancy, *n* is the degree of freedom (DOF), and *m* denotes the dimension of

**Citation:** Ruggiu, M.; Müller, A. Investigation of Cyclicity of Kinematic Resolution Methods for Serial and Parallel Planar Manipulators. *Robotics* **2021**, *10*, 9. https://doi.org/10.3390/ robotics10010009

Received: 20 October 2020 Accepted: 30 December 2020 Published: 3 January 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

the task space. For a redundant manipulator, the Jacobian is not square, and the inverse kinematics has *r*-dimensional solution set.

A number of solution techniques for solving the geometric IKP problem for redundant manipulators have been introduced in the last three decades, almost exclusively focusing on SKM [1]. Most of them determine a solution of the geometric inverse kinematics problem from a solution of the velocity inverse kinematics problem with the aim to provide an online solution scheme that can be executed on the robot controller with a well-defined computation time. As a consequence, all these method yield a local solution rather than the global solution of the geometric inverse kinematics problem. The two most widely used local solution schemes are the following:


Most industrial use cases involve repetitive tasks. That is, the EE performs a cyclic motion. An obvious requirement for safe application of redundant robots is that the IK solution is well-determined in the sense that for given pose **p** of the EE the IK solution scheme always yields the same pose **q** of the robot. This property is called cyclicity, which means that for a closed EE-path the IK solution method gives a closed path in joint space. It is well-known that null-space methods cannot be guaranteed to be cyclic [2]. This is an important drawback of the null-space methods since the motion of the chain is unpredictable during the motion. The cyclicity problem was first analyzed in a differential geometric framework in [3] and then explored in [4,5] for the pseudoinverse solution. In [6] the authors were able to design a feedback control law to produce convergence of the joint motion towards a cyclic trajectory assuming the null space vector as a linear map of the EE trajectory.

While task augmentation methods are cyclic [6], they suffer from algorithmic singularities [7]. It must be mentioned that this topic was extensively used and studied for SKM (aiming to avoid collision and singularities). To the authors knowledge these same techniques have been applied very rarely to PKMs as in [8]. For the sake of the completeness, we must point out that recently these classical schemes have been applied to the hyper-redundant manipulators [9,10] as well as algorithms based on the machine learning have been applied to the redundancy resolution [11].

In this paper the projected gradient method (PGM) and the augmented Jacobian method (AJM) are compared by means of examples when applied to SKM and PKM. The methods are designed to minimize an objective function. In the following, and w.l.o.g., the objective function will be the sum of the squared joint variables. The comparison regards (1) the cyclicity of the joint trajectories, and (2) the dependency of the (cyclic) joint trajectory with respect to changes in the initial configuration. The sensitivity is an aspect that has not yet been addressed in the literature as well as the cyclicity of the solutions for PKMs. It will be shown that only the projected gradient method is capable of minimizing the objective function independently of the initial configurations. Results are reported for both SKM and PKM.

This is the first publication that addresses the cyclicity of PGM. The novel contribution is a numerical investigation and comparison of the cyclicity of PGM and AJM when applied to SKM as well as to PKM.

The paper is organized as follows. In Section 2 the two methods for redundancy resolution at the velocity level are presented. Numerical results are shown in Section 3. The paper concludes with a discussion in Section 4.

### **2. Redundancy Resolution Methods**

Redundancy resolution of the differential kinematics of a manipulator is expressed as:

$$
\dot{\mathbf{q}} = \mathbf{J}'' \dot{\mathbf{p}}.\tag{3}
$$

**<sup>J</sup>**# is the Moore–Penrose pseudoinverse of the manipulator Jacobian **<sup>J</sup>** <sup>∈</sup> <sup>R</sup>*m*×*n*. In order to avoid integration drift, an algorithm solution (CLIK: closed-loop inverse kinematics) of Equation (3) is widely adopted and the geometric inverse kinematic problem expressed as:

$$\dot{\mathbf{q}} = \mathbf{J}^{\#}(\dot{\mathbf{p}}\_d + \mathbf{G}\mathbf{e}),\tag{4}$$

where **e** = **p***<sup>d</sup>* − **p** represents the error between the desired and actual EE location and **G** is a positive–definite diagonal gain matrix. Numerical integration of the Equation (4) given the initial configuration **q**<sup>0</sup> followed by the solution of the forward position kinematic equations closes the feedback loop.

1. PGM

The PGM exploits the fact that a general solution of the differential kinematics can be substituted to Equation (4) when a desired joint rate vector **v** is projected into the null space of **J**:

$$\dot{\mathbf{q}} = \mathbf{J}^{\boldsymbol{\theta}}(\dot{\mathbf{p}}\_d + \mathbf{G}\mathbf{e}) - \mathbf{P}\mathbf{v},\tag{5}$$

with **<sup>P</sup>** = (**<sup>I</sup>** <sup>−</sup> **<sup>J</sup>**#**J**) defined as the projector matrix. The added term generates self– motion of the kinematic chain without affecting the EE velocity. Vector **v** can be chosen in order to make a scalar objective function *h*(**q**) stationary by using the gradient projection method, such that **v** = ( *<sup>∂</sup>h*(**q**) *<sup>∂</sup>***<sup>q</sup>** )*T*, with **<sup>v</sup>** <sup>∈</sup> <sup>R</sup>*n*. *<sup>h</sup>*(**q**) may be any analytical differentiable function expressed in terms of the joint variables **q** only.

2. AJM

A different approach is followed in the AJM. An additional constraint task is imposed to the original task of the EE. Following [7,12], the objective function *h*(**q**) is projected onto the null space of **J** and imposed to be zero. Formally we can write:

$$\mathbf{g}(\mathbf{q}) = \mathbf{Z}^{\mathrm{T}} \mathbf{v} = \mathbf{0}.\tag{6}$$

**<sup>Z</sup>** <sup>∈</sup> <sup>R</sup>*n*×*<sup>r</sup>* is an orthonormal basis for the null space of **<sup>J</sup>** and **<sup>v</sup>** is the gradient of *<sup>h</sup>*(**q**) with respect to the joint variables as in the previous method. Therefore, Equation (6) yields *r* independent constraints keeping *h*(**q**) at the extreme at each time of the trajectory starting from the initial configuration **q**0.

The added Jacobian **<sup>J</sup>***<sup>a</sup>* <sup>∈</sup> <sup>R</sup>*r*×*<sup>n</sup>* can simply obtained as

$$\mathbf{J}\_a = \frac{\partial \mathbf{g}(\mathbf{q})}{\partial \mathbf{q}},\tag{7}$$

that leads to the CLIK kinematic resolution expressed as:

$$\dot{\mathbf{q}} = \mathbf{J}\_{\text{aug}}^{-1}(\begin{bmatrix} \dot{\mathbf{p}}\_d \\ \mathbf{0} \end{bmatrix} + \begin{bmatrix} \mathbf{G} & \mathbf{0} \\ \mathbf{0} & \mathbf{0} \end{bmatrix} \begin{bmatrix} \mathbf{e} \\ \mathbf{0} \end{bmatrix}), \quad \mathbf{J}\_{\text{aug}} = \begin{bmatrix} \mathbf{J} \\ \mathbf{J}\_a \end{bmatrix}. \tag{8}$$

The two methods are, therefore, designed to reach the same goal in order to have a direct comparison between them.

As pointed out in [6] for SKM, the cyclicity of the joints trajectories for a closed EE path can be reached exactly, asymptotically or it cannot even be reached. Exact cyclicity is obtained whenever **v** = **L**(**q**)**p**˙ at the steady–state, i.e., **e** = **0**. It has been shown in [6] that the augmented Jacobian can be expressed in this form thus proving its cyclicity. Furthermore, the cyclicity property may be verified on–line by calculating that the Lie

brackets formed by any columns of the control matrix, **J**−<sup>1</sup> *aug*, are linear combinations of the columns (involutive property) [13]. This calculation was indeed performed in this paper.

On the contrary, PGM cannot guarantee to lead to exact or even asymptotic cyclicity *a priori*. Indeed, it has been proved that control scheme in Equation (4) can only be asymptotic cyclic with **v** = **K**(**q***cyclic* − **q**), with **q***cyclic* representing the joint trajectory to converge to in order to reach repeatability.

### **3. Numerical Simulations**

PGM and AJM are implemented for the serial manipulators 4R and the parallel 2RRP planar manipulators and for a 6R industrial UR10e robot. The planar manipulators have a degree of redundancy *r* = 2, while the 6DOF arm has *r* = 3. The end effectors of the manipulators can track a circle at a given constant speed.

For all the manipulators, the objective function to be minimized was chosen to be:

$$h(\mathbf{q}) = \sum\_{i=1}^{n} q\_{i\prime}^{2} \tag{9}$$

with *qi* <sup>≡</sup> *<sup>θ</sup><sup>i</sup>* for the SKMs and **<sup>q</sup>** <sup>=</sup> *ψ*1, *ρ*1, *ψ*2, *ρ*<sup>2</sup> for the PKM. The simple nature of *h*(**q**) strongly simplifies the computations. In the PGM method, indeed, only the analytical calculation of the gradient of the objective function is required. Instead, AJM requires the analytical calculation of the gradient of **g**(**q**) in Equation (6) to obtain the added rows of the Jacobian. First, an analytical calculation of **Z** is needed. This can be avoided using an alternative procedure proposed in [14] to calculate only **Z** numerically.

The geometrical data and the data of the EE trajectory are given in Tables 1–3.

**Table 1.** 4R. *a*: links length; (*x*0, *y*0), *R*: coordinates of the centre and radius of the EE trajectory.


**Table 2.** 2RRP. *a*: links length; *H*: base width, (*x*0, *y*0), *R*: coordinates of the centre and radius of the EE trajectory.


**Table 3.** 6R. (*x*0, *y*0, *z*0), *R*: coordinates of the centre and radius of the EE trajectory. EE trajectory on the plane *y* = −0.256 m. *d*1, *a*2, *a*3, *d*4, *d*5, *d*6: *Denavit-Hartenberg* geometrical parameters.


### *3.1. Serial 4R*

The Jacobian of the manipulator in Figure 1a is trivial and it is not shown in the paper. Instead, Equation (10) shows **g**(**q**) used in the AJM:

$$\mathbf{g(q)} = \frac{2}{(s\_{234} + s\_{23} + s\_2)} \left( \begin{bmatrix} (\theta\_3 - \theta\_2)(s\_{234} + s\_{23}) + (\theta\_1 - \theta\_2)(s\_{34} + s\_3) + \theta\_3 s\_2\\ (\theta\_4 - \theta\_2)s\_{234} + (\theta\_1 - \theta\_2)(s\_{34} + s\_4) + \theta\_4(s\_{23} + s\_2) \end{bmatrix} \right) \tag{10}$$

**<sup>g</sup>**(**q**) is then differentiated with respect to the joint angles to obtain **<sup>J</sup>***<sup>a</sup>* <sup>∈</sup> <sup>R</sup>2×4.

**Figure 1.** The test cases: (**a**) Planar 4R, (**b**) Planar 2RRP, (**c**) UR10 industrial 6R robot.

In Figure 2 we show *h*(**q**) from both the methods for the planar 4R. The initial configuration **q**<sup>0</sup> was calculated such to be a solution of the inverse position problem and that **g**(**q**0) = **0**.

From Figure 2 we see that *h*(**q**) from the methods are in very good agreement. Both can keep the objective function stationary. Accordingly, also the joint trajectories of the two methods, shown in Figure 3, only exhibit small differences.

From Figure 3 we also see that both methods produce cyclic joints trajectories. In the case of the PGM we can verify that only by inspection while for the AJM we performed the numerical calculations of the Lie brackets associated to the augmented Jacobian. To have cyclicity we need to have that for any two columns **c***<sup>i</sup>* and **c***<sup>j</sup>* of **J**−<sup>1</sup> *aug*, their Lie bracket <sup>L</sup>*<sup>k</sup>* is a linear combination of the columns of the *control matrix* **<sup>J</sup>**−<sup>1</sup> *aug*. In this case, for each step of calculation we obtain <sup>L</sup>*<sup>k</sup>* of the order of 10−<sup>15</sup> proving the cyclicity. Figure <sup>4</sup> shows the objective function whenever **g**(**q**∗ <sup>0</sup> ) = **0** as **q**<sup>∗</sup> <sup>0</sup> = **q**<sup>0</sup> + *<sup>v</sup>* with *<sup>v</sup>* = 1, −1, 1, −1 . As expected, the results show that AJM can only keep the constraint as it is, while PGM quickly reaches the minimum.

The inability of AJM to converge to the minimum is stressed in Figure 5 where *h*(**q**) is plotted with different initial conditions.

**Figure 2.** *h*(**q**) in 4R. PGM: − (red line), AJM: −− (blue line).

**Figure 3.** *θ<sup>i</sup>* in 4R. PGM: − (red line), AJM: −− (blue line).

**Figure 4.** *h*(**q**) in 4R with = 0.3. PGM: − (red line), AJM: −− (blue line).

**Figure 5.** AJM: *h*(**q**) in 4R. = 0: −− (black line), = 0.25: −.− (red line), = 0.5: − (blue line).

We then tested the robustness of PGM as the initial conditions varied. In Figure 6, *h*(**q**) is computed by PGM when starting from different initial configurations. The method

allows reaching the minimum even when = 1.5. In this case the method converges in *Nc* = 1.2570 cycles chosen Δ*h* = (*h* − *h*0) = 1/1000 as threshold. Similarly, the joints trajectories coincide quickly. This is a striking characteristic of this method that proves to be robust with respect to the initial configuration of the motion that can be inaccurate for some reasons. In this case the method can reach the minimum curve even when the gap from the correct initial configuration is about 90◦.

**Figure 6.** PGM: *h*(**q**) in 4R. = 0: −− (black line), = 0.5: −.− (red line), = 1.0: − (blue line), = 1.5: ··· (green line).

### *3.2. Parallel 2RRP*

The kinematic relation of the PKM is [15]:

$$\mathbf{J}\_p \dot{\mathbf{p}} = \mathbf{J}\_q \dot{\mathbf{q}}.\tag{11}$$

The forward kinematics Jacobian in (2) is thus **J** = **J**−<sup>1</sup> *<sup>p</sup>* **J***q*. For the planar PKM in Figure 1b), the Jacobians **J***p* and **J***q* are

$$\begin{array}{rcl} \mathbf{J}\_{p} &=& \begin{bmatrix} \mathbf{x} - ac\_{\psi\_{1}} & y - as\_{\psi\_{1}} \\ (\mathbf{x} - H) - ac\_{\psi\_{2}} & y - ac\_{\psi\_{2}} \end{bmatrix} \\ \mathbf{J}\_{q} &=& \begin{bmatrix} a(yc\_{\psi\_{1}} - xs\_{\psi\_{1}}) & \rho\_{1} & 0 & 0 \\ 0 & 0 & a(yc\_{\psi\_{2}} - (\mathbf{x} - H)s\_{\psi\_{2}}) & \rho\_{2} \end{bmatrix} . \end{array} \tag{12}$$

The term **g**(**q**) used in the AJM is

$$\mathbf{g(q)} = 2(\begin{bmatrix} \rho\_1 (1 - \frac{\psi\_1}{a(\mathcal{yc}\_{\varphi\_1} - \mathbf{x}s\_{\varphi\_1})}) \\ \rho\_2 (1 - \frac{\psi\_2}{a(\mathcal{yc}\_{\varphi\_2} - (\mathbf{x} - H)s\_{\Psi\_2})}) \end{bmatrix}). \tag{1.3}$$

As with the planar serial case, Figure 7 shows *h*(**q**) from both the methods for the PKM. The initial configuration **q**<sup>0</sup> was calculated such to be a solution of the inverse position problem and that **g**(**q**0) = **0**.

**Figure 7.** *h*(**q**) in 2RRP. PGM: − (red line), AJM: −− (blue line).

Also in this case the methods coincide.

We repeat the tests considering **g**(**q**∗ <sup>0</sup> ) = **0**. However, in this case because of the different types of joints we use **q**∗ <sup>0</sup> = **q**<sup>0</sup> + *<sup>v</sup>* with *<sup>v</sup>* = 1, 1 *m*, 1, 1 *m* . The methods work consistently for the PKM as well. Only PGM can quickly reach the minimum of the objective function as shown in Figure 8.

Therefore, AJM cannot reach the minimum as pointed out in Figure 9.

The robustness of PGM when varying the initial conditions for the PKM is essentially the same of that for the serial 4*R*. The results of the test are shown in Figure 10. The convergence rate in terms of number of cycles was computed of *Nc* = 0.9975.

**Figure 8.** *h*(**q**) in 2RRP with = 0.15. PGM: − (red line), AJM: −− (blue line).

**Figure 9.** AJM: *h*(**q**) in 2RRP. = 0: −− (black line), = 0.1: −.− (red line), = 0.15: − (blue line).

**Figure 10.** PGM: *h*(**q**) in 2RRP. = 0: −− (black line), = 0.5: −.− (red line), = 1.0: − (blue line), = 1.5: ··· (green line).

### *3.3. Serial 6R*

Finally the AJM and PGM methods were tested for a 6DOF industrial robotic arm UR10e. In this case, due to their sizes, it is not possible to show explicitly the terms involved into calculations. Figure 11 shows that also in this case the methods coincide when **g**(**q**0) = **0**, i.e., the objective function is at the minimum at the initial configuration.

Also, in the case of **q**∗ <sup>0</sup> = **q**<sup>0</sup> + *<sup>v</sup>* with *<sup>v</sup>* = 1, −1, 1, −1, 1, −1 , the methods work as in the other cases as shown in Figure 12.

The response of the AJM in terms of variation of the initial conditions is shown in Figure 13.

Sensitivity of the PGM method to the initial configuration is shown in Figure 14. The behavior is similar to the other cases.

**Figure 11.** *h*(**q**) in 6R. PGM: − (red line), AJM: −− (blue line).

**Figure 12.** *h*(**q**) in 6R with = 0.5. PGM: − (red line), AJM: −− (blue line).

**Figure 13.** AJM: *h*(**q**) in 6R. = 0: −− (black line), = 0.25: −.− (red line), = 0.5: − (blue line).

**Figure 14.** PGM: *h*(**q**) in 6R. = 0: −− (black line), = 0.5: −.− (red line), = 1.0: − (blue line), = 1.5: ··· (green line).

The convergence measure in terms of number of cycles is *Nc* = 2.1145. In order to check the robustness of the methods, two further test were performed. First, the circular trajectory was tracked with a variable EE velocity obtained as *v* = ˙ *θR* with

$$
\theta = \frac{2\pi}{T}(1 - \cos 2\pi \tau), \quad \tau = \frac{t}{T}. \tag{14}
$$

Then, a test with a trajectory represented by *Bernoulli's Lemniscate* (∞-like curve) tracked with an EE velocity *v* = ˙ *θr*. In this case, the angular velocity follows Equation (14) while the radial velocity is kept constant. The authors did not find discrepancies with respect to the case of the circle trajectory tracked at the constant velocity for all the manipulators. As an example, Figure 15 shows the robustness of the PGM for the PKM when the EE is tracking the ∞-like curve at the *v* velocity.

**Figure 15.** PGM with *v* and ∞-trajectory: *h*(**q**) in 2RRP. = 0: − (black line), = 0.75: −− (red line), = 1.5: −. (blue line).

### **4. Conclusions**

This paper is the combination of two research directions pursued independently: investigation of the trajectories cyclicity and of the gradient-based algortihms for redundancy resolution. In the past, the augmented Jacobian method was shown to lead to cyclic trajectories while the projected gradient method for the inverse kinematics was only used to maximize an objective function. It was never attempted to use the projected gradient method to generate cyclic inverse kinematics solution. In this paper, the cyclicity property of the gradient-based and of the augmented Jacobian method for the redundancy resolution is investigated by means of numerical examples. This paper is the first presented in the literature in which a such comparison is performed.

Briefly, the main results of this paper are: (*a*) The PGM leads to cyclic trajectories at which the objective function attains a minimum, and converges to such cyclic trajectory even if the initial configuration is far from a point where the objective is minimized. The tests on the robustness with respect to the initial configurations showed similar behavior for all the manipulators investigated. On the contrary, as expected, the AJM can only keep the objective function at the value at the initial configuration. (*b*) The results are identical for serial and parallel manipulators since both kinematically characterized by a forward/inverse kinematic mapping. Both methods are cyclic either for the SKMs or for the PKM investigated. The cyclicity property of the AJM was confirmed by numerical calculation of the relevant Lie brackets, while the asymptotic cyclicity of the PGM was confirmed by numerical calculation. These results are novel and it remains to further investigate the properties of the PGM. Moreover, the analyses carried out in this paper are numeric and further research needs to address this topic from an analytical point of view.

**Author Contributions:** Conceptualization, M.R. and A.M.; software, M.R. and A.M.; investigation, M.R. and A.M.; writing—original draft preparation, M.R. and A.M.; writing—review and editing, M.R. and A.M.; Both authors have equally contributed to this paper. All authors have read and agreed to the published version of the manuscript.

**Funding:** The second author acknowledges support from the LCM K2 Center for Symbiotic Mechatronics within the framework of the Austrian COMET-K2 program.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author.

**Conflicts of Interest:** The authors declare no conflict of interest.

### **References**


*Article*
