Dual quaternion-based inverse kinematics of dexterous finger

Jinbao Chen1 , Dong Han2 , Hong Nie3 , Mei Cheng4

1, 2, 3College of Astronautics, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016, China

4Shanghai Institute of Satellite Engineering, Shanghai, 200240, China

1Corresponding author

Journal of Vibroengineering, Vol. 16, Issue 6, 2014, p. 2813-2820.
Received 29 October 2013; received in revised form 1 August 2014; accepted 22 August 2014; published 30 September 2014

Copyright © 2014 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

The inverse kinematics solution of a dexterous robotic finger has a significant impact on the real-time control of the robotic hand. Therefore a rapid method for solving is needed. The classical homogeneous matrix transformation is the most popular method used in robot kinematics. However, for the multi degree-of-freedom (DOF) robotic finger, the matrix parameters cost much storage and the inverse matrix calculation requires a large amount of computational cost. So it is not conducive to the real-time control of the robotic hand. Therefore, a method based on dual quaternions is presented for analysing the kinematics of a multi-DOF (4-DOF) robotic thumb. Firstly, the kinematics equation is expressed by dual quaternions. Then the multivariate kinematic equations are converted to binary quadratic equations with methods of separating variables and variable substitution, which is relatively easy to obtain the closed-form solution of the inverse kinematics. Finally, it proves that the dual quaternions method has advantages over the homogeneous matrix transformation in storage and computational cost by the specific numbers for the robotic thumb, which is conducive to the real-time control of robotic hand.

Keywords: inverse kinematics, dexterous finger, real-time control, homogeneous matrix transformation, dual quaternion.

1. Introduction

Dexterous finger is a multi-DOF and multi-link manipulator and its inverse kinematics have an important influence on the real-time control of dexterous hand. The thumb usually has the most DOFs in the whole hand, so the inverse kinematics of the thumb is the most complicated.

There are several methods for the formulation of the kinematic equations of manipulators with rigid links, including the homogeneous transformation and dual quaternion. The 4×4 matrix or the homogeneous transformation is an almost classical formulation of kinematic equation, which is a point transformation. The dual quaternion based on screw theory is the most efficient and compact way to express a screw displacement and is a line transformation [1, 2]. And Nichoals Aprove that the dual quaternion requires less computational time than the homogeneous transformation for the manipulators with more than 3 DOFs. And the homogeneous matrix needs 12 parameters to express the translation and rotation of rigid body while the dual quaternion 8. Therefore, the dual quaternion is an elegant and useful tool for kinematic analysis in many areas such as navigation [3], computer vision [4], inverse kinematics [5-7] and so on.

The relationship between the terminal and reference coordinates of multi-DOF manipulators is highly nonlinear [8], which makes the inverse kinematics difficult to solve. The solution of inverse kinematics can be classified into two types: closed-form analytical solution and numerical solution. Though the closed-form solution mainly relies on the kinematical structure of the manipulator, it is faster than the numerical solution and can avoid the singular problem and all the possible solutions of inverse kinematics can be given out [9]. The closed-form solution is the most ideal solution. The numerical solution is more general, but it can only find one solution for one set of initial values. In addition, numerical solution requires iterative calculation and a lot of iterative calculation can affect the computational time, which is not conducive to the real-time control. Furthermore, when the numerical method fails to converge, the solution of the inverse kinematics can’t be given out [9-11].

Many researchers choice the dual quaternion to solve the inverse kinematics of the multi-DOF robot manipulators because the kinematic equation with dual quaternion is more concise than that with homogeneous matrix, which is easy to get the closed-form solution. For example, Dgan, et al. adopt the dual quaternion to express the kinematic equation of 7-link 7R mechanism and get the closed-form solution by constructing Dixon resultant [6]; Ni Zhensong, et al. adopt the same method to get the closed-form solution of spatial 6R manipulator [7].

In this paper, the dual quaternion is applied in the inverse kinematics of 4-DOF thumb and the closed-form solution is presented. The algorithm is validated by an example. The whole process is very simple and easy to program, which is conducive to the real-time control of dexterous hand.

2. Quaternions and dual quaternions

2.1. Quaternions

The quaternion is a hypercomplex number with four parts and can well solve the problem of three dimensional space. A general quaternion has the following form [12]:

(1)
q = q 0 + q 1 i + q 2 j + q 3 k = s + v ,

where q0,q1,q2,q3 are real numbers and i,j,k are imaginary unit with the property i2=j2=k2= –1, ij=k, jk=i, ki=j. The quaternion can also be divided into a scalar (s) and a 3×1 vector (v).

Conjugate of the quaternion can be expressed as:

(2)
q * = q 0 - q 1 i - q 2 j - q 3 k .

The multiplication of two quaternions can be expressed as:

(3)
q q ' = q 0 + q 1 i + q 2 j + q 3 k q 0 ' + q 1 ' i + q 2 ' j + q 3 ' k
              = q 0 q 0 ' - q 1 q 1 ' - q 2 q 2 ' - q 3 q 3 ' + ( q 1 q 0 ' + q 0 q 1 ' - q 3 q 2 ' + q 2 q 3 ' ) i
              + q 2 q 0 ' + q 3 q 1 ' + q 0 q 2 ' - q 1 q 3 ' j + q 3 q 0 ' - q 2 q 1 ' + q 1 q 2 ' + q 0 q 3 ' k .

2.2. Dual quaternions

Quaternions can only be used for analysing of the rotation of rigid body around a fixed point. However, a mathematical tool that can express both rotational and translational transformations in a spatial mechanism is needed [13]. Dual numbers were introduced by Clifford W. K. in 1873 and the dual quaternion is a dual number that has two quaternions. The dual quaternion has the following form:

(4)
q ^ = q + ε q ε ,

where q and qε are quaternions and ε is dual unit with the property ε2= 0.

Conjugate of the dual quaternion can be expressed as:

(5)
q ^ * = q * + ε q 3 * .

The multiplication of two dual quaternions can be expressed as:

(6)
q ^ Q ^ = q + ε q ε Q + ε Q ε = q Q + ε q Q ε + q ε Q .

The inverse of a dual quaternion can be expressed as:

(7)
q ^ - 1 = q ^ * q ^ 2 ,       q 0 0 ,

where q^2is the norm of the dual quaternion, which can be expressed as:

(8)
q ^ = q ^ * q ^ = q ^ q ^ * = q q * + ε ( q q ε * + q ε q * ) .

When the dual quaternion q^ is a unit dual quaternion, that is q^=1, q^ should meet:

(9)
q q ε * + q ε q * = 0 , q q * = 1 ,

then q^-1=q^* [14].

The rigid body motion in three-dimensional space can be decomposed into rotation and translation. Let the unit quaternion (R) express the rotation and d express the translation vector. Therefore, the rigid motion can be described in dual queternions like this:

(10)
Q ^ = R + ε d R 2 ,

where R=cos(θ/2)+sinθ/2nr and θ is rotation angle and nr is direction vector of rotation axis (Fig. 1).

Fig. 1. The rotation and translation expressed by dual quaternions

 The rotation and translation expressed by dual quaternions

3. Kinematics of dexterous finger

This paper uses the 4-DOF thumb as an example to illustrate the method that adopts dual queternions to solve the inverse kinematics of dexterous finger. The mechanical structures of the thumb are shown in Fig. 2. And the relationship between the coordinate systems is shown in Fig. 3.

Fig. 2. 4-DOF thumb

 4-DOF thumb

Fig. 3. Relationship between the coordinate systems

Relationship between the coordinate systems

The Denavit-Hartenberg (D-H) parameters of the thumb are listed in Table 1.

Table 1. D-H parameters of the thumb

i
α i - 1   / (°)
a i - 1   / mm
d i / mm
θ i / (°)
1
90.0
10
0
θ 1
2
–90.0
17
0
θ 2
3
40.0
37
7
θ 3
4
0
38
0
θ 4

In this paper, dn is defined as the distance from Xn-1 to Xn along the Zn-1 axis; the other parameters are defined according to the classical definitions.

According to Fig. 3 and Table 1, the transformations of each link lever are expressed by dual queternions as follow:

(11)
Q ^ 0 1 = P ^ 0 X a 0 R ^ 0 Y θ 1 = 1 + ε a 0 i 2 c o s θ 1 2 + s i n θ 1 2 j + ε * 0 ,
(12)
Q ^ 1 2 = P ^ 1 X a 1 R ^ 1 Z θ 2 = 1 + ε a 1 i 2 c o s θ 2 2 + s i n θ 2 2 k + ε * 0 ,
(13)
Q ^ 2 3 = P ^ 2 Z d 3 R ^ 2 X α 2 P ^ 2 X a 2 R ^ 2 Z θ 3
            = 1 + ε d 3 k 2 c o s α 2 2 + s i n α 2 2 i + ε * 0 1 + ε a 2 i 2
              c o s θ 3 2 + s i n θ 3 2 k + ε * 0 ,
(14)
Q ^ 3 4 = P ^ 3 X a 3 R ^ 3 Z θ 4 = 1 + ε a 3 i 2 c o s θ 4 2 + s i n θ 4 2 k + ε * 0 ,
(15)
Q ^ 4 5 = P ^ 4 X a 4 = 1 + ε a 4 i 2 ,

where a0, a1, a2, a3, a4, α2, d3 are known D-H parameters; θn represents the angle displacement of the n-joint, n=1,2,3,4. The position and orientation of the end of the robotic thumb in the Cartesian coordinates are expressed by a dual queternion as follow:

(16)
N ^ = r 0 + r i i + r j j + r k k + ε s 0 + s i i + s j j + s k k .

where r0, ri, rj, rk, s0, si, sj and sk represent the position-orientation parameters of the end.

So, the kinematics equation of the thumb is given by:

(17)
Q ^ 0 1 Q ^ 1 2 Q ^ 2 3 Q ^ 3 4 Q ^ 4 5 = N ^ .

Eq. (17) is separated variables. A new equation is given by:

(18)
Q ^ 0 1 Q ^ 1 2 = N ^ Q ^ 2 3 - 1 Q ^ 3 4 - 1 Q ^ 4 5 - 1 .

It is easy to prove that Q^23, Q^34 and Q^45 meet Eq. (9), so their inverse is equal to their conjugate. Therefore, Eq. (18) can also be expressed as:

(19)
Q ^ 0 1 Q ^ 1 2 = N ^ Q ^ 2 3 * Q ^ 3 4 * Q ^ 4 5 * .

Setting left and right sides of Eq. (19) as U^ and V^:

(20)
U ^ = V ^ ,

where:

(21)
U ^ = U 1 + ε U 2 = u 11 + u 12 i + u 13 j + u 14 k + ε u 21 + u 22 i + u 23 j + u 24 k ,
(22)
V ^ = V 1 + ε V 2 = v 11 + v 12 i + v 13 j + v 14 k + ε v 21 + v 22 i + v 23 j + v 24 k .

If two dual queternions are equal, each element should be equal. So, u1i=v1i, u2i=v2i (i= 1, 2, 3, 4) and we get:

(23)
D j c θ 1 c θ 2 s θ 1 c θ 2 c θ 1 s θ 2 s θ 1 s θ 2 = C j ,       j = 1 ,   2 ,

where cθk=cos(θk/2), sθk=sin(θk/2), (k= 1, 2, 3, 4) and Dj is a 4×4 matrix that is decided by the D-H parameters of the thumb and the position and orientation of end-effector; Cj is a 4×1 vector, whose elements are the linear expression of cθ3cθ4, sθ3cθ4, cθ3sθ4 and sθ3sθ4.

Form Eq. (23), we can get:

(24)
D 1 - 1 C 1 = D 2 - 1 C 2 .

The left and right sides of Eq. (24) are divided by cθ3cθ4 and we get equations as follow:

(25)
E q e t 3 , t 4 = a e + m e t 3 + n e t 4 + p e t 3 t 4 = 0 ,       e = 1 ,   2 ,   3 ,   4 ,

where t3=sθ3/cθ3, t4=sθ4/cθ4, and ae, me, ne, pe are known quantity and decided by the D-H parameters of the thumb and the position and orientation of end-effector. Eq. (25) are binary quadratic equations. Hence, it is easy to get the closed-form solution of Eq. (25), that is to say, we can get the closed-form solution of θ3 and θ4. The closed-form solution of θ1 and θ2 can be similarly solved through Eq. (23) with known θ3 and θ4.

4. Verification and comparison

4.1. Verification with an example

Firstly, a simple example is given to illustrate the validity of the proposed method. Let θ1= 0.353 rad, θ2= 0.434 rad, θ3= 0.625 rad, θ4= 0.764 rad, which are the desired values. Then the position and orientation of the end-effector is calculated by the desired angle displacements of each joint, which is expressed by the dual quaternion as follow:

(26)
N ^ = 0.5943 + 0.4294 i - 0.0536 j + 0.6779 k
        + ε - 0.0128 + 0.0487 i - 0.0064 j - 0.0123 k .

Substitute the D-H parameters and Eq. (26) in Eq. (25) and get:

(27)
0.5671   - 0.6969 t 3 +   1.1478 t 4 + 2.5524 t 3 t 4 = 0 ,
- 9.1484   + 7.8760 t 3 + 14.4637 t 4 + 6.0844 t 3 t 4 = 0 ,    
- 1.8693   + 1.9064 t 3 +   3.0650 t 4 + 0.1278 t 3 t 4   = 0 ,
- 2.3280   + 2.1035 t 3 +   4.0748 t 4 + 0.1538 t 3 t 4   = 0 .

The solution of Eq. (27) is t3= 0.3262, t4= 0.4002, that is to say, θ3'=0.6304, θ4'=0.7610. Substitute θ3' and θ4' in Eq. (23) and get θ1'= 0.3513,θ2'= 0.4315. The results of inverse kinematics testify that the algorithm is effective.

Secondly, another example based on trajectories is presented to illustrate the practicality of the method. The trajectories are planned in joint space. For the stability of the manipulator’s motion, linear interpolation taking parabolic transition is adopted to plan the trajectories of each joint. And the planned trajectories of the four joints are respectively shown in Fig. 4(a), Fig. 5(a), Fig. 6(a) and Fig. 7(a). According to the trajectories in joint space, the trajectories of the end of the robtic thumb in Cartesian space can be computed. So, the inverse kinematics solution can be gotten by using the proposed method with knowing the trajectories of the end, as shown in Fig. 4(b), Fig. 5(b), Fig. 6(b) and Fig. 7(b). The computed results of the inverse kinematics are almost equal to the planned trajectories in joint space by comparing the Fig. 4-7. Therefore, the dual quaternion can be used for solving the inverse kinematics of the multi-DOF manipulator in practice.

4.2. Comparison with homogeneous transformation

The above two examples show the feasibility of the dual quaternion for solving the inverse kinematics. Furthermore, the dual quaternion outperforms homogeneous transformation in three aspects for the multi-DOF (more than 3 DOFs) manipulator [1], as follow:

1) Storage cost: For one transformation, the dual quaternion requires eight memory locations, while homogeneous transformation needs12, since the elements of the fourth row are known constants. According to Eq. (17), we know that the dual quaternion method requires 48 memory locations at least for an inverse kinematic computation, while homogeneous transformation requires 72. It is obvious that the storage cost of the dual quaternion method is minimum.

2) Computational cost: We define that the sign+ represents addition and subtraction operations and the sign* represents multiplication operations. The multiplication of two 4×4 homogeneous matrices needs 48* and 36+ [1]. The inverse of homogeneous matrix requires 9* and 6+. The multiplication of two dual quaternions needs 48* and 36+. The inverse of dual quaternions needs no computational cost. Therefore, the computation of (18) needs 219* and 162+ for homogeneous transformation, while 192* and 144+ for the dual quaternion.

3) Computational time: The storage and computational cost can both affect the computational time. Fetching an operand from memory will cost computational time. And the computational time of an addition is less than one multiplication. It turns out that the computational time of homogeneous transformation is more than the dual quaternion for the inverse kinematics of 4-DOF thumb according to the conclusions of 1) and 2).

In one word, the homogeneous transformation method needs more storage and computational cost than the dual quaternion method, which leads to more computational time that the former requires than the latter.

5. Conclusions

The dual quaternion is applied in the inverse kinematics solution of 4-DOF dexterous thumb, which greatly simplifies the process of getting the closed-form solution and provides a new thinking for the inverse kinematics solution of multi-DOF dexterous finger.

Through comparing with the homogeneous transformation, we can conclude that the dual quaternion can simplify the kinematics equation, effectively reduce the storage and computational cost, greatly improve computational speed and satisfy the requirements of the real-time control of dexterous hand. However, the dual quaternion has no advantage over the homogeneous transformation for the finger with less than 4 DOF.

Acknowledgements

The authors would like to thank the anonymous reviewers for their critical and constructive review of the manuscript. This work was supported by National Natural Science Foundation of China (Grant No. 51105196), Natural Science Foundation of Jiangsu Province (Grant No. BK2011733), Science and Technology Innovation Fund of Shanghai Aerospace (Grant No. SAST201320), and Shanghai Key Laboratory of Deep Space Exploration Technology.

References

  1. Nichoals A. Aspragathos, John K. Dimitros A comparative study of three methods for robot kinematics. IEEE Transactions on Systems, Vol. 28, Issue 2, 1998, p. 135-145. [Search CrossRef]
  2. Sariyildiz E., Temeltas H. Solution of inverse kinematic problem for serial robot using dual quaterninons and plücker coordinates. International Conference on Advanced Intelligent Mechatronics, 2009, p. 338-343. [Search CrossRef]
  3. Li Jing, Wang Huinan, Liu Haiying Application of dual quaternion in spacecraft relative navigation. Journal of Applied Sciences, Vol. 30, Issue 3, 2012, p. 311-316, (in Chinese). [Search CrossRef]
  4. Konstantinos Daniilidis Hand-eye calibration using dual quaternions. The International Journal of Robotics Research, Vol. 18, 1999, p. 286-298. [Search CrossRef]
  5. Perez Alba, McCarthy J. M. Dual quaternion synthesis of constrained robotic systems. Journal of Mechanical Design, Vol. 126, 2004, p. 425-435. [Search CrossRef]
  6. Gan D., Liao Q., Wei S., et al. Dual quaternion-based inverse kinematics of the general spatial 7R mechanism. Journal of Mechanical Engineering Science, Vol. 222(c), 2008, p. 1593-1598. [Search CrossRef]
  7. Ni Zhensong, Liao Qizheng, Wei Shimin, et al. Dual four element method for inverse kinematics analysis of spatial 6R manipulator. Journal of Mechanical Engineering, Vol. 45, Issue 11, 2009, p. 25-29, (in Chinese). [Search CrossRef]
  8. Rodolphe J. Gentili, Hyuk Oh, Javier Molina, et al. Cortex inspired model for inverse kinematics computation for a humanoid robotic finger. 34th Annual International Conference of the IEEE EMBS, 2012, p. 3052-3055. [Search CrossRef]
  9. Yang Si, Qingxuan Jia, Gang Chen, et al. A complete solution to the inverse kinematics problem for 4-DOF manipulator robot. 8th Conference on Industrial Electronics and Applications, 2013, p. 1880-1884. [Search CrossRef]
  10. Andrew A. Goldenberg, Member B. Benhabib, et al. A complete generalized solution to the inverse kinematics of robots. IEEE Journal of Robotics and Automation, Vol. 1, Issue 1, 1985, p. 14-20. [Search CrossRef]
  11. Jing Huang, Xianlun Wang, Dongsheng Liu, et al. A new method for solving inverse kinematics of an industrial robot. International Conference on Computer Science and Electronics, 2012, p. 53-56. [Search CrossRef]
  12. Hoang-Lan Pham, Véronique Perdereau, Bruno Vilhena Adorno, et al. Position and orientation control of robot manipulators using dual quaternion feedback. The International Conference on Intelligent Robots and Systemsg, 2010, p. 658-663. [Search CrossRef]
  13. Liu Yanzhu The mathematical expression of rigid body attitude. Mechanics in Engineering, Vol. 33, Issue 1, 2008, p. 98-101, (in Chinese). [Search CrossRef]
  14. Nan Ji, Wang Jianguo, Chen Jianyu Rigid dual quaternions related problems. Fire Control and Command Control, Vol. 37, 2012, p. 15-17, (in Chinese). [Search CrossRef]