Abstract

Performance indexes usually provide global evaluations of robot performances mixing their translational and/or rotational capabilities. This paper proposes a definition of performance index, called direction-selective index (DSI), which has been specifically developed for parallel manipulators and can provide uncoupled evaluations of robot translational capabilities along relevant directions. The DSI formulation is first presented within a general framework, highlighting its relationship with traditional manipulability definitions, and then applied to a family of parallel manipulators (4-RUU) of industrial interest. The investigation is both numerical and experimental and allows highlighting the two chief advantages of the proposed DSIs over more conventional manipulability indexes: not only are DSIs more accurate in predicting the workspace regions where manipulators can best perform translational movements along specific directions, but also they allow foreseeing satisfactorily the dynamic performance variations within the workspace, though being purely kinematic indexes. The experiments have been carried out on an instrumented 4-RUU commercial robot.

1. Introduction

Performance evaluation is one of the most important issues in the analysis and design of manipulators. Indeed, performance indexes may provide useful hints in the design and optimization of robots. In 1982, Salisbury and Craig firstly proposed the condition number of the transpose of the Jacobian matrix as a measure of the workspace quality [1]. Some years later, Yoshikawa introduced the concept of manipulability, based on the Jacobian matrix, in order to find the best postures for manipulators [2]. Subsequently, in [3], Gosselin and Angeles suggested using the global conditioning index in order to evaluate the dexterity of a robot over a given workspace. Such an index is basically an integral mean of the condition number discussed in [1]. In [4], the definitions of manipulability, condition number, and dexterity index have been revisited in order to be applied to parallel robots. Novel formulations have been proposed, which are different from the ones introduced in [13].

Performance indexes based on the Jacobian matrix may be affected by the presence of inhomogeneities in the velocity ratios that may lead to ineffective results. Such a problem was found in [5] in the analysis of the Jacobian matrix of serial manipulators. It was solved by defining a characteristic length of the manipulator and by proposing a weighting positive matrix that allows obtaining a homogenous Jacobian. Unfortunately, such a solution cannot be easily extended to parallel manipulators. The first formulation of a dimensionally homogeneous Jacobian matrix overcoming this problem for planar parallel manipulators was introduced in [6]. Such a Jacobian relates the actuator velocities to the ones of the Cartesian coordinates of two points of the mobile platform. Following such a reasoning, in [7], a new dimensionally homogeneous Jacobian matrix was formulated by relating the velocities of three points (coplanar with the mobile platform joints) to the actuator ones. The same approach has been proposed in [8] for force transmission analysis, and then it has been reformulated in [9] for the dexterity analysis of parallel manipulators.

A simple performance index based on the Jacobian matrix has been also presented in [10]. The proposed index is derived by developing an upper bound on the norm of the rate of change of the Jacobian matrix. Such an upper bound has a significant role in a sufficiency condition for Jacobian rank preservation. It is proved that the index proposed in [9] provides information coherent with manipulability [2] and the condition number [1], while being easier to compute.

In [11], the dexterity index has been exploited together with a space utility ratio in order to compose a mixed performance index for optimizing translational parallel manipulators. The introduction of the space utility ratio index ensures avoiding disproportion between the reachable workspace volume and the structure of the robot by using the ratio between the total workspace volume and the physical size of the manipulator.

A performance investigation on translational parallel manipulators has been presented in [12]. The performance index is given by the minimum and the maximum eigenvalues of the stiffness matrix obtained from the Jacobian as proposed in [13]. The effectiveness of this index has been proved by applying it to an innovative parallel manipulator whose kinematics is analyzed in [14].

Successively, in [15, 16], a kinetostatic performance index, called power manipulability, has been introduced. Power manipulability is related to the power within the mechanism and is fully homogeneous, whether the manipulator contains active joints of different type, or the task space combines both translation and rotation motion.

A different approach for the computation of manipulator dexterity has been proposed in some recent studies [17, 18], where a generalized transmission index that can evaluate the motion/force transmissibility of fully parallel manipulators has been defined through virtual coefficients and not through Jacobian matrix. Contrary to the other indexes, such a transmission index is frame-free (i.e., its value does not depend on the reference frame in which it is computed). Moreover, it can be used to identify a good transmission workspace, where not only is the parallel manipulator effective in terms of motion/force transmission, but also far from its singularities.

All the aforementioned indexes provide global evaluations of robot performances mixing the translational capabilities of robots or the rotational ones. However, in [4], the concept of manipulability has been extended to evaluate the translational and rotational capabilities of parallel manipulators. In such a case, the Jacobian matrix has been split in two submatrices, and , evaluated separately. Following a similar reasoning, in [19], it is proved that it is useful to split the Jacobian matrix in order to extend the manipulability definition further and analyze the translational performance of a robot along a specific direction, which originated the concept of direction-selective performance index (DSI). In [19], DSI formulation has been adopted to study and compare the performances along different directions of a 3-UPU manipulator and a delta-like (4-RUU) manipulator. Generally speaking, when uncoupled DSIs are available for horizontal and vertical movements, they may be usefully adopted to foresee performances on those tasks (e.g., the pick-and-place task), which consist of continuous sequences of vertical and horizontal displacements. Additionally, DSIs may be usefully employed to optimize the location of target frames with respect to a robot, or of a robot within its workcell.

There are two main motivations for this paper:(i)experimentally proving the practical usefulness of DSIs in predicting the regions within the workspace where a manipulator can best perform a specific movement;(ii)showing that DSIs, though being purely kinematic indexes, allow foreseeing dynamic performance variations within the workspace satisfactorily.

In order to prove the practical usefulness of DSIs, in this work, their formulation is applied to foresee performance dissimilarities in the vertical and horizontal movements of an industrial parallel manipulator belonging to the family of four-leg delta-like (4-RUU) manipulators: the Adept Quattro. Additionally, a comparative study is performed between the results obtained by applying DSIs and manipulability. The DSI formulation adopted in this paper is more general than the one first introduced in [19], while several mathematical formulations of manipulability are considered. The analyses made are both numerical and experimental.

The experimental investigation has been done mainly because the numerical results highlight major differences in the best performance regions identified by the DSIs and the manipulability indexes. The experiments have been carried out by instrumenting the robot with an MEMS accelerometer and by making the robot perform extensive tests consisting in constant-length vertical and horizontal movements of the moving platform within the whole workspace: the expected and the real performances of the robot have then been compared. The results from the experimental investigation prove the validity of the hints inferred through the proposed DSI formulation and hence its practical usefulness. In particular, as theoretically expected, DSIs turn up to be far more descriptive than manipulability.

The rest of the paper is set out as follows. Section 2 provides the definition of DSI within a general framework. Then, in Section 3, the DSI definition is applied to the family of 4-RUU parallel manipulators. Section 4 shows the results of a numerical investigation aimed at assessing and comparing the shapes of the DSIs and of manipulability indexes within the workspace of the Quattro robot. Planar and spatial graphical representations of the shapes of the indexes are provided. The experimental campaign and its results are then presented and discussed in Section 5. Finally, Section 6 draws the conclusions.

2. Definition of Direction-Selective Index

The congruence equation of the position analysis of a parallel manipulator can be generally defined as , where denotes the joint variable vector, is the generalized coordinates of the end effector, is an -dimensional implicit function of and , and is the -dimensional zero vector. By differentiating the congruence equation with respect to time, the velocity relation is achieved, where and .

The two matrices and are the Jacobian matrices of, respectively, the end-effector coordinates and the joint variables. Such matrices obviously depend on the configuration of the manipulator. From these equations, there immediately follow the expressions of and and also the solutions of the forward () and of the inverse () velocity kinematic problems.

In parallel robotics, forward position kinematics is usually harder to solve than inverse kinematics. Sometimes the analytical solution of the forward problem is not even available. Hence, the analytical expression of the inverse Jacobian matrix is usually easier to achieve than .

One of the most popular performance indexes is the manipulability index defined by Yoshikawa in [2] as follows: Such an index was conceived for serial robots, hence it is based on the Jacobian matrix of the forward kinematic problem. However, as previously discussed, in parallel manipulators, the inverse Jacobian matrix is usually available, which suggests introducing some small changes to the index definition. By recalling that if are the eigenvalues of a generic square matrix, the eigenvalues of its inverse are , and that the determinant of a matrix is equal to the dot product of its eigenvalues, a definition of the manipulability more suitable to parallel manipulators can be inferred By making use of such a definition, the values of become comparable with those computed through Yoshikawa’s formulation and become identical to such values in case of square Jacobian matrices. Moreover, the values taken by are coherent with the singularity analysis for parallel robots defined in [20].

Other measure criteria of manipulability can be used in the performance analysis of parallel manipulators [4]; manipulability can indeed be also defined by means of(i)the Euclidean norm (or Frobenius norm) ,(ii)the infinity norm ,

where are the elements of the matrix By following a reasoning similar to the one leading to (2) from (1), the following further definitions of manipulability can be provided for parallel manipulators: In order to try discriminating the performances related to linear and rotational displacements of the end effector of a robot, the Jacobian matrix may be split into two parts, and (), which can be evaluated separately [4]. The submatrix contains the ratios between the end-effector translations and the joint displacements; the submatrix contains the ratios between the end-effector rotations and the joint displacements. Hence, different performance indexes can be computed making use of and . Of course, such indexes still provide a global evaluation of the performances mixing either the translational capabilities of robot in the three directions or the rotational ones.

When parallel robots are considered, and are more easily computed through inverse kinematics; consider(i)the matrix , submatrix of comprising the partial derivates with respect to the translational coordinates of the end effector,(ii)the matrix , submatrix of comprising the partial derivates with respect to the rotational coordinates of the end effector.

In the following (4), the matrices and are then the inverse matrices for the translational and the rotational inverse velocity analyses: The analyses based on the Jacobian submatrix and the values of the performance indexes that may be computed from such a matrix are usually affected by the units adopted to measure the end-effector velocities and the joint ones. Some “normalization” criteria hence need to be adopted to make the results of the analyses comparable regardless of the features of the manipulator active joints (see [7, 20]). As a matter of fact, the coefficients of are the ratios between the linear velocities of the end effector and the velocities of the active joints, which might be either translational or rotational. In case of fully translational joints, the analysis of the matrix can be directly performed and leads to dimensionless indexes. Alternatively, if active joints are rotational, in order to keep the performance indexes dimensionless, it may be appropriate to modify the elements of matrix so as to relate the end-effector translational velocities to the tangent velocities of the drive cranks, rather than to their angular velocities. A new dimensionless matrix can therefore be adopted when the active joints are rotational; is obtained by premultiplying by a diagonal matrix made up of the lengths of the drive cranks. By replacing (and its transpose matrix) in (2) and (3), with or (and its transpose matrix), according to the features of the active joints, it is possible to get the manipulability indexes referring to end-effector translational degrees of freedom. In particular, in case of fully rotational active joints, it holds The idea behind direction-selective indexes (DSIs) relies on splitting matrices or into three column vectors in order to extend the manipulability formulation further and make it suitable to analyze robot performances along specific directions. As an example, the following holds in case of fully rotational active joints: DSIs can then be straightforwardly obtained by replacing (and its transpose matrix) in (2) with , , or (and their transpose vectors) The DSIs , and are here proposed to evaluate independently the translational capabilities of a manipulator along the axes of its world reference frame. The performances along any other relevant direction can be straightforwardly evaluated by suitably rotating the world reference frame, so that one axis coincides with the selected direction.

3. DSI Formulation for 4-RUU Parallel Manipulators

In order to assess the capability of the proposed indexes to foresee manipulator performance variations along selected directions, the definitions just provided for , , and have been applied to a generic 4-RUU parallel manipulator.

4-RUU manipulators are four-leg and four-degree-of-freedom parallel manipulators conceived for performing high-speed and high-acceleration pick-and-place operations. The kinematic and dynamic models for this family of parallel robots are reported in [21]. Due to their architecture, these robots can produce Schöenflies motions (i.e., translations along the X, Y and Z direction and a rotation about the Z axis), which are required in most pick-and-place operations, including packaging, picking, packing, and palletizing tasks. An example of 4-RUU manipulator is sketched in Figure 1, which also shows some essential geometrical parameters of the robot, of its workspace, and the robot world reference frame. The robot moving platform takes the shape of a rhombus whose side has length 2d. The platform is connected to the fixed base by four identical RUU serial chain legs, that is, legs with one revolute (R) and two universal (U) joints (see Figure 1).

Let be the end-effector pose, and the active joint coordinates. In order to compute the performance indexes for the studied family of parallel manipulators, since the active joints are all rotational, the following formulation of the translational matrix is used: Matrix is the inverse of the Jacobian matrix and is defined as follows: with: Matrix is the translational part of the Jacobian matrix and is defined as follows: where , , , , and (k = 1,,4).

Finally, matrix takes the following form: , where is the identity matrix of size four and is the crank length.

The DSIs , , and can be computed through (8) by simply extracting the column vectors , , and from ,

4. Numerical Investigation on an Industrial 4-RUU Manipulator: DSI versus Manipulability

A numerical investigation has been carried out by applying the theory developed above to an industrial 4-RUU manipulator (the Adept Quattro). The goal of the investigation is to point out the different suggestions provided by the DSIs , , and and the manipulability indexes, in terms of regions of the workspace where the robot performances should be optimal. In order to achieve exhaustive predictions of the robot performances over its full workspace (shown in Figure 1), it has been necessary to discretize the workspace through a regular and thick grid. The grid has been constructed by selecting 12 equally spaced horizontal planes, located at zu (u = 1,,12) in the robot world reference frame. Half the planes belong to the upper cylindrical region of the workspace, and half to the lower truncated conic region. Each plane contains 11 concentric circles of radius rv (v = 1,…,11). Each circle is split into 36 points with polar angle θw (w = 1,…,36). Hence, the grind comprises 396 points spread through the entire volume of the workspace.

The cylindrical coordinates of each point of the grid have been transformed into Cartesian coordinates in the robot world frame using the following basic expressions: The performance indexes in all the grid points have then been computed by setting , , and . The rotation of the end effector () has instead been set equal to zero in all the grid points.

Table 1 collects the chief geometrical parameters of the Quattro robot and of its workspace (refer to Figure 1 for the meaning of each parameter).

Figures 2 through 7 show the values taken by the manipulability indexes (Figure 2), (Figure 3), and (Figure 4) and by the DSIs (Figure 5), (Figure 6), and (Figure 7) over the robot workspace.

In these figures, a spatial view of the shapes of the indexes is given. So as to simplify the comparison among the predictions of the indexes, regardless of each index values and range of variation, the best performance regions take dark red color in all the figures, while the worst performance regions take dark blue color.

It may be immediately recognized that the manipulability indexes , , and have similar shapes; these indexes depict spherical-like and concentric isomanipulability surfaces. The indexes always take their best values in the inner regions of the workspace, which are hence identified as the best performance regions. The vertical position of such regions depends on the formulation chosen. The index values then decrease in the regions closer to the boundaries of the workspace. The shape shown by the manipulability indexes is rather intuitive; the best manipulator performances seems to be achievable close to the center of the manipulator workspace. However, considerably different hints come from the DSIs , , and, that is, when considering the expected performances for translations along different directions.

Let us first consider and (Figures 5 and 6), whose shapes are nearly identical, apart from a rotation by π/2 radians about the axis. The isomanipulability regions are no longer spherical like and concentric, but rather half-oval in the best performance regions and elliptical like elsewhere. Additionally, the best performance regions are not located at the center of the workspace, but they are close to the boundaries and symmetrical across the axis orthogonal to the selected direction (e.g., orthogonal to axis in the case of ; see Figure 5). The worst performance regions are instead the ones which are furthest from the symmetry axis.

As far as is concerned, Figure 7 shows that the best values of the index are found in a ring like region surrounding the upper part of the workspace. Less satisfactory values are instead computed for the lower part of the workspace.

A better understanding of the shapes of these indexes can be obtained by plotting the values they take on some significant sample planes within the workspace. A few examples are shown in the Figures from 8 through 12, where isomanipulability curves are shown. These figures allow highlighting better the differences between DSIs and manipulability. Since the shapes of the three manipulability indexes are alike (see Figures 2, 3, and 4) and are very different from the ones of the DSIs, in the following analyses and figures only the manipulability index will be compared with the DSIs.

Figures 8, 9, and 10 show the values taken in the plane with coordinate by the manipulability index , and by the DSIs and . Such figures clarify the above-mentioned differences between the best performance region predictions. It is apparent that the use of seems to lead to conclusions that are misleading: for example, the existence of four best performance regions rotated π/2 degrees with respect to each other and is in contrast with the results provided by the DSIs, where two best performance regions rotated π degrees with respect to each other, appears close to the workspace boundaries and symmetrically displaced across the direction orthogonal to the selected direction of motion.

The values taken in the plane with coordinate by and by the DSI are instead plotted in Figures 11 and 12. The manipulability index (Figure 11) reaches its best values in a region located close to the middle of the axis and close to the bottom (i.e., to low values of Z) of the this planar slice of the workspace. The manipulability values decrease in the areas close to the boundaries of the workspace and take the less satisfactory values at the top of the workspace.

Conversely, Figure 12 shows that takes its best values in a curved region crossing the planar slice of the workspace and achieving the very best values at the upper right and left boundaries. The region close to the bottom of workspace is where the less satisfactory values are computed.

5. Experimental Validation

The numerical results shown demonstrate that there are significant discrepancies between the predictions made by manipulability indexes and the DSIs. In order to assess experimentally which indexes provide the best predictions, some relevant experimental tests have been carried out on the Adept Quattro manipulator.

The experiments have been carried out with two objectives:(i)evaluating the translational velocity performances of the robot moving platform within the workspace,(ii)comparing the measured performance variations with the predictions inferred by the indexes analyzed.

Such a comparison makes sense given the apparent dependence on the Jacobian matrix elements of both the indexes and the translational velocities.

The moving platform of the Quattro robot has been instrumented with an uniaxial PCB MEMS accelerometer model 3741 (see Figure 13) alternatively employed to measure accelerations along the X and the Z axes of the robot world reference frame (see Figure 1). The adopted accelerometer has a measuring range of ±30 G and a frequency range between 0 (DC) and 2 kHz. An HBM MGC signal conditioner with ML10 modules has been used to power the MEMS accelerometer and to filter its output signals (a low-pass filter having a Butterworth characteristic and a cut-off frequency of 200 Hz has been employed). The filtered data have been acquired and logged by means of an LMS Pimento analyzer. Finally, experimental data have been processed (high pass filtered and integrated numerically) and compared to the performance index values using Matlab-Simulink.

The experimental tests have been carried out performing several repetitions of either straight horizontal displacements or straight vertical displacements in all the grid points into which the workspace has been split (see the grid definition in Section 4), and where performance indexes have already been computed.

As far as the Adept V+ software program written to run the robot is concerned, it is worth highlighting that(i)the “continuous path” feature of the robot controller has been disabled to avoid modifications of the desired path during the manipulator motion;(ii)the “speed” command has been set equal to 100% (i.e., no speed limitation has been introduced in path planning);(iii)the sharpest acceleration profile (square-wave) has been adopted.

As far as the horizontal (along the X direction) displacements are concerned, around each grid point, the end-effector has been made to perform a sequence of five 30 mm straight forward displacements and five 30 mm straight backwards displacements. The peak-to-peak amplitude of the displacements is 60 mm. The midpoint of these displacements coincides with the grid point considered. The amplitude chosen for the displacements has been selected after some trials to attain two concurrent results: to let the robot achieve relevant velocities at the midpoints and to allow assessing the robot performances in an adequate number of grid points throughout the workspace. The repetition of identical displacements has been motivated by the need of checking the repeatability of the measured accelerations.

Similar tests have been carried out with vertical displacements (along the Z direction): at each grid point, the end effector has been made perform a sequence of six 30 mm straight upwards displacements and five 30 mm straight downward displacements. The peak-to-peak amplitude of the displacements is again 60 mm, and the midpoint of these displacements coincides with the grid point considered. The acceleration data collected during the tests have been processed in order to gather the maximum horizontal and vertical velocity achieved by the end effector at each grid point. Maximum velocities are necessarily reached at the midpoints of the displacements (i.e., at the grid points of interests). As previously mentioned, velocities have been computed from accelerations by numerical integration. Signal detrending has been obtained by filtering the acceleration signals through a second-order Butterworth high-pass filter with passband frequency at 0.2 Hz.

A sample plot of the horizontal velocities computed is shown in Figure 14. The same figure also shows the displacements computed through double integration; the mentioned 60 mm peak-to-peak amplitude of the displacements can be recognized. By comparing the two plots, it can be also verified that maximum velocities are reached at about the midpoints of the paths. The same considerations hold for the vertical velocities (see Figure 15).

For each direction of motion (horizontal or vertical) and for each sequence of displacements (i.e., for each grid point), the very maximum absolute value of velocity achieved during the repetition of the displacements has then been detected and related to the grid point. For example, the subplots (b) of the Figures from 16 through 21 show the velocity values computed for the points belonging to three horizontal half-planes XY and three vertical half-planes YZ. The horizontal half-planes are located at  mm (Figure 16),  mm (Figure 17), and  mm (Figure 18). The vertical half-planes are instead located at  mm (Figure 19),  mm (Figure 20), and  mm (Figure 21). Velocity values have been normalized in order to take 1 as the maximum value.

The Figures from 16 through 21 also show the values taken by the studied performance indexes in the same planes; manipulability index is shown in the subplots (a) of the figures, while the DSIs and are shown in the subplots (c).

It should be observed that in these figures the comparison among numerical and experimental data is restricted to half-planes in consequence of the already proved symmetrical shapes of all the indexes.

It can be immediately recognized that in all cases the proposed DSIs and show shapes very similar to the normalized maximum velocities. Not only do the best performance regions predicted by the DSIs fit well with the experimentally recorded data, but the performance variations foreseen by DSIs seem to be in good agreement with the recorded velocity variations. What is more important, the predictions provided by DSIs are considerably more accurate than the predictions provided by the manipulability index . For example, in the case of horizontal displacements at  mm (Figure 17), suggests an inexistent best performance region at the center left of the half-plane (across the axis), while in the case of vertical displacements, regardless of the plane, it completely misses the best performance region locations.

As final evidence of the usefulness of the suggested DSIs, in the Figures from 22 through 25,  , , and have been related to the normalized maximum velocities achieved at each grid point of the vertical and horizontal half-planes shown in Figures 17 and 19. Generally speaking, though the definitions of the adopted indexes are purely kinematic, it would be informative for good performance indexes to provide satisfactory estimates of dynamic performance variations within the workspace. In other words, a performance index would be practically useful if, for example, an increase of the robot velocity could be predicted through a corresponding increase of the performance index. Figures 22 and 24 show that this is not the case when manipulability is considered; there seems to be no meaningful correlation between and the end-effector velocities in the horizontal and vertical directions. On the contrary, Figures 23 and 25 prove that data are less scattered when the DSIs and are considered; a satisfactory correlation appears between the DSIs and the velocity values.

6. Conclusions

Direction-selective indexes (DSIs) have been proposed to provide uncoupled evaluations of robot translational performances along specific directions. Such indexes, which might be considered an extension of the traditional manipulability index, overcome the limitations of manipulability (but also of other popular indexes such as the minimum singular value and the condition number of the Jacobian matrix) which just provide global evaluations of the robot capabilities that are often of limited practical usefulness.

The experimental investigation carried out on a suitably instrumented Adept Quattro commercial robot has proved that the proposed DSI formulation can provide reliable predictions of the robot performances in making movements along specific directions. In particular, the comparison among the values of DSIs ( and ), manipulability , and the normalized maximum velocities experimentally recorded performing horizontal and vertical movements in a large number of points has demonstrated that DSIs can better foresee the best performance regions within the workspace and that the performance variation predictions made through the DSIs are in good agreement with the recorded velocity variations. Overall the predictions made through the DSIs are considerably more accurate than the predictions provided by manipulability indexes.

So far, the effectiveness of DSIs has been assessed on a single family of parallel manipulators (4-RUU), but the results achieved are thought to have a general relevance; DSI definition just recourses to vectors extracted by the inverse Jacobian matrix and includes no constraints accounting for a specific parallel robot architecture. The proposed indexes can therefore be applied to any parallel robot architecture, as long as an inverse Jacobian matrix can be computed.

In conclusion, DSIs, though being purely kinematic indexes, may provide useful hints in foreseeing the robot dynamic performances along relevant directions of motion within the workspace. At robot design, installation, or programming stages, such information could be usefully employed to optimize the robot geometrical features, the robot location within a workcell, the location of the target frames with respect to the robot, and also the robot end-effector paths. All these considerations make the DSIs a useful tool for robot designers, manufacturers, and programmers.

Acknowledgment

The authors acknowledge the financial support of the FSU foundation (“Fondazione Studi Universitari”—Vicenza) through a Ph. D. scholarship and significant equipment grants.