Dynamics Modeling of a Delta Robot with Telescopic Rod for Torque Feedforward Control

: This paper presents dynamics modeling of a Delta robot with three revolute legs and a telescopic rod. Firstly, two generalized coordinate systems are established to describe the relationship between the movement of the telescopic rod and the position of the moving platform, and the telescopic rod system kinematics are established through singularity analysis. Secondly, taking the telescopic rod as the research object, the corresponding dynamics model is established using the Euler–Lagrange method. Moreover, this paper proposes a method to convert the force exerted by the telescopic rod motion on the moving platform into actuator torques. Thirdly, the dynamics model of the Delta robot with a telescopic rod is established, and numerical simulations are performed to demonstrate this approach. Finally, the inﬂuence of the telescopic rod on the actuator torques is veriﬁed using an experiment. A comparison is drawn between the two dynamics models used in torque feedforward control to validate the proposed dynamics model.


Introduction
Parallel robots are mechanical structures different from a traditional serial manipulator, which generally consists of multiple closed motion loops. Reymond developed the parallel prototype of the manipulator in 1988, which is called the Delta robot [1]. The Delta Robot is a three-degrees-of-freedom space mechanism with three revolute legs, i.e., each leg is connected to the fixed platform and the moving platform by a set of spherical joints. Moreover, the motion of the Delta robot is controlled by three motors that are mounted on the base. The three motors actuate the three legs to realize the movement of the moving platform.
Dynamics modeling of the parallel robotic manipulators is important for the machine's motion control [2][3][4][5] and its corresponding structure design [6,7]. Many researchers have studied the dynamics modeling methods of parallel manipulators. The following are the main methods: (1) The Newton-Euler method establishes the dynamics model of the parallel manipulator. This method was applied by Kourosh E. Zanganeh to a six-degree-of-freedom parallel robot [8]. The parallel robot consists of six legs connected to a moving platform and a fixed platform. Each leg passes through a universal joint and a spherical joint to form two revolving joints; Chunxia Zhu used Newton's second law to derive the dynamics model of a 3-TPT flexible parallel robot [9]; Swain proposed a general dynamics model derivation method for manipulators [10]; George H. Pfreundschuha proposed a dynamics model for a pneumatic parallel robot with three degrees of freedom [11].
(2) The Lagrangian method establishes the dynamics model of the parallel manipulator G. Lebret usedthe Lagrangian method to analyze the dynamics of the Stewart paralle platform [12]. (3) Kane's method establishes the dynamics model of a parallel manipulator. This method is proposed by Kane [13]. Liu used the Huston form of Kane's equation to derive the positive dynamics equation of the Gough-Stewart manipulator [14]. (4) The virtual work principle establishes the dynamics model of a parallel manipulator Tsai used the virtual work principle and the Jacobian matrix to derive the inverse dynamics model of the Gough-Stewart manipulator [15]. F. Caccavale proposed a dynamics modeling of Tricept robot. Its parallel structure has three DOFs which are described by the axial translation of the radial link, and by two rotations about two axes orthogonal to the link itself [16].
In this paper, the object of study is a parallel robot with four degrees of freedom which is based on the structure of the Delta robot [17][18][19]. With reference to Figure 1, a telescopic rod connecting the fixed platform and the moving platform allows the flange on the moving platform to rotate around the Z-axis. In recent years, extensive research has been reported on dynamics analysis of the Delta robot [20][21][22][23]. However, these works focus on the traditional Delta robot and do no consider the effect of the telescopic rod on the actuator torques. In this paper, we use the Lagrangian method to establish a dynamics model of the telescopic rod. Moreover, these equations are implemented in an algorithm that is used to study the telescopic rod influ ence on actuator torques. The application of this algorithm is illustrated through a numer ical example that simulates the dynamic behavior of the telescopic rod and actuators while following a sample trajectory.

Kinematics
Manipulator kinematics studies the movement of connecting rods under geometric constraints. In general, manipulator kinematics includes two parts: forward kinematics and inverse kinematics. The purpose of the forward kinematics solution is to determine the position of the moving platform using the input of the actuators. The purpose of the inverse kinematics solution is to determine the input of the actuators by determining the position of the moving platform. In recent years, extensive research has been reported on dynamics analysis of the Delta robot [20][21][22][23]. However, these works focus on the traditional Delta robot and do not consider the effect of the telescopic rod on the actuator torques. In this paper, we use the Lagrangian method to establish a dynamics model of the telescopic rod. Moreover, these equations are implemented in an algorithm that is used to study the telescopic rod influence on actuator torques. The application of this algorithm is illustrated through a numerical example that simulates the dynamic behavior of the telescopic rod and actuators while following a sample trajectory.

Kinematics
Manipulator kinematics studies the movement of connecting rods under geometric constraints. In general, manipulator kinematics includes two parts: forward kinematics and inverse kinematics. The purpose of the forward kinematics solution is to determine the position of the moving platform using the input of the actuators. The purpose of the inverse kinematics solution is to determine the input of the actuators by determining the position of the moving platform.
Because this paper aims to build a complete dynamics model of the Delta robot with a telescopic rod, it is necessary to study the relationship between the motion of the telescopic rod and position of the moving platform. Therefore, we divide the kinematics into two parts. The first part is the Delta robot kinematics, and the second part is to obtain the telescopic rod kinematics under the condition that the telescopic rod is regarded as a subsystem. The telescopic rod kinematics model and Jacobian matrix are proposed, and are used to study the relationships among the movement of the telescopic rod, the movement of the moving platform, and the actuator torques.

Robot Structure
The structure diagram of the Delta robot with a telescopic rod is shown in Figure 2, where numbers 1 and 2 represent the fixed platform and the moving platform, respectively. Three identical mechanical arms are connected to the fixed platform and the moving platform. Each mechanical arm consists of an upper arm and a lower arm. The upper arms are labeled 3, 4 and 5, and the lower arms are labeled 6, 7 and 8. The upper arm is connected to the drive motor to form rotating joints of A 1 , A 2 and A 3 . The upper arms and the lower arms are connected by ball hinge joints of B 1 , B 2 and B 3 . The lower arms have a parallelogram structure and are connected with the moving platform through ball hinge joints of C 1 , C 2 and C 3 . The telescopic rod is labeled 9, and is directly connected to the moving platform and the fixed platform, which realizes the rotation of the end effector.
The coordinate system XYZ is attached to the center of the fixed platform at point O. The X-axis and Y-axis lie in the same plane as defined by the joints of A 1 , A 2 and A 3 . α i defines the angular orientation of the leg relative to the XYZ frame on the fixed platform.
Because this paper aims to build a complete dynamics model of the Delta robot with a telescopic rod, it is necessary to study the relationship between the motion of the telescopic rod and position of the moving platform. Therefore, we divide the kinematics into two parts. The first part is the Delta robot kinematics, and the second part is to obtain the telescopic rod kinematics under the condition that the telescopic rod is regarded as a subsystem. The telescopic rod kinematics model and Jacobian matrix are proposed, and are used to study the relationships among the movement of the telescopic rod, the movement of the moving platform, and the actuator torques.

Robot Structure
The structure diagram of the Delta robot with a telescopic rod is shown in Figure 2, where numbers 1 and 2 represent the fixed platform and the moving platform, respectively. Three identical mechanical arms are connected to the fixed platform and the moving platform. Each mechanical arm consists of an upper arm and a lower arm. The upper arms are labeled 3, 4 and 5, and the lower arms are labeled 6, 7 and 8. The upper arm is connected to the drive motor to form rotating joints of 1  Figure 2. Schematic of the Delta robot with a telescopic rod: (1) Fixed platform; (2) moved platform; (3)(4)(5) upper arms; (6-8) lower arms; (9) telescopic rod.

Forward Kinematics
The purpose of studying forward kinematics is to be able to obtain the position of the end effector through a given rotation angle. For this type of robot, the rotation angle of

Forward Kinematics
The purpose of studying forward kinematics is to be able to obtain the position of the end effector through a given rotation angle. For this type of robot, the rotation angle of the upper arm is known. The position of the center of the moving platform in the XYZ coordinate system is derived by the rotation angle. Moreover, the movement of the telescopic rod can obtain the position of the movable platform by establishing the generalized coordinate system of the telescopic rod.
where 1 k , 2 k and 3 k are defined as: The positions x, y and z of the center of the moving platform at angles 11 θ , 12 θ , 13 θ are solved by Equation (1). With reference to Figure 3, a coordinate system U i V i W i is established at the joint A i for each leg, such that the U i -axis is perpendicular to the axis of rotation of the joint A i . The V i -axis is along the joint axis V i . The angle θ 1i is measured from U i to AB. The angle θ 2i is defined from U i to BC. The angle θ 3i is measured from V i to BC. The length from the center of the fixed platform to the joint is R, and the length from the center of the moving platform to the joint is r. The lengths of the upper arms and the lower arms are denoted as l a and l b , respectively.
We can obtain the equation as follows: where k 1 , k 2 and k 3 are defined as: The positions x, y and z of the center of the moving platform at angles θ 11 , θ 12 , θ 13 are solved by Equation (1).

The Forward Kinematics of the Telescopic Rod System
The movement of the telescopic rod does not affect the position of the moving platform, but it will exert a force on the moving platform, which will affect the actuator torques. The kinematics of the telescopic rod is the basis for studying the influence of the telescopic rod motion on the actuator torques.
The telescopic rod is regarded as a subsystem of the robot whose base is a fixed platform and whose end is a moving platform. The base and the telescopic rod are connected by ball joints. The Grubler formula is used to analyze the degrees of freedom of the telescopic rod: The n is rigid body degrees of freedom. N represents the number of mechanical components. K is the number of joints and f i is the degree of freedom of the Kth joint. For the telescopic rod, n = 3, N = 1, K = 1. The ball joint has three degrees of freedom, so the telescopic rod system has three degrees of freedom.
(1) The model of the telescopic rod is simplified, which is shown in Figure 4. The l is the length of the telescopic rod. θ is the angle between the telescopic rod and the negative direction of the z-axis. φ is the angle between the projection of the telescopic rod on the XY plane and the positive direction of the x-axis. l, θ and φ are used as variables to describe the motion of the telescopic rod.

The forward kinematics of the telescopic rod system
The movement of the telescopic rod does not affect the position of the moving platform, but it will exert a force on the moving platform, which will affect the actuator torques. The kinematics of the telescopic rod is the basis for studying the influence of the telescopic rod motion on the actuator torques.
The telescopic rod is regarded as a subsystem of the robot whose base is a fixed platform and whose end is a moving platform. The base and the telescopic rod are connected by ball joints. The Grubler formula is used to analyze the degrees of freedom of the telescopic rod: The n is rigid body degrees of freedom. N represents the number of mechanical components. K is the number of joints and i f is the degree of freedom of the Kth joint. For the telescopic rod, n = 3, N = 1, K = 1. The ball joint has three degrees of freedom, so the telescopic rod system has three degrees of freedom.
(1) The model of the telescopic rod is simplified, which is shown in Figure 4. The l is the length of the telescopic rod.  is the angle between the telescopic rod and the negative direction of the z-axis.  is the angle between the projection of the telescopic rod on the XY plane and the positive direction of the x-axis. l, and  are used as variables to describe the motion of the telescopic rod.  Defining l, θ,and φ as the axis variable. The axis variables and the position of the moving platform have the following relationship: (2) Generalized coordinate system 2 is shown in Figure 5. l is the length of the telescopic rod. θ is the angle between the projection of the rod in the XZ plane and the negative direction of the z-axis. φ is the angle between the projection of the rod in the YZ plane and the negative direction of the Z-axis.
(2) Generalized coordinate system 2 is shown in Figure 5. l is the length of the telescopic rod.  is the angle between the projection of the rod in the XZ plane and the negative direction of the z-axis.  is the angle between the projection of the rod in the YZ plane and the negative direction of the Z-axis. Similarly to generalized coordinate system 1, the axis variables and the position of moving platform have the following relationship: We can obtain the forward kinematics of the telescopic rod from Equation (4): Similarly to generalized coordinate system 1, the axis variables and the position of moving platform have the following relationship: We can obtain the forward kinematics of the telescopic rod from Equation (4):

Inverse Kinematics
A. The inverse kinematics of the Delta robot The first part of inverse kinematics is to solve the angle of the upper arm given the target position. By analyzing the structure of the leg i, we can readily derive the expression below [24]: p ui , p vi and p wi represent the position of P in the U i V i W i coordinate frame. The expressions for these variables are given by: Therefore, θ 3i can be obtained by Equation (8): We can obtain Equation (11) by summing the square of formula (7) and formula (9): Moreover, using Equation (11), we can write: where: Half-angle tangent is defined as: where h i can be obtained by solving a polynomial Equation (12). Then, θ 1i can be solved using Equation (14): B. The inverse kinematics of the telescopic rod system The second part of the inverse kinematics solution is to obtain the description of the motion of the telescopic rod in the generalized coordinate system through the position of the moving platform.
The inverse kinematics in generalized coordinate system 1 is given as: The inverse kinematics in generalized coordinate system 2 is given as:

Jacobian
A. The Jacobian of the Delta robot The first part of the Jacobian matrix provides a transformation from the velocity of the end-effector in Cartesian space to the actuated joint velocities, as shown in Equation (17): where v= v x v y v z T is the speed vector of the moving platform, and T is the rotation speed vector of the upper arms. The Jacobian matrix of this part is defined as where: Ja=   Ja 11 Ja 12 Ja 13 Ja 21 Ja 22 Ja 23 Ja 31 Ja 32 Ja 33   The Jacobian of the telescopic rod system Differentiating both sides of Equation (2) with respect to time, we obtain: Equation (20) can be derived by converting Equation (19) into matrix form: where: Similarly, differentiating both sides of Equation (5) with respect to time, we obtain: where: The kinematics and Jacobian matrices of telescopic rods in two generalized coordinate systems are established, respectively. In order to determine the kinematics of the telescopic rod in that generalized coordinate system, a singularity analysis is carried out on it. The singularities occur when the following condition is satisfied: The determinant of J t1 is zero when θ is equal to zero. At the same time, this also means that when the position of the moving platform passes (0, 0, z), the value φ will change suddenly. We have generated a linear trajectory in the Cartesian space from point (100, −100, −800) to (100, 100, −800), which is shown in Figure 6a. We can see from Figure 6b that the value of φ has a big mutation, which will cause the Jacobian matrix at this time to be ill-conditioned. However, the representation of the same trajectory in generalized coordinate system 2, which has no singularity, is shown in Figure 6c.
Therefore, generalized coordinate system 2 is used to describe the motion of the telescopic rod, because the determinant of J −1 t2 will not appear to be zero matrices. The Jacobian matrix of the telescopic rod and the Jacobian matrix of the Delta robot without a telescopic rod are derived through kinematics. This is important for establishing the dynamic model with the telescopic rod.
(100, −100, −800) to (100, 100, −800), which is shown in Figure 6a. We can see from Figure   6b that the value of φ has a big mutation, which will cause the Jacobian matrix at this time to be ill-conditioned. However, the representation of the same trajectory in generalized coordinate system 2, which has no singularity, is shown in Figure 6c.
Therefore, generalized coordinate system 2 is used to describe the motion of the telescopic rod, because the determinant of −1 2 t J will not appear to be zero matrices. The Jacobian matrix of the telescopic rod and the Jacobian matrix of the Delta robot without a telescopic rod are derived through kinematics. This is important for establishing the dynamic model with the telescopic rod.

Dynamics
The actuator torques of traditional Delta robot are mainly used to drive the upper arms, lower arms and moving platform. However, the Delta robot with a telescopic rod also needs to drive the telescopic rod, which means that the original dynamics model can no longer accurately represent the relationship between the motion of robot and the actuator torques.
The movement of the telescopic rod is produced by the moving platform, and the actuator torque of the moving platform is produced by the servo motor. The force of the telescopic rod acting on the moving platform and the force of the moving platform acting on the telescopic rod are a pair of interaction forces. The movement of the telescopic rod is analyzed to obtain the force on the moving platform, and then to obtain the torque for the three actuators.

Dynamics of Telescopic Rod
In the above section, the relationship between the axis variables Φ and the position of the moving platform is derived. Further, the Euler-Lagrange method is used to establish the dynamics equation of the telescopic rod, which is used to describe the relationship between the motion of the telescopic rod and the torque applied by the moving platform to the telescopic rod.
The Lagrangian function of the telescopic rod can be expressed as: where Γ and Ε represent the total kinetic energy and potential energy of the telescopic rod, respectively.

Dynamics
The actuator torques of traditional Delta robot are mainly used to drive the upper arms, lower arms and moving platform. However, the Delta robot with a telescopic rod also needs to drive the telescopic rod, which means that the original dynamics model can no longer accurately represent the relationship between the motion of robot and the actuator torques.
The movement of the telescopic rod is produced by the moving platform, and the actuator torque of the moving platform is produced by the servo motor. The force of the telescopic rod acting on the moving platform and the force of the moving platform acting on the telescopic rod are a pair of interaction forces. The movement of the telescopic rod is analyzed to obtain the force on the moving platform, and then to obtain the torque for the three actuators.

Dynamics of Telescopic Rod
In the above section, the relationship between the axis variables Φ and the position of the moving platform is derived. Further, the Euler-Lagrange method is used to establish the dynamics equation of the telescopic rod, which is used to describe the relationship between the motion of the telescopic rod and the torque applied by the moving platform to the telescopic rod.
The Lagrangian function of the telescopic rod can be expressed as: where Γ and E represent the total kinetic energy and potential energy of the telescopic rod, respectively. The Euler-Lagrange equation of the telescopic rod is expressed as follows: The structure of the telescopic rod is simplified into two parts, which are divided into an upper rod and a lower rod. The moment of inertia of the upper rod is: The moment of inertia of the lower rod is: where m l1 and m l2 are the masses of the upper rod and lower rod, respectively. The kinetic energy of the upper rod and lower rod can be given as: Similarly, we can obtain the potential energy equation of the upper rod and lower rod as follows: We can readily obtain the telescopic rod Lagrangian function and Lagrangian equation through Equations (23) and (24).

Dynamics Modeling of Delta Robot with Telescopic Rod
On the one hand, the driving torques generated by the actuators make the telescopic rod move; on the other hand, they make the moving platform move by driving the upper arms. In order to obtain the driving torque through the motion of the moving platform, a complete dynamics model of the four-DOF parallel robot needs to be established.
Newton-Euler equations of motion are used to determine the torques applied by actuators without a telescopic rod [24,25]. This dynamics model assumes that the mass of each lower rod is concentrated at the joints B i . Based on these assumptions, the equation of motion is written by summing the moments about the actuated joint for the ith leg: where ∑ τ Ai is the sum of the moments applied at joint A i ; I u is the mass moment of inertia of the upper rod and the lower rod; τ pi is the ith element of τ p , which is an array of the inertial loads at joint A i due to the acceleration of the moving platform; and k is the viscous damping coefficient of the actuators. The expressions for I u and τ p are given by Equations (36) and (37): where I m is the mass moment of inertia of the motor rotor, and a p is the acceleration of the moving platform. Taking into account the motor torques and the gravitational force, ∑ τ Ai can also be expressed as: Then, we can obtain actuator torques using Equations (35) and (38).
Furthermore, the total moment of the three actuators is given as: where F = F 1 F 2 F 3 can be obtained using Equation (33).

Simulations
To clearly present the proposed method, a numerically computed example will be discussed in this section. A trajectory for the moving platform and a set of manipulator design parameters are assumed. The gate-shaped trajectory is selected to represent a typical motion that might be used during a practical application of the Delta robot.
The Delta robot design parameters that are used for the simulations are given as: can be obtained using Equation (33).

Simulations
To clearly present the proposed method, a numerically computed example will be discussed in this section. A trajectory for the moving platform and a set of manipulator design parameters are assumed. The gate-shaped trajectory is selected to represent a typical motion that might be used during a practical application of the Delta robot.
The Delta robot design parameters that are used for the simulations are given as:  The axis variables of the telescopic rod can be obtained by the movement of the moving platform. The corresponding time-histories of the components of the position, velocity, and acceleration of the axis variables of the telescopic rod are shown in Figure 8a-c. The moment of the telescopic rod on the moving platform is shown in Figure 9d. Given this sample trajectory and set of manipulator design parameters, the actuator torques of the Delta robot are calculated without the telescopic rod using the proposed approach. Next, using same trajectory and Delta robot design parameters, the actuator torques of the Delta robot with a telescopic rod are calculated. Friction is ignored when calculating the actuator torques. The results from these two simulations are shown in Figures 9 and 10. Given this sample trajectory and set of manipulator design parameters, the actuator torques of the Delta robot are calculated without the telescopic rod using the proposed approach. Next, using same trajectory and Delta robot design parameters, the actuator torques of the Delta robot with a telescopic rod are calculated. Friction is ignored when calculating the actuator torques. The results from these two simulations are shown in Figures 9 and 10.
From the simulation results, it can be seen that the actuator torques with the telescopic rod and without the telescopic rod under the same trajectory are different, and the actuator torques with the telescopic rod are larger than those without the telescopic rod. This is because when there is a rod, the drive shaft not only drives the movement of the moving platform but also drives the movement of the telescopic rod.  From the simulation results, it can be seen that the actuator torques with the telescopic rod and without the telescopic rod under the same trajectory are different, and the actuator torques with the telescopic rod are larger than those without the telescopic rod. This is because when there is a rod, the drive shaft not only drives the movement of the moving platform but also drives the movement of the telescopic rod.

Experimental Application
An experimental study is now presented to test the dynamics model of the Delta robot with a telescopic rod. First, the actuator torques with and without the telescopic rod are compared under the same trajectory, which is used to verify the effect of the telescopic rod on the actuator torques. Second, a robot control system based on the torque feedforward control strategy is built, and the effects of the traditional Delta robot dynamics  From the simulation results, it can be seen that the actuator torques with the telescopic rod and without the telescopic rod under the same trajectory are different, and the actuator torques with the telescopic rod are larger than those without the telescopic rod. This is because when there is a rod, the drive shaft not only drives the movement of the moving platform but also drives the movement of the telescopic rod.

Experimental Application
An experimental study is now presented to test the dynamics model of the Delta robot with a telescopic rod. First, the actuator torques with and without the telescopic rod are compared under the same trajectory, which is used to verify the effect of the telescopic rod on the actuator torques. Second, a robot control system based on the torque feedforward control strategy is built, and the effects of the traditional Delta robot dynamics

Experimental Application
An experimental study is now presented to test the dynamics model of the Delta robot with a telescopic rod. First, the actuator torques with and without the telescopic rod are compared under the same trajectory, which is used to verify the effect of the telescopic rod on the actuator torques. Second, a robot control system based on the torque feedforward control strategy is built, and the effects of the traditional Delta robot dynamics model and the dynamics model proposed in this paper on trajectory tracking are compared.
The experimental installation constitutes of a Delta robot with a telescopic rod, in which the original control unit has been replaced by an open control system developed by ourselves. The mechanical parameters can be seen in Table 1. The experimental installation is shown in Figure 11. model and the dynamics model proposed in this paper on trajectory tracking are compared. The experimental installation constitutes of a Delta robot with a telescopic rod, in which the original control unit has been replaced by an open control system developed by ourselves. The mechanical parameters can be seen in Table 1. The experimental installation is shown in Figure 11.  Figure 11. The Delta robot for experimentation.

Comparing Actuator Torques
The above section discusses the effect of the telescopic rod motion on the actuator torques, and to verify this, an experiment is carried out. The robot is allowed move along a standard gate shape and collect the actuator torques. The actuator torques are collected again, along the same trajectory, without the telescopic rod.
Two different sets of actuator torques, obtained from data measured from the Delta robot, are compared-the first one with a telescopic rod and the second one without a telescopic rod. A comparison between the actuator torques with telescopic rod and actuator torques without telescopic rod is shown in Figure 12. It can be seen in the figure that the actuator torques are different in the two cases, and the actuator torques of the Delta robots with the telescopic rod are larger than those of the Delta robot without the telescopic rod when performing the same movement. The result indicates that the telescopic rod does have an effect on the actuator torques, and it is meaningful to study the effect of the torque produced by the telescopic rod.

Comparing Actuator Torques
The above section discusses the effect of the telescopic rod motion on the actuator torques, and to verify this, an experiment is carried out. The robot is allowed move along a standard gate shape and collect the actuator torques. The actuator torques are collected again, along the same trajectory, without the telescopic rod.
Two different sets of actuator torques, obtained from data measured from the Delta robot, are compared-the first one with a telescopic rod and the second one without a telescopic rod. A comparison between the actuator torques with telescopic rod and actuator torques without telescopic rod is shown in Figure 12. It can be seen in the figure that the actuator torques are different in the two cases, and the actuator torques of the Delta robots with the telescopic rod are larger than those of the Delta robot without the telescopic rod when performing the same movement. The result indicates that the telescopic rod does have an effect on the actuator torques, and it is meaningful to study the effect of the torque produced by the telescopic rod. (c) Figure 12. The actuator torques with telescopic rod and without telescopic rod (a) 1: the first joint actuator torque with telescopic rod, 2: the first joint actuator torque without telescopic rod; (b) 1: the second joint actuator torque with telescopic rod, 2: the second joint actuator torque without telescopic rod; (c) 1: the third joint actuator torque with telescopic rod, 2: the third joint actuator torque without telescopic rod.

Dynamics Parameter Identification
An important application of dynamics is to achieve more precise control, so the proposed dynamics are applied to torque feedforward control to compare the tracking errors with traditional robot dynamics. The torque feedforward control strategy needs to calculate the actuator torques in real-time according to the motion of the robot, which requires dynamics identification of the robot. The dynamics parameter identification of the Delta robot with a telescopic rod is divided into two parts. The first part is to identify the inertial parameters of the three legs. The second part is to identify the inertial parameters of the telescopic rod according to Equation (40).
The optimized trajectory used to perform the parameter identification is parameterized as a finite Fourier series according to the method described in Ref. [26]. The identification algorithm can be divided into two steps: (1) Dynamics parameter identification for Delta robot without telescopic rod The friction force is a non-negligible item in the actuator torques, so the friction force is considered based on Equation (39). The actuator torques without the telescopic rod can be expressed as: where f F is the friction force, and can be described as a coulomb viscous friction model. Figure 12. The actuator torques with telescopic rod and without telescopic rod (a) 1: the first joint actuator torque with telescopic rod, 2: the first joint actuator torque without telescopic rod; (b) 1: the second joint actuator torque with telescopic rod, 2: the second joint actuator torque without telescopic rod; (c) 1: the third joint actuator torque with telescopic rod, 2: the third joint actuator torque without telescopic rod.

Dynamics Parameter Identification
An important application of dynamics is to achieve more precise control, so the proposed dynamics are applied to torque feedforward control to compare the tracking errors with traditional robot dynamics. The torque feedforward control strategy needs to calculate the actuator torques in real-time according to the motion of the robot, which requires dynamics identification of the robot. The dynamics parameter identification of the Delta robot with a telescopic rod is divided into two parts. The first part is to identify the inertial parameters of the three legs. The second part is to identify the inertial parameters of the telescopic rod according to Equation (40).
The optimized trajectory used to perform the parameter identification is parameterized as a finite Fourier series according to the method described in Ref. [26]. The identification algorithm can be divided into two steps: (1) Dynamics parameter identification for Delta robot without telescopic rod The friction force is a non-negligible item in the actuator torques, so the friction force is considered based on Equation (39). The actuator torques without the telescopic rod can be expressed as: where F f is the friction force, and can be described as a coulomb viscous friction model.
Linearizing Equation (41), we obtain: where I is the parameter to be identified. The identification parameters can be obtained using the least squares method.
The identified inertial parameters are shown in Table 2. To verify the effect of inertial parameter identification, a verification experiment is carried out. The dynamics model is used to predict the actuator torques in real time and compare them with the real actuator torques. The experimental results are shown in Figure 13. (c) Figure 13. Comparison between the predicted and real actuator torques (a) 1: predicted actuator torque of first joint, 2: real actuator torque of first joint; (b) 1: predicted actuator torque of second joint, 2: actuator torque of second joint; (c) 1: predicted actuator torque of third joint, 2: actuator torque of third joint.

Experiment of Torque Feedforward Control
The torque feedforward strategy was implemented successfully on the Delta robot and is sufficient for the implementation of all desired control-algorithm elements [2]. The scheme used for the control of the Delta robot is shown in Figure 14. It is composed of a feedforward block of the dynamics model of the Delta robot. Two dynamics models are used, based on Equation (39) and Equation (40), respectively. The control system takes the position of the moving platform (Xd) as input. The desired joint position (qd) is obtained Figure 13. Comparison between the predicted and real actuator torques (a) 1: predicted actuator torque of first joint, 2: real actuator torque of first joint; (b) 1: predicted actuator torque of second joint, 2: actuator torque of second joint; (c) 1: predicted actuator torque of third joint, 2: actuator torque of third joint.

Experiment of Torque Feedforward Control
The torque feedforward strategy was implemented successfully on the Delta robot and is sufficient for the implementation of all desired control-algorithm elements [2]. The scheme used for the control of the Delta robot is shown in Figure 14. It is composed of a feedforward block of the dynamics model of the Delta robot. Two dynamics models are used, based on Equation (39) and Equation (40), respectively. The control system takes the position of the moving platform (X d ) as input. The desired joint position (q d ) is obtained using inverse kinematics, and the accelerations are obtained using numerical differentiation of the desired joint angles and moving platform position. The controller is a standard proportional-derivative.
To be able to evaluate the performance of the torque feedforward control algorithm, which uses the dynamics model of the Delta robot with a telescopic rod proposed above, it may be meaningful to compare it with the dynamics model of the Delta robot without a telescopic rod. The tracking errors of the three joints obtained using two different dynamics models are shown in Figure 15. The tracking error is roughly reduced by half using the dynamics model proposed in this paper. To be able to evaluate the performance of the torque feedforward control algorithm, which uses the dynamics model of the Delta robot with a telescopic rod proposed above, it may be meaningful to compare it with the dynamics model of the Delta robot without a telescopic rod. The tracking errors of the three joints obtained using two different dynamics models are shown in Figure 15. The tracking error is roughly reduced by half using the dynamics model proposed in this paper.
These results show that the dynamics model of the Delta robot proposed in this paper can better describe the actuator torques, which is beneficial for realizing the precise control of the robot. To be able to evaluate the performance of the torque feedforward control algorithm, which uses the dynamics model of the Delta robot with a telescopic rod proposed above, it may be meaningful to compare it with the dynamics model of the Delta robot without a telescopic rod. The tracking errors of the three joints obtained using two different dynamics models are shown in Figure 15. The tracking error is roughly reduced by half using the dynamics model proposed in this paper.
These results show that the dynamics model of the Delta robot proposed in this paper can better describe the actuator torques, which is beneficial for realizing the precise control of the robot.

Conclusions
In this article, dynamics modeling of a Delta robot with a telescopic rod was presented and the proposed dynamics model was applied to the dynamics feedforward control, to verify that the model can describe the actuator torques more accurately. The telescopic rod was regarded as a subsystem of the robot, its degrees of freedom were analyzed, and the kinematics model was deduced under the established generalized coordi- These results show that the dynamics model of the Delta robot proposed in this paper can better describe the actuator torques, which is beneficial for realizing the precise control of the robot.

Conclusions
In this article, dynamics modeling of a Delta robot with a telescopic rod was presented and the proposed dynamics model was applied to the dynamics feedforward control, to verify that the model can describe the actuator torques more accurately. The telescopic rod was regarded as a subsystem of the robot, its degrees of freedom were analyzed, and the kinematics model was deduced under the established generalized coordinate system. The dynamics modeling of the telescopic rod was then evaluated, based on the Euler-Lagrange method and Jacobian matrix of the telescopic rod, which projects forces acting on telescopic rod onto the actuator motors. This dynamics modeling was applied to inertial parameter identification experiments of the Delta robot with a telescopic rod, and the results of the inertial parameter identification experiments were used in the torque feedforward control experiment. The torque feedforward control experiment shows that using the dynamics modeling proposed in this paper can achieve more precise control of the Delta robot with a telescopic rod. Thus, the present study provides a framework for future research on the design and control of this type of parallel robot.