InTech uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Robotics » Robotics and Automation » "Parallel Manipulators, towards New Applications", book edited by Huapeng Wu, ISBN 978-3-902613-40-0, Published: April 1, 2008 under CC BY-NC-SA 3.0 license. © The Author(s).

Chapter 19

A Novel 4-DOF Parallel Manipulator H4

By Jinbo Wu and Zhouping Yin
DOI: 10.5772/5441

Article top

Overview

The Delta structure proposed by Prof. Clavel and one of its industrial version, the CE33 (courtesy of SIG Pack Systems)
Figure 1. The Delta structure proposed by Prof. Clavel and one of its industrial version, the CE33 (courtesy of SIG Pack Systems)
Joint-and-loop graph of the Delta robot and one of its equivalent structure
Figure 2. Joint-and-loop graph of the Delta robot and one of its equivalent structure
The structure of a Delta robot with no internal DOF.
Figure 3. The structure of a Delta robot with no internal DOF.
Architecture scheme of H4.
Figure 4. Architecture scheme of H4.
A simple scheme of H4 structure
Figure 5. A simple scheme of H4 structure
Symmetrical H4 with actuated prismatic joints and rotational joints
Figure 6. Symmetrical H4 with actuated prismatic joints and rotational joints
Asymmetrical design of H4.
Figure 7. Asymmetrical design of H4.
I4 prototype (courtesy of LIRMM)
Figure 8. I4 prototype (courtesy of LIRMM)
Part4 prototype (courtesy of LIRMM)
Figure 9. Part4 prototype (courtesy of LIRMM)
H4 manipulator with four lines drives and (S-S)2 chains
Figure 10. H4 manipulator with four lines drives and (S-S)2 chains
Design parameters of H4
Figure 11. Design parameters of H4
Geometrical conditions that characterize the varieties of dimension 1, 2 and 3 (Merlet, 2006)
Figure 12. Geometrical conditions that characterize the varieties of dimension 1, 2 and 3 (Merlet, 2006)
Geometrical conditions that characterize the varieties of dimension 4 and 5 (Merlet, 2006)
Figure 13. Geometrical conditions that characterize the varieties of dimension 4 and 5 (Merlet, 2006)
Forces and moments transmitted to the nacelle
Figure 14. Forces and moments transmitted to the nacelle
Forces and moments transmitted to C1C2
Figure 15. Forces and moments transmitted to C1C2
Parallelogram static analysis
Figure 16. Parallelogram static analysis
Serial singularity of H4 manipulator
Figure 17. Serial singularity of H4 manipulator
H4 robot with linear complexes (a): 3D view (b): top view
Figure 18. H4 robot with linear complexes (a): 3D view (b): top view
H4 robot with linear complexes (a): 3D view (b): top view
Figure 19. H4 robot with linear complexes (a): 3D view (b): top view
The constraint singular surface
Figure 20. The constraint singular surface
Minimum σ in the vicinity of the architecture singular configuration
Figure 21. Minimum σ in the vicinity of the architecture singular configuration
Pitch of linear complex with minimum σ in the vicinity of the architecture singular configuration
Figure 22. Pitch of linear complex with minimum σ in the vicinity of the architecture singular configuration

A Novel 4-DOF Parallel Manipulator H4

Jinbo Wu1 and Zhouping Yin2

1. Introduction

Parallel manipulators have the advantages of high stiffness and low inertia compared to serial mechanisms. Based on the Steward-Gough platform architecture, a lot of 6-DOF mechanical devices have been proposed. The 6-DOF parallel manipulators suffer from a small workspace, complex mechanical design, and difficult motion generation and control due to their complex kinematic analysis. To overcome these shortcomings, the limited-DOF manipulator, which has fewer than 6 DOFs, can be found in many production lines. It is clear today that most attention has been paid to 3-DOF family among the limited-DOF parallel manipulators (Carretero, 2000). However, in many industrial situations, there is a need for equipment providing more than 3-DOFs. For example, for most pick-and-place applications in semiconductor manufacturing, at least 4 DOFs are required (3 translation to move the carried die from one point to the other, 1 rotation to adjust the orientation in its final location). A new family of 4-dof parallel manipulators called H4 that could be useful for high-speed pick-and-place applications is proposed by Pierrot and Company (Pierrot, 1999). The H4 manipulator offers 3 DOFs in translation and 1 DOF in rotation about a given axis. The H4 manipulator is useful for high-speed handling in robotics and milling in machine-tool industry since it is a fully-parallel mechanism with no passive chain and able to provide high performance in terms of speed and acceleration.

This chapter discusses the kinematic analysis of the H4 manipulator. In section 2, synthesis methods for designing H4 are presented, and various possible mechanical architectures of the parallel manipulator are exposed. Section 3 discusses the inverse and forward kinematics problem of H4. Section 4 deals with singularity analysis of H4 utilizing line geometry tools and screw theory. Section 5 concludes this chapter by providing the development tendency of the parallel manipulators.

2. Structural synthesis and architectures

2.1. General concept of H4

Parallel manipulators are constituted of a moving platform that is connected to a fixed base by several chains (limbs). Generally, the number of limbs is equal to the degrees of freedom (DOF) of the moving platform such that each limb is driven by no more than one actuator and all actuators can be mounted on or near the fixed base. By acting on the limbs the platform pose (position and orientation) is controlled. Moreover, if the actuators are locked, the manipulator will become an isostatic structure in which all the legs carry the external loads applied to the platform. This feature makes the parallel manipulators with high stiffness possible throughout the whole space. From the first ideas proposed by Gough (Gough, 1956) or Steward (Steward, 1965), a lot of interesting mechanical devices or design methods have been extensively studied. In the beginning, many structures were based on the ingenuity of the researchers and not on a systematic approach. Subsequently, a new research domain called structure (or type) synthesis was proposed, in which vaious methodologies were tried to generate all the structures that have a desired kinematic performance. The most widely used synthesis approaches (and their variants) are graph theory, group theory and screw theory (Merlet, 2006). In this section, we will not discuss these theoretical problems, but focus on the structure generation of H4.

The 6-DOF parallel manipulators generally suffer from a small workspace, complex mechanical design and difficult motion generation and control due to their complex kinematic anslysis. To overcome these shortcomings, new structures for parallel manipulators having less than 6-DOF (which are called limited-DOF manipulators in this chapter, although many of these new structures were known well before) are explored. There is an overriding motivation behind such efforts: limited-DOF manipulators may be needed for many applications. For example, parallel wrists need only three rotational DOFs. It is clear today that most attention has been paid to 3-DOF family among the limited-DOF parallel manipulators (Gosselin & Angeles, 1989, Agrawal et al., 1995, Gosselin & St-Pierre, 1997, Fattah & Kasaei, 2000, Gregorio, 2001). However, in many industrial situations, there is a need for equipment providing more than 3 DOFs arranged in parallel and based on simpler arrangements than 6-DOF strctures. The reference (Cheung et al., 2002) developed a 4-DOF parallel manipulator E4 which could be used in a semiconductor packaging system. It is in the early 80’s when Reymond Clavel comes up with the brilliant idea of using parallelograms to build a parallel robot with three translational and one rotational degree of freedom (Bonev, 2001). Contrary to opinions published elsewhere, his inspiration was truly original and does not come from any parallel mechanism patents. Professor Clavel called his creation the Delta robot. The new family of 4-DOF parallel manipulators called H4 that could be useful for high-speed pick-and-place applications is just based on the idea of Delta structure (Pierrot & Company, 1999, Company & Pierrot, 1999). The prototype built in the Robotics Department of LIRMM can reach 10g accelerations and velocities higher than 5m/s (Robotics Department of LIRMM).

2.1.1. Delta structure

Even if an incredibly large number of different structures have been proposed by academic researchers in the last 30 years, most of these that are used widely in industry can be classfied into two basic types: the Delta structure and the so-called “hexapod“ with 6 U-P-S chains in parallel (U-P-S: Universal-Prismatic-Spherical). This may be a result of either the exceptional simplicity of the Delta 3-DOF solution, or the enormous research effort dedicated to “hexapod“ (Company, 1999).

The basic idea behind the Delta parallel robot design is the use of parallelograms. A parallelogram allows an output link to remain at a fixed orientation with respect to an input link. The use of three such parallelograms restrain completely the orientation of the mobile platform which remains only with three purely translational degrees of freedom. The input links of the three parallelograms are mounted on rotating levers via revolute joints. The revolute joints of the rotating levers are actuated in two different ways: with rotational (DC or AC servo) motors or with linear actuators. Finally, a fouth leg is used to transmit rotary motion from the base to an end-effector mounted on the mobile platform. Fig.1 shows the most famous Delta robot with three translation degrees of freedom, which was initially developed at École Polytechnique from Lausanne by Clavel (Clavel, 1988). All the kinematic chains of this robot are of the RRPaR type: a motor makes a revolute joint rotate about an axis w. On this joint is a lever, at the end of which another joint of the R type is set, with axis parallel to w. A parallelogram Pa is fixed to this joint, and allows translation in the direction parallel to w. At the end of this parallelogram is a joint of the R type, with axis parallel to w, and which is linked to the end effector.

media/image1.png

Figure 1.

The Delta structure proposed by Prof. Clavel and one of its industrial version, the CE33 (courtesy of SIG Pack Systems)

The Delta robot is firstly marked by the two Swiss brothers Marc-Olivier and Pascal Demaurex who created the Demaurex company (Bonev, 2001). The joint-and-loop graph of the Delta robot and one of its equivalent structures are shown in Fig.2, where P, R and S represent prismatic, revolute and spherical joint respectively. The displacement of the end-effector of the Delta robot is the result of the movement of the three articulated arms mounted on the base, each of which are connected to a pair of parallel rods. The three orientations are eliminated by joining the rods in a common termination and the three parallelograms ensure the stability of the end-effector. It is noted that the rotary actuator and lever part of the Delta could be replaced by a linear actuator, as suggested by Clavel itself and Zobel (Zobel, 1996). This type of Delta is sometime called a Linapod or a linear Delta.

media/image2.png

Figure 2.

Joint-and-loop graph of the Delta robot and one of its equivalent structure

2.1.2. Structure derivation of H4

H4 is based on 4 independent kinematic chains between the base and the nacelle, each chain is actuated where each actuator is fixed on the base (Pierrot & Company, 1999, Company & Pierrot, 1999). Such technological and conceptual ideas have already proven their efficiency on high-speed equipment for the Delta robot. So, knowing the advantages of this family of mechanisms, it is interesting to recall a few important design features of the Delta robot. The Delta robot is based on three actuated linear or rotational joints, and three pairs of rods equipped with ball joints at each end. As a matter of fact, two ball joints on each rod introduce an internal DOF for the rod which can rotate about its own axis. An arrangement such as the one depicted in Fig.3 suppresses these internal DOF and keeps the same global behaviour. To go a little further, it is possible to consider each pair of rods (that is: two U-S chains parallel one to the other) as equivalent to a unique rod equipped with universal joints at each extremity. The admissible motions for both chains are not exactly the same, but they can be regarded as “equivalent” in terms of number and type of degrees of motion.

media/image3.png

Figure 3.

The structure of a Delta robot with no internal DOF.

In order to define the basic principle of a 4-DOF mechanism, the classical kinematics formulas can be used. For a mechanism with a “closed loop” (a chain going from ground to the nacelle, and then back to ground), its mobility is given by:

m=i=1Nldofi6L
(1)

where L is the number of loops; dofi is the number of DOF of the ith link; Nl is the number of links.

For a fully-parallel 4-DOF mechanisms involving no passive chain, L=3 and m=4 . Thus:

i=1Nldofi=4+(6×3)=22
(2)

The four P-U-U chains provide: 4×5=20dof . Thus each additional joint must be only 1 DOF joint. Two rotational joints can be used because such joints are extremely easy to manufacture with good accuracy at low cost. So the arrangement shown in Fig.3 are created, each pair is connected to the nacelle by an R joint. Note that the actuated prismatic joints can be replaced by actuated rotational joints, and that the U-U chains can be replaced by (U-S)2 chains as well (shown in Fig.3(b)). To date, the conditions for a 4-DOF machine have been set up. The conditions to be fulfilled for obtaining 3 translations and 1 rotation about a given axis should be provided.

media/image12.png

Figure 4.

Architecture scheme of H4.

The demonstration of “quasi-equivalence” between P-U-U and P-(U-S)2 is proposed in (Company et al., 2003). The following discusses the possible motion of H4 (Pierrot & Company, 1999, Company & Pierrot, 1999).

Firstly, assume that the two rods of the (U-S)2 chains are parallel to each other. Fig.5 shows a simple scheme of a possible H4 structure. A i and B i are the centers of segments Ai1Ai2 and Bi1Bi2, respectively. Note Ai1Bi1=Ai2Bi2=AiBi and ui=Bi1Bi2. With such notations, the impossible motion for the tip of each P-(U-S)2 chain is the rotation about: AiBi×ui.

media/image13.png

Figure 5.

A simple scheme of H4 structure

The link B1B2 is connected to the ground by P-(U-S)2 chains. Each chain provides 3 translations; thus the link B1B2 can undergo the same translations. Moreover B1B2 cannot rotate about A1B1×u1, neither about A2B2×u2 (because the motion that B1B2 cannot make is the sum of the impossible motions for each chain). Consequently, B1B2 can only rotate about the vector normal to the two previous vectors, namely:

(A1B1×u1)×(A2B2×u2)
(3)

The link B3B4 has 3 DOFs in translation, and can only rotates about:

(A3B3×u3)×(A4B4×u4)
(4)

The link C1C2 is connected, on one side, to B1B2 by a revolute joint whose direction is represented by vector v1, on the other side, to B3B4 by a revolute joint whose direction is represented by vector v2. Thus, on each side of the nacelle, there are two “meta-chains” which have 5 DOFs. The first chain has 3 translations and 2 rotations about v1 and (A1B1×u1)×(A2B2×u2) respectively. The rotation about (A1B1×u1)×(A2B2×u2)×v1 is impossible. The second “mecha-chain” has 3 translations and 2 rotations about v2 and (A3B3×u3)×(A4B4×u4) respectively. The rotation about (A3B3×u3)×(A4B4×u4)×v2 is impossible. So the possible motions of the nacelle are 3 translations and 1 rotation about [(A1B1×u1)×(A2B2×u2)×v1]×[(A3B3×u3)×(A4B4×u4)×v2] . When the two revolute joints placed on the nacelle have the same direction represented by a unit vector v, the existence of a rotational motion about v can be written as follows:

α,α0[(w1×w2)×v]×[(w3×w4)×v]=αv
(5)

Where wi=AiBi×ui . Dot metrix the vector v on each side of the equation (4):

{[(w1×w2)×v]×[(w3×w4)×v]}v=αvv
(6)

Since v is a unit vector, then we can obtain

{[(w1×w2)×v]×[(w3×w4)×v]}=α0
(7)

Thus:

(w3×w4)[(w1×w2)×v]=[(w1×w2)×(w3×w4)]v=α0
(8)

So, the necessary condition that the nacelle can rotate about a given axis is:

[(w1×w2)×(w3×w4)]v0
(9)

It should be noted that the previous condition is only a necessary condition and any singular configuration is not taken into account. The following should demonstrate that the two rods of the P-(U-S)2 chain are parallel to each other when the nacelle moves, so the necessary condition (7) is also satisfied.

The rods of the P-(U-S)2 chain will stay parallel to each other if both links B1B2 and B3B4 move only in translation with respect to the ground. As a matter of fact, the possible motions of the nacelle imply that the two pivots fixed on the nacelle along the direction of vector v will keep this direction. Thus the only possible rotation for B1B2 and B3B4 is about vector v. But equation (2) indicates that B1B2 can only rotate about w1×w2 . So, as long as w1×w2 and v are not collinear, the link B1B2 has no possible rotation. A similar remark can be made for B3B4, leading to the two following conditions:

w1×w2vw3×w4v
(10)

In conclusion, if conditions (6) are fulfilled, the links B1B2 and B3B4 move only in translation, and the rods in a pair will stay parallel to each other.

2.1.3. H4 mechanism with symmetrical design

In this chapter, the symmetrical design of H4 refers to the structure scheme based on four identical P-U-U chains or P-(U-S)2 chains. The four identical chains provide 4×5=20 DOFs. Thus, two additional 1-DOF joints are needed. H4 uses two rotational joints because they are easy to manufacture with good accuracy at low cost. Two types of typical architecture scheme of H4 are shown in Fig.4. Note that the U-U chains can be replaced by (U-S)2 chains. Moreover, the actuated prismatic joints could be replaced by actuated rotational joints as well. A symmetrical mechanism with prismatic actuators is show in Fig.6 (a), while an equivalent mechanism with rotational actuators is shown in Fig.6 (b).

media/image28.png

Figure 6.

Symmetrical H4 with actuated prismatic joints and rotational joints

2.1.3. H4 mechanism with asymmetrical design

The link B1B2 has three translations and a possible rotation about a vector given by equation (2). It is possible to create more advanced mechanisms based on a particular case of equation (2). For example, if u1=u2, then equation (2) becomes:

(A1B1×u1)×(A2B2×u1)=βu1
(11)

As long as A1B1A2B2 , then β0 . So the only possible rotation of this link is about vector u 1. Thus, the first “meta-chain” already fulfilled our requirements: 3 translations and 1 rotation about a given axis. However, this “meta-chain” is only equipped with 2 actuators; thus 2 other chains (with one actuator each) are needed in parallel to lead to a fully-parallel manipulator. If the first “meta-chain” is constituted of two P-U-U chains, it provides 2×5=10 DOFs. So the two additional chains that constitute the second “meta-chain” should be 6

media/image32.png

Figure 7.

Asymmetrical design of H4.

DOFs, which are obviously compatible with the motions offered by the first “meta-chain”. Such mechanical structures are named asymmetrical H4 (Pierrot & Company, 1999, Company & Pierrot, 1999). An asymmetrical design of H4 using two P-U-Us and two P-U-S’s is shown in Fig.7 (a). The P-U-U’s can be replaced by P-(U-S)2’s, as shown in Fig.7 (b).

2.2. Evolution of H4

A new trend in the researches on parallel robotics is the development of lower mobility manipulators. Indeed, a lot of applications do not need six DOFs. A classification proposed by Brogardh (Brogardh, 2002) gives the necessary number of DOF for different industrial applications. He shows that applications such as pick-and-place, assembly, cutting, measurement, etc. just need from three up to five DOFs.

The parallel manipulator H4 can provide three translations and one rotation about a fixed axis, which are also called Schonflies motions (Hervé, 1999). SCARA robots were the first manipulators developed to produce these movements. Due their serial architecture, these robots involve high moving masses which are not suitable for high dynamics. Parallel architecture can overcome this problem. Delta robot is first introduced by Clavel to execute Schonflies motions (Clavel, 1988). This architecture is able to produce three translations and one rotational motion is obtained using a central “telescopic” leg built with universal and prismatic joints. This RUPUR chain suffers from a short service life, and involves a bad stiffness of the rotation motion. Orthoglide (Chablat, 2002) is another parallel robot that can provide Schonflies motions using linear actuators. Its particularity is to be isotropic at the center of its workspace. The machine tool HITA STT has been proposed to produce Schonflies motions (Thurneysen et al, 2002). The particularity of this architecture is to use additional parts in the traveling plate in order to amplify the rotational motion. Other Schonflies motion generators include Gross Manipulator (Angeles et al., 2000), SMG (Angeles, 2005), Kanuk and Manta (Rolland et al., 1999) robots.

In order to avoid the central telescopic leg of the Delta robot, in 1999, the Montpellier Laboratory of Computer Science, Robotics, and Microelectronics (LIRMM) in French invented the parallel manipulator H4 which used the concept of articulated traveling plate. The rotational movement is obtained using an internal mobility on the traveling plate, and a transforming device gives the desired range of rotational motion. Because placing actuators in a symmetrical way (i.e. at 90 one relatively to each other) involves “internal singularities” (Pierrot & Company, 1999), the robot has to be built using a particular arrangement of motors. This non-symmetrical arrangement entails a non-homogeneous behavior in the workspace and a limited stiffness of the robot (Company et al., 2005). Later, based on the concept of H4, an improved mechanical structure called I4 is proposed by LIRMM (Krut et al., 2004, Krut et al., 2003). The internal mobility of I4 is obtained with a prismatic joint (Fig. 8). The advantage of this architecture is to authorize a symmetrical arrangement of the actuators. As demonstrates by Krut (Krut et al., 2004), it is possible to place the actuators at 90 one relatively to each other. However, this architecture is more adapted to machine-tool application than to high speed pick-and-place. Indeed, commercial prismatic joints are not suitable for high speed and high accelerations, and have a short service life under such conditions. This inconvenient is due to the high pressure exerted on the balls of these elements at high acceleration conditions. Based on H4 and I4 architecture, an evolution of these mechanisms, named Part 4, has been developed with the wish of reaching high speeds and accelerations (Nabat et al., 2005). A paper written by Nabat (Nabat et al., 2006) presented experimental results showing that Part 4 was able to reach an acceleration of 15G. Part 4 is a parallel manipulator composed of four closed kinematic chains and an articulated traveling plate. The kinematic chains are similar for Delta, H4 and I4, each of which is composed of an arm and a spatial parallelogram (forearm) linked with spherical joints. The traveling plate is composed of four parts: two main parts (1, 2) linked by two bars (3, 4) with revolute joints (see Fig.9). Thus, its shape is a planar parallelogram and the internal mobility of this traveling plate is a circular translation. Generally, the “natural “range of the rotational operational motion is small, in order to make a complete turn: [-π; π], an amplification system has to be added on the traveling plate. Several options are available for this amplification such as gears or belt/pulleys. The prototype shown in Fig.9 (b) has been built using belt/pulleys system, with the first pulley fixed on one half traveling plate, and the second one is linked with a revolute joint to the second half traveling plate. Due to its traveling plate having the shape of a planar parallelogram, Part 4 has all the advantages of the previous robots, without their drawbacks. The traveling plate of this robot is articulated and is exclusively realized with revolute joints. Nabat (Nabat et al., 2006) has demonstrated that it is possible to have the same arrangements of the actuators as I4. Thus, this robot is well suited to reach high dynamics and, at the same time, to have a good stiffness and a homogenous behavior in the workspace.

media/image33.png

Figure 8.

I4 prototype (courtesy of LIRMM)

The use of base-mounted actuators and low-mass links allows the traveling plate to achieve great accelerations. This makes this type of parallel manipulator a perfect candidate for pick and place operations of light objects. Especially for semiconductor end-package equipments, which need high-speed and high-precision pick-and-place movement, these parallel manipulators provide a novel solution scheme.

In this chapter, only H4 is considered, because it firstly provides the most important concept of traveling plate. Based on the analytical method of H4, other types of evolution structure can be analyzed conveniently.

media/image34.png

Figure 9.

Part4 prototype (courtesy of LIRMM)

3. Kinematic analysis

This section will examine the relations between the actuated joint coordinates of H4 and the nacelle (or traveling plate) pose.

3.1. Inverse kinematics

The relation giving the actuated joint coordinates for a given pose of the end-effector is called the inverse kinematics, which is simple for parallel robots. The inverse kinematics consists in establishing the value of the joint coordinates corresponding to the end-effector configuration. Establishing the inverse kinematics is essential for the position control of parallel robots. There are multiple ways to represent the pose of a rigid body through a set of parameters X. The most classical way is to use the coordinates in a reference frame of a given point C of the body, and three angles to represent its orientation. But there are other ways such as kinematic mapping which maps the displacement to a 6-dimensional hyperquadric, the Study quadric, in a seventh-dimensional projective space. The kinematic mapping may have an interest as equations involving displacement are algebraic (and the structure of algebraic varieties is better understood than other non-linear structures) and may have interesting properties, for example, stating that a point submitted to a displacement has to lie on a given sphere is easily written as a quadric equation using Study coordinates (Merlet, 2006).

3.1.1. Analytic method

For a fully-parallel mechanism, each of the chains link the base to the moving platform. If A represents the end of the chain that is linked to the base, and B the end of the chain that is linked to the moving platform. By construction the coordinates of A are known in a fixed reference frame, while the coordinates of B may be determined from the moving platform position and orientation. Hence the vector AB is fundamental data for the inverse kinematic problem, this is why it plays a crucial role in the solution.

In order to write the position relationship of H4, consider the robotic structure (see Fig.10, chain #4 is not plotted for sake of simplicity). The geometrical parameters of the manipulator at hand are defined as follows:

Two frames are defined, namely {A}: a reference frame fixed on the base; {B}: a coordinate frame fixed on the nacelle (C 1 C 2);

  • The actuators slide along guide-ways oriented along a unitary vector, kz ( kz is the unity vector along z axis in the reference frame {A}), and the origin is the point Pi , so the position of each point Ai is given by: Ai=Pi+qizi qi is the moving distance of the actuator. The number i=1,2,3,4 represents each pair of kinematic chains;

  • The parameters Mi,d and h are the length of the rods, the offset of the revolute-joint from the ball-joint, and the offset of each ball-joint from the center of traveling plate, respectively;

  • The position of the end effector, namely point B , is defined by a position vector B=[xyz] , and a scalar θ , representing the orientation angle.

Without losing generality, in this section we consider Pi=(ai,bi,0) for simplicity, Qz is the offset of {A} from {B} along the direction of kz .

media/image48.png

Figure 10.

H4 manipulator with four lines drives and (S-S)2 chains

The position of points C1 and C2 are given by:

C1=B+Rot(v,θ)(BC1)C2=B+Rot(v,θ)(BC2)
(12)

Where Rot(v,θ) denotes the rotation matrix of nacelle with respect to reference frame. The position of each point B i is given by:

B1=C1+C1B1B2=C1+C1B2B3=C2+C2B3B4=C2+C2B4
(13)

The position relationship can then be written as:

AiBi2=Mi2
(14)

Then for the first axis:

A1B1=[B+Rot(v,θ)(BC1)+C1B1P1]q1z1(A1B1)2=q122q1d1z1+d12
(15)

where d1=B+Rot(v,θ)(BC1)+C1B1P1 . Finally, the two solutions are given by:

q1=d1z1±(d1z1)2+M12d12
(16)

Similar derivations give the solution for q 2, q 3 and q 4.

3.1.2. Geometrical method

The geometrical approach to the inverse kinematics problem is to consider that the extremities A, B of each chain have a known position in 3D space. The leg can be cut at a point M and two different mechanisms M A, M B constituted of the chain between A, M and the chain between B, M can be gotten. The free motion of the joints in these two chains will be such that point M, considered as a member of M A, will lie on a variety V A, while considered as a member of M B it will lie on a variety V B. If assume the mechanism have only classical lower pairs, these varieties will be algebraic with dimensions d A, d B. In the 3D space, a variety of dimension d is defined through a set of 3-d independent equations, and hence V A, V B will be defined by 3-d A and 3-d B equations. The solutions of the inverse kinematic problem lie at the intersection of these varieties. As the number of of solutions must be finite (otherwise the robot cannot be controlled), the rank of the intersection variety must be 0. In other words, in order to determine the 3 coordinates of the points, 3-d A+3-d B should equal to 3 or d A+d B equal to 3 (Merlet, 2006).

The key problem about the geometrical mthod rests with the choice of the cutting point. For the H4 structure shown in Fig.10, B i are chosen as the cutting points. Points B i has to lie on a sphere centered at A i with radius M i, while for the nacelle, the coordinates of Points B i can be described as (8). Hence to obtain the intersection of these 2 varieties the known distance between A, B of should be equal to M i: this equation will give the joint coordinates q i. The same results can be obtained as described as (11).

3.2. Direct kinematics

Direct kinematics addresses the problem of determining the pose of the end-effector of a parallel manipulator from its actuated joint coordinates. This relation has a clear practical interest for the control of the pose of the manipulator, but also for the velocity control of the end-effector.

In general, the solutions for determining the pose of the end-effector from measurements of the minimal set of joint coordinates that are necessary for control purposes is not unique. There are several ways of assembling a parallel manipulator with given actuated joint coordinates, and generally the direct kinematic relationship cannot be expressed in an analytical manner. Although various methods have been presented for finding all the solutions for this problem and their computation times are decreasing (Faugère, 1995, Hesselbach & Kerle, 1995, Gosselin, 1996, Merlet, 2004), it is still difficult to use these algorithms in a real time context. Furthermore, there is no known algorithm that allows the determination of the current pose of the platform among the set of solutions. The numerical methods using a-priori information on the current pose are more compatible with a real-time context, while their convergence and robustness is an important issue. As direct kinematics is an important issue for control, it may be necessary and interesting to investigate how to improve the computation time. Besides the progress in algorithms and processor speed, one possible approach to solve these problem is to add sensors (i.e. to have more than n sensors for a n-DOF manipulator) to obtain information, allowing a faster calculation of the current pose of the platform, at the cost of more complex hardware (Bonev et al., 2001; Chui & Perng; 2001; Parenti-Castelli, 2001). Merlet has pointed out several unsovled problems about the parllel robot’s direct kinematics (Merlet, 2000). Especially the algebraic geometry and computational kinematics should bring about enough attention.

The first step of solving direct kinematics is to determine a bound on its maximal number of solution. Then, the equations are reduced for the system to obtain the solution of a univariate polynomial whose degree should be determined in the previous analysis. Many researchers have tried to obtained closed-form solutions of parallel robots in a univariate polynomial equation. Especially, the 3- and 6- DOFs parallel robots, e.g., planar, Delta, Steward platform and Hexapod robots, have been mainly considered. This section provides a univariate polynomial equation for H4. It is shown that the solutions of the forward kinematics yield a 16th degree polynomial in a single variable (Choi et al., 2003).

3.2.1. Forward kinematics formulation

The forward kinematics of the H4 is the problem of computing the position and orientation of the nacelle (traveling plate) form the motor angles or the moving distance of the actuator. To get a closed-form solution, a univariate polynomial equation needs to be solved. The following method is extracted from the paper written by Choi (Choi et al., 2003).

The kinematic structure of H4 is shown in Fig.10 and design parameters are shown in Fig.11. The nacelle is composed of three parts (two lateral bars and one central bar). The four points B 1, B 2, B 3 and B 4 form a parallelogram. Let BAi and AAi represent the homogeneous coordinates of the points B i in {A} and A i in {A} respectively. Then the homogeneous coordinates of BAi and AAi can be written as:

BAi=[x+hE1icosθy+hE1isinθ+dE2iz1]=[Piix+xPiiy+y+Qiz1]AAi=[aibiqi1]T
(17)

where

ix=cosθ iy=sinθ Pi=hE1i Qi=dE2i E11=E12=1 E13=E14=1 E21=E24=1 ,
media/image67.png

Figure 11.

Design parameters of H4

The kinematic closure of each elementary chain can be written as:

AAiBAi2=Mi2
(19)

Consequently, the each chain provides us with the following equation:

(Piix+xai)2+(Piiy+y+Qibi)2+(zqi)2=Mi2
(20)

Expanding and rearranging equation (14), the following equation is obtained:

A˜ix+B˜iy+C˜iz+D˜i+E˜i=0
(21)

where, A˜i=2(Piixai) B˜i=2(Piiy+Qibi) C˜i=2qi D˜i=ai2+bi2+qi22Pi(aiix+biiyQiiy)2Qibi+Pi2+Qi2Mi2

E˜i=x2+y2+z2
(22)

From equation (15), the four equations become:

A˜1x+B˜1y+C˜1z+D˜1+E˜1=0
(23)
A˜2x+B˜2y+C˜2z+D˜2+E˜2=0
(24)
A˜3x+B˜3y+C˜3z+D˜3+E˜3=0
(25)
A˜4x+B˜4y+C˜4z+D˜4+E˜4=0
(26)

Choosing any three from the above four equations, we can eliminate x 2y 2 and z 2 simultaneously and obtain an equation with variables x, y and z. Subtracting equation (18) from equation (17) and equation (19) from equation (17) respectively, the following two equations are obtained, which are used to eliminate x and y to obtain a polynomial in a single variable z.

ΔA23x+ΔB23y+ΔC23z+ΔD23=0
(27)
ΔA24x+ΔB24y+ΔC24z+ΔD24=0
(28)

where, ΔA23=(A˜2A˜3) ΔB23=(B˜2B˜3) ΔC23=(C˜2C˜3) ΔD23=(D˜2D˜3) ΔA24=(A˜2A˜4) ΔB24=(B˜2B˜4) ΔC24=(C˜2C˜4)

ΔD24=(D˜2D˜4)
(29)

From equation (20) and (21), x and y are solved as:

x=Δ11z+Δ12Δ0=e1z+e2
(30)
y=Δ21z+Δ22Δ0=e3z+e4
(31)

where, e1=Δ11Δ0 e2=Δ12Δ0 e3=Δ21Δ0 e4=Δ22Δ0 , with Δ0=ΔA23ΔB24ΔA24ΔB23 Δ11=ΔB23ΔC24ΔB24ΔC23 Δ12=ΔB23ΔD24ΔB24ΔD23 Δ21=ΔA24ΔC23ΔA23ΔC24 ,

Δ22=ΔA24ΔD23ΔA23ΔD24
(32)

Substituting equations (22) and (23) into equation (17), the following quadratic equation is obtained:

λ0z2+λ1z+λ2=0
(33)

where, λ0=e12+e32+1 λ1=2e1e2+2e3e4+A˜2e1+B˜2e3+C˜2

λ2=e22+e42+A˜2e2+B˜2e4+D˜2
(34)

From equation (24), z can be calculated as:

z=λ1+ρ2λ0
(35)

where,

ρ=±λ124λ0λ2
(36)

Substituting equations (22), (23) and (25) into equation (16), the following equation is obtained:

ρ22(λ1λ1)ρ+λ122λ1λ1+4λ0λ2=0
(37)

where, λ1=2e1e2+2e3e4+A˜1e1+B˜1e3+C˜1

λ2=e22+e42+A˜1e2+B˜1e4+D˜1
(38)

The variable ρ can be calculated from equation (27):

ρ=(λ1λ1)±λ124λ0λ2
(39)

So the right parts of equations (26) and (28) are equivalent:

(λ1λ1)±λ124λ0λ2=±λ124λ0λ2
(40)

By taking the second power of the both sides of equation (29), a polynomial equation with variables ix and iy can be obtained as follows:

λ0[Δλ1(λ1λ2λ1λ2)+λ0Δλ22]=0
(41)

where, Δλ1=λ1λ1 Δλ2=λ2λ2 . Equation (20) can be rewritten as (Choi et al., 2003):

f(ix8,ix7,ix7iy,ix7iy,ix6,ix6iy,ix6iy2,ix5,ix5iy,ix5iy2,ix5iy3,ix4,ix4iy,ix4iy2,ix4iy3,ix4iy4,ix3,ix3iy,ix3iy2,ix3iy3,ix3iy4,ix3iy5
(42)
ix2,ix2iy,ix2iy2,ix2iy3,ix2iy4,ix2iy5,ix2iy6,ix,ixiy,ixiy2,ixiy3,ixiy4,ixiy5,ixiy6,ix0iy0,ix0iy,ix0iy2,ix0iy3,ix0iy4)=0
(43)

where f(x1,x2,xn) represents a linear combination of x1,x2,xn . The coefficients of equation (30) are all constant quantities which depend on the geometric parameters of the H4 robot.

Using the following well-known trigonometric function:

ix=1T21+T2,ix=2T1+T2
(44)

where T=tan(θ/2) , then equation (31) becomes a polynomial equation in a single variable T which contains terms up to the 16th power. This means that the mechanism may have up to 16 different configurations for a given set of actuator’s position. Substituting equation (32) into equations (26), (22), (23) and (25), the values of ρ,x,y and z can be gotten respectively.

3.2.2. Numerical example

In this section, a numerical example for the H4 robot shown in Fig.10 and Fig.11 is presented to illustrate the application of the method proposed in the previous section. The design parameters of H4 is chosen as: h=d=6, Q z=42, v=[0 0 1]T, M 1=M 2= 3 Qz, M3=M4= 2 Qz, the origin coordinates of points Ai and Bi are set as following:

A1=[48480] B1=[6642] A2=[48480]
B2=[6642]
(45)
A3=[6480] B3=[6642] A4=[4860]
B4=[6642]
(46)
u=[u1u2u3u4]=[111011010000] , where

When the position and orientation of the nacelle are set x=5, y=-5, z=-30, θ=π/6(rad), four actuators’ position which is calculated by inverse kinematics are given by:

q1=13.021,q2=7.488,q3=9.679,q4=15.770
(48)

For the actuators’ position given by equation (33), Matlab Symbolic Math Toolbox is used to solve the polynomial equation with variable T. The results indicate that there are only twelve solutions, the author believes that there are several repeated roots. This result needs further study. The two real roots are T=0.268 and T=-2.882 respectively. When set T=0.268, the configuration of the nacelle is calculated as: x=5mm, y=-5mm, z=-30mm, θ=0.524(rad), which is the expected solutions. When set T=-2.882, the configuration of the nacelle is calculated as: x=3mm, y=-7.95mm, z=-14.58mm, θ=-2.47(rad). This configuration cannot be realized in practice because of the mutual interference of the parallelograms as shown in Fig.11.

3.3. Jacobian matrix

Jacobian matrix relates the actuated joint velocities to the end-effector cartesian velocities, and is essential for the velocity and trajectory control of parallel robots. For parallel manipulators, their inverse Jacobian matrix can be established without a very high complexity, but their Jacobian matrix cannot be obtained directly, even with the help of symbolic computation, except in some particular cases (Bruyninckx, 1997; Pennock & Kassner, 1990). Theoretical analytic formulations of jacobians have been proposed, but require complicated matrix inversions (Dutré et al., 1997; Kim et al., 2000). Generally, the difficulty of the inversion does not lie in the complexity of the algorithm but in the sheer size of the result (Merlet, 2006).

The velocity of the nacelle can be defined by resorting to a velocity for the translation, vB=[x˙y˙z˙] and a scalar for the rotation about v, θ˙ . Thus, the velocity of points C1 and C2 can be written as follows (Pierrot & Company, 1999):

vC1=vB+θ˙v×BC1
vC2=vB+θ˙v×BC2
(49)

Moreover, since the links B1B2 and B3B4 move only in translation, the following relations hold:

vB1=vB2=vC1

On the other hand, velocity of points Ai is given by:

where zi is the unit vector along the direction of prismatic. The velocity relationship can then be written thanks to the classical property:

vAiAiBi=vBiAiBi
(52)

Equation (34) can be written for i=1,…,4 and the results grouped in a matrix form, such as:

where

Jq=[(A1B1z1)0000(A2B2z2)0000(A3B3z3)0000(A4B4z4)]
(54)
q˙=[q˙1q˙2q˙3q˙4]T
(55)
Jx=[(A1B1)x(A1B1)y(A1B1)z(A1B1×BC1)v(A2B2)x(A2B2)y(A2B2)z(A2B2×BC1)v(A3B3)x(A3B3)y(A3B3)z(A3B3×BC2)v(A4B4)x(A4B4)y(A4B4)z(A4B4×BC2)v]
(56)
x˙=[x˙y˙z˙θ˙]
(57)

If the mechanism is not in a singular configuration, the inverse Jacobian matrix is described as:

J1=Jq1Jx
(58)

And the Jacobian matrix is:

When the determinant of inverse Jacobian matrix equals to zero, the parallel manipulator is in a singular configuration, which is a general method for identifying the singular configurations of parallel manipulators. But for a manipulator with n<6DOF, it should be noted that it may not be efficient to identify all the singular configurations by determining only the n×n inverse kinematic Jacobian that relates the actuated joint velocities to the possible DOF velocities. Next section will discuss an inverse kinematic Jacobian that involves actuated joints velocities and the full twist of the end-effector, which may be essential for singularity analysis. This type of matrix is coined as overall jacobian by Joshi (Joshi & Tsai, 2002), while Merlet called it a full inverse kinematic jacobian (Merlet, 2006).

3.4. Determination of the joint velocities and twist

The purpose of this section is to calculate the actuated joint (or end-effector) velocities, being given the end-effector (or actuated joint) velocities.

In the previous section, equation (35) can be applied to compute the actuated joint velocities directly. For example, the structure and design parameters of H4 is described in Fig.10 and Fig.11, where the following parameters are chosen as described in section 3.2.2. So, according to equation (35), the inverse Jacobian matrix can be expressed as:

J1=[1116111601161010]
(60)

This inverse Jacobian matrix can be used to calculate the actuated joint velocity directly when the pose of the nacelle is shown in Fig.11.

It is usually difficult, even for robots with less than 6-DOFs robots, to invert J-1 analytically. So a numerical procedure will be needed to calculate the twist of the nacelle from the joint velocities. For a given pose of the end-effector, there are two methods to determine the Jacobian matrix, namely, a numerical inversion algorithm and the quasi-Newton scheme (Reboulet, 1985, Merlet, 2006).

Jk+1=Jk+J0(q˙J1Jk)
(61)

where q˙ is the actuated joint velocity; J0 denotes the kinematic jacobian matrix calculated for a given pose. The algorithm stops when the differences between the joint velocities and those that are calculated from the twist are lower than a fixed threshold.

For the H4 structure and design parameters described above, when the actuated joing velocity is q=[q1q2q3q4]=[1111] , the speed of the nacelle can reach x˙=[x˙y˙z˙θ˙]=[0.50.250.50.125]

3.5. Accelerations analysis

This section will determine what are the relations between the actuated joint accelerations and the cartesian and angular accelerations of the end-effector. H4 presents excellent characteristics as to acceleration which can reach 10g. For parallel robots it is generally easy to obtain these relations directly. From the equation q˙=J1x˙ , the following can be obtained by differentiation

q¨=J1x¨+J˙1x˙
(62)

For the various categories of parallel manipulators, the determination of the acceleration equations thus amounts to the determination of the derivative of the inverse kinematic jacobian matrix (Merlet, 2006).

But for H4, the method described above is not intuitionistic. So we derive the relations between the actuated joint accelerations and the accelerations of the end-effector using geometric method. From equation (12), accelerations of the point Bi and Ai in {A} can be obtained as:

aBi=[x¨hE1icosθθ¨y¨hE1isinθθ¨z¨]TaAi=[00q¨i]T
(63)

The acceleration projection of points Bi and Ai on the line AiBi is equal. So the following equation comes into existence:

aAiAiBi=aBiAiBi
(64)

Equation (40) can be written for i=1,…,4 and the results grouped in a matrix form, such as:

where

q¨=[q¨1q¨2q¨3q¨4]T
(66)
Ja=[(A1B1)x(A1B1)y(A1B1)z((A1B1)xcosθ+(A1B1)ysinθ)hE11(A2B2)x(A2B2)y(A2B2)z((A2B2)xcosθ+(A2B2)ysinθ)hE12(A3B3)x(A3B3)y(A3B3)z((A3B3)xcosθ+(A3B3)ysinθ)hE13(A4B4)x(A4B4)y(A4B4)z((A4B4)xcosθ+(A4B4)ysinθ)hE14]
(67)
x¨=[x¨y¨z¨θ¨]
(68)

For the H4 structure and design parameters described above, when the nacelle is in the origin pose and the actuated joint accelerations are q¨=[10101010] m/s2, the linear acceleration of the nacelle can reach 23m/s2 and the angular acceleration can reach 2.5rad/s2.

3.6. Kinetostatic performance indices

3.6.1. Manipulability

The inverse kinematic jacobian matrix reflects a linear relation between the manipulator accuracy Δx and the measurement errors Δq on q . If the measurement errors are bounded, then the hyper-sphere in the joint error space can be mapped into an ellipsoid in the generalized Cartesian error space with using the Euclidean norm. This ellipsoid is usually called the manipulability ellipsoid. If using the infinity norm, the joint errors are restricted to lie in a hyper-cube in the joint error space. The hyper-cube in the joint error space can be mapped into the kinematic polyhedron in the positioning errors space (Merlet, 2006). The shape and volume of the manipulability ellipsoid and the kinematic polyhedron characterize the manipulator dexterity. It is necessary to set up some kinetostatic performance indices that are used to quantify the dexterity of a robot. The famous Yoshikawa’s manipularity index JJT is used for serial robots for a long time. Another index is the condition number that is often used for parallel robots.

κ=J1J
(69)

The condition number expresses how a relative error in joint space gets multiplied and leads to a relative error in work space. The condition number is dependent on the choice of the metric norm and the most used norms are the 2-norm and the Euclidean norm. The condition number is mentioned as the main index for characterizing the accuracy of parallel robots. But there is major drawback to the condition number for H4, since the elements of the matrix corresponding to translations are dimensionless, whereas those corresponding to the rotations are lengths. A direct consequence is that the condition number has no clear physical meaning, as the rotations are transformed arbitrarily into “equivalent” translations (Merlet, 2006). There are various proposals that have been made to avoid this drawback (Gosselin, 1990; Kim & Ryu, 2003). But all these proposals have their own special drawbacks. How to design a general law for parallel robots’ manipulability is still an open problem.

3.6.2. Isotropy

Poses with a condition number of 1 are called isotropic poses, and robots having only such type of poses are called isotropic robots (Merlet, 2006). At the given pose as shown in Fig.11, the condition number of H4 is about 8.56 and there is no any isotropic poses in the workspace of H4 (Company et al., 2005). While Part4 robots which are the evolutional structure of H4 are isotropic at their original poses.

There are some other manipulability and accuracy indices, such as global conditioning index, uniformity of manipulability and maximal positioning error index et al. All the above definitions of the kinetostatic indices do not take into account other factor affecting the accuracy of parallel robots, such as manufacturing tolerances, clearance and friction in the joints. In order to include these factors in accuracy analysis, Monte Carlo statistical simulation technique can be used to evaluate the accuracy of parallel manipulators (Ryu & Cha, 2003).

4. Singularity analysis

This section will identify all the singular configurations of the H4 and analyze the manipulator’s self-motion when in or closed to a singular configuration. Generally, singular configurations refer to particular poses of the end-effector, for which parallel robots lose their inherent infinite rigidity, and in which the end-effector will have uncontrollable degrees of freedom. But for this new family of parallem manipulator H4, singularities are associated with either loss or gain of DOF. This section utilizes line geometry tools and screw theory to deal with singularity analysis of H4. Firstly, the basic theory including in the process of singularity analysis is introduced briefly. Then the static equilibrium condition of the end-effector is derived to obtain the full inverse kinematic jacobian 6×6 matrix, which is set of governing lines of the manipulator. Based on linear dependency of these lines, the singular configurations of the manipulator can be identified. Moreover, in order to deal with singularities associated with loss of DOFs (serial singularity), the static equilibrium of the actuators is also defined. Secondly, architecture and constraint singularities associated with gain of DOFs (parallel singularity) are defined and analyzed using linear complex approximation algorithm (LCAA), which is employed to obtain the closest linear complex, presented by its screw coordinates, to the set of governing lines. The linear complex axis and pitch provide additional information and a better physical understanding of the manipulator’s self-motion when in or closed to a singular configuration. Lastly, various singularities of an example H4 manipulator are presented and analyzed using the proposed methods (Wu et al., 2006).

4.1. Basic theory

In the context of designing a parallel manipulator, understanding the intrinsic nature of singularities and their relations with the kinematic parameters and the configuration spaces is of great importance. The phenomena of singularity in parallel manipulators have been approached from different points of view. One way to introduce singular configurations is to examine the relations obtained for inverse kinematics. Singularities correspond to the roots of the parallel manipulator Jacobian’s determinant. Based on the rank deficiency associated with the Jacobian matrices, singularities of parallel manipulators have been explained (Gosselin &Angeles, 1990, Zlatanov et al., 1994). Using the linear decomposition method and co-factor expansion, St-Onge and Gosselin (St-Onge & Gosselin, 2000) studied the singularity loci of the Gough-Steward platform, and obtained a graphical representation of these loci in the manipulator workspace. Due to the complexity of the kinematic model, several authors proposed numerical methods to analyze the singularities (Feng-Cheng &Haug, 1994, Funabashi & Takeda, 1995). However, even with symbolic computation software, the calculation of the determinant of Jacobian matrix is a complicated task. Moreover, identifying kinematic singularity through matrices does not provide physical insight into the nature of singular configuration of parallel manipulators.

A parallel manipulator is naturally associated with a set of constraint functions defined by its closure constraints (geometric relations of the closed-chain mechanism). The differential forms arising from these constraint functions completely characterize the geometric properties of the manipulator. So Liu et al. (Liu et al., 2003) used differential geometric tools to study singularities of parallel mechanisms and provided a finer classification of singularities. In their works, they classified singularity into configuration space singularities and parametrization singularities which include actuator and end-effector singularities as their special cases. But the method is too abstract and complicated when dealing with singularities of spatial parallel manipulators.

Another approach to analyze the parallel manipulators’ singularity is based on line geometry and screw theory. Merlet (Merlet, 1992) considered the Jacobian matrix of the Gough-Steward platform as the Plücker line coordinates of the robot’s actuators and identified the singular configurations of the robot. In (Romdhane et al., 2002), a mixed geometric and vector formulation is used to investigate the singularities of a 3-translational-DOF parallel manipulator. In (Collins, 1995), singularities of an in-parallel hand controller for force-reflected teleoperation were analyzed, the six columns of the Jacobian matrix were viewed as zero-pitch wrenches (lines) acting on the top platform, then line geometry and rank determining geometric constructions were used to obtain all configuration singularities. Basu and Ghosal (Basu & Ghosal, 1997) presented a geometric condition to deal with the singularity analysis associated with gain of DOF in a class of platform-type, multi-loop spatial manipulators. Joshi and Tsai (Joshi & Tsai, 2002) developed a methodology for the Jacobian analysis of limited-DOF parallel manipulators by making use of the theory of reciprocal screws. A 6×6 Jacobian matrix so derived provided information about both architecture and constraint singularities.

The present investigation identifies all the singular configurations of the H4 robot and provides physical insight into the nature of these singular configurations from the view point of geometry. Moreover, the behavior in these singular configurations or in the neighboring ones is also determined using linear complex approximation. All these analyses need some basic theories, include Grassmann geometry, Line geometry and Screw theory etc. The following gives a brief presentation of these theories.

4.1.1. Line geometry and Plücker coordinates

Line geometry investigates the set of lines in three-space. The ambient space can be a real projective, affine, Euclidean or a non-Euclidean space. Line geometry possesses a close relation to spatial kinematics (Bottema & Roth, 1990) and has therefore found applications in mechanism design and robot kinematics. A parallel manipulator has a singular position if the axes of its hydraulic cylindrical legs lie in a linear complex, a special three-parameter set of lines. In practice, several sources for errors (manufacturing, material properties, computing, etc.) can hardly be avoided. Thus, the question is whether the lines on the objects near their realization are close – within some tolerance – to the lines of a linear complex. This is an approximation or regression problem in line space (Helmut et al., 1999).

In real Euclidean space E3, a Cartesian coordinate system is used and the points are represented by their coordinate vectors M. A straight line L can be determined by a point pL and a normalized direction vector l of L, i.e. l=1 . To obtain coordinates for L, one forms the moment vector

with respect to the origin. If p is substituted by any point q=p+λl on L, equation (43) implies that l¯ is independent of p on L. The six coordinates (l,l¯) with

l=(l1,l2,l3) , and
l¯=(l4,l5l6)
(71)

are called the normalized Plücker coordinates of L. The set of lines in E3 is a four-dimensional manifold and accordingly the six Plücker coordinates satisfy two relations. One is the normalization and the other is, by equation (43), the so-called Plücker relation

which expresses the orthogonality between l and l¯ . Conversely, any six-tuple (l,l¯) with l=1 , which satisfies the Plücker relation ll¯=0 represents a line in E3. As the orientation are not concerned, (l,l¯) and (l,l¯) describe the same line L.

The topic about the Klein mapping and special sets of lines can refer to the paper written by Pottmann et al. (Pottmann et al., 1999).

4.1.2. Grassmann geometry

As shown in the previous section, two lines with Plücker coordinates l1=[p1,q1] l2=[p2,q2] intersect if and only if p1q1+q2p2=0 . Plücker vectors with p=0 do not represent real lines and are associated with a line at infinity. All lines at infinity belong to a plane, the plane at infinity. A point may also be represented by the Plücker coordinates (α, r) so that its coordinates are r/α. If α=0, then the point is at infinity, and a point (0, r) at infinity is on the line at infinity (0, q) if and only if r q=0. Consequently a point at infinity that belongs to the two lines at infinity (0, s1), (0, s2) has coordinates (0, s1×s2).

The columns of the full inverse kinematic jacobian matrices of most parallel robots are constructed from the Plücker vectors of lines associated with links of the manipulator. The singularity of this matrix therefore means that there will be a linear dependence between these vectors. Grassmann showed that linear dependence of Plücker vectors induced geometric relations between the associated lines, so that a set of n Plücker vectors creates a variety with dimension m<n. A thorough introduction to Grassmann geometry please refers to (Pottmann et al., 1999). The applications of Grassmann geometry can be found in references (Collins & Long, 1995, Hao & McCarthy, 1998, Wu et al., 2006).

Using the Grassmann’s geometrical conditions, an algorithm for finding the singular configurations of any type of parallel robot whose full inverse kinematic jacobian consists of Plücker vectors should be designed. Firstly consider all sets of n lines that are associated with the Plücker vectors constructed from full inverse kinematic jacobian matrices, and then the pose of the moving platform is determined so that the n lines satisfy one of the geometrical conditions which ensure that they span a variety of dimension n-1, thereby leading to a singularity of the robot. The following list the geometric conditions that ensure that the dimension of the variety spanned by a set of n+1 Plücker vectors is n, for each possible dimension n of the variety (Merlet, 2006). It should be noted that the notation of intersecting lines also include parallel lines which intersect at infinity. To consider the reading convenience, the following results is excerpted from the monograph written by Merlet (Merlet, 2006).

media/image187.png

Figure 12.

Geometrical conditions that characterize the varieties of dimension 1, 2 and 3 (Merlet, 2006)

For the variety of dimension 1 (called a point) there is just one Plücker vector and one line. A variety of dimension 2, called a line, may be constituted either by two Plücker vectors for which the associated lines are skew, i.e. they do not intersect and they are not parallel, or be spanned by more than two Plücker vectors if the lines that are associated with the vectors form a planar pencil of lines, i.e. they are coplanar and possess a common point (possibly at infinity, to cover the case of coplanar parallel lines). A variety of dimension 3, called a plane, is the set of lines F that are dependent on 3 lines F1, F 2, F 3. It is possible to show that all the points belonging to the lines F lie on a quadric surface Q. This quadric degenerates to a pair of planes P1, P2 if any two of the three lines intersect.

  • condition 3d: all the lines are coplanar, but do not constitute a planar pencil of lines; F1, F 2, F 3 are coplanar and P1, P2 are coincident.

  • condition 3c: all the lines possess a common point, but they are not coplanar; F1, F 2, F 3 intersect at the same point, possibly at infinity (this covers the case of parallel lines).

  • condition 3b: all the lines belong to the union of two planar pencils of non coplanar lines that have a line L in common; two of the lines intersect at a point p, and L intersect the last line at a. Two different cases may occur:

  • P1, P2 are distinct and intersect along the line L. The set of dependent lines are the lines in P1 that go through a, and the lines in P2 that go through p

  • P1, P2 are distinct and parallel. This occurs if two of the lines Fi are parallel; L is a line at infinity, and the set of dependent lines are two planes of parallel lines

  • condition 3a: all the lines belong to a regulus; F1, F 2, F 3 are skew.

  • A variety if dimension 4, called a congruence, corresponds to a set of lines which satisfies one of the following 4 conditions:

  • condition 4d: all the lines lie in a plane or meet a common point that lies within this plane. This is a degenerate congruence.

  • condition 4c: all the lines belong to the union of three planar pencils of lines, in different planes, but which have a common line. This is a parabolic congruence.

  • condition 4b: all the lines intersect two given skew lines. This is a hyperbolic congruence.

  • condition 4a: the variety is spanned by 4 skew lines such that none of these lines intersects the regulus that is generated by the other three. This is an elliptic congruence.

A variety of dimension 5, called a linear complex, is defined by two 3-dimensional vectors (c,c¯) as the set of lines L with Plücker coordinates (l,l¯) such that c¯l+cl¯=0 . The complex may be:

– condition 5b: all the lines of the complex intersect the line with Plücker coordinates (c,c¯) , namely, cc¯=0 . This is a singular configuration.

– condition 5a: cc¯0 . This is a general or non singular configuration.

The degree of freedom associated with a linear complex is a screw motion with axis defined by the line with Plücker vector (cc¯pc)/c , where p=cc¯/c is the pitch of the motion.

media/image195.png

Figure 13.

Geometrical conditions that characterize the varieties of dimension 4 and 5 (Merlet, 2006)

4.1.3. Screw theory

According to screw theory, the instantaneous velocities of a rigid body can be described by a 6-dimensional vector of the form (ω,v) , where ω is the angular velocity of the rigid body, and v is its translational velocity. These elements are called velocity twists or screws. Forces and torques exerted on the rigid body are important for motion and may be represented as a couple of 3-dimensional vectors (F,M) called a wrench. A twist and a wrench will be said to be reciprocal if

Fv+Mω=0
(73)

When a kinematic chain is connected to a rigid body the key point is that the possible instantaneous twists for the rigid body are reciprocal to the wrenches imposed by the kinematic chains (called the constraint wrenches). In other words, the DOF of the rigid body are determined by the constraint wrenches. So based on linear dependency of the constraint wrenches, the singular configurations can be identified. Another function of the screw theory is the mobility analysis (Li & Huang, 2003).

4.2. Static analysis of the H4

When deriving the Jacobian matrix of the H4 using the velocity-equation formulation described in section 3.3, the result is a 4×4 Jacobian matrix because of the 4-DOF of the manipulator. However, it is not clear as to this is the best way to express the Jacobian of this type of limited-DOF parallel manipulator when singular analysis is processed. Because this approach is valid for general-purpose planar or spatial parallel manipulators, for which the connectivity of each serial chain (limb) is equal to the mobility of the end effector, it is not necessary true for parallel manipulators with less than 6-DOF (Joshi & Tsai, 2002). For example, this approach leads to a 3×3 Jacobian matrix for the 3-UPU parallel manipulator assembled for pure translation (Tsai, 1996, , Tsai & Joshi 2000). However, such a 3×3 Jacobian matrix cannot predict all possible singularities. Di Gregorio and Parenti-Castelli (Di Gregorio & Parenti-Castelli, 1999) analyzed the singularities of the 3-UPU translational platform. They derived the conditions for which the actuators cannot control the linear velocity of the moving platform, generally known as architecture singularities. Bonev and Zlatanov also found that under certain configurations the 3-UPU translational platform would gain rotational DOFs, which was called constraint singularities (Bonev and Zlatanov, 2001). Constraint singularities occur when the limbs of a limited-DOF parallel manipulator lose their ability to constrain the moving platform to the intended motion.

The 4×4 Jacobian matrix of the H4 cannot predict all possible singularities. Under certain configurations the chains of the H4 parallel manipulator lose their ability to constrain the moving platform to the intended motion, and the nacelle will gain additional DOF(s). So it is necessary to derive a full 6×6 matrix of the robot. In order to obtain the full 6×6 matrix, including the moments of constraints, one can express the set of static equilibrium conditions of the nacelle. These expressions result in a 6×8 matrix that maps the external wrench acting on the nacelle to the moments of the internal forces that are generated by the (S-S)2 chains. Taking into account the constraint relations between the (S-S)2 chains, we can derive a 6×6 matrix. This matrix’s rows describe the six governing lines, from which the manipulator’s singularities are obtained.

The structure and the geometrical parameters of the manipulator are shown in Fig.10. The static analysis of H4 manipulator can be divided into two parts. One regards the static equilibrium of the nacelle with B1B2 and B3B4 . The second deals with the static equilibrium conditions of the parallelograms which are composed of P-(S-S)2 chains. For a P-(S-S)2 chain, Ai1Bi1=Ai2Bi2=AiBi Ai and Bi are centers of segments Ai1Bi1 and Ai2Bi2 . Every chain has 5 DOFs, and provides the same type: 3 translations and 2 rotations. Since both rods Ai1Bi1 and Ai2Bi2 are considered as plain solids, the impossible motion is the rotation about the vector: AiBi×ui , where ui=Bi1Bi2 . Observing the structure of the parallel part of the (S-S)2 chain in Fig.2, one can detect that due to the spherical joint, there are moments exerted by the rods to B1 (the direction is along AiBi×ui ) and forces to B1 (the direction is along BiAi ), the static equilibrium for the ith chain is given by (see Fig.14):

media/image214.png

Figure 14.

Forces and moments transmitted to the nacelle

i=14fiSiFe=0i=14RABri×Sifi+i=14miSniMe=0
(74)

where RAB is a rotation matrix, which expresses the orientation of the nacelle with respect to the fixed reference frame; ri is a vector connecting B to Bi Si and fi are the unit direction vectors and the internal forces acting on Bi in the direction of BiAi respectively; Sni and mi are the unit direction vectors and the internal moments acting on the Bi in the direction of AiBi×ui . Writing equation (46) in a matrix form yields:

[S1S2S3S40000RABr1×S1RABr2×S2RABr3×S3RABr4×S4Sn1Sn2Sn3Sn4][f1f2f3f4m1m2m3m4]=[FeMe]
(75)

the matrix in the left of equation (47) is a 6×8 matrix. In order to obtain the 6×6 matrix, the constraint relations between fi,mi(i=1,2,3,4) should be analyzed.

Link C1C2 is connected, on one side, to B1B2 by a revolute joint whose direction is represented by vector v1 , on the other side, to B3B4 by a revolute joint whose direction is represented by vector v2 . Thus, we can consider two 5-DOF “meta-chains” on each side of the nacelle. Those “meta-chains” have the following properties (Pierrot & Company, 1999):

    So, the moment’s direction exerted by B1B2 to C1C2 is along (Sn1×Sn2)×v1 , and the moment’s direction exerted by B3B4 to is along (Sn3×Sn4)×v2 . The static equilibrium of C1C2 is given by (see Fig.15):

    media/image236.png

    Figure 15.

    Forces and moments transmitted to C1C2

    FC1+FC2=FeBC1×FC1+BC2×FC2+mC1SC1+mC2SC2=Me
    (76)

    where

    FC1=f1S1+f2S2FC2=f3S3+f4S4
    (77)

    are the internal force vectors acting on the points C1 and C2 mC1 and mC2 are moments transmitted to C1C2 SC1 and SC2 are the unit vectors along the direction of (Sn1×Sn2)×v1 and (Sn3×Sn4)×v2 respectively.

    Introducing expressions (49) into equation (48) becomes as follows:

    i=14fiSiFe=0RABBC1×(i=12Sifi)+RABBC2×(i=34Sifi)+i=12mCiSCiMe=0
    (78)

    Writing equation (50) in a matrix form yields:

    [S1S2S3S400RABl1×S1RABl1×S2RABl2×S3RABl2×S4SC1SC2][f1f2f3f4mC1mC2]=[FeMe]
    (79)

    where l1=BC1,l2=BC2 . The forces and moments at the robot joints due to external load are given by:

    JS1[FeMe]=[f1f2f3f4mC1mC2]
    (80)

    Investigating the transformation matrix JS given in equation (52), one can analyze the singularity of the parallel part. The columns of JS are lines lying on the Klein quadric M24 , as they satisfy Klein’s equation (Hunt, 1978):

    l1l4+l2l5+l3l6=0
    (81)

    where each column of equation (53) is represented by its Plücker line coordinates given by (l1,l2,l3,l4,l5,l6)T

    Moreover, considering the forces and moments exerted on the actuator, one can obtain (see Fig.16).

    τikz+σikpi=fiSi
    (82)
    cosθi=kzSi
    (85)

    where kpi is the unit vector along the direction (kz×ui) . Equation (54.2) can be written for i=1,2,3,4, and the results grouped in a matrix form, τi can be written as:

    JP[f1f2f3f4]=[τ1τ2τ3τ4]
    (86)

    where

    JP=[cosθ10000cosθ20000cosθ30000cosθ4]
    (87)
    media/image264.png

    Figure 16.

    Parallelogram static analysis

    JS (equation (52)) and JP (equation (56)) provide the information about the singularities of the H4 manipulator. The analysis of JS provides the singularity of the parallel part of H4 robot known as the “parallel singularity” and JP provides the singularity due to the parallelogram structure, known as the “serial singularity”. In (Liu et al., 2003), the singularity caused by JP is also called as “actuator singularity”.
    media/image266.png

    Figure 17.

    Serial singularity of H4 manipulator

    Observing JP one can see that singularity occurs whenever θi=π/2+πk (k=1,…,n, i=1,2,3,4). In this case one or more columns of JP are zero and the nacelle can’t move in the Si direction, namely, the manipulator losses a DOF. Fig.17 shows the serial singularity of H4 manipulator.

    4.3. The behavior of the H4 manipulator in parallel singular configurations

    The parallel singular configurations of the H4 are identified by investigating matrix JS . The first four columns of JS define the four internal forces of parallelogram which are along the direction BiAi (see Fig.53 and 55). Utilizing Grassmann geometry, we can obtain all singularities of JS .

    Based on H4’s kinematic architecture, the first four columns of JS can form only the following linear varieties:

    • Type I: a line bundle where all lines intersect a common point, but they are not coplanar. This is the condition 3c descirbed in section 4.1.2.

    • Type II: all the lines belong to the union of two planar pencils of non coplanar lines that have a line in common. This is the condition 3b descirbed in section 4.1.2.

    Type I can be considered as the special condition of Type II. When the first four columns of JS are linear dependent, the actuators can’t control the linear velocity of the nacelle, generally known as “architecture singularity” (Joshi & Tsai, 2002).

    The last two columns of JS define the two moments of constraint which are along the direction SCi . We make here the assumption that the two revolute joints placed on the nacelle have the same direction represented by vector v . This choice may be not the only one, but it worth noting that placing two revolute joints parallel to each other on the nacelle is a good choice for practical matters. The existence of a rotational motion about v can be written as follows (Pierrot & Company, 1999):

    When SC1×SC2=0 , the vectors SC1 and SC2 are linearly dependent. The constraint moments exerted on the nacelle are along the same direction, and the robot would gain additional DOF(s). This singularity arises from deficiency of constraints, so we call this singularity as “constraint singularity”.

    In order to investigate the instantaneous motion the robot tends to perform while in singular configurations, the LCAA which was presented at (Pottmann et al., 1999) and further used for robotic analysis in (Wolf & Shoham, 2003) is adopted. The LCAA is based on determining the linear complex C=(c,c¯) which is the closest one to the given set of lines Li=(li,l¯i) . The moment of Li with respect to a linear complex C is given by:

    m(L,C)=|c¯li+cl¯i|c
    (89)

    Hence, given a set of lines Li , the closest linear complex among all linear complexes χ can be given by the minimization of:

    i=1km(Li,χ)
    (90)

    where χ is given by χ=(x,x¯)6 .

    Equation (59) is equivalent to minimizing the positive semi-definite quadratic form given by (Pottmann et al., 1999):

    F(X)=i=1k(x¯li+xl¯i)2=χTMχ
    (91)

    under the normalization condition 1=x2=χTDχ , where D=diag(1,1,1,0,0,0) χ presents the set of all linear complexes given by χ=(x,x¯)6 , and M is the Gramian matrix of JS . Equation (60) is a general eigenvalue problem and the solution, C , is the eigenvector χi corresponding to the smallest eigenvalue. Moreover, given the closest linear complex found by the algorithm, the axis, A , of the linear complex and its pitch can also be revealed. In fact, the minimization procedure shown in above minimizes the work of the wrench applied on the moving platform and the twists. When the smallest eigenvalue is zero, there is no work generated by the set of wrenches when the platform moves in the twist direction given by χ . So we can use LCAA to investigate the behavior of H4 manipulator in singular configurations.

    4.4. Simulation and results

    In this section, the parallel singularities of an example H4 manipulator are presented and the LCAA method is used to investigate the instantaneous motion the manipulator tends to perform while in these singular configurations. For simulation, the design parameters we have selected is described in Fig.11, where the parameters have been chosen as described in section 3.2.2.

    Architecture singularities

    Architecture singular configurations of the H4 manipulator are identified using line geometry (Merlet, 2006). Note the first four columns of JS as:

    [S1S2S3S4RBDl1×S1RBDl1×S2RBDl2×S3RBDl2×S4]=[a1a2a3a4a¯1a¯2a¯3a¯4]=[A1A2A3A4]
    (92)

    where ai=1 and aia¯i=0 , so the six coordinates with (ai,a¯i) are called the normalized Plücker coordinates of Ai (i=1,2,3,4). Deriving from the geometric conditions (Type I and Type II) in section 4.3, we can find that when the coordinates of point D satisfy the following condition:

    x=9.5812y=9.4130θ=1.5080
    (93)

    singularity occurs. The closest linear complex’s axis found by the algorithm is:

    A=[0015.29370.33310]
    (94)

    The standard deviation of the lines form C is given by: σ=λik5=0 and the linear complex’s pitch is given by: p=cc¯c2=9.4047 Fig.18 demonstrated the result of the LCAA.

    media/image296.png

    Figure 18.

    H4 robot with linear complexes (a): 3D view (b): top view

    Corresponding to the non-zero pitch value of the linear complex, at this singularity the robot performs a screw motion of rotation and translation around the linear complex axis A

    Constraint singularities

    From the discussion in section 4.3, we know that when SC1×SC2=0 the robot is in a constraint singular configuration. If we select the structure parameters as Fig.11, when:

    x=3.6683y=3.3709θ=2.3905
    (95)

    the robot reaches singular configuration.

    Constraint singularity occurs in limited-DOF manipulator when the screw system, formed by the constraint wrenches in all chains, loses rank. In this section, the geometric explanation of constraint singularity is given, then LCAA is applied to investigate the behavior of H4 manipulator in constraint singular configuration. From the viewpoint in (Zlatanov et al., 2002), the system of the constraining wrenches of H4, W, contains as a subspace the 2-system of all pure moments. Then its reciprocal system, the freedoms system T, must contain 3 translations and 1 rotation. In a H4 manipulator, the two “meta-chain” constrain system Wp contains two pure moments, one’s direction is along SC1 and the other’s direction is along SC2 . When the two moments of the two “meta-chain” are linearly independent, two pure moments will be in W.

    However, in the constraint singular configuration, SC1 and SC2 are parallel. The constraining moments of two “meta-chain” are identical and the system Wp are the same one-system. Then, W consists of only one screw. Hence, the twist system of the robot is a five-system and the mechanism gains new additional freedom. Applying the method of the closest linear complex on the H4 robot in its constraint singular configuration, we can find the closest linear complex’s axis is: A=[0.91370.35240.20240.32080.20641.0887] and the linear complex’s pitch is: p=cc¯c2=1.0469 . SC1=SC2=[0.35990.93300]T means that this configuration is singular and the robot is able to execute a screw motion of rotation and translation around the linear complex axis. The direction of the linear complex axis is perpendicular to SC1 (and SC2 ) because aSC1=0 (and aSC2=0 ). This corroborates the results deduced from the viewpoint of constraint singularities. Fig.19 demonstrates the result of the LCAA in constraint singular configurations.

    media/image304.png

    Figure 19.

    H4 robot with linear complexes (a): 3D view (b): top view

    Discussion: Simulation Results of the H4 robot

    From above discussion, only one architecture singular configuration can be found. In this configuration, the intersection point of A1,A2 and the intersection point of A3,A4 (see equation (61)) are superposed, which satisfies the geometric condition of Type I (in section 4.3). However, for constraint singularities, as long as SC1×SC2=0 is satisfied, the manipulator would reach singular configurations. Four parameters x,y,z,θ are used to describe the configuration of a H4 manipulator. z is independent of the singularities. The singular equation SC1×SC2=0 provides only one constraint equation, so the three parameters x,y,θ constitute a three dimension curved surface. All the points in the curved surface satisfy the equation SC1×SC2=0 . We call this surface as constraint singular surface. The equation of the constraint singular surface can be written as:

    (SC1×SC2)v=0
    (96)

    Fig.20 depicts the constraint singular surface.

    media/image313.png

    Figure 20.

    The constraint singular surface

    As we can see from the reference (Wolf & Shoham, 2003):

    F(X)=i=1k(x¯li+xl¯i)2=χTMχ=λ
    (97)

    λ has the meaning of the sum of the square of the virtual work generated by the structure to the nacelle, when the latter is moving instantaneously in the twist direction given by χ . Hence, if λ=0 , the chains would lose the constraints on the nacelle and the nacelle could move instantaneously in the twist direction given by χ

    The standard deviation of the lines Li from C, may be written as:

    σ=λk5=i=1km(Li,c)2k5
    (98)

    which is the square of error from Li to C . In case of a sufficiently small deviation σ , the given lines Li can be well approximated by lines of the linear complex C , which means that the current configuration is close to the singular configuration. When σ=0 , the lines Li lie on the linear complex C , which means that the manipulator is in a singular configuration and the nacelle can move in the twist direction of the linear complex axis.

    In order to investigate the robot’s behavior in the vicinity of the singular configuration, a simulation is performed in the neighborhood of the architecture singular configuration (Fig.21-Fig.22).

    media/image319.png

    Figure 21.

    Minimum σ in the vicinity of the architecture singular configuration

    media/image320.png

    Figure 22.

    Pitch of linear complex with minimum σ in the vicinity of the architecture singular configuration

    All the simulations are in the adjacent range of the singular configuration: x[10.19.1] y[9.98.9] z=0

    It can be seen that the manipulator is no longer in a singular position. However, when manufacturing tolerances are greater than the value σ , there is still possibility of an uncontrolled motion of the nacelle according to the interpretation of σ

    5. Conclusions

    This chapter mainly discusses the architectures and some kinematic problems about the 4-DOF parallel manipulator H4 which is suitable for pick and place operations. The geometric methods used to analyze the singularity poses are introduced in detail. As this type of parallel manipulator generally performs tasks implying significant accuracy, velocities and accelerations, there are still numerous topics need further study.

    It is the complexity of the parallel manipulator’s implantation, of their control and of their design that hinders their widely use. Control efficiency is related to kinematics calculation. Although the analytical formulation of the direct kinematics is deduced in section 3.2.1, solving the 16th degree univariate polynomial equation is not trivial and efficient control still requires very fast algorithms which give all of the solutions.

    The motions of H4 are mainly restricted by mechanical limits on passive joints, interference links and limitations due to the actuators. Because of the coupling between translations and rotations, it is generally complex to determine the workspace of H4. There are various methods can be used to calculate the workspace of H4, such as the discretisation approach, inequality constraints approach (Jo, 1989) and algebraic approach (Husty, 1996). Merlet (Merlet, 2006) has pointed out that the geometrical approaches are very efficient for the determination of various types of workspace and algebraic geometry should also play an important role in some of the algorithms which are used to calculate the workspace of parallel manipulators.

    Singularity analysis of parallel manipulators is crucial for practical applications. Besides the geometric method suggested in section 4, some new ideas are provided to analyze singularities of parallel manipulators by using different mathematic tools. Differential geometric is a powerful tool for general analysis of parallel manipulators. From the geometric properties of the parallel manipulator, Liu et al. presented a framework for analyzing the singularities of a parallel manipulator (Liu et al., 2003). Topological structure of these singularities and their relations with the kinematic parameters of the system were investigated systematically using Morse function theory. An intrinsic definition of parameterization singularities was given in this chapter. Due to the profound explanation of the parallel manipulators’ characteristics, although the differential geometric tool is abstract, its use for general analysis and extension to geometric control of parallel manipulator is still worthy of study. In some situations, it is only necessary to seek singular configurations in the workspace of parallel manipulators. Su et al. (Su et al., 2003) proposed a new singularity analysis method for Steward platform using genetic algorithm (GA). The square of determinant of the Jacobian matrix is selected as the objective function, and the minimal of this objective function is searched in the workspace of Steward platform by the GA. The singularity of Steward platform depends on this minimal objective function: if this value is zero, the singularity of Steward platform will take place, otherwise, the Sterward platform is singularity-free. This method can be used to determine the dimensioning of a H4 so that it is singularity-free, at least in a given workspace. Of course, the effectiveness and speediness of the GA are the chief problems which should be considered.

    The study of the novel parallel manipulator H4 and its various reformative structures has become increasingly important during recent years for its excellent performance of velocity and acceleration. But numerous problems still remain open, especially the positioning accuracy of the end-effector. From the paper written by Renaud et al. (Renaud et al., 2003), the positioning accuracy of a H4 can only reach lower than 0.5mm, which cannot meet the requirements (generally lower than 0.05mm) of the semiconductor end-package equipments. We hope that this chapter will arise some attention to extend this type of parallel manipulators to semiconductor industry.

    6. Acknowledgements

    The work is supported by the Natural Science Fund of China (NSFC) (Project No. 50625516) and the National Fundamental Research Program (973) (Project Nos. 2003CB716207 and 2007CB714000).

    References

    1 - S. K. Agrawal, G. Desmier, S. Li, 1995 Fabrication and analysis of a novel 3 DOF parallel wrist mechanism. ASME Journal of Mechanical Design, 117 June, 343 345
    2 - J. Angeles, A. Morozov, O. Navarro, 2000 A novel manipulator architecture for the production of scara motions, Proceedings of IEEE International Conference on Robotics and Automation, 2370 2375 , San Francisco, USA, April 2000
    3 - J. Angeles, 2005 The degree of freedom of parallel robots: a group-theoretic approach, Proceedings of IEEE International Conference on Robotics and Automation, 1017 1024 , Barcelona, Spain, April 2005
    4 - D. Basu, A. Ghosal, 1997 Singularity analysis of platform-type multi-loop spatial mechanisms. Mechanism and Machine Theory, 32 3 375 389
    5 - I. Benov, 2001 Delta parallel robot- the story of success. May, 6, http://www.parallemic.org/Reviews/Review002.html
    6 - L. Bonev, D. Zlatanov, 2001 The mystery of the singular SUN translational parallel robot. June, 12, 2001, http://www.parallemic.org/Reviews/Review004.html
    7 - I. A. Bonev, et al. 2001 A closed-form solution to the direct kinematics of nearly general parallel manipulators with optimally located three linear extra sensors. IEEE Transaction on Robotics and Automation, 17 2 148 156
    8 - O. Bottema, B. Roth, 1990 Theoretical kinematics. New York: Dover Publications
    9 - T. Brogardh, 2002 2002. PKM research- important issues, as seen from a product development perspective at ABB robotics, Workshop on Fundamental Issues and Future Research Directions for parallel Mechanisms and Manipulators, Quebec, Canada, 2002
    10 - H. Bruyninckx, 1997 The 321-hexa: a fully parallel manipulator with closed-form position and velocity kinematics, Proceedings of IEEE Conference on Robotics and Automation, 2657 2662 , Albuquerque, April, 1997
    11 - J. A. Carretero, R. P. Podhorodeski, M. A. Nahon, C. M. Gosselin, 2000 Kinematic analysis and optimization of a new three degree-of-freedom spatial parallel manipulator. ASME Journal of Mechanical Design, 122 17 24
    12 - D. D. Chablet, P. Wenger, 2002 Design of a three-axis isotropic parallel manipulator for machining applications: the orthoglide, Workshop on Fundamental Issues and Future Research Directions for Parallel Mechanisms and Manipulators, Quebec, Canada, Oct., 2002
    13 - J. W. F. Cheung, Y. S. Hung, G. P. Widdowson, 2002 Design and analysis of a novel 4-dof parallel manipulator for semiconductor applications, Proceedings of the 8th Mechatronics Forum International Conference, 1358 1366 , Twente, Netherlands: Drebbel Institute for Mechatronics
    14 - Y. J. Chiu, M. H. Perng, 2001 Forward kinematics of a general fully parallel manipulator with auxiliary sensors. International Journal of Robotics Research, 20 5 401 414
    15 - H. B. Choi, A. Konno, M. Uchiyama, 2003 Closed-form solutions for the forward kinematics of a 4-DOFs parallel robot H4, Proceedings of IEEE/RSJ International Conference on Intelligent Robots and Systems, 3312 3317 , Las Vegas, USA, October, 2003
    16 - R. Clavel, 1988 DELTA, a fast robot with parallel geometry, Proceedings of the 18th International Symposium on Industrial Robot, 91 100 , April, 1988, Lausanne
    17 - C. L. Collins, G. L. Long, 1995 The singularity analysis of an in-parallel hand controller for force-reflected teleoperation. IEEE Transactions on Robotics and Automation, 11 5 661 669
    18 - O. Company, F. Pierrot, 1999 A new 3T-1R parallel robot, Proceedings of the 9th International Conference on Advanced Robotics, 557 562 , October, 1999, Tokyo, Japan
    19 - O. Company, F. Marquet, F. Pierrot, 2003 A new high-speed 4-DOF parallel robot synthesis and modeling issues. IEEE Transactions on Robotics and Automation, 19 3 411 420
    20 - O. Company, F. Pierrot, J. C. Fouroux, 2005 A method for modeling analytical stiffness of a lower mability parallel manipulator, Proceedings of IEEE International Conference on Robotics and Automation, 18 22 , Barcelona, Spain, April, 2005, Milan, September, 1995
    21 - R. Di Gregorio, V. Parenti-Castelli, 1999 Mobility analysis of the 3-UPU parallel mechanism assembled for a pure translational motion, Proceedings of the IEEE/ASME International Conference on Advanced Intelligent Mechatronics, 520 525 , Atlanta, GA, September, 1999
    22 - A. Fattah, G. Kasaei, 2000 Kinematics and dynamics of a parallel manipulator with a new architecture. Robotica, 18 535 543
    23 - J. C. Faugère, D. Lazard, 1995 The combinatorial classes of parallel manipulators. Mechanism and Machine Theory, 30 6 765 776
    24 - Y. Feng-Cheng, E. J. Haug, 1994 Numerical analysis of the kinematic working capability of mechanism. ASME Journal of Mechanical Design, 116 111 117
    25 - H. Funabashi, Y. Takeda, 1995 Determination of singular points and their vicinity in parallel manipulators based on the transmission index, Proceedings of the Ninth World Congress on the Theory of Machine and Mechanism, 1977 1981 ,
    26 - C. M. Gosselin, J. Angeles, 1989 The optimum kinematic design of a spherical three-degree-of-freedom parallel manipulator. ASME Journal of Mechanisms, Transmission and Automation in Design, 111 2 202 207
    27 - C. Gosselin, J. Angeles, 1990 Singularity analysis of closed-loop kinematic chains. IEEE Transaction on Robotics and Automation, 6 3 281 290
    28 - C. Gosselin, 1990 Dexterity indices for planar and spatial robotic manipulators. Proceedings of the IEEE International Conference on Robotics and Automation, 650 655 , Cincinnati, May, 1990
    29 - C. Gosselin, 1996 Parallel computational algorithms for the kinematics and dynamics of planar and spatial parallel manipulators. ASME Journal of Dynamic Systems, Measurement and Control, 118 1 22 28
    30 - C. M. Gosselin, É. St-Pierre, 1997 Development and experimentation of a fast 3-DOF camera-orienting device. The International Journal of Robotics Research, 16 5 619 629
    31 - V. E. Gough, 1956 Contribution to discussion of papers on research in automobile stability, control and tyre performance, Proceedings of the Automotive Division of the Institution of Mechanical Engineers, Part D, 392 394
    32 - R. D. Gregorio, 2001 A new parallel wrist using only revolute pairs: the 3-RUU wrist. Robotica, 19 305 309
    33 - F. Hao, J. M. Mc Carthy, 1998 Conditions for line-based singularities in spatial platform manipulators. Journal of Robotic Systems, 15 1 43 55
    34 - J. Hesselbach, H. Kerle, 1995 Structurally adapted kinematic algorithms for parallel robots up to six degrees if freedom (dof), Proceedings of the 9th World Congress on the Theory of Machines and Mechanisms, 1935 1939 , Milan, Italy, August, 1995
    35 - K. H. Hunt, 1978 Kinematic geometry of mechanism. Department of Mechanical Engineering, Monash University, Clayton, Victoria, Australia
    36 - M. L. Husty, 1996 On the workspace of planar three-legged platforms. Proceedings of the World Automation Congress, 3 339 344 , Montpellier, May, 1996
    37 - D. Y. Jo, E. J. Haug, 1989 Workspace analysis of closed loop mechanism with unilateral constraints. Proceedings of the ASME Design Automation Conference, 53 60 , Montreal, September, 1989
    38 - S. A. Joshi, L. Tsai, W. , 2002 Jacobian analysis of limited-dof parallel manipulators. ASME Journal of Mechanical Design, 124 2 254 258
    39 - D. Kim, W. Chung, Y. Youm, 2000 Analytic jacobian of in-parallel manipulators. Proceedings of IEEE International Conference on Robotics and Automation, 2376 2381 , San Francisco, April, 2000
    40 - S. G. Kim, J. Ryu, 2003 New dimensionally homogeneous jacobian matrix formulation by three end-effector points for optimal design of parallel manipulators. IEEE Transactions on Robotics and Automation, 19 4 731 736
    41 - S. Krut, O. . Company, M. Benoit, et al. 2003 I4: a new parallel mechanism for scara motions, Proceedings of IEEE International Conference on Robotics and Automation, 1875 1880 , Taipei, Taiwan, September, 2003
    42 - S. Krut, V. Nabat, O. Company, F. Pierrot, 2004 A high-speed parallel robot for scara motions, Proceedings of IEEE International Conference on Robotics and Automation, 4109 4115 , New Orleans, USA, April, 2004
    43 - Q. Li, Z. Huang, 2003 Mobility analysis of lower-mobility parallel manipulators based on screw theory, Proceeding of IEEE International Conference on Robotics and Automation, 1179 1184 , Taipei, Taiwan, September, 2003
    44 - G. F. Liu, Y. J. Lou, Z. X. Li, 2003 Singularities of parallel manipulators: a geometric treatment. IEEE Transactions on Robotics and Automation, 19 4 579 594
    45 - J. P. Merlet, 1992 On the infinitesimal motion of a parallel manipulator in singular configurations, Proceedings of IEEE International Conference on Robotics and Automation, 320 325 , Nice, France, May, 1992
    46 - J. P. Merlet, 2000 Kinematics’ not dead, Proceedings of IEEE International Conference on Robotics & Automation, 1 6 , San Francisco, CA, April, 2000
    47 - J. P. Merlet, D. Daney, 2005 Dimensional synthesis of parallel robots with a guaranteed given accuracy over a specific workspace, Proceedings of IEEE International Conference on Robotics and Automation, 19 22 , Barcelona, Spain, April, 2005
    48 - J. P. Merlet, 2006 Parallel Robots (Second Edition), Springer, Dordrecht, The Netherlands
    49 - V. Nabat, O. Company, M. Krut, et al. 2005 Part 4: very high speed parallel robot for pick-and-place, Proceedings of IEEE International Conference on Intelligent Robots and Systems, 553 558 , Edmonton, Canada, August, 2005
    50 - V. Nabat, O. Company, F. Pierrot, 2006 Dynamic modeling and identification of Part 4, a very high speed parallel manipulator, Proceedings of the International Conference on Intelligent Robots and Systems, 496 501 , Beijing, China, October, 2006
    51 - P. Parenti-Castelli, P. Di Gregorio, 20010 Real-time actual pose determination of the general fully parallel spherical wrist, using only one extra sensor. Journal of Robotic Systems, 18 12 723 729
    52 - G. R. Pennock, D. J. Kassner, 1990 Kinematic analysis of a planar eight-bar linkage: application to a platform-type robot, ASME Proceedings of the 21th Biennial Mechanism Conference, 37 43 , Chicago, September, 1990
    53 - F. Pierrot, O. Company, 1999 H4: a new family of 4-dof parallel robots, Proceedings of the 1999 IEEE/ASME International Conference on Advanced Intelligent Mechatronics, 508 513 , Atlanta, USA, September, 1999
    54 - H. Pottmann, M. Peternell, B. Ravani, 1999 An introduction to line geometry with applications. Compter Aided Design, 31 3 16
    55 - C. Reboulet, A. Robert, 1985 Hybrid control of a manipulator with an active compliant wrist, The 3rd International Symposium of Robotics Research, 76 80 , Gouvieux, France, October, 1985
    56 - P. Renaud, N. Andreff, F. Marquet, et al. 2003 Vision-based kinematic calibration of a H4 parallel mechanism. Proceedings of the IEEE International Conference on Robotics & Automation, 1191 1196 , Taipei, Taiwan, September, 2003
    57 - Robotics Department of LIRMM. H4 with rotary drives: a member of H4 family, dedicated to high speed applications. http://www.lirmm.fr/w3rob/SiteWeb/detail_resultat.php?num_resultat=7&num_topic=1& num_projet=1&num_activite=3
    58 - L. Rolland, 1999 The manta and the kanuk: novel 4-dof parallel mechanisms for industrial handling. ASME Dynamic Systems and Control Division, IMECE’99 Conference, 831 844 , Nashville, USA, November 1999
    59 - L. Romdhane, Z. Affi, M. Fayet, 2002 Design and singularity analysis of a 3-translational-DOF in-parallel manipulator. ASME Journal of Mechanical Design, 124 419 426
    60 - J. Ryu, J. Cha, 2003 Volumetric error analysis and architecture optimization for accuracy of HexaSlide type parallel manipulators. Mechanism and Machine Theory, 38 227 240
    61 - D. Steward, 1965 A platform with 6 degrees of f reedom. Proceedings of the Institution of Mechanical Engineers, 180 Part 1, 15 371 386 , British
    62 - B. M. St-Onge, C. M. Gosselin, 2000 Singularity analysis and representation of the general Gaugh-Steward platform. International Journal of Robot Research, 19 3 271 288
    63 - Y. X. Su, B. Y. Duan, B. Peng, et al. 2003 Singularity analysis of fine-tuning stewart platform for large radio telescope using genetic algorthm. Mechatronics, 13 413 425
    64 - L. W. Tsai, 1996 Kinematics of a three-DOF platform with three extensible limbs. Recent Advances in Robot Kinematics, Kluwer Academic Publishers, 401 410
    65 - L. W. Tsai, S. A. Joshi, 2000 Kinematics and optimization of a spatial 3-UPU parallel manipulators. ASME Journal of Mechanical Design, 122 4 439 446
    66 - A. Wolf, M. Shoham, 2003 Investigation of parallel manipulators using linear complex approximation. ASME Journal of Mechanical Design, 125 3 564 572
    67 - J. Wu, Z. Yin, Y. Xiong, 2006 Singularity analysis of a novel 4-dof parallel manipulator H4. International Journal of Advanced Manufacturing and Technology, 29 794 802
    68 - D. Zlatanov, R. G. Fenton, B. Benhabit, 1994 Singularity analysis of mechanism and robots via a motion-space model of the instantaneous kinematics. Proceedings of IEEE International Conference on Robotics and Automation, 980 991 , San Diego, CA, May, 1994
    69 - D. Zlatanov, I. A. Benov, C. M. Gosselin, 2002 Constraint singularities of parallel mechanisms, Proceedings of the IEEE International Conference on Robotics and Automation, 496 502 , Washington DC, May, 2002
    70 - P. B. Zobel, P. Di Stefano, T. Raparelli, 1996 The design of a 3 dof parallel robot with pneumatic dirves, Proceedings of the 27th International Symposium on Industrial Robot, 707 710 , October, 1996, Milan