Next Article in Journal
Analysis of Natural Heat Dissipation Capacity of Hydraulic Tank and Relevant Influencing Factors
Next Article in Special Issue
New Frontiers in Parallel Robots
Previous Article in Journal
Damping Characteristics of Cantilever Beam with Obstacle Grid Particle Dampers
Previous Article in Special Issue
Consistent Solution Strategy for Static Equilibrium Workspace and Trajectory Planning of Under-Constrained Cable-Driven Parallel and Planar Hybrid Robots
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geometric Error Analysis of a 2UPR-RPU Over-Constrained Parallel Manipulator

1
School of Mechanical Engineering, Zhejiang Sci.-Tech. University, Hangzhou 310018, China
2
School of Mechanical Engineering, Hangzhou Dianzi University, Hangzhou 310018, China
*
Author to whom correspondence should be addressed.
Machines 2022, 10(11), 990; https://doi.org/10.3390/machines10110990
Submission received: 18 September 2022 / Revised: 25 October 2022 / Accepted: 27 October 2022 / Published: 29 October 2022
(This article belongs to the Special Issue New Frontiers in Parallel Robots)

Abstract

:
For a 2UPR-RPU over-constrained parallel manipulator, some geometric errors result in internal forces and deformations, which limit the improvement of the pose accuracy of the moving platform and shorten the service life of the manipulator. Analysis of these geometric errors is important for restricting them. In this study, an evaluation model is established to analyse the influence of geometric errors on the limbs’ comprehensive deformations for this manipulator. Firstly, the nominal inverse and actual forward kinematics are analysed according to the vector theory and the local product of the exponential formula. Secondly, the evaluation model of the limbs’ comprehensive deformations is established based on kinematics. Thirdly, 41 geometric errors causing internal forces and deformations are identified and the results are verified through simulations based on the evaluation model. Next, two global sensitivity indices are proposed and a sensitivity analysis is conducted using the Monte Carlo method throughout the reachable workspace of the manipulator. The results of the sensitivity analysis indicate that 10 geometric errors have no effects on the average angular comprehensive deformation and that the identified geometric errors have greater effects on the average linear comprehensive deformation. Therefore, the distribution of the global sensitivity index of the average linear comprehensive deformation is more meaningful for accuracy synthesis. Finally, simulations are performed to verify the results of sensitivity analysis.

1. Introduction

Parallel mechanisms with three DOFs have been successfully applied to hybrid serial–parallel machine tools, such as the well-known Eco-speed series, Tricept, and Exechon [1,2,3,4,5,6], owing to their high stiffness, large payload, and good dynamics. To achieve a simpler structure, Li et al. [3] designed a 2R1T (R denotes a rotational DOF, and T denotes a translational DOF) parallel mechanism named 2UPR-RPU. This mechanism is not only easier to control but also suitable for many operations along the surfaces. However, it is an over-constrained parallel mechanism with common constraints and over-constraints [7,8]. Some geometric errors in a manipulator based on this mechanism break the common constraints and over-constraints, resulting in internal forces and deformations. The internal forces and deformations not only limit the further improvement of the pose accuracy of the moving platform but also shorten the service life of the manipulator [9,10]. Therefore, it is necessary to restrict the internal-force-and-deformation-related geometric errors in the 2UPR-RPU parallel manipulator.
The accuracy design [11,12,13] can be applied to restrict geometric errors by determining the tolerances of the fabrication and assembly of machines. It consists of three components: error modelling [14,15,16], sensitivity analysis [17,18,19], and accuracy synthesis [20,21,22], where error modelling is the basis of sensitivity analysis and accuracy synthesis. Zhang et al. [13] applied the closed-loop vector and first-order perturbation methods to establish a geometric error model for a 2UPR-RPS over-constrained manipulator, and they identified the geometric errors that affected the pose errors of the moving platform. Zhang et al. [15] utilised the screw theory to establish a geometric error model for a 4RSR-SS over-constrained parallel tracking machine. With the use of the geometric error model, 53 geometric errors that had a significant influence on the pose errors of the moving platform were identified after sensitivity analysis. However, neither of the above two methods considers the deformations caused by internal forces in over-constrained parallel manipulators. Taking parameter uncertainties into account, Tang et al. [23] built a general interval kinetostatic model for a 2UPR-SPR over-constrained parallel machine to perform sensitivity analysis and tolerance allocation. To predict the pose errors of an over-constrained extendible support structure, Yu et al. [24] proposed a comprehensive model that simultaneously considered geometric errors, joint gaps, and link flexibility. In spite of good accuracy, these two models are complicated for the stiffness matrix needs to be derived and the stiffness coefficients of parts need to be obtained via finite element software.
Affected by geometric errors, the end poses of different limbs of a parallel manipulator should be theoretically inconsistent. However, they can be consistent in non-overconstrained parallel manipulators due to the existence of the moving platform and the motion deviations of passive joints. On this basis, a numerical iterative algorithm [25,26] was proposed to analyse the kinematics of non-overconstrained parallel manipulators with kinematic errors. Inspired by this algorithm, this study aims to establish an evaluation model based on kinematics to analyse the influence of geometric errors on the limbs’ comprehensive deformations for the 2UPR-RPU over-constrained parallel manipulator.
Based on the established evaluation model, sensitivity analysis can help reveal the influence of different internal-force-and-deformation-related geometric errors on the limbs’ comprehensive deformations. The interval analysis method and probabilistic method have been commonly used for sensitivity analysis of the moving platform’s pose error in literature. The interval analysis method treats geometric errors as interval variables and can get a balance between calculation speed and accuracy [11,18]. Treating geometric errors as random variables with a normal distribution, the probabilistic method can be divided into the Monte Carlo method and the probability modelling method. The Monte Carlo method calculates the moving platform’s pose errors according to the geometric error model and lots of random values of a geometric error [22,27]. It has good accuracy and low computational efficiency. The probability modelling method establishes an analytical model between the standard deviation of each geometric error and that of the moving platform’s pose error based on the geometric error model [28]. In spite of high computational efficiency, this method needs prior knowledge about probability distributions. Considering that the interval analysis method and probability modelling method are not suitable when the geometric error model is iterative, the Monte Carlo method is utilised to analyse the influence of geometric errors on the limbs’ comprehensive deformations in this paper.
The remainder of this paper is organised as follows. In Section 2, the 2UPR-RPU parallel mechanism is briefly introduced. Section 3 presents an analysis of the nominal inverse kinematics and actual forward kinematics. Section 4 establishes an evaluation model of the limbs’ comprehensive deformations caused by geometric errors. Based on the evaluation model, the internal-force-and-deformation-related geometric errors are identified and the results are verified through simulations in Section 5. In Section 6, two global sensitivity indices are proposed and sensitivity analysis is conducted. Simulations are also performed to verify the results of sensitivity analysis. Finally, the conclusions are drawn in Section 7.

2. 2UPR-RPU Parallel Mechanism

As shown in Figure 1, the 2UPR-RPU parallel mechanism mainly consists of a moving platform, two UPR limbs, one RPU limb, and one fixed base, where the moving platform and fixed base are represented by the isosceles right triangles ΔA1A2A3 and ΔB1B2B3. U, P, and R denote universal, prismatic, and revolute joints, respectively. B1, B2 and A3 are the centres of U, and A1, A2 and B3 are the centres of R. Because each universal joint is equivalent to two mutually perpendicular revolute joints, the UPR limb is equivalent to the RRPR limb, and the RPU limb is equivalent to the RPRR limb. The axis of the jth joint of the ith limb is denoted by si,j. A fixed coordinate system {oB; x, y, z} is established at the midpoint between B1 and B2, where x points from B2 to B1 and y points from oB to B3. Similarly, a moving coordinate system {oA; u, v, w} is also established, where u points from A2 to A1 and v points from oA to A3. The coordinate axes z and w are determined using the right-hand rule. For the 2UPR-RPU parallel mechanism, each limb exerts a force and a couple on the moving platform [8], where the two forces from the UPR limbs are parallel to v, and the three couples from the UPR and RPU limbs rotate around w. It is worth mentioning that the two forces parallel to v will lead to over-constraint, and the three couples rotating around w will lead to common constraints. Thus, the 2UPR-RPU parallel mechanism is an over-constrained parallel mechanism.

3. Kinematics

Inverse kinematics aims to calculate the displacements of all joints relative to their initial positions or angles according to a given target pose of the moving platform. Forward kinematics is the reverse operation of inverse kinematics. Inverse kinematics without considering geometric errors is called nominal inverse kinematics. In this section, the nominal inverse kinematics of actuated joints and passive joints is first introduced. Then, the actual forward kinematics of the limbs is derived.

3.1. Nominal Inverse Kinematics

The position and orientation of the moving platform shown in Figure 1 can be described by [ x y z ] T and [ α β γ ] T , respectively, where [ x y z ] T denotes the position coordinates of oA with respect to {oB; x, y, z} and [ α β γ ] T denotes the Euler angle with respect to z-x-v. Because only the translation motion along oBoA and the rotations around x and v can be achieved by the moving platform [8], [ z β γ ] T is sufficient to represent the poses. For a given target pose of the moving platform, the nominal displacements of actuated P-joints can be derived using the closed-loop vector method [10] as follows:
{ q 1 , 3 = B 1 A 1 B 1 A ˜ 1 q 2 , 3 = B 2 A 2 B 2 A ˜ 2 q 3 , 2 = B 3 A 3 B 3 A ˜ 3
where represents the Euclidean norm. A ˜ i denotes the initial position of Ai, which is determined by
{ B 1 A 1 = [ l A cos γ l B l A sin β sin γ z tan β l A cos β sin γ + z ] T B 2 A 2 = [ l A cos γ + l B l A sin β sin γ z tan β l A cos β sin γ + z ] T B 3 A 3 = [ 0 l A cos β z tan β l B l A s i n β + z ] T
where l A = A 1 A 2 / 2 and l B = B 1 B 2 / 2 .
For the first UPR limb, the first, second, and fourth joints are passive. The nominal displacement of the first joint can be expressed as
q 1 , 1 = β
The nominal displacement of the second joint can be expressed as
q 1 , 2 = arccos ( e 1 T B 1 A 1 B 1 A 1 ) arccos ( e 1 T B 1 A ˜ 1 B 1 A ˜ 1 )
where e1 is the unit vector along x.
The nominal displacement of the fourth joint can be expressed as
q 1 , 4 = arccos ( ( A 1 A ˜ 2 ) T ( B 1 A ˜ 1 ) A 1 A ˜ 2 B 1 A ˜ 1 ) arccos ( ( A 1 A 2 ) T ( B 1 A 1 ) A 1 A 2 B 1 A 1 )
Because the two UPR limbs are symmetrically distributed with respect to oAoB, we have
q 2 , 1 = β
q 2 , 2 = arccos ( e 1 T B 2 A 2 B 2 A 2 ) arccos ( e 1 T B 2 A ˜ 2 B 2 A ˜ 2 )
q 2 , 4 = arccos ( ( A 1 A 2 ) T ( B 2 A 2 ) A 1 A 2 B 2 A 2 ) arccos ( ( A 1 A ˜ 2 ) T ( B 2 A ˜ 2 ) A 1 A ˜ 2 B 2 A ˜ 2 )
Similarly, the nominal displacements of the first, third, and fourth joints of the RPU limb can be expressed as
q 3 , 1 = arccos ( e 2 T B 3 A 3 B 3 A 3 ) arccos ( e 2 T B 3 A ˜ 3 B 3 A ˜ 3 )
q 3 , 3 = arccos ( ( R e 2 ) T ( B 3 A 3 ) B 3 A 3 ) arccos ( ( R ˜ e 2 ) T ( B 3 A ˜ 3 ) B 3 A ˜ 3 )
q 3 , 4 = γ
Here, e2 is the unit vector along y, and R ˜ denotes the initial state of R, which is given as follows:
R = [ cos γ 0 sin γ sin β sin γ cos β sin β cos γ cos β sin γ sin β cos β cos γ ]

3.2. Actual Forward Kinematics

The nominal inverse kinematics described above does not consider geometric errors. However, geometric errors exist in the 2UPR-RPU parallel manipulator. In this section, the actual forward kinematics of the limbs in the manipulator is derived in detail.
As shown in Figure 2, four local coordinate systems {Fi,j; xi,j, yi,j, zi,j} are assigned to each limb to describe the geometric errors of the 2UPR-RPU parallel manipulator, where the initial pose of the moving platform is [ z 0 β 0 γ 0 ] T = [ 0.2 m 0 0 ] T under the home configuration. The coordinate systems {oB; x, y, z} and {oA; u, v, w} are identical to those in Figure 1. For brevity, we use {Fi,j} instead of {Fi,j; xi,j, yi,j, zi,j}. It is worth mentioning that this figure only shows xi,j and zi,j of the local coordinate systems, and yi,j can be determined according to the right-hand rule, which will not be illustrated in detail here. The definitions of the local coordinate systems for the two UPR limbs and the RPU limb are listed in Table 1, Table 2 and Table 3.
The end poses of the ith limb can be obtained from the local product of the exponential formula [25] as
g i ( q i ) = g i , 0 e ζ ^ i , 1 q i , 1 g i , 1 e ζ ^ i , 2 q i , 2 g i , 2 e ζ ^ i , 3 q i , 3 g i , 3 e ζ ^ i , 4 q i , 4 g i , 4 , i = 1 , 2 , 3
where gi denotes the homogeneous transformation matrix (HTM) of {oA; u, v, w} with respect to {oB; x, y, z} calculated using the ith limb. ζi,j denotes the screw coordinates of si,j with respect to {Fi,j}, which can be written as [25,26]
{ ζ i , j = [ 0 0 1 0 0 0 ] T   for   R   jo int ζ i , j = [ 0 0 0 0 0 1 ] T   for   P   joint
Here, e ζ ^ i , j q i , j denotes the exponential map from the Lie algebra se(3) to the special Euclidean group SE(3), which can be obtained using (A1)–(A4) in Appendix A. gi,j is the HTM between adjacent coordinate systems when the parallel manipulator is under the home configuration. To be more specific, gi,0 denotes the HTM of {Fi,1} with respect to {oB; x, y, z}; gi,4 denotes the HTM of {oA; u, v, w} with respect to {Fi,4}; when j ≠ 0 and j ≠ 4, gi,j is the HTM of {Fi,j+1} with respect to {Fi,j}. gi,j can be written as
{ g 1 , 0 = T r a n s ( x , l B + d ) R o t ( y , π / 2 ) g 1 , 1 = T r a n s ( z , d ) R o t ( x , π / 2 ) R o t ( z , q ˜ 1 , 2 π / 2 ) g 1 , 2 = R o t ( y , π / 2 ) g 1 , 3 = T r a n s ( z , q ˜ 1 , 3 ) R o t ( y , π / 2 ) R o t ( z , q ˜ 1 , 2 + π / 2 ) g 1 , 4 = T r a n s ( y , l A ) R o t ( y , π / 2 ) R o t ( z , π / 2 )
{ g 2 , 0 = T r a n s ( x , l B d ) R o t ( y , π / 2 ) g 2 , 1 = T r a n s ( z , d ) R o t ( x , π / 2 ) R o t ( z , q ˜ 2 , 2 π / 2 ) g 2 , 2 = R o t ( y , π / 2 ) g 2 , 3 = T r a n s ( z , q ˜ 2 , 3 ) R o t ( y , π / 2 ) R o t ( z , q ˜ 2 , 2 + π / 2 ) g 2 , 4 = T r a n s ( y , l A ) R o t ( y , π / 2 ) R o t ( z , π / 2 )
{ g 3 , 0 = T r a n s ( y , l B ) R o t ( y , π / 2 ) R o t ( z , q ˜ 3 , 1 π / 2 ) g 3 , 1 = R o t ( y , π / 2 ) g 3 , 2 = T r a n s ( z , q ˜ 3 , 2 ) R o t ( y , π / 2 ) R o t ( z , q ˜ 3 , 1 + π / 2 ) g 3 , 3 = T r a n s ( y , c ) R o t ( x , π / 2 ) g 3 , 4 = T r a n s ( z , c l A ) R o t ( y , π / 2 ) R o t ( z , π / 2 )
where Trans(x, lB) denotes the HTM that translates by lB along x, and Rot(y, π/2) denotes the HTM that rotates by π/2 around y. q ˜ 1 , 2 is the initial angle between x and B1A1, q ˜ 2 , 2 is the initial angle between B2B1 and B2A2, and q ˜ 3 , 1 is the initial angle between B3oB and B3A3, which can be expressed as
q ˜ 1 , 2 = arccos ( e 1 T B 1 A ˜ 1 B 1 A ˜ 1 )
q ˜ 2 , 2 = arccos ( e 1 T B 2 A ˜ 2 B 2 A ˜ 2 )
q ˜ 3 , 1 = arccos ( e 2 T B 3 A ˜ 3 B 3 A ˜ 3 )
In contrast to q ˜ 1 , 2 , q ˜ 2 , 2 , and q ˜ 3 , 1 , q ˜ 1 , 3 , q ˜ 2 , 3 , and q ˜ 3 , 2 are the initial positions of the actuated P-joints, and we have
q ˜ 1 , 3 = B 1 A ˜ 1
q ˜ 2 , 3 = B 2 A ˜ 2
q ˜ 3 , 2 = B 3 A ˜ 3
The linear errors δi,j of {Fi,j+1} along xi,j, yi,j, and zi,j can be expressed as follows:
δ i , j = [ δ i , j x δ i , j y δ i , j z ] T , i = 1 ,   2 ,   3   and   j = 0 ,   ,   4
In addition to linear errors, angular errors also exist. The angular errors εi,j of {Fi,j+1} around xi,j, yi,j, and zi,j can be expressed as follows:
ε i , j = [ ε i , j x ε i , j y ε i , j z ] T , i = 1 ,   2 ,   3   and   j = 0 ,   ,   4
where δi,0 and εi,0 denote the linear and angular errors of {Fi,1} with respect to {oB; x, y, z}, respectively. δi,4 and εi,4 denote the linear and angular errors of {Fi,4} with respect to {oA; u, v, w}, respectively. Among the 90 error parameters, ε 1 , 0 x , ε 1 , 1 y , δ 1 , 3 z , ε 1 , 3 x , ε 2 , 0 x , ε 2 , 1 y , δ 2 , 3 z , ε 2 , 3 x , ε 3 , 0 x , δ 3 , 2 z , ε 3 , 2 x , and ε 3 , 3 y represent the initial displacement errors of the 12 joints. In addition, the values of δ 1 , 2 x , δ 2 , 2 x , ε 1 , 4 y , ε 2 , 4 y , δ 3 , 1 x , ε 3 , 3 z , and ε 3 , 4 y are zeros since the definitions of local coordinate systems. Therefore, the rest 71 error parameters represent the linear and angular geometric errors.
Setting the values of error parameters other than geometric errors to zeros, the HTM of the geometric errors between adjacent coordinate systems can be written as
Δ g i , j = [ e ε ^ i , j δ i , j 0 1 × 3 1 ] , i = 1 ,   2 ,   3   and   j = 0 ,   ,   4
where e ε ^ i , j denotes the exponential map from the Lie algebra so(3) to the special orthogonal group SO(3), which can be determined using (A3) and (A5) in Appendix A.
The end poses of the ith limb that include sthe linear and angular geometric errors can then be obtained as follows:
g i g e ( q i ) = Δ g i , 0 g i , 0 e ζ ^ i , 1 q i , 1 Δ g i , 1 g i , 1 e ζ ^ i , 2 q i , 2 Δ g i , 2 g i , 2 e ζ ^ i , 3 q i , 3 Δ g i , 3 g i , 3 e ζ ^ i , 4 q i , 4 g i , 4 Δ g i , 4 1 , i = 1 , 2 , 3
which can be rewritten as
g i g e ( q i ) = g i , 0 g e e ζ ^ i , 1 q i , 1 g i , 1 g e e ζ ^ i , 2 q i , 2 g i , 2 g e e ζ ^ i , 3 q i , 3 g i , 3 g e e ζ ^ i , 4 q i , 4 g i , 4 g e , i = 1 , 2 , 3

4. Evaluation Model of Deformations

As mentioned previously, the 2UPR-RPU parallel manipulator is over-constrained. Theoretically, the end poses of any two limbs can also be consistent with each other through the motion deviations of passive joints when the internal-force-and-deformation-related geometric errors are zero, which can be expressed as
g i g e + Δ g i g e = g k g e + Δ g k g e
where Δ g i g e and Δ g k g e denote the end-pose deviations of the ith and kth limbs caused by the motion deviations of passive joints, respectively. The end-pose deviation between the ith and kth limbs can be written as [25,26]
Δ μ k , i = { log [ g k g e ( g i g e ) 1 ] }
where log [ ] stands for the logarithmic operation from SE(3) to se(3), and it can be obtained using (A6) and (A7) in Appendix A. represents the reverse operation of (A1). The end-pose deviation can be rewritten in screw form as follows:
Δ μ k , i = Δ μ i Δ μ k
where the screws Δμi and Δμk denote the end-pose deviations of the ith and kth limbs originating from the motion deviations of passive joints, respectively. Take Δμi as an example. Taking the partial differential of (28) with respect to the displacements of the passive joints, Δμi can be expressed as follows:
Δ μ i = Ψ i Φ i Δ q i , i = 1 , 2 , 3
where Δqi denotes the motion deviation of the passive joints of the ith limb.
When i = 1 and i = 2, we have
Δ q i = [ Δ q i , 1 Δ q i , 2 Δ q i , 4 ] T
and when i = 3, we have
Δ q i = [ Δ q i , 1 Δ q i , 3 Δ q i , 4 ] T
For the coefficient matrices Ψi and Φi, when i = 1 and i = 2, we obtain
Ψ i = [ I 6 Ad ( e ξ ^ i , 1 q i , 1 ) Ad ( e ξ ^ i , 1 q i , 1 e ξ ^ i , 2 q i , 2 e ξ ^ i , 3 q i , 3 ) ] 6 × 18
Φ i = Blkdiag ( ξ i , 1 , ξ i , 2 , ξ i , 4 ) 18 × 3
When i = 3, we obtain
Ψ i = [ I 6 Ad ( e ξ ^ i , 1 q i , 1 e ξ ^ i , 2 q i , 2 ) Ad ( e ξ ^ i , 1 q i , 1 e ξ ^ i , 2 q i , 2 e ξ ^ i , 3 q i , 3 ) ] 6 × 18
Φ i = Blkdiag ( ξ i , 1 , ξ i , 3 , ξ i , 4 ) 18 × 3
where I6 is an identity matrix of order six. Ad ( ) is an adjoint representation of SE(3) and is given in (A8) in Appendix A. Blkdiag ( ) denotes a block-diagonal matrix. ξi,j denotes the screw coordinates of si,j with respect to {oB; x, y, z}, which can be written as follows [25]:
ξ i , j = Ad ( g i , 0 g e , g i , 1 g e , , g i , j 1 g e ) ζ i , j
Combining (31) with (32) yields
[ Δ μ 2 , 1 Δ μ 3 , 2 Δ μ 1 , 3 ] = [ Ψ 1 Φ 1 Ψ 2 Φ 2 0 0 Ψ 2 Φ 2 Ψ 3 Φ 3 Ψ 1 Φ 1 0 Ψ 3 Φ 3 ] [ Δ q 1 Δ q 2 Δ q 3 ]
Let
Δ μ = [ Δ μ 2 , 1 T Δ μ 3 , 2 T Δ μ 1 , 3 T ] T 18 × 1
Δ q = [ Δ q 1 T Δ q 2 T Δ q 3 T ] T 9 × 1
J = [ Ψ 1 Φ 1 Ψ 2 Φ 2 0 0 Ψ 2 Φ 2 Ψ 3 Φ 3 Ψ 1 Φ 1 0 Ψ 3 Φ 3 ] 18 × 9
Note that when ∆μ is obtained using (30) and (41), the motion deviations of passive joints can be calculated as
Δ q = ( J T J ) 1 J T Δ μ = J c Δ μ
Based on the above work, an iterative model can be proposed to evaluate the deformations caused by geometric errors of the 2UPR-RPU over-constrained manipulator. The detailed processes are described below.
As shown in Figure 3, the proposed evaluation model mainly includes the following steps. Firstly, a target pose of the moving platform is input, and specified values are assigned to some of the 71 geometric errors; secondly, the nominal displacements of all joints are calculated based on the inverse kinematics; thirdly, the displacements of the passive joints are iteratively updated starting with the nominal values and the end condition is given as the maximum number of iterations or the target value of the infinity norm Δ μ j Δ μ j 1 ; finally, the latest end-pose deviation Δ μ j for the target pose is output. When the internal-force-and-deformation-related geometric errors are not all zeros, the end poses of the limbs cannot be consistent without deformations. Therefore, the latest Δ μ j and indices based on it can be used to indirectly evaluate the limbs’ comprehensive deformations caused by geometric errors of the 2UPR-RPU over-constrained manipulator.
Considering that a large amount of matrix calculation is included in the proposed evaluation model, MATLAB is used for programming in Section 5 and Section 6.

5. Geometric Error Identification

Finding the internal-force-and-deformation-related geometric errors is the basis of sensitivity analysis. In this section, the reachable workspace of the 2UPR-RPU parallel manipulator is described. For geometric error identification and verification, 692 and 1738 target poses are selected in the reachable workspace. Subsequently, internal-force-and-deformation-related geometric errors in the manipulator are identified based on the proposed evaluation model and an evaluation index. Finally, simulations are conducted to verify the correctness of the identification results.

5.1. Identification Analysis

The structural parameters of the 2UPR-RPU parallel manipulator are presented in Table 4. Using the space search method [29], the reachable workspace of the manipulator can be obtained. The search results are shown in Figure 4. Because the end poses at the boundaries of the reachable workspace are more sensitive to geometric errors, the 692 target poses shown in Figure 5 are uniformly selected for geometric error identification. To identify the internal-force-and-deformation-related geometric errors, the evaluation index of the maximum comprehensive deformation of a limb can be written as
Δ μ max = max ( Δ μ 2 , 1 j , Δ μ 3 , 2 j , Δ μ 1 , 3 j )
Based on Scheme I and (45), 692 Δ μ max s can be calculated for the selected target poses of the moving platform. If the 692 Δ μ max s are not close to zero, it means that there are internal-force-and-deformation-related geometric errors among the geometric errors that were assigned specified values. Without loss of generality, three groups of specified values for geometric errors are given, as listed in Table 5. The maximum iteration number λ and specified tolerance τ in Scheme I are set to 50 and 10−15, respectively.
Taking group 1 as an example, the detailed processes are described as follows: (1) δ 1 , 0 x is set to 0.005 m, and the remaining geometric errors are set to 0. (2) 692 Δ μ max s are calculated according to Scheme I and (45). (3) If the number of Δ μ max s that are smaller than 10−15 is less than 657 (≈95% of 692), then δ 1 , 0 x is referred to as an internal-force-and-deformation-related geometric error. After repeating the above steps for the 71 geometric errors, 39 internal-force-and-deformation-related geometric errors were initially identified and are listed in Table 6. In the table, “✓” denotes the internal-force-and-deformation-related geometric error; “–” denotes the error parameter that is not a geometric error.
Some geometric errors between any two adjacent coordinate systems in a limb may be linearly dependent. Therefore, it is necessary to analyse geometric errors simultaneously. Based on the results in Table 6, the set of the six error parameters, [ δ i , j x , δ i , j y , δ i , j z , ε i , j x , ε i , j y , ε i , j z ] , are regarded as one unit. Take [ δ 1 , 1 x , δ 1 , 1 y , δ 1 , 1 z , ε 1 , 1 x , ε 1 , 1 y , ε 1 , 1 z ] as an example. [ δ 1 , 1 x , δ 1 , 1 y , δ 1 , 1 z , ε 1 , 1 x , ε 1 , 1 y , ε 1 , 1 z ] is set to [0.005 m, 0, 0.005 m, 0, 0, 0.005 rad], and the remaining units are set to [0,0,0,0,0,0]. Then, 692 Δ μ max s are calculated according to Scheme I and (45). If the number of Δ μ max s that are smaller than 10−15 is less than 657, the internal-force-and-deformation-related geometric errors are included in δ 1 , 1 x , δ 1 , 1 z , and ε 1 , 1 z . Then, δ 1 , 1 x , δ 1 , 1 z , and ε 1 , 1 z are set to 0 in turn, and the remaining units are unchanged. The Δ μ max s are recalculated. If the number of Δ μ max s that decrease significantly is greater than 656, it is determined that the geometric error, which is set as 0, will cause internal forces and deformations. These steps were repeated for each error unit and the results are listed in Table 7. The identification results for groups 2 and 3 in Table 5 are the same as those shown in Table 7. The results demonstrate that there are 41 internal-force-and-deformation-related geometric errors, where the number of angular geometric errors is greater than that of linear geometric errors. In addition, the internal-force-and-deformation-related geometric errors of the first UPR limb are the same as those of the second UPR limb because of the symmetric distribution of the two limbs. For the RPU limb, the geometric errors that cause internal forces and deformations are angular geometric errors.

5.2. Simulation Analysis

To validate the correctness of the identified results listed in Table 7, three groups of numerical simulations were conducted using 1738 target poses of the 2UPR-RPU parallel manipulator, as shown in Figure 6. It is assumed that geometric errors are normally distributed with zero means [19,22]. Three groups of standard deviations are listed in Table 8. In the simulation, the internal-force-and-deformation-related geometric errors identified in Table 7 were set to 0, and the remaining 30 geometric errors were assigned random values generated by randn function using the standard deviations of δi,j and εi,j listed in Table 8. Then, according to Scheme I and (45), 1738 Δ μ max s were calculated for each group. The simulation results are shown in Figure 7. It can be seen that Δ μ max s of Group 1, Group 2, and Group 3, are all smaller than 10−15. This demonstrates that the internal-force-and-deformation-related geometric errors identified in Section 5.1 are correct.

6. Sensitivity Analysis

Sensitivity analysis can help reveal the influence of different internal-force-and-deformation-related geometric errors on the limbs’ comprehensive deformations. Since Δ μ j is calculated iteratively in Scheme I, the Monte Carlo method [22] is utilised to conduct sensitivity analysis in this section. Two global sensitivity indices are proposed and the results of sensitivity analysis are verified through simulations.

6.1. Sensitivity Indices

According to (41), Δ μ j consists of Δ μ 2 , 1 j , Δ μ 3 , 2 j , and Δ μ 1 , 3 j , which can be written as
Δ μ k , i j = [ ω k , i , 1 j ω k , i , 2 j ω k , i , 3 j v k , i , 1 j v k , i , 2 j v k , i , 3 j ] T
The end-orientation and end-position volumetric deviations between any two limbs are
Δ ω k , i j = ( ω k , i , 1 j ) 2 + ( ω k , i , 2 j ) 2 + ( ω k , i , 3 j ) 2
and
Δ v k , i j = ( v k , i , 1 j ) 2 + ( v k , i , 2 j ) 2 + ( v k , i , 3 j ) 2
Then, the evaluation indices of the average angular and linear comprehensive deformations of the three limbs can be written as
Δ ω a j = Δ ω 2 , 1 j + Δ ω 3 , 2 j + Δ ω 1 , 3 j 3
and
Δ v a j = Δ v 2 , 1 j + Δ v 3 , 2 j + Δ v 1 , 3 j 3
Under the condition that geometric errors are normally distributed with zero means, the sensitivity indices of the average angular and linear comprehensive deformations with respect to a geometric error can be written as
μ ω , p = σ ( Δ ω a , p j ) σ ( G e p ) ,   p = 1 ,   2 ,   ,   25  
and
μ v , p = σ ( Δ v a , p j ) σ ( G e p ) ,   p = 1 ,   2 ,   ,   25  
where σ ( ) denotes the standard deviation and can be calculated by std function. Geps are the internal-force-and-deformation-related geometric errors of the first and third limbs. Because of the symmetric distribution of the first and second limbs, the internal-force-and-deformation-related geometric errors of the second limb are not considered. Generally, the values of sensitivity indices vary with different target poses of the moving platform. Hence, m target poses should be chosen and the global sensitivity indices can be written as [16]
μ ω , p g = i = 1 m μ ω , p , i m + σ ( μ ω , p )
and
μ v , p g = i = 1 m μ v , p , i m + σ ( μ v , p )

6.2. Sensitivity Analysis

Based on the equations in Section 6.1 and Scheme I, the detailed processes to calculate the two global sensitivity indices with respect to Gep are described as follows: (1) Set σ ( G e p ) to 1 mm (0.001 m) or 1° (π/180 rad) for linear or angular geometric error. And the other 40 internal-force-and-deformation-related geometric errors are set to 0. In addition, the rest 30 linear or angular geometric errors are set to 1 mm or 1°. (2) Assign 1000 random values that obey the normal distribution to Gep and calculate 1000 Δ ω a , p j s and Δ v a , p j s. (3) Calculate σ ( Δ ω a , p j ) , μ ω , p , σ ( Δ v a , p j ) , and μ v , p for a target pose of the moving platform. (4) Repeat the above steps for m target poses and calculate the global sensitivity indices μ ω , p g and μ v , p g .
In order to improve the computational efficiency, 158 of the 1738 target poses shown in Figure 6 were selected uniformly to perform the above steps for each Gep. The global sensitivity indices of the average angular and linear comprehensive deformations with respect to Geps are shown in Figure 8 and Figure 9, respectively. It can be seen that the values of μ ω , p g with respect to Ge5( δ 1 , 1 y ), Ge7( ε 1 , 1 z ), Ge8( δ 1 , 2 z ), Ge11( δ 1 , 3 x ), and Ge14( δ 1 , 4 y ), are zero. This indicates that the corresponding geometric errors have no effects on the average angular comprehensive deformation. It is worth mentioning that δ 2 , 1 y , ε 2 , 1 z , δ 2 , 2 z , δ 2 , 3 x , and δ 2 , 4 y , have also no effects on the average angular comprehensive deformation due to the symmetric distribution of the first and second limbs. Comparing Figure 8 with Figure 9, it can also be found that the value of μ v , p g is larger than that of μ ω , p g for each Gep. This demonstrates that the internal-force-and-deformation-related geometric errors have greater effects on the average linear comprehensive deformation. Thus, the distribution of the global sensitivity index μ v , p g is more useful for accuracy synthesis. According to Figure 9, Geps can be sorted in descending order as follows: Ge16( ε 1 , 4 z ), Ge25( ε 3 , 4 z ), Ge15( ε 1 , 4 x ), Ge4( ε 1 , 0 z ), Ge24( ε 3 , 4 x ), Ge12( ε 1 , 3 y ), Ge6( ε 1 , 1 x ), Ge3( ε 1 , 0 y ), Ge13( ε 1 , 3 z ), Ge9( ε 1 , 2 x ), Ge18( ε 3 , 0 z ), Ge23( ε 3 , 3 x ), Ge17( ε 3 , 0 y ), Ge19( ε 3 , 1 x ), Ge22( ε 3 , 2 z ), Ge10( ε 1 , 2 y ), Ge20( ε 3 , 1 y ), Ge21( ε 3 , 2 y ), Ge8( δ 1 , 2 z ), Ge14( δ 1 , 4 y ), Ge5( δ 1 , 1 y ), Ge11( δ 1 , 3 x ), Ge1( δ 1 , 0 y ), Ge2( δ 1 , 0 z ), and Ge7( ε 1 , 1 z ). In order to lower the cost of fabrication and assembly, the allowable range of geometric errors should be larger and larger from Ge16 to Ge7.

6.3. Verification

6.3.1. Average Angular Comprehensive Deformation

As shown in Table 9, three groups of specified values for geometric errors are given. For each group, 1738 Δ ω a j s were calculated according to Scheme I and using the target poses shown in Figure 6. The maximum and average values of Δ ω a j are listed in Table 9. It can be seen that both the maximum and average values of Δ ω a j do not change from Group 1 to Group 3. This indicates that Ge5, Ge7, Ge8, Ge11, Ge14, δ 2 , 1 y , ε 2 , 1 z , δ 2 , 2 z , δ 2 , 3 x , and δ 2 , 4 y , have no effects on the average angular comprehensive deformation.

6.3.2. Average Linear Comprehensive Deformation

As shown in Table 10, Ge7( ε 1 , 1 z ), which has the smallest effect on the average linear comprehensive deformation, is set to 1°. Then, the other Geps are set to μ v , 7 g / μ v , p g mm or °. It is worth mentioning that the corresponding internal-force-and-deformation-related geometric errors of the second limb are assigned the same values as the first limb due to the symmetric distribution of the two UPR limbs. The remaining 30 geometric errors are set to 0.1 mm or °. According to Scheme I and using the target poses shown in Figure 6, 1738 Δ ω a j s and Δ v a j s were calculated. For comparison, the internal-force-and-deformation-related linear and angular geometric errors are set to their average values, 0.0209 mm and 0.0743 °, respectively, while the values of the remaining 30 geometric errors are unchanged. After recalculation, the maximum and average values of Δ ω a j and Δ v a j are listed in Table 11. It can be seen that both the maximum and average values of Δ v a j are larger than that of Δ ω a j for each group. This indicates that the internal-force-and-deformation-related geometric errors have greater effects on the average linear comprehensive deformation. It can also be found that from Group 2 to Group 1, the maximum and average values of Δ ω a j and Δ v a j decreased by 84%, 83%, 91%, and 89%, respectively. This demonstrates that at the same cost, restricting the values of geometric errors according to the sensitivity analysis results of the average linear comprehensive deformation can dramatically decrease the average angular and linear comprehensive deformations. Furthermore, it indirectly verifies the sensitivity analysis results of the average linear comprehensive deformation.

7. Conclusions

This paper deals with error modelling and sensitivity analysis of geometric errors that cause internal forces and deformations in the 2UPR-RPU over-constrained parallel manipulator. Conclusions are drawn as follows:
(1) The nominal inverse kinematics and actual forward kinematics of the over-constrained parallel manipulator are analysed according to the vector theory and the local product of the exponential formula. On this basis, an iterative model is established to indirectly evaluate the limbs’ comprehensive deformations caused by geometric errors.
(2) Based on the iterative evaluation model, the maximum Euclidean norm of the end-pose deviations of limbs is defined as an evaluation index of the maximum comprehensive deformation of a limb. Programming with MATLAB, 41 internal-force-and-deformation-related geometric errors are identified. Among the 41 geometric errors, the number of angular geometric errors is greater than that of linear geometric errors; the geometric errors of the first UPR limb are the same as those of the second UPR limb; the geometric errors of the RPU limb are all angular geometric errors. The correctness of the identification results is verified through simulations under the condition that geometric errors are normally distributed with zero means.
(3) The global sensitivity indices of the average angular and linear comprehensive deformations with respect to internal-force-and-deformation-related geometric errors are proposed and calculated based on the Monte Carlo method. The results of sensitivity analysis demonstrate that δ 1 , 1 y , ε 1 , 1 z , δ 1 , 2 z , δ 1 , 3 x , δ 1 , 4 y , δ 2 , 1 y , ε 2 , 1 z , δ 2 , 2 z , δ 2 , 3 x , and δ 2 , 4 y , have no effects on the average angular comprehensive deformation. Furthermore, the internal-force-and-deformation-related geometric errors have greater effects on the average linear comprehensive deformation. Therefore, the distribution of the global sensitivity index of the average linear comprehensive deformation with respect to geometric errors is more meaningful for accuracy synthesis. Finally, the results of sensitivity analysis are verified through simulations.
Based on the work presented in this paper, we will establish a model for accuracy synthesis and determine the tolerances of the fabrication and assembly of the manipulator in the future.

Author Contributions

Conceptualization, X.D. and J.Z.; methodology, X.D.; software, B.W.; validation, X.D., B.W. and J.Z.; formal analysis, B.W.; investigation, X.D.; resources, X.D.; data curation, B.W.; writing—original draft preparation, X.D.; writing—review and editing, J.Z.; visualization, X.D.; supervision, J.Z.; project administration, J.Z.; funding acquisition, J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by [National Natural Science Foundation of China] grant number [52275469] and [Open Fund of State Key Laboratory of Robotics and System] grant number [SKLRS-2021-KF-08].

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that there are no conflict of interest.

Appendix A. Lie Groups and Lie Algebras

Some equations about Lie groups and Lie algebras [25,30,31] are introduced here so that this work can be clearly understood. For a screw ζ = [ ω T ν T ] T , the operation denotes
ζ ^ = [ ω ^ ν 0 1 × 3 0 ] s e ( 3 )
The exponential map from the Lie algebra se(3) to the special Euclidean group SE(3) can be determined by
e ζ ^ q = [ I 3 + sin q ω ^ + ( 1 cos q ) ω ^ 2 V ν 0 1 × 3 1 ] S E ( 3 )
where I3 is an identity matrix of order three. ω ^ and V are expressed as
ω ^ = [ 0 ω 3 ω 2 ω 3 0 ω 1 ω 2 ω 1 0 ] s o ( 3 )
V = q I 3 + ( 1 cos q ) ω ^ + ( q sin q ) ω ^ 2
The exponential map from the Lie algebra so(3) to the special orthogonal group SO(3) can be determined by
e ω ^ = I 3 + sin ω ω ω ^ + 1 cos ω ω 2 ω ^ 2
For a HTM g S E ( 3 ) , the Lie algebra se(3) can be obtained as
log ( g ) = 1 8 csc 3 θ 2 sec θ 2 [ θ cos 2 θ sin θ θ cos θ 2 θ cos 2 θ + sin θ + sin 2 θ 2 θ cos θ + θ cos 2 θ sin θ sin 2 θ θ cos θ + sin θ ] T [ I 4 g g 2 g 3 ]
where
θ = arccos ( Tr ( g ) 2 2 ) , θ ( π , π )
The adjoint representation of g can be written as
Ad ( g ) = Ad ( [ R t 0 1 × 3 1 ] ) = [ R 0 3 × 3 t ^ R R ] 6 × 6

References

  1. Tang, T.F.; Fang, H.L.; Zhang, J. Hierarchical design, laboratory prototype fabrication and machining tests of a novel 5-axis hybrid serial-parallel kinematic machine tool. Robot. Comput.-Integr. Manuf. 2020, 64, 101944. [Google Scholar] [CrossRef]
  2. Xie, F.G.; Liu, X.J.; Li, T.M. Type synthesis and typical application of 1T2R-type parallel robotic mechanisms. Math. Probl. Eng. 2013, 2013, 206181. [Google Scholar] [CrossRef] [Green Version]
  3. Li, Q.C.; Hervé, J.M. Type synthesis of 3-DOF RPR-equivalent parallel mechanisms. IEEE Trans. Robot. 2014, 30, 1333–1343. [Google Scholar] [CrossRef]
  4. Li, Q.; Xu, L.; Chen, Q.; Ye, W. New family of RPR-equivalent parallel mechanisms: Design and application. Chin. J. Mech. Eng. 2017, 30, 217–221. [Google Scholar] [CrossRef]
  5. Huang, T.; Dong, C.; Liu, H.; Sun, T.; Chetwynd, D.G. A simple and visually orientated approach for type synthesis of overconstrained 1T2R parallel mechanisms. Robotica 2019, 37, 1161–1173. [Google Scholar] [CrossRef] [Green Version]
  6. Tang, T.; Fang, H.; Luo, H.; Song, Y.; Zhang, J. Type synthesis, unified kinematic analysis and prototype validation of a family of Exechon inspired parallel mechanisms for 5-axis hybrid kinematic machine tools. Robot. Comput.-Integr. Manuf. 2021, 72, 102181. [Google Scholar] [CrossRef]
  7. Zhao, Y.; Xu, Y.; Yao, J.; Jin, L. A force analysis method for overconstrained parallel mechanisms. China Mech. Eng. 2014, 25, 711–717. [Google Scholar] [CrossRef]
  8. Chai, X.X.; Xiang, J.N.; Li, Q.C. Singularity analysis of a 2-UPR-RPU parallel mechanism. J. Mech. Eng. Chin. Ed. 2015, 51, 144–151. [Google Scholar] [CrossRef] [Green Version]
  9. Jiang, Y.; Li, T.; Wang, L.; Chen, F. Kinematic error modeling and identification of the over-constrained parallel kinematic machine. Robot. Comput.-Integr. Manuf. 2018, 49, 105–119. [Google Scholar] [CrossRef]
  10. He, L.; Li, Q.; Zhu, X.; Wu, C.-Y. Kinematic calibration of a three degrees-of-freedom parallel manipulator with a laser tracker. J. Dyn. Syst. Meas. Control 2019, 141, 031009. [Google Scholar] [CrossRef]
  11. Song, Y.M.; Zhai, Y.P.; Sun, T. Interval analysis based accuracy design of a 3-DOF rotational parallel mechanism. J. Beijing Univ. Technol. 2015, 41, 1620–1626,1755. [Google Scholar]
  12. Ni, Y.; Zhang, B.; Sun, Y.; Zhang, Y. Accuracy analysis and design of A3 parallel spindle head. Chin. J. Mech. Eng. 2016, 29, 239–249. [Google Scholar] [CrossRef]
  13. Zhang, J.; Chi, C.C.; Jiang, S.J. Accuracy design of 2UPR-RPS parallel mechanism. Trans. Chin. Soc. Agric. Eng. 2021, 52, 411–420. [Google Scholar] [CrossRef]
  14. Huang, T.; Whitehouse, D.J.; Chetwynd, D.G. A unified error model for tolerance design, assembly and error compensation of 3-DOF parallel kinematic machines with parallelogram struts. CIRP Ann. 2002, 51, 297–301. [Google Scholar] [CrossRef]
  15. Zhang, J.T.; Lian, B.B.; Song, Y.M. Geometric error analysis of an over-constrained parallel tracking mechanism using the screw theory. Chin. J. Aeronaut. 2019, 32, 1541–1554. [Google Scholar] [CrossRef]
  16. Luo, X.; Xie, F.; Liu, X.-J.; Li, J. Error modeling and sensitivity analysis of a novel 5-degree-of-freedom parallel kinematic machine tool. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 2019, 233, 1637–1652. [Google Scholar] [CrossRef]
  17. Wang, W.; Yun, C. Orthogonal experimental design to synthesize the accuracy of robotic mechanism. J. Mech. Eng. Chin. Ed. 2009, 45, 18–24. [Google Scholar] [CrossRef]
  18. Yao, R.; Zhu, W.B.; Huang, P. Accuracy analysis of stewart platform based on interval analysis method. Chin. J. Mech. Eng. 2013, 26, 29–34. [Google Scholar] [CrossRef]
  19. Wu, M.; Yue, X.; Chen, W.; Nie, Q.; Zhang, Y. Accuracy analysis and synthesis of asymmetric parallel mechanism based on Sobol-QMC. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2020, 234, 4200–4214. [Google Scholar] [CrossRef]
  20. Li, X.Y.; Chen, W.Y.; Han, X.G. Accuracy analysis and synthesis of 3-RPS parallel machine based on orthogonal design. J. Beijing Univ. Aeronaut. Astronaut. 2011, 37, 979–984. [Google Scholar] [CrossRef]
  21. Tai-Ke, Y.; Xi, Z.; Feng, Z.; Li-Min, Z.; Yong, W. Accuracy synthesis of a 3-RPS parallel robot based on manufacturing costs. In Proceedings of the 31st Chinese Control Conference, Hefei, China, 25–27 July 2012; pp. 5168–5172. [Google Scholar]
  22. Liu, H.T.; Pan, Q.; Yin, F.W.; Dong, C.L. Accuracy synthesis of the TriMule hybrid robot. J. Tianjin Univ. Sci. Technol. 2019, 52, 1245–1254. [Google Scholar]
  23. Tang, T.; Luo, H.; Song, Y.; Fang, H.; Zhang, J. Chebyshev inclusion function based interval kinetostatic modeling and parameter sensitivity analysis for Exechon-like parallel kinematic machines with parameter uncertainties. Mech. Mach. Theory 2021, 157, 104209. [Google Scholar] [CrossRef]
  24. Yu, D.; Zhao, Q.; Guo, J.; Chen, F.; Hong, J. Accuracy analysis of spatial overconstrained extendible support structures considering geometric errors, joint clearances and link flexibility. Aerosp. Sci. Technol. 2021, 119, 107098. [Google Scholar] [CrossRef]
  25. Chen, G.; Kong, L.; Li, Q.; Wang, H.; Lin, Z. Complete, minimal and continuous error models for the kinematic calibration of parallel manipulators based on POE formula. Mech. Mach. Theory 2018, 121, 844–856. [Google Scholar] [CrossRef]
  26. Kong, L.; Chen, G.; Wang, H.; Huang, G.; Zhang, D. Kinematic calibration of a 3-PRRU parallel manipulator based on the complete, minimal and continuous error model. Robot. Comput.-Integr. Manuf. 2021, 71, 102158. [Google Scholar] [CrossRef]
  27. Sun, T.; Song, Y.; Li, Y.; Xu, L. Separation of comprehensive geometrical errors of a 3-DOF parallel manipulator based on Jacobian matrix and its sensitivity analysis with Monte-Carlo method. Chin. J. Mech. Eng. 2011, 24, 406–413. [Google Scholar] [CrossRef]
  28. Chen, Y.; Xie, F.; Liu, X.; Zhou, Y. Error modeling and sensitivity analysis of a parallel robot with SCARA (selective compliance assembly robot arm) motions. Chin. J. Mech. Eng. 2014, 27, 693–702. [Google Scholar] [CrossRef]
  29. Ma, J.G. Kinematics Analysis and Parameter Optimization of 2UPR-RPU Parallel Mechanism. Master’s Thesis, Zhejiang Sci-Tech University, Hangzhou, China, 2014. [Google Scholar] [CrossRef]
  30. Selig, J.M. Geometric Fundamentals of Robotics; Springer: New York, NY, USA, 2005. [Google Scholar] [CrossRef]
  31. Dai, J.S. Screw Algebra and Lie Groups and Lie Algebras; Higher Education Press: Beijing, China, 2014. [Google Scholar]
Figure 1. Schematic diagram of the 2UPR-RPU parallel mechanism.
Figure 1. Schematic diagram of the 2UPR-RPU parallel mechanism.
Machines 10 00990 g001
Figure 2. 2UPR-RPU parallel manipulator and its local coordinate systems.
Figure 2. 2UPR-RPU parallel manipulator and its local coordinate systems.
Machines 10 00990 g002
Figure 3. Scheme I: Evaluation of the limbs’ comprehensive deformations caused by geometric errors of the 2UPR-RPU over-constrained manipulator.
Figure 3. Scheme I: Evaluation of the limbs’ comprehensive deformations caused by geometric errors of the 2UPR-RPU over-constrained manipulator.
Machines 10 00990 g003
Figure 4. Reachable workspace of the 2UPR-RPU parallel manipulator.
Figure 4. Reachable workspace of the 2UPR-RPU parallel manipulator.
Machines 10 00990 g004
Figure 5. 692 target poses of the 2UPR-RPU parallel manipulator.
Figure 5. 692 target poses of the 2UPR-RPU parallel manipulator.
Machines 10 00990 g005
Figure 6. 1738 target poses of the 2UPR-RPU parallel manipulator.
Figure 6. 1738 target poses of the 2UPR-RPU parallel manipulator.
Machines 10 00990 g006
Figure 7. Simulation results using the standard deviations listed in Table 8. (a) Group 1; (b) Group 2; (c) Group 3.
Figure 7. Simulation results using the standard deviations listed in Table 8. (a) Group 1; (b) Group 2; (c) Group 3.
Machines 10 00990 g007
Figure 8. Global sensitivity of the average angular comprehensive deformation with respect to Geps.
Figure 8. Global sensitivity of the average angular comprehensive deformation with respect to Geps.
Machines 10 00990 g008
Figure 9. Global sensitivity of the average linear comprehensive deformation with respect to Geps.
Figure 9. Global sensitivity of the average linear comprehensive deformation with respect to Geps.
Machines 10 00990 g009
Table 1. Definitions of local coordinate systems for the first UPR limb.
Table 1. Definitions of local coordinate systems for the first UPR limb.
{Fi,j}The LocationFi,jxi,jzi,j
{F1,1}On the revolute shelfThe intersection of the right hole axis of the revolute shelf and the right end face of the revolute shelfParallel to the intersection of the front and rear symmetry plane of the right hole of the revolute shelf and the vertical plane of the right hole axisCoincide with the right hole axis of the revolute shelf
Point downPoint outwards
{F1,2}On the slider seatThe midpoint of the hole axis of the slider seatParallel to the intersection of the slider mounting plane and the vertical plane of the hole axis of the slider seatCoincide with the hole axis of the slider seat
Point to the moving platformPoint to the RPU limb
{F1,3}On the lead screwThe intersection of the lead screw axis and the plane passing through z1,2 and perpendicular to the slider mounting planeParallel to the intersection of the guide rail plane and the vertical plane of the lead screw axisCoincide with the lead screw axis
Point in the direction opposite to the RPU limbPoint to the moving platform
{F1,4}On the moving platformThe midpoint of the right hole axis of the moving platformParallel to the intersection of the vertical plane of the right hole axis of the moving platform and the plane constructed with v and wCoincide with the right hole axis of the moving platform
Point downPoint to the RPU limb
Table 2. Definitions of local coordinate systems for the second UPR limb.
Table 2. Definitions of local coordinate systems for the second UPR limb.
{Fi,j}The LocationFi,jxi,jzi,j
{F2,1}On the revolute shelfThe intersection of the left hole axis of the revolute shelf and the left end face of the revolute shelfParallel to the intersection of the front and rear symmetry plane of the left hole of the revolute shelf and the vertical plane of the left hole axisCoincide with the left hole axis of the revolute shelf
Point downPoint inwards
{F2,2}On the slider seatThe midpoint of the hole axis of the slider seatParallel to the intersection of the slider mounting plane and the vertical plane of the hole axis of the slider seatCoincide with the hole axis of the slider seat
Point to the moving platformPoint to the RPU limb
{F2,3}On the lead screwThe intersection of the lead screw axis and the plane passing through z2,2 and perpendicular to the slider mounting planeParallel to the intersection of the guide rail plane and the vertical plane of the lead screw axisCoincide with the lead screw axis
Point in the direction opposite to the RPU limbPoint to the moving platform
{F2,4}On the moving platformThe midpoint of the left hole axis of the moving platformParallel to the intersection of the vertical plane of the left hole axis of the moving platform and the plane constructed with v and wCoincide with the left hole axis of the moving platform
Point downPoint to the RPU limb
Table 3. Definitions of local coordinate systems for the RPU limb.
Table 3. Definitions of local coordinate systems for the RPU limb.
{Fi,j}The LocationFi,jxi,jzi,j
{F3,1}On the slider seatThe midpoint of the hole axis of the slider seatParallel to the intersection of the slider mounting plane and the vertical plane of the hole axis of the slider seatCoincide with the hole axis of the slider seat
Point to the moving platformPoint to the first UPR limb
{F3,2}On the lead screwThe intersection of the lead screw axis and the plane passing through z3,1 and perpendicular to the slider mounting planeParallel to the intersection of the guide rail plane and the vertical plane of the lead screw axisCoincide with the lead screw axis
Point to the second UPR limbPoint to the moving platform
{F3,3}On the U jointThe midpoint of the hole axis of the U jointParallel to the intersection of the vertical planes of the two hole axes of the U jointCoincide with the hole axis of the U joint
Point downPoint to the first UPR limb
{F3,4}On the moving platformThe intersection of the rear hole axis of the moving platform and the rear end face of the moving platformParallel to the intersection of the vertical plane of the rear hole axis of the moving platform and the plane constructed with v and wCoincide with the rear hole axis of the moving platform
Point downPoint to the RPU limb
Table 4. Structural parameters of the 2UPR-RPU parallel manipulator.
Table 4. Structural parameters of the 2UPR-RPU parallel manipulator.
SymbolsValuesUnits
lA0.06m
lB0.15m
c0.025m
d0.115m
Table 5. Specified geometric errors [13,26] for geometric error identification.
Table 5. Specified geometric errors [13,26] for geometric error identification.
SymbolsGroup 1Group 2Group 3Units
δi,j0.0050.0015 × 10−5m
εi,j0.005π/180π/7200rad
Table 6. Initially identified internal-force-and-deformation-related geometric errors.
Table 6. Initially identified internal-force-and-deformation-related geometric errors.
ij δ i , j x δ i , j y δ i , j z ε i , j x ε i , j y ε i , j z
1, 20
1, 21
1, 22
1, 23
1, 24
30
31
32
33
34
Table 7. Identified internal-force-and-deformation-related geometric errors.
Table 7. Identified internal-force-and-deformation-related geometric errors.
ij δ i , j x δ i , j y δ i , j z ε i , j x ε i , j y ε i , j z
1, 20
1, 21
1, 22
1, 23
1, 24
30
31
32
33
34
Table 8. Standard deviations of the geometric errors for the numerical simulations.
Table 8. Standard deviations of the geometric errors for the numerical simulations.
SymbolsGroup 1Group 2Group 3Units
The standard deviations of δi,j1.6667 × 10−33.3333 × 10−51.6667 × 10−5m
The standard deviations of εi,j1.6667 × 10−3π/540π/21,600rad
Table 9. Sensitivity analysis results of the average angular comprehensive deformation.
Table 9. Sensitivity analysis results of the average angular comprehensive deformation.
Group NumberGe5 ,   G e 7 ,   G e 8 ,   G e 11 ,   G e 14 ,   δ 2 , 1 y , ε 2 , 1 z ,   δ 2 , 2 z ,   δ 2 , 3 x ,   δ 2 , 4 y [mm or °]Other Geometric Errors [mm or °] The   Maximum   Value   of   Δ ω a j [°] The   Average   Value   of   Δ ω a j [°]
Group 10.10.10.14300.0961
Group 20.010.10.14300.0961
Group 30.0010.10.14300.0961
Table 10. Specified geometric errors for verification.
Table 10. Specified geometric errors for verification.
ij δ i , j x [mm] δ i , j y [mm] δ i , j z [mm] ε i , j x [°] ε i , j y [°] ε i , j z [°]
1, 200.10.01770.03810.00540.0033
1, 210.10.01740.10.00531
1, 220.10.01730.00560.00920.1
1, 230.01740.10.00370.0056
1, 240.10.01740.10.00280.0022
300.10.10.10.00790.0057
310.10.10.00830.01350.1
320.10.10.01370.0084
330.10.10.10.0061
340.10.10.10.00370.0024
Table 11. Sensitivity analysis results of the average linear comprehensive deformation.
Table 11. Sensitivity analysis results of the average linear comprehensive deformation.
Group Number The   Maximum   Value   of   Δ ω a j [°] The   Average   Value   of   Δ ω a j [°] The   Maximum   Value   of   Δ v a j [mm] The   Average   Value   of   Δ v a j [mm]
Group 10.01650.01180.06960.0390
Group 20.10610.07140.83740.3581
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Du, X.; Wang, B.; Zheng, J. Geometric Error Analysis of a 2UPR-RPU Over-Constrained Parallel Manipulator. Machines 2022, 10, 990. https://doi.org/10.3390/machines10110990

AMA Style

Du X, Wang B, Zheng J. Geometric Error Analysis of a 2UPR-RPU Over-Constrained Parallel Manipulator. Machines. 2022; 10(11):990. https://doi.org/10.3390/machines10110990

Chicago/Turabian Style

Du, Xu, Bin Wang, and Junqiang Zheng. 2022. "Geometric Error Analysis of a 2UPR-RPU Over-Constrained Parallel Manipulator" Machines 10, no. 11: 990. https://doi.org/10.3390/machines10110990

APA Style

Du, X., Wang, B., & Zheng, J. (2022). Geometric Error Analysis of a 2UPR-RPU Over-Constrained Parallel Manipulator. Machines, 10(11), 990. https://doi.org/10.3390/machines10110990

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop