Open access peer-reviewed chapter

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

Written By

Bongsoo Kang and James K. Mills

Published: April 1st, 2008

DOI: 10.5772/5380

Chapter metrics overview

3,278 Chapter Downloads

View Full Metrics

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.

Advertisement

2. Rigid-body analysis of a planar parallel manipulator

2.1. Architecture of a planar parallel manipulator

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

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.Aiis the origin of the ithlinear guide andBiis the position of the ithslider.Ciis the position of the revolute joint of the platform facing with the ithlinkage andPis 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φ]TE1

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

ρ¯:=[ρ1ρ2ρ3]E2
βiis defined as the angle at Bibetween the X-axis of the fixed frame and the ithintermediate link andαiis the angle at Aibetween the X-axis of the fixed frame and the ithlinear guide. The position and orientation of the ithlink at its mass center can also be written as
X¯i:=[xiyiβi]TE3

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,3E4

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 Ci, is written as

xci=xai+ρicosαi+lcosβiE5
yci=yai+ρisinαi+lsinβiE6

wherexaiandyaiare coordinates of point Airespectively and lis length of the linkage.

Figure 3.

Schematic diagram of loop close equation.

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

xci=xp+xcicosφycisinφE7
yci=yp+xcisinφ+ycicosφE8
xciandyciare x and y coordinates of Ci respectively measured from the mass center of the platform, P, whenφis zero, as shown in Fig. 4.

Figure 4.

Platform coordinates.

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

ρi=Mi±l2Si2      i=1,2,3E9

where:

Mi =(xcixai)cosαi+(yciyai)sinαiSi =(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).

Advertisement

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, wi(l), at the distal end of the ithlinkage, Ci, 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).

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βiE10
yci=yai+ρisinαi+lsinβi+wi(l)cosβiE11

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

xci=xp+xcicosφycisinφE12
yci=yp+xcisinφ+ycicosφE13

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,3E14

where:

Mi=(xcixai)cosαi+(yciyai)sinαiSi=(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)E15

wherei¯,j¯are unit vectors along the reference X-Y frame respectively ande¯iis shown in Fig. 5. Dot-multiplication of equation (15) byb¯ileads to

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

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) byb¯igives

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

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}E18
β¨i=1l2{[biybixeixbix+eiybiy]X¯¨P(b¯i×e¯i)φ˙2(b¯i×a¯i)ρ¨i}w¨i(l)/lE19

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,3E20

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 Biand Ci, their behavior is close to a pin (Bi)-free (Ci) motion. Normalized shape functions, satisfying this boundary condition, are selected as:

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

where:

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

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

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

X:=[ρ¯β¯X¯Pη¯]TR9+3rE22

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ρ˙i2E23

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)]dxE24

The kinetic energy of the platform is expressed as

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

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φ˙2E26

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)2dxE27

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+3rE28

where:Qi:=generalized force

λk:= kthLagrange multiplierΓk:= kthconstrained 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,3E29
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,3E30
ddt(TX˙P)(TV)XP=[mP000mP000IP][x¨Py¨Pφ¨]E31
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,..,rE32

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,3E33

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+φ)=0E34
Γ2i:=ρisinαi+lsinβi+i=1rηijcosβiyPrsin(φi+φ)=0E35

where:

rcos(φi):=xci,rsin(φi):=ycii=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,3E36
k=16λkΓkβi=λ2i1(lsinβij=1rηijcosβi)+λ2i(lcosβij=1rηijsinβi)    i=1,2,3E37
Fext+k=16λkΓkXP=Fext+[101010010101s31c31s32c32s33c33][λ1λ6]E38

where:

s3i:=eix,c3i:=eiy

Fextis an external force.

k=16λkΓkη1j=λ1sinβ1+λ2cosβ1     j=1,2,..,rE39
k=16λkΓkη2j=λ3sinβ2+λ4cosβ2    j=1,2,..,rE40
k=16λkΓkη3j=λ5sinβ3+λ6cosβ3    j=1,2,..,rE41

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

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

where:

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

R3×3r

M22=ml23[100010001]R3×3M33=[mp000mp000IP]R3×3M24=ml[ψ1ξdξψrξdξ000000ψ1ξdξψrξdξ000000ψ1ξdξψrξdξ]R3×3rM44=m[M^000M^000M^]R3r×3rM^=[ψ12dξ00ψr2dξ]Rr×rK=EIl3[K^000K^000K^]R3r×3rK^=[ψ12dξ00ψr2dξ]Rr×rV1=[0.5mlc1β˙12+j=1rmη˙1jβ˙1c1ψjdξ0.5mlc2β˙22+j=1rmη˙2jβ˙2c2ψjdξ0.5mlc3β˙32+j=1rmη˙3jβ˙3c3ψjdξ]R3V2=m[j=1rη˙1jρ˙1c1ψjdξj=1rη˙2jρ˙2c2ψjdξj=1rη˙3jρ˙3c3ψjdξ]R3V4=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)andci:=cos(αiβi)

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

where:s2i:=lsinβicosβij=1rηijandc2i:=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
Advertisement

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,3E43

where kpand kdare a proportional and a derivative feedback gain respectively.ρdiandρ˙diare desired displacement and velocity of the ithslider 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.

Figure 6.

Intermediate link with PVDF layer.

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, vi, is applied to the PVDF layer, the virtual work done by the PVDF layer is expressed as

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

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

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

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,3E46

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, Ci. The shape function,ψj(ξ=1), and its derivative,ψj(ξ=1), have same trend of variation at the distal end of the linkages, Ci, in all vibration modes, as shown in Fig. 6 ;

ψj(ξ=1)ψj(ξ=1)0j=1,2,,rE47

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

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)]δηijE48
a1anda2denote the positions of the two ends of the PZT actuator measured from Bialong 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,3E49

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))0E50

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.

Advertisement

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)E51

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(xf) within 10 msec(tf). Feedback gains of the control system for the slider actuators are listed in Table 2. The feedback gain for piezoelectric actuators, kI, 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 a1=0.66, a2=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 Ci, 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 a1=0.4, a2=0.65, only the first mode satisfies the stability condition, which leads to divergence of vibration modes, as shown in Figure 17.

Advertisement

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.

Figure 8.

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

Figure 9.

Flexible deformation of each 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).

Figure 11.

Control output for the first link.

Figure 12.

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

Figure 13.

Flexible deformation of each 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).

Figure 15.

Control output for the first link.

Figure 16.

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

References

  1. 1. AraiT.ClearyK.HommaK.AdachiH.NakamuraT.1991Development of parallel link manipulator for underground excavation task,Proceedings of International Symposium on Advanced Robot Technology,541548, Tokyo, Japan, March 1991
  2. 2. BaileyT.HubbardJ.1985Distributed piezoelectric-polymer active vibration control of a cantilever beam, Journal of Guidance, Control and Dynamics,85605611, 0731-5090
  3. 3. BellezzaF.LanariL.UliviG.1990Exact modeling of the flexible slewing link,Proceedings of IEEE International Conference on Robotics and Automation,734739, 0-8186-9061-5, Cincinnati, USA, May 1990, IEEE
  4. 4. CraigJ. J.2003Introduction to Robotics,Prentice-Hall,100387985069
  5. 5. FattahA.AngelesJ.MisraK. A.1995Dynamic of a 3-DOF spatial parallel manipulator with flexible links,Proceedings of IEEE International Conference on Robotics and Automation,627632, 0-7803-1965-6, Nagoya, Japan, May 1995, IEEE
  6. 6. GentaG.1993Vibration of structures and machines, Springer,100387985069
  7. 7. GosselinC.AngelesJ.1990Singularity analysis of closed-loop kinematic chains,IEEE Transactions on Robotics and Automation,63281290, 1552-3098
  8. 8. GosselinC.LemieuxS.MerletJ. P.1996A new architecture of planar three- degree-of-freedom parallel manipulator,Proceedings of IEEE International Conference on Robotics and Automation,37383743, 0-7803-2988, Minneapolis, USA, April 1996, IEEE
  9. 9. HeerahI.KangB.MillsJ. K.BenhabibB.2002Architecture 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,16, 791836037, Montreal, Canada, October 2002, ASME
  10. 10. KangB.YeungB.MillsJ. K.2002Two-time scale controller design for a high speed planar parallel manipulator with structural flexibility,Robotica,205519528, 0263-5747
  11. 11. KozakK.Ebert-UphoffI.SinghoseW. E.2004Locally linearized dynamic analysis of parallel manipulators and application of input shaping to reduce vibrations,ASME Journal of Mechanical Design,1261156168, 1050-0472
  12. 12. LowK. H.VidyasagarM.1988A Lagrangian formulation of the dynamic model for flexible manipulator systems,ASME Journal of Dynamic Systems, Measurement, and Control,110175181, 0022-0434
  13. 13. MerletJ. P.1996Direct kinematics of planar parallel manipulators,Proceedings of IEEE International Conference on Robotics and Automation,37443749, 0-7803-2988, Minneapolis, USA, April 1996, IEEE
  14. 14. SunD.MillsJ. K.1999PZT actuator placement for structural vibration damping of high speed manufacturing equipment,Proceedings of the American Control Conference,11071111, 0780349903, San Diego, USA, June 1999, IEEE
  15. 15. ToyamaT.ShibukawaT.HattoriK.OtuboK.TsutsumiM.2001Vibration analysis of parallel mechanism platform with tilting linear motion actuators,Journal of the Japan Society of Precision Engineering,6794581462, 0916-78X

Written By

Bongsoo Kang and James K. Mills

Published: April 1st, 2008