- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Author Guidelines ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Table of Contents
Advances in Mechanical Engineering
Volume 2010 (2010), Article ID 186203, 9 pages
Kinetostatic and Inertial Conditioning of the McGill Schönflies-Motion Generator
1Dipartimento Ingegneria Industriale e Meccanica, Università di Catania, Viale A. Doria 5, 95123 Catania, Italy
2Department of Mechanical Engineering & Centre for Intelligent Machines, McGill University, 817 Sherbrooke Street W., Montreal, QC, Canada H3A 2K6
Received 1 April 2009; Revised 31 August 2009; Accepted 30 October 2009
Academic Editor: Tian Huang
Copyright © 2010 Alessandro Cammarata et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This paper focuses on the optimization of the McGill Schönflies Motion Generator. Recent trends on optimum design of parallel robots led us to investigate the advantages and disadvantages derived from an optimization based on performance indices. Particularly, we optimize here two different indices: the kinematic conditioning and the inertial conditioning, pertaining to the condition number of the Jacobian matrix and to that of the generalized inertia matrix of the robot, respectively. The problem of finding the characteristic length for the robot is first investigated by means of a constrained optimization problem; then plots of the kinetostatic and the inertial conditioning indices are provided for a particular trajectory to be tracked by the moving platform of the SMG. Deep connections appear between the two indices, reflecting a correlation between kinematics and dynamics.
In recent years researchers and industry have been paying more attention to parallel-kinematics machines (PKMs) with reduced mobility [1–4]. The growing search for high performance in terms of accuracy, precision, maximum acceleration, and ease in control revealed some limits of PKMs with full mobility to accomplish certain types of tasks . The main characteristic of a PKM with reduced mobility is a reduced number of limbs with respect to a six-degrees-of-freedom (dof) PKM. This aspect yields a reduced weight of the global structure, which leads to an enhanced performance of the moving platform (MP). Besides, this feature implies a lower complexity of the architecture, which leads a lower risk of limb collisions, while increasing the workspace. Obviously, robots of this type are designed for specific applications and cannot perform the general tasks of their six-dof counterparts. Then again, industrial tasks usually call for specific types of motions; that is, not all dof of the MP are necessary. Some applications, as pick-and-place, drilling, riveting, assembling, painting, and so on, call for three to five dof. If a particular task required more dof than those provided by a single PKM, two or more PKMs might be assembled to cooperate. Unfortunately, the simpler architecture of these robot does not translate into a simpler analysis. The main issues lie in type synthesis, which excludes the use of the intuition at the conceptual design stage; a more systematic synthesis procedure of this new class of robots is needed .
Achieving high performance involves optimization in terms of kinematics and dynamics. Essentially, a performance index is a scalar function related to the kinematic or the dynamic performance of a robot. Kinematic-performance indices are mainly based on the robot Jacobian matrix. Various performance indices based on the Jacobian matrix have been proposed: we cite the conditioning of robotic manipulators [7–9], or the concept of manipulability, as proposed by Yoshikawa , is an attempt to measure the distance to singularity of the Jacobian matrix of a robot. One of the most frequently used performance indices, the kinetostatic conditioning index, is based on the condition number of the Jacobian matrix . We refer to [12–14] for a summary of the aforementioned kinematic-performance indices.
We shall devote our paper to only the kinetostatic conditioning index, the one studied most intensively.
Most of the dynamic performance indices have been derived through analogy from the kinetostatic indices. As a matter of fact, the dynamic conditioning index is obtained similar to the kinetostatic conditioning index by replacing the Jacobian matrix with the generalized inertia matrix (GIM). By the same token, Yoshikawa introduced the dynamic manipulability measure as the counterpart to his kinematic manipulability . Wiens et al.  proposed some indices for the measure of the nonlinear inertia forces while Khatib and Burdik  defined the isotropic acceleration as the greatest magnitude that the acceleration of the end-effector can reach in any direction starting from a given manipulator configuration. Ma and Angeles introduced the dynamic isotropy of an -dof serial manipulator and the “dynamic conditioning index’’ as a measure of the latter . Di Gregorio and Parenti-Castelli proposed the use of three homogeneous coordinates, chosen among the operational-space coordinates, for the definition of a three-dof-manipulator GIM to be used in the determination of the dynamically isotropic configurations meaningful for three-dof-manipulator design . Di Gregorio studied two-dof mechanisms and defined some indices that characterize the dynamic isotropy of these mechanisms .
In this paper the case of the McGill Schönflies-Motion Generator (SMG) is analyzed to apply the concepts of kinetostatic and inertial conditioning.
First, a brief description of the robot is provided. Then, the concept of characteristic length for positioning and orienting robots is recalled and adapted to the case of the McGill SMG. Further, a constrained optimization problem is set to find the posture of the robot with a Jacobian matrix and a GIM of minimum condition number. Finally, the same indices are analyzed while the robot moves along a standard trajectory, the test-cycle trajectory, to show that it would be possible to find an optimum location for the said trajectory where the indices attain minimum values.
2. Robot Description
Here a brief description of the McGill SMG is given, further details being available in . The McGill SMG is a parallel robot for SCARA-type motions: its moving platform can undergo three independent translations and one rotation about a vertical axis. A CAD model is shown in Figure 1, where one can distinguish the different parts composing the robot. Starting from the top of the fixed frame, each limb is composed of two motors, an epicyclic gear train (EGT), a right-angled gear box (RGB), and an RR kinematic chain. Then, the two chains, one for each limb, are coupled to a moving platform (MP). The EGTs and the RGBs are made of steel, while the two brackets, referred to as the distal and the proximal brackets, links of the RR chain, are fabricated of aluminum. Finally, the joints are slender plate elements fabricated of carbon-fiber to make these joints lighter.
Figure 2 shows the kinematic chain of the SMG. Hereafter, we will refer to the set of angles as the actuated-joint variables, while the remaining angles, , are the passive-joint variables. The DKP gives four different solutions for the actuated-joint variables when the generic configuration of the robot is defined. It is noteworthy that, during the optimization problem, we will refer only to a single manifold of solutions by imposing suitable initial values for the said angles. Furthermore, the plane of the joint of limb is normal to the unit vector , .
In Table 1 the notation used throughout the paper is described. All symbols with double subscript, comprising one Roman and one Arabic numeral, refer to the limb number, or , and to the corresponding item of the respective arm. In Table 1, and .
We will not dwell on the derivation of the kinematics and dynamics model of the SMG, described in detail in .
3. The Characteristic Length
Before starting with the conditioning-index definition, we recall the kinematics equation of a parallel robot: which relates the array of actuated joint rates with the twist array of the MP at the operation point, located on the moving platform, by means of two matrices, and , referred to as the direct-kinematics and the inverse-kinematics Jacobian matrices, respectively. For the case at hand, in light of the reduced mobility of the MP of the SMG, we have as derived in . It is noteworthy that and are Jacobians, (2) thus entailing six equations linear in the four independent components of either or . We resorted to this redundancy to add robustness to the kinematics model, so as to prevent algorithmic singularities—not related to the actual kinematics of the robot—which occur whenever the planes of the two limbs coincide. Here, solving (1) for the twist array , let us introduce the definition of the SMG Jacobian as where is the left Moore-Penrose generalized inverse of .
For robots accomplishing positioning and orienting tasks, like every SCARA, and the SMG is a parallel version thereof, the Jacobian bears different units, thus preventing the computation of the Jacobian condition number for optimization tasks. In order to remove the dimensional inhomogeneity, a characteristic length is introduced. The characteristic length is defined as that normalizing length that renders the condition number of the Jacobian matrix a minimum . As we will explain in detail in this section the condition number of the Jacobian matrix is bounded from below by unity but unbounded from above. Two situations can occur: robots whose Jacobian condition number can attain the minimum value of unity at one or more postures, termed isotropic robots, and robots whose Jacobian condition number attains a minimum larger than unity, at certain postures. Isotropy is not the rule but depends on the robot architecture.
It has been demonstrated that for spatial serial manipulators the characteristic length is the root mean square of the distances of the revolute axes to the operation point when the robot finds itself at a posture of minimum condition number . For parallel robots this geometric condition has not been found, although the original definition still stands. Finding a formula for the characteristic length is possible only for simple cases; more generally, the characteristic length is calculated via an optimization procedure.
The problem we solve here is formulated as follows: given a manipulator with a prescribed architecture, find its characteristic length as that which renders the Jacobian dimensionless and its condition number a minimum at an optimum posture. As the characteristic length is strictly connected to the concept of condition number of a matrix, that is, the Jacobian matrix, we briefly recall some pertinent concepts. It is noteworthy that the concept of condition number of a matrix is meaningful if and only if all the matrix entries bear the same physical units. The simplest definition of condition number is based on the matrix 2-norm, namely, , where and are the maximum and the minimum singular values of the matrix M. Apparently, we can order all the singular values of a matrix if and only if all its entries bear the same units. This definition, based on the 2-norm, is simple to state, but rather cumbersome to work with because is not an analytic function of the matrix entries everywhere, which brings about some computational difficulties. For this reason we use the condition number based on the Frobenius norm , namely, where, for SMGs, takes the form which is the weighted Frobenius norm of a matrix, as the Jacobian at hand has four nonzero singular values. Using (6), becomes We recall further that , like any form of condition number, is bounded from below but unbounded from above, that is,
Once again, a robot is termed isotropic when the condition number of its Jacobian reaches its minimum value of unity. It should be underlined that the condition number is configuration-dependent, and so, a manipulator can attain isotropic configurations only at certain points of its workspace if its design so permits.
For the McGill SMG we will set a constrained optimization problem to find its characteristic length L and its minimum condition number. Firstly, let be the design-variable vector whose components are the joint angles and the characteristic length L. We start by noticing that the angles , , have no effect on the Jacobian of (3); hence, only the remaining six angles will be taken into account. It should be mentioned that it is possible to chose a minimum set made of only actuated angles and the characteristic length L, but this choice, while reducing the number of design variables, makes the calculation more cumbersome, and so, we will not attempt it. The design-variable vector is then defined as Now, after introducing the vectors , , the optimum value of is found by solving the constrained optimization problem: subject to the constraints: where is the -component of , and is the radius of the MP. Notice that (11a)-(11b) implicitly take into account the dependency of the passive-joint angles from their actuated counterparts. The first constraint in (11a) must be set to preserve the parallel architecture, while the second in (11b) guarantees that the moving platform undergoes Schönflies displacements only.
It is noteworthy that the minimization problem is set for all configurations inside the workspace of the robot. As the moving platform has four dof, that is, three translations and one rotation, we might give a complete representation of the workspace only in a four-dimensional space. To overcome this issue we will recur to the constant orientation workspace representation, that is, the reachable volume attained by the operation point of the MP, its centroid in our case. In Figure 3 the said workspace is plotted for different values of the angle of rotation of the MP.
The results of the minimization problem are reported in Table 2 and Figure 4. As we can see, two postures, symmetric with respect to the YZ-plane, in which the robot Jacobian matrix attains minimum condition number , exist. The values of and are reported in Table 2.
Apparently, the condition number is very close to unity but is not unity; hence, the McGill SMG is not an isotropic robot.
4. Inertial Conditioning
Germane to the kinematic conditioning, the concept of dynamic isotropy relates to the GIM of a robotic mechanical system. It has been recognized that it is convenient to have an isotropic GIM, that is, of the form where is the identity matrix and a real, positive number, the unique eigenvalue of , of algebraic multiplicity . An isotropic GIM is convenient because this matrix must be inverted when solving for joint accelerations in the direct-dynamics problem. The inertia matrix being positive-definite, its eigenvalues are identical to its singular values.
Usually, the units of the GIM depend on the generalized coordinates chosen to describe a mechanical system. In robotics, for control purposes, it is often convenient to make these coincident with the actuated-joint variables. Therefore, the entries of the GIM of a robot can bear disparate units. The McGill SMG is actuated by two pairs of motors, all entries of thus bearing units of kg , and hence, the GIM is dimensionally homogeneous. We have chosen to consider in our analysis also the EGTs and RGBs because most of the inertia of the SMG is concentrated on them.
The global GIM is derived upon adding the three inertia matrices, , , and , of the three subsystems making up the SMG: EGTs, RGBs, and aluminum parts. Hence, the global GIM becomes All parameters used to obtain these matrices are recorded in , the corresponding expressions derived in the same reference.
We formulate now an optimization problem to find the posture of the SMG at which its GIM attains its minimum condition number . Now, is the objective function to minimize, as in (10), with respect to a new design-variable vector , defined as subject to the same constraints introduced in (11a)-(11b). It is noteworthy that the global GIM is a matrix whose entries bear all units of kg m. Figure 5 shows one of the two symmetrical minimum-condition-number postures, while Table 3 records the values of and for the said postures.
The GIM at the minimum-condition-number posture is close to isotropic, namely, Notice that matrix comes only two distinct blocks, which comes as no surprise, given the symmetries, geometric, and mechanical, of the robot.
5. Test-Cycle Trajectory
In this section we apply the concepts of kinetostatic and inertial conditioning to optimize the pose—position and orientation—of a test-cycle trajectory in the - plane. The test-cycle trajectory is a particular trajectory commonly adopted by the industry to represent the standard pick-and-place task performed by a SCARA robot. In our case, the test-cycle trajectory lies in the medium vertical plane of the robot. The SMG should pick an object, release it, and return to the initial MP pose in a given cycle of duration . Shown in Figure 6 is the trajectory, projected into the - plane, as generated using a 4-5-6-7 polynomial displacement program, which is reproduced here for quick reference : where is a nondimensional displacement function, with as a nondimensional time parameter, is time, and is the time elapsed between the ends of the trajectory. A characteristic of this polynomial is that it provides zero jerk at the start and the end of the motion, which means that the polynomial is -continuous . The fourth dimension of the MP trajectory involves a rotation of about a vertical axis, as the operation point of Figure 2 moves horizontally, going back to the original MP orientation on the return part of this motion. In Figure 7 the X and Z coordinates of point P and the angle of rotation of the moving plate are separately shown over the whole trajectory.
5.1. Kinetostatic Conditioning along the Test-Cycle Trajectory
The kinetostatic index is used here to compare different locations, on the vertical plane of symmetry of the robot, of the same test-cycle trajectory. Hence, the kinetostatic conditioning index has been evaluated on the test-cycle trajectory (in solid line) and on two other trajectories shifted along the -axis by 0.05 m and 0.10 m, but lying in the same XZ-plane, as shown in Figure 8. We can notice that the condition number worsens when the test-cycle trajectory is translated closer to the base.
5.2. Inertial Conditioning along the Test-Cycle Trajectory
Figure 9 shows how varies on the test-trajectory plane, increasing when the operation point of the moving platform approaches the base. Comparing Figures 8 and 9 we realize that the inertial conditioning index follows the same trend of its kinetostatic counterpart, but the values seem to be amplified and distorted in their scale and maximum peaks.
Finally, in Figure 10 we plot the condition number of on the -- space. A volume of size (0.2 m, 0.25 m, 2.4 rad) has been considered, while the coordinate has been held fixed to 0.3 m, the value at the XZ-plane of symmetry of the robot. Some slicing planes have been introduced to better visualize results. Following the said figure we can observe how worsens rapidly as the robot moves towards lower values of z. This means that it will be more convenient to move the test-cycle trajectory, onto the XZ-plane of symmetry of the robot, as far as possible from the base.
Summarizing results we can conclude that the kinetostatic and inertial conditioning indices seem to have similar trends, the latter being a stretched and scaled version of the former. Therefore, one might consider only one of the two optimization problems to observe the robot behaviour on a quality level. Sometimes, as in the case of the McGill SMG, it might be easier to work with the inertial conditioning index rather than with its kinetostatic counterpart, for the latter does not require a characteristic length as a further design-variable. On the contrary, the kinetostatic index is usually simpler to calculate than the inertial index, for the latter needs the complete dynamics model.
The kinematic and dynamic conditioning of the McGill SMG were studied. Firstly, a constrained optimization problem was set to find the robot characteristic length and the minimum condition number of the robot Jacobian. Then, a similar constrained optimization problem was set to find the optimum robot posture of minimum condition number of the generalized inertia matrix. The kinetostatic and inertial conditioning indices were evaluated on a particular trajectory, the test-cycle, to be tracked by the MP. The evaluation has revealed similar behaviors between kinetostatic and dynamic indices. Particularly, the inertial effects seem to amplify and distort the values of the kinetostatic index, thereby modifying both the scale and the maximum peaks. Finally, the inertial conditioning index was plotted over a constant- volume to better understand the effects of a displacement of the test-cycle trajectory on the XZ-plane of symmetry of the robot.
- M. Callegari, M. Palpacelli, and M. Scarponi, “Kinematics of the 3-CPU parallel manipulator assembled for motions of pure translation,” in Proceedings of the IEEE International Conference on Robotics and Automation, pp. 4020–4025, Barcelona, Spain, April 2005.
- X.-J. Liu, J. Wang, F. Gao, and L.-P. Wang, “On the analysis of a new spatial three-degrees-of-freedom parallel manipulator,” IEEE Transactions on Robotics and Automation, vol. 17, no. 6, pp. 959–968, 2001.
- W. Li, F. Gao, and J. Zhang, “R-CUBE, a decoupled parallel manipulator only with revolute joints,” Mechanism and Machine Theory, vol. 40, no. 4, pp. 467–473, 2005.
- M. Carricato and V. Parenti-Castelli, “On the topological and geometrical synthesis and classification of translational parallel mechanisms,” in Proceeding of the 11th World Congress in Mechanism and Machine Science, pp. 1624–1628, Tianjin, China, April 2004.
- X.-J. Liu, J. Wang, and G. Pritschow, “A new family of spatial 3-DoF fully-parallel manipulators with high rotational capability,” Mechanism and Machine Theory, vol. 40, no. 4, pp. 475–494, 2005.
- X. Kong and C. Gosselin, Type Synthesis of Parallel Mechanisms, vol. 33 of Springer Tracs in Advanced Robotics, Springer, Berlin, Germany, 2007.
- D. C. Yang and Z. C. Lai, “On the conditioning of robotic manipulators-service angle,” ASME Journal of Mechanisms, Transmissions, and Automation in Design, vol. 107, pp. 262–270, 1985.
- X.-J. Liu, Z.-L. Jin, and F. Gao, “Optimum design of 3-DOF spherical parallel manipulators with respect to the conditioning and stiffness indices,” Mechanism and Machine Theory, vol. 35, no. 9, pp. 1257–1267, 2000.
- T. Huang, Z. Li, M. Li, D. G. Chetwynd, and C. M. Gosselin, “Conceptual design and dimensional synthesis of a novel 2-DOF translational parallel robot for pick-and-place operations,” ASME Journal of Mechanical Design, vol. 126, no. 3, pp. 449–455, 2004.
- T. Yoshikawa, “Manipulability of robotic mechanisms,” International Journal of Robotics Research, vol. 4, no. 2, pp. 3–9, 1985.
- J. K. Salisbury and J. J. Craig, “Articulated hands: force and kinematic issues,” International Journal of Robotics Research, vol. 1, no. 1, pp. 4–17, 1982.
- J. Angeles, Fundamentals of Robotic Mechanical Systems, Springer, New York, NY, USA, 3rd edition, 2007.
- J. P. Merlet, “Jacobian, manipulability, condition number, and accuracy of parallel robots,” ASME Journal of Mechanical Design, vol. 128, no. 1, pp. 199–206, 2006.
- G. Gogu, Structural Synthesis of Parallel Robots—Part 1: Methodology, Springer, Berlin, Germany, 2007.
- T. Yoshikawa, “Dynamic manipulabilty of robot manipulators,” in Proceedings of the International Conference on Robotics and Automation, pp. 1033–1038, St. Louis, Mo, USA, 1985.
- G. J. Wiens, R. A. Scott, and M. Y. Zarrugh, “The role of inertia sensitivity in the evaluation of manipulator performance,” ASME Journal of Dynamic Systems, Measurement and Control, vol. 111, no. 2, pp. 194–199, 1989.
- O. Khatib and J. Burdik, “Optimization of dynamics in manipulator design: the operational space formulation,” The International Journal of Robotics and Automation, vol. 2, pp. 90–98, 1987.
- O. Ma and J. Angeles, “The concept of dynamic isotropy and its applications to inverse kinematics and trajectory planning,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '90), pp. 481–486, Cincinnati, Ohio, USA, 1990.
- R. Di Gregorio and V. Parenti-Castelli, “Dynamic performance indices for 3-DOF parallel manipulators,” in Advances in Robot Kinematics: Theory and Applications, J. Lenarcic and F. Thomas, Eds., pp. 11–20, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2002.
- R. Di Gregorio, “Dynamic model and performances of 2-dof mechanisms,” in Proceedings of the ASME Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Salt Lake City, Utah, USA, September-October 2004, paper no. DETC2004-57129.
- A. Cammarata, J. Angeles, and R. Sinatra, “The dynamics of parallel Schönflies motion generators: the case of a two-limb system,” Proceedings of the Institution of Mechanical Engineers I, vol. 223, no. 1, pp. 29–52, 2009.
- G. H. Golub and C. F. Van Loan, Matrix Computations, The Johns Hopkins University Press, Baltimore, Md, USA, 1989.
- W. A. Khan and J. Angeles, “The kinetostatic optimization of robotic manipulators: the inverse and the direct problems,” ASME Journal of Mechanical Design, vol. 128, no. 1, pp. 168–178, 2006.