Next Article in Journal
On Neural Observer in Dynamic Sliding Mode Control of Permanent Magnet Synchronous Wind Generator
Next Article in Special Issue
Robust Leader–Follower Formation Control Using Neural Adaptive Prescribed Performance Strategies
Previous Article in Journal
A Novel Coupled Memristive Izhikevich Neuron Model and Its Complex Dynamics
Previous Article in Special Issue
L1 Adaptive Control for Marine Structures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Collision-Free Trajectory Planning Optimization Algorithms for Two-Arm Cascade Combination System

1
Institute of Advanced Manufacturing and Intelligent Technology, Beijing University of Technology, Beijing 100124, China
2
Genertec Machine Tool Engineering Research Institute Co., Ltd., Beijing 100102, China
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(14), 2245; https://doi.org/10.3390/math12142245
Submission received: 27 June 2024 / Revised: 11 July 2024 / Accepted: 15 July 2024 / Published: 18 July 2024

Abstract

:
As a kind of space robot, the two-arm cascade combination system (TACCS) has been applied to perform auxiliary operations at different locations outside space cabins. The motion coupling relation of two arms and complex surrounding obstacles make the collision-free trajectory planning optimization of TACCS more difficult, which has become an urgent problem to be solved. For the above problem, this paper proposed collision-free and time–energy–minimum trajectory planning optimization algorithms, considering the motion coupling of two arms. In this method, the screw-based inverse kinematics (IK) model of TACCS is established to provide the basis for the motion planning in joint space by decoupling the whole IK problem into two IK sub-problems of two arms; the minimum distance calculation model is established based on the hybrid geometric enveloping way and basic distance functions, which can provide the efficient and accurate data basis for the obstacle-avoidance constraint condition of the trajectory optimization. Moreover, the single and bi-layer optimization algorithms are presented by taking motion time and energy consumption as objectives and considering obstacle-avoidance and kinematics constraints. Finally, through example cases, the results indicate that the bi-layer optimization has higher convergence efficiency under the premise of ensuring the optimization effect by separating variables and constraint terms. This work can provide theoretical and methodological support for the efficient and intelligent applications of TACCS in the space arena.

1. Introduction

1.1. Background

In the field of modern aerospace, space robots have been used to help people complete certain tasks inside and outside space capsules. With the development of aerospace technology, different kinds of space robots have been studied and applied for on-orbit servicing tasks, such as a 6-DOFs manipulator [1], 7-DOFs manipulator [2], 14-DOFs super redundant manipulator [3], and a dual-arm humanoid robot [4]. The safe and efficient operation of space robots depends on the motion planning, which has always been the focus of researchers in recent studies [5,6]. For example, Shrivastava et al. proposed a jerk-optimized motion planning method for the redundant space robot which used the gray wolf optimization algorithm to search for the optimal result [7]; Xie et al. presented the trajectory planning method for dual-arm free-floating space robot to minimize the final base attitude change and ensure the trajectory smoothness by using an enhanced bidirectional approach [8]. In recent years, the two-arm cascade combination system (TACCS) shown in Figure 1 has been developed and applied to handling and maintenance outside capsules, and it is obtained through the docking combination of two manipulators (which can largely increase the workable area and the robotic dexterity) and better adapted to the operational tasks on the outer surface of the cylindrical cabin. Two manipulators in TACCS can be regarded as the fixed manipulator (first-level arm) and the mobile manipulator (second-level arm), respectively, and the base coordinate frame (BCF) of the latter can be changed by the joint movement of the former.

1.2. Literature Review of Motion Planning

In existing studies, the motion planning of manipulators has been the focus of research aimed at finding efficient, accurate and safe modes of application. Actually, the motion planning of TACCS has the same problem with that of the mobile manipulator (composed of a mobile car and a manipulator), which is how to consider the influence of the motion of the first-level arm or the mobile car to the motion performance of the second-level arm or the manipulator. There exist two solving methods for the motion planning of the mobile manipulator: one of them is that the motion planning algorithms of two subsystems are carried separately (called SMP in the following), the other is that the motion of the high-DOFs system is planned as a whole (called WMP in the following).
For SMP, researchers mainly focused on the improvement of the whole planning efficiency of the robotic system in previous studies. For example, Li et al. presented a hierarchical motion planning method for the mobile manipulator, which includes two planning stages: firstly the path of the two-dimensional mobile base was optimized using a hybrid sampling strategy that combined a bridge test and uniform sampling, aiming to improve the planning efficiency; secondly the secure configuration and trajectory of the manipulator was searched [9]. Chen et al. also proposed a hierarchical motion planning method which used the optimization-based A* algorithm and the sampling-based heuristic algorithm for the collision-free motion plannings of the mobile base and the manipulator, respectively, which can effectively improve the planning efficiency and success rate [10]. Rastegarpanah et al. used an A* algorithm and rapidly exploring random trees (RRTs) to separately plan the paths for the mobile platform and manipulator and completed the task of picking up an object from one table to one box [11]. In order to consider the coupling relation of configurations of two subsystems, criteria like reachability or manipulability of the manipulator’s end effector has been considered to decouple the planning of the two subsystems. For example, Colucci et al. developed a motion planning algorithm based on a manipulability index which considers the mobile base and manipulator separately and provides the closed inverse kinematics algorithm for the redundant manipulator [12,13]. Actually, the motion planning for only the mobile car or manipulator has been popular for quite a long time. Some mature technologies have been proposed. For the motion safety, some distance calculation or collision detection algorithms between the robot and surrounding obstacles based on different envelope ways were developed by scholars. For example, basic geometries like spheres and cylinders [14,15], and hierarchical bounding boxes [16] were used to envelope the robotic arms or obstacles and then applied to solve the collision-free motion planning problem of the manipulator. For the motion stability or trajectory-tracking accuracy, a cubic/quintic polynomial curve [17,18], B-spline curve [19,20] or segmented function curve [21,22] were utilized to interpolate the motion trajectory, which can eliminate vibrations caused by speed or acceleration mutations when the robot starts and stops. Moreover, the motion time [23,24], energy consumption [25], angular jerk [26,27] or their combination [28] were always taken as the objectives of the motion planning optimization of the manipulator.
For WMP, the grid-based [29], sampling-based [30] and optimization-based [31] methods were proposed in previous studies. The grid-based method are generally proposed based on the free space represented by a swept volume or a graph-creating algorithm, and since the dimension of space related to the number of robotic joints has a major impact on the planning efficiency, it is more applicable for the path planning of the mobile car or the low-DOFs robot. For example, Liu et al. considered the mobile manipulator (composed of a 2-DOFs manipulator and 3-DOFs mobile device) as a whole system and presented a path planning method based on representation space [29]. The sampling-based methods have been developed well to solve the path planning of mobile manipulators with high DOFs, including RRTs, probabilistic roadmap (PRM) and their modifications. For example, Dai et al. presented a path planning method based on a novel potential bidirectional RRT* for redundant manipulators in joint space [32] which can effectively improve the obstacle-avoidance ability of the robot; Chen et al. proposed an improved PRM based on a new sampling strategy using the virtual force field and the design of a new connection strategy to improve the planning efficiency [33]. However, the above methods are mostly used for the path planning of the robot, which refers to determining a pure geometric description of a series of positions of the end effector or mobile base; additionally, it does not consider the system kinematics, including the change in velocities or accelerations over time, which means the motion stability or motion time cannot be planned or optimized. The optimization-based methods are mostly presented by the modeling and solving of the motion planning optimization problem, wherein constraints for the obstacle avoidance, accuracy, kinematics and dynamics can be considered, and multiple optimization objectives, including the motion time, energy consumption and joint jerks, can be taken to ensure the working efficiency, energy conservation and motion stability [34,35,36]. In existing works, many intelligent optimization algorithms (like the particle swarm optimization (PSO) algorithm [35], genetic algorithm [36], etc.) have been introduced to solve the above optimization problem. But for the mobile manipulators with high DOFs, the high optimization performance (including the ability of searching for the global optimal result, the convergence efficiency) is difficult to obtain. Therefore, many scholars focus on the modification of the current optimization algorithms or the proposal of the novel optimization algorithm. For example, Jin et al. used the chaotic PSO to solve the motion trajectory planning optimization problem, which improved the premature phenomenon of the traditional PSO algorithm [37]; Li et al. proposed a novel hybrid heuristic algorithm, which combined PSO and a whale optimization algorithm (PSO-WOA) to solve the multi-objective optimization problem of the trajectory planning of the space robot, considering the end-effector pose, base disturbance, motion time and manipulability as objectives [38]; Cao et al. proposed a modified multi-objective PSO algorithm (GMOPSO) by combining a mutation operator, annealing factor and feedback mechanism for the trajectory planning of fruit picking manipulator [39].
The comparison of the existing representative works is shown in Table 1. Since only the path planning was implemented without optimization in many related works, the robot type, DOFs, planning methods, operational performance (like safety, manipulability, stability, efficiency, energy consumption) and the consideration of motion coupling of all works are listed in the comparison table in order to compare them more comprehensively. From the comparison, we can observe that the optimization-based planning method can directly obtain an optimal motion trajectory with the minimum of efficiency or energy consumption as the objective for the robot, not only a collision-free operational path, and it can better control the motion performance of the robot in terms of obstacle avoidance, motion stability, tracking accuracy, efficiency, or energy consumption, etc. However, the low planning optimization efficiency has always been a problem and the focus of the current research.

1.3. Problem Formulation

As shown in Figure 1, when one manipulator (DOFs ≥ 6) is mounted as the real end effector of the other manipulator (DOFs ≥ 6), the whole robot will be a super redundant system with DOFs of more than 12 which has been applied to perform on-orbit servicing tasks in the aerospace field, such as assembly, maintenance, the capture of cooperative and non-cooperative objects, deployment and debris removal [40]. Moreover, there exist many protruding structures outside the cylindrical cabin surface used for measurement, external docking, entrance, etc., which means the robotic motion planning needs to consider the obstacle avoidance. Therefore, the high redundancy, the large cylindrical working surface and the complex surrounding obstacles make the collision-free motion planning of TACCS quite difficult. In addition, from the summary of the current motion planning methods, the optimization-based planning method can better control the motion performance of the robot. However, since TACCS has the higher DOFs, as well as faces the special and complex working environment, the motion planning optimization becomes more difficult. The difficulties are mainly reflected in several aspects: (1) the accurate analytical solution of the inverse kinematics of the TACCS is the main basis for its motion planning in joint space, which can better ensure the motion stability of robotic joints, but it has always been a problem due to the high redundancy; (2) the accurate minimum distance between manipulator and its workspace is difficult to calculate quickly, which is very important in obtaining an efficient and reliable motion planning algorithm; (3) the number of variables and constraints increases over the robotic DOFs, which makes the trajectory optimization less efficient and more easily trapped into the local optimum.
In this work, a collision-free and time–energy–minimum (CF-TEM) trajectory planning optimization method is proposed for TACCS considering the dexterity of the second-level arm under the complex and specific working environment outside the capsule. Three problems are solved based on the following contributions: (1) the IK problem of TACCS is decoupled into two IK sub-problems of two arms and then the screw-based IK model is established, wherein the IK solution of second-level arm is related to the end-effector pose of the first-level arm; (2) the minimum distance calculation model is established based on the hybrid geometric enveloping way and basic distance functions, wherein TACCS is enveloped by a convex polyhedrons, slices-based circular/polygonal edges and planes, and the the surrounding obstacles are enveloped by convex polygonal planes; (3) the single and bi-layer optimization algorithms are both proposed by taking motion time and energy consumption as objectives and considering constraints of obstacle avoidance and robotic kinematics, and they are finally compared in terms of efficiency and effect of optimization, aiming to find an efficient and reliable optimization algorithm for the motion planning problem of the system.
The remainder of this work is organized as the following three sections: Section 2 presents the screw-based inverse kinematics of TACCS, which gives the explicit mathematical mapping from the target pose of the end effector in a base coordinate frame (BCF) of the system to joint angles of two arms. Section 3 gives the minimum distance calculation model based on the hybrid geometric envelope way, considering the characteristics of operation task and obstacles outside the capsule; Section 4 gives the single and bi-layer optimization models and algorithms of the system; Section 5 discusses the results obtained by two optimization algorithms from the efficiency and effect of optimization through a simulation case; Section 6 gives the main conclusions and future works.

2. Screw-Based Inverse Kinematics Model of TACCS

For the inverse kinematics (IK) analysis, the kinematics parameters and three coordinate frames of TACCS are shown in Figure 2, wherein the frame 1 ( o x y z ) refers to the BCF of the whole system and also the BCF of the first-level arm; the frame 2 ( o E 1 x E 1 y E 1 z E 1 ) refers to the end-effector coordinate frame (EECF) of the first-level arm and also the BCF of the second-level arm; the frame 3 ( o E 2 x E 2 y E 2 z E 2 ) refers to the EECF of the second-level arm, which is directly related to the operational task. Then, the IK solution problem can be described as the solving of joint angles of two arms by taking the end-effector target pose of the second-level arm in frame 1 as a known condition. Moreover, in Figure 2, q i represents one point taken from the rotational axis of the i t h joint of one arm of TACCS, while q e is the real end effector of each arm and l i means the kinematics parameters that will be used in the following IK model. Since the pose g s t E 1 of the end effector of the first-level arm has a major impact on the operation effect of the second-level arm, it will be optimized in the motion planning optimization of the system, which means the IK problem of TACCS can be decoupling into the IK sub-problems of two arms. Through the above decoupling processing, the IK sub-problem 1 aims to solve the joint angles of the first-level arm with the given pose g s t E 1 in frame 1; IK sub-problem 2 aims to solve the joint angles of the second-level arm with the given pose g s t E 2 in frame 1. In the traditional IK problem, the end-effector pose of manipulator is given in its own base coordinate frame, that is, for IK sub-problem 1, the inverse solution can be modeled directly using the screw-based method, but for IK sub-problem 2, the pose g s t E 2 in frame 1 must be transferred in frame 2, as then the solution can be modeled based on the new end-effector pose g s t E 2 , N .
It is noted that, for the whole IK problem of TACCS, the operation task determines the pose g s t E 2 4 × 4 , which can be expressed by a homogeneous coordinate matrix as follows:
g s t E 2 = [ R e q e 0 1 ] = [ n x o x a x x n y o y a y y n z o z a z z 0 0 0 1 ]
wherein R e 3 × 3 and q e 3 denote the orientation and the 3D coordinates of the end effector, respectively.
Moreover, in the following optimization of the pose g s t E 1 , three linear displacements Δ x , Δ y , Δ z along three axes and three rotational angles α , β , γ around three axes of frame 1 are taken as optimization variables. Then, the optimal pose g s t E 1 4 × 4 can be expressed by the following matrix:
g s t E 1 = [ c α c β c α s β s γ s α c γ c α s β c γ + s α s γ Δ x s α c β s α s β s γ + c α c γ s α s β c γ c α s γ Δ y s β c β s γ c β c γ Δ z 0 0 0 1 ]
wherein c α = c o s α ,   s α = s i n α ; c β = c o s β ,   s β = s i n β ; c γ = c o s γ ,   s γ = s i n γ .
Since the poses g s t E 2 and g s t E 2 , N have the transfer relation g s t E 2 = g s t E 1 · g s t E 2 , N , the pose g s t E 2 , N 4 × 4 can be expressed by
g s t E 2 , N = [ g s t E 1 ] 1 · g s t E 2
Based on the Chasles principle, the motion screw of the rotational joint is defined by ξ = [ ω θ ] T , and the exponential form of the screw coordinate can be written as follows,
e x p ( ξ ^ θ ) = [ e x p ( ω ^ θ ) [ I e x p ( ω ^ θ ) ] ( ω × v ) + θ ω ω T v 0 1 ]
wherein ω is the unit vector of the rotational axis; ω ^ is the anti-symmetric matrix of ω , while setting ω = ( ω 1 , ω 2 , ω 3 ) , ω ^ = [ 0 ω 3 ω 2 ω 3 0 ω 1 ω 2 ω 1 0 ] ; θ represents the rotational angle; v means the linear velocity of the rotational motion, v = q × ω , where q denotes the 3D coordinate of one point located on the rotational axis; I is a 3 × 3 identity matrix; e x p ( ω ^ θ ) can be calculated by e x p ( ω ^ θ ) = I + ω ^ s i n θ + ω ^ 2 ( 1 c o s θ ) .
Based on the above analysis, two IK sub-problems can be transferred into the general IK problems. Taking the robotic model shown in Figure 2 as the example case, two arms are manipulators with the same configuration but different values to the kinematics parameters listed in Table 2. According to our previous work [41], the IK problem of each arm can be firstly decomposed into three sub-problems: one is for the rotation about three non-intersecting axes; the second is for the rotation about two intersecting axes; and the third is for the rotation about one single axis. For the first two IK sub-problems, the rotational motions are shown as Figure 3, respectively, wherein c 1 , c 2 , c 3 are the three transition points of two motions, and their 3D coordinates c 1 = [ x 1 , y 1 , z 1 ] , c 2 = [ x 2 , y 2 , z 2 ] , c 3 = [ x 3 , y 3 , z 3 ] can be obtained by solving a hexagonal quadratic equation system and a ternary quadratic equation system. Based on the above calculations, the IK solution of the arm (the angular displacements of all joints of the arm) can be explicitly modeled and can be expressed as Equation (5). The detailed solving process can be referred to through reference [41].
{ θ 1 = atan 2 ( x 2 y y 2 x , x 2 x y 2 y ) θ 2 = atan 2 [ x 2 ( z 1 l 1 ) x 1 ( z 2 l 1 ) , x 1 x 2 ( z 1 l 1 ) ( z 2 l 1 ) ] θ 3 = atan 2 [ x 1 l 4 , l 4 ( z 1 l 1 l 3 ) ] θ 4 = atan2 ( x 3 y e 1 y 3 x e 1 , x 3 x e 1 y 3 y e 1 ) θ 5 = atan2 [ x 3 l 2 , l 2 ( z 3 l 1 l 3 l 4 ) ] θ 6 = atan 2 ( y e 2 , x e 2 )
where, [ x e 1 , y e 1 , z e 1 , 1 ] T = [ e x p ( ξ ^ 1 θ 1 ) e x p ( ξ ^ 2 θ 2 ) e x p ( ξ ^ 3 θ 3 ) ] 1 [ x s 1 , y s 1 , z s 1 , 1 ] T . In this equation, [ x s 1 , y s 1 , z s 1 ] is the initial 3D coordinates of one point located on the sixth joint axis but not on the fourth and fifth axes; [ x e 2 , y e 2 , z e 2 , 1 ] T = [ e x p ( ξ ^ 1 θ 1 ) e x p ( ξ ^ 5 θ 5 ) ] 1 [ x s 2 , y s 2 , z s 2 , 1 ] T , in this equation, [ x s 2 , y s 2 , z s 2 ] is the initial 3D coordinates of one point located out of the sixth joint axis.
Then, IK models can be written as two functions, θ A 1 = f u n _ I K _ A 1 ( g s t E 1 ) and θ A 2 = f u n _ I K _ A 2 ( g s t E 2 , N ) , which together form an important basis for the motion planning optimization of the TACCS.

3. Minimum Distance Calculation Based on Hybrid Geometric Envelope

The minimum distance calculation is an important basis for the spatial constraint determination of the motion planning optimization. In real applications, the TACCS will be applied to perform works with targets at different places on the capsule body. Moreover, the outer surface of the capsule may have some raised structures for docking or measuring, which makes the accurate distance calculation more difficult. In order to solve the above problem, this work uses the circular slices and convex polygonal planes to envelop the robotic links and connection structure between two arms, respectively, also using the convex polygonal planes to envelop the capsule and other obstacle structures located at the surface of the capsule, respectively. The envelope structures are shown in Figure 4. From Figure 4, the environmental obstacles include convex obstacles and cylindrical obstacles. Since the distance calculation between the circular slice and the cylinder is difficult, the cylinder is also enveloped by convex polygonal planes, which needs to firstly envelope the end circle using its outer joined regular polygon. For the envelope of a circle, the largest error between the original and new geometry can be calculated by ϵ = R [ 1 / c o s ( π / N e ) 1 ] (as shown in Figure 5); therefore; if we hope ϵ 5   mm , the number of edges must meet the condition as N e π / a c o s [ R / ( 5 + R ) ] (based on this condition, when R = 50   mm , N e 7.3 ).
In our previous study [42], the basic functions of the minimum distance between the circular edge and the convex polygonal plane and that between two convex polygonal planes have been established, which can be represented by f u n _ C P ( ρ , ς ) and f u n _ P P ( ς 1 , ς 2 ) , wherein ρ represents one circular edge and is determined by the coordinate of the circular center and ς represents one convex polygonal plane and is determined by the coordinates of polygonal vertices. Based on the envelope ways shown in Figure 3, the structure of the TACCS is finally expressed as a set of geometries, including circular edges and convex polygonal planes in mathematics, and the structure of its workspace obstacles are expressed as a set of convex polygonal planes. Therefore, the minimum distance between TACCS and its workspace obstacles can be written as follows:
D I = min { f u n _ C P ( ρ R , i , τ O , k ) ,   f u n _ P P ( ρ R , j , τ O , k )   1 i I ,   1 j J ,     1 k K }
wherein ρ R , i and ρ R , j denote the i t h circular edge and the j t h polygonal plane; τ O , k means the k t h polygonal plane of obstacles; I , J are the numbers of circular edges and polygonal planes of the TACCS, and K is the number of polygonal planes of obstacles.
As shown in Figure 5, the envelope of circle will lead to the largest error, ϵ . In order to further analyze the influence of ϵ on the error of the shortest distance calculation, the distance between the circular plane with different radii ( R = 50 ,   100 ,   150 ,   200   mm ) and an circular edge ( R e = 150   mm ) are taken as objectives. For the analysis, the center locations of two circles are set by C = ( 0 , 0 , 0 ) , C e = ( 200 , 300 , 400 ) , and their normal vectors are set by n = ( 0 , 0 , 1 ) , n e = 3 3 · ( 1 , 1 , 1 ) ,   respectively . During the calculation process, taking different values for ϵ , and determining the minimum number of edges by N e = ceil { π / a c o s [ R / ( ϵ + R ) ] } ( ceil { } represents the function that returns the minimum integer larger than the given real number), the shortest distance between the circular plane and the circular edge can be predicted by the distance between the convex polygonal plane ς and the circular edge ρ can be calculated by the function d p = f u n _ C P ( ρ , ς ) ; the real shortest distance d r can be calculated by taking N e = 1000 , which can make the convex polygonal plane close enough to the circular outline; then the error of the shortest distance calculation is determined by E = d r d p .
Based on the above strategy, the errors of the shortest distance under different values of ϵ are shown in Figure 6, from which the following conclusions can be obtained: (1) The error basically increases with the increase of the value of ϵ ; this is because the increase of ϵ leads to the lower number of edges of the enveloped polygon (for example, when ϵ 5 , the errors for four circular planes are all less than 5 mm, but when R = 15 mm, ϵ = 17 , the error reaches 11.51 mm). (2) In some cases, even where ϵ takes the large value, the error may be very small. The reasons for that are as follows: ϵ represents the largest error between the original circle and its enveloped polygon, which means some points on their edges coincide or are close together; moreover, the shortest distance between the enveloped polygonal plane and the edge must be the distance between two point elements located on two geometries, and so when the point element at the polygonal plane is located just at the position where the envelope error is small, the predicted shortest distance will be very close to its real value. (3) In order to both ensure the accuracy and efficiency, the reasonable value for ϵ should be taken, since, the smaller the value is, the lower the calculation efficiency and the higher the calculation accuracy. In real applications, the value can be taken according to the safety threshold value, which is set for the absolute safety of the robotic motion.
Moreover, the pose of the TACCS determines the locations of feature points of the envelope geometries, which can be calculated based on screw-based conversion matrices as follows:
q e = [ i = 1 N e x p ( ξ ^ i θ i ) ] q s
wherein q s and q e represents the start and final homogeneous coordinates of one point in frame 1, and they are 4D vectors extended from the corresponding 3D coordinate q s by adding a unit column, i.e., q s = ( q s , 1 ) T 4 ; N means the number of joints that affects the pose of this point.
Therefore, since the angular displacement θ i of the i t h joint changes over time during the operational motion, locations of feature points are related to time, which means the minimum distance during motion can be described as a function over time and expressed by D I ( t ) . The above function is a major basis for the following motion planning optimization.

4. Single and Bi-Layer Optimization Models and Algorithms

In the motion planning problem of TACCS, the variables of optimization include the following six motion variables: Δ x , Δ y , Δ z , α , β , γ , T A 1 and T A 2 . T A 1 and T A 2 denote the motion times of the first-level arm and the second-level arm, respectively. In order to ensure the operational stability, it is supposed that two arms do not move simultaneously in real applications, which means although the motion dexterity of the second-level arm will be influenced by the pose of the first-level arm, the coupling influence is not dynamic. For this problem, considering both the total motion time and the motion energy consumption as the optimization objective, the objective of the single optimization can be modeled as follows:
f S = w 1 ( T A 1 + T A 2 ) + w 2 w { i = 1 D A 1 t = 0 T A 1 S A 1 , i ( t ) + i = 1 D A 2 t = 0 T A 2 S A 2 , i ( t ) }
wherein w 1 and w 2 are weighting factors for items of the motion and energy consumption, w 1 + w 2 = 1 ; w denotes the coefficient used to eliminate the influence of the magnitude difference of items; D A 1 and D A 2 mean the numbers of DOFs of two arms; S A 1 , i ( t ) and S A 2 , i ( t ) are time functions of the angular displacement of the i t h joint of two arms, respectively, which can be obtained through motion trajectory interpolation. In this work, the quintic polynomial interpolation method is used, and the function of one joint angle can be expressed by
S ( t ) = a 1 t 5 + a 2 t 4 + a 3 t 3 + a 4 t 2 + a 5 t + a 6
wherein the coefficients of the function can be calculated by a 1 = 6 ( θ E θ S ) / T 5 , a 2 = 15 ( θ E θ S ) / T 4 , a 3 = 10 ( θ E θ S ) / T 3 , a 4 = a 5 = 0 and a 6 = θ S , respectively, which are obtained by the boundary conditions according to requirements of movement stability as S ( 0 ) = θ S , S ( T ) = θ E and S ˙ ( 0 ) = S ˙ ( T ) = S ¨ ( 0 ) = S ¨ ( T ) = 0 ; S ˙ ( t ) and S ¨ ( t ) are the functions of angular velocity and acceleration of the joint and can be obtained through one derivative and two derivatives of S ( t ) ,   respectively .
Then, considering kinematics and obstacle-avoidance constraints, the single optimization problem can be modeled as follows:
{ minimize : f = f S ( Δ x , Δ y , Δ z , α , β , γ , T A 1 , T A 2 ) subject to : Δ x [ Δ x m i n , Δ x m a x ] Δ y [ Δ y m i n , Δ y m a x ] Δ z [ Δ z m i n , Δ z m a x ] α [ α m i n , α m a x ] β [ β m i n , β m a x ] γ [ γ m i n , γ m a x ] S A 1 , i ( t ) D A 1 , i ; S A 2 , i ( t ) D A 2 , i max [ | S ˙ A 1 , i ( t ) | ] V A 1 , i , m a x ; max [ | S ¨ A 1 , i ( t ) | ] A A 1 , i , m a x max [ | S ˙ A 2 , i ( t ) | ] V A 2 , i , m a x ; max [ | S ¨ A 2 , i ( t ) | ] A A 2 , i , m a x min [ D I ( t ) ] D I thred T A 1 + T A 2 T m a x
where in the subscripts m i n and m a x mean the upper and lower bounds of the parameter, for example Δ x m i n and Δ x m a x , are the upper and lower bounds of the linear displacement of the end effector of the first-level arm, respectively; D A 1 , i and D A 2 , i are the allowable range of the angular displacement of the i t h joints of two arms, respectively; the functions max [ g ] and min [ g ] refer to taking the maximum value and the minimum value of the function g respectively; V A 1 , i , m a x and A A 1 , i , m a x are the maximum values of the angular velocity and acceleration allowed by the i t h joints of the first-level arm, respectively; D I ( t ) is the minimum distance function of TACCS; and D I t h r e d denotes the threshold value taken for the reliable obstacle avoidance.
Therefore, even if the allowable ranges or maximum values of the angular displacements, velocity and acceleration of all joints of each arm are the same, there are 14 constraint conditions in the above model. In order to increase the optimization efficiency, the upper and lower values of variables Δ x , Δ y , Δ z , α , β and γ need to be limited as much as possible. For this purpose, this work proposed a feasible space analysis method based on the IK functions of two arms. In this method, the Monte Carlo method is firstly used to randomly take values for variables Δ x , Δ y , Δ z , α , β , γ ; then the target pose of the end effector of the first-level arm and the base frame of the second-level arm can be determined, based on which the IK functions of two arms can be solved; finally, values of variables with valid solutions in two IK problems can be recorded, based on which the feasible space of frame 2 (as Figure 2) is obtained, and can help to determine the upper and lower bounds of the corresponding variable. The framework of this method is shown in Figure 7.
Since the number of variables and constraint conditions have a major impact on the convergence time and effect of the optimization, the bi-layer optimization strategy is also proposed in this work. In this strategy, the optimization problem is decoupled into two layers of sub-optimizations: the first layer of sub-optimization is the end-effector pose optimization of the first-level arm considering the kinematics and obstacle-avoidance constraints, as well as its influence on the operational dexterity of the second-level arm; the second layer includes two sub-optimizations, which are the motion trajectory optimizations of two arms, respectively, based on the optimal result of the first layer. Moreover, the bi-layer optimization problem also considers the motion time and energy consumption of the arms. Therefore, the first sub-optimization takes the sum of Euclidean distances of all joints of two arms as the object, which refers to the equivalent index of the energy consumption; the other sub-optimizations take the motion time and the length of motion path of the arm as the object. Then, the objectives of three sub-optimizations can be expressed by
f B 1 = i = 1 D A 1 θ A 1 , i 2 + i = 1 D A 2 θ A 2 , i 2
f B 2 = w 1 T A 1 + w 2 w i = 1 D A 1 t = 0 T A 1 S A 1 , i ( t )
f B 3 = w 1 T A 2 + w 2 w i = 1 D A 2 t = 0 T A 2 S A 2 , i ( t )
wherein θ A 1 , i and θ A 2 , i are the inverse solutions of the i t h joint at the target pose of the end effector for two arms, respectively.
The bi-layer optimization framework is shown in Figure 8, wherein there exists a global judgement for obstacle avoidance of the whole robotic motion before the second-layer sub-optimizations. From this framework, the optimal values of the first-layer optimization is the input for sub-optimization of the second-layer optimization, which means two sub-optimizations can perform parallel computing and output optimal trajectories. In actual applications, the parallel computing can largely improve the execution efficiency of the bi-layer optimization algorithm.
Moreover, in the global judgement, the motion times of two arms are firstly taken reasonably to determine the overall robotic poses at different motion times, since their values have a minimal influence on the robotic pose at a certain point of the trajectory but largely affect the motion speed and acceleration of each joint; then, the smallest distance between the TACCS and its surrounding obstacles during motion can be calculated based on the model of the minimum distance of Section 3; finally, the smallest distance is compared with its threshold value (taken to ensure the absolutely safety of the system) to complete the judgement, which can detect the possible collisions that may happen during motion accurately and reliably. In the bi-layer optimization framework, if the global judgement condition is not satisfied, it then returns to the first-layer optimization, or else it continues the second-layer optimizations. The above strategy can further help improve the efficiency of the whole optimization, which can be modeled as follows:
{ minimize : f = f B 1 ( Δ x , Δ y , Δ z , α , β , γ ) subject to: : Δ x [ Δ x min , Δ x max ] Δ y [ Δ y min , Δ y max ] Δ z [ Δ z min , Δ z max ] α [ α min , α max ] β [ β min , β max ] γ [ γ min , γ max ] θ A 1 , i D A 1 , i ; θ A 2 , i D A 2 , i D I E D I thred Judgement : min [ D I ( t ) ] D I thred minimize : f = f B 2 ( T A 1 ) subject to : S A 1 , i ( t ) D A 1 , i max [ | S ˙ A 1 , i ( t ) | ] V A 1 , i , max max [ | S ¨ A 1 , i ( t ) | ] A A 1 , i , max T A 1 T A 1 , max minimize : f = f B 3 ( T A 2 ) subject to : S A 2 , i ( t ) D A 2 , i max [ | S ˙ A 2 , i ( t ) | ] V A 1 , i , max max [ | S ¨ A 2 , i ( t ) | ] A A 2 , i , max T A 2 T A 2 , max
wherein D I E means the minimum distance when the TACCS is at its target pose.
In this work, the particle swarm optimization algorithm (PSO) is adopted to solve the optimization models proposed above; the PSO is an intelligent algorithm that searches for the optimal solution in the solution space by a certain number of particles and has advantages of easy implementation and fewer algorithm parameters. Based on the calculation functions of IK and the minimum distance of TACCS, pseudo codes of algorithms based on the single and bi-layer optimizations can be written as Algorithm A1 and Algorithm A2 in Appendix A. For the PSO, parameters include inertia weight ω , maximum flying velocity v m a x of particles and the number of particles N p , and they have major impacts on the convergence speed and optimization effect. Among them, the value of v m a x for each variable can be taken by 10% to 20% of its searching range, and if the ranges for the variables are different, the ratio φ of the value to its corresponding range can be taken as a uniform parameter for values of the maximum velocities of all variables, 0.1 φ 0.2 . Therefore, for the motion planning optimization, reasonable values of PSO parameters ( ω , φ and N p ) should be taken to obtain the optimal results according to several trails of each optimization. Moreover, in order to ensure obtain the global optimal result, the stopping criterion for iteration is defined by max [ | f k τ 1 f k τ | ,   0 τ 14 ] < ε , which means that the change in the objective value during fourteen successive iterations is less than a threshold value ε = 0.001 .

5. Verification Based on Simulation Cases

In order to show the performance of two proposed optimization algorithms, this section designed two example cases. The principle of determining initial conditions (including the target pose and locations of obstacles) for the simulations involves ensuring the solvability of the optimization problem, i.e., there is at least one collision-free path during the conversion from the initial pose to the target pose. The operational environment of case 2 is shown in Figure 4 (the target structure is located on one side of capsule), and for that of case 1, the target structure and its surrounding obstacles are located on the top of the capsule. The target pose and the initial 3D locations of the feature points of the TACCS and obstacles in frame 1 (Table A1) are given for the following simulation analysis based on the above basic principle.
From the proposed method, there also exists one important step before the optimization, which is to determine the upper and lower values of variables Δ x , Δ y , Δ z , α , β and γ according to the framework shown in Figure 8. Figure 9 and Figure 10 show the feasible spaces of variables in two example cases, wherein valid solutions of variables Δ x , Δ y and Δ z and that of variables α , β and γ are plotted as points in a 3D coordinate system, respectively, which finally form point clouds to show the feasible spaces.
In each sub-figure of Figure 9 and Figure 10, the red box is used to show the final space size chosen to further decrease the allowable ranges of variables that can help improve the optimization efficiency. The box selection principles include (1) trying to select regions with dense point clouds, which refers to the point clouds with higher density; (2) ensuring that the proportion of selected points to all points is greater than a reasonable value; and (3) taking the limits of the red box based on the above proportion requirement. From Figure 9 and Figure 10, since the density of the points is uniform, the regions for variables α , β and γ are taken as the whole; the distribution of point clouds for variables Δ x , Δ y and Δ z features sparse edges and a dense middle, and so a reasonable value of the proportion can largely improve the optimization efficiency under the premise of ensuring the optimization effect. In order to show the influence of the proportion of variables Δ x , Δ y and Δ z on the optimization effect, the iterative curves of the sub-optimization 1 (since the size of the allowable range of variables Δ x , Δ y and Δ z has a major impact on sub-optimization 1 of the bi-layer optimization) and the single optimization under different values of the proportion (50%, 70%, 90% and 100%) are shown in Figure 11, respectively, from which, when the proportion is taken by 90%, the optimization efficiency can be largely improved with only a small loss of the objective value. From two figures, the allowable ranges of six variables in two example cases can be finally obtained.
Moreover, as presented in Section 4, PSO parameters (including the inertia weight ω , the velocity ratio φ and the particle number N p ) have a major influence on the optimization effect, and they are always related to ranges of variables and the complexity of the optimization problem. Therefore, different values for these parameters should be taken to perform the optimizations before cases analysis, aiming to determine the reasonable values for PSO parameters and ensure the optimization effect. Iterative curves for the single optimization and three sub-optimizations of the bi-layer optimization with different values of PSO parameters are shown in Figure 12. From this figure, the reasonable values for the PSO parameters can be determined as [ ω = 0.8 ,   φ = 0.15 ,   N p = 15 ] , [ ω = 1.2 ,   φ = 0.15 ,   N p = 10 ] and [ ω = 0.8 ,   φ = 0.10 ,   N p = 10 ] for the first-layer optimization (sub-optimization 1), the second-layer sub-optimizations of the bi-layer optimization (sub-optimizations 2 and 3) and the single optimization, respectively, which can ensure the fast convergence speed and global optimization capability.
Since the randomness of PSO iterations has a major impact on the convergence time of optimization, in order to analysis the influence degrees of the randomness for two proposed algorithms, the single and bi-layer optimizations for two example cases are performed 30 times. The comparison of convergence times of two optimization algorithms is shown in Table 3, wherein μ , σ and c v mean the average convergence time, the standard deviation and the coefficient of variation of 30 optimizations, respectively, while μ is used to evaluate the average optimization efficiency, σ and c v are used to evaluate the robustness of the algorithms, σ can reflect the dispersion of convergence time, σ = 1 30 i = 1 30 ( t i μ ) 2 and c v is the dimensionless representation of σ , c v = σ / μ × 100 % , (and it can be used to compare the convergence performances of different optimizations). From the result, we can obtain the following conclusions: (1) The average convergence times of the bi-layer optimization algorithm in two cases are reduced by 72.43% and 63.64%, respectively, compared with that of the single optimization algorithm, which are calculated by ( μ s μ b ) / μ s × 1 00%, wherein μ s and μ b are the average convergence time of the single and bi-layer algorithms. The reason is that reducing the number of variables can largely increase the optimization speed; additionally, the collision judgement during the whole motion is placed outside the optimization iteration process. (2) Values of the dispersion of the bi-layer optimization algorithm in two cases are reduced by 79.52% and 82.53%, respectively, and values of the coefficient of variation in two cases are reduced by 25.68% and 51.90%, respectively, compared with that of the single optimization algorithm. That is, the randomness of PSO iterations has a greater impact on the convergence time of the single optimization algorithm. (3) Although the initial conditions have a major impact on the convergence time of optimizations, the average convergence times of the single optimization in two cases are obviously larger than that of the bi-layer optimization, which indicates that the advantage of bi-layer optimization is confidently certain under different initial conditions.
The constraints and the optimal results of the main parameters based on three algorithms, including the bi-layer optimization, single optimization and traditional optimization, are listed in Table 4, wherein V m a x and A m a x mean the maximum values of angular velocity and acceleration of all robotic joints and the values of difference 1 and difference 2 are calculated by | p s p b | / p s × 100 % and | p t p b | / p t × 100 % , while p s , p b and p t represent the parameter values based on single optimization, bi-layer optimization and traditional optimization, respectively. The results show that values for parameters Δ x , Δ y , Δ z , α , β , γ , V m a x and A m a x , based on three algorithms, are all less than their corresponding limits given in constraints; compared with the bi-layer optimization, the overall optimization objectives in two cases obtained based on the traditional optimization increase 6.94% and 4.36%, respectively. In particular, the energy consumption increases 11.21% in case 1, and the motion time increases 22.75% in case 2; the differences of the overall optimization objectives based on the bi-layer optimization and the single optimization in two cases are only 3.21% and 1.33%, respectively. The main reasons for the differences include (1) the hierarchical processing of the optimization problem taking into account the end-effector pose of the first-level arm on the motion of the second-level arm, which can rapidly find the optimal result and finally improve the optimization efficiency, but it also sacrifices a small amount of the valid solution space, which leads to some objectives based on single optimization being potentially less than that based on bi-layer optimization (for example, the overall objective based on single optimization is 3.21% smaller than that based on bi-layer optimization in case 1); (2) in spite of the above reasons, it is very difficult for single optimization to find the real global optimal solution, so its objectives may be higher than that of the bi-layer optimization (for example, the overall objective based on single optimization is 1.33% larger than that based on bi-layer optimization). For the constraints based on bi-layer optimization, the maximum velocity and acceleration of robotic joints are 77.71 d e g / s and 103.80 d e g / s 2 , and they are obviously less than their allowable values 120 d e g / s and 200 d e g / s 2 , respectively. Also, Figure 13 gives the variation in the shortest distance between the TACCS and its obstacles over motion time based on two algorithms, which indicates that the optimal trajectories of the TACCS based on two algorithms can ensure the motion safety of the system; additionally, in order to clearly show the details of the shortest distance between the TACCS and obstacles, the corresponding structure no. and obstacle no. are also shown in Figure 10 (for example, R2-O1 means the distance between the second robotic structure and the first obstacle is the shortest distance, wherein the no. information can be found in Table A1).
Moreover, taking example case 1 as the objective, curves of angular parameters of each joint of TACCS and a series of motion poses relative to capsule based on the bi-layer optimization algorithm are shown in Figure 14 and Figure 15, and those based on the single optimization algorithm are shown in Figure A1 and Figure A2, which can further intuitively demonstrate the optimization effect of two algorithms, including the guarantee of the motion stability and safety. Moreover, taking example case 2 as the objective, a series of motion poses relative to the capsule based on the bi-layer and single optimization algorithms are shown in Figure A3 and Figure A4.
From the comprehensive analysis results of convergence speed and optimization effect, it can be concluded that the bi-layer optimization has a higher execution efficiency under the premise of ensuring the optimization effect.

6. Conclusions and Future Works

6.1. Conclusions

This work proposed the single and bi-layer optimization algorithms for the collision-free motion planning problem of the TACCS, aiming to optimize the operating time and energy consumption and considering both the kinematics and obstacle-avoidance constraints under the complex operational environment. For the implementation of optimizations, the functions for the inverse kinematics of TACCS and the shortest distance between the TACCS and its obstacles are established based on screw theory and a hybrid geometric envelope, respectively. During the optimization problem modeling, a single optimization problem with more than eight variables is converted into three optimization sub-problems with fewer variables of bi-layer optimization, which aims to improve the convergence efficiency of optimizations. Finally, two algorithms are applied into two example cases. Through experimental simulations, three main conclusions are obtained as follows:
(1)
The reasonable PSO parameters are selected to perform optimizations of two algorithms, both considering their influences on convergence time and optimization effect. Moreover, a box selection principle is designed to determine the allowable ranges of motion variables of the first level arm, which can further improve the execution speed of algorithms by sacrificing very little workspace. The above strategies can help ensure and improve the optimization performances of two algorithms.
(2)
The performance indexes are taken as the average value and dispersion of convergence times of 30 executions, which are obtained in two example cases and based on two algorithms. The comparison of performance indexes indicates that the bi-layer optimization has a higher convergence speed than the single optimization, which can be increased by more than 60%, and the randomness of the PSO iterations has less of an impact on the convergence speed of the bi-layer optimization.
(3)
The optimal results based on two algorithms are compared to show the significance of this work. The results indicate that values of motion parameters and shortest distances during the whole motion are less than their corresponding limits given in constraints; additionally, the overall optimization objective based on the bi-layer optimization is smaller than that based on the single optimization. Combining with the comparison result of convergence speed, it can be concluded that the bi-layer optimization has a higher execution efficiency under the premise of ensuring the optimization effect.

6.2. Future Works

In this work, the proposed algorithms can be used for the collision-free motion planning of the TACCS under the complex structured environment, for which the locations of obstacles and targets are given information, i.e., they are fixed or regularly changed. Therefore, this work cannot solve the real-time collision-free motion planning of the TACCS in randomly dynamic scenarios, which will be one focus of future works. Moreover, the real application of the proposed algorithms will be the other focus of future works, for which we will actively seek cooperation opportunities with aerospace enterprises and apply our result in actual scenarios.

Author Contributions

Conceptualization and methodology, J.X.; validation and writing—original draft preparation, L.T. and Y.P.; writing—review and editing, all authors; supervision, Q.C. and H.C.; project administration, J.X. and T.Z.; funding acquisition, Q.C., H.C. and T.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key R&D Program of China (No. 2023YFB3406400), the Joint Funds of the National Natural Science Foundation of China (No. U23B20104), the 2022 industrial technology basic public service platform (No. 2022-232-223-02), and the R&D Program of the Beijing Municipal Education Commission (No. KM202210005031).

Data Availability Statement

The original contributions presented in the study are included in the article; further inquiries can be directed to the corresponding author.

Conflicts of Interest

Author Yanhu Pei was employed by the Genertec Machine Tool Engineering Research Institute Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The Genertec Machine Tool Engineering Research Institute Co., Ltd. had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Nomenclature

TACCSTwo-arm cascade combination system: a system that one manipulator is mounted at the end effector of the other manipulator.
IKInverse kinematics: a mapping from the end-effector pose to the angular displacements of joints for one manipulator.
DOFsDegrees of freedom.
BCFBase coordinate frame of the manipulator.
SMPMotion planning algorithms of two subsystems are carried separately for the mobile manipulator.
WMPMotion planning of the high-DOFs system is carried in a whole.
A* algorithmA-star algorithm: a heuristic search algorithm to find the shortest path from the starting point to the target point in a graph.
RRTRapidly exploring random trees: build a tree by random sampling and gradually expand the tree to approach the target point, eventually finding a viable path from the starting point to the end point.
PRMProbabilistic roadmap: the possible motion path of the robot is constructed by probabilistically sampling points in the configuration space and establishing edges between the connectable points.
PSOParticle swarm optimization.
PSO-WOAA novel hybrid heuristic algorithm, which combined PSO and whale optimization algorithm.
GMOPSOA modified multi-objective PSO algorithm.
CF-TEMA collision-free and time–energy–minimum trajectory planning optimization method.
EECFEnd-effector coordinate frame of the manipulator.
3D coordinateThree-dimensional coordinate.
4D vectorFour-dimensional vector.

Appendix A. Optimization Algorithms

Based on the calculation functions of IK and minimum distance of TACCS, pseudo codes of algorithms based on the single and bi-layer optimizations can be written as Algorithm A1 and Algorithm A2 as follows.
Algorithm A1 Pseudo Codes of the Single Optimization Algorithm
1.Input:  g s t
2.PSO parameters setting N p , ω , c 1 , c 2 , φ ;
3.Initialize:   [ Δ x ,   Δ y ,   Δ z ,   α ,   β ,   γ ,   T A 1 ,   T A 2 ]
4.while count 15 do
5.    for  n = 1   to   N p  do
6.     Calculate θ A 1 f u n _ I K A 1 ( [ Δ x ,   Δ y ,   Δ z ,   α ,   β ,   γ ] ) ,
7.     Calculate θ A 2 f u n _ I K A 2 ( g s t , [ Δ x ,   Δ y ,   Δ z ,   α ,   β ,   γ ] ) ,
8.     if θ A 1   & &   θ A 2  then
9.       Determine S A 1 , i ( t A 1 ) , S ˙ A 1 , i ( t A 1 ) , S ¨ A 1 , i ( t A 1 ) , D I ( t A 1 ) ;
10.       Determine S A 2 , i ( t A 2 ) , S ˙ A 2 , i ( t A 2 ) , S ¨ A 2 , i ( t A 2 ) , D I ( t A 2 ) ;
11.       Calculate θ A 1 , i , m i n = m i n S A 1 , i ( t A 1 ) and θ A 1 , i , m a x = m a x S A 1 , i ( t A 1 ) ,   0 t A 1 T A 1 .
12.       Calculate θ A 2 , i , m i n = m i n S A 2 , i ( t A 2 ) and θ A 2 , i , m a x = m a x S A 2 , i ( t A 2 ) ,   0 t A 2 T A 2
13.       Calculate θ ˙ A 1 , i , m i n = m i n S ˙ A 1 , i ( t A 1 ) and θ ˙ A 1 , i , m a x = m a x S ˙ A 1 , i ( t A 1 ) ,   0 t A 1 T A 1 .
14.       Calculate θ ˙ A 2 , i , m i n = m i n S ˙ A 2 , i ( t A 2 ) and θ ˙ A 2 , i , m a x = m a x S ˙ A 2 , i ( t A 2 ) ,   0 t A 2 T A 2
15.       Calculate θ ¨ A 1 , i , m i n = m i n S ¨ A 1 , i ( t A 1 ) and θ ¨ A 1 , i , m a x = m a x S ¨ A 1 , i ( t A 1 ) ,   0 t A 1 T A 1 .
16.       Calculate θ ¨ A 2 , i , m i n = m i n S ¨ A 2 , i ( t A 2 ) and θ ¨ A 2 , i , m a x = m a x S ¨ A 2 , i ( t A 2 ) ,   0 t A 2 T A 2
17.       Calculate D I A 1 , m i n = m i n D I ( t A 1 ) ,   0 t A 1 T A 1
18.       Calculate D I A 2 , m i n = m i n D I ( t A 2 ) ,   0 t A 2 T A 2
19.       if  θ A 1 / A 2 , i , m i n , θ A 1 / A 2 , i , m a x D A 1 / A 2 , i   & &   θ ˙ A 1 / A 2 , i , m i n , θ ˙ A 1 / A 2 , i , m a x D ˙ A 1 / A 2 , i   & &  
θ ¨ A 1 / A 2 , i , m i n ,   θ ¨ A 1 / A 2 , i , m a x D ¨ ˙ A 1 / A 2 , i   & &   D I A I / A 2 , m i n D I t h then
20.             Calculate f i t ,
21.                 if  f i t n < f p  then
22.                        f p , n = f i t n ; x p , n = x n .
23.                  end if
24.         end if
25.       end if
26.     end for
27.     if m i n ( f p ) < f g  then
28.        [ f g , n u m ] = m i n ( f p )   ; x g = x p , n u m .
29.    end if
30.    if f g , k 1 f g , k     0.001  then
31.        count = count + 1
32.     else
33.         count = 1
34.end if
35. Update the current flying velocities and positions of all particles;
36.end while
Algorithm A2 Pseudo Codes of the Bi-Layer Optimization Algorithm
1.Input:  g s t
2.PSO parameters setting N p , ω , c 1 , c 2 , φ ;
3.Initialize:   [ Δ x ,   Δ y ,   Δ z ,   α ,   β ,   γ ]
4.while condition 0 do
5.   while count1 15 do
6.     for  n = 1   to   N p  do
7.      Calculate θ A 1 f u n _ I K A 1 ( [ Δ x ,   Δ y ,   Δ z ,   α ,   β ,   γ ] ) ,
8.      Calculate θ A 2 f u n _ I K A 2 ( g s t , [ Δ x ,   Δ y ,   Δ z ,   α ,   β ,   γ ] ) ,
9.      if θ A 1   & &   θ A 1   & &   θ A 1 / A 2 D A 1 / A 2  then
10.       Calculate D I θ A 1 , D I θ A 2
11.       if    D I θ A 1 / θ A 2 D I t h  then
12.    Calculate f i t 1 ,
13.      if  f i t 1 n < f 1 p  then
14.         f 1 p , n = f i t 1 n ; x 1 p , n = x 1 n .
15.       end if
16.        end if
17.       end if
18.end for
19.if   m i n ( f 1 p ) < f 1 g  then
20.        [ f 1 g , n u m ] = m i n ( f 1 p )   ; x 1 g = x 1 p , n u m .
21.end if
22.if f 1 g , k 1 f 1 g , k     0.001  then
23.       count1 = count1 + 1
24.else
25.      count1 = 1
26.end if
27.  Update the current flying velocities and positions of all particles;
28. end while
29. Calculate D I A 1 , m i n = m i n D I ( t A 1 , C 1 ) ,   0 t A 1 , C 1 T A 1 , C 1
30. Calculate D I A 2 , m i n = m i n D I ( t A 2 , C 2 ) ,   0 t A 2 , C 2 T A 2 , C 2
31. if    D I A I / A 2 , m i n D I t h  then
32. condition = 1
33.  end if
34.end while
35.Initialize: T A 1
36.while count2 15 do
37.  for  n = 1   to   N p  do
38. Determine S A 1 , i ( t A 1 ) , S ˙ A 1 , i ( t A 1 ) , S ¨ A 1 , i ( t A 1 ) , D I ( t A 1 ) ;
39. Calculate θ A 1 , i , m i n = m i n S A 1 , i ( t A 1 ) and θ A 1 , i , m a x = m a x S A 1 , i ( t A 1 ) ,   0 t A 1 T A 1 .
40. Calculate θ ˙ A 1 , i , m i n = m i n S ˙ A 1 , i ( t A 1 ) and θ ˙ A 1 , i , m a x = m a x S ˙ A 1 , i ( t A 1 ) ,   0 t A 1 T A 1 .
41. Calculate θ ¨ A 1 , i , m i n = m i n S ¨ A 1 , i ( t A 1 ) and θ ¨ A 1 , i , m a x = m a x S ¨ A 1 , i ( t A 1 ) ,   0 t A 1 T A 1 .
42. Calculate D I A 1 , m i n = m i n D I ( t A 1 ) ,   0 t A 1 T A 1
43.if  θ A 1 , i , m i n , θ A 1 , i , m a x D A 1 , i   & &   θ ˙ A 1 , i , m i n , θ ˙ A 1 , i , m a x D ˙ A 1 , i   & &  
θ ¨ A 1 , i , m i n ,   θ ¨ A 1 , i , m a x D ¨ ˙ A 1 , i   & &   D I A I , m i n D I t h  then
44.  Calculate f i t 2 ,
45.  if  f i t 2 n < f 2 p  then
46.        f 2 p , n = f i t 2 n ; x 2 p , n = x 2 n .
47.   end if
48.end if
49.  end for
50. if m i n ( f 2 p ) < f 2 g  then
51.   [ f 2 g , n u m ] = m i n ( f 2 p )   ; x 2 g = x 2 p , n u m .
52.end if
53. if f 2 g , k 1 f 2 g , k     0.001  then
54.count2 = count2 + 1
55.  else
56. count2 = 1
57.  end if
58. Update the current flying velocities and positions of all particles;
59.end while
60.Initialize: T A 2
61.while count3   15 do
62.  for  n = 1   to   N p  do
63. Determine S A 2 , i ( t A 2 ) , S ˙ A 2 , i ( t A 2 ) , S ¨ A 2 , i ( t A 2 ) , D I ( t A 2 ) ;
64. Calculate θ A 2 , i , m i n = m i n S A 2 , i ( t A 2 ) and θ A 2 , i , m a x = m a x S A 2 , i ( t A 2 ) ,   0 t A 2 T A 2
65. Calculate θ ˙ A 2 , i , m i n = m i n S ˙ A 2 , i ( t A 2 ) and θ ˙ A 2 , i , m a x = m a x S ˙ A 2 , i ( t A 2 ) ,   0 t A 2 T A 2
66. Calculate θ ¨ A 2 , i , m i n = m i n S ¨ A 2 , i ( t A 2 ) and θ ¨ A 2 , i , m a x = m a x S ¨ A 2 , i ( t A 2 ) ,   0 t A 2 T A 2
67. Calculate D I A 2 , m i n = m i n D I ( t A 2 ) ,   0 t A 2 T A 2
68.if  θ A 2 , i , m i n , θ A 2 , i , m a x D A 2 , i   & &   θ ˙ A 2 , i , m i n , θ ˙ A 2 , i , m a x D ˙ A 2 , i   & &  
θ ¨ A 2 , i , m i n ,   θ ¨ A 2 , i , m a x D ¨ ˙ A 2 , i   & &   D I A 2 , m i n D I t h  then
69.  Calculate f i t 3 ,
70.  if  f i t 3 n < f 3 p  then
71.        f 3 p , n = f i t 3 n ; x 3 p , n = x 3 n .
72.   end if
73.end if
74.  end for
75. if m i n ( f 3 p ) < f 3 g  then
76.   [ f 3 g , n u m ] = m i n ( f 3 p )   ; x 3 g = x 3 p , n u m .
77.end if
78. if f 3 g , k 1 f 3 g , k     0.001  then
79.count3 = count3 + 1
80.  else
81. count3 = 1
82.  end if
83. Update the current flying velocities and positions of all particles;
84.end while

Appendix B. Experimental Settings

For two example cases, the target position of the end effector (the target orientation is the same as the initial state) and the initial 3D locations of feature points of the TACCS and obstacles in frame 1 are given for the simulation implementation, and they are listed in Table A1 as follows.
Table A1. Target pose and the initial 3D locations of feature points of TACCS and obstacles in frame 1.
Table A1. Target pose and the initial 3D locations of feature points of TACCS and obstacles in frame 1.
StructuresCharacteristic PointsCoordinateRadius/mm (Number of Edges)
Structure TypeNumberX/mmY/mmZ/mm
RobotPolygon 11−200−200400-
2200−200400-
3200200400-
4−200200400-
5−200−200470-
6200−200470/
7200200470/
8−200200470/
Cylinder 2100470180
2001010180
Cylinder 31001010180
204051010180
Cylinder 4104051010198
207401010198
Cylinder 5105401010198
205402504198
Cylinder 6104052504198
207402504198
Cylinder 7104052504180
2002504180
Cylinder 81002504180
2003764180
Cylinder 910−2703764180
202703764180
Cylinder 101003764153
2004034153
Cylinder 111004034153
2004366153
Polygon 121−200−2004366/
2200−2004366/
32002004366/
4−2002004366/
5−200−2004516/
6200−2004516/
72002004516/
8−2002004516/
Cylinder 131004516100
2004816100
Cylinder 141004816100
202254816100
Cylinder 15104154816110
202254816110
Cylinder 16103004816110
203005646110
Cylinder 17104155646110
202255646110
Cylinder 181005646100
202255646100
Cylinder 191005646100
2006346100
Cylinder 2010−1506346100
201506346100
Cylinder 21100634685
200649685
Cylinder 22100649685
200668085
ObstaclesPolygonal obstacles 1 (Cases 1,2)1850450500/
21150450500/
31150750500/
4850750500/
58504501300/
611504501300/
711507501300/
88507501300/
Polygonal obstacles 2 (Cases 1,2)11840−660500/
22160−660500/
32160−340500/
41840−340500/
51840−6601500/
62160−6601500/
72160−3401500/
81840−3401500/
Cylindrical obstacles 3 (Cases 1,2)1800−500550200 ( N e   = 12)
2800−5001300200 ( N e   = 12)
Cylindrical obstacles 4 (Cases 1,2)11800500550200 ( N e   = 12)
218005001500200 ( N e   = 12)
Polygonal obstacles 5 (Case 1)12850−150700/
23150−150700/
33150150700/
42850150700/
52850−150850/
63150−150850/
73150150850/
82850150850/
Polygonal obstacles 6 (Case 1)13580−120700/
23820−120700/
33820120700/
43580120700/
53580−1201100/
63820−1201100/
738201201100/
835801201100/
Cylindrical obstacles 7 (Case 1)13000−550550150 ( N e   = 12)
23000−5501000150 ( N e   = 12)
Cylindrical obstacles 8 (Case 1)13000550550150 ( N e   = 12)
230005501000150 ( N e   = 12)
Polygonal obstacles 9 (Case 2)12850388.91601.04/
23150388.91601.04/
33150601.04388.91/
42850601.04388.91/
52850494.97707.11/
63150494.97707.11/
73150707.11494.97/
82850707.11494.97/
Polygonal obstacles 10 (Case 2)13580410.12579.83/
23820410.12579.83/
33820579.83410.12/
43580579.83410.12/
53580692.97862.67/
63820692.97862.67/
73580862.67692.97/
83580862.67692.97/
Cylindrical obstacles 11 (Case 2)130000777.82150 ( N e   = 12)
23000318.21096.02150 ( N e   = 12)
Cylindrical obstacles 12 (Case 2)13000777.820150 ( N e   = 12)
230001096.02318.2150 ( N e   = 12)
CaseTarget point130000950/
23000636.4676.4/

Appendix C. Experimental Results

Taking example case 1 as the objective, curves of angular parameters of each joint of TACCS and a series of motion poses relative to capsule based on the single optimization algorithm are shown in Figure A1 and Figure A2 as follows.
Figure A1. Curves of angular parameters of each joint of TACCS based on single optimization in example case 1.
Figure A1. Curves of angular parameters of each joint of TACCS based on single optimization in example case 1.
Mathematics 12 02245 g0a1
Figure A2. Series of motion poses relative to capsule based on single optimization in example case 1.
Figure A2. Series of motion poses relative to capsule based on single optimization in example case 1.
Mathematics 12 02245 g0a2
Moreover, taking example case 2 as the objective, a series of motion poses relative to capsule based on the bi-layer and single optimization algorithms are shown in Figure A3 and Figure A4 as follows.
Figure A3. Series of motion poses relative to capsule based on bi-layer optimization in example case 2.
Figure A3. Series of motion poses relative to capsule based on bi-layer optimization in example case 2.
Mathematics 12 02245 g0a3
Figure A4. Series of motion poses relative to capsule based on single optimization in example case 2.
Figure A4. Series of motion poses relative to capsule based on single optimization in example case 2.
Mathematics 12 02245 g0a4

References

  1. Kulakov, F.M. Some Russian research on robotics. Robot. Auton. Syst. 1996, 18, 365–372. [Google Scholar] [CrossRef]
  2. Thronson, H.; Akin, D.; Grunsfeld, J.; Lester, D. The Evolution and Promise of Robotic In-Space Servicing. In Proceedings of the AIAA Space 2009 Conference & Exposition, Pasadena, CA, USA, 14–17 September 2009. [Google Scholar]
  3. Jing, Z.; Qiao, L.; Pan, H.; Yang, Y.; Chen, W. An overview of the configuration and manipulation of soft robotics for on-orbit servicing. Sci. China Inf. Sci. 2017, 60, 050201. [Google Scholar] [CrossRef]
  4. Bluethmann, W.; Ambrose, R.; Diftler, M.; Askew, S.; Magruder, D. Robonaut: A Robot Designed to Work with Humans in Space. Auton. Robot. 2003, 14, 179–197. [Google Scholar] [CrossRef] [PubMed]
  5. Al-Kamil, S.J.; Szabolcsi, R. Enhancing Mobile Robot Navigation: Optimization of Trajectories through Machine Learning Techniques for Improved Path Planning Efficiency. Mathematics 2024, 12, 1787. [Google Scholar] [CrossRef]
  6. Dai, Y.; Xiang, C.F.; Zhang, Y.; Jiang, Y.P.; Qu, W.Y.; Zhang, Q.H. A Review of Spatial Robotic Arm Trajectory Planning. Aerospace 2022, 9, 361. [Google Scholar] [CrossRef]
  7. Shrivastava, A.; Dalla, V.K. Jerk Optimized Motion Planning of Redundant Space Robot Based on Grey-Wolf Optimization Approach. Arab. J. Sci. Eng. 2022, 48, 2687–2699. [Google Scholar] [CrossRef]
  8. Xie, Z.; Zhao, X.; Jiang, Z.; Yang, H.; Chongyang, L.I. Trajectory planning and base attitude restoration of dual-arm free-floating space robot by enhanced bidirectional approach. Front. Mech. Eng. Engl. Ed. 2022, 17, 16. [Google Scholar] [CrossRef]
  9. Li, Q.H.; Mu, Y.Q.; You, Y.; Zhang, Z.; Feng, C. A Hierarchical Motion Planning for Mobile Manipulator. IEEJ Trans. Electr. Electron. Eng. 2020, 15, 1390–1399. [Google Scholar] [CrossRef]
  10. Chen, H.L.; Zang, X.Z.; Liu, Y.B.; Zhang, X.H.; Zhao, J. A Hierarchical Motion Planning Method for Mobile Manipulator. Sensors 2023, 23, 6952. [Google Scholar] [CrossRef]
  11. Rastegarpanah, A.; Gonzalez, H.C.; Stolkin, R. Semi-Autonomous Behaviour Tree-Based Framework for Sorting Electric Vehicle Batteries Components. Robotics 2021, 10, 82. [Google Scholar] [CrossRef]
  12. Colucci, G.; Botta, A.; Tagliavini, L.; Cavallone, P.; Baglieri, L.; Quaglia, G. Kinematic Modeling and Motion Planning of the Mobile Manipulator Agri.Q for Precision Agriculture. Machines 2022, 10, 321. [Google Scholar] [CrossRef]
  13. Colucci, G.; Tagliavini, L.; Botta, A.; Baglieri, L.; Quaglia, G. Decoupled motion planning of a mobile manipulator for precision agriculture. Robotica 2023, 41, 1872–1887. [Google Scholar] [CrossRef]
  14. NithyaP, K.; Priya, P.S.; Benjamin, G.E.; Venkateswaran, J. Optimal Path Planning and Static Obstacle Avoidance for a Dual Arm Manipulator Used in On-Orbit Satellite Servicing. IFAC-Pap. 2020, 53, 189–194. [Google Scholar]
  15. Yang, H.; Li, D.; Xu, X.R.; Zhang, H. An Obstacle Avoidance and Trajectory Tracking Algorithm for Redundant Manipulator End. IEEE Access 2022, 10, 52912–52921. [Google Scholar] [CrossRef]
  16. Zhang, J.; Zhang, J.; Zhang, Q.; Wei, X. Obstacle Avoidance Path Planning of Space Robot Based on Improved Particle Swarm Optimization. Symmetry 2022, 14, 938. [Google Scholar] [CrossRef]
  17. Cui, L.L.; Wang, H.C.; Chen, W.D. Trajectory planning of a spatial flexible manipulator for vibration suppression. Robot. Auton. Syst. 2020, 123, 103316. [Google Scholar] [CrossRef]
  18. Song, Q.; Li, S.; Bai, Q.; Yang, J.; Zhang, A.; Zhang, X.; Zhe, L. Trajectory Planning of Robot Manipulator Based on RBF Neural Network. Entropy 2021, 23, 1207. [Google Scholar] [CrossRef] [PubMed]
  19. Liu, X.M.; Qiu, C.R.; Zeng, Q.F.; Li, A.P.; Xie, N. Time-energy Optimal Trajectory Planning for Collaborative Welding Robot with Multiple Manipulators. Procedia Manuf. 2020, 43, 527–534. [Google Scholar] [CrossRef]
  20. Cheng, Q.; Hao, X.L.; Wang, Y.; Xu, W.X.; Li, S.J. Trajectory planning of transcranial magnetic stimulation manipulator based on time-safety collision optimization. Robot. Auton. Syst. 2022, 152, 104039. [Google Scholar] [CrossRef]
  21. Chen, S.; Zhang, C.; Yi, J. Time-Optimal Trajectory Planning for Woodworking Manipulators Using an Improved PSO Algorithm. Appl. Sci. 2023, 13, 10482. [Google Scholar] [CrossRef]
  22. Sun, J.G.; Han, X.Y.; Zuo, Y.M.; Tian, S.Q.; Song, J.W.; Li, S.H. Trajectory Planning in Joint Space for a Pointing Mechanism Based on a Novel Hybrid Interpolation Algorithm and NSGA-II Algorithm. IEEE Access 2020, 8, 228628–228638. [Google Scholar] [CrossRef]
  23. Zhang, T.; Cheng, J.; Zou, Y.B. Time-optimal and Smooth Trajectory Planning for Multi-axis Motion Systems Based on ISC Similarity. Int. J. Control. Autom. Syst. 2024, 22, 1238–1251. [Google Scholar] [CrossRef]
  24. Ma, J.F.; Gao, S.; Yan, H.Y.; Lv, Q.; Hu, G.Q. A new approach to time-optimal trajectory planning with torque and jerk limits for robot. Robotics Auton. Syst. 2021, 140, 103744. [Google Scholar] [CrossRef]
  25. Liu, Z.F.; Xu, J.J.; Yang, C.B.; Zhao, Y.S.; Zhang, T. A TE-E Optimal Planning Technique Based on Screw Theory for Robot Trajectory in Workspace. J. Intell. Robot. Syst. 2018, 91, 363–375. [Google Scholar] [CrossRef]
  26. Zhou, Y.; Han, G.; Wei, Z.; Huang, Z.; Chen, X. A Time-Optimal Continuous Jerk Trajectory Planning Algorithm for Manipulators. Appl. Sci. 2023, 13, 11479. [Google Scholar] [CrossRef]
  27. Pu, Q.C.; Xu, X.R.; Li, Q.Q.; Zhang, H. Robotic arm time–jerk optimal trajectory based on improved dingo optimization. J. Braz. Soc. Mech. Sci. Eng. 2024, 46, 198. [Google Scholar] [CrossRef]
  28. Kashima, T.; Isurugi, Y. Trajectory Planning of Manipulators Based on a Minimum-Energy Criterion and Operating Time. J. Robot. Soc. Jpn. 2010, 15, 1012–1018. [Google Scholar] [CrossRef]
  29. Liu, K.P.; Sui, J.L.; Yue, N.; Liu, S.S. Path planning method of mobile manipulator based on the representation space. In Proceedings of the 2016 IEEE International Conference on Mechatronics and Automation, Harbin, China, 7–10 August 2016; pp. 322–326. [Google Scholar]
  30. Luna, R.; Moll, M.; Badger, J.; Kavraki, L. A scalable motion planner for high-dimensional kinematic systems. Int. J. Robot. Res. 2020, 39, 361–388. [Google Scholar] [CrossRef]
  31. Zhang, W.M.; Fu, S.X. Time-optimal Trajectory Planning of Dulcimer Music Robot Based on PSO Algorithm. In Proceedings of the Chinese Control and Decision Conference (CCDC), Hefei, China, 22–24 August 2020; pp. 4769–4774. [Google Scholar]
  32. Dai, J.; Zhang, Y.; Deng, H. Novel Potential Guided Bidirectional RRT* with Direct Connection Strategy for Path Planning of Redundant Robot Manipulators in Joint Space. IEEE Trans. Ind. Electron. 2024, 71, 2737–2747. [Google Scholar] [CrossRef]
  33. Chen, G.; Luo, N.; Liu, D.; Zhao, Z.H.; Liang, C.C. Path planning for manipulators based on an improved probabilistic roadmap method. Robot. Comput. Integr. Manuf. 2021, 72, 102196. [Google Scholar] [CrossRef]
  34. Niu, P.; Cheng, Q.; Chen, C.; Yang, C.; Liu, Z. An approach for crucial geometric error analysis and accuracy enhancement of gantry milling machines based on generalized correlation sensitivity. J. Manuf. Process. 2024, 119, 401–413. [Google Scholar] [CrossRef]
  35. Shi, B.H.; Xu, J.X. Time-Optimal Trajectory Planning of Industrial Robot based on Improved Particle Swarm Optimization Algorithm. In Proceedings of the 39th Chinese Control Conference (CCC), Shenyang, China, 27–29 July 2020; pp. 3683–3688. [Google Scholar]
  36. Yu, X.L.; Dong, M.S.; Yin, W. Time-optimal trajectory planning of manipulator with simultaneously searching the optimal path. Comput. Commun. 2021, 181, 446–453. [Google Scholar] [CrossRef]
  37. Jin, R.Y.; Rocco, P.; Geng, Y. Cartesian trajectory planning of space robots using a multi-objective optimization. Aerosp. Sci. Technol. 2021, 108, 106360. [Google Scholar] [CrossRef]
  38. Li, R.; Liu, M.; Teutsch, J.; Wollherr, D. Constraint trajectory planning for redundant space robot. Neural Comput. Appl. 2023, 35, 24243–24258. [Google Scholar] [CrossRef]
  39. Cao, X.; Yan, H.; Huang, Z.; Ai, S.; Xu, Y.; Fu, R.; Zou, X. A Multi-Objective Particle Swarm Optimization for Trajectory Planning of Fruit Picking Manipulator. Agronomy 2021, 11, 2286. [Google Scholar] [CrossRef]
  40. Ding, X.L.; Wang, Y.C.; Wang, Y.B.; Xu, K. A review of structures, verification, and calibration technologies of space robotic systems for on-orbit servicing. Sci. China Technol. Sci. 2021, 64, 462–480. [Google Scholar] [CrossRef]
  41. Xu, J.J.; Liu, Z.F.; Cheng, Q.; Zhao, Y.S.; Pei, T.H.; Yang, C.B. Models for Three New Screw-based IK Sub-problems Using Geometric Descriptions and Their Applications. Appl. Math. Model. 2018, 67, 399–412. [Google Scholar] [CrossRef]
  42. Xu, J.J.; Liu, Z.F.; Zhang, C.X.; Yang, C.B.; Pei, Y.H. Minimal distance calculation between the industrial robot and its workspace based on circle/polygon-slices representation. Appl. Math. Model. 2020, 87, 691–710. [Google Scholar] [CrossRef]
Figure 1. Classic structure of TACCS.
Figure 1. Classic structure of TACCS.
Mathematics 12 02245 g001
Figure 2. Kinematics parameters of each arm of TACCS.
Figure 2. Kinematics parameters of each arm of TACCS.
Mathematics 12 02245 g002
Figure 3. Motions of two IK sub-problems.
Figure 3. Motions of two IK sub-problems.
Mathematics 12 02245 g003
Figure 4. Envelope structures of TACCS and capsule.
Figure 4. Envelope structures of TACCS and capsule.
Mathematics 12 02245 g004
Figure 5. Envelop of a circle.
Figure 5. Envelop of a circle.
Mathematics 12 02245 g005
Figure 6. Errors of the shortest distance under different values of ϵ .
Figure 6. Errors of the shortest distance under different values of ϵ .
Mathematics 12 02245 g006
Figure 7. Framework of the feasible space analysis based on IK functions of two arms.
Figure 7. Framework of the feasible space analysis based on IK functions of two arms.
Mathematics 12 02245 g007
Figure 8. Bi-layer optimization framework.
Figure 8. Bi-layer optimization framework.
Mathematics 12 02245 g008
Figure 9. Feasible spaces in example case 1.
Figure 9. Feasible spaces in example case 1.
Mathematics 12 02245 g009
Figure 10. Feasible spaces in example case 2.
Figure 10. Feasible spaces in example case 2.
Mathematics 12 02245 g010
Figure 11. Iterative curves of the sub-optimization 1 in bi-layer optimization and the single optimization under different values of the proportion (50%, 70%, 90% and 100%).
Figure 11. Iterative curves of the sub-optimization 1 in bi-layer optimization and the single optimization under different values of the proportion (50%, 70%, 90% and 100%).
Mathematics 12 02245 g011
Figure 12. Iterative curves for the single optimization and three sub-optimizations of the bi-layer optimization with different values of PSO parameters.
Figure 12. Iterative curves for the single optimization and three sub-optimizations of the bi-layer optimization with different values of PSO parameters.
Mathematics 12 02245 g012
Figure 13. Variation in the shortest distance between TACCS and its obstacles over motion time based on two algorithms.
Figure 13. Variation in the shortest distance between TACCS and its obstacles over motion time based on two algorithms.
Mathematics 12 02245 g013
Figure 14. Curves of angular parameters of each joint of TACCS based on bi-layer optimization.
Figure 14. Curves of angular parameters of each joint of TACCS based on bi-layer optimization.
Mathematics 12 02245 g014
Figure 15. Series of motion poses relative to capsule based on bi-layer optimization.
Figure 15. Series of motion poses relative to capsule based on bi-layer optimization.
Mathematics 12 02245 g015
Table 1. Comparison of some representative works wherein the superscript #F means the robot type is free-floating space robot.
Table 1. Comparison of some representative works wherein the superscript #F means the robot type is free-floating space robot.
WorksRobot TypeDOFsPlanning MethodsSafetyStabilityEfficiencyEnergy
Consumption
Consideration
of Motion Coupling
Ref. [9]Mobile
manipulator
3 + 6A* algorithm; Uniform
sampling-based algorithm
YesNoNoNoNo
Ref. [10]3 + 6A* algorithm; Sampling-based
heuristic algorithm
YesNoNoNoNo
Ref. [11]3 + 6A* algorithm; RRTYesNoNoNoNo
Refs. [12,13]5 + 6A* algorithm; RRTYesNoNoNoYes
Ref. [29]3 + 2 A* algorithmYesNoNoNoYes
Ref. [39]3 + 6PSO-based optimizationYesYesYesYesNo
Ref. [30]High-DOFs
robot
14 or 21Sampling-based algorithmYesNoNoNo/
Ref. [31]5PSO-based optimizationNoYesYesNo/
Ref. [32]AnyBidirectional RRT*YesNoNoNo/
Ref. [33]AnyPRMYesNoNoNo/
Ref. [37]6 #FCPSO-based optimizationNoYesYesNo/
Ref. [38]7 #FPSO-WOA-based optimizationYesYesYesNo/
The proposed methodTACCS6 + 6PSO-based optimizationYesYesYesYesYes
Table 2. Kinematics parameters of two arms.
Table 2. Kinematics parameters of two arms.
Parameters l 1   (mm) l 2   (mm) l 3   (mm) l 4   (mm) l 5   (mm) l 6   (mm)
First-level arm54054014941260270332
Second-level arm300300830700150184.5
Table 3. Comparison of convergence times of two optimization algorithms.
Table 3. Comparison of convergence times of two optimization algorithms.
Example CaseOptimization μ / s σ / s c v / %
Case 1Single145.5451.6735.51
Bi-layer40.1210.5826.39
Difference72.43%79.52%25.68%
Case 2Single52.4818.6635.55
Bi-layer19.083.2617.10
Difference63.64%82.53%51.90%
Table 4. Constraints, optimal results of main parameters based on three algorithms, including the bi-layer opt. (bi-layer optimization), single opt. (single optimization) and trad. opt. (traditional optimization).
Table 4. Constraints, optimal results of main parameters based on three algorithms, including the bi-layer opt. (bi-layer optimization), single opt. (single optimization) and trad. opt. (traditional optimization).
Example CaseParameterConstraintsBi-Layer Opt.Single Opt.Difference 1Trad. Opt.Difference 2
Case 1 Δ x / ( mm ) [1500, 3000]2016.081628.35-2500.00-
Δ y / ( mm ) [−1000, 1000]604.94675.01-815.81-
Δ z / ( mm ) [820, 2200]2148.841797.81-2097.21-
α / ( deg ) [−180, 180]92.66−24.45-104.74-
β / ( deg ) [−90, 90]76.6564.71-5.24-
γ / ( deg ) [−180, 180]86.732.64-90.00-
V m a x / ( deg / s ) ≤12077.7164.51-55.59-
A m a x / ( deg / s 2 ) ≤200103.8095.67-54.13-
T   ( s ) ≤204.674.986.22%4.513.55%
P a t h   ( deg ) -9.829.068.39%11.0611.21%
F-14.4914.043.21%15.576.94%
Case 2 Δ x / ( mm ) [1500, 2800]1871.381652.65-1653.14-
Δ y / ( mm ) [−200, 1800]1102.911246.74-1214.89-
Δ z / ( mm ) [800, 1880]1394.951521.54-1619.18-
α / ( deg ) [−180, 180]22.37−77.82-67.21-
β / ( deg ) [−90, 90]72.6125.74-45.72-
γ / ( deg ) [−180, 180]17.02−72.95-59.02-
V m a x / ( deg / s ) ≤12095.9796.94-81.29-
A m a x / ( deg / s 2 ) ≤200136.82138.06-115.61-
T   ( s ) ≤204.214.373.66%5.4522.75%
P a t h   ( deg ) -11.3711.420.44%10.844.89%
F-15.5815.791.33%16.294.36%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xu, J.; Tao, L.; Pei, Y.; Cheng, Q.; Chu, H.; Zhang, T. Collision-Free Trajectory Planning Optimization Algorithms for Two-Arm Cascade Combination System. Mathematics 2024, 12, 2245. https://doi.org/10.3390/math12142245

AMA Style

Xu J, Tao L, Pei Y, Cheng Q, Chu H, Zhang T. Collision-Free Trajectory Planning Optimization Algorithms for Two-Arm Cascade Combination System. Mathematics. 2024; 12(14):2245. https://doi.org/10.3390/math12142245

Chicago/Turabian Style

Xu, Jingjing, Long Tao, Yanhu Pei, Qiang Cheng, Hongyan Chu, and Tao Zhang. 2024. "Collision-Free Trajectory Planning Optimization Algorithms for Two-Arm Cascade Combination System" Mathematics 12, no. 14: 2245. https://doi.org/10.3390/math12142245

APA Style

Xu, J., Tao, L., Pei, Y., Cheng, Q., Chu, H., & Zhang, T. (2024). Collision-Free Trajectory Planning Optimization Algorithms for Two-Arm Cascade Combination System. Mathematics, 12(14), 2245. https://doi.org/10.3390/math12142245

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop