Abstract
Aiming at reducing joint velocity jumps caused by an unexpected joint-locked failure during space manipulator on-orbit operations without shutting down manipulator, trajectory optimization strategy considering the unexpectedness characteristics of joint-locked failure is proposed in the paper, which can achieve velocity jumps reduction in both operation space and joint space simultaneously. In the strategy, velocity in operation space concerning task completion directly is treated as equality constraints, and velocity in joint space concerning motion performance is treated as objective function. Global compensation vector which consists of coefficient, gradient of manipulability, and orthogonal matrix of null space is constructed to minimize the objective function. For each particular failure time, unique optimal coefficient can be obtained when the objective function is minimal. As a basis, a method for optimal coefficient function fitting is proposed based on a priori failure information (possible failure time and the corresponding optimal coefficient) to guarantee the unexpectedness characteristics of joint-locked failure. Simulations are implemented to validate the efficiency of trajectory optimization strategy in reducing velocity jumps in both joint space and operation space. And the feasibility of coefficient function is also verified in reducing velocity jump no matter when joint-locked failure occurs.
1. Introduction
Space manipulator [1, 2] is a kind of complex electromechanical system, which plays an important role in space station maintenance, rendezvous and docking, and space exploration. During its long-time service on orbit, affected by radiation, energetic particle, and high-low temperature alternating of space environment, many kinds of malfunctions are inevitable, such as link deformation, performance degradation, and joint failure [3, 4]. Joint as a precision unit, integrating with machinery, communications, and automatic control technology, is the core component of space manipulator to achieve safe and smooth operation. Joint failure may bring about performance degradation of space manipulator and lead to bad effect on task completion.
For general manipulator on the ground, when joint fails during manipulator motion, the manipulator can be shut down to maintain or replace the failure joint. However, if shutting down space manipulator once joint failure occurs, the inertia force acting on the manipulator will cause serious impact force on the end-effector and joints, as well as cause instability of the base. When joint fails especially during load operations such as space cabin docking [5] and astronauts assistant [6], the bigger mass and inertia will cause huger impact, which may exceed the rated torque and cause serious hurt on the manipulator parts and operating objects. On the other hand, when trying to execute task continuously after joint failure, control strategy should be adjusted in real time. As a result, velocity jumps will be caused in joint space and operation space, which will also generate impact on manipulator structure. Thus, it is of great significance to smooth motion parameters when joint failure occurs during on-orbit operations no matter whether space manipulator will be shut down or execute task continuously.
Joint failure can be divided into joint-locked [7] and free swing [8] according to whether joint can be operated or not after failure occurs. In free swing state, failure joint moves freely, making the manipulator out of control. In joint-locked state, failure joint is locked at the angle of failure time. Although the degree of freedom (DOF) of space manipulator is reduced, the healthy joints can move normally, and the task can be executed continuously when joint-locked failure occurs. Thus in order to avoid the impact force caused by shutting down space manipulator after joint-locked failure, the paper is devoted into trajectory optimization research to achieve task completion and motion parameters smoothness (velocity jumps reduction especially) without shutting down manipulator.
In order to achieve task completion after manipulator failure, scholars concentrate on the following researches. First, in order to obtain the failure information such as failure time and failure joint, fault detection, fault isolation, and fault identification [9–12] are researched. Second, abilities of manipulator are analyzed in failure condition such as fault-tolerant workspace [13, 14] and dexterity [15, 16]. Based on these researches above, pseudoinverse method [17], self-motion manifold method [18], quadratic programming algorithm [19], and configuration optimization [20, 21] are applied to achieve task completion after failure occurs. Meanwhile, constraints are analyzed during fault tolerance to guarantee the practical application. Ralph and Pai [22] analyzed the minimal constraints necessary in fault tolerance, Jamisola Jr. et al. [23] took environment obstacle constraints into account, and Xie et al. [24] introduced the singularity constraints into the fault-tolerant method.
In the mentioned fault-tolerant researches, some of them shut down manipulator after failures; then fault-tolerant method is introduced to achieve task completion. However, it is not suitable for space manipulator fault-tolerance. The others execute task via fault-tolerant method in real time. Unfortunately, when guaranteeing velocity and position in operation space for task completion, velocity jumps are brought about into joint space. Scholars try to reduce the velocity jumps; Abdi et al. [25] concentrated on the minimum velocity jump in operation space and proposed an optimization strategy for joint velocity redistribution. Meanwhile, the velocity jump in joint space is reduced with the application of least square method and matrix perturbation method [26]. Jing and Cheng [27] achieved velocity jump reduction in joint space by optimizing the initial configuration of a dual-arm manipulator. In the authors’ former work [28], velocity jumps in joint space were reduced by constructing compensation vectors based on the gradient of manipulability. Further, the coefficient form of compensation vector was discussed to achieve better reduction performance.
However, simultaneous reduction for velocity jumps in operation space and joint space is not achieved. What is more, velocity jumps reduction in the operation space may even increase the velocity jumps in joint space. Since velocity jumps in operation space are concerned with task completion and velocity jumps in joint space are concerned with joint structure safety, it is of equal importance to achieve velocity jumps reduction in the two spaces to guarantee task completion and motion safety. According to the self-motion characteristics of space manipulator, a set of velocities in operation space corresponds to multiple sets of joint velocities. So when the velocity in operation space is fixed, joint velocity can be adjusted to find an optimal value minimizing the jumps. On the basis, the trajectory optimization strategy can be formed by treating velocity in operation space as constraints and treating velocity in joint space as optimization objective.
Besides, the mentioned fault-tolerant methods depend on failure information; only after the failure condition is confirmed can they begin to work. However, in practical application, joint-locked failure occurs unexpectedly. It needs time to obtain failure information with fault detection method, so time delay is inventible between failure discovery and fault-tolerant method application, which makes fault-tolerant performance worse than the ideal condition. Although controllers are designed based on neural networks to achieve real-time fault tolerance [29–31], the adaptation of controllers for rapid response and stable operation is not enough for on-orbit application. Therefore, in order to avoid time delay problem of fault detection method and to guarantee the unexpectedness characteristics of joint-locked failure, global compensation for joint velocity during the entire path planning is proposed in the paper, based on which the trajectory optimization is not relevant to the failure information and can handle an unexpected joint-locked failure.
In conclusion, the trajectory optimization strategy proposed in the paper aims at dealing with three problems. First, the task completion on-orbit is guaranteed without shutting down space manipulator, and the parameters of space manipulator particularly the velocity jumps should be smoothed after joint-locked failure occurs. Second, the reduction for velocity jumps in both operation space and joint space should be achieved simultaneously. It can be achieved by considering velocity in operation space as equality constraints and considering velocity in joint space as optimization objective in the strategy. Third, the avoidance of time delay problem and the unexpectedness characteristics of joint-locked failure should be handled. For this purpose, a global compensation vector is constructed for velocity compensation over the entire path planning task, which makes the optimization strategy irrelevant with the real-time failure information acquisition. In addition, with global compensation for joint velocity, a unique coefficient which makes objective function minimal can be obtained related to a particular failure time. Taking it as a priori failure information, the law of optimal coefficient related to arbitrary failure time is concluded. A coefficient function is obtained to deal with velocity jump reduction caused by unexpected joint-locked failure.
The presented paper is organized as follows: in Section 2, velocity jumps in operation space and joint space caused by joint-locked failure are analyzed. The mathematical expression of trajectory optimization strategy is derived, during which velocity jumps in operation space are expressed as equality constraint and velocity jumps in joint space are expressed as objective function. In Section 3, global compensation vector for joint velocity is constructed, and the solving method is proposed to obtain the optimal coefficient of global compensation when joint-locked failure occurs at a particular time. Then considering the unexpectedness characteristics of joint-locked failure, coefficient of global compensation vector is fitted in function form based on a priori failure information. In Section 4, simulations are implemented to validate the effectiveness of trajectory optimization strategy in reducing the velocity jump in both joint space and operation space when joint-locked failure occurs at arbitrary time. The last part is conclusion.
2. Mathematical Expression of Trajectory Optimization Strategy
Joints are the key parts of space manipulator operations. When joint-locked failure occurs, it will bring about velocity jumps in operation space and joint space, which would tremendously reduce the operability of space manipulator, especially the ability of pose-reaching and task-completing. In this section, trajectory optimization is proposed considering simultaneous velocity jumps reduction in two spaces. Based on the analysis of velocity jumps, the constraints and objective of trajectory optimization are derived.
Trajectory optimization is discussed on a typical space manipulator path planning task, which starts at time and ends at time . The control interval is represented as , and the execution time of path planning can be represented as , , .
2.1. Analysis on Velocity Jumps in Operation Space and Joint Space
For a degree space manipulator, when joints are healthy during path planning task, the Jacobian matrix and joint velocity vector at time can be expressed as
The Jacobian matrix is the function of joint angle vector at current time . And represents the dimension of operation space. represents the th column of Jacobian matrix.
When joint-locked failure occurs in joint at time , the dimension of joint space degenerates to . Correspondingly, the Jacobian matrix degenerates to , which is named as reduced Jacobian matrix [32]. The reduced Jacobian matrix and joint velocity vector at time can be expressed as
Correspondingly, represents the th column of reduced Jacobian matrix. And represents the reduced joint angle vector. After joint fails at time , the velocity of failure joints turns to 0, and velocity of healthy joints keeps constant. Thus, at time , it has
Then the velocity jump vector in joint space at time can be expressed as
represents the nominal velocity of th joint when the joint is healthy. On the basis, the velocity jump in operation space caused by joint-locked failure is expressed as
According to (5), it will lead to accumulative pose (position and orientation) deviation in operation space:
When joint failure occurs unexpectedly during path planning task, an average velocity deviation in the operation space can be obtained:
Equation (7) can be used to evaluate the influence of velocity deviation caused by an unexpected joint-locked failure, while, for a general path planning task, the pose accuracy should fulfill accuracy threshold at each time to guarantee task completion. Set the threshold value as ; when , the task cannot be completed. Therefore, the velocity jump in operation space should be eliminated preferentially to guarantee task completion.
Generally, velocity jump in operation space can be reduced by compensation for healthy joint. When compensation vector is introduced, (5) is turned into
When the velocity deviation in operation space is completely eliminated, namely, , the compensation vector for joint velocity can be expressed as
represents the pseudoinverse of the reduced Jacobian matrix. By given the compensation as (9), the velocity jump in operation space can be completely eliminated. The contribution of failure joint is redistributed to healthy joints. However, additional velocity jumps are introduced into healthy joints. The velocity compensation may be too large for safety threshold of joint velocity and the impact caused by compensation may damage the joint structure, leading to terribly concatenate failure of joints. Thus, it is of equal importance to decrease the velocity jumps in operation space and joint space.
2.2. Constraints and Objective of Trajectory Optimization
In order to achieve simultaneous velocity jumps reduction in both operation space and joint space, the velocity in operation space is considered as equality constraint, which is preferentially guaranteed for task completion. And joint velocity is used to construct the objective function, which is adjusted to achieve joint velocity jumps minimal. The constraints and objective function of trajectory optimization can be expressed as follows.
2.2.1. Constraints
Assume only position accuracy in operation space is considered during path planning task. The actual value of linear velocity in operation space should be in accordance with the nominal value; namely,
represent the actual and nominal value of the linear velocity in the operation space. Correspondingly, the position in operation space should be tracked as the nominal value:
represent the actual and nominal position in operation space, respectively. Equations (10) and (11) are named as task constraints. When task constraints are guaranteed, it has (; ).
Besides, the parameters of space manipulator should fulfill the thresholds in path planning. Constraints of joint velocity and joint angular acceleration are taken into account:wherein represents the number of joints and and represent the threshold value of joint velocity and angular acceleration of the th joint, respectively. When (12) is fulfilled, it has .
2.2.2. Objective Function
The objective function is constructed as the joint velocity deviation at failure time . Since the dimension of joint velocity is reduced after joint failure, in order to guarantee the dimension, is defined based on :
Then the objective function can be expressed as
Based on the analysis above, the mathematical expression of trajectory optimization strategy is derived as
3. Solution for Trajectory Optimization considering Unexpectedness Characteristics of Joint-Locked Failure
In this section, trajectory optimization when joint fails at a particular time is solved firstly. A global compensation vector for joint velocity is constructed based on the orthogonal basis of Jacobian matrix null space. Then the optimal coefficient of compensation vector is obtained when the objective function is minimal. Since global compensation vector compensates joint velocity in each control interval ignoring the failure time, time delay problem can be efficiently avoided. Taking the failure time and corresponding optimal coefficient as a priori failure information, the law of optimal coefficient related to arbitrary failure time is concluded, and a method for coefficient function fitting is proposed based on a priori failure information to guarantee the unexpectedness characteristics of joint-locked failure.
3.1. Global Compensation Vector for Joint Velocity
In order to construct the global compensation vector, the orthogonal basis of Jacobian matrix null space is used to represent the self-motion characteristics of space manipulator. Manipulability [33] is introduced into the vector, whose gradient can be used to represent the velocity compensation weight of each joint. When joint-locked failure occurs, it is reduced as . Besides, a constant coefficient is used to adjust the velocity compensation rate for each joint. Then for each time, it has
In (16a) and (16b), and are used to compensate joint velocity at each time in path planning. Before failure occurs, precompensation is implemented with the vector of (16a). After failure occurs, joint velocity is compensated according to (16b). Namely, joint velocity at each time corresponds to a vector or . Then the global compensation vector for joint velocity over the entire task is constructed:
With the compensation of , joint velocity can be expressed as , which is detailed in
Joint angle at time can be obtained as
By substituting and into (16a) and (16b), and can be obtained. The recurrence relation between , , and is generated, which is shown in Figure 1.
Since the dimension of velocity is reduced after joint failure, similar to (13), is defined as the same dimension as , whose th element equals 0. Then, according to (18), the objective function can be turned into
For a particular path planning task, when failure time is certain, the global compensation vector will change with coefficient . Thus the objective function can be optimized to achieve joint velocity jump minimal by finding the optimal .
3.2. Calculation for the Optimal Coefficient
In order to find the optimal , the gradient of reduced manipulability in global compensation vector should be calculated firstly. It haswherein the partial derivative of the reduced manipulability to the th joint angle can be expressed as
Make , and expresses the th row of . Thenwherein represents the th row of and is the dimension of operation space. Then
In (24), . Substitute (24) into (23); the partial derivative of reduced manipulability to joint angle can be obtained:
Substitute (25) into (21); the gradient of the reduced manipulability can be obtained. The gradient of manipulability can be obtained in a similar way. Then, by given , the value of global compensation vector can be calculated.
Before searching for the optimal , its feasible region should be preferentially determined. According to constraints (10) and (11), the feasible region for can be defined aswherein the subscripts represent the component of vector .
Equation (23) means the component of deviation between actual position and desired position should fulfill the threshold value of position accuracy simultaneously, where is expressed as . In the feasible region, the optimal can be obtained as follows:(a)According to the constraints, obtain the feasible region of .(b)Set the initial value of and interval value , and calculate the global compensation vector .(c)Calculate the joint velocity and acceleration compensated with .(d)Confirm whether joint velocity and acceleration guarantee the constraints (12). If the constraints are not guaranteed, the corresponding is abandoned. Else the value of objective function is obtained . is stored in set .(e). If , turn to step (c). Else select the minimal value of in set , and the corresponding is the optimal value.
The steps are also detailed as shown in Figure 2.
3.3. The Method of Fitting Optimal Coefficient Function
With the analysis above, the optimal can be obtained to minimize the objective function. However, during practical application, the accurate failure time is not easy to be obtained. In this section, the trajectory optimization is developed considering the unexpectedness characteristics of joint-locked failure when failure information is unknown. And a method is proposed to derive the law of optimal coefficient for unexpected joint failure.
Since failure time is unexpected in practical application, the objective function is changed as the maximal value of joint velocity deviation between adjacent control cycles:wherein represents the joint velocity deviation, which is expressed as
According to the analysis in Sections 3.1 and 3.2, a unique can be obtained to minimize the objective function when joint-locked failure occurs at particular time . Thus, when the failure time is unknown, a function of operation time can be constructed to guarantee the arbitrary failure time. Then (18) can be turned into
As a basis, the trajectory optimization strategy considering the unexpectedness characteristics of joint-locked failure can be expressed as
In order to solve the optimization strategy and obtain the coefficient function , a method of fitting the coefficient function based on polynomial fitting using a priori failure information is proposed. Define the polynomial function as follows:
represents the coefficient of polynomial, represents the order of polynomial, and represents the task execution time. Generally, joint velocity at the initial time and the end time should be zero during path planning; namely, and . According to , function should guarantee the boundary condition:
In order to obtain the polynomial coefficient , a priori failure information is used to fit function (31). The a priori failure information includes failure time, failure joint angle, and the corresponding optimal coefficient . When failure time is certain, the failure joint angle can be theoretically calculated via manipulator path planning algorithm, and the optimal can be obtained based on the analysis in Section 3.2. In this way, a priori failure information can be obtained by supposing possible failure time in path planning. Then numbers of failure time and the corresponding optimal coefficient can be obtained as samples for function fitting.
The samples consist of failure time and the corresponding optimal coefficient . The th sample can be represented as . and are the th element of and , which can also be represented as and , respectively. Thus, the th sample can be represented as .
Define the evaluation criterion of function fitting as the sum of square for deviations between samples and fitting function value, which has the expression of
The threshold of evaluation criterion is defined as . When , function has a high fitting quality, and the fitting result is acceptable, while means more samples should be selected or polynomial order should be increased for better fitting. In conclusion, the process of fitting function is shown in Figure 3.
By considering the unexpectedness characteristics of joint-locked failure, a method is proposed based on a priori failure information for optimal coefficient function fitting. So when accurate joint failure information cannot be acquired in practical application, global compensation vector based on can be introduced to achieve trajectory optimization, with which velocity jumps can be efficiently reduced no matter when joint-locked failure occurs.
4. Simulation Experiments
4.1. Research Object
An SSRMS [2] type space manipulator who has 7 DOF is used to carry out simulation experiments. The coordinate systems are shown in Figure 4, and the DH parameters are listed in Table 1. Path planning task mentioned in the paper concerns the position accuracy in the operation space of space manipulator only.
4.2. Velocity Jump Reduction When Joint-Locked Failure Occurs at a Particular Time
In a particular path planning task, the initial configuration of space manipulator is set as and the target pose is set as . The other parameters of path planning task are listed in Table 2.
According to Sections 3.1 and 3.2, under the position threshold of , the feasible region of is calculated as . Then set the initial value of as 0, and carry out search algorithm with interval in the feasible region. The value of objective function changing with is shown in Figure 5. When , the objective function is minimal; .
When , the joint velocities are shown in Figure 6. Figure 6(a) represents the velocity curve of joint 2 (failure joint), and Figure 6(b) represents the velocity curve of joint 3 (healthy joint). It is obvious that the joint velocity jump is significantly reduced after global compensation. In detail, the value of velocity jumps, acceleration, and jump reduction rate of each joint is listed in Table 3.
(a) Joint 2
(b) Joint 3
In Table 3, velocity jump of each joint is reduced significantly except joint 5. Since objective function reflects the accumulation of velocity jump of all the joints, when objective function is minimal, it does not mean the velocity jump of each joint is minimal. The velocity jumps of some joints may increase, such as joint 5 in Table 3. In a former paper of the author [28], the coefficient is changed into diagonal matrix form to achieve velocity jump minimal for each joint. However, it brings about higher computation cost and more complex solution process. With constant coefficient of , the norm of velocity jumps is reduced , which shows great efficiency of the global compensation vector constructed with constant coefficient in reducing joint velocity jumps. Although jumps of some joints (joint 5) increase, the velocity jump is very small and can be accepted in practical application.
Besides, the position deviation in operation space when is obtained in Figure 7. With the compensation in joint velocity at each control interval, the position deviation in operation space is accumulated. The maximum deviation is less than , which is much smaller than the length of manipulator ( order of magnitude). During the 20 s’ path planning task, the velocity jump in operation space is about order of magnitude, which can be basically ignored. In conclusion, the velocity jumps in operation space can be mostly eliminated, while the velocity jump in joint space can be reduced to minimal when joint-locked failure occurs at a particular time.
4.3. Coefficient Function Fitting for Unexpected Joint-Locked Failure
Define the initial configuration of space manipulator path planning task as , and set the target pose as , , , , , . The velocity in operation space is planned by trapezoidal planning with parabola transition. Task cycle is 12 s, acceleration time is 2 s, and control interval is 50 ms. Joint-locked failure is assumed to occur in joint 2.
Since joint-locked failure time is unexpected, a priori failure information is used to fit the coefficient function. A priori failure information includes possible failure time, failure joint angle, and the corresponding optimal . During the mentioned path planning task, 40 moments are selected uniformly as possible failure time. According to (18) and (19), the joint angle at failure time can be obtained. Correspondingly, the optimal related to each moment can be obtained, which is shown in Figure 8. Regardless of the samples , the other samples can be used to fit a five-order polynomial function when meeting the threshold of fitting :
Considering the samples , the coefficient function for an unexpected joint-locked failure can be expressed as
and are segment points. The criterion to determine the segment points is defined as follows: the values of objective function compensated with and are the same; namely,
represents infinitesimal. However, for each possible failure moment , there is a feasible region for the related coefficient . If selected according to (36) guarantees , it should be abandoned. Instead of it, the segment points are set as .
Accordingly, two segment points can be obtained as , . Then (35) can be expressed as
Based on (37), the compensation vector of each control interval can be obtained, with which the velocity jump can be reduced no matter when joint-locked failure occurs.
4.4. Verification of the Coefficient Function Related to Unexpected Joint-Locked Failure
In order to verify the efficiency of the optimal coefficient function (37), assume joint-locked failure occurs at each control interval; namely, failure occurs at . Via the velocity compensation vector based on (37), the velocity jumps related to each failure time are shown in Figure 9.
At most failure time, the joint velocity jump can be significantly reduced. The average velocity jump equals after compensation. Compared with the average value of before compensation, velocity jump is reduced averagely. Especially in period of , the performance of jump reduction is remarkable, and the velocity jump is reduced below .
Before velocity jump reduction, during a control interval , joint acceleration needs to change averagely when joint-locked failure occurs. In order to afford the acceleration, huge joint torque is needed. When the safe threshold of joint acceleration is set as to protect joint torque not exceeding the rated torque, the change of velocity cannot be achieved and huge impact will be affected on joint motor. After velocity jump reduction, the requirement for joint acceleration is reduced to averagely, which fulfills the safe threshold and the impact on the joint can be eliminated. In this way, the proposed trajectory strategy is of great efficiency in smoothing the velocity parameter and guaranteeing motion safety.
During periods of and , the reduction performance of velocity jump is worse than other periods. According to (37), the compensation coefficient in this period equals 0, which means no compensation is introduced. So when joint-locked failure occurs in this period, the velocity jump cannot be reduced with the proposed method. However, the velocity jump in this period is below , which is acceptable for space manipulator in most applications.
Periods of and are almost consistent with the acceleration stage and deceleration stage of path planning. Velocity in operation space changes quickly in the two stages, making the velocity compensation insufficient for jump reduction. As a result, the reduction performance is not ideal in these periods. Period corresponds to the uniform motion stage of path planning, the velocity in operation space during this stage is stable, and the velocity jump reduction is significant. It can be concluded that velocity changing in operation space at the failure time can influence the performance of velocity jump reduction.
In conclusion, the coefficient function obtained based on a priori failure information is verified efficiently in reducing the velocity jumps caused by an unexpected joint-locked failure. And velocity changing in operation space at the failure time can influence the reduction performance. The limitation of the method is that when joint-locked failure occurs at the initial and end period of path planning, the velocity jump cannot be reduced. However, velocity jumps in these periods are always small, which can be ignored in space manipulator applications.
5. Conclusion
With the proposed trajectory optimization strategy in the paper, the on-orbit operations can be executed to completion without shutting down the space manipulator when joint-locked failure occurs unexpectedly. The velocity jumps in both joint space and operation space caused by joint-locked failure are simultaneously reduced, and the impact caused on space manipulator structure can be eliminated. As a result, the safety of space manipulator and related parts can be guaranteed. The core of the strategy is the construction of global compensation vector for joint velocity, which consists of the gradient of reduced manipulability, orthogonal matrix of null space, and coefficient. The strategy can deal with velocity jumps reduction problem caused by joint-locked failure at either a particular time or an unexpected time. For joint-locked failure at a particular time, the optimal coefficient of compensation vector can be obtained to achieve velocity jump minimal. For joint-locked failure at an unexpected time, a priori failure information is obtained by assuming the possible failure time and calculating the corresponding optimal coefficient. Then the law of coefficient is concluded for trajectory optimization. As a result, coefficient is turned into function form and a function fitting method based on a priori failure information is proposed. With simulations, the efficiency of the method is verified in reducing the velocity jumps no matter when joint-locked failure occurs.
For further study, the effect of the initial configuration of manipulator on the velocity jump reduction performance should be discussed in-depth. Besides, the velocity in operation space is planned by the trapezoidal planning with parabola transition in the paper. However, the reduction performance in the acceleration and deceleration stage is not ideal compared with the uniform motion stage, which is considered to be the limitation of the velocity planning method in operation space. By changing the planning way of velocity in operation space, the reduction performance for velocity jump is worthy to be researched. In addition, the paper concentrates on space manipulator and single joint failure. For arbitrary serial manipulator, the trajectory optimization strategy is adaptable, which can also be developed for multiple joint failure tolerance.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is supported by the National Natural Science Foundation of China (61403038), the National Key Basic Research Program of China (2013CB733000), the National Natural Science Foundation of China for the Major Equipment Development (no. 61327806), and Specialized Research Fund for the Doctoral Program of Higher Education (20120005120004).