NUMERICAL METHOD FOR OBTAINING THE INVERSE KINEMATICS OF SIX-DEGREE-OF-FREEDOM SERIAL ROBOT WITH AN OFFSET WRIST

Abstract

Disclosed is a method for obtaining the inverse kinematics of a six-degree-of-freedom serial robot with an offset wrist. The method uses the analytical solution of the inverse kinematics of the six-degree-of-freedom serial robot with a non-offset wrist as an approximate solution of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist and an initial point for iteration, and obtains a numerical solution of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist meeting the accuracy through continuously iterative approaching. The present disclosure has faster convergence and less calculation amount relative to the traditional calculation method, reduces the computation burden for a robot controller, and improves the real-time performance.

Claims

1. A numerical method for obtaining the inverse kinematics of a six-degree-of-freedom serial robot with an offset wrist, comprising the following steps: A. use the pose of the terminal coordinate frame of the six-degree-of-freedom serial robot with an offset wrist as the pose of the terminal coordinate frame of a six-degree-of-freedom serial robot with a non-offset wrist to obtain an analytical solution .sub.1.sub.6 of the inverse kinematics of the six-degree-of-freedom serial robot with a non-offset wrist; B. use the analytical solution .sub.1.sub.6 as an approximate solution of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist and an initial point for iterative calculation using an iteration method; use .sub.1.sub.6 to obtain the approximate pose of the terminal coordinate frame according to a formula for obtaining the forward kinematics of the six-degree-of-freedom serial robot with an offset wrist; C. obtain a pose increment dX of the expected pose of the terminal coordinate frame with respect to the approximate pose of the terminal coordinate frame through calculation, and obtain the synthetic error e of the pose of the terminal coordinate frame by using dX; D. judge whether the synthetic error e is within a reasonable error range or not; if the synthetic error is within the reasonable error range, stop and return the current numerical solution of the joint variable as an iterative numerical solution of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist; if the synthetic error e is not within the reasonable error range, obtain the joint variable through a jacobian matrix of the six-degree-of-freedom serial robot with an offset wrist, and put the new joint variable into step B and step C to conduct iterative calculation until the synthetic error e is within the reasonable error range or the maximum number of iterations is reached and then return a current numerical solution of the joint variable as the iterative numerical solution of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist.

2. The numerical method for obtaining the inverse kinematics of a six-degree-of-freedom serial robot with an offset wrist according to claim 1, further comprising before step A: establish a base coordinate frame by using an intersection point of a rotation axis of a first joint and a rotation axis of a second joint of the six-degree-of-freedom serial robot with a non-offset wrist as the original point, using a straight line where the rotation axis of the first joint is located as a Z axis of the base coordinate frame, and using a straight line where the rotation axis of the second joint is located as a Y axis; establish the coordinate frame {0} {1} {2} {3} {4} {5} and coordinate frame {6}, from the base coordinate frame towards the terminal coordinate frame of the robot in sequence according to the six-degree-of-freedom, wherein the coordinate frame {0} is the base coordinate frame, and the coordinate frame {6} is the terminal coordinate frame.

3. The numerical method for obtaining the inverse kinematics of a six-degree-of-freedom serial robot with an offset wrist according to claim 1, wherein step A comprises the following steps: A1. set the expected original point of the terminal coordinate frame {6} of the six-degree-of-freedom serial robot with an offset wrist as (X.sub.d, Y.sub.d, Z.sub.d), and set the expected pose matrix of the terminal coordinate frame {6} as R.sub.XYZ(.sub.d, .sub.d, .sub.d), and set the base coordinate frame {0} as a reference coordinate frame; use the pose of the terminal coordinate frame of the six-degree-of-freedom serial robot with an offset wrist as the pose of the terminal coordinate frame of the six-degree-of-freedom serial robot with a non-offset wrist, and use a method to obtain the analytical solution (.sub.1.sub.6) of the inverse kinematics of the six-degree-of-freedom serial robot with a non-offset wrist; calculate former three angles .sub.1, .sub.2 and .sub.3 of the inverse kinematics firstly: as the non-offset wrist meets Pieper criterion, and according to homogeneous transformation matrixes between coordinate frames: ${}_{}{}^{0}P_{4.Math..Math.\mathrm{org}}^{}=\left[\begin{array}{c}{g}_1.Math.\cos .Math..Math.{}_1^{}+g_2.Math.\sin .Math..Math.{}_1^{}\\ g_1.Math.\sin .Math..Math.{}_1^{}-g_2.Math.\sin .Math..Math.{}_1^{}\\ g_3\\ 1\end{array}\right]=\left[\begin{array}{c}{X}_d\\ Y_d\\ Z_d\\ 1\end{array}\right],.Math.\mathrm{wherein}.Math..Math.\{\begin{array}{c}{g}_1=d_4.Math.\sin \left({}_2^{}+{}_3^{}\right)+a_2.Math.\cos .Math..Math.{}_2^{}\\ g_2=d_2+d_3\\ g_3=a_2.Math.\sin .Math..Math.{}_2^{}-d_4.Math.\cos \left({}_2^{}+.Math.{}_3^{}\right)\end{array},$ obtain the analytical solution of .sub.1, .sub.2 and .sub.3 through solving the equations above, wherein ${}_{}{}^{0}P_{4.Math..Math.\mathrm{org}}^{}=\left[\begin{array}{c}{X}_d\\ Y_d\\ Z_d\\ 1\end{array}\right]$ is a vector of the non-offset wrist point, a.sub.2 is the distance from the Z axis of the coordinate frame {2} to the Z axis of the coordinate frame {3} along the X axis of the coordinate frame {2}, d.sub.2 is the distance from the X axis of the coordinate frame {1} to the X axis of the coordinate frame {2} along the Z axis of the coordinate frame {2}, d.sub.3 is the distance from the X axis of the coordinate frame {2} to the X axis of the coordinate frame {3} along the Z axis of the coordinate frame {3}, and d.sub.4 is the distance from the X axis of the coordinate frame {3} to the X axis of the coordinate frame {4} along the Z axis of the coordinate frame {4}; A2.obtain the last three angles .sub.4, .sub.5 and .sub.6 of the inverse kinematics of the six-degree-of-freedom serial robot with a non-offset wrist by using a transformation matrix R.sub.XYZ(.sub.d, .sub.d, .sub.d) which is based on the Euler angles XYZ(.sub.d, .sub.d, .sub.d):
R.sub.XYZ(.sub.d, .sub.d, .sub.d)=.sub.4.sup.0R|.sub..sub.4.sub.=0 6.sup.4R|.sub..sub.4.sub.=.sub.5.sub.=.sub.6.sub.=0 E.sup.6R.sub.XYZ(.sub.4, .sub.5, .sub.6), wherein .sub.4.sup.0R is a rotation matrix of the coordinate frame {4} with respect to the base coordinate frame {0}, .sub.6.sup.4R is a rotation matrix of the coordinate frame {6} with respect to the coordinate frame {4}, and .sub.E.sup.6R.sub.xyz is a rotation matrix of the coordinate frame {6} obtained after rotation of last three rotation axes; and obtain the matrix of .sub.E.sup.6R.sub.XYZ(.sub.4, .sub.5, .sub.6) through calculating the formula above, and when cos .sub.50, the analytic solutions of .sub.4, .sub.5 and .sub.6 can be obtained.

4. The numerical method for obtaining the inverse kinematics of a six-degree-of-freedom serial robot with an offset wrist according to claim 1, wherein step B comprises the following steps: B1. use the analytical solution obtained in step A .sub.1.sub.6 of the inverse kinematics of the six-degree-of-freedom serial robot with a non-offset wrist as the approximate solution of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist and the initial point for iteration; according to the matrix expression formula of the forward kinematics of the six-degree-of-freedom serial robot with an offset wrist: ${}_6^0.Math.T={}_1^0.Math.T\left({}_1^{}\right).Math.{}_2^1.Math.T\left({}_2^{}\right).Math.{}_3^2.Math.T\left({}_3^{}\right).Math.{}_4^3.Math.T\left({}_4^{}\right).Math.{}_5^4.Math.T\left({}_5^{}\right).Math.{}_6^5.Math.T\left({}_6^{}\right)=\left[\begin{array}{cccc}{n}_x& o_x& a_x& X_c\\ n_y& o_y& a_y& Y_c\\ n_z& o_z& a_z& Z_c\\ 0& 0& 0& 1\end{array}\right],$ wherein, the Euler angles $R_{\mathrm{XYZ}}^{}\left({}_c,{}_c,{}_c\right)=\left[\begin{array}{ccc}{n}_x& o_x& a_x\\ n_y& o_y& a_y\\ n_z& o_z& a_z\end{array}\right],$ obtain the approximate position (X.sub.c, Y.sub.c, Z.sub.c) of the original point of the terminal coordinate frame {6} and the approximate pose matrix R.sub.XYZ(.sub.c, .sub.c, Y.sub.c) of the terminal coordinate frame {6 }by solving the formula above.

5. The numerical method for obtaining the inverse kinematics of a six-degree-of-freedom serial robot with an offset wrist according to claim 1, wherein step C comprises the following steps: C1. obtain a position increment dP of the original point of the terminal coordinate frame {6} according to the expected position (X.sub.d, Y.sub.d, Z.sub.d) of the original point of the terminal coordinate frame {6} and the approximate position (X.sub.c, Y.sub.c, Z.sub.c) of the original point of the terminal coordinate frame {6}; C2. obtain a rotation increment dAng of the terminal coordinate frame {6} according to the rotation relation between the expected pose matrix R.sub.XYZ(.sub.d, .sub.d, .sub.d) of the terminal coordinate frame {6} and the approximate pose matrix R.sub.XYZ(.sub.c, .sub.c, .sub.c) of the terminal coordinate frame {6}; C3. combine the position increment dP of the original point of the terminal coordinate frame and the rotation increment dAng of the pose of the terminal coordinate frame into a differential motion dX of Cartesian space, which is represented as $\mathrm{dX}=\left[\begin{array}{c}\mathrm{dP}\\ \mathrm{dAng}\end{array}\right],$ and define the synthetic error of the pose of the terminal coordinate frame as e=idP.sup.2+dAng.sup.2.

6. The numerical method for obtaining the inverse kinematics of a six-degree-of-freedom serial robot with an offset wrist according to claim 1, wherein step D comprises the following steps: D 1. judge whether the synthetic error e is within a reasonable error range or not by the step C above, if the synthetic error e is within the reasonable error range, stop the iteration and return a current numerical value of the joint variable as the iterative numerical solution of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist; if the synthetic error is not within the reasonable error range, proceeding with step D2; D2.calculate d by using a jacobian matrix J of the six-degree-of-freedom serial robot with an offset wrist through dX=Jd, which is as follows:
d=J.sup.1dX, obtain the joint variable =+d, and put the new joint variable into step B1 to conduct iterative calculation until e is within the reasonable error range or the maximum number of iterations is reached, and then return the current numerical solution of the joint variable as the iterative numerical solution of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0030] The present disclosure is further described hereinafter with reference to the drawings and the embodiments.

[0031] FIG. 1 is a front view of a six-degree-of-freedom serial robot with an offset wrist according to the present disclosure;

[0032] FIG. 2 is a side view of the six-degree-of-freedom serial robot with an offset wrist according to the present disclosure;

[0033] FIG. 3 is a schematic diagram of a spatial coordinate frame established according to the six-degree-of-freedom serial robot with an offset wrist related in the present disclosure;

[0034] FIG. 4 is a flow chart of an algorithm of the present disclosure.

DETAILED DESCRIPTION

[0035] As shown in FIG. 1 and FIG. 2, a six-degree-of-freedom serial robot with an offset wrist according to the present disclosure comprises a base 0, a first joint 11 that connects to the base 0, a second joint 21 by a linkage 1, a third joint 31 that connects to the second joint 21 by a linkage 2, a fourth joint 41 that connects to the third joint 31 by a linkage 3 and linkage 4, a fifth joint 51 that connects to the fourth joint 41 by a linkage 5, a sixth joint 61 that connects to the fifth joint 51 by a linkage 6.

[0036] FIG. 3 shows a schematic diagram of spatial coordinates established according to the six-degree-of-freedom serial robot with an offset wrist of FIGS. 1 and 2. A base coordinate frame is established by utilizing the intersection point of a rotation axis of the first joint 11 and the rotation axis of the second joint 21 as an original point, a straight line where the rotation axis of the first joint 11 is located as Z axis of the base coordinate frame, a straight line where the rotation axis of the second joint 21 is located as a Y axis of the base coordinate frame. The coordinate frames are established from the base towards the terminal coordinate frame of the robot in sequence, i.e., coordinate frame {0} to coordinate frame {6}, the coordinate frame {0} is the base coordinate frame, and the coordinate frame {6} is the terminal coordinate frame.

[0037] As shown in FIG. 4, the inverse kinematics calculation method of a six-degree-of-freedom serial robot with an offset wrist comprises the following steps:

[0038] set a point where the tail end of a wrist of the six-degree-of-freedom serial robot with an offset wrist is located as (X.sub.d, Y.sub.d, Z.sub.d), i.e., the original point of the terminal coordinate frame {6},set the base coordinate frame {0} as a reference frame, set the pose matrix as R.sub.XYZ(.sub.d, .sub.d, .sub.d), and use the pose of the terminal coordinate frame of the six-degree-of-freedom serial robot with an offset wrist as the pose of the terminal coordinate frame of the six-degree-of-freedom serial robot with a non-offset wrist, and use a method for obtaining the inverse kinematics of the six-degree-of-freedom serial robot with a non-offset wrist, to obtain the analytical solution (.sub.1.sub.6) of the inverse kinematics of the six-degree-of-freedom serial robot with a non-offset wrist. Calculate the former three angles .sub.1, .sub.2 and .sub.3 of the inverse kinematics firstly; as the non-offset wrist meets the Pieper criterion, and according to transformation matrix of the coordinate frame {1 } with respect to the coordinate frame {0}, transformation matrix of the coordinate frame {2}with respect to the coordinate frame {1}, and transformation matrix of the coordinate frame {3 } with respect to the coordinate frame {2}:

[00006]${}^{0}P_{4.Math.\mathrm{org}}={}_1^0.Math.T\left({}_1^{}\right).Math.{}_2^1.Math.T\left({}_2^{}\right).Math.{}_3^2.Math.T\left({}_3^{}\right).Math.{}^{3}P_{4.Math.\mathrm{org}}=\left[\begin{array}{c}{X}_d\\ Y_d\\ Z_d\\ 1\end{array}\right],$ [0039] by arranging the formula above, it can obtain that

[00007]${}_{}{}^{0}P_{4.Math..Math.\mathrm{org}}^{}=\left[\begin{array}{c}{g}_1.Math.\cos .Math..Math.{}_1^{}+g_2.Math.\sin .Math..Math.{}_1^{}\\ g_1.Math.\sin .Math..Math.{}_1^{}-g_2.Math.\sin .Math..Math.{}_1^{}\\ g_3\\ 1\end{array}\right],.Math.\mathrm{wherein}.Math..Math.\{\begin{array}{c}{g}_1=d_4.Math.\sin \left({}_2^{}+{}_3^{}\right)+a_2.Math.\cos .Math..Math.{}_2^{}\\ g_2=d_2+d_3\\ g_3=a_2.Math.\sin .Math..Math.{}_2^{}-d_4.Math.\cos \left({}_2^{}+.Math.{}_3^{}\right)\end{array},\mathrm{and}$ [0040] .sub.1, .sub.2 and .sub.3 can be obtained through solving the equations above wherein,

[00008]${}_{}{}^{0}P_{4.Math..Math.\mathrm{org}}^{}=\left[\begin{array}{c}{X}_d\\ Y_d\\ Z_d\\ 1\end{array}\right]$

is a vector of the point where the tail end of the wrist is located; [0041] a.sub.2 is the distance between the Z axis of the coordinate frame {2} and the Z axis of the coordinate frame {3} along the X axis of the coordinate frame {2}; [0042] d.sub.2 is the distance between the X axis of the coordinate frame {1} and the X axis of the coordinate frame {2} along the Z axis of the coordinate frame {2}; [0043] d.sub.3 is the distance between the X axis of the coordinate frame {2} and the X axis of the coordinate frame {3 } along the Z axis of the coordinate frame {3}; [0044] d.sub.4 is the distance between the X axis of the coordinate frame {3 } and the X axis of the coordinate frame{4} along the Z axis of the coordinate frame {4}; use the transformation matrix R.sub.XYZ (.sub.d, .sub.d, .sub.d) with Euler angles XYZ(.sub.d, .sub.d, .sub.d) to obtain the last three angles .sub.4, .sub.5 and .sub.6 of the inverse kinematics of the six-degree-of-freedom serial robot with a non-offset wrist:

R.sub.XYZ(.sub.d, l.sub.d, .sub.d)=.sub.4.sup.0R|.sub..sub.4.sub.=0 6.sup.4R|.sub..sub.4.sub.=.sub.5.sub.=.sub.6.sub.=0 E.sup.6R.sub.XYZ(.sub.4, .sub.5, .sub.6), [0045] wherein .sub.4.sup.0R is a rotation matrix of the coordinate frame {4} with respect to the base coordinate frame {0}, .sub.6.sup.4R is a rotation matrix of the coordinate frame {6 } with respect to the base coordinate frame { 4}, and .sub.E.sup.6R.sub.xyz is a rotation matrix of the coordinate frame {6} obtained after rotation of the last three rotation axes; [0046] the solution process of the formula above is as follows:

.sub.E.sup.6R.sub.XYZ(.sub.4, .sub.5, .sub.6)=.sub.6.sup.4R|.sup.1.sub..sub.4.sub.=.sub.5.sub.=.sub.6.sub.=0 4.sup.0R|.sup.1.sub..sub.4.sub.=0R.sub.XYZ(.sub.d, .sub.d, .sub.d), [0047] from which it can be obtained that:

[00009]${}_E{}^{6}R_{\mathrm{XYZ}}^{}\left({}_4^{},{}_5^{},{}_6^{}\right)=\left[\begin{array}{ccc}{c}_5.Math.c_6& -c_5.Math.c_6& s_5\\ s_4.Math.s_5.Math.s_6& -s_4.Math.s_5.Math.s_6+c_4.Math.c_6& -s_4.Math.c_5\\ -c_4.Math.s_5.Math.s_6+s_4.Math.s_6& c_4.Math.s_5.Math.s_6+s_4.Math.c_6& c_4.Math.c_5\end{array}\right],$ [0048] wherein: c.sub.4=cos .sub.4, s.sub.4=sin .sub.4, c.sub.5=cos .sub.5, s.sub.5=sin .sub.5, and c.sub.6=cos .sub.6, s.sub.6=sin .sub.6,

[0049] AS .sub.1, .sub.2, and .sub.3 are obtained, .sub.4.sup.0R|.sub..sub.4.sub.=0 is obtained if .sub.4=0; therefore, .sub.4, .sub.5 and .sub.6 can be obtained through the transformation matrix R.sub.XYZ(.sub.d, .sub.d, .sub.d) with Euler angles XYZ(.sub.d, .sub.d, .sub.d) , and .sub.E.sup.6R.sub.XYZ(.sub.4, .sub.5, .sub.6)can be represented as the following form:

[00010]${}_E{}^{6}R_{\mathrm{XYZ}}^{}\left({}_4^{},{}_5^{},{}_6^{}\right)=\left[\begin{array}{ccc}{r}_{11}& r_{12}& r_{13}\\ r_{21}& r_{22}& r_{23}\\ r_{31}& r_{32}& r_{33}\end{array}\right],$ [0050] when cos .sub.50, it can be obtained that the analytical expressions of .sub.4, .sub.5 and .sub.6 are as follows:

[00011]$\{\begin{array}{c}{}_5^{}=\arctan .Math..Math.2.Math.\left(r_{13},\sqrt{r_{11}^2+r_{12}^2}\right)\\ {}_4^{}=\arctan .Math..Math.2.Math.\left(-\frac{r_{23}}{\cos .Math..Math.{}_5^{}},\frac{r_{33}}{\cos .Math..Math.{}_5^{}}\right)\\ {}_6^{}=\arctan .Math..Math.2.Math.\left(-\frac{r_{12}}{\cos .Math..Math.{}_5^{}},\frac{r_{11}}{\cos .Math..Math.{}_5^{}}\right)\end{array},$

[0051] The analytical solution .sub.1.sub.6 of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist can be obtained from the steps above, and the analytical solution is used as the approximate solution of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist and the initial point for iteration.

[0052] Use the approximate solution .sub.1.sub.6 of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist obtained through the calculating steps above and according to the matrix of the forward kinematics of the six-degree-of-freedom serial robot with an offset wrist:

[00012]${}_6^0.Math.T={}_1^0.Math.T\left({}_1^{}\right).Math.{}_2^1.Math.T\left({}_2^{}\right).Math.{}_3^2.Math.T\left({}_3^{}\right).Math.{}_4^3.Math.T\left({}_4^{}\right).Math.{}_5^4.Math.T\left({}_5^{}\right).Math.{}_6^5.Math.T\left({}_6^{}\right)=\left[\begin{array}{cccc}{n}_x& o_x& a_x& X_c\\ n_y& o_y& a_y& Y_c\\ n_z& o_z& a_z& Z_c\\ 0& 0& 0& 1\end{array}\right],.Math..Math.\mathrm{wherein}.Math.\text{:}$$.Math.{}_1^0.Math.T\left({}_1^{}\right)=\left[\begin{array}{cccc}\cos .Math..Math.{}_1^{}& -\sin .Math..Math.{}_1^{}& 0& 0\\ \sin .Math..Math.{}_1^{}& \cos .Math..Math.{}_1^{}& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 1\end{array}\right],.Math..Math.{}_2^1.Math.T\left({}_2^{}\right)=\left[\begin{array}{cccc}\cos .Math..Math.{}_2^{}& -\sin .Math..Math.{}_2^{}& 0& 0\\ 0& 0& -1& -d_2\\ \sin .Math..Math.{}_2^{}& \cos .Math..Math.{}_2^{}& 0& 0\\ 0& 0& 0& 1\end{array}\right],.Math..Math.{}_3^2.Math.T\left({}_3^{}\right)=\left[\begin{array}{cccc}\cos .Math..Math.{}_3^{}& -\sin .Math..Math.{}_3^{}& 0& a_2\\ \sin .Math..Math.{}_3^{}& \cos .Math..Math.{}_3^{}& 0& 0\\ 0& 0& 1& d_3\\ 0& 0& 0& 1\end{array}\right],.Math..Math.{}_4^3.Math.T\left({}_4^{}\right)=\left[\begin{array}{cccc}\cos .Math..Math.{}_4^{}& -\sin .Math..Math.{}_4^{}& 0& 0\\ 0& 0& -1& -d_4\\ \sin .Math..Math.{}_4^{}& \cos .Math..Math.{}_4^{}& 0& 0\\ 0& 0& 0& 1\end{array}\right],.Math..Math.{}_5^4.Math.T\left({}_5^{}\right)=\left[\begin{array}{cccc}\cos .Math..Math.{}_5^{}& -\sin .Math..Math.{}_5^{}& 0& 0\\ 0& 0& 1& d_5\\ -\sin .Math..Math.{}_5^{}& -\cos .Math..Math.{}_5^{}& 0& 0\\ 0& 0& 0& 1\end{array}\right],.Math..Math.{}_6^5.Math.T\left({}_6^{}\right)=\left[\begin{array}{cccc}\cos .Math..Math.{}_6^{}& -\sin .Math..Math.{}_6^{}& 0& 0\\ 0& 0& -1& -d_6\\ \sin .Math..Math.{}_6^{}& \cos .Math..Math.{}_6^{}& 0& 0\\ 0& 0& 0& 1\end{array}\right],.Math..Math.R_{\mathrm{XYZ}}^{}\left({}_c,{}_c,{}_c\right)=\left[\begin{array}{ccc}{n}_x& o_x& a_x\\ n_y& o_y& a_y\\ n_z& o_z& a_z\end{array}\right],$

[0053] wherein:

[0054] d.sub.5 is the distance between the X axis of the coordinate frame {4} and the X axis of the coordinate frame {5} along the Z axis of the coordinate frame {5};

[0055] d.sub.6 is the distance between the X axis of the coordinate frame {5} and the X axis of the coordinate frame {6} along the Z axis of the coordinate frame {6}.

[0056] The approximate position (X.sub.c, Y.sub.c, Z.sub.c) of the original point of the terminal coordinate frame and the approximate pose matrix R.sub.XYZ(.sub.c, .sub.c, .sub.c) are obtained through solving the formula above.

[0057] (X.sub.d, Y.sub.d, Z.sub.d) is the expected position of the original point of the terminal coordinate frame and R.sub.XYZ(.sub.d, .sub.d, .sub.d) is the expected pose matrix of the terminal coordinate frame. (X.sub.c, Y.sub.c, Z.sub.c) is the approximate position of the original point of the terminal coordinate frame and R.sub.XYZ(.sub.c, .sub.c, .sub.c) is the approximate pose matrix of the terminal coordinate frame, which can be obtained from the steps above. Assuming that the position increment dP of the original point of the terminal coordinate frame and the rotation increment dAng of the pose of the terminal coordinate frame are replaced with a differential motion dX of Cartesian space, which can be represented as:

[00013]$\mathrm{dX}=\left[\begin{array}{c}\mathrm{dP}\\ \mathrm{dAng}\end{array}\right],$

[0058] the position increment dP=[X.sub.dX.sub.c, Y.sub.dY.sub.c, Z.sub.dZ.sub.c], and

[0059] the rotation increment dAng is calculated as follows:

[0060] Let vector K(k.sub.x, k.sub.y, k.sub.z) be the vector respectively passing through the original point of the terminal coordinate frame; assuming that the terminal coordinate frame rotates by an angle around K axis at a constant speed, which rotates from the Euler angles XYZ(, , ) of the approximate pose to the expected pose angle XYZ(.sub.d, .sub.d, .sub.d) of the terminal coordinate frame, then:

[00014]$R\left(K,\right)=R_{\mathrm{XYZ}}^{.Math.-1}\left({}_c,{}_c,{}_c\right).Math.R_{\mathrm{XYZ}}^{}\left({}_d,{}_d,{}_d\right)=\left[\begin{array}{ccc}{k}_x.Math.k_x.Math.v+\cos .Math..Math.& k_x.Math.k_y.Math.v-k_z.Math..Math.\sin .Math..Math.& k_x.Math.k_z.Math.v+k_y.Math.\sin .Math..Math.\\ k_x.Math.k_y.Math.v+k_z.Math..Math.\sin .Math..Math.& k_y.Math.k_y.Math.v+\cos .Math..Math.& k_y.Math.k_z.Math.v-k_x.Math.\sin .Math..Math.\\ k_x.Math.k_z.Math.v-k_y.Math.\sin .Math..Math.& k_y.Math.k_z.Math.v+k_x.Math.\sin .Math..Math.& k_z.Math.k_z.Math.v+\cos .Math..Math.\end{array}\right]=\left[\begin{array}{ccc}{r}_{11}^{}& r_{12}^{}& r_{13}^{}\\ r_{21}^{}& r_{22}^{}& r_{23}^{}\\ r_{31}^{}& r_{32}^{}& r_{33}^{}\end{array}\right],.Math..Math.\mathrm{wherein}.Math..Math.v=1-\cos .Math..Math.,$

[0061] through calculating the formula above, dAng is obtained:

[00015]$\mathrm{dAng}=\frac{1}{2.Math.\sin .Math..Math.}\left[\begin{array}{c}{r}_{32}^{}-r_{23}^{}\\ r_{13}^{}-r_{31}^{}\\ r_{21}^{}-r_{12}^{}\end{array}\right].$

[0062] In order to describe the deviation(dP, dAng) between the current approximate pose of the terminal coordinate frame and the expected pose of the terminal coordinate frame, the synthetic error of the pose of the terminal coordinate frame is defined as follows:

e=dP.sup.2+dAng.sup.2.

[0063] According to the formula above and dP, dAng obtained above, judge whether the synthetic error e is within a reasonable range or not; if e is within a reasonable range, stop the iterative calculation, and return the current numerical solution of the joint variable as an iterative numerical solution of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist; and if the synthetic error e is not within the reasonable range, then the iterative calculation continues.

[0064] Use the jacobian matrix dX=Jd of the six-degree-of-freedom serial robot with an offset wrist, wherein J is represented by the following formula:

[00016]$J=\left(\begin{array}{c}{}_{}^0.Math.Z_i^{}\left({}_{}{}^{0}P_n-{}_{}{}^{0}P_i\right)\\ {}_{}{}^{0}Z_i^{}\end{array}\right),i=1,2,3,\mathrm{.Math.}.Math.,6,n=6,$

[0065] wherein, .sup.0Z.sub.i is the vector representation of the i.sup.th rotation axis Z in the base coordinate frame, and .sup.0P.sub.i is the vector representation of the original point of the coordinate frame {i} with respect to the base coordinate frame.

[0066] Through reversibility of the jacobian matrix, d=J.sup.1dX is obtained, so that the new joint variable =+d can be calculated, and the new joint variable is put into the matrix expression formula for obtaining the forward kinematics of the six-degree-of-freedom serial robot with an offset wrist, and the approximate position (X.sub.c, Y.sub.c, Z.sub.c) of the original point of the terminal coordinate frame and the approximate pose matrix R.sub.XYZ(.sub.c, .sub.c, .sub.c) of the terminal coordinate frame are updated:

.sub.6.sup.0T=.sub.1.sup.0T(.sub.1).sub.2.sup.1T(.sub.2).sub.3.sup.2T(.sub.3).sub.4.sup.3T(.sub.4).sub.5.sup.4T(.sub.5).sub.6.sup.5T(.sub.6).

[0067] Execute iterative calculation by repeating the steps to obtain the approximate pose and the synthetic error e above, until e is within the reasonable error range or the maximum number of iterations is reached, then return the current numerical solution of the joint variable as the iterative numerical solution of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist.

[0068] The algorithm skillfully uses the analytical solution of the inverse kinematics of the six-degree-of-freedom serial robot with a non-offset wrist as the approximate solution of the inverse kinematics of the six-degree-of-freedom serial robot with an offset wrist to obtain the approximate pose; the error between the approximate pose and the expected pose is obtained and an equivalent axial angle is used to represent the pose rotation increment between the approximate pose and the expected pose of the terminal coordinate frame, and the joint variable increment d is obtained according to the transformation relation between the differential motion dX of Cartesian space and the joint variable space motion d, i.e., the jacobian matrix J, so as to obtain the new joint variable =+d, which is used as a new iterative point. Through the iterative calculation of a formula for obtaining the forward kinematics of the six-degree-of-freedom serial robot with an offset wrist, the iterative point is continuously updated to approach the expected pose, so as to obtain the numerical solution of the inverse kinematics meeting the actual accuracy requirement finally. The algorithm has less calculation amount, faster convergence, higher efficiency, less burden to a robot controller and this method provides better real-time performance and can improve the efficiency.

[0069] The contents above are only preferred embodiments of the present disclosure, and the present disclosure is not limited to the embodiments. Any solution with the same technical effects by the same means shall all fall within the protection range of the present disclosure.

[0070] Without further elaboration, it is believed that one skilled in the art can, using the preceding description, utilize the present invention to its fullest extent. The preceding preferred specific embodiments are, therefore, to be construed as merely illustrative, and not limitative of the remainder of the disclosure in any way whatsoever.

[0071] In the foregoing and in the examples, all temperatures are set forth uncorrected in degrees Celsius and, all parts and percentages are by weight, unless otherwise indicated.

[0072] The entire disclosures of all applications, patents and publications, cited herein and of corresponding Chinese application No. 201710969628.6, filed Oct. 18, 2017are incorporated by reference herein.

[0073] The preceding examples can be repeated with similar success by substituting the generically or specifically described reactants and/or operating conditions of this invention for those used in the preceding examples.

[0074] From the foregoing description, one skilled in the art can easily ascertain the essential characteristics of this invention and, without departing from the spirit and scope thereof, can make various changes and modifications of the invention to adapt it to various usages and conditions.