Stiffness Model of a 3-DOF Parallel Manipulator with Two ...cdn.intechopen.com/pdfs/47697.pdf · of a 3-DOF parallel manipulator with two additional legs. ... Compared to serial manipulators,

Stiffness Model of a 3-DOF Parallel Manipulator with Two Additional Legs

Regular Paper

Guang Yu1,2, Jun Wu1,2,* and Liping Wang1,2 1 State Key Laboratory of Tribology and Institute of Manufacturing Engineering, Department of Mechanical Engineering, Tsinghua University, Beijing, China 2 Beijing Key Lab of Precision/Ultra-precision Manufacturing Equipment and Control, Beijing, China * Corresponding author E-mail: [email protected] Received 11 Jun 2014; Accepted 30 Aug 2014 DOI: 10.5772/59306 © 2014 The Author(s). Licensee InTech. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract This paper investigates the stiffness modelling of a 3-DOF parallel manipulator with two additional legs. The stiffness model in six directions of the 3-DOF parallel manipulator with two additional legs is derived by performing condensation of DOFs for the joint connection and treatment of the fixed-end connections. Moreover, this modelling method is used to derive the stiffness model of the manipulator with zero/one additional legs. Two performance indices are given to compare the stiffness of the parallel manipulators with two additional legs with those of the manipulators with zero/one additional legs. The method not only can be used to derive the stiffness model of a redundant parallel manipulator, but also to model the stiffness of non-redundant parallel manipulators. Keywords Stiffness, Actuation Redundancy, Parallel Manipulator, Comparison Study

1. Introduction Parallel manipulators have attracted much attention from the academic and industrial communities due to potential

applications not only as robot manipulators but also as machine tools. Compared to serial manipulators, parallel manipulators bring the benefits of much higher payload-mass ratios, higher stiffness and better dynamic performance [1-3]. However, relatively small workspace and the abundance of singularities within the workspace partly nullify the aforementioned advantages. It has been shown by some researchers that redundant actuation is an effective method to tackle these shortcomings [4-5]. Redundant actuation eliminates singularities, increases the usable workspace and partially controls the internal forces [6-8]. Thus, redundantly actuated parallel manipulators are rapidly becoming a focus of research due to the multitude of potential advantages they provide [9-11]. Stiffness performances are of great importance for parallel manipulators, particularly those which are used as machine tools, since higher stiffness allows higher machining speeds and feeds while providing the desired precision, surface finish and tool life. The stiffness of a parallel kinematic machine (PKM) in some configurations in the usable workspace might be lower than desirable. Therefore, it is important to carry out the stiffness modelling and to evaluate the stiffness of a PKM in the

1Guang Yu, Jun Wu and Liping Wang: Stiffness Model of a 3-DOF Parallel Manipulator with Two Additional Legs

ARTICLE

Int J Adv Robot Syst, 2014, 11:173 | doi: 10.5772/59306

International Journal of Advanced Robotic Systems

usable workspace. Due to the complex geometry together with the changing rigidity throughout the workspace, it is a complex task to achieve the stiffness of a parallel manipulator. It can be accomplished with the aid of finite element analysis (FEA) software such as ANSYS, ABAQUS and NASTRAN [12]. For a parallel manipulator with fixed length and linear drives, the FEA model has to be remeshed over and over again to make the nodes on the machine frame match with those on the sliders, which vary with the changing configurations of the manipulator. This results in a very tedious and time-consuming routine [13]. An accurate stiffness matrix of a parallel manipulator in stiffness modelling is critical for stiffness evaluation. Some publications have concentrated on the stiffness modelling of parallel manipulators [14-16]. Gosselin [17] and Svinin [18] investigated the stiffness of the Stewart platform and the mapping between the driving force and the platform deformation based on the Jacobian matrix. Li et al. [19] derived intuitively the stiffness matrix of a 3-PUU PKM based on an alternative approach considering actuations and constraints. However, there is little work on the stiffness modelling of redundant parallel manipulators and even less work on the stiffness comparison of redundant parallel manipulators with their non-redundant counterparts. In this paper, a stiffness modelling method is proposed to derive the stiffness of a parallel manipulator in six directions. Two performance indices are given to compare the stiffness of the manipulator without additional legs with the stiffnesses with one/two additional legs. 2. Structure description and kinematic analysis 2.1. Structure description The parallel manipulator studied here, with two translational DOFs and one rotational DOF, is shown in Figure 1. It is composed of sliders 1 3A A and 2 4A A , the moving platform 1 3B B , constant-length links 1 1A B and

2 2A B and extendable links 3 3A B , 4 4A B and 5 5A B . One end of the links 1 1A B and 2 2A B is fixed to sliders

1 3A A and 2 4A A , and the other two are fixed to the moving platform 1 3B B , respectively. When sliders 1 3A A and 2 4A A slide along the vertical guide rails, the

moving platform can move in Y and Z directions. Extendable links 3 3A B and 4 4A B are driven by inner active actuators, with one end fixed to the sliders 1 3A A ,

2 4A A and the other two fixed to the moving platform

1 3B B respectively. Extendable link 5 5A B is also driven by an inner active actuator, with one end fixed on the beam and the other to the moving platform 1 3B B . It can be concluded that the manipulator has a 3-DOF motional

capability even if links 4 4A B and 5 5A B are removed. In this paper, parallel manipulators with zero/one legs are as shown in Figure 1: without legs 4 4A B and 5 5A B , or without leg 5 5A B , respectively.

A1

A3

B1(B2)

A5

A2

A4

C1 C2

l3

e

l1

l2

l4

B3(B4B5)

l5

O

Z

Y

Y'

Z'

O'

d

a

β1

β2

γ

α1

α2

h

Figure 1. Kinematic model of the parallel manipulator 2.2. Kinematics analysis As shown in Figure 1, two coordinate systems O XYZ− and O X Y Z′ ′ ′ ′− are created. The global coordinate system, O XYZ− , is attached to the middle of the frame bottom, and the moving coordinate system, O X Y Z′ ′ ′ ′− , is positioned at point 1 2( )B B of the moving platform,

where its Z′ -axis is along the direction of the constant-length link 1 3B B . The angle a is the orientation (pose) of the moving platform. O[ ]i ′A and

O[ ]i ′B are defined as the position vectors of the points

iA and iB with respect to O X Y Z′ ′ ′ ′− , respectively. The position vector of the origin point O′ with respect to O XYZ− is [ , ]Ty z=r . Let the position vector of the point iA with respect to O XYZ− be [ ]i OA . The position vector of the point iB with respect to O XYZ− can be expressed as

[ ] [ ]i O i O′= +B R B r ( 1,2,3,4,5)i = (1)

where R is the rotation matrix of the platform with

respect to O XYZ− , and cos sinsin cos

a aR

a a−

=

. Then, the

constrained equation of the kinematic chain can be written as

[ ] [ ]i O i O il− =A B ( 1,2,3,4,5)i = (2)

2 Int J Adv Robot Syst, 2014, 11:173 | doi: 10.5772/59306

where 1l and 2l are the lengths of two constant-length links 1 1A B and 2 2A B , and 3l , 4l and 5l are the lengths of the extendable links 3 3A B , 4 4A B and 5 5A B . The inverse kinematic solutions of the parallel manipulator can be written as

2 21

12 2

222 2

3 1

2 242

5 2 2

( / 2)

( / 2 )

( sin / 2) ( cos )

( sin / 2) ( cos )

( sin ) ( cos )

z e l y dzz e l d yz

l y e a d z e a zl y e a d z e a zl

y e a z e a h

+ + − + + + − − = − + + + − − − + + − − + + −

(3)

where e is the length of sliders 1 3A A and 2 4A A and the moving platform 1 3B B , h is the height of the column, d is the width between two columns, and 1z and 2z are the Z-directional position of points 3A and

4A on the sliders, respectively. 3. Stiffness modelling In stiffness modelling, the parallel manipulator is simplified as many spatial beam elements. Based on structural mechanics, the element stiffness matrix in the element coordinate system is built. Then, constraint processing according to the connection forms between two elements is performed for the element stiffness matrix. Finally, the global stiffness matrix is obtained by using the element integration method to integrate all the element stiffness matrices. 3.1 Element stiffness matrix Let the end nodes of each element be m and n , the displacements of each node in the element coordinate system be ( ) [ , , , , , ]e T

m m m m m m mx y z α β γ=δ and ( ) [ , , , , , ]e Tn n n n n n nx y z α β γ=δ , and the force vectors of each

node be ( ) [ , , , , , ]e Tm m m m xm ym zmX Y Z M M M=F and

( ) [ , , , , , ]e Tn n n n xn yn znX Y Z M M M=F [20]. Based on the

relationship between the force and the corresponding displacement, the stiffness equation of element ( e ) can be expressed as

( )( )

( )

emem

en n

=

δFk

F δ (4)

where ( ) 12 12e ×∈k R is the element stiffness matrix.

In the global coordinate system, the displacements of nodes m and n are ( ) [ , , , , , ]e T

m m m m m m mx y z α β γ=δ and ( ) [ , , , , , ]e Tn n n n n n nx y z α β γ=δ , and the force vectors of each

node are ( ) [ , , , , , ]e Tm m m m xm ym zmX Y Z M M M=F and

( ) [ , , , , , ]e Tn n n n xn yn znX Y Z M M M=F .

The displacements and force vectors of each node in the global coordinate system can be expressed as

( )( )

( ) ( )

eemTm

e en n

=

δδT

δ δ(5)

( ) ( )

( ) ( )

e eTm m

e en n

=

F FT

F F(6)

where T is the coordinate transformation matrix from the element coordinate system to the global coordinate system, and 12 12×∈T R . By substituting Eqs. (5) and (6) into Eq. (4), the relationship between node displacements and force vectors in the global coordinate system can be written as

( ) ( )( )

( ) ( )

e eem m

e en n

=

F δk

F δ (7)

where ( )ek is the element stiffness matrix in the global coordinate system and ( ) ( )e T e=k T k T . 3.2 DOF condensation Forms of connection between the elements can be divided into three types: rigid connection (e.g., the connection between the beam and the column), joint connection and insertion end. Elements with rigid connection can be integrated directly with the element matrices by the integration method. Elements with joint connection require treatment of DOF condensation. Elements with fixed ends require treatment for the insertion process. Before using the element integration method to obtain the stiffness model of the parallel manipulator, treatment of DOF condensation and insertion process should be performed for the element stiffness matrix, except for the rigid connection elements. In the element coordinate system, the element stiffness matrix is divided into blocks to describe the process of DOF condensation. The stiffness equation of element ( e ) can be expressed as:

( ) ( ) ( )e e ess sn st s s

ns nn nt n n

tts tn tt t

=

k k k δ Fk k k δ F

Fk k k δ (8)

where ssk , snk , stk , nsk , ntk , nnk , tsk , tnk and ttk are

the corresponding block-stiffness matrices. nδ and nF are the displacement and the corresponding nodal force of the DOFs that need to be condensed. sδ and tδ are

the displacements and sF and tF are the corresponding nodal forces of the DOFs that retain and


apply to the element integration method for each element. Based on Eq. (8), nδ can be expressed as

1( )n nn n ns s nt t−= − −δ k F k δ k δ (9)

By substituting Eq. (9) into Eq. (8), the element stiffness matrix can be written as

1 1 * 1

1 1 * 1

00 0 0 0 0

0

ss sn nn ns st sn nn nt s s sn nn n

ts tn nn ns tt tn nn nt t t tn nn n

− − −

− − −

− − −

= − − −

k k k k k k k k δ F k k F

k k k k k k k k δ F k k F

(10)

3.3 Insertion process The lower end of the column is fixed to the base and cannot move in any direction. Thus, fixed-end treatment should be carried out with the insertion process before integrating all the element stiffness matrices to obtain the global stiffness matrix. Since node n is fixed, =0nδ . In the element coordinate

system, the total displacement ( )eδ can be expressed as

( ) ( )( )

( )

e eme m

en

= =

δ δδ0δ

(11)

The element stiffness matrix in the element coordinate system can be expressed as

( )

11 12( )

21 22

e

e

=

k kk

k k(12)

Since =0nδ , considering the relationship between the displacements and force vectors the elements in the element stiffness matrix can be written as

12 21 22= = =k k k 0 . 3.4 Integration method All element stiffness matrices in the element coordinate system can be transferred to the global coordinate system. Then, the global stiffness matrix can be obtained by integrating all the element stiffness matrices in the global coordinate system. The following equation can be obtained:

( )

( )

e

e=K k (13)

where K is the global stiffness matrix of the parallel manipulator. The stiffness model of the parallel manipulator in the global coordinate system can be written as

=F KΔ (14)

where F is the applied force vector of all nodes, and Δ is the displacement vector of all nodes. The displacement vector of all nodes is derived as

-1=Δ K F (15)

Let the stiffness values of the i th stressed point in each direction be xik , yik , zik , mxik , myik and mzik . Applying the

force xif in the X direction of the i th stressed point, the applied force vector of all nodes can be expressed as

6( 1)

0,0, ,0, ,0,0, ,0xi

i

f−

=

F (16)

Therefore, the stiffness of the i th stressed point in X direction can be expressed as

6( 1) 1

xixi

i

fkδ − +

= 1,2,i = (17)

where 6( 1) 1iδ − + is the displacement of the i th stressed

point in X direction. Similar methods can be used to obtain yik , zik , mxik , myik

and mzik . 4. Numerical simulation 4.1 Stiffness distribution The proposed modelling method is used to derive the stiffness model of the 3-DOF parallel manipulator. The parallel manipulator is simplified as a structural system composed of a plurality of spatial beam elements, as shown in Figure 2. The 3-DOF parallel manipulator includes 14 nodes and 17 elements.

5

1

42

8

7

109

3

12

13 14

6

11

Figure 2. Nodes and elements of the 3-DOF parallel manipulator


The elastic modular is 210GPa and the Poisson’s ratio is 0.3. The geometrical parameters of the 3-DOF parallel manipulator are 1 2 1150 ,= =l l mm 1170 ,=d mm

3160h mm= and 250e mm= . The cross-sections of the constant-length link, beam and columns are rectangles and the extendable links are annular, as shown in Figures 3 (a) and 3 (b), respectively. The cross-sectional parameters of constant-length link, beam, columns and moving platform are shown in Table 1. The cross-sectional parameters are 53oR mm= , 40iR mm= for the upper part of the extendable link, and 40oR mm= and 0iR mm= for the lower part [21].

(a) Constant-length link, beams and columns

(b) Extendable link

Figure 3. Shapes of cross-section

Parameters Ld Wd Lh WhBeams and columns 520 560 18 18

Constant-length links

550 114 5.3 5.3

Moving platform 120 25 0 0

Table 1. Rectangular cross-sectional parameters (mm) In order to verify the accuracy of the stiffness model, finite-element simulation software (Ansys APDL) is also used to find the stiffness of the parallel manipulator. The displacement of the moving platform can be obtained by APDL with the force imposing on the moving platform, and the stiffness of parallel manipulator can be easily obtained by the relationship between the displacement and the force of the moving platform. The stiffness of the 3-DOF parallel manipulator in the workspace with

[ 400 ,400 ]y mm mm∈ − , [700 ,1150 ]z mm mm∈ , 0oa = is simulated by theoretical method and finite-element simulation, respectively.

Figure 4 shows the stiffness distribution with two additional legs obtained by using the method proposed in this paper and using the finite-element method. For both the method proposed in this paper and the finite-element method, the Y-directional and Z-directional position stiffness distributions are symmetrical with respect to the Z axis, due to the symmetrical architecture of the parallel manipulator with two additional legs. When the Z coordinate increases, the stiffness in each direction, except for the Z-directional position stiffness, goes down; Z-directional position stiffness increases due to the effect of an additional leg

5 5A B near the centre of symmetry.

-0.4-0.2

00.2

0.4

0.8

1

1.2

0.98

1

1.02

1.04

x 108

y(m)z(m)

X-d

irect

iona

l pos

ition

stif

fnes

s(N

/m)

simulation resulttheoretical calculations

-0.4-0.2

00.2

0.4

0.8

1

1.2

2.5

3

3.5

4

x 108

y(m)z(m)

Y-d

irect

iona

l pos

ition

stif

fnes

s(N

/m)


-0.4-0.2

00.2

0.4

0.8

1

1.2

1.3

1.4

1.5

1.6

1.7

x 109

y(m)z(m)

Z-di

rect

iona

l pos

ition

stif

fnes

s(N

/m)



-0.4-0.2

00.2

0.4

0.8

1

6.5

7

7.5

8

x 106

y(m)z(m)

X-d

irect

iona

l rot

atio

nal s

tiffn

ess(

N*m

/rad) simulation result

theoretical calculations

-0.4-0.2

00.2

0.4

0.8

1

1.2

4.7

4.8

4.9

5

5.1

x 107

y(m)z(m)

Y-d

irect

iona

l rot

atio

nal s

tiffn

ess(

N*m



-0.4-0.2

00.2

0.4

0.8

1

1.2

7

7.2

7.4

7.6

x 107

y(m)z(m)

Z-d

irect

iona

l rot

atio

nal s

tiffn

ess(

N*m



Figure 4. Comparison of theoretical calculations and simulation results In order to validate the approach proposed in this paper with a finite-element approach, the mean stiffness values in each direction and the error rate between the values obtained by using the two approaches are computed and shown in Table 2. From Table 2, one may see that the theoretical values of the proposed method are consistent with the simulation values. The stiffnesses in each direction are also in close agreement. The overall trend is similar with very small fluctuation.

Mean stiffness

Theoretical values 106(N/m, N*m/rad)

Simulation values 106(N/m, N*m/rad)

Error rate (%)

X position 100.21 102.29 2.03 Y position 278.75 289.10 3.57 Z position 1579.4 1524.4 3.68 X rotation 7.4540 6.8842 7.6 Y rotation 47.769 49.056 2.62 Z rotation 71.952 75.107 4.2

Table 2. Theoretical values and simulation values 4.2 Stiffness comparison In order to evaluate the stiffness performance of the manipulator with two additional legs, the stiffness of the manipulators with zero/one additional legs is also investigated in the workspace, as shown in Figure 5. Compared with those of the manipulators with zero/one additional legs, the X-directional and Z-directional position stiffnesses and the Z-directional rotational stiffness of the manipulator with two additional legs are much higher. The distribution of the X-directional rotational stiffness of the manipulator without additional legs is asymmetrical in the workspace. However, the distribution of the X-directional rotational stiffness of the manipulator with one/two additional legs is symmetrical due to the symmetrical architecture. The manipulator without additional legs has the lowest stiffness among the three parallel manipulators in each direction. To explicitly compare the stiffness of the 3-DOF parallel manipulator with zero/one/two additional legs, two evaluating indices are defined: the mean stiffness in the workspace (denoted as 1kh ) and the ratio of the difference between the maximum stiffness and the minimum stiffness to the maximum stiffness (denoted as

2kh ). 2kh can be expressed as

max min2

maxk

k khk

−= (18)

where maxk and mink are the maximum stiffness and the minimum stiffness in the workspace, respectively[22]. Table 3 shows 1kh and 2kh of the planar parallel manipulators with zero/one/two additional legs with

[ 400 ,400 ]y mm mm∈ − , [700 ,1150 ]z mm mm∈ , and 0oa = . From Table 3, one may see that the Z-directional average position stiffness of the manipulator with zero/one/two additional legs is higher than the X-directional and Y-directional average position stiffness, and the X-directional average position stiffness is the lowest. The ratio of the difference according to the maximum stiffness and the minimum stiffness in the Y-directional and Z-directional position stiffness is greater than that in the X-directional position stiffness.


-0.4-0.2

00.2

0.4

0.8

1

1.2

9.2

9.4

9.6

9.8

10

10.2

x 107

y(m)z(m)

X-d

irect

iona

l pos

ition

stif

fnes

s(N

/m)

without additional legswith one additional legwith two additional legs

-0.4-0.2

00.2

0.4

0.8

1

2.5

3

3.5

x 108

y(m)z(m)

Y-d

irect

iona

l pos

ition

stif

fnes

s(N

/m)


-0.4-0.2

00.2

0.4

0.8

1

1.2

1

1.2

1.4

1.6

x 109

y(m)z(m)

Z-di

rect

iona

l pos

ition

stif

fnes

s(N

/m)


-0.4-0.2

00.2

0.4

0.8

1

2

4

6

8

x 106

y(m)z(m)

X-d

irect

iona

l rot

atio

nal s

tiffn

ess(

N*m

/rad

)


-0.4-0.2

00.2

0.4

0.8

1

1.2

4.65

4.7

4.75

4.8

4.85

4.9

4.95

x 107

y(m)z(m)

Y-d

irect

iona

l rot

atio

nal s

tiffn

ess(

N*m

/rad

)


-0.4-0.2

00.2

0.4

0.8

1

1.26.4

6.6

6.8

7

7.2

7.4

x 107

y(m)z(m)

Z-di

rect

iona

l rot

atio

nal s

tiffn

ess(

N*m

/rad

)


Figure 5. Comparison result with zero/one/two additional legs

Stiffness

Without additional legs

With one additional leg

With two additional legs

hk1

(N/μm) hk2

(%)hk1

(N/μm) hk2

(%) hk1

(N/μm)hk2

(%) X position 94.5 6.6 94.6 6.5 100.2 6.5 Y position 275.3 41.7 278.3 42.9 278.7 43.1 Z position 941.8 23.5 1342.3 23.5 1579.4 17.3 X rotation 4.1 93.6 7.3 18.3 7.4 19.8 Y rotation 47.3 6.4 47.5 6.4 47.7 7.1 Z rotation 66.7 8.0 67.0 8.0 71.9 6.9

Table 3. 1kh and 2kh of planar parallel manipulators

According to the value of hk1 from Table 3, the mean rotational stiffness in the X direction is the lowest among the X-directional, Y-directional and Z–directional rotational and position stiffnesses. Moreover, the ratio of the difference according to the maximum stiffness and the minimum stiffness in the X-directional rotational stiffness is the greatest. Thus, the X-directional rotation stiffness of the 3-DOF parallel manipulator with zero/one/two additional legs reflects the stiffness performance of the parallel manipulator, since the stiffness of a manipulator is determined by the lowest stiffness. Compared with the parallel manipulator without additional legs, the X-directional rotation stiffness of the 3-DOF parallel


manipulator with one additional leg approximately doubled due to the one additional leg; the X-directional rotation stiffnesses of the 3-DOF parallel manipulator with one and two additional legs were similar. This means that the stiffness performances of the parallel manipulator with one/two additional legs are better than the manipulator without additional legs. 5. Conclusions The stiffness model of a 3-DOF parallel manipulator with two additional legs is investigated. The condensation of degrees of freedom for the joint connection and the treatment of fixed ends for the fixed-end connections are taken into account. The stiffness of the manipulator with two additional legs is compared with that of the manipulator with zero/one additional legs. The X-directional and Z-directional position stiffnesses and Z-directional rotational stiffnesses of the manipulator with two additional legs are much higher than those of the manipulator with zero/one additional legs. The manipulator without any additional leg has the lowest stiffness among the three parallel manipulators. This modelling method can derive six-directional stiffness of a parallel manipulator. Moreover, the method is suitable for both redundantly actuated parallel manipulators and non-redundant parallel manipulators. 6. Acknowledgements This work is supported by the National Natural Science Foundation of China (Grant No. 51105225, 51225503, 51375210), as well as a national grant for an excellent doctoral dissertation of PR China (201137), the National Science and Technology Specific Projects (2014ZX04002051, 2013ZX04000030), and the Fund of State Key Laboratory of Tribology (No. SKLT11C01). 7. References [1] Wu H.P (2008) Parallel manipulators towards new

applications. Vienna, Austria: I-Tech Education and Publishing.

[2] Tian Y.L, Shirinzadeh B, Zhang D, Zhong Y (2011) Modeling and Analysis of a three-revolute-revolute- revolute parallel micro-positioning mechanism. Journal of Mechanical Engineering Science. IMechE Part C. 255(5): 1027-1043.

[3] Zhang D, Lei J (2011) Kinematic analysis of a novel 3-DOF actuation redundant parallel manipulator using artificial intelligence approach. Robotics and Computer-Integrated Manufacturing. 27(1): 157-163.

[4] Zhao Y, Gao F (2009) Dynamic performance comparison of the 8PSS redundant parallel manipulator and its non-redundant counterpart- the 6PSS parallel manipulator. Mechanism and Machine Theory. 44(5): 991-1008.

[5] Muller A, Hufnagel T (2012) Model-based control of redundantly actuated parallel manipulators in redundant coordinates. Robotics and Autonomous Systems. 60(1) 563-571.

[6] Wu J, Wang J.S, Wang L.P, Li T.M (2009) Dynamics and control of a planar 3-DOF parallel manipulator with actuation redundancy. Mechanism and Machine Theory. 44(4): 835-849.

[7] Wu J, Li T.M, Wang J.S, Wang L.P (2013) Performance analysis and comparison of planar 3-DOF parallel manipulators with one and two additional branches. Journal of Intelligent and Robotic Systems. 72(1): 73-82.

[8] Wu J, Chen X.M, Li T.M, Wang L.P (2013) Optimal design of a 2-DOF parallel manipulator with actuation redundancy considering kinematics and natural frequency. Robotics and Computer Integrated Manufacturing. 29(1): 80-85.

[9] Mohamed M.G, Gosselin C.M (2005) Design and analysis of kinematically redundant parallel manipulators with configurable platforms. IEEE Transactions on Robotics. 21(3): 277-287.

[10] Kim J, Park F.C, Ryu S.J, Kim J, Hwang J.C, Park C, Iurascu C.C (2001) Design and analysis of a redundantly actuated parallel mechanism for rapid machining. IEEE Transactions on Robotics and Automation. 17(4): 423-434.

[11] Nokleby S.B, Fisher R, Podhorodeski R.P, Firmani F (2005) Force capabilities of redundantly-actuated parallel manipulators. Mechanism and Machine Theory. 40(5): 578-599.

[12] Huang D.T.-Y, Lee J.-J (2001) On obtaining machine tool stiffness by CAE techniques. International Journal of Machine Tools and Manufacture. 41(8): 1149-1163.

[13] Huang T, Zhao X, Whitehouse D.J (2002) Stiffness estimation of a tripod-based parallel kinematic machine. IEEE Transactions on Robotics and Automation. 18(1): 50-58.

[14] Lu Y, Yu J.J, Chen L.W, Hu B, Zhang X.L, Han J.D, Sui C.P (2014) Stiffness and elastic deformation of a 3-leg 5-dof parallel manipulator with one composite leg. International Journal of Robotics and Automation. 29(1): 23-31.

[15] Hu B, Mao B, Yu J, Lu Y (2014) Unified stiffness model of lower mobility parallel manipulators with linear active legs. International Journal of Robotics and Automation. 29(1): 58-66.

[16] Rezaei A, Akbarzadeh A (2013) Position and stiffness analysis of a new asymmetric 2 PRR-PPR parallel CNC machine. Advanced Robotics. 27(2): 133-145.

[17] Gosselin C.M (1990) Stiffness mapping for parallel manipulators, IEEE Transactions on Robotics and Automation. 6(3): 377-382.

[18] Svinin M.M, Hosoe S, Uchiyana M (2001) On the stiffness and stability of Gough-Stewart platforms. In:


Proceedings of the 2001 IEEE International Conference on Robotics and Automation, May 21-26, 2001. Seoul, Korea. IEEE. pp. 3268-3273.

[19] Li Y.M, Xu Q.S (2008) Stiffness analysis for a 3-PUU parallel kinematic machine. Mechanism and Machine Theory. 43(00): 186-200.

[20] Zhang H (2004) Analysis and experiments of stiffness and dynamics for hybrid machine tool. Beijing: Tsinghua University.

[21] Zhang H, Li Y.W, Wang L.P, Tang X.Q (2004) Static stiffness analysis of the gantry hybrid machine tool, Journal of Tsinghua University (Science and Technology). 02: 182-185.

[22] Chen J, Wang L.P, Zhang H, Tang X.Q (2005) Static stiffness analysis and design of a 2-DoF Planar parallel manipulator. Chinese Journal of Mechanical Engineering. 07: 158-163.


Stiffness Model of a 3-DOF Parallel Manipulator with Two ...cdn.intechopen.com/pdfs/47697.pdf · of a 3-DOF parallel manipulator with two additional legs. ... Compared to serial manipulators,

Documents