Next Article in Journal
An Integrated Application of Motion Sensing and Eye Movement Tracking Techniques in Perceiving User Behaviors in a Large Display Interaction
Previous Article in Journal
Lean Manufacturing in Industry 4.0: A Smart and Sustainable Manufacturing System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Constraints and Parasitic Motions of a Tripod Parallel Continuum Manipulator

1
Department of Mechanical Engineering, University of the Basque Country UPV/EHU, 48013 Bilbao, Spain
2
Department of Mechanical and Aerospace Engineering, Politecnico di Torino, 10129 Torino, Italy
3
Mechanical Systems Design Laboratory, Tokyo Institute of Technology, Tokyo 152-8550, Japan
*
Author to whom correspondence should be addressed.
Machines 2023, 11(1), 71; https://doi.org/10.3390/machines11010071
Submission received: 17 November 2022 / Revised: 26 December 2022 / Accepted: 28 December 2022 / Published: 6 January 2023
(This article belongs to the Section Machine Design and Theory)

Abstract

:
A parallel continuum manipulator (PCM) is a mechanism of closed-loop morphology with flexible elements such that their deformation contributes to its mobility. Flexible hexapods are six-degrees-of-freedom (DoF) fully parallel continuum mechanisms already presented in the literature. Devices of reduced mobility, i.e., lower mobility than six DoF, have not been studied so far. An essential characteristic of lower mobility mechanisms is that reduced mobility is due to kinematic constraints generated by mechanical arrangements and passive joints. In rigid-link parallel manipulators, those constraints are expressed as a set of equations relating to the parameters representing the end effector’s pose. As a consequence, independent output pose variables are controllable with the position equations, while dependent output variables undergo parasitic motions. In this paper, the performance of a tripod-type parallel continuum manipulator, 3 P F ̲ S , is compared with the operation of its rigid counterpart 3 P ̲ R S . We will show that in PCMs there are no such geometric constraints expressible with algebraic equations, but it is difficult to perform some types of motion in the end effector with the input torques. Another goal of this paper is to evaluate such limitation of motion in a tripod-like PCM and compare it with the constraints of the rigid 3 P ̲ R S . Finally, the paper shows that there are strong similarities in the reduced mobility of both mechanisms.

1. Introduction

Parallel kinematic mechanisms (PKMs) are closed-loop mechanisms with rigid links and kinematic joints that connect an end effector to a fixed frame via several kinematic chains with actuators distributed on them to achieve better stiffness and accuracy [1]. The hexapod is a typical mechanism of this type with full mobility (six DoFs). Lower mobility PKMs are those that have less than six DoFs. Such mechanisms are subjected to a set of permanent geometrical constraints that condition some output motions. In some cases, those geometric constraints impose a certain constant value on several of the output parameters, as in a Delta robot with translational three-DoF motion. In other cases, the geometric constraints generate a set of equations that relate some output parameters to others, as in tripod parallel kinematic machines of the type 3 P ̲ R S chosen in this paper to illustrate the subject, where two translational output coordinates and one orientation angle are dependent on the other two orientation angles and the other one translational coordinate. The type of motion resultant is a 2 R 1 T mixed-freedom motion [2].
The constraint analysis is a key task of the kinematic analysis of a mechanism because it determines the type of motion allowed at the end effector, and then, it shows the suitability of the morphology chosen to accomplish the objective motions of the machine. However, it is also fundamental in order to define the inverse kinematic (IK) position problem. This is because we need to know which are the independent output pose variables to be used to control the IK problem and which are the functions that generate the values of the dependent output variables. With that information, the output pose is completely defined and the IK problem can be solved using the loop-closure equations of each limb. In the case of the tripod chosen for this paper, 3 P ̲ R S , these are two angles and one coordinate, while the dependent output parameters acquire values along motions that are called parasitic motions [3,4]. The permanent geometrical constraints imposed are due to passive kinematic pairs and rigid links not controlled by the actuator system.
Compliant mechanisms are a group of mechanisms that acquire mobility thanks to the relative flexibility of some of their parts [5,6,7]. A subset of them is the so-called continuum mechanisms, these are slender systems with a flexible backbone actuated using parallel assemblies of tendons. They are motivated by the problem of manipulation in confined, hard-to-reach workspaces [8], and suitable for minimally invasive procedures due to their dexterity and ease of miniaturization [9,10]. An evolution of the aforementioned systems is parallel continuum mechanisms (PCMs). These are flexible mechanical devices with a rigid end effector that is connected to a fixed frame using flexible slender links whose nonlinear deformation is the cause of its mobility [11]. An appropriate model for the nonlinear deformation of the rods used as elements is the Cosserat rod model [12,13]. This model is used extensively, and it can be mathematically expressed in different ways, a nonlinear system of differential equations (ODEs) being the most accurate one. The closed-loop system of flexible rods and rigid elements generates a coupling of large deformation modeling equations and force equilibrium ones to accommodate the deformation’s internal forces and moments to the load. A sustained line of research on hexapod-like flexible devices has been performed, solving a variety of problems such as real-time position control and elastic analysis [14,15,16,17,18,19].
In the case of parallel continuum mechanisms, permanent constraints to reduce the full mobility as in the above-mentioned tripod are not feasible because the flexible links can deform under actuation and load reaching any configuration in the space. However, some mechanical arrangements can introduce a much higher limitation of deformation in some directions, producing a similar constraining effect. Nevertheless, the resultant type of motion can no longer be obtained with an algebraic manipulation of some geometric constraints, and the expected relationships between output parameters cannot be obtained in closed form.
In this paper, we are aiming to show that PCMs of lower mobility exist, they have strong similarities to their rigid-links counterparts, their basic kinematic principles apply equally, and whether numerical procedures have to be devised to perform their kinematic constraint analysis. In rigid lower-mobility parallel manipulators, the kinematic constraints that reduce their mobility are defined by a set of equations. In the case of parallel continuum mechanisms, these geometric constraints expressible with algebraic equations do not exist. In this paper, we will try to propose an analytical expression that related the output independent parameters to the parasitic motion. To better understand this and other aspects of the kinematics of this type of mechanism, we propose the analysis of a flexible 3 P F ̲ S parallel continuum tripod model, and introduce some changes to the above-mentioned flexible hexapod already studied in the literature. Then, a numerical analysis will be used to solve the forward kinematics (FK) position problem following the procedure in [14] with some modifications. By introducing any set of inputs in that FK problem, it is possible to obtain the possible output motion, to compare it to the rigid-link tripod 3 P ̲ R S , and to deduce from those numerical results the type of motion and the parasitic motions.

2. Method

The goal of this research is to verify if tripod PCMs have a type of motion, constraint equations, and parasitic motions that match the results in their rigid-links PKM counterparts, and if the constraint analysis on the latter applies to the former. In order to achieve the goal of this research we will follow the following procedure.
First, we choose a classical lower mobility PKM with mixed freedoms, the rigid-link 3 P ̲ R S tripod (see Figure 1), as a way of example to demonstrate the relevance of the problem of constraint analysis in closed-loop spatial mechanisms.
The family of tripod-like parallel mechanisms is a paradigmatic example because it fulfills several conditions: it has a high number of constraints, the constraints are coupled so that output motions are related between them, it has symmetry, it has many different variants, and it has been extensively studied in the literature using many different kinematic methods. This type of spatial mechanism has reduced mobility to 3 degrees of freedom, involving three geometrical constraints; and very often the resultant type of motion on the end effector is a 2 R 1 T (2 rotational and 1 translational). This means that the constraint equations impose a set of equations among the output parameters of the end effector such that the other 2 translations and the 1 rotation parameters are dependent on the independent ones, and hence, as motion is produced, we acquire values that are called parasitic motions. For our purpose, we have chosen the symmetric geometric arrangement shown in Figure 1. The three limbs are placed in vertical planes π i 120 apart. The linear actuators are vertical and fixed to the base at points A i . The end effector is an equilateral triangle with barycenter at P and circumradius r and spherical joints at its vertices B i join with rigid links of length L attached to revolute pairs at C i . These are oriented with vectors u i normal to π i .
The constraint analysis of this PKM is done using Euler–Rodrigues parameters to express the orientation of the end effector. This redundant parametrization is an alternative to the use of any set of three angles (Euler angles for example). Redundant parametrizations avoid the use of trigonometric functions in the analysis and generate polynomial equations for easier analysis. By using the geometric conditions imposed by passive revolute joints on the feasible output motion of the end effector, it is possible to get a closed-form solution of the equations that relate output parameters between them. That result shows the type of output motion of the mechanism, and the constrained motion of the dependent output parameters, i.e., the parasitic motion. These results serve as reference comparisons for the kinematic analysis of the corresponding flexible mechanism.
Second, we devise a lower mobility PCM with 3 DoF looking for an analogy to the rigid-link tripod. The 6-degrees-of-freedom flexible hexapod shown in Figure 2, where the end effector is connected to 6 flexible rods whose length is controlled by actuators, has been studied in [14]. There, the FK and IK problems are solved using the Cosserat rod model for the nonlinear deformation of flexible rods, and the minimization of a set of geometric and force equilibrium conditions. It requires a known home position configuration so that contiguous configurations are solved with a boundary value problem upon slight changes of either input or output variables.
The hexapod mechanism in Figure 2 can be modified to get the tripod-like PCM shown in Figure 3. The end of the rods ( 2 j 1 ) and 2 j ( j = 1 , 2 , 3 ) are joined to the end effector at coincident points B ( 2 j 1 ) and B 2 j with a single spherical joint that introduces no restriction of the rod’s self-rotation. The same rods ( 2 j 1 ) and 2 j are connected to the base at fixed points A ( 2 j 1 ) and A ( 2 j ) distributed symmetrically at circumradius r A and with a fixed vertical orientation, again with no restriction of intrinsic rotation so that torsion effects on the rods are avoided. At the connection to the fixed base, rods are conducted through said guiding holes A i to the 3 linear actuators j below the base that control the length of the rods ( 2 j 1 ) and 2 j together ( j = 1 , 2 , 3 ). Rods are free to undergo deformations between the base and the end effector, being subjected only to bending, not torsion, shear, or axial deformations. Due to the arrangements of the flexible rods, the bending of each pair from A ( 2 j 1 ) and A ( 2 j ) to B ( 2 j 1 ) = B ( 2 j ) is easier about an horizontal axis of direction A ( 2 j 1 ) A ( 2 j ) , being analogous morphology to the P ̲ R S kinematic chain used in the above-mentioned tripod PKM. The resultant PCM is denoted 3 P F ̲ S , where F indicates the bending flexibility of the limb.
The following step is to solve the position problem, i.e., to evaluate the relationship between the location of the end effector platform under some load and the effective length L i of each pair of flexible rods. By effective length we mean the longitude of the rod that can deform, i.e., the one between the attachments to fixed frame A i and the join to end effector at B i . The location of the end effector platform is depicted through the combination of the position vector of reference point P in the fixed frame, p , and the orientation of the moving frame with respect to the fixed frame given by the rotation matrix R E E . The FK problem determines the pose of the end effector, i.e., the position vector p and the rotation matrix R E E when each rod i takes a certain value of its effective length L i and a given load is imposed. Its solution is straightforward using the procedure in [14] because the tripod PCM works, in this case, as a particular type of hexapod morphology actuated with 6 input lengths that are equal two by two L ( 2 j 1 ) = L ( 2 j ) ( j = 1 , 2 , 3 ).
However, the IK problem cannot be solved straightforwardly because to start the process we need to introduce a desired output pose of the end effector. For that, we need to know which are the independent parameters that can be used and which is the constrained value of the dependent ones. Again, as in classical PKMs we need to solve a constraint analysis first. In this paper, we will check if the constraint analysis performed for the rigid-link PKM is valid for the PCM, and if we can use the same independent output variables and the same constraint equations.
For that purpose, the final step in this study is to introduce a variety of input values in the FK solver of the tripod PCM, obtain the output pose parameters, and verify the degree of accomplishment of the constraint equations obtained for the tripod PKM.

3. Results

In this section we will elaborate on the above-mentioned methodology: first, to deduce the constraint equations applicable and the parasitic motions produced in the rigid-link parallel mechanism 3 P ̲ R S ; second, to perform the kinematic analysis of the analogous flexible parallel continuum mechanism; third, to compare and discuss the fulfillment of 3 P ̲ R S constraint equations with the simulated output motions of the 3 P F ̲ S .

3.1. Constraint Analysis of a 3 P ̲ R S Tripod using Euler–Rodrigues Parameters

The rigid-link parallel kinematic mechanism chosen for the discussion is a tripod-like mechanism with the 3 P ̲ R S morphology, as shown in Figure 1. The mechanism is formed by a rigid end effector joined by spherical joints to rigid links of constant length L, that are joined to vertical linear actuators through revolute joints of horizontal axes. A moving frame is attached to the end effector at reference point P, whose position in a fixed frame is p = [ x , y , z ] T . The orientation of the moving frame with respect to the fixed frame, i.e., orientation of the end effector, is defined either using a rotation matrix R E E or Euler–Rodrigues parameters (a unit quaternion e ˜ ). The position vector of spherical joints’ attachment points B i expressed in the end effector’s frame is r i . The position vector of the actuators’ attachment point A i expressed in the fixed frame attached to the base at origin O is a i . Revolute joints at points C i have horizontal axis u i .
Three loop-closure equations at B i can be stated through the three limbs of the parallel kinematic mechanism as:
p + R E E r i = a i + ( c i a i ) + ( b i c i ) i = 1 . . . 3
As the passive revolute joint of each limb constrains its corresponding point B i to be on a vertical plane π i perpendicular to u i through C i , which also contains points O and A i , we get the set of three constraint equations by dot-multiplying Equation (1) by vectors u i , namely:
( p + R E E r i ) · u i = 0 i = 1 . . . 3
Using Euler–Rodrigues parameters ( e 0 , e 1 , e 2 , e 3 ) to express the rotation matrix allows a polynomial representation of output orientation, namely:
R E E = e 0 2 + e 1 2 e 2 2 e 3 2 2 ( e 1 e 2 e 0 e 3 ) 2 ( e 1 e 3 + e 0 e 2 ) 2 ( e 1 e 2 + e 0 e 3 ) e 0 2 e 1 2 + e 2 2 e 3 2 2 ( e 2 e 3 e 0 e 1 ) 2 ( e 1 e 3 e 0 e 2 ) 2 ( e 2 e 3 + e 0 e 1 ) e 0 2 e 1 2 e 2 2 + e 3 2
subject to:
e 0 2 + e 1 2 + e 2 2 + e 3 2 = 1
where the Euler–Rodrigues parameters are related to the Euler pole a and the angle rotated ϕ in the following way:
e 0 = cos ( ϕ / 2 )
e 1 e 2 e 3 = sin ( ϕ / 2 ) · a
If the symmetrical arrangement of the system shown in Figure 1 is considered, we have an equilateral triangular end effector with a circumradius r so that in the moving frame we get:
r 1 = 0 r 0 r 2 = 3 r / 2 r / 2 0 r 3 = 3 r / 2 r / 2 0
and revolute axes oriented in the fixed frame as:
u 1 = 1 0 0 u 2 = 1 / 2 3 / 2 0 u 3 = 1 / 2 3 / 2 0
Upon substitution of Equations (3), (7), and (8) into Equation (2), and considering Equation (4), we get the following constraint equations:
e 0 e 3 = 0 x + 2 · r · e 1 e 2 = 0 y + r · ( e 1 2 e 2 2 ) = 0 e 0 2 + e 1 2 + e 2 2 + e 3 2 = 1
As a consequence, for the kinematic analysis of this tripod mechanism we can choose z , e 1 , e 2 as the three independent output pose parameters, with dependent parameters x , y , e 0 , e 3 determined from the constraint equations Equation (9). Moreover, as e 0 = cos ( ϕ / 2 ) is only null when there is no output rotation the constraint equations can be expressed as:
e 3 = 0 x + 2 · r · e 1 e 2 = 0 y + r · ( e 1 2 e 2 2 ) = 0 e 0 2 + e 1 2 + e 2 2 = 1
Dependent parameters x , y , e 0 acquire non-null values as the end effector moves; these are called parasitic motions. We can plot the workspace of the mechanism in the space z , e 1 , e 2 . Moreover, as dependent parameters are expressed solely in terms of e 1 , e 2 in constraint equations Equation (10), parasitic motions for any possible motion of the mechanism can be plotted in the space [ e 1 , e 2 ] where e 1 2 + e 2 2 1 as shown in Figure 4 for r = 1 . This result of the constraint analysis is proof of its reduced mobility of three degrees of freedom. Moreover, the analysis of the effect of the three geometrical constraints justifies that the type of motion on the end effector is a 2 R 1 T (two rotational and one translational).
The kinematic analysis of the relationship between the actuators’ input lengths ρ i and the independent output parameters z , e 1 , e 2 , using the loop-closure equations (Equation (1)), provides the solution to either the forward (FK) or the inverse (IK) position problem. First, we get rid of passive variables by finding the Euclidean norm of ( b i c i ) in Equation (1), namely:
( b i c i ) = p + R E E r i a i ( c i a i ) i = 1 . . . 3 L = p + R E E r i a i ( c i a i ) i = 1 . . . 3
where a symmetrical arrangement of actuators with A i points at a circumradius r A defines:
a 1 = 0 r A 0 a 2 = 3 r A / 2 r A / 2 0 a 3 = 3 r A / 2 r A / 2 0
and actuators’ inputs are given by:
c 1 a 1 = 0 0 ρ 1 c 2 a 2 = 0 0 ρ 2 c 3 a 3 = 0 0 ρ 3
Finally, we introduce constraint equations from Equation (10) into Equation (11). The resultant system has three equations in three unknowns: z , e 1 , e 2 for the FK and ρ 1 , ρ 2 , ρ 3 for the IK.
Velocity analysis is performed on the derivative of the loop-closure equations Equation (1), namely
p ˙ + ω × r i = ρ ˙ i k + ω i × ( b i c i ) i = 1 3
where ω is the angular velocity of the end effector, and ω i is the angular velocity of the rigid connecting link ( C i B i ) .
Dot-multiplying Equation (14) by ( b i c i ) , we get rid of ω i , obtaining a relationship between velocity inputs and velocity outputs:
( b i c i ) · p ˙ + ( r i × ( b i c i ) ) · ω = ρ ˙ i [ ( b i c i ) · k ] i = 1 3
and dot-multiplying Equation (14) by u i , we get the relationship between derivatives of output independent and dependent variables, i.e., the derivative of the constraint equations:
u i · p ˙ + ( r i × u i ) · ω = 0 i = 1 3
Upon assembly of Equations (15) and (16), we get the Jacobian equation:
( b i c i ) r i × ( b i c i ) u i r i × u i p ˙ ω = ( b 1 c 1 ) · k 0 0 0 ( b 2 c 2 ) · k 0 0 0 ( b 3 c 3 ) · k 0 0 0 0 0 0 0 0 0 ρ ˙ 1 ρ ˙ 2 ρ ˙ 3
With those Jacobians, we can perform the singularity analysis to determine the limits of the position analysis. IK singularity is determined by the conditions for rank deficiency in:
J I K = ( b 1 c 1 ) · k 0 0 0 ( b 2 c 2 ) · k 0 0 0 ( b 3 c 3 ) · k
which developed in terms of output variables result:
1 3 e 1 2 + e 2 2 = r A ± L r 1 + 2 3 e 1 e 2 2 e 2 2 = r A ± L r 1 2 3 e 1 e 2 2 e 2 2 = r A ± L r
As a result, the space of parasitic motions shown in Figure 4 will be limited by IK singularity loci (see Figure 5 with r A = 0.25 , L = 0.35 , and r = 0.125 ).
The aforementioned Jacobians can be analyzed on this workspace to obtain dexterity indices, as shown in Figure 6 and Figure 7.
Nevertheless, the parasitic motion functions do not change, and the mobility analysis can be fully defined with the constraint equation (Equation (10)) and the plots in Figure 4. These will be used in the comparison with the motion of the flexible tripod mechanism.

3.2. Forward Kinematic Analysis of a 3 P F ̲ S Parallel Continuum Tripod

The second step in the method is to analyze the flexible 3 P F ̲ S parallel continuum tripod shown in Figure 3 and to check the type of motion performed as inputs are introduced. As before, a moving frame is attached to the end effector at reference point P, whose position in a fixed frame is p . The orientation of the moving frame with respect to the fixed frame, i.e., the orientation of the end effector, is defined using a rotation matrix R E E in terms of Euler–Rodrigues parameters (i.e., a unit quaternion e ˜ ). The position vector of rod attachment point B i expressed in the end effector frame is r i . The position vector of rod attachment point A i expressed in the fixed frame attached to the base at origin O is a i .
The forward kinematics problem consists of determining the pose of the end effector, i.e., the position vector of reference point P, p , and the orientation given by the rotation matrix R E E , when each rod i takes a certain value of its effective length L i and a given load F e x t ; M e x t is imposed.
The nonlinear large deformation of the flexible rods is analyzed with the Kirchhoff model [12]. Each cross-section of the rod i along the arc length s is located with a vector pointing its centroid p i ( s ) and a local frame oriented with a unit quaternion q ˜ i ( s ) (see Figure 8). The rods’ internal moments m ( s ) are related to the rods’ curvature u ( s ) along length s, or vice-versa, with the material constitutive law given by u = K B T 1 R T m , where K B T is a stiffness matrix for bending and torsion; no extension–compression or shear effects are considered in this model to affect internal forces n ( s ) .
The change of shape of the flexible rod and the equilibrium of internal forces and moments with the load along s are related through a system of differential equations d y d s = f that can be stated with the vector y of variables
y = p q ˜ n u
and the vector of functions f :
f = R e 3 1 2 q ˜ u ˜ 0 K B T 1 u ^ K B T u + e ^ 3 R T n
Considering the expression of the rotation matrix that defines the orientation of each cross section along the arc length s in terms of a unit quaternion:
R ( s ) = q 0 2 + q 1 2 q 2 2 q 3 2 2 ( q 1 q 2 q 0 q 3 ) 2 ( q 1 q 3 + q 0 q 2 ) 2 ( q 1 q 2 + q 0 q 3 ) q 0 2 q 1 2 + q 2 2 q 3 2 2 ( q 2 q 3 q 0 q 1 ) 2 ( q 1 q 3 q 0 q 2 ) 2 ( q 2 q 3 + q 0 q 1 ) q 0 2 q 1 2 q 2 2 + q 3 2
and the change of the quaternion along s with u
d q ˜ d s = 1 2 q ˜ u ˜ = 1 2 0 u x u y u z u x 0 u z u y u y u z 0 u x u z u y u x 0 q 0 q 1 q 2 q 3
where u z = 0 everywhere because the flexible rods used in the mechanism are free of torsional effects.
We get the following system of differential equations:
d d s x y z q 0 q 1 q 2 q 3 n x n y n z u x u y = 2 ( q 1 q 3 + q 0 q 2 ) 2 ( q 2 q 3 q 0 q 1 ) q 0 2 q 1 2 q 2 2 + q 3 2 1 2 ( u x q 1 u y q 2 ) 1 2 ( u x q 0 u y q 3 ) 1 2 ( u y q 0 + u x q 3 ) 1 2 ( u y q 1 u x q 2 ) 0 0 0 1 E I ( 2 n x ( q 1 q 2 q 0 q 3 ) + n y ( q 0 2 q 1 2 + q 2 2 q 3 2 ) + 2 n z ( q 2 q 3 + q 0 q 1 ) ) 1 E I ( n x ( q 0 2 + q 1 2 q 2 2 q 3 2 ) 2 n y ( q 1 q 2 + q 0 q 3 ) 2 n z ( q 1 q 3 q 0 q 2 ) )
As it can be observed, with no load applied along the arc length of the rod, the internal force n is constant along the rod, while the evolution of p i ( s ) , q ˜ i ( s ) , and u x i ( s ) , u y i ( s ) can be obtained upon integration of the system of equations using, for example, the Runge–Kutta method, from s = 0 to s = L i .
Hence, solving the position problem of the mechanism is started with the integration of each rod numerically by Runge–Kutta from known data values at s = 0 (the position of the base tip of the rod p i ( s = 0 ) and orientation of the base tip q ˜ i ( s = 0 ) , and given inputs L i ), and guess values of u x i ( s = 0 ) , u y i ( s = 0 ) , n i ( 0 ) . To reach convergence with certain security, starting from a known home pose helps to initialize the solving of the problem because those values serve as guess values to find the next solution with the new input values of actuators chosen close enough to secure convergence. The step-by-step solution will lead from that home pose to the required one.
As a result of this iteration, we get for each rod its end tip location p i ( s = L i ) and curvature components u x i ( s = L i ) , u y i ( s = L i ) , to verify the conditions imposed by the mechanism assembly and load. This is done by expressing those conditions in terms of residuals that have to be minimized following an iterative procedure ruled by a Newton scheme.
There are two main conditions that must be satisfied so that a pose of the flexible mechanism is achieved with the data introduced: first, that it can be assembled geometrically; and second, that the configuration is in static equilibrium.
Regarding the first condition, we have the coordinates of the distal end B i as a function of the output pose of the end effector, namely
b i = p + R E E r i
to be related to the coordinates of the end tip of the flexible rod p i ( s = L i ) obtained from the deformation model (as explained before) and the known and fixed position a i of the proximal end A i of each rod:
b i = a i + p i ( s = L i )
so the following set of conditions ( i = 1 6 ) must be fulfilled:
0 = a i + p i ( s = L i ) p R E E r i
Additionally, as distal ends are joined to the end effector through spherical joints, no moment must appear there. Moreover, no torsional effects are present anywhere on the rods due to the type of joint to the base. Bearing in mind that we can easily change from the rods’ internal moments m to the rods’ curvature u along length s, or vice-versa, with the material constitutive law given by u = K B T 1 R T m , the null moment conditions at distal ends ( i = 1 . . . 6 ) to be verified are simply expressed in terms of the components of curvature in the cross-section of the distal ends as:
u x i ( s = L i ) u y i ( s = L i ) = 0 0
Static equilibrium conditions, unlike geometric ones, are stated for the whole device at the same time. This means that the end effector and all rods’ values are highly coupled in the set of equations. The end effector may be subjected to external load F e x t at P and M e x t , and to reaction forces and moments at the attachment to rods due to their deformation, namely internal forces n i and moments m i at ( s = L i ) . So, the following static equilibrium can be stated, taking moments about O:
i = 1 6 [ n i ( L i ) ] F e x t = 0 i = 1 6 [ a i + p i ( L i ) × n i ( L i ) + m i ( L i ) ] p × F e x t M e x t = 0
The problem to be solved is made up of the geometric conditions on the position of distal ends (Equation (27)), i.e., 3 scalar equations at each rod for a total of 18; the geometric conditions on the attachment of distal ends to end effector (Equation (28)), i.e., 2 equations for each rod, 12 in total; and the static equilibrium conditions on the end effector (Equation (29)), i.e., 6 more, to make a total of 36 conditions.
The unknowns are due to the state variables of each rod necessary to integrate the six subsystems of differential equations for each rod, i.e., the unknowns in n and m (alternatively u ) being five because there is no torsion; plus the six unknowns that define the end effector pose p and R E E , a total of thirty-six unknowns.
As the orientation of the end effector with respect to the fixed frame is more effectively expressed with a unit quaternion e ˜ = [ e 0 e 1 e 2 e 3 ] , the corresponding rotation matrix R E E has four components to be used as unknowns instead of three, making a total of thirty-seven unknowns and another condition has to be added, namely the unity condition of the quaternion:
e 0 2 + e 1 2 + e 2 2 + e 3 2 = 1
Hence resulting in a problem of 37 unknowns in 37 conditions.
Due to the integration involved in each rod, the solution must be approached as a boundary value problem: guess values of unknowns at the lower end s = 0 for each rod (namely u i | x y ( 0 ) , n i ( 0 ) ) and the output pose ( p and e ˜ ) are introduced, and a shooting method iterates until residuals of the boundary conditions are below a tolerance.
The chosen order of those residuals will provide an adequate distribution of terms to lower the computational cost. First, the geometric constraints for the rod i = 1 , i.e., Equations (27) and (28), then the same geometric constraints but for the rod i = 2 , and so on until i = 6 . The last elements of the residue vector are formed by the equations of equilibrium (Equation (29)), and the normalization condition of the end effector quaternion (Equation (30)), namely
g r e s = p 1 ( L 1 ) p R E E r 1 + a 1 u x 1 ( L 1 ) u y 1 ( L 1 ) i n e i n e p 6 ( L 6 ) p R E E r 6 + a 6 u x 6 ( L 6 ) u y 6 ( L 6 ) i n e i = 1 6 [ n i ( L i ) ] F e x t i = 1 6 [ a i + p i ( L i ) × n i ( L i ) + m i ( L i ) ] p × F e x t M e x t i n e e 0 2 + e 1 2 + e 2 2 + e 3 2 1
The Newton method used in the shooting method requires the evaluation of a Jacobian of the residue vector function with respect to the variables of the problem, i.e., the guess values u i | x y ( 0 ) , n i ( 0 ) and output pose p and e ˜ = [ e 0 e 1 e 2 e 3 ] , in order to update the guess values accordingly. A special order for the guess values is also considered to generate a Jacobian sparse matrix. They are grouped so that the guess values associated with a rod are together:
v g u e s s = u x 1 g u e s s u y 1 g u e s s n 1 g u e s s i n e i n e u x 6 g u e s s u y 6 g u e s s n 6 g u e s s i n e p g u e s s e ˜ g u e s s
The Jacobian is then obtained as:
J ( v g u e s s ) = g r e s ( v g u e s s ) v g u e s s
The Jacobian related to the forward kinematic problem is a sparse matrix J R 37 x 37 of the form
J = A 1 0 0 0 0 0 C 1 D 1 0 A 2 0 0 0 0 C 2 D 2 0 0 A 3 0 0 0 C 3 D 3 0 0 0 A 4 0 0 C 4 D 4 0 0 0 0 A 5 0 C 5 D 5 0 0 0 0 0 A 6 C 6 D 6 B 1 B 2 B 3 B 4 B 5 B 6 C 7 0 0 0 0 0 0 0 0 D 7
Each sub-matrix A i R 5 x 5 is calculated considering the residue associated with the geometric conditions for the rod i and the kinematic variables used as guess values for the rod i. Similarly, matrix B i R 6 x 5 is defined considering the residue associated with the static conditions that affect the whole mechanism and the kinematic variables used as guess values for rod i. Matrix C i R 5 x 3 is calculated considering the residue associated with the geometric conditions for the rod i and the coordinates of P are used as guess values. Matrix D i R 5 x 4 is calculated considering the residue associated with the geometric conditions for the rod i and the kinematic variables used as guess values for the orientation of the end effector. Finally, matrix C 7 R 6 x 3 is defined considering the residue associated with the static conditions and the output position of P of the end effector; while matrix D 7 R 1 x 4 is defined considering the residue associated to the quaternion normalization condition and the orientation of the end effector.
After each iteration j, guess values are updated with:
v g u e s s j + 1 = v g u e s s j J 1 g r e s j
and the process starts again until g r e s is below a given tolerance.
Starting from a home pose where the end effector is placed horizontally at a certain height, introducing a variation in inputs following the distribution on the jointspace, as shown in Figure 9, and solving the FK position problem using the procedure above, we get the set of output poses represented in the workspace e 1 , e 2 , z P .
As the geometrical arrangement of flexible rods and actuators produces a workspace that is independent of z P as in the 3 P ̲ R S , it is convenient to get the jointspace subset that produces the output poses on a horizontal plane, as shown in Figure 10.
From the solutions found in the above-mentioned analysis, it is possible to do a comparison with those obtained for the rigid link parallel mechanism analyzed in the first part of this work. In the next section some of these preliminary results will be discussed.

4. Discussion

In order to compare the motion capabilities of the flexible tripod 3 P F ̲ S with the rigid link counterpart 3 P ̲ R S , we have chosen to plot together the results of parasitic motions of the dependent output variables e 3 , x P , y P against the independent output parameters e 1 , e 2 , for different z P = h values of the end effector’s initial position. The same jointspace represented in Figure 10a has been explored in this analysis. For the study presented in this section, the data of the flexible rods used in the 3 P F ̲ S manipulator are rod material, N i t i n o l ; Young modulus, 83 GPa; rod diameter, 3 mm.
The output-dependent parameter e 3 is the one better defining the type of motion of the tripod. Its null value in the rigid-link mechanism implies that only rotations about horizontal axes are possible. The comparison plot corresponding to it is shown in Figure 11. The e 3 parameter in the flexible mechanisms has small changes from the null value, so the motion is fairly similar to the rigid tripod. To quantitatively compare the result obtained for the 3 P F ̲ S with the 3 P R S rigid counterpart, the root mean square error ( R M S E ) between the two datasets has been computed for each value h. The results corresponding to the initial heights h = 0.2 m, h = 0.4 m and h = 0.6 m are respectively R M S E e 3 = 0.0132 , R M S E e 3 = 0.0095 , and R M S E e 3 = 0.0065 . Therefore, due to the low value of the error, the motion of the 3 P F ̲ S is similar to the rigid tripod. Moreover, for the flexible mechanism there seems to be some dependence of results with h contrary to the rigid case. As expected, the magnitude of the parasitic motions increases as the length of the rods decreases because of the increase in stiffness.
Output dependent parameters x P a n d y P acquire non-null values for the rigid and flexible tripods with some slight differences, as shown in Figure 12 and Figure 13. Such difference is due to the fact that the mechanical arrangement of flexible rods does not constrain rigidly the spherical joints on the corresponding vertical planes π i . However, these differences do not compromise the type of motion of the mechanism.
To obtain an analytical expression that relates the output independent parameters e 1 and e 2 with the parasitic motions e 3 , x p , and y p , a fitting method based on polynomial approximation has been used. The resulting Equation (36) fits with low error to the blue plots of Figure 11, Figure 12 and Figure 13.
e 3 = 3 c 0 e 1 2 e 2 + c 0 e 2 3 x p = c 3 e 1 4 e 2 + 2 c 2 e 1 3 e 2 4 c 4 e 1 2 e 2 3 + c 3 e 1 2 e 2 + 2 c 2 e 1 e 2 c 4 e 1 e 2 3 + c 1 e 2 + c 3 e 2 3 c 5 e 2 5 y p = 2 c 2 e 1 5 + c 6 5 e 1 4 2 c 4 e 1 3 e 2 2 + c 3 e 1 3 + c 6 e 1 2 e 2 2 + c 2 e 1 2 + 2 c 4 e 1 e 2 4 + c 3 e 1 e 2 2 c 1 e 1 c 5 e 2 4 c 2 e 2 2 e 0 2 + e 1 2 + e 2 2 + e 3 2 = 0
where c 0 = 1.6; c 1 = 0.077 m; c 2 = 0.24 m; c 3 = 0.64 m; c 4 = 0.60 m; c 5 = 0.11 m; c 6 = 0.80 m
The expression (36) represents for the 3 P F ̲ S what Equation (10) does for the 3 P ̲ R S . These relations will be fundamental to evaluating feasible poses of the end effector to be used as inputs of the inverse kinematic problem that will be investigated in future works.
This study has been done for a null applied load. Figure 14, Figure 15 and Figure 16 show the simulation results in terms of output parameters e 1 , e 2 , and z p , and parasitic motions e 3 , x p , and y p under three different load configurations, respectively. Figure 14 is related to a constant centred vertical load F e x t = [ 0 0 F e x t , z ] T , Figure 15 to a constant torque around the x axis of the fixed frame M e x t = [ M e x t , x 0 0 ] T , and Figure 16 to a constant torque around the y axis of the fixed frame M e x t = [ 0 M e x t , y 0 ] T .
In the case of a vertical load, a little contraction of the output parameters is shown. Higher differences from the null load results are observed as the end effector rotates. This behavior indicates higher stiffness of the robot when its orientation is parallel to the home configuration. The parasitic motions under vertical loads are characterized by a scaling effect while their shapes over the e 1 and e 2 workspace remain the same.
Similar considerations can be made in the cases of non-null external torques applied to the end effector. Moreover, in these cases an additional shifting effect is observed. Due to the nonsymmetric arrangement of the robot over the x and y axes of the fixed frame, the effects of torques around these two axes are not equal, as expected.
To quantitatively describe the influence of external load on the end effector’s position, the root mean square error ( R M S E ) between each load configuration and the unloaded one has been computed. To make it easier to read, the RMSE has been normalized over the size of the range of each variable; therefore, the results presented in Table 1 are relative errors between the loaded and unloaded case.
The same analysis has been performed again with a rod diameter of 4 mm to study the effect of the stiffness on the 3 P F ̲ S response to external loads. The results are presented in Table 2. As expected, an increase of stiffness in the system results in lower differences between loaded and unloaded configuration.

5. Conclusions

In this paper, the performance of a tripod-type parallel continuum manipulator, 3 P F ̲ S , is compared with the operation of its rigid counterpart 3 P ̲ R S . On the one hand, in rigid lower mobility parallel manipulators their reduced mobility is expressed as a set of equations relating the parameters representing the pose of the end effector. On the other hand, in continuum architectures, such geometric constraints expressible with algebraic equations are not feasible because flexible links can deform under actuation and load reaching any configuration in space. However, some mechanical arrangements can introduce a much greater limitation of deformation in some directions, producing a similar constraint effect.
The flexible 3 P F ̲ S parallel continuum tripod is modeled by introducing some changes to the flexible hexapod already studied in the literature and the forward kinematics position problem is solved numerically. Introducing a variety of input values in the FK solver of the tripod PCM, we verify the degree of accomplishment of the constraint equations obtained for the tripod PKM. Thus, for the flexible mechanism there seems to be some dependence of the output-dependent parameters on the end effector, contrary to the rigid case. This difference is due to the fact that the mechanical arrangement of the flexible rods does not rigidly constrain the spherical joints in the corresponding vertical planes. However, these differences do not compromise the type of motion of the mechanism. From this analysis, an analytical expression relating the output independent parameters to the parasitic motions is obtained by means of a fitting method. This expression represents for the 3 P F ̲ S what the constraint equations do for the 3 P ̲ R S .
In conclusion, in this paper we show that PCMs of lower mobility exist and that they have strong similarities with their rigid-link counterparts and that their basic kinematic principles are applicable. The obtained analytical expression referred to above will be fundamental to evaluate the feasible end effector poses to be used as inputs for the inverse kinematic problem to be investigated in future work.

Author Contributions

Conceptualization, O.A.; methodology, O.A.; software, O.A., L.T., J.L.R.-E. and Y.L.; validation, O.A., L.T. and V.P.; formal analysis, O.A.; investigation, O.A.; resources, Y.L.; writing—original draft preparation, O.A. and L.T.; writing—review and editing, all authors; supervision, O.A.; funding acquisition, O.A. and V.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish Government through the Ministerio de Ciencia e Innovación (Project PID2020-116176GB-I00) financed by MCIN/AEI/10.13039/501100011033, and the Regional Basque Government through the Departamento de Educación via Research Groups Project IT1480-22.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
PCMParallel Continuum Mechanism
DoFDegree of Freedom
PKMParallel Kinematic Mechanism
2R1TTwo Rotational plus 1 Translational
IKInverse Kinematics position problem
FKForward Kinematics position problem
3PRSThree limbed mechanism with chain Prismatic + Rotational + Spherical
ODEOrdinary Differential Equation
3PFSThree limbed flexible mechanism with chain Prismatic + Flexible rod + Spherical

References

  1. Merlet, J. Parallel Robots; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar] [CrossRef]
  2. Wang, D.; Wu, J.; Wang, L. Research on the error transfer characteristics of a 3-DOF parallel tool head. Robot. Comput.-Integr. Manuf. 2018, 50, 266–275. [Google Scholar] [CrossRef]
  3. Maraje, S.; Nurahmi, L.; Caro, S. Operation Modes Comparison of a Reconfigurable 3-PRS Parallel Manipulator Based on Kinematic Performance. In Proceedings of the International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. ASME, 40th Mechanisms and Robotics Conference, Charlotte, NC, USA, 21–24 August 2016; Volume 5B. [Google Scholar] [CrossRef]
  4. Nurahmi, L.; Caro, S.; Wenger, P. Operation Modes and Singularities of 3-PRS Parallel Manipulators With Different Arrangements of P-Joints. In Proceedings of the International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. ASME, 39th Mechanisms and Robotics Conference, Boston, MA, USA, 2–5 August 2015; Volume 5C. [Google Scholar] [CrossRef]
  5. Howell, L.L. Compliant Mechanisms; John Wiley & Sons: Hoboken, NJ, USA, 2002. [Google Scholar]
  6. McClintock, H.; Temel, F.Z.; Doshi, N.; Koh, J.s.; Wood, R.J. The milliDelta: A high-bandwidth, high-precision, millimeter-scale Delta robot. Sci. Robot. 2018, 3, aar3018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Ruiz, A.; Campa, F.J.; Roldán-Paraponiaris, C.; Altuzarra, O.; Pinto, C. Experimental validation of the kinematic design of 3-PRS compliant parallel mechanisms. Mechatronics 2016, 39, 77–88. [Google Scholar] [CrossRef]
  8. Anderson, P.L.; Mahoney, A.W.; Webster, R.J. Continuum Reconfigurable Parallel Robots for Surgery: Shape Sensing and State Estimation With Uncertainty. IEEE Robot. Autom. Lett. 2017, 2, 1617–1624. [Google Scholar] [CrossRef] [PubMed]
  9. Burgner-Kahrs, J.; Rucker, D.C.; Choset, H. Continuum Robots for Medical Applications: A Survey. IEEE Trans. Robot. 2015, 31, 1261–1280. [Google Scholar] [CrossRef]
  10. Kota, S.; Lu, K.J.; Kreiner, Z.; Trease, B.; Arenas, J.; Geiger, J. Design and Application of Compliant Mechanisms for Surgical Tools. J. Biomech. Eng. 2005, 127, 981–989. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Bryson, C.E.; Rucker, D.C. Toward parallel continuum manipulators. In Proceedings of the Robotics and Automation (ICRA), 2014 IEEE International Conference on, IEEE, Hong Kong, China, 31 May–7 June 2014; pp. 778–785. [Google Scholar]
  12. Antman, S.S. Nonlinear Problems of Elasticity; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  13. Bretl, T.; McCarthy, Z. Quasi-static manipulation of a Kirchhoff elastic rod based on a geometric analysis of equilibrium configurations. Int. J. Robot. Res. 2014, 33, 48–68. [Google Scholar] [CrossRef]
  14. Black, C.B.; Till, J.; Rucker, D.C. Parallel Continuum Robots: Modeling, Analysis, and Actuation-Based Force Sensing. IEEE Trans. Robot. 2018, 34, 29–47. [Google Scholar] [CrossRef]
  15. Rucker, D.C.; III, R.J.W. Statics and Dynamics of Continuum Robots With General Tendon Routing and External Loading. IEEE Trans. Robot. 2011, 27, 1033–1044. [Google Scholar] [CrossRef]
  16. Till, J.; Bryson, C.E.; Chung, S.; Orekhov, A.; Rucker, D.C. Efficient computation of multiple coupled Cosserat rod models for real-time simulation and control of parallel continuum manipulators. In Proceedings of the Robotics and Automation (ICRA), 2015 IEEE International Conference on IEEE, Seattle, WA, USA, 26–30 May 2015; pp. 5067–5074. [Google Scholar] [CrossRef]
  17. Orekhov, A.L.; Aloi, V.A.; Rucker, D.C. Modeling parallel continuum robots with general intermediate constraints. In Proceedings of the 2017 IEEE International Conference on Robotics and Automation (ICRA), Singapore, 29 May–3 June 2017; pp. 6142–6149. [Google Scholar] [CrossRef]
  18. Till, J.; Rucker, D.C. Elastic Stability of Cosserat Rods and Parallel Continuum Robots. IEEE Trans. Robot. 2017, 33, 718–733. [Google Scholar] [CrossRef]
  19. Till, J.; Rucker, D.C. Elastic rod dynamics: Validation of a real-time implicit approach. In Proceedings of the 2017 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Vancouver, BC, Canada, 24–28 September 2017; pp. 3013–3019. [Google Scholar] [CrossRef]
Figure 1. Rigid tripod: 3 P ̲ R S parallel kinematics manipulator.
Figure 1. Rigid tripod: 3 P ̲ R S parallel kinematics manipulator.
Machines 11 00071 g001
Figure 2. Flexible hexapod: 6 P ̲ F R .
Figure 2. Flexible hexapod: 6 P ̲ F R .
Machines 11 00071 g002
Figure 3. Flexible tripod: 3 P F ̲ S .
Figure 3. Flexible tripod: 3 P F ̲ S .
Machines 11 00071 g003
Figure 4. 3 P ̲ R S parasitic motions: (a) x P , (b) y P , (c) e 0 .
Figure 4. 3 P ̲ R S parasitic motions: (a) x P , (b) y P , (c) e 0 .
Machines 11 00071 g004
Figure 5. IK singularity locus for 3 P ̲ R S .
Figure 5. IK singularity locus for 3 P ̲ R S .
Machines 11 00071 g005
Figure 6. IK Jacobians’ norm.
Figure 6. IK Jacobians’ norm.
Machines 11 00071 g006
Figure 7. FK Jacobians’ norm.
Figure 7. FK Jacobians’ norm.
Machines 11 00071 g007
Figure 8. Flexible rod deformed in space.
Figure 8. Flexible rod deformed in space.
Machines 11 00071 g008
Figure 9. 3 P F ̲ S Jointspace transformation to workspace.
Figure 9. 3 P F ̲ S Jointspace transformation to workspace.
Machines 11 00071 g009
Figure 10. 3 P F ̲ S Jointspace transformation to workspace for a constant z P plane.
Figure 10. 3 P F ̲ S Jointspace transformation to workspace for a constant z P plane.
Machines 11 00071 g010
Figure 11. Parasitic motion e 3 comparison.
Figure 11. Parasitic motion e 3 comparison.
Machines 11 00071 g011
Figure 12. Parasitic motion x P comparison.
Figure 12. Parasitic motion x P comparison.
Machines 11 00071 g012
Figure 13. Parasitic motion y P comparison.
Figure 13. Parasitic motion y P comparison.
Machines 11 00071 g013
Figure 14. Influence of vertical centred loads on the 3 P F ̲ S output parameters e 1 , e 2 , and z p (a), and parasitic motions e 3 (b), x p (c), and y p (d).
Figure 14. Influence of vertical centred loads on the 3 P F ̲ S output parameters e 1 , e 2 , and z p (a), and parasitic motions e 3 (b), x p (c), and y p (d).
Machines 11 00071 g014
Figure 15. Influence of constant torque around the x axis of the fixed frame on the 3 P ̲ F S output parameters e 1 , e 2 , and z p (a), and parasitic motions e 3 (b), x p (c), and y p (d).
Figure 15. Influence of constant torque around the x axis of the fixed frame on the 3 P ̲ F S output parameters e 1 , e 2 , and z p (a), and parasitic motions e 3 (b), x p (c), and y p (d).
Machines 11 00071 g015
Figure 16. Influence of constant torques around the y axis of the fixed frame on the 3 P ̲ F S output parameters e 1 , e 2 , and z p (a), and parasitic motions e 3 (b), x p (c), and y p (d).
Figure 16. Influence of constant torques around the y axis of the fixed frame on the 3 P ̲ F S output parameters e 1 , e 2 , and z p (a), and parasitic motions e 3 (b), x p (c), and y p (d).
Machines 11 00071 g016
Table 1. Relative R M S error of variables z p , e 3 , x p , and y p between loaded and unloaded case under three loading configurations: F e x t = [ 0 0 F e x t , z ] T , M e x t = [ M e x t , x 0 0 ] T , and M e x t = [ 0 M e x t , y 0 ] T . Rod diameter 3 mm.
Table 1. Relative R M S error of variables z p , e 3 , x p , and y p between loaded and unloaded case under three loading configurations: F e x t = [ 0 0 F e x t , z ] T , M e x t = [ M e x t , x 0 0 ] T , and M e x t = [ 0 M e x t , y 0 ] T . Rod diameter 3 mm.
F ext , z , (N) RMS z p , (%) RMS e 3 , (%) RMS x p , (%) RMS y p , (%)
−2.500.29940.33160.05870.1847
−5.000.65670.72740.12500.3962
−7.501.09031.20760.20190.6465
−10.001.62761.80270.29500.9569
M e x t , x , (Nm) R M S z p , (%) R M S e 3 , (%) R M S x p , (%) R M S y p , (%)
0.251.08012.91730.08842.1147
0.502.13336.00130.18524.2585
0.753.16619.29410.33306.5019
1.004.186812.84310.55998.9256
M e x t , y , (Nm) R M S z p , (%) R M S e 3 , (%) R M S x p , (%) R M S y p , (%)
0.252.57731.21630.10802.1356
0.505.19382.45090.26224.3486
0.757.89553.72450.50556.7339
1.0010.75125.06380.90969.4405
Table 2. Relative R M S error of variables z p , e 3 , x p , and y p between loaded and unloaded case under three loading configurations: F e x t = [ 0 0 F e x t , z ] T , M e x t = [ M e x t , x 0 0 ] T , and M e x t = [ 0 M e x t , y 0 ] T . Rod diameter 4 mm.
Table 2. Relative R M S error of variables z p , e 3 , x p , and y p between loaded and unloaded case under three loading configurations: F e x t = [ 0 0 F e x t , z ] T , M e x t = [ M e x t , x 0 0 ] T , and M e x t = [ 0 M e x t , y 0 ] T . Rod diameter 4 mm.
F ext , z , (N) RMS z p , (%) RMS e 3 , (%) RMS x p , (%) RMS y p , (%)
−2.500.03480.03850.00700.0219
−5.000.07020.07780.01410.0441
−7.500.10650.11790.02130.0667
−10.000.14340.15890.02860.0896
M e x t , x , (Nm) R M S z p , (%) R M S e 3 , (%) R M S x p , (%) R M S y p , (%)
0.250.13670.35710.01220.2652
0.500.27300.71620.02390.5300
0.750.40871.07750.03520.7943
1.000.54401.44100.04611.0585
M e x t , y , (Nm) R M S z p , (%) R M S e 3 , (%) R M S x p , (%) R M S y p , (%)
0.250.32150.15170.01250.2655
0.500.64310.30350.02510.5312
0.750.96490.45540.03790.7972
1.001.28680.60730.05101.0636
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Altuzarra, O.; Tagliavini, L.; Lei, Y.; Petuya, V.; Ruiz-Erezuma, J.L. On Constraints and Parasitic Motions of a Tripod Parallel Continuum Manipulator. Machines 2023, 11, 71. https://doi.org/10.3390/machines11010071

AMA Style

Altuzarra O, Tagliavini L, Lei Y, Petuya V, Ruiz-Erezuma JL. On Constraints and Parasitic Motions of a Tripod Parallel Continuum Manipulator. Machines. 2023; 11(1):71. https://doi.org/10.3390/machines11010071

Chicago/Turabian Style

Altuzarra, Oscar, Luigi Tagliavini, Yuhang Lei, Victor Petuya, and Jose Luis Ruiz-Erezuma. 2023. "On Constraints and Parasitic Motions of a Tripod Parallel Continuum Manipulator" Machines 11, no. 1: 71. https://doi.org/10.3390/machines11010071

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