Trajectory Optimization Algorithm for a 4-DOF Redundant Parallel Robot Based on 12-Phase Sine Jerk Motion Proﬁle

: To improve high motion accuracy and efﬁciency in the high-speed operation of a 4-DOF (4 degrees of freedom) redundant parallel robot, this paper introduces a trajectory planning of the parallel robot in joint space based on the twelve-phase sine jerk motion proﬁle. The 12-phase sine jerk motion proﬁle utilizes the characteristics of a sine function. Furthermore, the penalty function is used to optimize the trajectory energy consumption under the constraint condition. The simulation and experimental results show that the energy consumption of joint space is slightly higher than that of the three-phase sine jerk motion proﬁle, but the overall operation is more accurate and stable. Speciﬁcally, the sudden change of force and velocity in each joint is eliminated, which is the cause of mechanism oscillation. Moreover, the force of each joint is more average. The results indicate that each movement is closer to the maximum allowable limit and the running efﬁciency is higher.


Introduction
The 4-DOF redundant parallel robot is widely used in industrial production lines. It has the advantages of good rigidity, high movement precision, and compact structure [1][2][3][4][5][6]. Thus, it is especially suitable for the operation of lightweight objects, sorting and packing in high-speed handling, etc. [7,8]. The inertia of the robot arm, however, causes vibration and jerk when the robot moves at high speed, leading to the low accuracy and stability of the robot. Scholars have carried out in-depth research from the point of trajectory planning. For example, in Reference [9], the three-phase sine jerk motion profile was used as trajectory planning to reduce the jerk caused by the high-speed point-to-point motion of the robot. Fang et al. [10] proposed a method of expanding the three-phase sine jerk motion profile into a fifteen-phase sine jerk motion profile. Based on the research of Li and Fang, with the limit of each joint of the robot motion, Valente et al. [11] proposed a multivariable time optimization method to plan the optimal trajectory of joint space on the three-phase sine jerk motion profile and made the running time of the trajectory shortest. To reduce the computation load of trajectory planning for online calculation of trajectory planning, some researchers adopt the method of piecewise interpolation. Wang et al. [12], based on the three-phase sine jerk motion profile, constructed a new compound sine trajectory to realize the online fast interpolation feature of joint space motion trajectory, through using polynomial transition interpolation at the beginning and end of joint motion.
However, with the acceleration of the speed of the robot, the energy loss in the running process has become a non-ignorable problem. Therefore, some scholars began to discuss the possibility of robot trajectory optimization from the point of energy and to obtain many valuable research results on the optimization algorithm of robot energy   Figure 1b, the base coordinate is xoy, and the coordinates of the joints are x i o i y i (i = 1, 2, 3). Three plinths A1 (Xa1, Ya1), A2 (Xa2, Ya2), and A3 (Xa3, Ya3) are used to connect the 3 driving manipulators. The driving and driven angle is q 1 (θ 1 , α 1 , β 1 ) and q 2 (θ 2 , α 2 , β 2 ), respectively. The length of the Actuators 2021, 10, 80 4 of 18 3 driving arms and the 3 driven arms is L. The end-effector is C(X, Y). The distance Aic from the base to the end-effector is as follows: Using geometry and the Law of Cosines, the driving angle q 1 can be denoted as follows: /(X C −xai)) + arccos(Aic/(2L))]180/π (2) The follower angle q 2 is as follows: where the q, q 1 , and q 2 represent each driving joint variable, driven joint variable, and each joint variable, respectively. The coordinate vector of the end-effector is q C = X C Y C T .
From the kinematic equation, the velocity of the end-effector output terminal point C can be expressed as follows: . .
In Equation (5), J is the velocity Jacobian Matrix, which maps the joint space onto the operating space. The generalized velocity of the end-effector is Acceleration is a(t) = ..
we have the following: From Equations (1)-(8), we can see that the velocity, acceleration, and acceleration rate in the joint space and the operating space are all sine functions of time. They are obtained by integrating the acceleration rate curve in joint space several times. These sine functions can reduce the calculation and improve the stability of the system. The trajectory optimization equation of system joint space is constructed according to the motion requirement of joint space, the indices of time, and the stability of system operation.

Construction Principle of 12-Phase Sine Jerk Motion Profile
During the process from initial status to final status, the robot component has to go through several stages, such as acceleration, constant speed, and deceleration, which form a typical jerk curve. To prevent the jerk caused by acceleration and deceleration of the mechanism, the sine curve is used to describe the acceleration rate in the acceleration and deceleration phase in Reference [7]. The three-phase sine jerk motion profile, a previous method to the 12-phase sine jerk motion profile, is shown in Figure 2.
In Figure 2, jerk(t), the function of the jerk curve to time, can be expressed as follows: acceleration section, constant speed section, deceleration section.
In Figure 2, ( ), the function of the jerk curve to time, can be expressed as follows: According to Figure 2 and the jerk(t) function, the three-phase sine jerk motion profile, though the machine can move smoothly, fails to move with the highest speed or efficiency, causing each part cannot reach the movement limit of all key quantities, which can result in low running efficiency. Therefore, for reducing the running time of the component, the method, adding the middle constant velocity section, is adopted to make the component reach the peak value quickly and then gradually reduce to zero. Meanwhile, to ensure the component is in the motion saturated state as long as possible and to reduce the jerk frequency, the acceleration rate of the three-phase sine jerk motion profile is revised to the 12-phase sine jerk motion profile, as shown in Figure 3.
According to Figure 2 and the ( ) function, the three-phase sine jerk motion profile, though the machine can move smoothly, fails to move with the highest speed or efficiency, causing each part cannot reach the movement limit of all key quantities, which can result in low running efficiency. Therefore, for reducing the running time of the component, the method, adding the middle constant velocity section, is adopted to make the component reach the peak value quickly and then gradually reduce to zero. Meanwhile, to ensure the component is in the motion saturated state as long as possible and to reduce the jerk frequency, the acceleration rate of the three-phase sine jerk motion profile is revised to the 12-phase sine jerk motion profile, as shown in Figure 3.  As is shown in Figure 3a, the angular acceleration rate curve ... q (t) of the 12-phase sine jerk motion profile is a piecewise function of constant j peak motion at each positive Actuators 2021, 10, 80 7 of 18 and negative peak of a continuous sine curve. The analytical expression can be described as follows: ..
Furthermore, the angular acceleration function .. q(t), the angular velocity curve function . q(t), and the angular displacement curve function q(t) are derived as follows, respectively: .. . .
According to the time interval shown in Figure 3, the piecewise function of the angular velocity . q(t) of the 12-phase sine jerk motion profile function in each time interval can be expressed as Table 1.
Note: ... q (t) is triple-integrated to obtain the angular displacement q(t) curve. Similarly, the angular acceleration function and the angular displacement function can be treated by Equations (9) and (11).

Total Mechanical Energy Consumption Model
According to the interval segment of the 12-phase sine jerk motion trajectory, considering the time of the robot in each interval, the maximum torque, and rotational speed, the energy consumption of the mechanism will be reflected. Therefore, to minimize energy consumption, the 12-phase sine jerk motion profile is optimized, and the interval time is determined accordingly, which provides the basis for the follow-up control of the system.

Establishment of Total Mechanical Energy Consumption Model
According to Figure 1, the potential energy of a 4-DOF redundant parallel robot is zero due to its planar motion, and the effect of friction in each joint is ignored. It is supposed that the energy of the system is converted to kinetic energy completely under ideal conditions.
The angular velocity of the driving arm is . β 2 (t) . The mechanical energy consumption E(t) of the whole robot is as follows: In Equation (12), K i1 is the mechanical energy consumption of the driving arm Aibi. K i2 is the mechanical energy consumption of the driven arm BiC (i = 1, 2, 3). K i1 and K i2 correspond to the kinetic energies in the parallel robot.
Thus, the mechanical energy consumption of the drive shaft is: The mechanical energy consumption of the driven shaft is: In Equation (13), p i1 is the distance from the centroid of the driving shaft to the center of the base. In Equation (14), p i2 is the distance from the center of mass of the driven shaft to the center of the joint, as shown in Figure 1b.
According to Equation (3), q 2 (t) is a function of q 1 (t). q 1 (t) is obtained by integration of . q 1 (t). The Equation (12) is substituted into the Equations (13) and (14), to obtain the total mechanical energy consumption function, E(t), in joint space.

Time-Energy Optimal Solution
In order to find the minimum value of E(t), the total mechanical energy consumption function, and to determine the time of each interval, E(t) is taken as the objective function. The energy and the time of each interval are solved by subdivisional calculation method under the constraint of kinematics and dynamics equations, for achieving a comprehensive optimization. Therefore, according to Equation (10), the angular velocity expression . q(t) of the robot driving arm is obtained. The rotational speed, n, of the driving arm input is found to be as follows: According to the principle of electromechanics0, the motor maximum output power is P max , and the minimal speed is n min . The transmission ratio between the motor is i, and the maximum working efficiency of the driving arm is η. Then, the maximum input torque, M max , of the driving arm can be obtained: Finally, the torque relationship between the driving torque and the joints on the arm is obtained by utilizing the Lagrange dynamic Equation (17) as follows: In Equation (17), M is the joint input torque: According to the objective function E(t), taking rapidity and stability as the constraint conditions, the solution model is established as follows: Objective function: min E(t) Constraint condition: In the equation, q 0 is the initial angle of the driving joint; q 12 is the endpoint angle of the driving joint. The variable is X = T 1 T 2 T . Moreover, g u (X)(u = 1, 2, · · · , 6) is an inequality constraint, where g 1 (X), g 2 (X), and g 3 (X) denote that the values of T 1 , T 2 , and j peak are greater than zero. Moreover, g 4 (X) is the torque constraint, g 5 (X) is the angular velocity constraint, and g 6 (X) is the angular acceleration constraint. An equality constraint h z (X)(Z = 1, 2) is established from the starting and ending positions of the locus, where h 1 (X) represents the angular displacement starting angle of the 12-phase sine jerk motion profile, and h 2 (X) represents the angular displacement ending angle of the 12-phase sine jerk motion profile. The optimization variables is T 1 , T 2 , and j peak . The optimization algorithm penalty function is used to optimize multiple variables of the objective function under constraints. According to the objective function, equality constraint condition, and inequality constraint condition, the penalty function φT 1 , T 2 , r (k) is constructed as follows: In Equation (18), r (k) is the penalty factor, a decreasing sequence of positive numbers. According to Equation (18), the parameters T 1 , T 2 , and j peak in the motion trajectory are calculated, and the penalty function φT 1 , T 2 , r (k) is minimized to obtain the optimal trajectory value.

Simulation and Experiment Results Analysis
To verify the correctness of the above theoretical analysis, Matlab was used in the simulation. Servo studio, GPM II, Mct 2008, Matlab, and other software were used in the experiment. Taking the 4-DOF redundant parallel robot shown in Figure 4a as the foundation of experiment and simulation analysis, the joints and coordinate positions are defined as shown in Figure 4b.

Simulation and Experiment Results Analysis
To verify the correctness of the above theoretical analysis, Matlab was used in the simulation. Servo studio, GPM II, Mct 2008, Matlab, and other software were used in the experiment. Taking the 4-DOF redundant parallel robot shown in Figure 4a as the foundation of experiment and simulation analysis, the joints and coordinate positions are defined as shown in Figure 4b.  The structure and dynamic parameters of the robot in Figure 4b are shown in Table 2. In Figure 4a, the 4-DOF parallel robot adopts Tamogawa AC servo motor TS4603 N7185 E200 with rated speed = 3000 / , max speed = 5000 / , rated power = 400 , and reduction ratio i = 80. The number of pulses in the encoder per turn The structure and dynamic parameters of the robot in Figure 4b are shown in Table 2. In Figure 4a, the 4-DOF parallel robot adopts Tamogawa AC servo motor TS4603 N7185 E200 with rated speed n e = 3000 r/min, max speed n max = 5000 r/min, rated power P e = 400 w, and reduction ratio i = 80. The number of pulses in the encoder per turn is 2 17 . The motion parameters, maximum allowable limits, of each joint in the joint space are shown in Table 3. Then, the analysis object, the linear motion of the end effector from point A to point B in the operating space in Figures 1 and 4b, is calculated by the mathematical software, using the 12-phase sine jerk motion profile function and penalty function equation. The obtained characteristic parameters and time segmentation values of the joint space are taken as the basis for the optimal trajectory simulation and experimental verification of the specific process, as follows.

Simulation Analysis of Stability and Accuracy of 12-Phase Sine Jerk Motion Profile
According to Table 1, an optimization algorithm is used to calculate the time T 1 = T 2 = 0.1 s, and the total motion time is T = 8T 1 + 4T 2 = 1.2 s, j peak = 346.261 deg/s 3 . Figure 4, the three plinths of the parallel mechanism form an equilateral triangle, and the end-effector of the robot moves from point A to point B, just on the midline of the equilateral triangle. Thus, the kinematics curves of Joint 1 are equal to Joint 2, with Joint 3 equaling to Joint 6, and Joint 4 equaling to Joint 5. Based on the obtained curves of angular displacement, angular velocity, angular acceleration, and angular acceleration rate of a joint, the motion profile of each joint can be calculated by the relationship between the joint and other joints. Setting the rotation of the anticlockwise direction as the positive direction, the Joints 4 and 5 are in a positive direction, while Joints 1, 2, 3, and 6 are in a negative direction. The angular displacement, angular velocity, angular acceleration, and rate of change of angular acceleration of each joint are shown in Figure 5.

Combined with
taken as the basis for the optimal trajectory simulation and experimental verification of the specific process, as follows.

Simulation Analysis of Stability and Accuracy of 12-Phase Sine Jerk Motion Profile
According to Table 1, an optimization algorithm is used to calculate the time 1 = 2 = 0.1 , and the total motion time is = 8 1 + 4 2 = 1.2 , = 346.261 / 3 . Combined with Figure 4, the three plinths of the parallel mechanism form an equilateral triangle, and the end-effector of the robot moves from point A to point B, just on the midline of the equilateral triangle. Thus, the kinematics curves of Joint 1 are equal to Joint 2, with Joint 3 equaling to Joint 6, and Joint 4 equaling to Joint 5. Based on the obtained curves of angular displacement, angular velocity, angular acceleration, and angular acceleration rate of a joint, the motion profile of each joint can be calculated by the relationship between the joint and other joints. Setting the rotation of the anticlockwise direction as the positive direction, the Joints 4 and 5 are in a positive direction, while Joints 1, 2, 3, and 6 are in a negative direction. The angular displacement, angular velocity, angular acceleration, and rate of change of angular acceleration of each joint are shown in Figure 5. According to Figure 5, the above curves of each joint of the 4-DOF parallel robot have no abrupt changes during the operation of the mechanism, which shows that the whole process of the mechanism runs smoothly. There is a little transient jerk, a little periodic or aperiodic oscillation phenomenon. When the robot moves from the initial position to the terminal position, the curves are continuous from zero, within the allowable limit of the mechanism joints, ensuring sufficient precision and smoothness of the optimal trajectory. The peak value of each movement is very close to the limit value of constraint so that the operating efficiency of the machine can be fully reflected. To make the above process more intuitive, the key parameters are presented in Table 4.  According to Figure 5, the above curves of each joint of the 4-DOF parallel robot have no abrupt changes during the operation of the mechanism, which shows that the whole process of the mechanism runs smoothly. There is a little transient jerk, a little periodic or aperiodic oscillation phenomenon. When the robot moves from the initial position to the terminal position, the curves are continuous from zero, within the allowable limit of the mechanism joints, ensuring sufficient precision and smoothness of the optimal trajectory. The peak value of each movement is very close to the limit value of constraint so that the operating efficiency of the machine can be fully reflected. To make the above process more intuitive, the key parameters are presented in Table 4.  Tables 3 and 4, the coincidence degree of the position of a starting point and endpoint between the initial set value and simulation result is 100%. The coincidence degree of the peak value of the limit angular velocity, the limit angular acceleration, the limit angular acceleration rate, and the maximum allowable limit value in each joint are all 91.050%. In theory, it is proved that the 12-phase sine jerk motion profile can meet the requirements of stability and accuracy of joint space trajectory of the 4-DOF parallel robot.

Comparative Analysis of Characteristic Parameters of 12-Phase, 3-Phase, and 15-Phase Sine Jerk Motion Profile
Experimentally, taking the three-phase and 15-phase sine jerk motion profile as the comparative basis to prove the superiority of the 12-phase sine jerk motion profile in trajectory planning, the experimental verification research is carried out around the characteristic parameters such as jerk amount, joint angular velocity, and joint energy consumption during a joint operation.

Comparison of Operational Stability
In the experimental platform shown in Figure 4a, the 3-phase, 12-phase, and 15-phase sine jerk motion profiles are sequentially input by offline programming in the robot motion program control language instruction, and the optimized time segments and j peak values are substituted.
For the straight-line trajectory from A to B in the operating space, it is set that the three-phase and 12-phase sine jerk motion profiles have the same starting and stopping point, and equaling motion time. In total, 101 sampling experimental data on Joint 1 were obtained. The experiment and simulation curves of each characteristic parameter on Joint 1 are drawn together, as shown in Figure 6, for showing the difference between simulation and experimental values of the robot under three-phase and 12-phase sine jerk motion profiles visually.
It can be seen from Figure 6 that the experimental and simulation curve of angular displacement are in good agreement under the three-phase or 12-phase sine jerk motion pro-file of Joint 1, which shows that the previous theoretical analysis is correct. The angular displacement of the 12-phase sine jerk motion profile is more smooth, which is a typical "S" type sine displacement curve. Compared with the angular dis-placement curve of the three-phase sine jerk motion profile, which is approximately a straight line, it can effectively alleviate the jerk during the operation of the mechanism and make the mechanism more flexible. Moreover, Figure 6b-d further shows that the angular velocity variation range, the magnitude of the joint's torque, and the speed of the torque direction change of the robot joint are more conducive to the stable operation of the robot and more conducive to reducing the mutual jerk between the joints when running 12-phase sine jerk motion than the three-phase sine jerk motion during the operation of the robot.
For the curve trajectory from C to D in the operating space, it is set that the three-phase and 15-phase sine jerk motion profiles have the same starting and stopping point and equaling motion time. A total of 402 sampling experimental data on Joint 1 were obtained. The experiment and simulation curves of each characteristic parameter on Joint 1 are drawn together as shown in Figures 7 and 8, for showing the difference between simulation and experimental values of the robot under 15-phase and 12-phase sine jerk motion profiles visually. Figure 7 shows the simulation experimental results of 12-phase sine jerk motion profile of Joint 1 as follows.
In the experimental platform shown in Figure 4a, the 3-phase, 12-phase, and 15-phase sine jerk motion profiles are sequentially input by offline programming in the robot motion program control language instruction, and the optimized time segments and values are substituted.
For the straight-line trajectory from A to B in the operating space, it is set that the three-phase and 12-phase sine jerk motion profiles have the same starting and stopping point, and equaling motion time. In total, 101 sampling experimental data on Joint 1 were obtained. The experiment and simulation curves of each characteristic parameter on Joint 1 are drawn together, as shown in Figure 6, for showing the difference between simulation and experimental values of the robot under three-phase and 12-phase sine jerk motion profiles visually. It can be seen from Figure 6 that the experimental and simulation curve of angular displacement are in good agreement under the three-phase or 12-phase sine jerk motion pro-file of Joint 1, which shows that the previous theoretical analysis is correct. The angular displacement of the 12-phase sine jerk motion profile is more smooth, which is a typical "S" type sine displacement curve. Compared with the angular dis-placement curve of the three-phase sine jerk motion profile, which is approximately a straight line, it can effectively alleviate the jerk during the operation of the mechanism and make the mechanism more flexible. Moreover, Figure 6b-d further shows that the angular velocity variation range, the magnitude of the joint's torque, and the speed of the torque direction change of the robot joint are more conducive to the stable operation of the robot and more conducive to reducing the mutual jerk between the joints when running 12-phase sine jerk motion than the three-phase sine jerk motion during the operation of the robot. For the curve trajectory from C to D in the operating space, it is set that the threephase and 15-phase sine jerk motion profiles have the same starting and stopping point and equaling motion time. A total of 402 sampling experimental data on Joint 1 were obtained. The experiment and simulation curves of each characteristic parameter on Joint 1 are drawn together as shown in Figures 7 and 8, for showing the difference between simulation and experimental values of the robot under 15-phase and 12-phase sine jerk motion profiles visually. Figure 7 shows the simulation experimental results of 12-phase sine jerk motion profile of Joint 1 as follows. To use the filter to process the data with a certain delay, the delay in Figure 7b, and c is not obvious, and the delay in Figure 7d is about 0.25 s. In the case of removing the effect of delay, it can be seen from Figure 7 that the experimental and simulation curve of angular displacement are in agreement under the 12-phase sine jerk motion profile of Joint 1, which shows that the previous theoretical analysis is correct. Figure 8 shows the simulation experimental results of 15-phase sine jerk motion pro-  In the initial stage of the movement, the error between the actual position curve and the planned trajectory is small, and then the error gradually increases. Finally, the actual position curve has two sudden changes. The reasons for this phenomenon are as follows. The 4-DOF redundant parallel robot has a small working space. It means shorter trajectory movement, acceleration, and velocity peaks are both constrained by displacement, and there is no need for a constant acceleration phase or a constant velocity phase in the trajectory. The 15-phase sine jerk motion profile contains a constant acceleration phase and a constant speed phase. Therefore, in a small working space, the 12-phase sine jerk motion profile is more applicable than the 15-phase curve. Therefore, it shows that the 15-phase sine jerk motion profile is not fit for the small operation motion.

Comparison of the Amount of Jerk on Each Joint
To further quantify the above analysis process and expand the jerk on each joint, based on the experimental data and the method given in Reference [26], assume that the angular displacement, angular velocity, angular acceleration, and angular acceleration rate of the joint are all restricted. The jerk on the whole moving process is measured by the integral of the square of the ( ) function, and then the calculation formula of the amount of jerk I of each joint is obtained as follows:  To use the filter to process the data with a certain delay, the delay in Figure 7b,c is not obvious, and the delay in Figure 7d is about 0.25 s. In the case of removing the effect of delay, it can be seen from Figure 7 that the experimental and simulation curve of angular displacement are in agreement under the 12-phase sine jerk motion profile of Joint 1, which shows that the previous theoretical analysis is correct. Figure 8 shows the simulation experimental results of 15-phase sine jerk motion profile of Joint 1 as follows. Figure 8 shows that the experimental and simulation curve of angular displacement is in a bad agreement under the 15-phase sine jerk motion profile of Joint 1. In the initial stage of the movement, the error between the actual position curve and the planned trajectory is small, and then the error gradually increases. Finally, the actual position curve has two sudden changes. The reasons for this phenomenon are as follows. The 4-DOF redundant parallel robot has a small working space. It means shorter trajectory movement, acceleration, and velocity peaks are both constrained by displacement, and there is no need for a constant acceleration phase or a constant velocity phase in the trajectory. The 15-phase sine jerk motion profile contains a constant acceleration phase and a constant speed phase. Therefore, in a small working space, the 12-phase sine jerk motion profile is more applicable than the 15-phase curve. Therefore, it shows that the 15-phase sine jerk motion profile is not fit for the small operation motion.

Comparison of the Amount of Jerk on Each Joint
To further quantify the above analysis process and expand the jerk on each joint, based on the experimental data and the method given in Reference [26], assume that the angular displacement, angular velocity, angular acceleration, and angular acceleration rate of the joint are all restricted. The jerk on the whole moving process is measured by the integral of the square of the jerk(t) function, and then the calculation formula of the amount of jerk I of each joint is obtained as follows: According to Equation (19), the jerk capacity of each joint is calculated as shown in Table 5 when running on the three-phase and 12-phase sine jerk motion profile respectively: The data in Table 5 directly show that the jerk of each joint, especially Joint 1, Joint 2, Joint 4, and Joint 5, decrease greatly, the jerk of the 12-phase sine jerk motion profile is greatly smaller than the three-phase sine jerk motion profile, and the force on each joint is more uniform. It proves, when the abovementioned 12-phase sine jerk motion profile is used, that the robot acceleration and deceleration process is more stable, with a better stability conclusion (lower jerk means more stable).

Speed Comparison
In order to respectively verify the rapidity of three-phase and 12-phase sine jerk motion profiles, under the process of end-effector motion, the peak angular velocity of each joint was obtained for quantitative comparison. The peak angular velocity of each joint is shown in Table 6. Table 6. Comparison of average peak angular velocity of joints' 3-phase and 12-phase sine jerk motion profile in the experiment. When comparing the peak angular velocity of each joint under the three-phase and 12-phase sine jerk motion profiles in Table 6, we see that the peak value of the angular velocity of the 12-phase sine jerk curve is 34.86% higher than the three-phase sine jerk curve. The robot can move at a higher speed by the 12-phase sine jerk motion profile, and each joint can fully reach the operation limit of key parameters. Moreover, compared with the limit angular velocity of each joint in Table 2, the 12-phase sine jerk motion profile can accelerate each joint to a higher speed; that is, each joint is closer to the limit of motion, so the operation efficiency is higher.

Comparison of Energy Consumption
The energy consumption is taken as the optimization objective in the trajectory optimization of the 12-phase sine jerk function curve. The energy consumption of the whole process is measured by the integral of torque square, according to the driving torque function M(t) in Lagrange dynamic Equation (20), combining with the calculation method of point-to-point track energy optimization in Reference [27]. The relationship between the energy consumption function P and the moment M(t) is as follows: According to Equation (20), the energy consumption of each joint of the three-phase and 12-phase sine jerk motion profile is calculated as shown in Table 7. Table 7. Comparison of energy consumption of joints' 3-phase and 12-phase sine jerk motion profile in the experiment.  According to Table 7, the total energy consumption of each joint of the 12-phase sine jerk motion profile is increased by 66.67% compared with that of the three-phase sine jerk motion profile. However, compared with Tables 4 and 5, the energy consumption of the 12-phase sine jerk motion profile can effectively reduce the transient jerk on each joint and improve the operation efficiency of each joint. Therefore, from the comprehensive comparison of operation accuracy, stability, and rapidity, the 12-phase sine jerk motion profile has more advantages and engineering application value.

Conclusions
To sum up, the high-speed motion trajectory optimized planning of the 4-DOF parallel robot was studied in this paper. By comparing the key characteristic parameters of three-phase and 12-phase sine jerk curves in simulation and experiment, the following conclusions are obtained: (1) When the joint space of the robot runs the 12-phase sine jerk motion profile, the overall operation is accurate and stable, and no sudden change of force and velocity exists in each joint, so the oscillation phenomenon is avoided. (2) Compared with the three-phase sine jerk curve, when the robot operates the 12-phase sine jerk motion profile, the force of each joint is more uniform, and the movement amount is closer to the maximum allowable limit value, so the operation efficiency is higher. (3) The energy consumption of the 12-phase sine jerk motion profile is higher than that of the three-phase sine jerk motion profile, but the former is far better than the latter in the stress uniformity and operation efficiency of each joint, which makes up for the defect of high energy consumption. As a result, the 12-phase sine jerk motion profile is still better than the three-phase sine jerk motion profile on the whole. (4) The 15-phase sine jerk motion profile is not suitable for trajectory planning in a small working space.
In future research, I will establish a dynamic equation closer to the motion model that includes friction force, based on the planned 12-phase sine impact motion profile. Then, I will design a controller to improve the accuracy of robot motion from the control level.