A Novel Point-to-Point Trajectory Planning Algorithm for Industrial Robots Based on a Locally Asymmetrical Jerk Motion Profile

: Suitable trajectories with minimum execution time are essential for an industrial robot to enhance productivity in pick and place operations. A novel point-to-point trajectory planning algo-rithm (PTPA) is proposed to improve the motion efficiency of industrial robots. The jerk profile for a trajectory model is determined by five intervals and the jerk constraint. According to the kinematic constraints and two shape coefficients, a velocity threshold and three displacement thresholds are calculated for an individual joint to transfer the proposed jerk motion profile into four specific profiles. The optimal trajectory model of the joint is developed for the minimum-time and jerk-contin-uous trajectory via the performance evaluation with the input displacement and three displacement thresholds. Moreover, time-based motion synchronization for all joints is taken into account in PTPA to decrease unnecessary burdens on the actuators. The simulations illustrate that the execution time by PTPA is more efficient than that by other techniques. The experiments of a point-to-point application on a real six-axis industrial robot show that the absolute errors at the end of the motion for all joints are within 0.04°. These results prove that PTPA can be an effective point-to-point trajectory planner for industrial robots


Introduction
Industrial robots (IRs) have become one of the most important devices in the automotive, electronics, metal and other industries to replace the manual resource, due to their outstanding performances in terms of accuracy, repeatability and efficiency [1][2][3][4][5].Trajectory planning for IRs, a fundamental issue in these automated industrial activities [6], generates reference inputs for the controller to implement certain tasks [7], such as pick and place operations [8], and assembly operations [9].A pick and place operation, in essence, is a point-to-point (PTP) trajectory planning problem for a suitable and smooth trajectory from an initial position to a goal position.Otherwise, the possible excitation of mechanical resonance mode may result in noise and decrease of the lifetime of IRs [10][11].
In the literature, some research efforts have attempted to employ elementary functions for a smooth trajectory, including the polynomial, trigonometric, exponential function, etc. [12].Jond et al. employed the polynomials of degrees 3, 4 and 5 to perform trajectory motions with continuous velocity and acceleration, where the acceleration values at the beginning and end motion are non-zero [13].Boryga et al. planned the acceleration profile with the polynomials of degrees 5, 7 and 9 [14].Fang  Wang and Zhigang Jiang cyclic point-to-point tasks [15].Third-order exponential function [16] and hyperbolic tangent function [17] were also adopted to generate smooth trajectories for PTP motions.These approaches result in jerk-continuous trajectories in the case of kinematic constraints, and the jerk values at the beginning and end movement are zero [14][15][16][17].However, the trajectory based on the elementary functions causes low-efficiency motion of a robot due to no acceleration phase or constant velocity phase.
Substantial research has been focused on the composition of elementary functions to improve the motion efficiency of trajectory models.The quadrinomial polynomial was utilized to connect the initial point and intermediate point (a virtual knot), while the quintic polynomial was employed to link the intermediate point and final point [18][19][20].Similarly, Bureerat et al. adopted the quintic polynomial to connect the initial point and intermediate point (a virtual knot) and then to the final point [21].Other trajectory models based on polynomial compositions have also been utilized to generate smooth PTP trajectories, such as cubic spline [22], fifth-order B-spline [23], trapezoidal velocity profile [24] and polynomial S-curve models (i.e., motion profiles) [25][26][27][28].In addition, Perumaal and Jawahar employed a trigonometric S-curve model with a polynomial of degree 2 to generate jerk-bound trajectories with synchronized acceleration, constant velocity and deceleration phases [29].Motived by their work, Valente et al. adopted the generated trajectories with asynchronous acceleration and deceleration phases, and found they greatly extended the number of feasible solutions [30].To better utilize the kinematic constraints, one or both of the constant acceleration and jerk phases was added to the motion profiles for trajectory models [31][32][33].Among these models, the trajectory model based on a sinusoidal jerk motion profile with fifteen segments has shown its superiority in motion efficiency for pick and place operations [31].Nevertheless, the trajectory resulting from the jerk profile with a high ramp coefficient leads to the low-efficiency motion of a robot due to the similar ramp-up to the ramp-down in the jerk profile.Some researchers also paid abundant attention to pursuing optimal solutions based on various objectives, such as minimum-time [14][15][16][30][31][32][33], minimum-jerk [21][22][23], minimum-energy [34][35], etc.The numerical and analysis methods are two optimal solutions for the solver type.Valente et al. adopted an active-set algorithm to globally minimize the total motion time [30].Chettibi et al. employed the sequential quadratic programming (SQP) method to obtain the minimum cost with transfer time and the quadratic average of both actuator efforts and corresponding power [22].Intelligent algorithms, such as genetic algorithm (GA) [13,18,19], artificial bee colony (ABC) algorithm [20] and differential evolution (DE) algorithm [21] can solve the multi-objective optimization problems for trajectory planning.The numerical method may fail to obtain a solution or become trapped in local convergence, due to the unsuitable selection of an initial solution and high computational complexity.Then, some studies are focused on the closed-form solutions for the global minimum and high computational efficiency by using the kinematic constraints.The properties of root multiplicity were employed to generate the high-order polynomial trajectory model by only one polynomial coefficient which can be computed from the kinematic limits [14,27].The execution time can be minimized from the kinematic limits and input displacements [15,16,25,[31][32][33].
The goal of this work is to provide a closed-form solution for the time-optimal and smooth trajectory planning of IRs based on an improved jerk motion profile.Since the motion with high efficiency for a robot can be produced by the kinematic peak values of a trajectory model, the locally asymmetrical jerk motion profile with low jerk ramp-up is proposed for the trajectory model to obtain rapidly the jerk peak value.The ramp coefficient and local asymmetry coefficient are adopted to reduce the number of determined intervals for the jerk profile.Furthermore, the relationship between the execution time and the two shape coefficients is investigated quantitatively to improve the motion efficiency with lower jerk ramp-up.According to the two coefficients and the kinematic constraints, four thresholds are calculated for an individual joint to transfer the proposed jerk motion profile into four specific profiles.Then, the optimal trajectory model of the joint is ascertained for the minimum-time and jerk-continuous trajectory via the performance evaluation with the input displacement and three displacement thresholds.Moreover, time-based motion synchronization for all joints is also taken into account in the proposed algorithm to decrease unnecessary burdens on the actuators.
The structure of this paper is organized as follows.Section 2 introduces the proposed trajectory model based on a locally asymmetrical sine jerk profile and spreads out the corresponding mathematical formulas.This section also presents the impact of two shape coefficients on the execution time.Section 3 presents a formal definition of the time-optimal trajectory planning problem.PTPA is proposed in Section 4. Section 5 presents the simulation results compared with other methods and experimental results on a real industrial robot.Concluding remarks and future extensions are given in Section 6.

Locally Asymmetrical Jerk Motion Profile
The conventional trajectory model with constant acceleration and velocity phases has similar ramp-up and ramp-down segments in the jerk motion profile (red line in Figure 1), which may reduce the motion efficiency in the case of a high ramp coefficient [31].Thus, to overcome this deficiency, we propose the locally asymmetrical jerk motion profile with low ramp-up, as shown (blue line in Figure 1).We define the terms similarly in [36] to describe the trajectory segments.D f denotes the displacement at the end of the sustained pulse of the velocity.V p , A p and J p represent the peak values of the velocity, acceleration and jerk during the motion, respectively.ti represents the time at the end of the ith segment, where i ∈ {0, 1, …, 15}.The trajectory model comprises three phases: the acceleration phase (AP) when t∈ [t0, t7], the constant velocity phase (CVP) when t ∈ [t7, t8], and the deceleration phase (DP) when t ∈ [t8, t15].Let Ti denote the interval of the ith segment, i.e., Ti = ti-ti-1, where i ∈ {1, …, 15}.We assume that T1 = T7 = T9 = T15, T2 = T6 = T10 = T14, T3 = T5 = T11 = T13 and T4 = T12.Then, as shown in Figure 1, the ramp-up and ramp-down of the velocity and acceleration are all symmetrical all the time.Note that a locally asymmetrical jerk motion profile exists in the case of T1 ≠ T3.Then, the execution time Tt can be expressed as where T1, T2 and T3 are the intervals of the ramp-up, cruise and ramp-down of the jerk, respectively; T4 and T8 are the intervals of the acceleration cruise and the velocity cruise, respectively.The locally asymmetrical jerk motion profile in Figure 1 is defined as Let Ai, Vi and Di denote the motion profiles of the acceleration, velocity and displacement for the ith segment, respectively.Given the initial conditions of the ith segment, Ai, Vi and Di can be computed by sequentially integrating their corresponding jerk equations, and are given by where Ai−1 e , Vi−1 e and Di−1 e (i = 1, 2, …, 15) denote the values of the acceleration, velocity and displacement at the end of the (i−1)th segment, respectively.For example, A0 e , V0 e and D0 e , represent the initial conditions of the acceleration, velocity and displacement of the first segment, respectively.The symmetry of the proposed jerk profile can simplify the integration process.The detailed equations of each segment are expressed in the Appendix A. Here, several important equations are utilized to compute the determined intervals.
A p and V3 e are the acceleration and velocity reached at the end of the first acceleration ramp-up from zero acceleration to A p , respectively, and are given by V p and D7 e are the velocity and displacement reached at the end of the velocity rampup from zero speed to V p , respectively, and are described as

T T T T T T T T T T T T D J T T T T T T T T T T T T T T T T T
D f is the displacement reached at the end of the movement and is given by If J p is specified, the proposed trajectory model can be fully determined by the five intervals T1, T2, T3, T4 and T8.However, the number of the determined intervals is still too many to form a closed-form solution for a time-optimal trajectory.

Shape Coefficients
To decrease the number of the determined intervals and enhance the computational efficiency, the relationships between T1 and T2, T1 and T3 are established, respectively.T4 and T8 are both set to zero in this subsection, because they may not exist in some circumstances, e.g., a case with large kinematic limits and low displacement.Moreover, to avoid the impact of the constraints on the execution time, we assume that the kinematics limits cannot be reached.Then, two shape coefficients are adopted to describe the shape of the jerk profile.Let α and β denote the ramp coefficient and the local asymmetry coefficient, respectively, which are expressed as Nonzero jerk can occur at the initial and goal positions during the motion when α or β is set to 0 (i.e., T1 = 0), which will increase the wear of the robot structure and should be avoided to enhance the reliability of IRs.Thus, α > 0 and β > 0, i.e., T1 > 0. In addition, to generate a jerk-continuous trajectory, β ≠ 1 (i.e., T3 ≠ 0).Note that the local jerk profile is symmetrical when β = 0.5 (i.e., T1 = T3).Therefore, the domain of α is set to (0, 1] and β to (0, 1), to ensure that the jerk of the generated trajectory is continuous and its values are zero at the initial and goal positions.The jerk motion profiles with different α and β are shown in Figure 2. As T4 = 0 and T8 = 0, substituting Equations ( 8)- (10) into Equation ( 1), Tt is related to α and β by α β αβ αβ where Thus, D f and J p can be considered as the proportionality coefficients of Tt, and the variation trend of Tt is only determined by α and β.For example, when D f = 140 deg and J p = 70 rad/s 3 , the impact of different α and β on Tt is depicted in Figure 3.In Figure 3, Tt declines significantly and steadily with the decreasing β when α ≥ 0.5.The reason is that the jerk can rapidly reach the peak as T1 decreases with smaller β.However, the impact is different when α < 0.5.There exists a threshold with the smallest execution time Tt.Let βth denote the threshold.βth increases with the decreasing α, and approaches 0.5 when α is fairly close to 0, because the jerk profile considerably approaches the bang-bang jerk profile.Therefore, for a given α, there is a corresponding β with the smallest execution time.

Optimal Problem Formulation
The optimal 14 of PTP trajectory planning for IRs is generally presented in the Cartesian space or the joint space.Planning trajectory in the joint space is difficult to predict the motion of the end effector of IRs, but it can avoid kinematic singularities and reduce the computational complexity, compared with that in the Cartesian space [37].Therefore, the goal of this paper is to obtain time-optimal trajectories in the joint space.
Let ℕ>0 be the set of positive integers and ℝ the set of reals.Let ℝ≥0, ℝ>0 and ℝ<0 represent the sets of non-negative, positive and negative reals, respectively.Then, let n ∈ ℕ>0 be the number of the joints of an industrial robot.Clearly, |V p |, A p , J p , Tt ∈ ℝ>0 n .Let m ∈ ℕ>0 denote the number of the trajectory segments.Then, let Γ ∈ ℝ≥0 m×n denote the matrix of the interval, i.e., Γ = {Ti,g}, where i, g ∈ ℕ>0, i ≤ m and g ≤ n.The 1-norm is denoted by ||||1.Note that ||Γ||1 = max {Tt}.Thus, the objective function in this paper can be described as Let V max , A max , J max ∈ ℝ>0 n denote the kinematic limits of the velocity, acceleration and jerk of an industrial robot, respectively.We assume that V min = −V max , A min = −A max and J min = −J max , where V min , A min , J min ∈ ℝ<0 n are the constraint minimum of the velocity, acceleration and jerk, respectively.Then, the kinematic constraints of the optimality problem are given by Let θinitial, θgoal ∈ ℝ n denote the initial and goal positions, respectively.Given a, b ∈ ℝ, a closed interval between a and b is denoted by [a, b].Consequently, the optimality problem of the PTP trajectory planning in this paper is to find feasible trajectories with minimal objective function while respecting the kinematic constraints, formalized as follows.
Given an initial position θinitial, a goal position θgoal, and the kinematic limits of the velocity V max , acceleration A max and jerk J max , find trajectories Ω:

Point-to-Point Trajectory Planning Algorithm
A PTPA is proposed to provide a closed-form solution for the time-optimal PTP trajectory.PTPA is divided into four procedures: data input, initialization, minimum-time calculation and motion synchronization.Pseudo-code of PTPA is shown in Algorithm 1. Compute the thresholds V amax , D amax , D vmax1 , D vmax2 using Equations ( 18), ( 19), ( 23) and ( 26) Compute T1,g, T2,g, T3,g, T4, g and T8,g according to Type I Compute T1,g, T2,g, T3,g and T4, g according to Type II T4, g = 0 Compute T1,g, T2,g, T3,g and T8,g according to Type III else T4, g = 0, T8,g = 0 Compute T1,g, T2,g and T3,g according to Type IV end if end for Generate the execution time τ using Equation (28) Compute the scaling coefficient k using Equation (29) Synchronize the motion profiles of each joint trajectory end

Minimum Time
According to T4 and T8, the jerk motion profiles are classified into four specific ones to improve the computational efficiency, as shown in Figure 4.For a given PTP task of an industrial robot, the type of jerk motion profile of a joint is mainly determined by a velocity threshold, three displacement thresholds and an input displacement.Let Vg amax ∈ ℝ be the velocity threshold, and Dg amax , Dg vmax1 , Dg vmax2 ∈ ℝ the three displacement thresholds.Note that the subscript g refers to the gth joint.Vg amax : the velocity reached at the end of the acceleration pulse from zero acceleration to Ag max and back down to zero acceleration, without taking the velocity limit into account.
Dg amax : the displacement reached at the end of the velocity pulse from zero velocity to Vg amax and back down to zero velocity (i.e., Vg p = Vg amax ).
Dg vmax1 : the displacement reached at the end of the velocity pulse from zero velocity to Vg max and back down to zero velocity (i.e., Vg p = Vg max ), while Vg amax is smaller than Vg max (i.e., Vg amax < Vg max = Vg p ).
Dg vmax2 : the displacement reached at the end of the velocity pulse from zero velocity to Vg max and back down zero velocity (i.e., Vg p = Vg max ), while Vg max is smaller than or equal to Vg amax (i.e., Vg p = Vg max ≤ Vg amax ).
The shortest execution time can be achieved by utilizing the kinematic limits [14-16,31-33].However, the velocity and acceleration limits may be unavailable in some tasks, such as a small displacement with large kinematic limits.Therefore, we assume that J p = J max in this part.Note that only the definitions of T1,g, T4,g and T8,g are given, since T2,g and T3,g can be obtained from T1,g by using Equations ( 9)-( 10) if α and β are specified.
T4,g >0 implies that the acceleration limit Ag max must be reached at the end of the first acceleration ramp-up (i.e., Ag p = Ag max ), and Vg amax must be smaller than Vg p (i.e., Vg p > Vg amax ).Thus, T1,g is obtained by substituting Equations ( 9) and (10) into Equation (4): T8,g > 0 represents that the velocity limit Vg max must be reached at the end of the velocity ramp-up (i.e., Vg p = Vg max ).Then, T4,g is obtained by using Equation ( 6): Furthermore, the displacement Dg f must be greater than Dg vmax1 (i.e., Dg f > Dg vmax1 ).Then, T8,g is obtained by using Equation ( 8): Consequently, this type can be selected when Vg max > Vg amax and Dg f > Dg vmax1 .Then, Vg amax can be obtained by using Equation ( 6) when T4,g = 0: Moreover, Dg vmax1 can be obtained by using Equation ( 8) when T8,g = 0: Type II: T4,g > 0 , T8,g = 0.
T4,g = 0 implies that Vg p must be equal to or smaller than Vg amax (i.e., Vg p ≤ Vg amax ).T8,g > 0 represents that Vg p = Vg max as mentioned previously.Additionally, the displacement Dg f must be greater than Dg vmax2 (i.e., Dg f > Dg vmax2 ).Consequently, this type can be selected when Vg max ≤ Vg amax and Dg f > Dg vmax2 .Note that T8,g can be computed from Equation (17), if T1,g is given.
Except for the three types above, T1,g can be obtained by substituting Equations ( 6), ( 9) and ( 10) into Equation ( 7): Four thresholds can be obtained, when α and β are specified for an industrial robot.Then, the jerk motion profile for the time-optimal trajectory model can be transferred into a specific one by comparing the input displacement and three displacement thresholds in a given PTP task.Thus, we can obtain the matrix of the interval Γ and the execution time Tt generated based on the specific profile.

Synchronization of Motion Profiles
Since the task and the kinematic limits are different for each joint of an industrial robot, their execution time varies in lengths.This may result in unnecessary burdens on the actuators without the synchronization of all joints.Time-based synchronization is taken into account to decrease the burden and ensure that all joints can start and stop at the designated position simultaneously.Due to the utilization of the kinematic limits, Tt can be considered as the minimum set of the execution time for each joint, such that Let k ∈ ℝ>0 n denote the set of the synchronization coefficients, i.e., k = {kg}, where g ∈ ℕ>0, g ≤ n, Let t ⸦ ℝ≥0 denote the time variable, and Dg(kgt), Vg(kgt), Ag(kgt), Jg(kgt) ⸦ ℝ n the trajectory and its derivatives.Then, the following holds: For the joints with kg > 1 in a given PTP task, the peak values of the velocity, acceleration and jerk decrease after synchronization, which can protect the drivers in the joints to extend their service life.
Herein, we summarized the procedures of PTPA.The initial and goal positions of each joint with their kinematic limits are defined in advance, as well as two shape coefficients.According to the kinematic constraints and the shape coefficients, the displacements and the thresholds are computed as the selection criterion for the jerk profile.A specific jerk profile of the optimal trajectory model for a joint is ascertained for the minimum-time and jerk-continuous trajectory via the performance evaluation with the input displacement and three displacement thresholds.Finally, the time-based motion synchronization for all joints is carried out to decrease unnecessary burdens on the actuators.

Results
To demonstrate the feasibility of the proposed algorithm, the numerical simulation and experiments tested on a real industrial robot were carried out.We ran the numerical simulations using the MATLAB R2018b on a Windows 10 personal computer with a Core i5-1035G1 CPU and 16 GB of RAM.An experimental platform was built up to mimic a pick and place operation, and validate the numerical simulations.This platform mainly included a six-axis industrial robot Universal Robots UR5, a robot controller CB3 and a personal computer with Linux operating system, as shown in Figure 5.The robot controller was controlled over TCP/IP protocol by a script implemented in Python using the Real-Time Data Exchange (RTDE) interface ur-rtde 1.4.3.The script fed the joint position targets to the robot controller to simulate a pick and place operation.This script was implemented on a personal computer with a Core-i3 4130 processor and 4 GB of RAM.

Numerical Simulation
A typical PTP task for a six-axis industrial robot was chosen to compare the simulated execution time by five different types of trajectories [22,[29][30][31][32].The PTP task to be tested is listed in Table 1, where the initial and goal positions of each joint are presented, as well as the kinematic limits of the velocity, acceleration and jerk.α was set to 0.1, 0.5 and 1 [31].β was set to 0.3, 0.1 and 0.1 according to the Equations ( 11) and (12).The motion profiles of joint 4 which determined the execution time, generated by the benchmark methods and this work, are shown in Figure 6.The cubic polynomial trajectory model failed to produce jerk-continuous motion, and the values of the acceleration and jerk were both non-zero at the initial and goal positions, which can lead to the vibration and the wear of the joint.The motion profiles generated by [22,29] show that the jerk limits failed to be reached, which resulted in the longer execution time compared with other algorithms.In contrast, the jerk limit was utilized from Refs.[30][31][32], which resulted in higher efficiency trajectories compared with Refs.[22,29].The proposed algorithm yielded the best velocity in the shortest time, but it failed to obtain better acceleration than the benchmark methods.Table 2 shows the measured performances including the execution time and maximum jerk by the benchmark algorithms and the proposed one for the test task.The results indicate that the algorithm presented in this paper yielded the smallest execution time compared with the benchmark ones.The obtained execution time from the trajectory model with locally asymmetrical jerk motion profiles was smaller than that with symmetrical jerk motion profiles [31] in the case of the appropriate value of β.For α = 0.5, the resulting execution time was decreased by 2.1% compared with the symmetrical 15-phase sine jerk trajectory model [31], and decreased by 40% compared with the cubic polynomial trajectory model [22].The max jerk appeared simultaneously at joint 1 and joint 4 in Refs.[29,31] and this work, although the constraints of these two joints were different, because the execution time by these three trajectory models [22,29,31] was mainly determined by the input displacements and the jerk constraints.The execution time by the proposed algorithm with α = 0.01 and β = 0.5 was close to that achieved by Ref. [32] with Smax = 4000rad/s 4 , because the bang-bang jerk profiles were considerably approached by these two trajectory models in these cases.

Work
Trajectory model Performance Measure Ex.Time (s) Max. jerk (rad/s3 ) Chettibi et al. [22] Cubic polynomial 2.675 30.56 (joint 2) Perumaal et al. [29] 3-segment sine jerk (sync acc/dec) 2.5133 6.63 (joint 1, 4) Valente et al. [30] 3-segment sine jerk (async acc/dec) 1.7395 20 (joint 4) Fang et al. [32] 15-segment sigmoid jerk 1.8760 (Smax=150rad/s 4 ) 20.3 (joint 1) Fang et al. [31] Symmetrical 15-segment sine jerk 1.5309 (α=0.1) The superiority in efficiency will be more significant, when the kinematic constraints with the same acceleration limit have a larger velocity limit or a smaller jerk limit.This feature enables the proper management of the degrading patterns for the specific joints.For instance, the contrastive task was calculated by this work and Ref. [30], in which the displacement and constraints of the joint 4 were adjusted.The results in Table 3 indicate that the percentage reduction (PR) of the execution time increased from 8.05% to 9.63% due to the longer CVP.Table 4 indicates that the proposed method had better performance on reducing the execution time, due to the jerk cruise in the profile and the better utilization of the jerk limits.When α = 0.5 and β = 0.1, the motion profiles of the joint position, velocity, acceleration and jerk generated by the proposed algorithm for the test task are illustrated in Figure 7.The optimized trajectories were smooth, and the motion curves under the kinematic constraints were all bounded and continuous all the time.The values of velocity, acceleration and jerk were zero at the initial and goal positions.In addition, all the joints started and stopped at the designated positions synchronously.

Experimental Results
To satisfy the operating conditions of the platform and for the industrial robot to avoid collisions with the workbench, the test task needed to be adjusted.The initial and goal positions of joint 2 and joint 3 were altered while maintaining their displacements.Specifically, the initial and goal positions of joint 2 were -3π/4rad and -5π/12rad, respectively.Similarly, the initial and goal positions of joint 3 were π/4rad and π/2rad, respectively.Moreover, the velocity limits of all joints were 3rad/s within the maximum speed range.Note that α and β were 0.5 and 0.1, respectively.
Twenty trials were carried out to compare the results of the experiments and simulations for reducing the random errors, and to obtain the measurement uncertainty expressed in terms of the variance.The results of a trial with the corresponding tracking errors in the joint space are depicted in Figure 8.The experimental curves were consistent with the simulated ones in general, except for some tiny deviations resulting from the parameters configured in PID.Additionally, several abnormal data, circled in Figure 8, appeared randomly, which may have been caused by the interference in the data transmission from the sensors to PC. Ignoring these abnormal data, the mean values of maximum absolute tracking errors and terminal errors with their variances of all joints for twenty trials are listed in Table 5.The results showed that the maximum absolute tracking error of a joint was proportional to the assigned displacement of the joint.Moreover, the absolute terminal errors at the goal positions were all within 0.04° and their variances were close to zero.Therefore, the experiments on a real industrial robot verify the effectiveness of PTPA.

Conclusions
A novel point-to-point trajectory planning algorithm (PTPA) based on a locally asymmetrical jerk motion profile is proposed for the time-optimal and smooth joint trajectories of IRs.The ramp coefficient and the local asymmetry coefficient determine the shape of a jerk motion profile, and affect the execution time.The results show that the locally asymmetrical jerk motion profile with lower ramp-up can generate a higher efficiency motion.
The execution time by PTPA was the smallest compared with the existing techniques, while the jerk, acceleration and velocity of the generated trajectories were all bounded and continuous, with the zero values at the initial and goal positions.The experimental results of angular profiles on a real industrial robot were consistent with the simulated results.Both the simulations and experiments prove that the PTPA can be enacted as an effective tool for the trajectories planning of IRs.
Future work will be devoted to applying the present method for multi-point trajectory planning to broaden real-time performance.Furthermore, the collision avoidance strategy in the operation space with a complex working environment should be investigated.Finally, energy optimization and dynamic constraints are also worth exploring in the proposed algorithm.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
et al. utilized the trigonometric function to build up a frequency central pattern generator (CPG) trajectory model for Citation: Wu, Z.; Chen, J.; Bao, T.; Wang, J.; Zhang, L.; Xu, F. A Novel Point-to-Point Trajectory Planning Algorithm for Industrial Robots Based on a Locally Asymmetrical Jerk Motion Profile.Processes 2022, 10, 728.https://doi.org/10.3390/pr10040728 Academic Editors: Wei Cai, Yan

Figure 1 .
Figure 1.Jerk, acceleration, velocity and displacement profiles of the proposed trajectory model.

Figure 3 .
Figure 3. Impact on Tt with different values of α and β.

Figure 4 .
Figure 4. Four types of the locally asymmetrical jerk motion profile.

Figure 6 .
Figure 6.Motion profiles of joint 4 for different trajectory models.

Figure 8 .
Figure 8. Measured joint trajectories for the adjusted task.

Author Contributions:
Conceptualization, Z.W. and J.C.; methodology, Z.W.; writing-original draft preparation, Z.W., J.C. and T.B.; writing-review and editing, J.C. and J.W.; supervision, L.Z. and F.X.All authors have read and agreed to the published version of the manuscript.Funding: This research was funded by [National Key Research and Development Program of China] grant number [2018YFB1309404], [Research on Public Welfare Technology Application Projects of Zhejiang Province] grant number [LGG18E050023] and [General Scientific Research Project of Education Department of Zhejiang Province] grant number [Y202147343].Data Availability Statement: Not applicable.Acknowledgments: The authors acknowledge the support and inspiration of the national key research and development program of China, research on public welfare technology application projects of Zhejiang province and the general scientific research project of the education department of Zhejiang province.

Table 1 .
Desired positions and kinematic limits of the test task.

Table 2 .
Measured performances from various trajectory models for the test task.

Table 3 .
The elements of the execution time for two tasks.

Table 4 .
Measured performances in different degradation modes.

Table 5 .
Mean values of maximum absolute tracking errors and terminal errors with their variances of all joints.