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, New Developments", book edited by Jee-Hwan Ryu, ISBN 978-3-902613-20-2, Published: April 1, 2008 under CC BY-NC-SA 3.0 license. © The Author(s).

Chapter 21

Dynamic Modelling and Vibration Control of a Planar Parallel Manipulator with Structurally Flexible Linkages

By Bongsoo Kang and James K. Mills
DOI: 10.5772/5380

Article top

Overview

Configuration of a PRR type planar parallel manipulator.
Figure 1. Configuration of a PRR type planar parallel manipulator.
Coordinate system of rigid-body model.
Figure 2. Coordinate system of rigid-body model.
Schematic diagram of loop close equation.
Figure 3. Schematic diagram of loop close equation.
Platform coordinates.
Figure 4. Platform coordinates.
Coordinate system of structurally-flexible model.
Figure 5. Coordinate system of structurally-flexible model.
Intermediate link with PVDF layer.
Figure 6. Intermediate link with PVDF layer.
Intermediate link with PZT actuator.
Figure 7. Intermediate link with PZT actuator.
Error profile of the platform (dotted: no damping, solid: with PVDF layer).
Figure 8. Error profile of the platform (dotted: no damping, solid: with PVDF layer).
Flexible deformation of each link (dotted: no damping, solid: with PVDF layer).
Figure 9. Flexible deformation of each link (dotted: no damping, solid: with PVDF layer).
The first three vibration modes of the first link (dotted: no damping, solid: with PVDF layer).
Figure 10. The first three vibration modes of the first link (dotted: no damping, solid: with PVDF layer).
Control output for the first link.
Figure 11. Control output for the first link.
Error profile of platform (dotted: no damping, solid: with PZT actuator).
Figure 12. Error profile of platform (dotted: no damping, solid: with PZT actuator).
Flexible deformation of each link (dotted: no damping, solid: with PZT actuator).
Figure 13. Flexible deformation of each link (dotted: no damping, solid: with PZT actuator).
The first three vibration modes of the first link(dotted: no damping, solid: with PZT actuator).
Figure 14. The first three vibration modes of the first link(dotted: no damping, solid: with PZT actuator).
Control output for the first link.
Figure 15. Control output for the first link.
The first three vibration modes of the first link with PZT actuator located on inappropriate place.
Figure 16. The first three vibration modes of the first link with PZT actuator located on inappropriate place.

Dynamic Modelling and Vibration Control of a Planar Parallel Manipulator with Structurally Flexible Linkages

Bongsoo Kang1 and James K. Mills2

1. Introduction

A parallel manipulator provides an alternative design to serial manipulators, and can be found in many applications such as mining machines (Arai et al., 1991). Through the design of active joints such that actuators are fixed to the manipulator base, the mass of moving components of the parallel manipulator is greatly reduced, and high speed and high acceleration performance may be achieved. Parallel manipulators, comprised of closed-loop chains due to multiple linkages of the parallel structure, also provide high mechanical rigidity, but adversely exhibit smaller workspace and associated singularities. Considerable research has focused on kinematic analysis and singularity characterization of these devices (Gosselin & Angeles, 1990, Merlet, 1996). Planar parallel manipulators typically consist of three closed chains and a moving platform. According to the arrangement of their joints in a chain, these mechanisms are classified as PRR, RRR etc. where P denotes a prismatic joint and R denotes a revolute joint respectively.

An assembly industry, such as the electronic fabrication, demands high-speed, high acceleration placement manipulators, with corresponding lightweight linkages, hence these linkages deform under high inertia forces leading to unwanted vibrations. Moreover, such multiple flexible linkages of a parallel manipulator propagate their oscillatory motions to the moving platform where a working gripper is located. Therefore, such vibration must be damped quickly to reduce settling time of the manipulator platform position and orientation. A number of approaches to develop the dynamic model of parallel manipulators with structural flexibility have been presented in the literature(Fattah et al., 1995, Toyama et al. 2001), but relatively few works related to vibration reduction of a parallel manipulator have been published. Kozak (Kozak et al., 2004) linearized the dynamic equations of a two-degree-of-freedom parallel manipulator locally, and applied an input shaping technique to reduce residual vibrations through modification of the reference command given to the system. Kang (Kang et al., 2002) modeled a planar parallel manipulator using the assumed modes method, and presented a two-time scale controller for linkage vibration attenuation of the planar parallel manipulator. Since both the input shaping technique and the two-time scale control scheme, applied to parallel manipulators, can only control command inputs to joint actuators of a manipulator, their performance on vibration reduction of flexible linkage are limited.

This chapter introduces a methodology for the dynamic analysis of a planar parallel mechanism beginning with the rigid-body model of a planar parallel manipulator. Flexible deformations of each linkage are expressed by the product of time-dependant functions and position-dependant functions, i.e. an assumed modes model. Overall dynamic equations of the motion for a planar parallel manipulator are formulated by Lagrangian equations. Then, an active damping approach using piezoelectric material actuators is presented to damp out oscillation of linkages of a planar parallel manipulator. Also an integrated control scheme is designed to permit the platform of the parallel manipulator to follow a given trajectory while simultaneously damping structural vibration of flexible linkages. Attached directly to the surface of flexible linkages, piezoelectric materials deform under an applied control voltage, producing shear forces, which counteract shear stresses that occur due to deformation of the linkages. Transducers are often developed from either of two types of material: polyvinylidene fluoride (PVDF) or lead zirconium titanate (PZT). PVDF is lightweight and is mainly used as a measurement device for detecting vibration, although the PVDF can also be used as an actuator. PZT has been used as actuators for micro mechanisms and has a higher strain constant than PVDF.

In simulations, both the PZT and PVDF piezoelectric actuators are applied to a planar parallel manipulator with flexible linkages and the respective performance of two actuators are compared. The dynamics of the planar parallel platform are selected such that the linkages have considerable flexibility, to better exhibit the effects of the vibration damping control system proposed. Simulation results show that the PZT actuator can give better performance in vibration attenuation than the PVDF layer.

2. Rigid-body analysis of a planar parallel manipulator

2.1. Architecture of a planar parallel manipulator

media/image1.jpeg

Figure 1.

Configuration of a PRR type planar parallel manipulator.

The architecture of a PRR type parallel mechanism is illustrated in Fig. 1. The underline of the first character of PRR means that the first Prismatic joint is an active joint derived by actuators. The moving platform, a regular triangular shape, exhibits translation and rotational motion in a plane. Three intermediate links between the moving platform and sliders play a role to convert the actuating force into movements of the platform. Both ends of the intermediate link are composed of non-actuated revolute joints. The sliders move along the linear guide and their motions can be achieved by a ball-screw mechanism. The proposed planar manipulator is categorized as a PRR type, because a closed-loop chain consists of a prismatic joint and two consecutive revolute joints. In contrast to well-known RPR type parallel manipulators, actuators of the proposed PRR configuration remain stationary that results in low inertia of moving parts. Workspace analysis of a planar parallel manipulator has been addressed (Gosselin et al., 1996, Heerah et al., 2002) and singularity analysis of a planar parallel manipulator has been studied (Gosselin & Angeles, 1990, Merlet, 1996).

2.2. Kinematics

media/image2.jpeg

Figure 2.

Coordinate system of rigid-body model.

We will begin with formulation of the rigid-body model of the proposed planar parallel manipulator, and then include structural flexibility of the linkage in Section 3. Prior to derivation of dynamics of the parallel manipulator, the inverse kinematic solution of the manipulator is formulated to define kinematic relations between active joints, and position and orientation of the platform. Then, based on these kinematic relations, equations of the motion for a parallel manipulator are formulated later.

Generalized coordinates for the PRR manipulator are defined, as shown in Fig. 2. The position of the reference X-Y frame is arbitrary. Ai is the origin of the i th linear guide and Bi is the position of the i th slider. Ci is the position of the revolute joint of the platform facing with the i th linkage and P is the position of the platform at its mass center. Three linkages, including associated coordinates, are numbered with a subscript, i, starting from the lower right link in a counterclockwise direction. The pose of the moving platform at its mass center can be written with respect to the reference X-Y frame as

X¯P:=[xPyPφ]T
(1)

The displacement of the sliders from their origin, Ai , to B i , are expressed as

ρ¯:=[ρ1ρ2ρ3]
(2)
βi is defined as the angle at B i between the X-axis of the fixed frame and the i th intermediate link and αi is the angle at A i between the X-axis of the fixed frame and the i th linear guide. The position and orientation of the i th link at its mass center can also be written as
X¯i:=[xiyiβi]T
(3)

For each chain, a loop close equation can be written using position vectors defined as

A¯iP¯+P¯C¯i=A¯iB¯i+B¯iC¯ii     =1,2,3
(4)

Fig. 3. shows the diagram of the loop close equation when the first link is considered.

The right-hand side of equation (4), the coordinates of C i , is written as

xci=xai+ρicosαi+lcosβi
(5)
yci=yai+ρisinαi+lsinβi
(6)

where xai and yai are coordinates of point A i respectively and l is length of the linkage.

media/image18.jpeg

Figure 3.

Schematic diagram of loop close equation.

Also, the x and y coordinates of point C i , i.e., the left-hand side of equations (4), can be formulated using platform coordinates as

xci=xp+xcicosφycisinφ
(7)
yci=yp+xcisinφ+ycicosφ
(8)
xci and yci are x and y coordinates of Ci respectively measured from the mass center of the platform, P, when φ is zero, as shown in Fig. 4.
media/image24.jpeg

Figure 4.

Platform coordinates.

From equations (5 -8), a closed-form solution is calculated as

ρi=Mi±l2Si2      i=1,2,3
(9)

where:

Mi =(xcixai)cosαi+(yciyai)sinαi Si =(xcixai)sinαi(yciyai)cosαi

Since there are two possible solutions for each chain, this manipulator can take on a maximum of eight configurations for a set of given coordinates of the platform. Note, only if the argument of the square root in equation (9) becomes zero, dose equation (9) have a unique solution. If the argument turns out to be negative, there is no solution to satisfy given kinematic requirements. Compared with the inverse kinematic solution described above, forward kinematic solutions of a planar parallel manipulator are much more difficult to solve (Merlet, 1996).

3. Dynamic analysis of flexible linkages

As industry demands high-speed machines, and hence lightweight linkages which deform under high inertial forces, we must consider structural flexibility of linkages in modeling a parallel manipulator. A single flexible link has been modeled in (Bellezza, 1990) and a serial type manipulator with both rigid and flexible links has been presented in (Low & Vidyasagar, 1988). Flexible models of a parallel manipulator have been studied in (Fattah et al., 1995).

A coordinate system for the flexible model is identical with the rigid-body model shown in Fig. 2. Only difference is the existence of the lateral deformation, w i (l), at the distal end of the i th linkage, C i , due to flexibility of the linkage, as shown in Fig. 5. Out of X-Y plane deformations are not considered here. This is the subject of another analysis. If length of the linkage, l, is much longer than thickness of the linkage, the linkage can be treated as an Euler-Bernoulli beam (Genta, 1993).

media/image28.jpeg

Figure 5.

Coordinate system of structurally-flexible model.

Coupled with flexible deformations, the kinematic equations (5 -6) are hence modified as follows;

xci=xai+ρicosαi+lcosβiwi(l)sinβi
(10)
yci=yai+ρisinαi+lsinβi+wi(l)cosβi
(11)

Since the left-hand side of equation (4) remains valid for a flexible model,

xci=xp+xcicosφycisinφ
(12)
yci=yp+xcisinφ+ycicosφ
(13)

Through Equations (10 -13), the inverse kinematic solution of a structurally flexible manipulator is formulated as

ρi=Mi±l2+wi(l)2Si2           i=1,2,3
(14)

where:

Mi=(xcixai)cosαi+(yciyai)sinαi Si=(xcixai)sinαi(yciyai)cosαi

Comparing with the inverse kinematic solution of the rigid-body model, equation(9), the linkage deformation is added in the right-hand side of equation (14). In addition, large linkage deformation may lead to no solution to ρi , i=1,2,3, because the argument of the square root in equation (14) has a negative value.

Evaluation of the derivative of equations (10 -13), with respect to time, gives

(x˙Pi¯+y˙Pj¯)+φ˙(k¯×e¯i)=ρ˙ia¯i+(β˙i+w˙i(l)/l)(k¯×b¯i)
(15)

where i¯ , j¯ are unit vectors along the reference X-Y frame respectively and e¯i is shown in Fig. 5. Dot-multiplication of equation (15) by b¯i leads to

ρ˙i=1a¯ib¯i[bixbiyeixbiyeiybix][x˙Py˙Pφ˙]T:=JPi X¯˙P
(16)

Subscripts of a position vector, named as x and y, represent X-directional and Y-directional components of the corresponding vector respectively. Cross-multiplication of equation (15) by b¯i gives

β˙i=1l2{[biybixeixbix+eiybiy](b¯i×a¯i)JPi}X¯˙Pw˙i(l)/l
(17)

Accelerations of the sliders and the links are given respectively by

ρ¨i=1a¯ib¯i{[bixbiyeixbiyeiybix]X¯¨Pb¯ie¯iφ˙2+β˙i2l2+β˙iw˙i(l)l}
(18)
β¨i=1l2{[biybixeixbix+eiybiy]X¯¨P(b¯i×e¯i)φ˙2(b¯i×a¯i)ρ¨i}w¨i(l)/l
(19)

Since three linkages in this analysis are assumed to have structural flexibility, the linkage deforms under high acceleration, as shown in Fig. 5. Flexible deformations can be expressed by the product of time-dependant functions and position-dependant functions, i.e. an assumed modes model (Genta, 1993);

wi(x,t):=j=1rη(t)ijψj(ξ)         i=1,2,3
(20)

where:

ξ:=x/l , r :=the number of assumed modes.

Functions η (t) can be considered the generalized coordinates expressing the deformation of the linkage and functions ψ(ξ) are referred to as assumed modes.

Considering boundary conditions of the linkage on B i and C i , their behavior is close to a pin (B i )-free (C i ) motion. Normalized shape functions, satisfying this boundary condition, are selected as:

ψj(ξ):=12sin(γj)[sin(γjξ)+sin(γj)sinh(γj)sinh(γjξ)]
(21)

where:

0ξ1 and γj:=(j+0.25)πl j=1,2,…,r

The first four shape functions are shown in Fig. 6 where the left end (B i ) exhibits zero deformation and the right end (C i ) exhibits a maximum deformation, as expected.

All generalized coordinates are collected to form of a single vector X defined as:

X:=[ρ¯β¯X¯Pη¯]TR9+3r
(22)

where:

β¯:=[β1β2β3]T η¯:=[η11η1rη21η2rη31η3r]T

Using inertia parameters of the manipulator and generalized coordinates, the kinetic energy of three sliders is written as

TS=i=1312msρ˙i2
(23)

The kinetic energy of the three links is expressed as

TL=i=1312ρA[ρ˙i2+(xβ˙i+w˙i)2+2ρ˙i(xβ˙i+w˙i)sin(αiβi)]dx
(24)

The kinetic energy of the platform is expressed as

TP=12mP(x˙P2+y˙P2)+12IPφ˙2
(25)

Therefore, collecting all kinetic energies, equations (23 -25), the total kinetic energy of the system is

T=i=1312msρ˙i2+i=1312ρA[ρi2+(xβ˙i+w˙i)2+2ρ˙i(xβ˙i+w˙i)sin(αiβi)]dx                                        +12mP(x˙P2+y˙P2)+12IPφ˙2
(26)

Since gravitational force is applied along Z-direction, perpendicular to the X-Y plane, potential energy due to gravitational force does not changed at all during any in-plane motions of the manipulator. Considering potential energy due to deformation of the linkage, total potential energy of the system is given as

V=i=13EI(wi)2dx
(27)

where: ρA :=mass per length of the linkage

E :=elastic modulus of the linkage

I :=area moment of inertia of the linkage

Evaluating Lagrangian equations of the first type given by

ddt(TX˙i)(TV)Xi=Qi+k=1mλkΓkXi;  i=1,2,, 9+3r
(28)

where: Qi :=generalized force

λk := k th Lagrange multiplier Γk := k th constrained equation

the left-hand side of equation (28) is formulated as follows:

ddt(Tρ˙i)(TV)ρi=(ms+m)ρ¨i+0.5mlsin(αiβi)β¨i+j=1rη¨ijsin(αiβi)ρAψjdx                0.5mlcos(αiβi)β˙i2j=1rη˙ijβ˙icos(αiβi)ρAψjdx    i=1,2,3
(29)
ddt(Tβ˙i)(TV)βi=0.5mlsin(αiβi)ρ¨i+ml2β¨i/3+j=1rη¨ijρAxψjdx                          +j=1rη˙ijρ˙icos(αiβi)ρAψjdx           i=1,2,3
(30)
ddt(TX˙P)(TV)XP=[mP000mP000IP][x¨Py¨Pφ¨]
(31)
ddt(Tη˙ij)(TV)ηij=sin(αiβi)ρ¨iρAψjdx+β¨iρAxψjdx+η¨ijρAψj2dx       ρAψjdxcos(αiβi)β˙iρ˙i+EI(ψj)2dx  i=1,2,3andj=1,2,..,r
(32)

Since the number of generalized coordinates excluding vibration modes is nine, greater than the number of the degrees-of-freedom of the manipulator, three, six constraint equations should be considered in equations of the motion. From the geometry of three closed-loop chains, equation (4), a fundamental constrained equation is given by

A¯iP¯+P¯C¯iA¯iB¯iB¯iC¯i=0  i=1,2,3
(33)

Dividing equations (33) into an X-axis’s component and a Y-axis’s component, six constraint equations are given by

Γ2i1:=ρicosαi+lcosβij=1rηijsinβixPrcos(φi+φ)=0
(34)
Γ2i:=ρisinαi+lsinβi+i=1rηijcosβiyPrsin(φi+φ)=0
(35)

where:

rcos(φi):=xci , rsin(φi):=yci i=1,2,3

From equation (34) and (35), the right-hand side of equation (28) is

Fai+k=16λkΓkρi=Fai+λ2i1cosαi+λ2isinαi    i=1,2,3
(36)
k=16λkΓkβi=λ2i1(lsinβij=1rηijcosβi)+λ2i(lcosβij=1rηijsinβi)    i=1,2,3
(37)
Fext+k=16λkΓkXP=Fext+[101010010101s31c31s32c32s33c33][λ1λ6]
(38)

where:

s3i:=eix , c3i:=eiy

F ext is an external force.

k=16λkΓkη1j=λ1sinβ1+λ2cosβ1     j=1,2,..,r
(39)
k=16λkΓkη2j=λ3sinβ2+λ4cosβ2    j=1,2,..,r
(40)
k=16λkΓkη3j=λ5sinβ3+λ6cosβ3    j=1,2,..,r
(41)

Putting equations (29 -32) and equations (36 -41) together, the equations of motion for the planar parallel manipulator are complete with a total of 9+3×r equations;

[M11M120M14M12TM220M2400M330M14TM24T0M44][ρ¯¨β¯¨X¯¨Pη¯¨]+[V1V20V4]+[000000000000000K][ρ¯β¯X¯Pη¯]=[Fa0Fext0]+[JΓ1JΓ2JΓ3JΓ4][λ1λ2λ3λ4λ5λ6]
(42)

where:

M11=(ms+m)[100010001] R3×3 M12=ml2[s1000s2000s3] R3×3

R3×3r

M22=ml23[100010001] R3×3 M33=[mp000mp000IP] R3×3 M24=ml[ψ1ξdξψrξdξ000000ψ1ξdξψrξdξ000000ψ1ξdξψrξdξ] R3×3r M44=m[M^000M^000M^] R3r×3r M^=[ψ12dξ00ψr2dξ] Rr×r K=EIl3[K^000K^000K^] R3r×3r K^=[ψ12dξ00ψr2dξ] Rr×r V1=[0.5mlc1β˙12+j=1rmη˙1jβ˙1c1ψjdξ0.5mlc2β˙22+j=1rmη˙2jβ˙2c2ψjdξ0.5mlc3β˙32+j=1rmη˙3jβ˙3c3ψjdξ] R3 V2=m[j=1rη˙1jρ˙1c1ψjdξj=1rη˙2jρ˙2c2ψjdξj=1rη˙3jρ˙3c3ψjdξ] R3 V4=m[c1ρ˙1β˙1ψ1dξc1ρ˙1β˙1ψrdξc2ρ˙2β˙2ψ1dξc2ρ˙2β˙2ψrdξc3ρ˙3β˙3ψ1dξc3ρ˙3β˙3ψrdξ]R3r

where : si:=sin(αiβi) and ci:=cos(αiβi)

JΓ1=[cosα1sinα1000000cosα2sinα2000000cosα3sinα3]R3×6 JΓ2=[s21c21000000s22c22000000s23c23]R3×6

where: s2i:=lsinβicosβij=1rηij and c2i:=lcosβisinβij=1rηij

JΓ3=[101010010101s31c31s32c32s33c33]R3×6

where: s3i:=eix , c3i:=eiy

Jr4=[sinβ1cosβ100sinβ1cosβ1sinβ2cosβ200sinβ2cosβ2sinβ3cosβ300sinβ3cosβ3]R3r×6

4. Active vibration control

If the intermediate linkages of the planar parallel manipulator are very stiff, an appropriate rigid body model based controller, such as a computed torque controller (Craig, 2003), can yield good trajectory tracking of the manipulator. However, structural flexibility of the linkages transfers unwanted vibration to the platform, and may even lead to instability of the whole system. Since control of linear motions of the sliders alone can not result in both precise tracking of the platform and vibration attenuation of the linkages simultaneously, an additional active damping method is proposed through the use of smart material. As discussed, the vibration damping controller proposed here is applied separately to a PVDF layer and PZT segments, and the performance of each actuator is then compared. Attached to the surface of the linkage, both of these piezoelectric materials generate shear force under applied control voltages, opposing shear stresses which arise due to elastic deformation of the linkages.

The integrated control system for the planar parallel manipulator proposed here consists of two components. The first component is a proportional and derivative (PD) feedback control scheme for the rigid body tracking of the platform as given below:

ui(t)=kP(ρdiρi)kd(ρ˙diρ˙i)           i=1,2,3
(43)

where k p and k d are a proportional and a derivative feedback gain respectively. ρdi and ρ˙di are desired displacement and velocity of the i th slider respectively. This signal is used as an input to electrical motors actuating ball-screw mechanisms for sliding motions. In the following, we introduce the second component of the integrated control system separately, for each of the piezoelectric materials examined, a PVDF layer and PZT segments, shown respectively in Fig. 7 and 8.

media/image116.jpeg

Figure 6.

Intermediate link with PVDF layer.

media/image117.jpeg

Figure 7.

Intermediate link with PZT actuator.

4.1. PVDF actuator control formulation

A PVDF layer can be bonded uniformly on the one side of the linkages of the planar parallel manipulator, as shown in Fig. 7. When a control voltage, v i , is applied to the PVDF layer, the virtual work done by the PVDF layer is expressed as

δWPVDF=cvi(t)j=1rψj(l)δηij
(44)

where c is a constant representing the bending moment per volt (Bailey & Hubbard, 1985) and l is the link length. ()' implies differentiation with respect to x. If the control voltage applied to the PVDF layer, v i , is formulated as

vi(t)=kIw˙i(l,t)        i=1,2,3
(45)

the slope velocity of the linkages, w˙(l,t) , converges to zero, assuming no exogenous disturbances applied to the manipulator, hence vibration of the linkages is damped out. Since the slope velocity, w˙(l,t) , is not easily measured or estimated by a conventional sensor system, an alternative scheme, referred to as the L-type method (Sun & Mills, 1999), is proposed as follows:

vi(t)=kIw˙i(l,t)       i=1,2,3
(46)

Instead of the slope velocity, w˙(l,t) , the linear velocity, w˙(l,t) , is employed in formulation of the control law. The linear velocity, w˙(l,t) , can be calculated through the integration of the linear acceleration measured by an accelerometer installed at the distal end of the linkages, C i . The shape function, ψj(ξ=1) , and its derivative, ψj(ξ=1) , have same trend of variation at the distal end of the linkages, C i , in all vibration modes, as shown in Fig. 6 ;

ψj(ξ=1)ψj(ξ=1)0j=1,2,,r
(47)

Therefore, the control system maintains stability when employing the L-type method to formulate the control voltages, v i .

4.2. PZT actuator control formulation

PZT actuators are manufactured in relatively small sizes, hence several PZT segments can be bonded together to a flexible linkage to damp unwanted vibrations. Assuming that only one PZT segment is attached to each intermediate linkage of the planar parallel manipulator, as shown in Fig. 8, the virtual work done by the PZT actuator is expressed as

δWPZT=cvi(t)j=1r[ψj(a2)ψj(a1)]δηij
(48)
a1 and a2 denote the positions of the two ends of the PZT actuator measured from B i along the intermediate linkage, as shown in Fig. 8. As the PVDF layer is, the PZT actuator is controlled using the L-type method as
vi(t)=kI[w˙i(a2,t)w˙i(a1,t)]  i=1,2,3
(49)

In contrast to the PVDF layer bonded uniformly to the manipulator linkages, the performance of the L-type scheme for the PZT actuator depends on the location of the PZT actuator. In order to achieve stable control performance, the PZT actuator should be placed in a region along the length of the linkage i.e. x[a1,a2] as discussed in (Sun & Mills, 1999), where ψj(x) and ψj(x) have the same trend of variation,

(ψj(a2)ψj(a1))(ψj(a2)ψj(a1))0
(50)

As the number of vibration modes increases, it is difficult to satisfy the stability condition, given in equation(50), for higher vibration modes, since the physical length of a PZT actuator is not sufficiently small.

5. Simulation results

Simulations are performed to investigate vibrations of the planar parallel manipulator linkages and damping performance of both piezoelectric actuators used in the manipulator with structurally-flexible linkages. Specifications of the manipulator for simulations are listed in Table 1. The first three modes are considered in the dynamic model, i.e. r=3. A sinusoidal function with smooth acceleration and deceleration is chosen as the desired input trajectory of the platform;

xP=xftftxf2πsin(2πtft)
(51)

Considering the target-performance in an electrical assembly process, such as wire bonding in integrated circuit fabrication, the goal for the platform is designed to move linearly 2 mm (x f ) within 10 msec (t f ). Feedback gains of the control system for the slider actuators are listed in Table 2. The feedback gain for piezoelectric actuators, k I , is selected so that the control voltage, applied to the PVDF layer, does not exceed 600 Volts. A fourth order Runge-Kutta method was used to integrate the ordinary differential equations, given by Equation (42) at a control update rate of 1 msec, using MATLABTM software. Parameters of piezoelectric materials, currently manufactured, are listed in Table 3. The placement position of the PZT actuator is adjusted to a 1 =0.66, a 2 =0.91, so that the first two vibration modes satisfy the stability condition given in equation (50).

Results of the PVDF layer are shown in Figures 9-12. Figure 9 shows that the error profile of the manipulator platform exhibits large oscillation at the initial acceleration, but continuously decreases due to the damping effect of the PVDF layer applied to the flexible linkages. The error profile of the platform without either of PVDF or PZT, labeled as “no damping” in Figure 9, shows typical characteristics of an undamped system. With Figure 10 showing deformation of the linkages on C i , it reveals that the PVDF layer can damp structural vibration of the linkages in a gradual way. The first three vibration modes are illustrated in Figure 11. The first mode has twenty times the amplitude than the second mode, and one hundred times the amplitude than the third mode. The control output for the first slider actuator is shown in the upper plot of Figure 12, and control voltage for the first PVDF layer is shown in the lower plot of Figure 12. The control voltage, applied to the PVDF layer, decreases as the amplitude of vibration does.

Results of the PZT actuator are shown in Figures 13-17. Comparing Figure 13 with Figure 9, the PZT actuator exhibits better damping performance than the PVDF layer. The error profile of the platform, with the PZT actuator activated, enters steady state quickly and does not exhibit any vibration in steady state. The structural vibrations of the linkages, illustrated in Figure 14, are completely damped after 60 msec. The first three vibration modes are shown in Figure 15. The first mode has ten times the amplitudes than the other modes. Since the PZT actuator has higher strain constant than the PVDF, the PZT actuator can generate large shear force with relatively small voltage applied. The maximum voltage of the lower plot of Figure 16 is about 200 Volts, while that of the Figure 12 reaches 600 Volts. Due to the length of the linkage and the PZT actuator applied to the linkage, only the first two modes satisfy the stability condition, given by equation (50). However, this has little effect on damping performance, as shown in Figure 14 since the first two modes play dominant roles in vibration. If the placement of the PZT actuator change to a 1 =0.4, a 2 =0.65, only the first mode satisfies the stability condition, which leads to divergence of vibration modes, as shown in Figure 17.

6. Conclusion

In this chapter, the equations of motion for the planar parallel manipulator are formulated by applying the Lagrangian equation of the first type. Introducing Lagrangian multipliers simplifies the complexities due to multiple closed loop chains of the parallel mechanism and the structurally flexible linkages. An active damping approach applied to two different piezoelectric materials, which are used as actuators to damp unwanted vibrations of flexible linkages of a planar parallel manipulator. The proposed control is applied to PVDF layer and PZT segments. An integrated control system, consisting of a PD feedback controller, applied to electrical motors for rigid body motion control of the manipulator platform, and a L-type controller applied to piezoelectric actuators to damp unwanted linkage vibrations, is developed to permit the manipulator platform to follow a given trajectory while damping vibration of the manipulator linkages. With an L-type control scheme determining a control voltage applied, the piezoelectric materials have been shown to provide good damping performance, and eventually reduce settling time of the platform of the planar parallel manipulator. Simulation results show that the planar parallel manipulator, with the lightweight linkages, during rigid body motion, undergoes persistent vibration due to high acceleration and deceleration. Additionally, the PZT actuator yields better performance in vibration attenuation than the PVDF layer, but may enter an unstable state if the position of the PZT actuator on the linkage violates the stability condition for the dominant vibration modes. In the near future, we will perform vibration experiments with a prototype planar parallel manipulator based on presented simulation results.

Platformside length mass0. 1 m 0.2 kg
Sliderm ass0.2 kg
Linear guide (Ball-screw)s troke incline angle0. 4 m 1 5 0 o , 270 o , 30 o
L inkl ength density modulus cross-section0. 2 m 2770 kg/m 3 73 GPa 0.0 25 m(W) * 0.0 15 m(H)

Table 1.

Specification of the planar parallel manipulator.

k p10,000 N/m
k d500 N-sec/m
k I4,000 V-sec/m for PVDF 1,500 V-sec/m for PZT

Table 2.

Feedback control gains.

PVDFPZT
modulus length thickness width density d 312 GPa 0.2 m 0.28 mm 0.0 25 m 1800 kg/m 3 22 * 10 -12 m/V63 GPa 0.0 5 m 0.75 mm 0.0 25 m 7600 kg/m 3 110 * 10 -12 m/V

Table 3.

Parameters of piezoelectric materials.

media/image136.png

Figure 8.

Error profile of the platform (dotted: no damping, solid: with PVDF layer).

media/image137.png

Figure 9.

Flexible deformation of each link (dotted: no damping, solid: with PVDF layer).

media/image138.png

Figure 10.

The first three vibration modes of the first link (dotted: no damping, solid: with PVDF layer).

media/image139.png

Figure 11.

Control output for the first link.

media/image140.png

Figure 12.

Error profile of platform (dotted: no damping, solid: with PZT actuator).

media/image141.png

Figure 13.

Flexible deformation of each link (dotted: no damping, solid: with PZT actuator).

media/image142.png

Figure 14.

The first three vibration modes of the first link(dotted: no damping, solid: with PZT actuator).

media/image143.png

Figure 15.

Control output for the first link.

media/image144.png

Figure 16.

The first three vibration modes of the first link with PZT actuator located on inappropriate place.

References

1 - T. Arai, K. Cleary, K. Homma, H. Adachi, T. Nakamura, 1991 Development of parallel link manipulator for underground excavation task, Proceedings of International Symposium on Advanced Robot Technology, 541 548 , Tokyo, Japan, March 1991
2 - T. Bailey, J. Hubbard, 1985 Distributed piezoelectric-polymer active vibration control of a cantilever beam, Journal of Guidance, Control and Dynamics, 8 5 605 611 , 0731-5090
3 - F. Bellezza, L. Lanari, G. Ulivi, 1990 Exact modeling of the flexible slewing link, Proceedings of IEEE International Conference on Robotics and Automation, 734 739 , 0-8186-9061-5, Cincinnati, USA, May 1990, IEEE
4 - J. J. Craig, 2003 Introduction to Robotics, Prentice-Hall, 100387985069
5 - A. Fattah, J. Angeles, K. A. Misra, 1995 Dynamic of a 3-DOF spatial parallel manipulator with flexible links, Proceedings of IEEE International Conference on Robotics and Automation, 627 632 , 0-7803-1965-6, Nagoya, Japan, May 1995, IEEE
6 - G. Genta, 1993 Vibration of structures and machines, Springer, 100387985069
7 - C. Gosselin, J. Angeles, 1990 Singularity analysis of closed-loop kinematic chains, IEEE Transactions on Robotics and Automation, 6 3 281 290 , 1552-3098
8 - C. Gosselin, S. Lemieux, J. P. Merlet, 1996 A new architecture of planar three- degree-of-freedom parallel manipulator, Proceedings of IEEE International Conference on Robotics and Automation, 3738 3743 , 0-7803-2988, Minneapolis, USA, April 1996, IEEE
9 - I. Heerah, B. Kang, J. K. Mills, B. Benhabib, 2002 Architecture selection and singularity analysis of a 3-degree-of-freedom planar parallel manipulator, Proceedings of ASME 2002 Design Engineering Technical Conferences and Computers and Information in Engineering Conference, 1 6 , 791836037, Montreal, Canada, October 2002, ASME
10 - B. Kang, B. Yeung, J. K. Mills, 2002 Two-time scale controller design for a high speed planar parallel manipulator with structural flexibility, Robotica, 20 5 519 528 , 0263-5747
11 - K. Kozak, I. Ebert-Uphoff, W. E. Singhose, 2004 Locally linearized dynamic analysis of parallel manipulators and application of input shaping to reduce vibrations, ASME Journal of Mechanical Design, 126 1 156 168 , 1050-0472
12 - K. H. Low, M. Vidyasagar, 1988 A Lagrangian formulation of the dynamic model for flexible manipulator systems, ASME Journal of Dynamic Systems, Measurement, and Control, 110 175 181 , 0022-0434
13 - J. P. Merlet, 1996 Direct kinematics of planar parallel manipulators, Proceedings of IEEE International Conference on Robotics and Automation, 3744 3749 , 0-7803-2988, Minneapolis, USA, April 1996, IEEE
14 - D. Sun, J. K. Mills, 1999 PZT actuator placement for structural vibration damping of high speed manufacturing equipment, Proceedings of the American Control Conference, 1107 1111 , 0780349903, San Diego, USA, June 1999, IEEE
15 - T. Toyama, T. Shibukawa, K. Hattori, K. Otubo, M. Tsutsumi, 2001 Vibration analysis of parallel mechanism platform with tilting linear motion actuators, Journal of the Japan Society of Precision Engineering, 67 9 458 1462 , 0916-78X