A Novel Analytical Inverse Kinematics Method for SSRMS-Type Space Manipulators Based on the POE Formula and the Paden- Kahan Subproblem

Space manipulators which have a similar symmetrical structure with seven revolute joints, such as the space station remote manipulator system (SSRMS), can be called SSRMS-type space manipulators. The analytical inverse kinematics of an SSRMStype manipulator can be solved by locking a single joint; the locked joint (joint 1, 2, 6, or 7) can be determined by configuration analysis. Although widely used in establishing the kinematics of SSRMS-type manipulators, the Denavit-Hartenberg (DH) method has a singular problem when two adjacent joint axes are nearly parallel. To avoid this problem, this paper proposes a novel analytical inverse kinematics method for SSRMS-type manipulators based on the product of exponentials (POE) formula and the Paden-Kahan subproblem. Because of the symmetrical structure, an SSRMS-type manipulator degrades to two kinds of 6-degree-of-freedom (DOF) manipulators when locking a single joint (joint 1, 2, 6, or 7). The analytical inverse kinematics of these two kinds of 6-DOF manipulators is solved by combining the Paden-Kahan subproblems and geometric and algebraic methods, respectively. The proposed approach is not only singularity free compared with the traditional DH-based methods but also more accurate than the POE-based numerical solution. The simulation results verify the efficiency of the proposed approach.


Introduction
Space manipulators, such as the space station remote manipulator system (SSRMS) and the special purpose dexterous manipulator (SPDM), have played vital roles in on-orbit servicing (OOS) missions [1,2] such as assembly, capture [3,4], refueling, repair [5], and debris removal [6,7]. Compared with nonredundant space manipulators, 7-degree-of-freedom (DOF) space manipulators have advantages in obstacle avoidance, joint limit avoidance, operability enhancement [8], joint torque optimization, singularity avoidance, and joint failure [9,10]. There are various kinds of possible configurations for 7-DOF space manipulators. Among these configurations, the spherical-revolute-spherical (SRS) type may be the optimal version [11]. The positioning and orienting problems of an SRS-type arm can be addressed, respectively, and the analytical inverse kinematics can be easily solved subsequently [12]. However, for transportation, storage, and increasing workspace purposes, the joints of 7-DOF space manipulators are always designed with broad joint movement ranges. To avoid the structural interference problem, many 7-DOF space manipulators are designed with joint offsets. At the same time, their positioning and orienting problems cannot be, respectively, resolved, and the analytical inverse kinematics cannot be easily addressed.
The SSRMS type is a typical configuration with joint offsets. The SSRMS and the Chinese space station remote manipulator system (CSSRMS), including the core module manipulator (CMM) and the experimental module manipulator [13], are typical SSRMS-type manipulators. The SPDM and the European robotic arm (ERA) have similar structures with joint offsets of different sizes. In general, an SSRMS-type space manipulator has a symmetrical design, and its joints are arranged as follows: a 3-DOF shoulder (roll-yaw-pitch), an elbow joint (pitch), and a 3-DOF wrist (pitch-yaw-roll), where the three axes of the shoulder/wrist do not intersect at a common point.
It is well known that the inverse kinematics of 7-DOF serial manipulators can be divided into two categories: numerical solution and analytical solution [14]. The numerical solution is usually achieved by calculating the Jacobian matrix, such as the extended Jacobian method [15], the gradient projection method [16,17], and the weight least-norm method [18]. These methods are suitable for manipulators with different configurations as long as the Jacobian matrix equations have solutions. However, they not only are computationally expensive but also have accumulated errors. The analytical solution can be calculated by locking a single joint [19] or fixing the arm angle [12], it can directly solve all possible joint angles without knowing the initial configuration, and it is faster and more accurate than the numerical solution. Its disadvantage is that it is heavily dependent on the configurations of the manipulators. In comparison, the analytical solution is more desirable when solving the inverse kinematics of a specific 7-DOF manipulator.
For the analytical inverse kinematics of SSRMS-type manipulators, scholars have carried out related research. Based on configuration analysis, Xu et al. [9,20] proposed a novel approach to solve the inverse kinematics of an SSRMS-type manipulator where analytical and semianalytical solutions were calculated by locking a single joint. Yu et al. [21] proposed a method to solve the analytical inverse kinematics of a 7-DOF space manipulator with an offset wrist. Based on constructing a virtual manipulator with a spherical wrist, the analytical inverse kinematics of a 7-DOF manipulator was solved by establishing the relationship between the virtual manipulator and the actual manipulator. According to the above studies, it can be found that the DH method is employed to establish the kinematics of SSRMStype manipulators most of the time. Although the DH method is quite popular, it has many disadvantages. For instance, the DH method has a singular problem when the axes of two adjacent joints are nearly parallel [22]; this will cause discontinuity when it is employed in kinematics calibration of an SSRMS-type manipulator. Besides, the joint coordinates cannot be arbitrarily arranged when the DH method is employed. To overcome the singular problem, scholars had proposed many new methods, for instance, the MDH method [22], the S method [23], and the complete and parametrically continuous (CPC) method [24,25], but these methods have not been widely used in modeling the kinematics of the SSRMS-type manipulators.
Another idea to establish the kinematics of a manipulator is to employ the product of exponentials (POE) formula. Due to the smooth mapping between Lie Algebras se(3) and Lie Groups SE(3), the POE formula is singularity free and has inherent advantages in modeling the kinematics of SSRMStype manipulators [26][27][28]; for instance, it is needless to establish the joint coordinate systems compared with the DH method. Besides, it has been widely used in the kinematics calibration of industrial manipulators [29][30][31]. Zhao et al. [32] discussed the singular problem of the DH and POE methods in combination with a 6-DOF industrial robot with adjacent parallel joints. In our previous work, we carried out the on-orbit kinematics calibration of an SSRMS-type manipulator based on the POE formula and achieved considerable results [33]. Owing to these advantages, it is essential to establish a POE-based analytical inverse kinematics method for SSRMS-type manipulators.
In this paper, to avoid the singular problem of the DH method, we propose a POE-based analytical inverse kinematics method for SSRMS-type manipulators. Our approach not only avoids the singularity of the DH method but is also more accurate than the POE-based numerical solution. The remainder of this paper is organized as follows: in Section 2, the kinematics model of a typical SSRMS-type manipulator is established based on the POE formula. To determine the locked joint (joint 1, 2, 6, or 7), the joints are classified into three categories by configuration analysis. Due to the symmetrical structure, the SSRMS-type manipulator degrades to two kinds of 6-DOF manipulators when locking a single joint (joint 1, 2, 6, or 7). In Section 3, we reconstruct the kinematics models of these two kinds of new 6-DOF manipulators. Analytical inverse kinematics of these two new 6-DOF manipulators is solved by combining the Paden-Kahan subproblems, geometric methods, and algebraic methods. Simulations are carried out to verify the efficiency of the proposed method in Section 4. Section 5 presents the conclusions of this paper.

The Kinematics and Model Analysis
2.1. Singularity Analysis of the DH Method. When establishing the kinematics of a manipulator, the DH method is the most popular, and it has four parameters: twist angle α i−1 , link length a i−1 , joint angle θ i , and translation distance d i , to represent the relationship between two adjacent joints. However, when the adjacent joint axes of a manipulator are parallel or nearly parallel, the corresponding DH parameters are discontinuous, which leads to the singularity of the DH method. Related studies have demonstrated the singularity of the DH method [22][23][24][25]. For instance, as shown in Figure 1  International Journal of Aerospace Engineering joint coordinate system rotates a small angle βðβ ⟶ 0Þ about axis y i , axes z i−1 and z i will intersect at a common point r; the corresponding DH parameters will change from nominal values to actual values: α i−1 = 0°, a i−1 = 0, θ i = 0°, and d i = 0. It can be found that when the adjacent joint axes are parallel, small deviations of the adjacent parallel joint axes lead to huge changes in the DH parameters. In fact, due to the manufacturing and assembly errors, there is no guarantee that the adjacent joint axes are completely parallel. Therefore, due to the singularity of the DH method, it is not suitable for establishing the kinematics of manipulators with adjacent parallel joints, especially SSRMS-type space manipulators with three adjacent parallel joints and large structural size.

The Kinematics Model Based on the POE Formula.
SSRMS-type space manipulators, such as the SSRMS, CSSRMS, SPDM, and ERA, play key roles in various OOS missions. As shown in Figure 2, they all have seven revolute joints and a similar joint arrangement: shoulder (roll-yaw-pitch), elbow (pitch), and wrist (pitch-yaw-roll). Although they differ in structural sizes and joint offset, the configuration of SSRMS-type manipulators can be expressed in a unified form, as shown in Figure 3. Due to the smooth mapping between Lie Algebras se(3) and Lie Groups SE(3), the POE formula is employed to establish the kinematics of a typical SSRMS-type manipulator. As shown in Figure 3, we attach a base coordinate system fSg to the base link and a tool coordinate system fTg to the end  3 International Journal of Aerospace Engineering effector. The pose of the tool frame fTg with respect to the base frame fSg can be expressed with a homogeneous transformation matrix g. All these matrixes form a rigid motion group, and an element of this group can be expressed as g ∈ SEð3Þ where the 3 × 3 matrix R ∈ SOð3Þ is a special orthogonal matrix describing the coordinate axes of the tool frame fTg in the base frame fSg, and the 3 × 1 vector p ∈ ℝ 3 is the position vector of the origin of the tool frame fTg in the base frame fSg. Besides, each joint is associated with a twist ξ i ði = 1, 2, ⋯, 7Þ, the twist coordinate of the ith joint ξ i in the base frame fSg: where ω i is the unit vector along the ith joint axis and r i is the position vector located at any point of the ith joint axis. Therefore, the twist coordinate of the ith revolute joint ξ i in the base frame fSg satisfies the following constraints: where k·k represents the 2-norm of a vector. The twist coordinate of the ith joint axis ξ i can also be expressed in a matrix form b ξ i ∈ seð3Þ: where the symbol ∧ represents a map from a 6-dimensional (or 3-dimensional) vector space to Lie Algebra space se (3) (soð3Þ) and b ω i is the antisymmetric matrix of ω i and can be expressed as The forward kinematics of the SSRMS-type manipulator can be given by using the POE formula: where g ST ðθÞ and g ST ð0Þ represent the actual and initial pose of the tool frame fTg in the base frame fSg, respectively, and θ i , i = 1, 2, ⋯, 7, represents the position of the ith joint. e ξ∧ i θ i can be expressed as The kinematics parameters of the SSRMS-type manipulator are listed in Table 1.

Configuration Analysis.
After establishing the forward kinematics of the SSRMS-type manipulator by using the POE formula, the configuration analysis is required to determine which joint can be locked to solve the analytical inverse kinematics of the SSRMS-type manipulator. It is Table 1: The kinematics parameters of an SSRMS-type manipulator.

No. i
Axis direction ω i Position r i The initial pose of the end effector Middle group (3,4,5) End group (6,7) By analyzing the configuration of the SSRMS-type manipulators, it can be found that there are three adjacent parallel joints (joints 3, 4, and 5). If all these three adjacent parallel joints are not locked in the new 6-DOF manipulator, there exist analytical solutions. In contrast, if any one of them is locked, the analytical solution condition for the new 6-DOF manipulator is not satisfied. To facilitate analysis, the joints of an SSRMS-type manipulator are classified into three groups [9]: the front group (joints 1 and 2), the middle group (joints 3, 4, and 5), and the end group (joints 6 and 7), as shown in Figure 4. Then, the analytical inverse kinematics of an SSRMS-type manipulator can be solved only when a joint in the front or end group (joint 1, 2, 6, or 7) is locked.

Analytical Inverse Kinematics for Locking a Single Joint
When locking a joint in the front or end group (joint 1, 2, 6, or 7), the SSRMS-type manipulator degrades to a new 6-DOF manipulator; the analytical inverse kinematics of the new 6-DOF manipulator can be solved by using the Paden-Kahan problems. The three basic Paden-Kahan subproblems [34] are listed as follows: (1) Rotation about a single axis (2) Rotation about two subsequent axes (3) Rotation to a given distance Due to the symmetrical structure of the SSRMS-type manipulators, the cases for locking a front or end joint (joints 1 and 7 or joints 2 and 6) are similar. Hence, only the analytical inverse kinematics for locking a joint in the front group (joints 1 and 2) are detailed.
3.1. Solving θ 2 When the First Joint Is Locked. When the first joint is locked at a fixed position θ L 1 , the SSRMS-type manipulator degrades to a new 6-DOF manipulator. As shown in Figure 5, the forward kinematics of the new 6-DOF manipulator can be written in the following form: where ξ i ′ represents the ith joint twist coordinate of the new 6-DOF manipulator and ξ i ′ = Ad e ξ∧ 1 θ 1 ξ i ði = 2, 3,⋯,7Þ. Ad T represents the adjoint transformation associated with which transforms a joint twist from one coordinate system to another, and it can be written in the following form: Right multiplying both sides of equation (8) by g −1 ST ð0Þ e −ξ∧ 1 θ L 1 gives the following relationship: By right multiplying both sides of equation (11) by p a , the screw motion of point O 7 in joints 2, 3, 4, and 5 can be obtained: where p a and p e are the initial and final positions of point O 7 , respectively. As shown in Figure 6, p b is obtained by rotating p a around the 5th joint axis by θ 5 , p c is obtained by rotating p b around the 4th joint axis by θ 4 , p d is obtained by rotating p c around the 3rd joint axis by θ 3 , and p e is obtained by rotating p d around the 2nd joint axis by θ 2 .
According to the geometrical relationship shown in Figure 6, we can obtain the following equation set: p d can be obtained using equation (14). The screw motion of point O 7 in the second joint can be expressed as Hence, θ 2 can be obtained using Paden-Kahan subproblem 1 [34].

3.2.
Solving θ 1 When the Second Joint Is Locked. When the second joint of the SSRMS-type manipulator is locked at a fixed position θ L 2 , the SSRMS-type manipulator degrades to another new 6-DOF manipulator. As shown in Figure 7,

2⨯
The middle group 2⨯ The end group 2⨯ The front group 2⨯ L L  International Journal of Aerospace Engineering the forward kinematics of the new 6-DOF manipulator can be written in the following form: where ξ i ′ represents the ith joint twist of the new 6-DOF manipulator; ξ i ′ = Ad Right multiplying both sides of equation (16) by g −1 ST ð0Þ e −ξ∧ 2 θ L 2 gives the following relationship: By right multiplying both sides of equation (11) by p a , the screw motion of point O 7 in joints 1, 3, 4, and 5 can be obtained: where p a and p e are the initial and final positions of point O 7 , respectively. As shown in Figure 8, p b is obtained by rotating p a around the 5th joint axis by θ 5 , p c is obtained by rotating p b around the 4th joint axis by θ 4 , p d is obtained by rotating p c around the 3rd joint axis by θ 3 , p e is obtained by rotating p d around the 1st joint axis by θ 1 , and q 1 and q 2 are two arbitrarily different points on the first joint axis. According to the geometrical relationship shown in Figure 8, we can obtain the following equation set: According to equations (19) and (20), we can get p d can be obtained by solving the above equation set. According to Figure 8, the screw motion of point O 7 in the second joint can be expressed as Hence, θ 1 can be obtained by using Paden-Kahan subproblem 1 [34].
3.3. Solving θ 3 , θ 4 , θ 5 , θ 6 , and θ 7 . According to Sections 3.1 and 3.2, both θ 1 and θ 2 can be solved. Left multiplying both 7 International Journal of Aerospace Engineering sides of equation (11) by e −ξ∧ 2 ′ θ 2 and equation (18) by e −ξ∧ 1 θ 1 gives the following relationships: It can be found that equations (23) and (24) are similar. Hence, the problems of resolving the analytical inverse kinematics of these two 6-DOF manipulators become the problem of resolving the analytical inverse kinematics of a 5-DOF manipulator.
Until now, the analytical inverse kinematics of an SSRMS-type manipulator has been solved by using the POE formula and the Paden-Kahan problem. The resolution pro-cess is as follows: resolving the front joints ⟶ resolving the end joints ⟶ resolving the middle joints. As shown in Figure 9, when locking a front joint (joint 1 or 2) at a fixed position, two sets of solutions for the other (joint 2 or 1) can be solved by using the geometrical relationships and the Paden-Kahan subproblems. Then, for each set of the front joints (joints 1 and 2), two sets of the end joints (joints 6 and 7) and θ 3 + θ 4 + θ 5 can be solved by using algebraic methods. Finally, for every set of the front and end joints (joints 1, 2, 6, and 7), two sets of the middle joints (joints 3, 4, and 5) can be solved by using algebraic methods. Because of the symmetrical structure of the SSRMS-type manipulators, the solution flow is similar when locking an end joint (joint 6 or 7).

Simulation
To verify the efficiency of the proposed method, simulations including two inverse kinematics cases, two motion planning cases, and an on-orbit kinematics calibration case are carried out on the CMM which will be deployed on the core module of the Chinese space station. The structural parameters of the CMM are set as follows:    (1) When the first joint is locked at θ L 1 = 30 ∘ , we set The pose of the tool frame g ST ðθÞ in the base frame can be calculated by using the forward kinematics equation (8): (2) When the second joint is locked at θ L 2 = 45 ∘ , we set Then, taking equations (29) and (31) as known conditions, the analytical inverse kinematics of the CMM can be solved by using the proposed method in Section 3. The corresponding simulation results are listed in Tables 2 and 3, respectively. It can be found that the proposed method can solve all 8 sets of possible solutions when locking joint 1 or 2. Next, we solve the forward kinematics again to calculate the poses of the tool frame corresponding to these 8 sets of   Table 2 or 3 and calculate their corresponding position and orientation errors to evaluate the maximum position and orientation errors. Finally, it can be concluded that when the first joint is locked, the maximum position and orientation errors of the new 6-DOF manipulator are 0.025 mm and 6:7 × 10 −5°, and when the second joint is locked, the maximum position and orientation errors of the new 6-DOF manipulator are 0.021 mm and 6:9 × 10 −5°.
To further evaluate the computational accuracy of the proposed method, a linear trajectory and a spline curve trajectory for the CMM are employed. As shown in Figure 10, these two trajectories have the same initial pose T ini , and their target poses are T tar lin andT tar spl , respectively.
15001 and 9001 discrete poses on the linear and spline curve trajectories are sampled to calculate the corresponding joint angles by using the proposed method. For comparison, the POE-based pseudo-inverse method is also employed in the motion planning of these two testing trajectories, and it can be expressed in the following form: where _ θ denotes the joint velocity column vector, J + is the pseudo-inverse of the Jacobian matrix, and _ X denotes the velocity column vector of the end effector.
As shown in Figure 11, our proposed method is with higher computational accuracy in solving the inverse kinematics of the CMM than the POE-based pseudo-inverse method. Our proposed method improves the mean position

10
International Journal of Aerospace Engineering and orientation accuracy by approximately 42% and 33%, respectively. Besides, the computational time of the proposed method is also verified. We implement the proposed algorithm to solve the inverse kinematics of the linear trajectory in Figure 10, in VC++ environment (Windows 10 64-bit, Intel Core i5 9300H 2.4 GHz CPU, 8 GB DDR4 RAM, 500 GB SSD), and record the computation time of every sampling pose 800 times. When locking the first joint, the mean and maximum computational time of the proposed method is 4.513 μs and 5.571 μs, respectively. When locking the second joint, the mean and maximum computational time of the proposed method is 4.735 μs and 5.657 μs, respectively. The results show that our proposed method has sufficient computational efficiency to solve the analytical solutions of the SSRMS-type manipulators in practical applications.

On-Orbit Kinematics Calibration
Case. SSRMS-type manipulators work in the space environment most of the time. In addition to machining and assembly errors, the huge impact during the launching process and the deformation after long-term use will inevitably cause changes in structural parameters. Therefore, it is essential to conduct on-orbit kinematics calibration for SSRMS-type space manipulators. Due to the singularity of the DH method when two adjacent joints are nearly parallel, the DH method is not suitable for establishing the forward and inverse kinematics of the SSRMS-type manipulators with three adjacent parallel joints.
Although the nonsingularity of the POE formula has been well proven by scholars, to verify that the proposed method based on POE formula has avoided the singularity of the DH method, it is employed in the on-orbit kinematics self-calibration of the CMM. Different from the kinematics calibration of traditional industrial robots, no external measuring equipment is available when carrying out on-orbit kinematics calibration of a space manipulator, the stereo camera mounted on the end effector of an SSRMS-type manipulator is always employed to measure the absolute poses of the end effector [35], as shown in Figure 12. Limited by the measurement range of the stereo camera and the huge structural size of the CMM, the measurement configurations cannot be randomly generated in the joint space but generated in the operational space, and then, the corresponding joint angles for calibration can be calculated by inverse kinematics at every iteration step.
The on-orbit kinematics calibration of the CMM is implemented in MATLAB environment, the simulation flow is shown in Figure 13, the nominal and actual kinematics parameters of the SSRMS-type manipulator are set as Table 4, pose measurements are generated within the measurement range of the cameras mounted on the CMM's end effector, each pose measurement is injected with uniformly distributed noise within [-2.4, 2.4] mm and [-0.3, 0.3]°to emulate the real situation, and the corresponding joint angles for calibration are calculated using the proposed method instead of the POE-based pseudo-inverse method. In the parameter identification and compensation phases, the adjoint error model [30] is employed to establish the relationship between the kinematics parameter errors Δx and the pose errors of the end effector y i : where A i is the Jacobian matrix. Given m end effector pose measurements y i and the corresponding Jacobian matrix With enough pose measurements that make A satisfy the full column rank condition, the least square solution of the parameter errors Δx can be solved.
At the end of each iteration step, the kinematics parameters ξ a i and g a ST0 are updated using the following two equations [30], and they are regarded as the new nominal values in the next iteration step.
where η i and η st are variable twists. The iteration is not terminated until the norm of Δx converges to a preset small threshold ε.
To verify the accuracy of the kinematics calibration, sixty random configurations are employed to calculate the mean position and orientation errors before and after calibration. As shown in Figure 14, compared with the POE-based pseudo-inverse method, the proposed method improves the mean position and attitude accuracy by 9% and 7%, respectively. Due to the singularity of the DH method when the adjacent joints are close to parallel, the DH method is not suitable for establishing the actual kinematics of the CMM, and the traditional DH-based calibration algorithm cannot converge to stable values. Therefore, the calibration simulation results show that both the POE-based pseudoinverse method and our proposed method have avoided the singularity of the DH method and can establish the actual kinematics of an SSRMS-type manipulator with 3 adjacent nearly parallel joints, which owes to the smooth mapping between Lie Algebras se(3) and Lie Groups SE(3).

Conclusion
In this paper, we propose a novel POE-based analytical inverse kinematics method for the SSRMS-type space manipulators. Owing to the smooth mapping between Lie Algebras se(3) and Lie Groups SE(3), the POE formula is employed to solve the analytical inverse kinematics of an SSRMS-type space manipulator to avoid the singular problem of the widely used DH method. To determine the locked joint (joint 1, 2, 6, or 7), the joints of an SSRMS-type manipulator are divided into three categories-the front group, the middle group, and the end group-by configuration analysis. Because of the symmetrical structure, an SSRMS-type space manipulator degrades to two kinds of 6-DOF manipulators when locking a front or end joint. The analytical inverse kinematics of these two 6-DOF manipulators is solved by combining the Paden-Kahan subproblems, geometric methods, and algebraic methods. Finally, to verify the efficiency of the proposed method, it is employed in two analytical inverse kinematics cases, two motion planning cases, and an on-orbit kinematics calibration case. The simulation results show that the proposed method can not only avoid the singularity of the traditional DH method but also are with higher calculation accuracy than the POE-based pseudo-inverse method.
In our future work, we will further examine the singularity configuration when solving the analytical inverse kinematics. Optimizing the motion planning algorithms for the SSRMS-type manipulators by determining the position of the locked joint will also be a research direction.

Data Availability
The data used to support the findings of this study are included within the article.

Conflicts of Interest
The authors declare that they have no conflicts of interest.

12
International Journal of Aerospace Engineering