Motion Planning of Ground Simulator for Space Instable Target Based on Energy Saving

: In order to achieve the high-precision motion trajectory in ground experiment of space instable target (SIT) while reducing the energy consumption of the motion simulator, a robot motion planning method based on energy saving is proposed. Observable-based ground robot motion experiment system for SIT is designed and motion planning process is illustrated. Discrete optimization mathematical model of energy consumption of motion simulator is established. The general motion form of the robot joints in ground test is given. The optimal joint path of motion simulator based on energy consumption under discontinuous singularity configuration is solved by constructing the complete energy consumption directed path and Dijkstra algorithm. An improved method by adding the global optimization algorithm is used to decouple the coupled robot joints to obtain the minimum energy consumption path under the continuous singularity configuration of the motion simulator. Simulations are carried out to verify the proposed solution. The simulation data show that total energy saving of motion simulator joints adopting the proposed method under the condition of non-singularity configuration, joints coupled motion with continuous singularity configuration, and coexistence of non-singularity path and continuous singularity path are, respec-tively, 72.67%, 28.24%, and 62.23%, which proves that the proposed method can meet the requirements of ground motion simulation for SIT and effectively save energy.


Introduction
As unstable targets in space continue to increase and become more threatening, maneuvering tasks, such as on-orbit removal, have become more and more urgent [1][2][3][4]. Space instable targets are objects to be manipulated in orbit, such as satellites and space debris, in space operation missions. The movement of SIT is complicated, and their manipulation task is also difficult [5]. In order to verify the feasibility of the space control method and optimize the control algorithm, it is necessary to conduct ground manipulation simulation experiments of SIT [6,7]. The primary task of the ground simulation test is to accurately simulate the relative spatial movement and relative attitude motion between SIT and space operation aircraft.
A serial robot with multi-DOF, widely used in areas such as grinding, polishing, and parts processing, has the advantages of a large range of pose motion and high repeat positioning accuracy. It is an effective method to achieve the given motion of SIT in ground test platform through using six-DOF industrial robot [8,9]. During simulating given space movement in the ground test, it is necessary to optimize the robot path on the premise that the simulation accuracy meets the requirements, and at the same time ensure the dynamic characteristics of the simulator. With the increasingly tense global energy situation, the optimization and control of robot energy consumption has gradually become an important goal in industrial systems [10].
Many motion planning methods and trajectory optimization algorithms have been proposed and studied to achieve smooth movement, obstacle avoidance, and flexible contact collisions in different application scenarios [11][12][13][14][15]. Interpolation, curve fitting the trajectory of the robot joint angular displacement or the cartesian pose of the robot end target position, generally including polynomial curve, Bézier spline, and NURBS curve, is adopted to smooth the path in robot motion planning to reduce acceleration shock and movement jitter [16,17]. In high-precision machining of parts, Pavel [18] adopts the analytic geometry method on the basis of the symmetries in the Euclidean plane and achieves a trajectory error of less than 5% in experiments, which proves the feasibility of generating complex processing motion for robot. Modelled predictive control approach, based on vision system and linear interpolation, is applied for the task of ball-catching to generate robot real-time trajectory [19]. Aiming at the multi-path optimization issue, R.A. [20] uses genetic algorithms to search for the shortest path without collision, which connects preset points in the robot's task space. Chaos enhanced accelerated particle swarm optimization (CAPSO) is developed by Taghavifar [21] to plan the optimal path of a wheeled robot, and ensures the shortest path from the starting point to the target position by avoiding collision with static and dynamic obstacles. The improved algorithm combining Dijkstra algorithm or A* algorithm with other algorithms can effectively perform the shortest path search, which has the advantages of a small calculation and fast convergence speed, and is used in varieties of scenarios of robot operations [22][23][24][25][26].
In terms of optimal motion planning combined with minimal energy consumption and robot operation tasks, many implementation methods have emerged. One DoF mechatronic system is studied by Giovanni to reduce the energy consumption and the minimum energy consumption conditions have been found in a closed form taking into account the chance of recovering the braking energy [27]. A robot configuration selection method based on the lowest power consumption is studied by Abdullah and verified through a straight-line scenario and a square-path scenario, which achieve the goal of optimizing the calculated joint configuration [28]. Luo uses Chebyshev interpolation points and direct iteration method to plan the trajectory of industry robotic manipulators for energy minimization, and the method is applied to the three-link plane mechanism, which obtained optimal joint angles under the premise that the start and end points of the movement are known [29]. A task-related analysis is used by Scalera to enhance the energy efficiency of a 4-DOF parallel robot [14]. Carabin adopts the trajectory planning mean, which is on the basis of electromagnetic field model and the derivation of the energy formulation, to realize the minimum-energy consumption of the Cartesian robot [30]. Bitar proposed a method of obtaining energy-optimized trajectory planning while managing obstacle constraints for ASVs [31]. The enhanced bacterial foraging optimization algorithm is researched by Abbas and applied to path planning in two-dimensional motion space [32]. An artificial bee colony algorithm and a genetic algorithm is studied to optimize the working path length and welding time of the welding robot, which minimizes the total amount of robot joint movement [33,34]. Sathiya [35] takes the execution time and execution tasks of mobile robots as the optimization goals, and proposes the application of evolutionary algorithms to achieve multi-objective optimization in trajectory planning. Under the conditions of kinematics and dynamics constraints, Amruta [36] presents the EMOTLBO approach to get the best joint motion that defines the weld path, which effectively enhances the efficiency of the robot in operation tasks. Wang [37] applies the improved whale optimization algorithm to obtain the constant velocity motion of the robot end effector, which improves the processing efficiency and quality of the grinding task. Under the working conditions of the power generation mode and the energy exchange, Christian [38] proposes an energy-based robot model based on B-spline functions and realizes the optimization of energy between different joints. Aiming at the robot PTP motion planning without collision, Francisco [39] uses a third-order polynomial method to fit the joint space motion, which obtains an approximate optimal time trajectory and considers the constraint of energy consumed.
Inspired by the above-mentioned research works and considering the particularity of ground simulation of space instability targets, there are still several issues that need to be studied and solved. Firstly, it is hard to completely avoid robot singularity configuration due to the large coverage of the motion range in ground simulation experiment of SIT [40,41], which makes the motion optimization with singular position particularly important. Secondly, although the intelligent optimization algorithm based on the global search is sufficiently robust to acquire the optimized path, the optimization result is more likely to converge prematurely in the process of processing multi-objective optimization [42], which will cause the problem of difficulty in motion reproduction and simulation validity when the robot exists with different error sources. Furthermore, there may be multiple optimal motion paths or cases where the sub-optimal path is better than the optimal path on certain evaluation functions and the path search algorithm needs to be improved accordingly. Finally, the spin motion of SIT on orbit is not restricted, that is, it can be rotated infinitely, in contrast to the limited simulation of rotational motion subjected to robot structure. How to achieve the simulation of spin motion as much as possible with minimal energy consumption through the combined motion and motion optimization of the robot joints is also an issue of motion planning of the SIT movement simulation robot, especially when singularity causes multi-axis coupling.
Focusing on the characteristics of ground motion simulation of SIT and the mentioned challenges above, the study proposes an optimization method of simulator motion path including a singularity configuration based on energy saving. Ground motion planning of SIT is performed through motion mapping and path optimization of motion simulation robots. The local optimal solution is avoided through path search algorithm on the basis of constructing a complete directional motion path of energy consumption. The versatility of ground motion simulation of SIT is achieved using the presented approach and it enhances the authenticity of the ground experiment and reduces the energy consumption of the ground system. The rest of the paper is organized as follows: ground motion planning method of SIT and motion optimization model of SITMSR are described in Section 2 and Section 3, respectively; solution to the optimal joint trajectory of SITMSR and results obtained by a series of simulations are presented in Section 4 and Section 5, respectively; and conclusions are drawn in Section 6.

Ground Motion Simulation State Description of SIT
The motion of SIT in space is in the free-floating state, which can be decomposed into attitude movement and centroid translation movement. The attitude general form of SIT is spin movement around the axis of maximum inertia and the nutation angle will change with the pose controllability in the close-range orbit and attitude tracking stage and the final approach stage according to space dynamics [4]. Position movement is the threedimensional space motion obtained after task-level trajectory planning. Therefore, the movement of the SIT is a motion of 6-DOF in space. The motion simulation state of SIT in the ground experiment is described in Figure 1. To ensure the effectiveness of the ground test, the angular velocity and translation velocity are required to meet the given precision. The expressions of ground simulation motion accuracy of SIT are as follows: Through the analysis of the simulation state of the ground motion of SIT, it can be seen that the state exists where projection of the centroid velocity and angular velocity of the ground motion in GS GS GS GS O X Y Z is not zero, which means 6-DOF motion of SIT should be simulated in the ground test. In other words, velocity motion and angular velocity motion of SIT in the ground test possess the generality of space motion.

Composition and Motion Planning Process of SITMSRS
Compared with parallel robot and cooperative manipulator, 6-DOF industrial robots have the characteristics of large motion range, high load capacity, and high motion accuracy. At the same time, they have the advantages of small envelope space, high technology maturity, and low cost over a truss robot. Therefore, the 6-DOF industrial robot system is chosen to simulate the motion of SIT in the ground experiment. The space instable target motion simulation robot system (SITMSRS), mainly composed of a 6-DOF industrial robot (SITMSR), the space instable target simulator (SITS), a motion control computer, and electrical system, is designed to achieve the planned positions and attitudes given above. SITS is installed on the end flange of SITMSR, which equipped with capture ring, on-board analogue elements, and six-axis force/torque sensor. The capture ring is the capture object for the space manipulator to carry out takeover tasks during the ground test. The capture ring and on-board analogue elements are used as recognition features for motion parameter estimation and 3D reconstruction. The six-axis force/torque sensor can measure the forces and moments applied on SITS. SITS need to have the same properties as the real SIT, including the width and height of the SIT, the geometric size and material of the capture ring, etc. The motion control computer is used to realize the planning and optimization of the motion trajectory and interact with the robot system. The composition of designed SITMSRS is shown in Figure 2. The spatial movement of the centroid of the SIT and attitude movement, including precession, nutation, and spin planned by the motion control computer, can be achieved by joint compound motion of SITMSR to form the general motion form of SIT. It is hard to control angular velocity of SITMSR joints directly. Nevertheless, the spatial position GS IT r and attitude GS IT P can be obtained by integrating the velocity and angular velocity of the SIT in the ground simulation experiment. Thus, the ground motion of the SIT is simulated by controlling the angular displacement motion of the SITMSR joints.
Since there are multiple solutions for the angular displacement of robot joints under the same position and pose of the SIT, it is necessary to optimize the joint movement of SITMSR to save energy consumption. The input to robot motion control is based on discrete point planning, and when the number of discrete points Ndis is constant, the motion planning of the SIT is transformed into the path planning and optimization problem of the SITMSR joint motions. The motion planning flowchart of SITMSRS is indicated in Figure Figure 3. The motion planning flowchart of the SITMSRS.
From Figure 3, the joint path optimization of the SITMSR based on energy saving is the crucial part of the above motion planning process and needs to be further studied in detail.

Mathematical Model of Joint Trajectory of the SITMSR
Define qi (i=1,2,···,6 as generalized coordinate of SITMSR and SITMSR kinematics is modeled by the D-H method. Suppose 0 0 0 OX Y Z is the base coordinate system and OX Y Z i i i is the joint coordinate system. The definition of DH parameters and joint coordinate systems of SITMSR are shown in Figure 4. = q ,q ,q ,q ,q ,q ,a ,a ,a ,d ,d , , , , where 0 6 T is the homogeneous transformation matrix of SITMSR and s T is the position and attitude matrix planned of SIT in the ground experiment. a1, a2, a3, d1, and d4 are the geometric parameters of the coordinate transformation and α1, α3, α4, and α5 are the rotation angles of the joint coordinate system using D-H method as shown in Figure 4. The value of qi, which is also the joint angular displacement of each axis of the SITMSR, can be obtained by solving Equations (4)- (7).
The motion control of the SITMSR is based on the displacement control of the discrete points in the ground experiment. The angular displacement is discretized based on the dispersing of the given trajectory. The approach of polynomial fitting is used to carry out the continuity of the discrete path of the SITMSRS.

Discrete Optimization Mathematical Model of Energy Consumption
SITMSR need to accurately simulate the given motion of SIT firstly, which meet the given motion simulation error. Secondly, it is necessary to ensure the total energy consumption of the SITMSR axes is low during the simulating trajectory. Additionally, in order to avoid violent measurement data jitter of the six-axis force/torque sensor measurement values, the smoothness of the movement of the SITS is also required. Based on the above considerations, the objective function in the mathematical model of motion optimization of SITMSR can be described as follows: where p t is the given motion simulation time, i, j q  is the joint angular velocity, Nq is the number of SITMSR joint axes, Ndis is the number of discrete points of a given trajectory of SITS, i, j q  is the angular velocity of qi at the discrete point j, and qi, j T is the joint drive torque when the angular displacement is i, j q .
The constraint conditions corresponding to above objective function are: where IT ITa M and IT IT M are the values of the six-axis force/torque sensor in SITS when the trajectory is given and when the trajectory is planned, respectively. M γ is the allowable value of relative error between IT ITa M and IT IT M , traj N is the number of discrete points from the joint path to the joint trajectory. imax q  and imax q  are the maximum allowable values of the joint angular velocity and the angular acceleration, respectively. SITMSR dynamics equations in discrete form using Newton-Euler algorithm is presented in Equation (10).
where i,j f is the interaction force between link i and link i + 1 at the discrete point j of the trajectory, i,ci r is the position vector of the gravity center of link i in the framework of OXiYiZi , i,j ω is the angular velocity of the link i, i I is the inertial tensor matrix and i,j a is the linear acceleration at the end point of the link i.
SITMSR is considered as rigid bodies in the process of dynamic modeling, and the effects of viscous force are ignored. The measurements for six-axis force/torque can be obtained by solving the dynamics model of SITS based on Newton-Euler algorithm. The expressions of i,j f , i,ci r , i,j ω , and i,j a are obtained in [28]. 7,j f and 7, j τ in the framework of OX6Y6Z6 can be written as: where GSITS is the gravity of SITS.

The General Form of Joint Motion Trajectories
For the planned motion of SITS is complex and diverse, it leads to the generality of the motion of SITMSR. As a result, the compound situations of singularity configuration and non-singularity configuration of SITMSR should be considered. The singularity can be judged by the full rank of Jacobian matrix, which can also be quickly determined by using three singular factors extracted from the joint expressions of the Jacobian matrix. Another method is to judge whether there is a problem of DOF degradation. The general forms of the SITMSR joint motion trajectories are seen in Figure 5. To facilitate the solution of the optimal trajectory, the joint movement above can be segmented into paths with non-singularity configuration (Stage 2 and Stage 4 in Figure 5) and paths with singularity configuration. Paths with singularity configuration are divided into trajectory with multipoint singular position (Stage 1 in Figure 5) and trajectory with continuous singular position (Stage 3 in Figure 5). Although joint qi and qk have constraints, respectively, even under the premise of coupled angular displacement constraints, there are infinite sets of solutions for their joint angular displacements. If the robot continues to run in this condition, then the robot is in the continuous singularity configuration as shown in Stage 3. We solved each segment of the optimal path separately, and, finally, merged the optimal solutions of the multiple segments to obtain the optimal trajectory of SITMSR under the general path form of SIT in the ground experiment.

Optimal Selection of Ndis
From the motion planning of SITMSRS and motion optimization model based on energy saving, it can be seen that Ndis is closely related to the simulation accuracy of a given motion trajectory. The ground motion simulation error of SIT is large when Ndis is small, which, obviously, does not meet the simulation requirements. This is because when the motion of SITS is determined, the fewer the number of discrete points, the greater the deviation of the motion trajectory between two adjacent discrete points, and the greater the error between the corresponding velocity and angular velocity and the ideal state. The trajectory simulation precision will continue to improve as Ndis continues to increase and it is paradoxical that the total number of paths will also increase exponentially leading to lower efficiency in solving the optimal trajectory. Motion planning of the given trajectory of SITS can be adopted to select the optimized Ndis, the approach of polynomial fitting is used to carry out the continuity of the discrete path of SITMSRS. Optimal selection process of Ndis is depicted in Figure 6. Nd0 represents iteration initial value of Ndis. Approximate solutions of joint angular displacements at finite intermittent singular points can be achieved by choosing the appropriate Ndis through the above approach, which reduces the complexity of path optimization at a certain level.

Complete Path Solution Space Generation of Joints Power
In the compound movement of each joint of SITMSR to simulate the planned motion of SIT, there are several sets of solutions at each discrete point to satisfy a given position and attitude. Complete energy consumption path solution space is needed to generate by expanding the solution of each joint at each discrete path point into a full solution matrix. Figure 8 illustrates the generation process of complete path solution space of joint power.
where Nisol is the number of solutions of i, j q and Nsol is the maximum value of Nisol at all path discrete points of SITMSR.

Directed Motion Paths Construction of Energy Consumption
Constructing a networked directed non-circular motion paths of joints energy consumption of SITMSR, consisting of path nodes, directed paths, and path weights, is carried out to implement the path optimization of SITMSR on the basis of generating complete path solution space of joints power. The path node of energy consumption takes each vector of the path solution space as a node, and all nodes in the path are numbered in order. In addition, the start node and the end node are added to form a complete path to be optimized. The average of the product of the path discrete time and the sum of the energy consumption is used as the weight of the path between the two nodes except the start node and the end node. The nodes in the first row have one path weight to the initial node and Nsol path weights to other nodes. The node in row Ndis has one path weight to the destination node and Nsol path weights to other nodes. Aside from the first and last two rows, the number of input path weights and output path weights of other nodes are both Nsol. Figure 9 constructs the directed energy consumption path without continuous singularity of SITMSR. Figure 9. Construction of directed energy consumption paths without continuous singularity of SITMSR.

Solution to Optimal Energy Consumption Path with Discontinuous Singularity
From the building of directed energy consumption path, solving the optimal path of SITMSR is searching the shortest path starting from the initial node to the end node, which meet the requirements of Equations (8) and (9). The node vectors and the corresponding weight vectors need to be established first, and then get the shortest path from the created directed motion paths through the Dijkstra algorithm. Other shortest path or secondshortest path can be achieved by repeatedly calling the Dijkstra algorithm on the basis of modifying the weight vector. The solution to the optimal path based on energy saving without continuous singularity is described as follows: Step 1: Calculate the starting node vector Ini S from node (NdisNsol+1) to the Nsol nodes in the first row, end node vector Ini t , and corresponding path weight vector Ini L .
Step 2: Compute the starting node vector Inik S between node k located from row one to row (Ndis-1) and the next row of nodes, end node vector Inik t , and weight vector k L . Here, The expression of Inik S , Inik t , and k L are given as follows: In Equations (18)- (25), the definition of 1k n is to obtain path discrete point sequence number 1k m where the node k is located. k1 n is introduced to get the row number and the column number of the solution vector where node k is located. Step Step 4: Form the complete start node vector Int S , end node vector Int t , and weight vector Int L of SITMSR motion path to be optimized. Step 6: Traverse the path weights of other nodes except the node paths related to the initial node and the end node of 1 π in Step 5 one by one, and individually assign a new weight of +∞ . Update the path weight vector Int L before solving the new directional motion path and remove duplicate path solutions to achieve the vector of the node solution number, the path weight vector, and the path value of the other shortest paths.
Step 7: Reverse coding to acquire SITMSR joint angle matrix

Solution to the Optimal Singular Path
Since the finite number of discontinuous singular points can be solved by trajectory planning, the main focus here is on the condition of the continuous singularity configuration. When SITMSR is in the singular motion path, theoretically, there are countless solutions because of the coupled motion of multi-axis motion. Singularity configuration makes it difficult to directly solve the optimal motion path based on energy saving by way of the directed node path. The reason is that when there is no singular position, the path weight between the nodes of the SITMSR path network is unique. However, there are countless path weights between the nodes at each discrete location in the singular path.
When solving the optimal energy consumption path of SITMSR with multipoint singular position and continuous singular position in given trajectory of SITS, the joint angular displacements of the singular position at the discrete point N can be set as unknown variables to be optimized. Assume that the joint angle of the coupled motion axis caused by the singularity configuration is ( ) where 1,N,mi,ni Q is the (10mi+ni)-th solution of Di,N q at path discrete point N.
It is hard to solve the above situation by directly using the path optimization method under the condition of non-continuous singularity position. An approximate method is to assume that the angular displacements of the coupled joints have motion constraints, and then decouple the joint motions at singular positions. The minimum energy consumption path is obtained by adjusting the coupling coefficient and using the solution algorithm for non-continuous singular configuration. The linear relationship of the coupled joints can be expressed as: where i λ is the coupling factors of Di,N q .
Another feasible method is to take the objective function, constraint function, and constraint function of the coupling joint as the fitness function of the intelligent optimization algorithm to solve the above problems. The problem with this method is that when Di and Ndis are small, the optimal motion path of SITMSR can be achieved. Nevertheless, the convergence process of the above optimization method is longer and the optimization fails in the case where Ndis takes a larger value. Therefore, it is necessary to improve the path optimization algorithm for the non-singular condition, so that it can handle path optimization at the SITMSR singularity configuration.
When generating the complete path solution space of energy consumption with singular positions, take the joint angle constraint relationship in equation (32) of the coupled joint as a whole, and combine with the other uncoupled joint angle displacements to form the whole path solution space. When SITMSR is in a path with singularity, dis N ,m,n q is defined as: The optimal path of the uncoupled joints and the optimal path of the constraint relationship of the coupled joints QOpt π can be acquired according to the first seven steps when optimizing with non-singular paths on the basis of generating the complete solution space of the singular path and the corresponding directional motion path on the basis of solving the optimal path of the coupling joints, which form Opt π together.
To solve the optimal motion path of each coupled joint angle, the decoupling mathematical optimization model needs to be additionally introduced for secondary optimization to gain the coupled joint angle displacements for any value of Di,N q to get the minimum value of objective function f2 under path constraints QOpt π . It is an optimization problem with linear constraints, which can be expressed as: where 1,N Q is the value at the discrete point N of QOpt π and Di,N P is the power of Di,N q .
The optimal path of each coupled joint under the singular path can be sought through applying the global optimization algorithm to decoupling the mathematical optimization model above. The optimal energy consumption trajectory of the uncoupled joint can be obtained according to Step 8 in the solution to the optimal path without singularity and the path fitting based on the polynomial difference.

Results
To verify the effectiveness of the proposed method, optimization solutions based on energy saving of SITMSR using KuKa KR60 programmed by Matlab introduced in previous sections are tested under different movements given in the SIT ground motion experiment. Three examples are shown in this section. The path in the first one is under the condition of discontinuous singularity configuration. Example 2 reveals the energy optimization results of coupled joints using three different solution methods under continuous singularity configuration of SITMSR. Example 3 shows the energy consumption of continuous singularity configuration and non-singularity configuration under the premise of given trajectory. The parameters and corresponding values required for examination are summarized in Table 1.  in this example, respectively. The initial value of nutation angle is −3°. In order to make full use of the robot's reachable range of motion and to simulate the relative motion between the SIT and the aircraft operation in space as realistically as possible, the motion of OIT in SITS is set as spiral movement. Other forms of positional movement can also be set, which do not affect the verification of the proposed method. On average, the total time taken to complete the motion planning is 3.95 s in Example 1. Figure 10 shows the trajectory of SIT in the ground experiment in Example 1. The above given trajectory is discretized by setting different values of Ndis, and the cubic polynomial is used to fit and discretize the joint paths of SITMSR. Figure 11a shows the relationship between velocity and angular velocity errors and Ndis of SITS in motion planning. Force and torque error measured by six-axis force/torque sensor between the given motion and the planned motion are given in Figure 11b. As can be seen from the above curves, compared with path πSL, although the optimal path πOpt rotates 23.43% more in the joint angular displacement, its energy consumption is saved by 72.67%.
There are also differences in joints driving the torque of SITMSR on the above paths. Figure 13a,b, respectively, present the performance of joints driving the torque in the planned motion path of πOpt and πSL. The maximum joint driving the torque when the SITMSR moves along the planning path πOpt is only 63.84% of that along path πSL. In terms of smoothness of motion, the jitter of the robot joints under the path πOpt and path πSL are shown in Figure 14a,b. The maximum jitter value of SITMSR when moving along path πSL is 0.034 rad/s 3 , which occurs at joint 4, while the maximum jitter value of path πOpt is 0.162 rad/s 3 , which also happens at joint 4. The joint jitter along the optimal path of energy consumption is, obviously, larger than that along the shortest path, which means that the joints are also subject to greater impact, but its value is within the range of the SITMSR. The comparison results of path πOpt and path πSL in Example 1 are summarized in Table 2. In this example, the condition where the wrist joint of SITMSR is in continuous singularity configuration, that is, the rotation axis of robot joint four and joint six is collinear. Under the precondition that the end trajectory of SITS is given, the trajectory of other joints is determined, except for q4 and q6. Therefore, optimizing the joint space path of the coupling joint is focused to ensure optimal energy consumption. The precession angular velocity and nutation angular velocity of SITS are set as zeros and the angular velocity of spin of SITS is IT2 ψ  . The constraints of the angular displacement and the angular velocity of the coupled joints are shown in Figure 15a. The trajectory of the centroid of SITS is shown in Figure 15b The proposed algorithm (method 1) in the previous section is applied to optimize the energy consumption of the coupling joints in the case of continuous singularity of SITMSR. The friction torque during the motion of the coupled joint is considered, and it is assumed that it satisfies a linear relationship with the angular velocity, and the friction coefficients of q4 and q6 are C4 and C6, respectively. The convergence curve of FD with iteration is shown in Figure 16a. The average running time of motion planning based on energy saving using method 1 is 36.67 s. After the solution is completed, the trajectory of each coupled joint can be obtained, as shown in Figure 16b. It can be seen from the optimization results that when the iteration runs to Step 68, the energy consumption converges to 2.49 J, and the maximum residual error is 1.02 × 10 −12°. The rotation angle of q4 is 76.89° more than that of q6. In order to verify the effectiveness of the above solving of the algorithm, a directional path search algorithm based on the linear decoupling method is used to solve the above issue. The relationship between energy consumption and its ratio to method 1 and λ1 is shown in Figure 17a. The angular displacements of the coupled joints and energy consumption when λ1 = 1.17 are shown in Figure 17b. In Figure 17, the energy consumption of the coupled joints decreases with the increase of λ1 until it reaches the minimum value of 3.47 J. Compared with method 1, the minimum relative increase energy consumption is 28.24%, and the value of λ1 reaches the maximum at this time. If λ1 continues to increase, it will cause the joint angular displacement to exceed the range of motion.
The genetic algorithm is adopted to solve the minimum energy consumption under the condition in Example 2. Set the initial population number of the genetic algorithm to 300, and the maximum generation to 800. Since the results of the genetic algorithm are different each time, the joint angular displacements take a set of solutions corresponding to the minimum energy consumption value. The convergence curve of energy consumption of SITMSR in the case of continuous singularity using genetic algorithm is shown in Figure 18a. The angular displacements of the coupled joint corresponding to the minimum energy consumption obtained by the genetic algorithm is shown in Figure 18b. It can be seen from the Figure 18 that the minimum energy consumption of the genetic algorithm solution is 2.78 J, which is similar to method 1, which also verifies the effectiveness of the proposed solution method. The comparison of the results of the above three methods under the condition of continuous singularity configuration of SITMSR are summarized in Table 3. Table 3. The comparison of the results of the above three methods under the condition of continuous singularity configuration of SITMSR.

Method
Energy In this example, the given centroid trajectory of SITS is shown in Figure 18b, and its precession and nutation motion settings are the same as in Example 2, but the spin angular velocity is 0.5 times that in Example 2. The optimization method of continuous singularity of SITMSR is used to optimize the coupled joint motion, and then the directional path map of energy consumption is constructed for path optimization. Here, the non-singularity path πNs is compared with the path πOpt, which is also the energy-efficient path. The relationship between path length ratio and energy ratio of πNs and πOpt with time are shown in Figure 19. Although the lengths of the two paths are approximately the same, the optimal path πOpt saves 6982.36 J compared to the path πNs. The comparison results of path πNs and path πOpt in Example 3 are shown in Table 4.

Conclusions
In this paper, a motion planning method of SITMSR based on energy saving and its corresponding joint trajectory solving algorithm are proposed. Simulations based on different given motions of SIT in the ground experiment are conducted to verify the effectiveness of presented approach. The above research mainly consists of three parts: first, the motion trajectory and accuracy requirements of SIT are mapped to the ground test system (SITMSRS); second, discrete optimization mathematical model and constraints of energy consumption of SITMSR are established; third, the Dijkstra algorithm and adding the global optimization algorithm are used to get the solutions of the SITMSR path with minimal energy consumption under the condition of discontinuous singularity configuration and continuous singularity configuration, and the trajectory of each joint is obtained through polynomial fitting and discretization of solved path.
The general form of SIT makes it difficult for ground simulation robots to avoid singular positions completely, which brings new challenges to motion planning based on energy optimization. The proposed method is based on the kinematics and dynamics of the robot, which helps to solve the problem of optimal energy consumption when manipulator exists singularity configuration under a given end trajectory. From the results in this paper, the robot consumes less energy in singularity configuration under the premise that the trajectory of the motion is satisfied. In addition, the motion planning method based on energy saving is also applicable to 7-DOF manipulators and robot systems with more DOF. The algorithm presented can obtain the global optimal solution. However, the dynamics of the robot becomes complicated with the increase of DOF, which will increase the computation amount of the proposed method and the cost of hardware implementation.
Future work will further improve the efficiency of the algorithm operation. The main work is speeding up the generation process of complete path solution space of joints power based on robot dynamics. The interface of cyber-human interaction will need to be designed and the motion planning method will need to be embedded in the real-time software of the ground semi-physical test system. High-performance computer hardware will be used to support the operation of the algorithm. In addition, joint stiffness should be considered when the proposed method is applied to a cooperative manipulator with 7-DOF or robots with more DOF. The evaluation will be carried out to quantify its impact on energy consumption in the process of robot motion planning.