- 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 2012 (2012), Article ID 171682, 10 pages

http://dx.doi.org/10.1155/2012/171682

## Choosing Actuation Scheme for Optimal Performance of 3-DOF Planar Parallel Manipulators

^{1}Department of Mechanical Engineering, G.V.P. College of Engineering, Visakhapatnam 530048, India^{2}Department of Mechanical Engineering, JNTU College of Engineering, Vizianagaram 535003, India^{3}Department of Mechanical Engineering, AU College of Engineering, Visakhapatnam 530003, India

Received 17 October 2011; Revised 6 March 2012; Accepted 30 March 2012

Academic Editor: Adib Becker

Copyright © 2012 S. Ramana Babu 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.

#### Abstract

This paper describes how the actuation scheme affects the size and shape of the workspace of 3RPR and 3PRR, which are the two popular planar parallel manipulators. This study helps in achieving the required singularity-free zones for a specific task by adopting a suitable actuation scheme without resorting to changing the architecture of the manipulator. In addition, the performance index GCI for each actuation scheme has been evaluated and compared with each other. Finally, the GCI for each actuation scheme has been optimized subjected to the geometric constraints for the 3PRR planar parallel manipulator.

#### 1. Introduction

Parallel manipulators are robots that consist of separate serial chains that connect the fixed link to the moving end-effectors. Even though parallel manipulators possess several advantages like, high stiffness, low inertia, and large payload capacity, they are suffered from limited workspace, complicated singularities, and design difficulties. Many researchers have analyzed kinematics, dynamics, workspace, and control of parallel manipulators. Mohammadi Daniali et al. [1] presented an in-depth study of velocity relationships and singular conditions for general planar parallel robots. Gosselin [2] presents computational algorithms for kinematics and dynamics of planar and spatial parallel robots. Merelt et al. [3] solved the forward pose kinematics for a broad class of planar parallel manipulators. Y. Jin et al. [4] described a 3-limb selectively actuated parallel manipulator (SA-PM) in which the end-effecter exhibits either a 3-DOF spherical, 3-DOF translational, or a complete 6-DOF spatial motion depending on the types of the actuation chosen for the actuators. They also presented singularity analysis of the SA-PM based on geometry for all actuation schemes to facilitate the kinematic design of the manipulator. Arakelian et al. [5] proposed the legs of variable structure for the increase of reachable workspace of a spatial manipulator. Su et al. [6] studied trajectory planning and actuation schemes, which are crucial for their control. Ebrahimi et al. [7] presented a model of an actuation scheme and its effects on the singularities of parallel manipulators for a given path in the workspace. Actuation schemes and actuator relocation for serial planar manipulator have also been studied in detail by Matone and Roth [8] and Ramji et al. [9] respectively. Li and Xu [10] designed a new XYZ novel micromanipulator, which was called a 3-PRC manipulator in order to generate a cuboids’ shape workspace by using the (PZT) piezoelectric actuator to drive the prismatic active joint due to its major advantages of high stiffness and compact size; however, stroke of the adopted PZT cannot meet the application requirements so in order to enlarge the workspace the stroke of the PZT is amplified by using a lever amplification mechanism and making the leg as an equivalent PRPR linkage. Xu and Li [11] proposed a new PRC translational parallel manipulator with orthogonally arranged fixed actuators for producing a singularity-free workspace through an extensive reduction in cost in building and controlling such a manipulator. Li and Xu [12] analyzed the mobility of the 3-PRC translational parallel manipulator via screw theory and also illustrated the effect of architectural parameters on the shape and size of the workspace; further more dexterity characteristics are investigated in the local and global sense for real machine design.

Actuation scheme depends on the Jacobian matrix that transforms the active input joint rates into end-effecter task space velocities. Earlier researchers on PKMs studied the kinematic analysis and performance of the manipulator for only one set of active joint rates (consists of only one joint from one leg) preferably the base joints. Where the kinematic and workspace analysis is not only confined to one set of active joints and also tried to investigate the changes in reachable workspace and performance of manipulators for the other set of joints, which are chosen as active joints. We have to enumerate the optimal geometric parameters like link lengths, joint variables, and orientation of the prismatic joints, which greatly influence the workspace and performance of the manipulator for each possible actuation scheme. Though actuation schemes have been proposed by various researchers, these have not been applied to any specific case except by Rakotomanga et al. [13]. This paper is an attempt to apply the concept of actuation schemes to two popular planar manipulators, namely, 3RPR and 3PRR.

The organization of the paper is as follows: Method explains the basic theory of architecture, kinematic modeling, and actuation schemes of 3RPR manipulator first, followed by a similar treatment for the 3PRR manipulator. Results describes through graphs and tables, comparison of the performance of a manipulator for all possible actuation schemes. Also included is the optimal geometric design of 3PRR planar parallel manipulator. These results are followed by Discussion.

#### 2. Method

##### 2.1. Mechanism Architecture of 3RPR Planar Manipulator

The manipulator considered here is symmetric and consists of three identical legs connecting the fixed base to the end-effecter as shown in Figure 1. The three revolute joints are located on the fixed base triangle at and another three revolute joints are located on the moving triangle at where . The prismatic joint variables correspond to the leg lengths . The moving frame is located at the centroid of the end-effecter, which is an equilateral triangle as shown in Figure 2. The fixed frame is chosen at the origin of one of the legs on the fixed base, which is also an equivalent triangle.

For the 3RPR manipulator, the workspace limitations are due to the limitations of prismatic actuators, and the maximum and minimum lengths of the prismatic actuators of the th chain are denoted by , , respectively, as the prismatic joint variables. The maximum possible geometric workspace of 3RPR manipulator is the intersection of three annular regions centered in whose radii are and as shown in Figure 3. This annular region will exist if and only if . The determination of the maximal workspace has been addressed by Kumar [14] and Pennock and Kassner [15] who pointed out that the boundary of this workspace is composed of circular arcs and portions of sextic curves.

###### 2.1.1. Kinematic Modeling of 3RPR Parallel Manipulator

The solutions are found using the distance norm constraint In the forward pose kinematics for the given prismatic joint lengths , the Cartesian pose coordinate can be calculated. Equation (1) can be expressed in analytical form as where is the orientation of the moving platform with reference -axis, is angle between vector and -axis of the moving platform, is length of vector , and is distance of the vertex of the fixed platform from the origin of the reference frame. The revolute joint angle, which is attached to the fixed platform, is expressed as It is expressed in analytical form as

###### 2.1.2. Actuation Schemes and Velocity Analysis

The major disadvantage of parallel kinematic chains is their limited workspace and in order to overcome this, recently some researchers have investigated hybrid serial-parallel manipulators that present a compromise between the high rigidity of parallel manipulators and the extended workspace of serial manipulators. However, the usefulness of this manipulator is also limited because the compliance of the manipulator is more, and so it leads to inaccuracy in controlling the manipulator. There is another trend to eliminate the singularities and increase the workspace by the redundant actuation with the use of more than one actuator in a leg. However, this redundant actuation is expensive due to the additional actuators and the complication in controlling the manipulator.

The actuation scheme is the set of active input joint rates of the manipulator in which only one joint from each leg is to be selected. Each actuation scheme has its own singularity free zones of the end-effecter. Since each scheme implies a different set of actuator motions, this leads to different movements of the end-effecter. A proper actuation scheme not only avoids singularities of the workspace but even enhances the stiffness and accuracy of the manipulator.

The 3RPR manipulator can be actuated either by revolute joint or by prismatic joint in each leg so to attain the required end-effecter output velocities . In order to change the end-effecter position, we have to actuate either the base revolute joint or the linear actuator , that is, only one input from each leg so as to avoid redundancy since it has its own drawbacks. Only six joints are available for input out of the nine, and among these, we have to choose three joints, one from each leg. Hence there exist eight possible combinations and each such set is called an actuation scheme. These actuation schemes are presented in Table 1. The underlined letters show the input pairs, “__R__” for revolute with input angle and “__P__” for prismatic with input displacement .

###### 2.1.3. Jacobian Matrix

The velocity vector of a point is given in two different directions of a loop closure. Each loop closure consists of a fixed base, moving platform, and all links of a leg. The unactuated joint rates in each leg are eliminated by taking dot product of the velocity vector loop equation with a vector that is normal to all vectors of unactuated joint rates. Finally the resulting equations are assembled into a Jacobian matrix. In order to obtain the Jacobian matrix for the actuating joint differentiate (2) with respect to time and is expressed as where
The Jacobian for the actuating joint is obtained by differentiating (4) with respect to time and is expressed as
Using the above equations (5) and (7), the Jacobian matrix for eight actuation schemes is expressed as
For actuation scheme 1 (R__P__R1-R__P__R2-R__P__R3), the joint variables are passive joints, while the active joints are the prismatic joints . The Jacobian matrix mapping joint rates into Cartesian rates is also expressed in two separate Jacobian matrices as
where

Which is analogous to the , where and are inverse and direct Jacobian matrices, respectively, the overall Jacobian matrix . For actuation scheme 2 (__R__PR1-__R__PR2-__R__PR3), the active joint rates into the Cartesian rates . The velocity relation is expressed in two separate matrix forms as

##### 2.2. 3PRR Manipulator Architecture

A 3-DOF planar parallel manipulator with three identical legs, which has been studied extensively by many researchers including Zhang et al. [16], is shown in Figure 4. Each kinematic chain is of PRR type and consists of one prismatic joint and two revolute joints. The manipulator is intended to position and orient the equilateral triangle-shaped platform . The geometric centre of platform denoted by is the operation point of the manipulator, with its position and orientation being . Three intermediate links for connect the sliders and the moving platform.

The base of the manipulator, also an equilateral triangle with vertices , , and , is fixed to the ground. The joint coordinate associated with the th actuator that actuates the prismatic joint is for . The variable is the angle at between the -axis of the fixed frame and th linear guide . The revolute variable is the angle made by the th intermediate link with the -axis of the fixed frame. The vector loop closure equation for the th leg is written as Writing the above equation in component form, Using (11) and (12), we can establish a quadratic equation in where It is clear that there are two possible solutions for each individual link, and the manipulator can take on a maximum of eight configurations for a given set of coordinates of the moving platform.

###### 2.2.1. Velocity Analysis of Actuation Schemes

The 3RPR manipulator can be actuated either by revolute joint or by prismatic joint in each leg of the manipulator in order to attain the required end-effecter output velocities. Thus there are only three input joint rates out of six . There exist eight possible combinations; each combinational set of actuators is called an actuation scheme and these are shown in Table 2. The underlined letters show the input pairs, “__R__” for revolute with input angle and “__P__” for prismatic with input displacement . The actuation scheme 1 is represented as __P__RR1-__P__RR2-__P__RR3, where the prismatic joint variables are active joints. To obtain the velocity equations for actuation scheme 1, we differentiate (11) with respect to time so that
where , is the positional vector from to , is the positional vector from to , and is unit vector orthogonal to both and . The left-hand side of (9) can be interpreted as summation of the velocity of the platform at its mass centre and a velocity term due to the rotation of the platform. The right-hand side of (9) is the summation of the velocity of the slider and the velocity term due to rotation of the link. To remove the velocity term related to the rotational velocity of the link, in (17), we apply vector dot product on both sides of (17) with
The above equation can be expressed as
Equation (19) presents the velocity relationship between the sliders and the platform with the Jacobian matrix . Using (19) for , we can write Jacobian matrix for actuation scheme 1 (__P__RP1-__P__RR2-__P__RR3).

Performing vector cross product to (17) by , From (19) and (20), the angular velocity of the link is given as

Using (21) we have to get Jacobian matrix for the actuation scheme 2 represented as (P__R__R1-P__R__R2-P__R__R3). Here is the passive joint variable and is the active revolute joint variable for . The Jacobian matrix of other actuation schemes is derived using (19) and (21).

###### 2.2.2. Dexterity Indices

Yoshikawa [17] proposed the concept of kinematics manipulability as where is the Jacobin matrix and depends on the instantaneous configuration of the manipulator defined by a joint vector . This index is treated as a quantitative measure of the easiness to arbitrarily change the position and orientation at a certain posture. The significance of manipulability is the measure of the fast recovery ability from the escapable singular point for the redundant manipulator. Manipulability has two problems: scale and order deficiencies. Those problems prevent the fair comparison among the manipulators with different dimension and make it impossible to derive the physical meaning of the manipulability. Trying to eliminate the order dependency, Kim and Khosla [18] gave another way to calculate manipulability where is the number of degrees of freedom of the manipulator.

Another usually used index is the condition number of the Jacobian matrix that is expressed as Condition number signifies the error amplification factor between the joint rate errors in to task space errors. The condition number depends on the manipulator configuration, link lengths. It varies from one at isotropic configurations to infinity at singular configurations so this is also defined as the measure of degree of ill-conditioning of the Jacobian matrix. Both indices, condition number and manipulability, are local measures of performance at a particular pose of the end-effecter within the workspace region; using these indices, we cannot compare the performance of any two different manipulators for their usefulness. In order to evaluate the global behavior of a manipulator over the workspace, a global index is proposed by Gosselin and Angeles [19] as where is the workspace and is the condition number. In other sense, is average value of over workspace region; it represents the uniformity of dexterity over the entire workspace so this index makes sense for the optimal design of manipulator for which the average value performance is an important design factor.

#### 3. Results and Discussion

##### 3.1. Performance of Actuation Schemes for 3RPR Manipulator

Symmetric architectures are considered here since they are widely used in practice. The values chosen for the sides of the fixed link are m, while m are the sides of the moving platform. The leg length range is chosen as m.

The reachable workspace for all eight actuation schemes using the Monte Carlo algorithm is plotted and is shown in Figure 5. The Global Condition Index for each actuation scheme within the maximum possible reachable workspace is computed, and the results are shown in Table 3. Actuation scheme 1 (R__P__R1-R__P__R2-R__P__R3) has the highest of 0.4351 as well as a reasonably large dexterous workspace as indicated by the Monte Carlo points (29328). Actuation scheme 2 (__R__PR1-__R__PR2-__R__PR3), on the other hand, showed very poor performance with a GCI of 0.0074 and also lowest dexterous workspace (27778). The potential of each actuation scheme is quantified by these two numbers. The shape and magnitude of the reachable workspace is varied from one actuation scheme to other actuation scheme. The reachable workspace for actuation scheme 1 (R__P__R1-R__P__R2-R__P__R3) is registered for almost all orientations of the end-effecter as shown in Figure 5(a). The reachable workspace for actuation scheme2 and for actuation scheme3 is shown in Figures 5(b) and 5(c), respectively, where the tool platform is unable to move in certain orientations . For actuation scheme 4 (__P__RR1-__P__RR2-P__R__R3) a continuously increasing reachable workspace for wide range of orientations of moving platform is observed.

###### 3.1.1. Effect of Geometric Parameters on GCI

The performance values shown above pertain to one set of dimensions, and obviously these need not necessarily be the best. Hence various dimensions have been tried out keeping in mind that generally linear actuators are available in standard sizes. Of all possible combinations, the optimum ones based on GCI have been determined and the results have been plotted as shown in Figure 6. The variables considered are maximum length of the linear actuator and the size of the end-effecter for the given m and the base platform m. The peak values of GCI for each actuation scheme are shown in Table 4.

It can be seen that the global maximum of all GCIs occurs for actuation scheme 1 (R__P__R1-R__P__R2-R__P__R3) with a value of 0.4487 at the joint variable m, m. Even though the is greatly influenced by the geometric parameters in all the actuation schemes as shown in Figures 6(a) to 6(h), this value is changed from one actuation scheme to other actuation scheme for a given set of geometric parameters.

##### 3.2. Workspace and Performance Evaluation of 3PRR

Precise calculation of workspace and its boundaries is important because they influence the dimensional design, manipulator’s positioning in the work environment, and its dexterity within the reachable workspace. The workspace is limited by the boundary obtained through solving inverse kinematics, reachable extent of drives and joints, and the occurrence of singularities. A numerical example of the force-unconstrained poses of the 3PRR planar manipulator is presented. Both fixed and moving platforms are equilateral triangles, whose sides are m for the fixed and m for the moving platform. The stroke of the prismatic joint is m for , the intermediate leg lengths m for . For the workspace calculation, Rastegar and Fardanesh [20] introduced the Monte Carlo random sampling numerical method to generate the workspace and the boundary surfaces. The Global Condition Index and the number of Monte Carlo points are shown in Table 5 for the eight actuation schemes, and the workspace plots are shown in Figure 7. The size and shape of the workspace of the manipulator varies with the actuation. Among all the eight actuation schemes, the actuation scheme 2 (P__R__R1-P__R__R2-P__R__R3) showed the best performance with a GCI of 0.3241 and the number of Monte Carlo workspace points 218622. The GCI itself represents the quality of the conditioned workspace, and the number of Monte Carlo points indicates the size of the reachable workspace. Actuation scheme 7 (P__R__R1-__P__RR2-P__R__R3) also shows good performance with GCI of 0.1699 along with larger workspace size next to the actuation scheme 2 among all eight actuation schemes for 3PRR manipulator. The reachable workspace for actuation scheme 2 (P__R__R1-P__R__R2-P__R__R3) is uniform and in regular shape for all the orientations of the moving platform as shown in Figure 7(b). The workspace for actuation scheme 2 (P__R__R1-P__R__R2-P__R__R3) as shown in Figure 7(b) is registered for only one orientation of the tool platform; the moving platform is unable to move in other orientations within the maximal reachable workspace due to the presence of complex singularities. Actuation scheme 4 (__P__RR1-__P__RR2-P__R__R3) and actuation scheme 8 (__P__RR1-P__R__R2-__P__RR3) have got very limited reachable workspaces, and also the moving platform failed to move many orientations as shown in Figures 7(d) and 7(h), respectively.

###### 3.2.1. Optimization of Geometric Structure for 3PRR Manipulator

Sequential quadratic programming (SQP) methods attempt to solve a nonlinear program directly rather than convert it to a sequence of unconstrained minimization problems. The basic idea is analogous to Newton’s method for unconstrained minimization. At each step, a local model of the optimization problem is constructed and solved, yielding a step toward the solution of the original problem. A quadratic approximation to this function is now constructed and so along with linear constraints forms a quadratic programming problem that is, the minimization of a function quadratic in the variables, subject to linear and nonlinear constraints. The solution of the original optimization problem, say , is now obtained from an initial estimate and solving a sequence of updated quadratic programs. Sequential quadratic programming was used by Attolico et al. [21] and Xu and Wang [22]. It is a successful approach to solving constrained optimization problems. The vector of design parameters is .

The objective of this optimization algorithm is to determine optimal dimensions and optimal performance index for the eight actuations for the three legged parallel manipulators for a specified task. The optimization algorithm solves a nonlinear constrained optimization problem defined as follows:

For the given , , , , , , , , m; m. The optimized results are shown in Table 6.

Actuation scheme 2 (P__R__R1-P__R__R2-P__R__R3) possesses maximum reachable workspace size with maximum GCI of 0.5489, the optimal design vector corresponding to this actuation mode is .

#### 4. Conclusions

Kinematic analysis of 3RPR and 3PRR parallel manipulators with symmetric properties is presented in this paper. In this problem, actuation scheme is identified as one of the design variables of the mechanism. It is shown that actuation scheme affects the performance of the manipulator and size and shape of the workspace significantly. The dimensional synthesis for optimization of Global Condition Index for each actuation scheme is presented.

Thus this paper offers a methodology to choose the ideal actuation scheme so as to attain the best performance of the manipulator with desired workspace; this is illustrated on two planar 3PRR and 3RPR manipulators as an example. The same technique can be applied to other planar PKMs.

#### References

- H. R. Mohammadi Daniali, P. J. Zsombor-Murray, and J. Angeles, “Singularity analysis of planar parallel manipulators,”
*Mechanism and Machine Theory*, vol. 30, no. 5, pp. 665–678, 1995. View at Scopus - C. M. Gosselin, “Parallel computational algorithms for the kinematics and dynamics of planar and spatial parallel manipulators,”
*Journal of Dynamic Systems, Measurement and Control, Transactions of the ASME*, vol. 118, no. 1, pp. 22–28, 1996. View at Scopus - J. P. Merlet, S. Lemieux, and C. M. Gosselin, “New architecture of planar three-degree-of-freedom parallel manipulator,” in
*Proceedings of the 13th IEEE International Conference on Robotics and Automation*, pp. 3738–3743, April 1996. View at Scopus - Y. Jin, I. M. Chen, and G. Yang, “Structure synthesis and singularity analysis of a parallel manipulator based on selective actuation,” in
*Proceedings of the IEEE International Conference on Robotics and Automation*, pp. 4533–4538, New Orleans, La, USA, May 2004. View at Scopus - V. Arakelian, S. Briot, and V. Glazunov, “Improvement of functional performance of spatial parallel manipulators using mechanisms of variable structure,” in
*Proceedings of the 12th IFTOMM World Congress*, Besancon, France, 2007. - H. J. Su, P. Dietmaier, and J. M. McCarthy, “Trajectory planning for constrained parallel manipulators,”
*Journal of Mechanical Design, Transactions of the ASME*, vol. 125, no. 4, pp. 709–716, 2003. View at Publisher · View at Google Scholar · View at Scopus - I. Ebrahimi, J. A. Carretero, and R. Boudreau, “Actuation scheme for a 6-dof kinematically redundant planar parallel manipulator,” in
*Proceedings of the 12th IFTOMM World Congress*, Besancon, France, 2007. - R. Matone and B. Roth, “The effects of actuation schemes on the kinematic performance of manipulators,”
*Journal of Mechanical Design, Transactions of the ASME*, vol. 119, no. 2, pp. 212–217, 1997. View at Scopus - K. Ramji, V. Ramachandra Raju, and S. Ramanababu, “Effect of actuator relocation on performance of industrial manipulators,”
*International Journal of Mechanics of Solids*, vol. 5, pp. 71–82, 2009. - Y. Li and Q. Xu, “Design and optimization of an XYZ parallel micromanipulator with flexure hinges,”
*Journal of Intelligent and Robotic Systems*, vol. 55, no. 4-5, pp. 377–402, 2009. View at Publisher · View at Google Scholar · View at Scopus - Q. Xu and Y. Li, “Design and analysis of a new singularity-free three-prismatic-revolute-cylindrical translational parallel manipulator,”
*Proceedings of the Institution of Mechanical Engineers Part C*, vol. 221, no. 5, pp. 565–577, 2007. View at Publisher · View at Google Scholar · View at Scopus - Y. Li and Q. Xu, “Kinematic analysis and design of a new 3-DOF translational parallel manipulator,”
*Journal of Mechanical Design, Transactions of the ASME*, vol. 128, no. 4, pp. 729–737, 2006. View at Publisher · View at Google Scholar · View at Scopus - N. Rakotomanga, D. Chablat, and S. Caro, “Kinetostatic performance of planar parallel mechanism with variable actuation,”
*Advances in Robot kinematics*, pp. 1–11, 2008. - V. Kumar, “Characterization of workspaces of parallel manipulators,” in
*Proceedings of the ASME 21th Biennial Mechanisms Conference*, pp. 321–329, Chicago, Ill, USA, 1990. - G. R. Pennock and D. J. Kassner, “Workspace of a general geometry planar three-degree-of-freedom platform-type manipulator,”
*Journal of Mechanical Design, Transactions Of the ASME*, vol. 115, no. 2, pp. 269–276, 1993. View at Scopus - X. Zhang, J. K. Mills, and W. L. Cleghorn, “Dynamic modeling and experimental validation of a 3-PRR parallel manipulator with flexible intermediate links,”
*Journal of Intelligent and Robotic Systems*, vol. 50, no. 4, pp. 323–340, 2007. View at Publisher · View at Google Scholar · View at Scopus - T. Yoshikawa, “Manipulability of robotic mechanisms,”
*International Journal of Robotics Research*, vol. 4, no. 2, pp. 439–446, 1984. - J. O. Kim and P. Khosla, “Dexterity measures for design and control of manipulators,” in
*Proceedings of the IEEE/RSJ International Workshop on Intelligent Robot and Systems (IRSO '91)*, pp. 758–763, November 1991. View at Scopus - C. Gosselin and J. Angeles, “Global performance index for the kinematic optimization of robotic manipulators,”
*Journal of Mechanisms, Transmissions, and Automation in Design*, vol. 113, no. 3, pp. 220–226, 1991. View at Scopus - J. Rastegar and B. Fardanesh, “Manipulation workspace analysis using the Monte Carlo Method,”
*Mechanism and Machine Theory*, vol. 25, no. 2, pp. 233–239, 1990. View at Scopus - M. Attolico, P. Masarati, and P. Mantegazza, “Trajectory optimization and real-time solution for robotics applications,” in
*Proceedings of the ECCOMAS Thematic Conference: Multibody Dynamics*, Madrid, Spain, 2005. - J. X. Xu and W. Wang, “Two optimization algorithm for solving robotics inverse kinematics with redundancy,” in
*Proceedings of the IEEE International Conference on Control and Automation (ICCA '07)*, pp. 3021–3028, June 2007. View at Publisher · View at Google Scholar · View at Scopus