Dynamic modeling and performance analysis of a lower-mobility parallel robot with a rotatable platform

: Recently, applications of high-speed, lightweight parallel robots have been gaining increasing interest. Studies have shown that their elastic that their elastic deformation during operation often affects the robot's dynamic performance. In this paper, we designed and studied a 3 DOF parallel robot with a rotatable working platform. We developed a rigid-flexible coupled dynamics model consisting of a fully flexible rod and a rigid platform by combining the Assumed Mode Method with the Augmented Lagrange Method. The driving moments under three different modes were used as feedforward in the model's numerical simulation and analysis. We conducted a comparative analysis demonstrating that the flexible rod's elastic deformation under a redundant drive is significantly smaller than that of a non-redundant one, leading to a better suppression effect on vibration. The system's dynamic performance under the redundant drive was significantly superior compared to that of the non-redundant one. Additionally, the motion accuracy was higher and the driving mode b was better than that of the driving mode c. Finally, the proposed dynamics model's correctness was verified by modeling it in Adams.


Introduction
Parallel robots possess features of high stiffness, high load-bearing capacity, low cumulative error, and high motion accuracy [1].They have been widely used in processing and manufacturing, aerospace, and medical fields [2][3][4][5][6][7].However, they also have some disadvantages, such as smaller working spaces, more singularities, and difficulty in dynamic solutions.These limitations have restricted the development of parallel robots to some extent.
Due to their structural features, parallel robots have several singularities.Near the singularity, the robot's performance sharply drops.In more critical cases, it may lead to control loss [1].Therefore, when designing and applying parallel robots, singular dislocations should be avoided, especially type II singularities [8].The redundancy method is a widely-used method that permits the avoidance of such singularities [9][10][11][12][13].This method is further divided into the kinematic [9,10] and the drive redundancy methods [11][12][13].Compared to the prior, the latter is more operational because it does not introduce a new DOF in the elimination process.
With the industry's rapid development, demand for the dynamic performance of parallel robots has been increasing.To obtain superior dynamic characteristics, parallel robots are often designed to be light.However, their flexible components are prone to elastic deformations during high-speed motion, potentially causing motion errors in the moving platform in addition to elastic vibration or resonance in the mechanism.This consequently results in the parallel mechanism's kinematic and dynamic performance degradation.Therefore, it is of the utmost importance to consider the elastic deformation of components during the modeling process, in addition to establishing a flexible model.To simplify this study, we discretized the components and used the following main methods: 1) the Finite Element Method (FEM) [14,15] ; 2) the Assumed Mode Method (AMM) [16,17] ; 3) the finite section method [18] ; and 4) the lumped parameter method [19].The assumed mode method could be further combined with the modal truncation technique to truncate the dominant low-order modes.After that, we can use the LaGrangian formulation and Kane's equation to establish the system's dynamic equations with lower dimensions.Therefore, the assumed mode method is suitable for the construction of active control strategies for flexible robots.
Academic research on the treatment of rigid-flexible kinematic relationships of flexible bodies has led to the development of various modeling methods.The most widely used of these are the following: 1) The Kinematic Elastic Dynamics (KED); 2) the Floating Frame of Reference (FFRF); and 3) the Absolute Nodal Coordinate Formulation (ANCF).The floating frame of reference proved to be a reliable method of establishing the rigid-flexible coupled dynamics model.
Recently, there have been some studies on the dynamic modeling of flexible parallel robots.For example, Liang Dong et al [20] studied a multi-drive 2 DOF parallel robot and established a systematic rigid-flexible coupling dynamics model by combining the assumed mode method with the Augmented Lagrange method.They also conducted a study on the active control strategy.İDer, S. K et al [21] established a dynamic model for flexible joints and studied the inverse dynamic control of flexible joint parallel robots.Sheng et al [22] developed a dynamics model for a 3-RRR parallel robot with a flexible intermediate linkage and investigated the vibration features.Cammarata et al [23] proposed an MSA/CMS framework to describe the complete and simplified model of fully flexible parallel robots' dynamic elasticity.
We present in this paper a novel 4-RRR redundantly-driven parallel robot.First, we discretized the flexibility of the system's arbitrary rod, j, based on the assumed mode method.Then, we established the flexible rod's dynamic equations.Using the Lagrangian formulation, we established the endeffector's rigid platform's dynamic equations.We also categorized the eight rods into four branch chains to establish the open-chain system's dynamic equations.We established the complete system's dynamic equations on the system's closed-loop vector equations and transformed them into the closed-loop system dynamic equations in the form of independent generalized system coordinates.
We regarded the robot's eight rods as flexible bodies while the end-effector platform as a rigid body.They thus formed a rigid-flexible coupled system and were numerically simulated.In these simulations, the dynamic performance was discussed under 3 types of drive-force distribution: a) a non-redundant drive; b) a redundant drive with minimum input torque; and c) a redundant drive with minimum input energy.Based on the inverted rigid body dynamics model, we respectively calculated three sets of corresponding driving moments.Taking the previously mentioned modes' driving moments as feedforward inputs, numerical calculations were conducted in MATLAB to obtain the system's dynamic response.To provide a reference for control and strategy, we then compared and analyzed the dynamic performance differences under these three modes.
We developed in this paper a dynamic model of a rigid-flexible coupled parallel manipulator with a rigid end-effector.Moreover, we analyzed the model's dynamic performance under three driving modes.The modular modeling method was conducive to the generalization to other models; additionally, it was helpful for the spatial structure's modeling.The established rigid-flexible coupling system dynamics model has proved to be beneficial for the optimization of light design and dynamic performance of parallel robots.Moreover, it is conducive to the development of subsequent control strategies.

Model
We present in this paper a planar 3 DOF parallel robot design as shown in Figure 1.It contains a total of four branch chains, each of which contains two rods and a square end-effector platform.The four branch chains are connected to the end-effector platform's four corners.Among all the branch chains, the rods labeled 1-4 are active while the 5-8 ones are driven.We placed the robot in the inertial coordinate system  −  − .Let  denote the solid joint coordinates between the i th branch chain and the inertial coordinate system, its coordinates denoted as (   ,    ) , and its angle as   ( = 1, ⋯ ,4) .Let  denote the position of the joint between the two rods of the i th branch chain, its coordinate denoted as (   ,    ), and its angle as   ( = 5, ⋯ ,8).Let   denote the position of the joint at the end of the follower rod of the i th branch chain, its coordinate denoted as (   ,    ).The operating platform's center point's position is denoted as P, having coordinates of (  ,   ).We denote the angle between the line of the center point P,   , and the  −  axis as  .

The discretization of an arbitrary flexible rod j
According to the floating coordinate frame method, the motion of an arbitrary flexible rod j can be decomposed into two parts.One is a large range of rigid bodies solidly connected to the floating coordinate system and the other is a small range of elastic deformation motions relative to it.Thus, we selected the arbitrary rod's near-frame endpoint as the origin of its floating coordinate system.The xdirection is the rod's axial direction while the y-direction is its tangential one.At this point, we simulated the deformation of the flexible rod j as the cantilever beam's deformation.Figure 2 shows this deformation while moving.In Figure 2,  −  −  represents the global coordinate system.The origin O is fixed in the plane of motion.  −   −   represents the floating coordinate system fixed on the flexible rod j.The position vector of the origin   in the overall coordinate system is    = [       ]  .Before the flexible rod's motion deformation,  is a point on it and its position vector in the floating coordinate system   −   −   is  ,  = [ ̄0 ]  .After motion deformation, the deformation vector diameter of the point   in the floating coordinate system becomes    = [       ]  .The position vector of the point   in the overall coordinate system is denoted as    .According to Fig. 2, the following position-vector relationship can be obtained as follows: In ( In the equation,    (̄, ) and    (, ) denote the axial and transverse elastic deformations on the flexible rod j at a distance ̄ from the origin under the floating coordinate system   −   −   at the moment of motion t, respectively.
In the equation， .
After discretizing the flexible rod j in the assumed mode method, we can express the position vector of the arbitrary point Pj after the rod's kinematic deformation in the overall coordinate system as: In the equation， By taking the first order derivative of equation ( 4) with respect to time, we can obtain the velocity vector at the global coordinates of point Pj as: In the equation ∈ ℝ (3+  +  )×1 are defined as a column matrix of the flexible rod's generalized coordinates.It contains both rigid generalized coordinates,    and   , for a wide range of variation and flexible generalized coordinates,    , for a small one.

The arbitrary flexible rod's kinetic energy
By integrating the flexible rod's microelements' kinetic energy, we obtained its kinetic energy at moment t as: In the equation， denotes the flexible rod j's density; A denotes its cross-sectional area;   denotes its length;  , (  ) = ∫      0   denotes its mass matrix, which is a function of   .

The arbitrarily flexible rod's potential energy
In this paper, the 4-RRR parallel robot always moves in the horizontal plane, and the gravitational potential energy is constant; thus, we can nullify the latter.This implies that we only need to calculate the mechanism's elastic potential energy.Neglecting the rod's shear deformation and considering only the deformation and bending ones, we can obtain its potential energy as: In the equation, E denotes the flexible rod's material's modulus of elasticity.I denotes the crosssectional's area moment of inertia, and  , denotes its stiffness matrix.

The flexible rod's kinetic equations
We can express the Lagrangian formulation in the following form: Where L is the difference between the kinetic and potential energies,  =  , −  , ,   is the generalized coordinate, ̇ is the generalized velocity, and   is the generalized force corresponding to the generalized coordinate.
In the equation,  , (  ) and  , denote the flexible rod's mass and stiffness matrices, respectively., denotes its generalized external moment column matrix., (  , ̇) denotes its centrifugal Force/Coriolis force matrix with the following expressions: where I denotes the unit matrix and ⨂ the Kronecker Product symbol.

The flexible rod's kinetic equations in the system's generalized coordinate form
For this paper's system, we can express the near-frame end joint coordinates    of the driven rod     of the arbitrary branch chain i in terms of the drive rod's angle of rotation and deformation modal coordinates,    .Therefore, we chose all the rods' angles of rotation and modal coordinates as the generalized coordinates.Since the end effector is a rigid body, its center point coordinate,   , and its rotation angle,  1 , are chosen as its generalized coordinates.By grouping all flexible rods with the end effector, a total of eight flexible rods and one rigid operator are included.The whole system's generalized coordinates are denoted as   (non-independent coordinates) and are specified below： To represent the flexible bar's generalized coordinates   in the system's generalized coordinates, it is necessary to express each quantity in   as   .
The arbitrary flexible rod's angular coordinates of rotation,   can be expressed as: In the equation,    denotes the transfer matrix between the flexible rod's angular coordinates of rotation and the system's generalized coordinates, which is a Boolean matrix; i.e.: its elements are either 0 or 1.We specify it as: Similarly, we can express the arbitrary flexible rod's array of modal coordinates,    , as: In the equation,     denotes the transfer matrix between the array of the arbitrary flexible rod's modal coordinates and the system's generalized coordinates.It is also a Boolean matrix specified as: Due to the mechanism's configuration, we cannot express the flexible rod's near-frame end joint coordinates,    , as a general formula.Thus, it is necessary to divide each rod's near-frame end joint coordinates into two categories expressed as a function of the system's generalized coordinates,   , respectively： ( ) ( ) ( 1, 2,3, 4) ( 5, 6, 7,8) 0 Taking each of these two equations first order derivative with respect to time t, we can obtain： According to equations ( 12), (14), and ( 17), we can express the rod's generalized coordinate column matrix's first-order derivative, ̇, in terms of the system's generalized coordinates' first-order derivative, ̇, as:
To organize equation (19)  force/Coriolis force, stiffness, and generalized external moment column matrices expressed in the form of the system's generalized coordinates.

The end effector's kinetic equations
The end effector has a square structure and high stiffness.We can consider its motion to be purely rigid.We denoted its position and posture vector by   = [      1 ]  , its position vector as   = [     ]  =  1   , and its rotation angle as  1 =  2   , where  1 and  2 are the transformation matrices expressed as: The end effector's kinetic energy, Tp , consists of both translational and rotational kinetic energies expressed as follows:  Where   denotes the end effector's mass,   its rotational inertia, and   its mass matrix, which is expressed as   =    1   1 +    2   2 .As shown in Figure 3, we theoretically divided the end manipulator's four jointed sections and analyzed them separately.Since the end effector constantly is in a pure rigid motion, we can set both its gravitational and elastic potential energies as zero.Moreover, its kinetic energy can be plugged into the Lagrangian formulation to obtain its kinetic equation: In the equation,   denotes the column matrix of the equivalent external force on the center of mass of the endeffector and its external torque.
To represent its kinetic equations as a system in a generalized coordinate form, it is necessary to map  ̈ onto ̈. ̈can represented as： In the equation， Substituting (23) into equation (22) allows us to write it as: Multiplying both sides of ( 24) by     gives the end effector system's dynamic equations in the generalized coordinate form: In the equation， Here,    and    denote the end effector's mass and generalized external moment column matrices in the form of the system's generalized coordinates, respectively.

The kinetic equations of open chain systems
For this paper's 4-RRR parallel robot system, we expressed all of the components' dynamic equations using the system's generalized coordinates.We can obtain the open chain system's dynamic equations by grouping all of the components:   (  )̈ +   (  , ̇)̇ +     =   . (26) In the equation，   (  ),   (  , ̇),   , and   respectively denote the open chain system's mass, centrifugal force/Coriolis force, stiffness, and generalized external moment column matrices.

The kinetic equations of closed chain systems
The system dynamic equations obtained above are in the form of an open chain; i.e.: the four branch chains are hypothetically severed from the end table.To obtain the complete closed chain system dynamic equations, we must add back the closed chain constraint to the equation.
We can write four sets of closed chain constraint equations (i=1,2,3,4) based on the closed chain vector relationship: In the equation,   denotes the angle between the line from P to   and the  −  axis;    denotes the length of the line from P to   .
Set  1 ,  2 ,  3 , and 4 in a matrix to obtain the closed chain constraint matrix equation for the system: The above equation does not explicitly include the time t for which the first-order derivative with respect to time is obtained.

The kinetic equations of closed chain systems written in the form of system independent generalized coordinates
The system dynamic equations obtained above contain the system generalized coordinate array,   , with a total of (11 + 8(  +   )) coordinate variables that are mutually coupled and nonindependent.To solve the system's dynamic equations, we must obtain them in the form of the system's independent generalized coordinates.Considering that this paper's parallel robot is in the horizontal plane and contains two translational and one rotational DOFs, we can utilize the position and posture array   = [      1 ]  at point P as the rigid system's independent coordinates.The mechanism consists of 8 flexible rods.
All the flexible rods' modal coordinates can be used as the flexible system's independent coordinates.Therefore, we can define the system's independent generalized coordinates as: (31) The first 3 terms in the equation are the end effector's rigid independent coordinates and the last 8 terms are the independent modal coordinates of the 8 flexible rods.
To represent the system's non-independent generalized coordinate array,   , as a systemindependent one, ̄, we divide all terms in   into 3 classes and rewrite them as: ] .
At this point, the system's non-independent generalized coordinate array's first-order derivative ̇ has been expressed by the first-order derivative ̄ of the system's independent generalized coordinate array.Plugging (36) in (29) gives: Since ̄ is a system independent generalized coordinate, the following consistently holds for the arbitrary ̄ in the system workspace: Hereby, we can define G as the system's closed-chain constrained Jacobi's zero-space matrix.

The dynamic response's numerical simulation
Based on the closed chain system's dynamics model in the form of system-independent generalized coordinates, we conducted a dynamics simulation for a 4-RRR parallel robot containing two parts when moving at high speed.One part is a large range of rigid motions while the other is a small range of flexible ones.Since the two influence each other, we utilized the ode23tb solver to solve the differential equations and simplify the model for efficiency.We took all flexible rods as the first order of modalities [25].To compare the robot's dynamic performance under different drive torque assignments, we first proposed different optimization schemes for the distribution of the rigid body model's active joint drive torques.We then inputted them into the rigid-flexible coupled 4-RRR parallel robot.The positive kinematic solution was then obtained to determine the system's dynamic response.
The 4-RRR parallel robot has 3 DOFs and 4 active joints.If a non-redundant drive scheme is used, we can select any 3 active joints as inputs, where the input torque solution is unique and deterministic.Otherwise, if a redundant drive scheme is used, then there is an infinite number of possibilities for the allocation of input moments.The latter can be broadly divided into three categories: the first is the minimum input moment, i.e.: the minimum Euclidean norm   = min√∑   for the four active moments.The second is the minimum input energy, i.e.: the minimum Euclidean norm   = min√∑ (  ) 2

𝑖=1
for the four active joints' input power.The third is the minimum power of some of the robot's internal forces.However, since the dynamics model's Lagrangian derivation only involves energy and not internal forces, we will not discuss this optimization scheme.
In this simulation, the planned path was the following: the end effector's center point P moves from point (0, 0) to the point (1.2mm, -17.5mm) in a spiral trajectory within 0-0.4s, while  1 = 2 3  remains constant.We chose the eccentric thread because the spiral trajectory has a gradually increasing velocity and acceleration on the one hand, which facilitates dynamic analysis.On the other hand, it permits us to study the system's response under such an eccentric trajectory as an eccentric thread.The trajectory, velocity, and acceleration of point P are shown in figures 4, 5, and 6, respectively.Table 1 demonstrates the parameters of a set of 4-RRR parallel robots.We assumed that all bars and end manipulators were of uniformly-textured aluminum and that all joints were lightweight ideal Volume 20, Issue 2, 3918-3943.bearings; i.e.: they are massless and frictionless, etc.Based on the previously-presented planning trajectory, we discuss the dynamic performance under 3 types of drive force distributions, namely mode a: a non-redundant drive (drive joints A1, A2, A3), mode b: a redundant drive with a minimum input moment, and mode c: a redundant drive with minimum input energy.According to the inverse rigid body dynamics model, their drive moments are calculated as shown in Figure 7.In MATLAB, the drive moments of these drive modes were used as prior inputs for the numerical calculations of the system's dynamic response.Figure shows the end elastic deformation of rods 1-8 in the redundant drive modes.Since the axial deformation is much smaller than the transverse one, we only calculated the prior one here.To improve the calculation efficiency, the first-order mode was intercepted for all flexible rods for calculation.(1) For the active rod, the elastic deformation of the flexible rod's end in the non-redundant driving mode had the largest vibration amplitude in each small range, which is about 2-3 times bigger compared to that of the other two redundant driving modes.(2) For the active rod, the flexible rod's end's elastic deformation was in the non-redundant drive mode.The vibration amplitude in the period between (0-0.35s) was slightly larger than that of the other two redundant drive modes; the vibration amplitude in the period between (0.35-0.4s) was much larger than that of the other two redundant drive modes.(3) For all three drive modes, the driven rod's deformation was much smaller than that of the drive rod.(4) For the middle 4 driven rods, the deformation magnitudes under the 3 drive force modes did not greatly differ.(5) The redundant drive with the smallest input torque had a slightly smaller maximum deformation for rods 2, 3, and 4 compared to the redundant drive with the smallest input energy.( 6) Since the robot's 1-4 rod was the active rod, its elastic deformation was affected by the driving force, demonstrating a trend fluctuation; the robot's 5-8 rod was the driven rod and its elastic deformation fluctuation was irregular.The deformation of the remaining flexible rods was equivalent.Overall, the deformation of the flexible rods in the two redundant drive modes was much smaller than that in the non-redundant drive mode.Compared to the two redundant drive modes, the redundant drive with the smallest input torque had a slightly smaller deformation of flexible rods compared to redundant drive with the smallest input energy.
Figure 9 shows point P's actual trajectory in different drive modes.Figure 10 shows its displacement in each direction, the desired displacement error, the end table's angle of rotation, and the desired angle error in different drive modes.The redundant drive with the lowest input torque and the one with the lowest energy input have relatively small deviations compared to the nonredundant drive.Figure 11 shows the error between the actual motion and desired angles of rods 1-8 in different drive modes.We can observe that the deformation of each flexible rod increases sharply from 0.2s in all three drive modes.Additionally, the degree of deformation is greater in the non-redundant drive mode compared to that in the remaining two redundant modes.
Figure 12 shows the actual kinematic angular velocities of rods 1-8 in different drive modes compared to the desired angular velocities.It is shown that for the active rod, the flexible rod in the non-redundant drive mode's angular velocity has the largest vibration amplitude in all small ranges, which is approximately 2-3 times greater than that of the other two redundant drive modes.Moreover, for the driven rod, the flexible rod in the non-redundant drive mode's angular velocities in some periods are slightly greater than the other two modes (for example, rod 5 at 0.35-0.4s).
Figure 13 shows the actual kinematic angular acceleration of rods 1-8 in different drive modes compared to the desired angular acceleration.The amplitude of vibration is slightly greater in some periods than in the other two modes (for example, rod 5 at 0.35-0.4s).In general, the flexible rod's angular acceleration vibration tendency is identical to that of the angular velocity, albeit at a higher frequency and with a greater amplitude.

Multi-Body software simulation validation
We based the simulation validation on a parametric model of a rigid-flexible coupled system for a 4-RRR parallel robot built on Adams.It was modeled using an approach similar to finite segment discretization, in which we discretized the eight rods by dividing each one into nine rigid body subsections, each linked by a massless torsion spring.As shown in Figure 14, the end-effector platform was treated as a rigid body, i.e.: it was not divided [18].We conducted the simulation in Adams using the numerical calculation's same trajectory and the simulation duration set to 0.4s.Since the workspace is in the horizontal plane, we set the gravity to 0.
Without loss of generality, the numerical calculation contained three torque distribution modes.The simulation was conducted with the lowest torque redundant drive.Figure 15 shows the trajectory displacement curve at the end console's center point P, which was obtained from the numerical calculation and the Adams simulation.Figure 16 shows the velocity profile at P. We can observe that the point's displacements and velocities obtained using the two methods are in good correspondence and have only some minor deviations.Figure 17 shows each rod's angular velocity profile obtained from the numerical calculation and the Adams simulation (taking rods 1-4 as an example).As observed, the angular velocities of each rod obtained by the two methods are in general agreement; their deviations are bound by 15%.These differences are due to the different means of discretization used by the two methods; the numerical calculation uses the assumed mode method, while the Adams simulation uses the finite section method, resulting in a slight variation in the number of cells.
(a) The x-direction velocity.
(b) The y-direction velocity.

Conclusion
In this paper, we designed and fabricated a planar 4-RRR parallel robot by combining the assumed mode and floating coordinate frame methods to discretize an arbitrary, flexible rod j.The modal truncation method was used to intercept the rod's first-order modal state.Then, assembling the two flexible rods of an arbitrary branch chain i, one can further obtain the dynamics model of a single chain.Additionally, we introduced the closed chain constraint to obtain its complete dynamic equations.We utilized the zero-space method that eliminates the generalized binding terms in the equations.There are seven functions in MATLAB that could be used to solve for the numerical solution of differential equations: ode45, ode23, ode113, ode15s, ode23s, ode23t, and ode23tb.Among them, ode23tb is realized by TR --BDF2 because it is based on the implicit Runge --Kutta formula, which is more adequate for solving rigid equations.Accordingly, we used the ode23td solver to calculate the established model and obtain the dynamic response.After that, we used the Adams software to verify the results.From that, the following could be concluded： (1) The mechanism has planar rotational degrees of freedom, allowing it to move in the desired position and posture without being limited to a translational movement like the conventional planar parallel robots.This allows for more complex movements.
(2) We proposed an effective modular rigid-flexible coupling modeling method based on the assumed mode one.We simplified it to enhance the accuracy and computational efficiency, which is of value for models with more flexible components and larger computational volumes.
(3) Different torque distribution schemes for the drive joints were used as feedforward inputs to obtain different dynamic performances.The results show the following: firstly, the redundant drive mode can suppress the rod's elastic deformation to a certain extent and effectively improve the robot end-effector's trajectory tracking accuracy, which is an important guideline for the practical control of future robots.Secondly, the minimum input torque redundant drive has slightly less flexible rod deformation compared to the minimum input energy redundant one.However, they are almost equivalent.
Furthermore, the dynamic model of the fully flexible link system proposed in this study is expected to aid in the optimization of lightweight design and dynamic performance of parallel robots, which are expected to be used in special processing scenarios and packaging industries.

Work outlook
(1) We did not consider the effect of joint gaps in this paper's modeling, which affects the actual model's accuracy.Thus, we need to consider this in our future research.
(2) In our future work, we should further investigate the expansion of the concept of parallel robots and their theoretical analysis method's applications to spatially parallel mechanisms.

Figure 1 .
Figure 1.A structural illustration of the 4-RRR parallel robot and its deformation under high-speed motion.

Figure 2 .
Figure 2. The deformation of the flexible rod j during ovements. ̇.

Figure 3 .
Figure 3.A schematic diagram of the 4-RRR parallel robot with the cut joint.

Figure 4 .
Figure 4.The planning trajectory at point P.

Figure 5 .
Figure 5. Planned velocity at point P.

Figure 6 .
Figure 6.Planned acceleration at point P.
(a) Non-redundant drives.(b) Minimum input torque redundant drive.(c)Minimum input power drive.

Figure 7 .
Figure 7.The drive torque for different drive modes.
(a) Elastic deformation at the end of rod 1.(b) Elastic deformation at the end of rod 2. (b) Elastic deformation at the of rod 3. (d) Elastic deformation at the end of rod 4.(e) Elastic deformation at the end of rod 5. (f) Elastic deformation at the end of rod 6.(g) Elastic deformation at end of rod 7.(h) Elastic deformation at end of rod 8.

Figure 8 .
Figure 8.The rod end's elastic deformation for different drive modes.

Figure 9 .
Figure 9. Point P's motion trajectory for different drive modes.

Figure 10 .
Figure 10.The displacement error in each direction at point P for different drive modes.

Figure 11 .
Figure 11.The angular error of 1-8 rods in different drive modes.

Figure 12 .
Figure 12.The angular velocity of rods 1-8 in different drive modes.

Figure 13 .
Figure 13.The angular acceleration of rods 1-8 in different drive modes.

Figure 16 .
Figure 16.The P-Point speed comparison.
−   −   , in which the flexible rod j is located with respect to the overall coordinate system  −  − , where   is the relative angle of rotation.The following equation can be obtained using the assumed mode method to discretize the flexible rod j: ,  (̄) denotes the   order axial modal shape function.,  (̄) denotes the transverse modal shape function of the   order.,  () denotes the axial modal time-varying coordinates.,  () denotes the transverse modal time-varying coordinates. and   denote axial and transverse modal orders, respectively.   () denotes the axial mode shape function matrices.   (̄) denotes the transverse mode shape function matrices.    () denotes the column matrix of the axial modal time-varying coordinates.

Table 1 .
Parameters of the 4-RRR parallel robot (in mm).