Abstract

This article proposes a new approach to evaluate the stiffness of a three degrees-of-freedom decoupled Cartesian parallel manipulator which uses revolute joints. Prior to the discussion on the new approach, the stiffness of an individual limb is discussed. This discussion focuses on the influence of kinematic parameters on the limb's stiffness. Taking advantage of the decoupled behavior of the Cartesian parallel manipulator, the stiffness evaluation is initially performed in multiple analyses at the individual limb level. Subsequently, the end-effector deformation is obtained by summing the individual limbs' deformations. This approach does not require the assembly of the complete manipulator. Accordingly, it offers a simpler and more intuitive way to analyze the stiffness of the manipulator. The approach was demonstrated in a 3PRRR decoupled parallel manipulator. The result of the proposed approach was compared with that of a standard approach that considers a completely assembled manipulator. The comparison shows an acceptable agreement, suggesting that the proposed approach can quickly estimate the stiffness of the entire manipulator. The stiffness evaluations in the article employed several methods, including analytical approximation, matrix structural analysis, and finite element analysis.

1 Introduction

The high stiffness-to-mass ratio is often seen as a major advantage of parallel manipulators. In the design stage of a parallel manipulator, one typically evaluates the manipulator's stiffness and subsequently optimizes it along with other important measures, such as mass and workspace, in a multi-objective optimization. The stiffness analysis evaluates the deformation of the manipulator with certain geometry, boundary conditions, and material, due to an applied wrench. The manipulator's stiffness can be expressed as a stiffness matrix or a resulting deformation due to a specified applied wrench. In the former case, the elements of the stiffness matrix indicate the amount of wrench elements required to cause a unit deformation in the constrained degrees-of-freedom (DOF). In the latter case, the amount of deformation due to the specified wrench indicates the stiffness. There are two common purposes for stiffness evaluation of a manipulator. The first purpose, which is the purpose of this work, is to be used for design optimization. This usually aims to optimize the dimension of the manipulator while assuming the use of a certain material. With this goal in mind, an overall optimization of the manipulator geometry typically starts with the optimization of the links' lengths. Once the optimum links' lengths are obtained, other geometric parameters, such as the cross-sectional size and the detailed links' topology, can be subsequently optimized. The second purpose of a stiffness evaluation is to quantify the absolute deformation of the manipulator given a certain wrench. When this is the goal, it is important to include all the stiffness components or at least all the significant stiffness components. Furthermore, the detailed geometry and material of the manipulator's parts should be considered when the absolute value of the manipulator's deformation is the concern. In the following section, several stiffness evaluation methods will be reviewed to discuss their advantages and disadvantages, and the methods will be used selectively in this article.

Based on the coupling condition, parallel manipulators can be classified into fully coupled, partially decoupled (group decoupled), and fully decoupled manipulators [1]. In fully decoupled manipulators, each end-effector (EE) coordinate depends on only one active joint coordinate independently [2]. In other words, the task space coordinates and the joint space coordinates have a one-to-one correspondence relationship. While the stiffness of a decoupled parallel manipulator can be analyzed in a similar manner to that of general parallel manipulators, the authors cannot find a work discussing a method to evaluate specifically the stiffness of a decoupled parallel manipulator except Ref. [3]. Nevertheless, the stiffness evaluation in that work has the following limitations:

  • The work only performed the stiffness evaluation in one Cartesian direction, and it did not evaluate the manipulator stiffness subject to a combined force, i.e., a force consisting of non-zero components working in all the Cartesian directions.

  • The work did not investigate the stiffness of the manipulator due to an applied moment.

  • The work did not explicitly present a systematic method to evaluate the stiffness of the entire manipulator subject to general force and moment applied to the end-effector.

  • The work did not study how the variation of the manipulator kinematic parameters affects the manipulator stiffness.

In this work, we overcome the aforementioned limitations by proposing a systematic approach and performing parametric studies to evaluate the stiffness of a decoupled parallel manipulator. This work was driven by our need to evaluate the stiffness of a 3DOF decoupled Cartesian parallel manipulator [4,5] composed of three serial RRR links in a quick and intuitive manner. We will start by evaluating the stiffness of an individual limb with two links connected by R joints with parallel joints' axes. The study on the individual two-link limb focuses on the behavior of the limb before it is assembled into the entire manipulator. Subsequently, we will evaluate the stiffness of the entire 3DOF decoupled Cartesian parallel manipulator using a systematic approach we propose. The idea of the approach is to break down the stiffness evaluation into multiple separate analyses corresponding to elementary Cartesian wrench components resolved from a general wrench applied to the manipulator's end-effector. The separate analyses are performed based on the decoupled behavior of the Cartesian parallel manipulator. The proposed method is simpler and more intuitive than the standard stiffness evaluation method, which requires an assembly procedure. However, it comprises multiple analyses and only applies to decoupled Cartesian parallel manipulators. Using the approach we propose, as will be shown in this article, the stiffness of the entire manipulator can be evaluated at the individual limbs in a simple manner and subsequently completed by superposing the results to obtain the stiffness of the entire manipulator. Using the proposed approach, a simpler and more intuitive stiffness evaluation can be performed in the manipulator design stage. This study complements the previous work [3] by the following: (1) proposing a systematic method to estimate the stiffness of the manipulator subject to combined force and moment, (2) comparing the result of the proposed approach with the standard approach that is based on an assembled manipulator model. Furthermore, this study also aims to investigate the influence of the manipulator's kinematic parameters on its stiffness. Since the aim of the stiffness evaluation is topology and dimensional optimization, the compliance of the manipulator in this work is assumed to be given by the compliance of its links. The manipulator's joints are assumed to be rigid.

The remainder of this article is organized as follows. Section 2 briefly reviews the methods used to evaluate the stiffness of a manipulator. Section 3 presents the analysis setup used in this work, including the model boundary conditions, material properties, and geometrical dimensions of the analyzed mechanism. Section 4 evaluates the out-of-plane stiffness of the individual two-link limb. The term “plane” here is defined as the plane in which the links of the limb perform their motion. After the discussion of the individual limb stiffness, Sec. 5 presents our simplified approach to evaluate the stiffness of an entire decoupled Cartesian parallel manipulator without a need to perform assembly. Finally, Sec. 6 concludes the article. Throughout the article, the revolute and prismatic joints are indicated by R and P, respectively. An underlined joint indicates that the joint is actuated.

2 Review of Manipulator Stiffness Evaluation

2.1 Stiffness Modeling.

The stiffness modeling of a manipulator can be broadly classified into lumped and distributed modeling methods. The most commonly used lumped model is the Jacobian-based virtual joint model (VJM). Using this method, it is quite common to select one or more stiffness components, such as the actuator stiffness [611], the axial stiffness of the links [812], the bending stiffness of the links [7,9,10], and the torsional stiffness [11] of the links. To incorporate the stiffness in all translational and rotational directions, Pashkevich et al. [13] proposed a VJM method that uses 6DOF springs to represent the stiffness of the manipulator's links. A VJM method based on the overall Jacobian [14] that consists of the actuation and constraint Jacobians defines stiffness as the relationship between the applied wrench and the displacements of the actuations and the constraints [11,15]. While the determination of the actuators' stiffness is quite straightforward, the determination of the constraints' stiffness is not always straightforward. Hoevenaars et al. [16] proposed a VJM method based on the generalized Jacobian [17] that is broader than the overall Jacobian. Compared with the VJM method proposed by Pashkevich et al. [13], this method results in a smaller dimension, i.e., the size of matrices. However, this advantage is most likely not quite significant, particularly with the current computing hardware technology, as any VJM method does not have a high dimension.

The most commonly used distributed models are the matrix structural analysis (MSA) and finite element analysis (FEA). The FEA is convenient to accommodate complex geometries with high accuracy by using suitable element types and meshing strategies. An experimental investigation by Pinto et al. [3] shows that the FEA model is only different by 0–2% from the experiment. However, the stiffness evaluation using the FEA should be performed for every manipulator's posture since an FEA user typically does not have the parametrical expressions of the finite element model of the manipulator. The MSA, sometimes called the direct stiffness method, can be seen as a special case of the FEA as it is basically FEA performed through matrix operations that are conveniently used for one-dimensional finite elements such as bars and beams. Since the MSA is only applied to simple one-dimensional finite elements, it is quite feasible to code the parametrical expressions of the model. Accordingly, the model can be used across the postures of the modeled manipulator without a need to re-setup for every posture. Klimchik et al. [18] proposed a new formulation of the MSA in which they introduced boundary condition equations (matrices) that eliminate the necessity to remove the rows and columns corresponding to the boundary conditions from the stiffness matrix as in the conventional MSA [19].

Pashkevich et al. [13] and Klimchik et al. [20] proposed a hybrid approach that combines VJM with the FEA through stiffness matrix identification using the FEA. This is achieved by performing virtual experiments in an fea software. Another hybrid approach is performing a model order reduction to an FEA model. Taghvaeipour et al. [21] obtained the stiffness matrices of complex anisotropic links through substructuring to obtain reduced-order stiffness matrices. Furthermore, Cammarata [22] employed several model order reduction methods, namely Guyan, Craig–Bampton, and Improved Reduced System (IRS) methods, to reduce the order of a finite element model of a multibody elastodynamic system.

2.2 Stiffness Indices.

Stiffness indices are typically used as a handy scalar indicator of the manipulator stiffness. The VJM provides a 6 × 6 stiffness Cartesian matrix (CSM). The stiffness magnitude is commonly measured by using the trace and the weighted trace of the CSM, the minimum, maximum, and average eigenvalues of the CSM, the mean value of the CSM, and sometimes the determinant of the CSM. The stiffness uniformity is typically measured by using the condition number of the CSM. Despite the convenience of extracting a stiffness index from the CSM, the stiffness index suffers from the non-uniformity problem when a manipulator mixes translational and rotational degrees-of-freedom. There are typically two strategies to deal with such a problem. The first strategy is separating the evaluation between the translational and rotational stiffness components. The second strategy is homogenizing the stiffness matrix. Some global indices such as the global compliance index [7] and the global stiffness index (GSI) [23] are commonly used to evaluate a manipulator's stiffness across the workspace.

The distributed stiffness models such as MSA and FEA typically do not provide a 6 × 6 Cartesian stiffness matrix. Using these models, the computed deformation of the end-effector is typically used as the base of a stiffness index. A global stiffness index is typically obtained by averaging the deformations across the workspace. To have an alternative stiffness index based on a distributed stiffness model, Raoofian et al. [24] proposed an enhanced MSA model, verified it with an FEA model, and proposed a method to obtain a 6 × 6 CSM from the enhanced MSA model. Subsequently, they proposed a stiffness index based on the obtained 6 × 6 CSM.

3 Analysis Setup

Figure 1 shows the two-link limb to be evaluated in this work. This two-link limb also serves as the main mechanism of the decoupled Cartesian parallel manipulator discussed at the end of this article. Point O is the fixed position of the proximal R joint. The proximal link makes an angle α1 with respect to the X-axis. This angle is called the limb angle across this article. In practice, the virtual segment connecting point O and the end-effector can be oriented at any angle. For convenience, this virtual segment will be illustrated as a horizontal line coincident with the X-axis, implying that the end-effector keeps its rigid body position along the X-axis despite the change of the angle θ and the links' lengths L1 and L2. Assuming that each R joint is passive, the moment about the joint axis is not constrained; hence, it cannot be modeled as a fixed joint. However, when the two-link limb with passive R joints is subject to only an out-of-plane load, the passive R joints can be modeled as fixed joints. In this work, the two-link limb is only subject to an out-of-plane load. Accordingly, the passive R joints are modeled as fixed joints.

Fig. 1
(a) Two-link limb with passive RR joints and (b) its simplification assuming rigid joints
Fig. 1
(a) Two-link limb with passive RR joints and (b) its simplification assuming rigid joints
Close modal

Table 1 shows the material properties and geometrical dimensions of the links used in this work. The links are made from an aluminum alloy with the mentioned properties. Unless stated otherwise, the links of the mechanism have identical lengths, i.e., L1=L2.

Table 1

Materials properties and geometrical dimensions of the links

ParameterUnitValue
Young’s modulus (E)GPa71.1
Shear modulus (G)GPa26.7
Link length L1=L2mm250
Link square cross sectionmm30 × 30
ParameterUnitValue
Young’s modulus (E)GPa71.1
Shear modulus (G)GPa26.7
Link length L1=L2mm250
Link square cross sectionmm30 × 30
Figure 2 shows two-link limbs with passive R joints at a general configuration. The two links create an angle θ between them. The angle α1 is measured between link 1 and a virtual line connecting O and B. The virtual line is not necessarily oriented horizontally. However, it is oriented horizontally here for convenience. The axis 1 is perpendicular to link 1, OA, whereas the axis 2 is perpendicular to link 2, AB. Figures 2(a) and 2(b) depict two-link limbs with L1L2 and with L1=L2, respectively. The distance between O and B is defined as the effective length of the limb, L. Applying the law of sines yields the following:
(1)
Fig. 2
Two-link limb (a) with L1≠L2 at an arbitrary angle θ and (b) with L1=L2 at an arbitrary angle θ
Fig. 2
Two-link limb (a) with L1≠L2 at an arbitrary angle θ and (b) with L1=L2 at an arbitrary angle θ
Close modal
Assuming constant link lengths and constrained motion of the end-effector B only along a line collinear with the virtual line OB, the angle θ is dependent on the angle α1
(2)

The EE deformation in this work is evaluated using the analytical, MSA, and FEA methods. We proposed some analytical approximation in this work to evaluate the end-effector deformation. We coded the MSA solver using the Euler–Bernoulli beams with six degrees-of-freedom at each node to model the links. Each link has only nodes at its ends. On the other hand, we obtained the FEA solution using ansys mechanical. The FEA model discretizes each link into 20 beam elements. Figure 3 shows the FEA models in ansys mechanical. Point A is fixed. Point B is the end-effector of the two-link limb. In Fig. 3(a), the out-of-plane force, working in a direction perpendicular to the paper, is applied at point B. In Fig. 3(b), the out-of-plane moment, with its rotation axis lying on the plane and perpendicular to the link, is applied at point B. The parametric evaluation in the ansys fea was performed for several discrete values of the parameters within the range of interest. To evaluate the stiffness across the workspace, we plotted stiffness maps based on the end-effector deformation computed using the MSA method.

Fig. 3
FEA model of two-link limb subject to (a) an out-of-plane force and (b) an out-of-plane moment at its tip
Fig. 3
FEA model of two-link limb subject to (a) an out-of-plane force and (b) an out-of-plane moment at its tip
Close modal

4 Out-of-Plane Stiffness of Two-Link Limb With Passive R Joints

4.1 Analytical Solution for General Configuration.

Given an out-of-plane force F applied to the end-effector, the out-of-plane end-effector deformation at an arbitrary angle θ can be approximated as [3]
(3)
where I and J are the area moment of inertia and the polar moment of inertia of the link cross section, respectively.
If an out-of-plane moment M is applied to the end-effector, we propose that the out-of-plane end-effector deformation due to the applied moment can be approximated as
(4)

4.2 Analytical Solution in Special Cases.

Now consider the two-link limb in three special cases, namely fully extended, fully folded, and right-angled configurations, as depicted in Fig. 4. Although the three configurations are typically feasible to achieve in a single limb, some of them, particularly the fully extended and fully folded configurations, can be impractical to achieve in multiple limbs composing a parallel manipulator, due to the kinematic constraints of the manipulator. The out-of-plane deformation at the end-effector due to an out-of-plane force Fz and an out-of-plane moment My in the three special cases, as shown in Table 2, can be derived using the textbook mechanics of materials. The analytical solutions shown in Table 2 are identical to those derived from Eqs. (3) and (4) by setting the angles to values corresponding to the special cases. The table also compares the analytical solution with the MSA and FEA solutions. Since the analytical and MSA solutions are identical, the percentage error shown in the table indicates the difference between either of them and the FEA solution.

Fig. 4
A two-link limb in its special cases: (a) fully extended, (b) right-angled, and (c) fully folded configurations
Fig. 4
A two-link limb in its special cases: (a) fully extended, (b) right-angled, and (c) fully folded configurations
Close modal
Table 2

Out-of-plane EE deformation due to out-of-plane load in special cases

Configurationθ
(deg)
α2
(deg)
LoadAnalytic solutionMSA solution
(mm)
ansys fea solution
(mm)
Percent error (%)
FormulaValue
(mm)
Fully extended1800FzδEE,z=Fz(L1+L2)33EI0.43470.43470.43590.3
MyδEE,z=My(L1+L2)22EI0.65210.65210.65210
Right-angled9090FzδEE,z=FzL133EI+FzL233EI+FzL1L22GJ0.32550.32550.363310.4
MyδEE,z=MyL222EI+MyL1L2GJ0.59670.59670.669710.8
Fully folded0180FzδEE,z=FzL133EI+FzL233EI+FzL1L22EIFzL12L2EI0.10870.10870.11091.9
MyδEE,z=MyL122EI+MyL222EI0.32600.32600.32770.5
Configurationθ
(deg)
α2
(deg)
LoadAnalytic solutionMSA solution
(mm)
ansys fea solution
(mm)
Percent error (%)
FormulaValue
(mm)
Fully extended1800FzδEE,z=Fz(L1+L2)33EI0.43470.43470.43590.3
MyδEE,z=My(L1+L2)22EI0.65210.65210.65210
Right-angled9090FzδEE,z=FzL133EI+FzL233EI+FzL1L22GJ0.32550.32550.363310.4
MyδEE,z=MyL222EI+MyL1L2GJ0.59670.59670.669710.8
Fully folded0180FzδEE,z=FzL133EI+FzL233EI+FzL1L22EIFzL12L2EI0.10870.10870.11091.9
MyδEE,z=MyL122EI+MyL222EI0.32600.32600.32770.5

4.3 Parametric Studies.

Figure 5 shows one-link and two-link limbs with an out-of-plane force, perpendicular to the paper, applied to the tip of the limbs. For simplicity, let us assume a one-link limb. In this simple case, the bending problem is similar to that in the bending of a straight cantilever beam. The smaller the limb length L, the stiffer the limb is, and vice versa. In the case of the two-link limb, the appearance of the limb is not only dictated by the effective length of the limb, L, but also W1 and W2. Assuming an identical length of the two links, W1 = W2 = W, L and W are dependent on each other if the link length is already given. The smaller L is, the larger W is, and vice versa.

Fig. 5
A limb seen from the bending force direction
Fig. 5
A limb seen from the bending force direction
Close modal

First of all, let us assume that the one-link and two-link limbs have identical link cross section and material. If both the limbs have identical L, we can hypothesize that they do not have an identical stiffness. Furthermore, if two-link limbs have identical L but different W, we also hypothesize that they do not have identical stiffness. In a single two-link limb, we are interested in evaluating its stiffness at its varied posture. In this case, a specific posture is defined by its limb angle, effective length L, and folding depth W, which depend on each other. To study the effects of the parameter values, parametric studies of a two-link limb were performed using analytical analysis based on the analytical, MSA, and ansys fea solutions.

4.3.1 The Effect of the Limb Angle.

Here, the out-of-plane stiffness of the limb is evaluated, and the angle of the limb is varied. Figure 6 shows the out-of-plane end-effector deformation of the limb due to a separate out-of-plane force of 50 N and an out-of-plane moment of 25 Nm. The figure shows that the out-of-plane stiffness is affected by the limb angle. Furthermore, the data show that the larger the limb angle α1, the stiffer the limb in the out-of-plane direction. This is an intuitive behavior since a larger limb angle α1 corresponds to a shorter effective length of the limb. The data also show that the analytical and MSA solutions are identical. Both solutions do not have a significant difference from the FEA solution at the fully extended and fully folded configurations. The closer the limb to the right-angled configuration, the larger the difference between the analytical/MSA and FEA solutions. As shown in Table 2, the largest difference is around 10%, occurring at the right-angled configuration.

Fig. 6
Out-of-plane EE deformation due to (a) out-of-plane force and (b) out-of-plane moment over limb angle, obtained using the analytical solution, MSA, and ansys fea. The EE deformation refers to the left vertical axes, whereas the effective length refers to the right vertical axes.
Fig. 6
Out-of-plane EE deformation due to (a) out-of-plane force and (b) out-of-plane moment over limb angle, obtained using the analytical solution, MSA, and ansys fea. The EE deformation refers to the left vertical axes, whereas the effective length refers to the right vertical axes.
Close modal

4.3.2 The Effects of the Links' Lengths.

First, the links' lengths vary while the total length of the two-link limb is maintained constant. The effective length is not maintained to be constant. This evaluation was performed at three representative limb angles, namely 30 deg, 50 deg, and 80 deg. Figure 7 shows the out-of-plane end-effector deformation in this evaluation. The figure shows that the end-effector deformation follows the effective length. The smaller the effective length, the smaller the end-effector deformation, and vice versa.

Fig. 7
Out-of-plane EE deformation of a two-link limb subject to an out-of-plane force applied to the end-effector, computed in MSA. The EE deformation refers to the left vertical axis, whereas the total length and effective length refer to the right vertical axis.
Fig. 7
Out-of-plane EE deformation of a two-link limb subject to an out-of-plane force applied to the end-effector, computed in MSA. The EE deformation refers to the left vertical axis, whereas the total length and effective length refer to the right vertical axis.
Close modal

To exclude the influence of the effective length, another evaluation was performed by varying the link lengths while keeping both the total length and effective length of the limb constant. Accordingly, each pair of L1 and L2 corresponds to a unique limb angle. This evaluation was performed in the MSA and ansys fea. In this evaluation, an out-of-plane force of 100 N and an out-of-plane moment of 25 Nm are separately applied to the end-effector. Figure 8(a) shows the ansys fea model, whereas Fig. 8(b) shows the limb postures at all the varied link lengths. The out-of-plane end-effector deformation of the limb is shown in Fig. 9. The figure shows that a larger link length ratio results in a smaller out-of-plane end-effector deformation.

Fig. 8
(a) The FEA model of the two-link limb subject to out-of-plane force and moment at its tip and (b) the plot of the limb postures upon varying the link lengths while keeping the total length and effective length constant
Fig. 8
(a) The FEA model of the two-link limb subject to out-of-plane force and moment at its tip and (b) the plot of the limb postures upon varying the link lengths while keeping the total length and effective length constant
Close modal
Fig. 9
Out-of-plane EE deformation of a two-link limb subject to an out-of-plane force and moment applied to the end-effector, computed in the MSA and ansys fea
Fig. 9
Out-of-plane EE deformation of a two-link limb subject to an out-of-plane force and moment applied to the end-effector, computed in the MSA and ansys fea
Close modal

5 Simplified Stiffness Evaluation of the 3PRRR Parallel Manipulator

After discussing the stiffness of a two-link limb with passive R joints, let us discuss the stiffness of a 3PRRR parallel manipulator, which has three RRR limbs with passive R joints and loads applied perpendicular to the motion plane of the limbs. An actuated P joint is added to the proximal end of each limb, transforming the limb from RRR to PRRR topology. Figure 10(a) shows the 3PRRR parallel manipulator. The underlined letter indicates the actuated joint. It is a 3DOF decoupled Cartesian parallel manipulator. The first link in each limb is typically very short, so its stiffness can be neglected. Accordingly, the manipulator topology can be simplified, as shown in Fig. 10(b). Since all the R joints in each limb are parallel, the limbs' motion planes are always parallel with the Cartesian planes, namely YZ, XZ, and XY planes.

Fig. 10
(a) Topology of 3PRRR parallel manipulator and (b) its simplified topology
Fig. 10
(a) Topology of 3PRRR parallel manipulator and (b) its simplified topology
Close modal

5.1 Stiffness Against Forces.

For convenience, let us see the parallel manipulator at one of its postures in a two-dimensional view, as depicted in Fig. 11. For simplicity, let us assume that the manipulator's end-effector is only subject to an external force. In other words, an external moment is not applied to the end-effector. With the P joints locked at their positions, namely x1, x2, and x3, let us apply an external force to the manipulator's end-effector. This force can be resolved into three components, namely Fx, Fy, and Fz, as depicted in Fig. 10.

Fig. 11
3PRRR parallel manipulator in two-dimensional view, with external forces Fx, Fy, and Fz (perpendicular to the paper) applied to its end-effector
Fig. 11
3PRRR parallel manipulator in two-dimensional view, with external forces Fx, Fy, and Fz (perpendicular to the paper) applied to its end-effector
Close modal

When an external force in the X-direction, namely Fx, is applied, the locked limb 1 prevents the end-effector from a rigid body translation in the X-direction. Accordingly, limb 1 is suffering from a transversal bending in the X-direction due to the force Fx. This leads to a deformation in the X-direction, namely dx1 (with the superscript 1 indicating the limb and the subscript x indicating the deformation direction). Meanwhile, limbs 2 and 3 cannot resist the deformation of limb 1 in the X-direction since all the R joints in both limbs are passive. This is a physical consequence of the decoupled characteristics of the manipulator. The out-of-plane deformation of a limb is only resisted by the structure of that particular limb. A similar behavior applies to limbs 2 and 3 when forces Fy and Fz are applied, respectively. This behavior is summarized in Table 3.

Table 3

Bending and its resulting deformations in the limbs due to an external force

External force componentLimb subject to the transversal bendingDeformation due to the external forceLimbs not resisting deformation in the external force direction
FxLimb 1dx1Limbs 2 and 3
FyLimb 2dy2Limbs 1 and 3
FzLimb 3dz3Limbs 1 and 2
External force componentLimb subject to the transversal bendingDeformation due to the external forceLimbs not resisting deformation in the external force direction
FxLimb 1dx1Limbs 2 and 3
FyLimb 2dy2Limbs 1 and 3
FzLimb 3dz3Limbs 1 and 2

Consequently, the end-effector's deformations dx1, dy2, and dz3 are only resisted by the individual limbs 1, 2, and 3, respectively. More specifically, the resistance comes from the structural constraint of the individual limb that behaves as a cantilever subject to bending. When an external force applied to the end-effector has a diagonal direction, i.e., not purely in X, Y, or Z direction, then two or more limbs are subject to bending. By resolving the external force into its Cartesian components, a separate analysis can be performed for each loaded limb.

On the other hand, rigid body rotations about the X, Y, and Z axes are constrained collectively by the three limbs. These rotational constraints are also handled by the structure of the limbs. Table 4 summarizes the constrained rotations in each limb. Since the constraints of the end-effector are the sum of the constraints of the limbs, it is easy to deduce from Table 4 that the end-effector is constrained by rotation about the X, Y, and Z axes. Accordingly, the manipulator has only 3T degrees-of-freedom.

Table 4

Constrained rotations in the limbs

LimbConstrained rotation
Limb 1About Y and Z
Limb 2About X and Z
Limb 3About X and Y
LimbConstrained rotation
Limb 1About Y and Z
Limb 2About X and Z
Limb 3About X and Y

It is worth mentioning that the rotational constraints exist at all times, regardless of whether the actuators are locked or not. In contrast, the translational constraints theoretically apply only when the corresponding actuators are locked. Nevertheless, in practice, if the translational actuation is implemented using lead screws, the translational constraints at all times can be achieved due to the translational resistance of the lead screw mechanism, regardless of whether the actuators are locked or not.

5.2 Stiffness Against Moments.

Consider the 3PRRR manipulator with a moment applied to its end-effector. The moment can be resolved into three components about the X-, Y-, and Z-axis, namely Mx, My, and Mz, respectively, as depicted in Fig. 10. Let us consider the application of a single moment component at a time. When Mx is applied to the end-effector, limbs 2 and 3 resist the rotation, but limb 1 does not resist the rotation since all the R joints in the limb are passive. Similar behavior applies when My and Mz are separately applied to the end-effector. Table 5 summarizes this behavior. Accordingly, a limb only resists two moment components as shown in Fig. 12 as it cannot resist the moment component working about an axis parallel with the R joints' axes. Table 6 summarizes the resisted moments in each limb. Comparing the behavior of the manipulator against forces with that against moments, it can be concluded that a force component is only resisted by a single limb that moves in a plane perpendicular to the force, whereas a moment component is resisted by two limbs whose R joints' axes are not parallel with the moment axis.

Fig. 12
The resisted moments in limb 1 of the 3PRRR parallel manipulator subject to an external moment
Fig. 12
The resisted moments in limb 1 of the 3PRRR parallel manipulator subject to an external moment
Close modal
Table 5

Moments applied to the end-effector and the limbs' resistance to the moments

Moment componentLimbs resisting the rotationLimbs not resisting the rotation
MxLimbs 2 and 3Limb 1
MyLimbs 1 and 3Limb 2
MzLimbs 1 and 2Limb 3
Moment componentLimbs resisting the rotationLimbs not resisting the rotation
MxLimbs 2 and 3Limb 1
MyLimbs 1 and 3Limb 2
MzLimbs 1 and 2Limb 3
Table 6

Moment components resisted in each limb

LimbMoment components resistedMoment component is not resisted
Limb 1My and MzMx
Limb 2Mx and MzMy
Limb 3Mx and MyMz
LimbMoment components resistedMoment component is not resisted
Limb 1My and MzMx
Limb 2Mx and MzMy
Limb 3Mx and MyMz

5.3 Stiffness of the Entire Manipulator.

Based on the behavior of the manipulator's limbs under a force and a moment applied to the manipulator's end-effector as discussed in Secs. 5.1 and 5.2, the stiffness of the entire manipulator can be evaluated through multiple separate analyses at the individual limb level. The evaluation is performed in the following steps:

Step 1: Given a general force and moment applied to the manipulator's end-effector, decompose the applied force and moment into their Cartesian components, namely Fx, Fy, Fz, Mx, My, and Mz.

Step 2: Solve the kinematics of the manipulator. This is required to obtain the posture of each limb.

Step 3: Analyze the deformation of the resisting limb(s) using any feasible method. Please note that the number of separate analyses depends on the number of force and moment components applied to the end-effector. For example, only one analysis is required if only one force component or one moment component exists. At maximum, nine analyses are required in the presence of three force components (Fx, Fy, and Fz) and three moment components (Mx, My, and Mz).

As shown in Table 3, the limb tip out-of-plane deformation evaluation due to forces only considers one force component. Similarly, only one moment component should be used in the evaluation of each limb, although Table 6 shows that two moment components are resisted in each limb. This is to avoid using a moment component more than once, because the deformation of the manipulator's end-effector is obtained by summing the end-effector's deformation obtained from all the limbs. Observing limb 1 in Fig. 12, we notice that Mz bends the limb whereas My rotates the whole limb about the Y-axis. It was observed from the simulation that the end-effector deformation due to the bending moment (Mz) is more dominant. Hence, it is used in the evaluation of limb 1. Similarly, Mx and My are used in the evaluation of limbs 2 and 3, respectively.

Step 4: Sum the deformations from all the analyses to obtain the deformation of the end-effector due to all the force and moment components working at the end-effector.

The four aforementioned steps only require the assembly of the links in each limb. The assembly of all the limbs into a complete manipulator is not required. Although this approach requires multiple individual analyses, each analysis is simpler, more intuitive, and easier to verify since it is performed in the individual limb that involves a small number of links. This approach also provides better knowledge of the behavior of the manipulator, such as its end-effector deformation, under each load component.

5.4 Example.

Consider a 3PRRR parallel manipulator with proximal and distal links in all the limbs that have identical lengths of 250 mm. All the links have a square cross section of 30 × 30 mm and are made of aluminum with the Young modulus E of 71.1 GPa and the shear modulus G of 26.7 GPa. The end-effector offset in the X, Y, and Z directions is equally 80 mm. With the manipulator at a pose of x=y=z=350mm, forces Fx=Fy=Fz=50N and moments Mx=My=Mz=25Nm are separately and altogether applied to its end-effector.

5.4.1 Deformation Against Forces.

Consider force components Fx=Fy=Fz=50N are applied to the end-effector. Table 7 shows the end-effector deformation due to separate force components, evaluated at the limb level. The last line of the table shows the end-effector deformation if all the force components are applied to the end-effector. It is obtained by simply summing the end-effector deformations due to the separate force components, without a need to assemble the complete manipulator.

Table 7

EE deformation in the 3PRRR parallel manipulator due to forces, computed with the proposed approach using the analytical/MSA solution

LimbEnd-effector wrench componentdx(mm)dy(mm)dz(mm)
Limb 1Fx=50N−0.343200
Limb 2Fy=50N0−0.34320
Limb 3Fz=50N00−0.3432
Combined (Fx=Fy=Fz=50N)−0.3432−0.3432−0.3432
LimbEnd-effector wrench componentdx(mm)dy(mm)dz(mm)
Limb 1Fx=50N−0.343200
Limb 2Fy=50N0−0.34320
Limb 3Fz=50N00−0.3432
Combined (Fx=Fy=Fz=50N)−0.3432−0.3432−0.3432

Table 8 shows the deformation of the end-effector computed using various methods, including the simplified approach proposed in this article. The table shows that the deformation computed using the proposed approach is close to that computed using the assembled manipulator model.

Table 8

EE deformation of the 3PRRR parallel manipulator due to forces

Deformation componentDeformation using the proposed approach
(mm)
Deformation using the assembled manipulator model
MSA
(mm)
ansys fea
(mm)
dx−0.3432−0.3388−0.3275
dy−0.3432−0.3388−0.3275
dz−0.3432−0.3388−0.3275
Deformation componentDeformation using the proposed approach
(mm)
Deformation using the assembled manipulator model
MSA
(mm)
ansys fea
(mm)
dx−0.3432−0.3388−0.3275
dy−0.3432−0.3388−0.3275
dz−0.3432−0.3388−0.3275

Figure 13 shows the deformation of the 3PRRR manipulator in the MSA model. Figure 14 shows the meshed manipulator in the ansys fea model.

Fig. 13
The undeformed (dashed line) and deformed (continuous line) 3PRRR manipulator solved using the MSA, with deformation shown at the amplifying factor of 50
Fig. 13
The undeformed (dashed line) and deformed (continuous line) 3PRRR manipulator solved using the MSA, with deformation shown at the amplifying factor of 50
Close modal
Fig. 14
(a) The mesh and (b) the deformation visualization of the FEA model of the 3PRRR manipulator in ANSYS
Fig. 14
(a) The mesh and (b) the deformation visualization of the FEA model of the 3PRRR manipulator in ANSYS
Close modal

Figure 15 shows the deformation of the manipulator's end-effector across its workspace, evaluated in the MSA model. Figure 16 shows a possible useful workspace with minimal deformation that represents the optimal stiffness. The stiffness map shows that the optimal stiffness is achieved at positions far enough from the workspace boundaries. In other words, the stiffness is poor at positions around the workspace boundaries. This can be related to the singularities. Since the singularities of the manipulator occur along the workspace boundaries, the stiffness of the manipulator is poor at singular positions or close to the singular positions.

Fig. 15
The EE deformation (in mm) across (a) the whole workspace from x=y=z=150mm to x=y=z=350mm and (b) a cross section of the workspace at x=350mm
Fig. 15
The EE deformation (in mm) across (a) the whole workspace from x=y=z=150mm to x=y=z=350mm and (b) a cross section of the workspace at x=350mm
Close modal
Fig. 16
The EE deformation map (in mm) at x=350mm, y from 175 mm to 375 mm, and z from 175 mm to 375 mm, with a square indicating the minimal deformations in mm
Fig. 16
The EE deformation map (in mm) at x=350mm, y from 175 mm to 375 mm, and z from 175 mm to 375 mm, with a square indicating the minimal deformations in mm
Close modal

5.4.2 Deformation Against Moments.

Now consider moments Mx=My=Mz=25Nm are applied to the manipulator's end-effector. The moment components Mz, Mx, and My are respectively applied to limbs 1, 2, and 3 in the limb evaluation. The end-effector of the complete manipulator due to all the moment components can be approximated by summing the deformations obtained from the limbs' evaluation, as shown in Table 9. A comparison between the solution using this approach and that using an assembled manipulator model is shown in Table 10.

Table 9

EE deformation in the 3PRRR parallel manipulator due to moments, computed with the proposed approach using the analytical/MSA solution

LimbEnd-effector wrench componentdx(mm)dy(mm)dz(mm)
Limb 1Mz=25Nm−0.328700
Limb 2Mx=25Nm0−0.32870
Limb 3My=25Nm00−0.3287
Combined (Mx=My=Mz=25Nm)−0.3287−0.3287−0.3287
LimbEnd-effector wrench componentdx(mm)dy(mm)dz(mm)
Limb 1Mz=25Nm−0.328700
Limb 2Mx=25Nm0−0.32870
Limb 3My=25Nm00−0.3287
Combined (Mx=My=Mz=25Nm)−0.3287−0.3287−0.3287
Table 10

EE deformation in the 3PRRR parallel manipulator due to moments

Deformation componentDeformation using the proposed approach (mm)Deformation using the assembled manipulator model
MSA
(mm)
ANSYS FEA
(mm)
dx−0.3287−0.4168−0.3046
dy−0.3287−0.4168−0.3046
dz−0.3287−0.4168−0.3046
Deformation componentDeformation using the proposed approach (mm)Deformation using the assembled manipulator model
MSA
(mm)
ANSYS FEA
(mm)
dx−0.3287−0.4168−0.3046
dy−0.3287−0.4168−0.3046
dz−0.3287−0.4168−0.3046

5.4.3 Deformation Against Forces and Moments.

Now consider forces Fx=Fy=Fz=50N and moments Mx=My=Mz=25Nm are applied to the manipulator's end-effector. The evaluation of each limb can be performed by combining the evaluations against forces and moments presented earlier. The end-effector of the complete manipulator due to all the force and moment components can be approximated by summing the deformations obtained from the limbs' evaluation, as shown in Table 11. A comparison between the solution using this approach and that using an assembled manipulator model is shown in Table 12.

Table 11

EE deformation in the 3PRRR parallel manipulator due to forces and moments, computed with the proposed approach using the analytical/MSA solution

LimbEnd-effector wrench componentdx(mm)dy(mm)dz(mm)
Limb 1Fx=50N−0.343200
Mz=25Nm−0.328700
Limb 2Fy=50N0−0.34320
Mx=25Nm0−0.32870
Limb 3Fz=50N00−0.3432
My=25Nm00−0.3287
Combined (Fx=Fy=Fz=50N, Mx=My=Mz=25Nm)−0.6719−0.6719−0.6719
LimbEnd-effector wrench componentdx(mm)dy(mm)dz(mm)
Limb 1Fx=50N−0.343200
Mz=25Nm−0.328700
Limb 2Fy=50N0−0.34320
Mx=25Nm0−0.32870
Limb 3Fz=50N00−0.3432
My=25Nm00−0.3287
Combined (Fx=Fy=Fz=50N, Mx=My=Mz=25Nm)−0.6719−0.6719−0.6719
Table 12

EE deformation in the 3PRRR parallel manipulator due to forces and moments

Deformation componentDeformation using the proposed approach (mm)Deformation using the assembled manipulator model
MSA
(mm)
ANSYS FEA
(mm)
dx−0.6719−0.6201−0.5479
dy−0.6719−0.6201−0.5479
dz−0.6719−0.6201−0.5479
Deformation componentDeformation using the proposed approach (mm)Deformation using the assembled manipulator model
MSA
(mm)
ANSYS FEA
(mm)
dx−0.6719−0.6201−0.5479
dy−0.6719−0.6201−0.5479
dz−0.6719−0.6201−0.5479

5.5 Parametric Studies.

To investigate the influence of the manipulator's kinematic parameters on the manipulator's stiffness, a global stiffness index will be used to measure the stiffness quality at various parameter values. The global stiffness index is the mean of the average end-effector deformations across a specified workspace. This index will be evaluated based on the MSA model of the manipulator.

5.5.1 The Influence of the Links' Lengths.

Figure 17 shows the means of averaged end-effector deformations over a range of links' lengths L1 and L2. The means of averaged end-effector deformations dEE,mean is obtained as
(5)
where dEE,x, dEE,y, and dEE,z indicate the end-effector deformation in the X, Y, and Z directions, respectively. The figure shows that the shorter the links the less the end-effector deformation and accordingly the stiffer the manipulator. The color map indicates that the best stiffness is achieved when both links are shorter.
Fig. 17
Mean of the averaged end-effector deformations (in mm) across the workspace over the links' lengths
Fig. 17
Mean of the averaged end-effector deformations (in mm) across the workspace over the links' lengths
Close modal

5.5.2 The Influence of the Lower and Upper Links' Lengths Ratio.

Figure 18 shows the end-effector deformation changes with the link length ratio L1/L2. The plotted curves indicate the means of averaged end-effector deformations over the changing ratios at some fixed values of L1 or L2 at a time. The dashed curves correspond to the cases in which L1 is fixed at some values while L2 varies from 150 mm to 300 mm. On the other hand, the solid curves correspond to the cases in which L2 is fixed at some values while L1 is varied from 150 mm to 300 mm. The plot shows that the value of the fixed link lengths affects the best ratio, i.e., the ratio resulting in the minimum end-effector deformation. An observation of the whole curves shows that a ratio of 1 is in general a tradeoff.

Fig. 18
Mean of the averaged end-effector deformations across the workspace over the links' lengths
Fig. 18
Mean of the averaged end-effector deformations across the workspace over the links' lengths
Close modal

5.5.3 The Influence of the End-Effector Offset Distance.

In practice, the end-effector cannot be zero. Figure 19 shows that shorter end-effector offset distance leads to less mean averaged end-effector deformation. The deformation shown in the figure is based on modeling the end-effector as three orthogonal beams connected by fixed joints at the end-effector position (x, y, z). All of the beams have an identical length that represents the end-effector offset distance.

Fig. 19
Mean of the averaged end-effector deformations across the workspace over the end-effector offset distance
Fig. 19
Mean of the averaged end-effector deformations across the workspace over the end-effector offset distance
Close modal

6 Conclusions

An investigation on the stiffness of a two-link limb subject to out-of-plane loads was conducted. We proposed an analytical solution of a two-link limb subject to an out-of-plane moment to complement a previously published analytical solution of a two-link limb subject to an out-of-plane force. We verified that the analytical solution is identical to the MSA solution, while both solutions have a slight difference from the FEA solution. Furthermore, this study investigated the effect of the limb angle and the links' lengths on the limb's out-of-plane stiffness. The following can be concluded from the parametric analysis of the two-link limb subject to an out-of-plane loading:

  • The closer the limb is to the folded configuration, the stiffer the limb in the out-of-plane direction. This corresponds to a shorter effective length of the limb.

  • The effect of the link length can be evaluated based on the resulting effective length. The shorter the effective length, the stiffer the limb, and vice versa.

  • A larger link length ratio L1/L2 with constant total length and effective length results in a higher out-of-plane stiffness.

The stiffness evaluation of the two-link limb serves as a base for the subsequent stiffness evaluation of a decoupled Cartesian parallel manipulator. A new approach to quickly estimate the stiffness of the entire manipulator without the need to assemble the limbs of the manipulator shows results close to those obtained by using the common approach that requires the assembly of the whole manipulator. The new approach only evaluates the deformation of each limb independently, followed by superposing the deformations of all the limbs. The new approach can be used utilizing any existing deformation-based stiffness evaluation method such as analytical method, MSA, and FEA. Not only the new approach provides a quick stiffness evaluation as it avoids an assembly process, but it also provides a direct insight into the contribution of the limb stiffness to the entire manipulator stiffness. Furthermore, a parametric analysis was performed on the decoupled parallel manipulator. The parametric analysis shows that shorter links' lengths and end-effector offset result in a stiffer manipulator while a link length ratio of 1 is the best tradeoff for the manipulator stiffness.

Funding Data

  • This work was supported by Khalifa University of Science and Technology under Grant No. RIG-2023-003.

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

References

1.
Jin
,
Y.
,
Chen
,
I.-M.
, and
Yang
,
G.
,
2009
, “
Kinematic Design of a Family of 6-DOF Partially Decoupled Parallel Manipulators
,”
Mech. Mach. Theory
,
44
(
5
), pp.
912
922
.
2.
Kong
,
X.
, and
Gosselin
,
C. M.
,
2004
, “
Type Synthesis of Input-Output Decoupled Parallel Manipulators
,”
Trans. Can. Soc. Mech. Eng.
,
28
(
2A
), pp.
185
196
.
3.
Pinto
,
C.
,
Corral
,
J.
,
Altuzarra
,
O.
, and
Hernández
,
A.
,
2010
, “
A Methodology for Static Stiffness Mapping in Lower Mobility Parallel Manipulators With Decoupled Motions
,”
Robotica
,
28
(
5
), pp.
719
735
.
4.
Rosyid
,
A.
,
Stefanini
,
C.
, and
El-Khasawneh
,
B.
,
2022
, “
A Reconfigurable Parallel Robot for On-Structure Machining of Large Structures
,”
Robotics
,
11
(
5
), p.
110
.
5.
Rosyid
,
A.
,
Stefanini
,
C.
, and
El-Khasawneh
,
B.
,
2023
, “
A Novel Walking Parallel Robot for On-Structure Three-Axis Machining of Large Structures
,”
ASME J. Mech. Rob.
,
15
(
6
), p.
061016
.
6.
Gosselin
,
C. M.
,
1990
, “
Stiffness Mapping for Parallel Manipulators
,”
Trans. Rob. Autom.
,
6
(
3
), pp.
377
382
.
7.
Xi
,
F.
,
Zhang
,
D.
,
Mechefske
,
C. M.
, and
Lang
,
S. Y. T.
,
2004
, “
Global Kinetostatic Modelling of Tripod-Based Parallel Kinematic Machine
,”
Mech. Mach. Theory
,
39
(
4
), pp.
357
377
.
8.
Ceccarelli
,
M.
, and
Carbone
,
G.
,
2002
, “
A Stiffness Analysis for CaPaMan (Cassino Parallel Manipulator)
,”
Mech. Mach. Theory
,
37
(
5
), pp.
427
439
.
9.
Zhang
,
D.
,
Xi
,
F.
,
Mechefske
,
C. M.
, and
Lang
,
S. Y. T.
,
2004
, “
Analysis of Parallel Kinematic Machine With Kinetostatic Modeling Method
,”
Rob. Comput.-Integr. Manuf.
,
20
(
2
), pp.
151
165
.
10.
Rosyid
,
A.
,
El-Khasawneh
,
B.
, and
Alazzam
,
A.
,
2018
, “
Genetic and Hybrid Algorithms for Optimization of Non-Singular 3PRR Planar Parallel Kinematics Mechanism for Machining Application
,”
Robotica
,
36
(
6
), pp.
839
864
.
11.
Li
,
Y.
, and
Xu
,
Q.
,
2008
, “
Stiffness Analysis for a 3-PUU Parallel Kinematic Machine
,”
Mech. Mach. Theory
,
43
(
2
), pp.
186
200
.
12.
Rosyid
,
A.
,
Stefanini
,
C.
, and
El-Khasawneh
,
B.
,
2024
, “
A Novel Reconfigurable 3-DOF Parallel Kinematics Machine
,”
ASME J. Mech. Rob.
,
16
(
2
), p.
021005
.
13.
Pashkevich
,
A.
,
Chablat
,
D.
, and
Wenger
,
P.
,
2009
, “
Stiffness Analysis of Overconstrained Parallel Manipulators
,”
Mech. Mach. Theory
,
44
(
5
), pp.
966
982
.
14.
Joshi
,
S. A.
, and
Tsai
,
L.
,
2002
, “
Jacobian Analysis of Limited-DOF Parallel Manipulators
,”
ASME J. Mech. Des.
,
124
(
2
), pp.
254
258
.
15.
Wang
,
Y.
,
Liu
,
H.
,
Huang
,
T.
, and
Chetwynd
,
D. G.
,
2009
, “
Stiffness Modeling of the Tricept Robot Using the Overall Jacobian Matrix
,”
ASME J. Mech. Rob.
,
1
(
2
), p.
021002
.
16.
Hoevenaars
,
A. G.
,
Lambert
,
P.
, and
Herder
,
J. L.
,
2016
, “
Jacobian-Based Stiffness Analysis Method for Parallel Manipulators With Non-Redundant Legs
,”
Proc. Inst. Mech. Eng., Part C: J. Mech. Eng. Sci.
,
230
(
3
), pp.
341
352
.
17.
Huang
,
T.
,
Liu
,
H. T.
, and
Chetwynd
,
D. G.
,
2011
, “
Generalized Jacobian Analysis of Lower Mobility Manipulators
,”
Mech. Mach. Theory
,
46
(
6
), pp.
831
844
.
18.
Klimchik
,
A.
,
Pashkevich
,
A.
, and
Chablat
,
D.
,
2019
, “
Fundamentals of Manipulator Stiffness Modeling Using Matrix Structural Analysis
,”
Mech. Mach. Theory
,
133
, pp.
365
394
.
19.
Soares Júnior
,
G. D. L.
,
Carvalho
,
J. C. M.
, and
Gonçalves
,
R. S.
,
2016
, “
Stiffness Analysis of Multibody Systems Using Matrix Structural Analysis—MSA
,”
Robotica
,
34
(
10
), pp.
2368
2385
.
20.
Klimchik
,
A.
,
Pashkevich
,
A.
, and
Chablat
,
D.
,
2013
, “
CAD-Based Approach for Identification of Elasto-Static Parameters of Robotic Manipulators
,”
Finite Elem. Anal. Des.
,
75
, pp.
19
30
.
21.
Taghvaeipour
,
A.
,
Angeles
,
J.
, and
Lessard
,
L.
,
2010
, “
Online Computation of the Stiffness Matrix in Robotic Structures Using Finite Element Analysis
,”
Department of Mechanical Engineering and Centre for Intelligent Machines, McGill University
,
Montreal, Canada
.
22.
Cammarata
,
A.
,
2020
, “
Full and Reduced Models for the Elastodynamics of Fully Flexible Parallel Robots
,”
Mech. Mach. Theory
,
151
, p.
103895
.
23.
Liu
,
X-J;
,
Jin
,
Z.-L.
, and
Gao
,
F.
,
2000
, “
Optimum Design of 3-DOF Spherical Parallel Manipulators With Respect to the Conditioning and Stiffness Indices
,”
Mech. Mach. Theory
,
35
(
9
), pp.
1257
1267
.
24.
Raoofian
,
A.
,
Taghvaeipour
,
A.
, and
Kamali
,
A.
,
2018
, “
On the Stiffness Analysis of Robotic Manipulators and Calculation of Stiffness Indices
,”
Mech. Mach. Theory
,
130
, pp.
382
402
.