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

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

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.$$A_i$$is the origin of the i^thlinear guide and$$B_i$$is the position of the i^thslider.$$C_i$$is the position of the revolute joint of the platform facing with the i^thlinkage 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

$${\overline{X}}_P:={\left[\begin{array}{ccc}{x}_P& y_P& \phi \end{array}\right]}^T$$E1

The displacement of the sliders from their origin,$$A_i$$, to B_i, are expressed as

$$\overline{\rho }:=\left[\begin{array}{ccc}{\rho }_1& {\rho }_2& {\rho }_3\end{array}\right]$$E2

$${\beta }_i$$is defined as the angle at B_ibetween the X-axis of the fixed frame and the i^thintermediate link and$${\alpha }_i$$is the angle at A_ibetween the X-axis of the fixed frame and the i^thlinear guide. The position and orientation of the i^thlink at its mass center can also be written as

$${\overline{X}}_i:={\left[\begin{array}{ccc}{x}_i& y_i& {\beta }_i\end{array}\right]}^T$$E3

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

$${\overline{A}}_i\overline{P}+\overline{P}{\overline{C}}_i={\overline{A}}_i{\overline{B}}_i+{\overline{B}}_i{\overline{C}}_{i}i =1,2,3$$E4

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

$$x_{ci}=x_{ai}+{\rho }_i\cos {\alpha }_i+l\cos {\beta }_i$$E5

$$y_{ci}=y_{ai}+{\rho }_i\sin {\alpha }_i+l\sin {\beta }_i$$E6

where$$x_{ai}$$and$$y_{ai}$$are coordinates of point A_irespectively and lis length of the linkage.

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

$$x_{ci}=x_p+{x^{\prime }}_{ci}\cos \phi -{y^{\prime }}_{ci}\sin \phi$$E7

$$y_{ci}=y_p+{x^{\prime }}_{ci}\sin \phi +{y^{\prime }}_{ci}\cos \phi$$E8

$${x^{\prime }}_{ci}$$and$${y^{\prime }}_{ci}$$are x and y coordinates of Ci respectively measured from the mass center of the platform, P, when$$\phi$$is zero, as shown in Fig. 4.

Figure 4.

Platform coordinates.

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

$${\rho }_i=M_i\pm \sqrt{l^2-S_i^2} i=1,2,3$$E9

where:

$$M_i =(x_{ci}-x_{ai})\cos {\alpha }_i+(y_{ci}-y_{ai})\sin {\alpha }_i$$$$S_{i }=(x_{ci}-x_{ai})\sin {\alpha }_i-(y_{ci}-y_{ai})\cos {\alpha }_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^thlinkage, 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).

Figure 5.

Coordinate system of structurally-flexible model.

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

$$x_{ci}=x_{ai}+{\rho }_i\cos {\alpha }_i+l\cos {\beta }_i-w_i(l)\sin {\beta }_i$$E10

$$y_{ci}=y_{ai}+{\rho }_i\sin {\alpha }_i+l\sin {\beta }_i+w_i(l)\cos {\beta }_i$$E11

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

$$x_{ci}=x_p+{x^{\prime }}_{ci}\cos \phi -{y^{\prime }}_{ci}\sin \phi$$E12

$$y_{ci}=y_p+{x^{\prime }}_{ci}\sin \phi +{y^{\prime }}_{ci}\cos \phi$$E13

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

$${\rho }_i=M_i\pm \sqrt{l^2+w_i{(l)}^2-S_i^2} i=1,2,3$$E14

where:

$$M_i=(x_{ci}-x_{ai})\cos {\alpha }_i+(y_{ci}-y_{ai})\sin {\alpha }_i$$$$S_i=(x_{ci}-x_{ai})\sin {\alpha }_i-(y_{ci}-y_{ai})\cos {\alpha }_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$${\rho }_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

$$({\dot{x}}_P\overline{i}+{\dot{y}}_P\overline{j})+\dot{\phi }(\overline{k}\times {\overline{e}}_i)={\dot{\rho }}_i{\overline{a}}_i+({\dot{\beta }}_i+{\dot{w}}_i(l)/l)(\overline{k}\times {\overline{b}}_i)$$E15

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

$${\dot{\rho }}_i=\frac{1}{{\overline{a}}_i\cdot {\overline{b}}_i}\left[\begin{array}{ccc}{b}_{ix}& b_{iy}& e_{ix}{b}_{iy}-e_{iy}{b}_{ix}\end{array}\right]{\left[\begin{array}{ccc}{\dot{x}}_P& {\dot{y}}_P& \dot{\phi }\end{array}\right]}^T:=J_{Pi} {\dot{\overline{X}}}_P$$E16

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$${\overline{b}}_i$$gives

$${\dot{\beta }}_i=\frac{1}{l^2}\{\left[-\begin{array}{ccc}{b}_{iy}& b_{ix}& e_{ix}{b}_{ix}+e_{iy}{b}_{iy}\end{array}\right]-({\overline{b}}_i\times {\overline{a}}_i)J_{Pi}\}{\dot{\overline{X}}}_P-{\dot{w}}_i(l)/l$$E17

Accelerations of the sliders and the links are given respectively by

$${\ddot{\rho }}_i=\frac{1}{{\overline{a}}_i\cdot {\overline{b}}_i}\{\left[\begin{array}{ccc}{b}_{ix}& b_{iy}& e_{ix}{b}_{iy}-e_{iy}{b}_{ix}\end{array}\right]{\ddot{\overline{X}}}_P-{\overline{b}}_i\cdot {\overline{e}}_i{\dot{\phi }}^2+{\dot{\beta }}_i^2{l}^2+{\dot{\beta }}_i{\dot{w}}_i(l)l\}$$E18

$${\ddot{\beta }}_i=\frac{1}{l^2}\{\left[\begin{array}{ccc}-b_{iy}& b_{ix}& e_{ix}{b}_{ix}+e_{iy}{b}_{iy}\end{array}\right]{\ddot{\overline{X}}}_P-({\overline{b}}_i\times {\overline{e}}_i){\dot{\phi }}^2-({\overline{b}}_i\times {\overline{a}}_i){\ddot{\rho }}_i\}-{\ddot{w}}_i(l)/l$$E19

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

$$w_i(x,t):={\displaystyle \sum _{j=1}^r\eta {(t)}_{ij}{\psi }_j(\xi )} i=1,2,3$$E20

where:

$$\xi :=x/l$$, r:=the number of assumed modes.

Functions$$\eta$$(t) can be considered the generalized coordinates expressing the deformation of the linkage and functions$$\psi \left(\xi \right)$$are referred to as assumed modes.

Considering boundary conditions of the linkage on B_iand 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:

$${\psi }_j(\xi ):=\frac{1}{2\sin ({\gamma }_j)}[\sin ({\gamma }_j\xi )+\frac{\sin ({\gamma }_j)}{\sinh ({\gamma }_j)}\sinh ({\gamma }_j\xi )]$$E21

where:

$$0\le \xi \le 1$$and$${\gamma }_j:=(j+0.25)\pi 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 Xdefined as:

$$X:={\left[\begin{array}{cccc}\overline{\rho }& \overline{\beta }& {\overline{X}}_P& \overline{\eta }\end{array}\right]}^T\in R^{9+3r}$$E22

where:

$$\overline{\beta }:={\left[\begin{array}{ccc}{\beta }_1& {\beta }_2& {\beta }_3\end{array}\right]}^T$$$$\overline{\eta }:={\left[\begin{array}{ccccccccc}{\eta }_{11}& \cdots & {\eta }_{1r}& {\eta }_{21}& \cdots & {\eta }_{2r}& {\eta }_{31}& \cdots & {\eta }_{3r}\end{array}\right]}^T$$

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

$$T_S={\displaystyle \sum _{i=1}^3\frac{1}{2}{m}_s{\dot{\rho }}_i{}^2}$$E23

The kinetic energy of the three links is expressed as

$$T_L={\displaystyle \sum _{i=1}^3\frac{1}{2}{\displaystyle \int {\rho }_A[}{\dot{\rho }}_i{}^2+{(x{\dot{\beta }}_i+{\dot{w}}_i)}^2+2{\dot{\rho }}_i(x{\dot{\beta }}_i+{\dot{w}}_i)\sin ({\alpha }_i-{\beta }_i)]dx}$$E24

The kinetic energy of the platform is expressed as

$$T_P=\frac{1}{2}{m}_P({\dot{x}}_P{}^2+{\dot{y}}_P{}^2)+\frac{1}{2}{I}_P{\dot{\phi }}^2$$E25

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

$$\begin{array}{l}T={\displaystyle \sum _{i=1}^3\frac{1}{2}{m}_s{\dot{\rho }}_i{}^2}+{\displaystyle \sum _{i=1}^3\frac{1}{2}{\displaystyle \int {\rho }_A[}{\rho }_i{}^2+{(x{\dot{\beta }}_i+{\dot{w}}_i)}^2+2{\dot{\rho }}_i(x{\dot{\beta }}_i+{\dot{w}}_i)\sin ({\alpha }_i-{\beta }_i)]dx}\\ +\frac{1}{2}{m}_P({\dot{x}}_P{}^2+{\dot{y}}_P{}^2)+\frac{1}{2}{I}_P{\dot{\phi }}^2\end{array}$$E26

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={\displaystyle \sum _{i=1}^3{\displaystyle \int EI{({w^{″}}_i)}^{2}dx}}$$E27

where:$${\rho }_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

$$\frac{d}{dt}(\frac{\partial T}{\partial {\dot{X}}_i})-\frac{\partial (T-V)}{\partial X_i}=Q_i+{\displaystyle \sum _{k=1}^m{\lambda }_k\frac{\partial {\Gamma }_k}{\partial X_i}}; i=1,2,\dots ,9+3r$$E28

where:$$Q_i$$:=generalized force

$${\lambda }_k$$:= k^thLagrange multiplier$${\Gamma }_k$$:= k^thconstrained equation

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

$$\begin{array}{l}\frac{d}{dt}(\frac{\partial T}{\partial {\dot{\rho }}_i})-\frac{\partial (T-V)}{\partial {\rho }_i}=(m_s+m){\ddot{\rho }}_i+0.5ml\sin ({\alpha }_i-{\beta }_i){\ddot{\beta }}_i+{\displaystyle \sum _{j=1}^r{\ddot{\eta }}_{ij}\sin ({\alpha }_i-{\beta }_i){\displaystyle \int {\rho }_A{\psi }_{j}dx}}\\ -0.5ml\cos ({\alpha }_i-{\beta }_i){\dot{\beta }}_i{}^2-{\displaystyle \sum _{j=1}^r{\dot{\eta }}_{ij}{\dot{\beta }}_i\cos ({\alpha }_i-{\beta }_i){\displaystyle \int {\rho }_A{\psi }_{j}dx}} i=1,2,3\end{array}$$E29

$$\begin{array}{l}\frac{d}{dt}(\frac{\partial T}{\partial {\dot{\beta }}_i})-\frac{\partial (T-V)}{\partial {\beta }_i}=0.5ml\sin ({\alpha }_i-{\beta }_i){\ddot{\rho }}_i+m{l}^2{\ddot{\beta }}_i/3+{\displaystyle \sum _{j=1}^r{\ddot{\eta }}_{ij}{\displaystyle \int {\rho }_{A}x{\psi }_{j}dx}}\\ +{\displaystyle \sum _{j=1}^r{\dot{\eta }}_{ij}{\dot{\rho }}_i\cos ({\alpha }_i-{\beta }_i){\displaystyle \int {\rho }_A{\psi }_{j}dx}} i=1,2,3\end{array}$$E30

$$\frac{d}{dt}(\frac{\partial T}{\partial {\dot{X}}_P})-\frac{\partial (T-V)}{\partial X_P}=\left[\begin{array}{ccc}{m}_P& 0& 0\\ 0& m_P& 0\\ 0& 0& I_P\end{array}\right]\left[\begin{array}{c}{\ddot{x}}_P\\ {\ddot{y}}_P\\ \ddot{\phi }\end{array}\right]$$E31

$$\begin{array}{l}\frac{d}{dt}(\frac{\partial T}{\partial {\dot{\eta }}_{ij}})-\frac{\partial (T-V)}{\partial {\eta }_{ij}}=\sin ({\alpha }_i-{\beta }_i){\ddot{\rho }}_i{\displaystyle \int {\rho }_A{\psi }_{j}dx+{\ddot{\beta }}_i}{\displaystyle \int {\rho }_{A}x{\psi }_{j}dx+{\ddot{\eta }}_{ij}{\displaystyle \int {\rho }_A{\psi }_j{}^2}}dx\\ -{\displaystyle \int {\rho }_A{\psi }_{j}dx\cos ({\alpha }_i-{\beta }_i){\dot{\beta }}_i{\dot{\rho }}_i}+{\displaystyle \int EI({{\psi }^{″}}_j}{)}^{2}dx i=1,2,3andj=1,2,\mathrm{..},r\end{array}$$E32

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

$${\overline{A}}_i\overline{P}+\overline{P}{\overline{C}}_i-{\overline{A}}_i{\overline{B}}_i-{\overline{B}}_i{\overline{C}}_i=0 i=1,2,3$$E33

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

$${\Gamma }_{2i-1}:={\rho }_i\cos {\alpha }_i+l\cos {\beta }_i-{\displaystyle \sum _{j=1}^r{\eta }_{ij}\sin {\beta }_i-}{x}_P-r\cos ({\phi }_i+\phi )=0$$E34

$${\Gamma }_{2i}:={\rho }_i\sin {\alpha }_i+l\sin {\beta }_i+{\displaystyle \sum _{i=1}^r{\eta }_{ij}\cos {\beta }_i-}{y}_P-r\sin ({\phi }_i+\phi )=0$$E35

where:

$$r\cos ({\phi }_i):={x^{\prime }}_{ci}$$,$$r\sin ({\phi }_i):={y^{\prime }}_{ci}$$i=1,2,3

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

$$F_{ai}+{\displaystyle \sum _{k=1}^6{\lambda }_k}\frac{\partial {\Gamma }_k}{\partial {\rho }_i}=F_{ai}+{\lambda }_{2i-1}\cos {\alpha }_i+{\lambda }_{2i}\sin {\alpha }_i i=1,2,3$$E36

$${\displaystyle \sum _{k=1}^6{\lambda }_k}\frac{\partial {\Gamma }_k}{\partial {\beta }_i}={\lambda }_{2i-1}(-l\sin {\beta }_i-{\displaystyle \sum _{j=1}^r{\eta }_{ij}}\cos {\beta }_i)+{\lambda }_{2i}(l\cos {\beta }_i-{\displaystyle \sum _{j=1}^r{\eta }_{ij}}\sin {\beta }_i) i=1,2,3$$E37

$$F_{ext}+{\displaystyle \sum _{k=1}^6{\lambda }_k}\frac{\partial {\Gamma }_k}{\partial X_P}=F_{ext}+\left[\begin{array}{cccccc}-1& 0& -1& 0& -1& 0\\ 0& -1& 0& -1& 0& -1\\ s{3}_1& c{3}_1& s{3}_2& c{3}_2& s{3}_3& c{3}_3\end{array}\right]\left[\begin{array}{c}{\lambda }_1\\ ⋮\\ {\lambda }_6\end{array}\right]$$E38

where:

$$s{3}_i:=e_{ix}$$,$$c{3}_i:=e_{iy}$$

F_extis an external force.

$${\displaystyle \sum _{k=1}^6{\lambda }_k}\frac{\partial {\Gamma }_k}{\partial {\eta }_{1j}}=-{\lambda }_1\sin {\beta }_1+{\lambda }_2\cos {\beta }_1 j=1,2,\mathrm{..},r$$E39

$${\displaystyle \sum _{k=1}^6{\lambda }_k}\frac{\partial {\Gamma }_k}{\partial {\eta }_{2j}}=-{\lambda }_3\sin {\beta }_2+{\lambda }_4\cos {\beta }_2 j=1,2,\mathrm{..},r$$E40

$${\displaystyle \sum _{k=1}^6{\lambda }_k}\frac{\partial {\Gamma }_k}{\partial {\eta }_{3j}}=-{\lambda }_5\sin {\beta }_3+{\lambda }_6\cos {\beta }_3 j=1,2,\mathrm{..},r$$E41

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\times r$$equations;

$$\left[\begin{array}{cccc}{M}_{11}& M_{12}& 0& M_{14}\\ M_{12}^T& M_{22}& 0& M_{24}\\ 0& 0& M_{33}& 0\\ M_{14}^T& M_{24}^T& 0& M_{44}\end{array}\right]\left[\begin{array}{c}\ddot{\overline{\rho }}\\ \ddot{\overline{\beta }}\\ {\ddot{\overline{X}}}_P\\ \ddot{\overline{\eta }}\end{array}\right]+\left[\begin{array}{c}{V}_1\\ V_2\\ 0\\ V_4\end{array}\right]+\left[\begin{array}{cccc}0& 0& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& K\end{array}\right]\left[\begin{array}{c}\overline{\rho }\\ \overline{\beta }\\ {\overline{X}}_P\\ \overline{\eta }\end{array}\right]=\left[\begin{array}{c}{F}_a\\ 0\\ F_{ext}\\ 0\end{array}\right]+\left[\begin{array}{c}{J}_{\Gamma 1}\\ J_{\Gamma 2}\\ J_{\Gamma 3}\\ J_{\Gamma 4}\end{array}\right]\left[\begin{array}{c}{\lambda }_1\\ {\lambda }_2\\ {\lambda }_3\\ {\lambda }_4\\ {\lambda }_5\\ {\lambda }_6\end{array}\right]$$E42

where:

$$M_{11}=(m_s+m)\left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right]$$$$\in R^{3\times 3}$$$$M_{12}=\frac{ml}{2}\left[\begin{array}{ccc}{s}_1& 0& 0\\ 0& s_2& 0\\ 0& 0& s_3\end{array}\right]$$$$\in R^{3\times 3}$$

$$\in R^{3\times 3r}$$

$$M_{22}=\frac{m{l}^2}{3}\left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right]$$$$\in R^{3\times 3}$$$$M_{33}=\left[\begin{array}{ccc}{m}_p& 0& 0\\ 0& m_p& 0\\ 0& 0& I_P\end{array}\right]$$$$\in R^{3\times 3}$$$$M_{24}=ml\left[\begin{array}{ccccccccc}{\displaystyle \int {\psi }_1\xi d\xi }& \cdots & {\displaystyle \int {\psi }_r\xi d\xi }& 0& \cdots & 0& 0& \cdots & 0\\ 0& \cdots & 0& {\displaystyle \int {\psi }_1\xi d\xi }& \cdots & {\displaystyle \int {\psi }_r\xi d\xi }& 0& \cdots & 0\\ 0& \cdots & 0& 0& \cdots & 0& {\displaystyle \int {\psi }_1\xi d\xi }& \cdots & {\displaystyle \int {\psi }_r\xi d\xi }\end{array}\right]$$$$\in R^{3\times 3r}$$$$M_{44}=m\left[\begin{array}{ccc}\widehat{M}& 0& 0\\ 0& \widehat{M}& 0\\ 0& 0& \widehat{M}\end{array}\right]$$$$\in R^{3r\times 3r}$$$$\widehat{M}=\left[\begin{array}{ccc}{\displaystyle \int {\psi }_1^{2}d\xi }& \cdots & 0\\ ⋮& \cdots & ⋮\\ 0& \cdots & {\displaystyle \int {\psi }_r^{2}d\xi }\end{array}\right]$$$$\in R^{r\times r}$$$$K=\frac{EI}{l^3}\left[\begin{array}{ccc}\widehat{K}& 0& 0\\ 0& \widehat{K}& 0\\ 0& 0& \widehat{K}\end{array}\right]$$$$\in R^{3r\times 3r}$$$$\widehat{K}=\left[\begin{array}{ccc}{\displaystyle \int {{\psi }^{″}}_1^{2}d\xi }& \cdots & 0\\ ⋮& \cdots & ⋮\\ 0& \cdots & {\displaystyle \int {{\psi }^{″}}_r^{2}d\xi }\end{array}\right]$$$$\in R^{r\times r}$$$$V_1=-\left[\begin{array}{c}0.5ml{c}_1{\dot{\beta }}_1^2+{\displaystyle \sum _{j=1}^{r}m{\dot{\eta }}_{1j}{\dot{\beta }}_1{c}_1{\displaystyle \int {\psi }_{j}d\xi }}\\ 0.5ml{c}_2{\dot{\beta }}_2^2+{\displaystyle \sum _{j=1}^{r}m{\dot{\eta }}_{2j}{\dot{\beta }}_2{c}_2{\displaystyle \int {\psi }_{j}d\xi }}\\ 0.5ml{c}_3{\dot{\beta }}_3^2+{\displaystyle \sum _{j=1}^{r}m{\dot{\eta }}_{3j}{\dot{\beta }}_3{c}_3{\displaystyle \int {\psi }_{j}d\xi }}\end{array}\right]$$$$\in R^3$$$$V_2=m\left[\begin{array}{c}{\displaystyle \sum _{j=1}^r{\dot{\eta }}_{1j}{\dot{\rho }}_1{c}_1{\displaystyle \int {\psi }_{j}d\xi }}\\ {\displaystyle \sum _{j=1}^r{\dot{\eta }}_{2j}{\dot{\rho }}_2{c}_2{\displaystyle \int {\psi }_{j}d\xi }}\\ {\displaystyle \sum _{j=1}^r{\dot{\eta }}_{3j}{\dot{\rho }}_3{c}_3{\displaystyle \int {\psi }_{j}d\xi }}\end{array}\right]$$$$\in R^3$$$$V_4=-m\left[\begin{array}{c}{c}_1{\dot{\rho }}_1{\dot{\beta }}_1{\displaystyle \int {\psi }_{1}d\xi }\\ ⋮\\ c_1{\dot{\rho }}_1{\dot{\beta }}_1{\displaystyle \int {\psi }_{r}d\xi }\\ c_2{\dot{\rho }}_2{\dot{\beta }}_2{\displaystyle \int {\psi }_{1}d\xi }\\ ⋮\\ c_2{\dot{\rho }}_2{\dot{\beta }}_2{\displaystyle \int {\psi }_{r}d\xi }\\ c_3{\dot{\rho }}_3{\dot{\beta }}_3{\displaystyle \int {\psi }_{1}d\xi }\\ ⋮\\ c_3{\dot{\rho }}_3{\dot{\beta }}_3{\displaystyle \int {\psi }_{r}d\xi }\end{array}\right]\in R^{3r}$$

where :$$s_i:=\sin ({\alpha }_i-{\beta }_i)$$and$$c_i:=\cos ({\alpha }_i-{\beta }_i)$$

$$J_{\Gamma 1}=\left[\begin{array}{cccccc}\cos {\alpha }_1& \sin {\alpha }_1& 0& 0& 0& 0\\ 0& 0& \cos {\alpha }_2& \sin {\alpha }_2& 0& 0\\ 0& 0& 0& 0& \cos {\alpha }_3& \sin {\alpha }_3\end{array}\right]\in R^{3\times 6}$$$$J_{\Gamma 2}=\left[\begin{array}{cccccc}s{2}_1& c{2}_1& 0& 0& 0& 0\\ 0& 0& s{2}_2& c{2}_2& 0& 0\\ 0& 0& 0& 0& s{2}_3& c{2}_3\end{array}\right]\in R^{3\times 6}$$

where:$$s{2}_i:=-l\sin {\beta }_i-\cos {\beta }_i{\displaystyle \sum _{j=1}^r{\eta }_{ij}}$$and$$c{2}_i:=l\cos {\beta }_i-\sin {\beta }_i{\displaystyle \sum _{j=1}^r{\eta }_{ij}}$$

$$J_{\Gamma 3}=\left[\begin{array}{cccccc}-1& 0& -1& 0& -1& 0\\ 0& -1& 0& -1& 0& -1\\ s{3}_1& -c{3}_1& s{3}_2& -c{3}_2& s{3}_3& -c{3}_3\end{array}\right]\in R^{3\times 6}$$

where:$$s{3}_i:=e_{ix}$$,$$c{3}_i:=e_{iy}$$

$$J_{r4}=\left[\begin{array}{ccccccccc}-\sin {\beta }_1& & \cos {\beta }_1& & & & & & \\ ⋮& & ⋮& & & 0& & 0& \\ -\sin {\beta }_1& & \cos {\beta }_1& & & & & & \\ & & & -\sin {\beta }_2& & \cos {\beta }_2& & & \\ & 0& & ⋮& & ⋮& & 0& \\ & & & -\sin {\beta }_2& & \cos {\beta }_2& & & \\ & & & & & & -\sin {\beta }_3& & \cos {\beta }_3\\ & 0& & & 0& & ⋮& & ⋮\\ & & & & & & -\sin {\beta }_3& & \cos {\beta }_3\end{array}\right]\in R^{3r\times 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:

$$u_i(t)=-k_P({\rho }_{di}-{\rho }_i)-k_d({\dot{\rho }}_{di}-{\dot{\rho }}_i) i=1,2,3$$E43

where k_pand k_dare a proportional and a derivative feedback gain respectively.$${\rho }_{di}$$and$${\dot{\rho }}_{di}$$are desired displacement and velocity of the i^thslider 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, v_i, is applied to the PVDF layer, the virtual work done by the PVDF layer is expressed as

$$\delta W_{PVDF}=c{v}_i(t){\displaystyle \sum _{j=1}^r{{\psi }^{\prime }}_j(l)}\delta {\eta }_{ij}$$E44

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

$$v_i(t)=-k_I{{\dot{w}}^{\prime }}_i(l,t) i=1,2,3$$E45

the slope velocity of the linkages,$${\dot{w}}^{\prime }(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,$${\dot{w}}^{\prime }(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:

$$v_i(t)=-k_I{\dot{w}}_i(l,t) i=1,2,3$$E46

Instead of the slope velocity,$${\dot{w}}^{\prime }(l,t)$$, the linear velocity,$$\dot{w}(l,t)$$, is employed in formulation of the control law. The linear velocity,$$\dot{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,$${\psi }_j(\xi =1)$$, and its derivative,$${{\psi }^{\prime }}_j(\xi =1)$$, have same trend of variation at the distal end of the linkages, C_i, in all vibration modes, as shown in Fig. 6 ;

$${\psi }_j(\xi =1){{\psi }^{\prime }}_j(\xi =1)\ge 0j=1,2,---,r$$E47

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

$$\delta W_{PZT}=c{v}_i(t){\displaystyle \sum _{j=1}^r[{{\psi }^{\prime }}_j(a_2)-{{\psi }^{\prime }}_j}(a_1)]\delta {\eta }_{ij}$$E48

$$a_1$$and$$a_2$$denote the positions of the two ends of the PZT actuator measured from B_ialong the intermediate linkage, as shown in Fig. 8. As the PVDF layer is, the PZT actuator is controlled using the L-type method as

$$v_i(t)=-k_I[{\dot{w}}_i(a_2,t)-{\dot{w}}_i(a_1,t)] i=1,2,3$$E49

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\in [a_1,a_2]$$as discussed in (Sun & Mills, 1999), where$${\psi }_j(x)$$and$${{\psi }^{\prime }}_j(x)$$have the same trend of variation,

$$({\psi }_j(a_2)-{\psi }_j(a_1))({{\psi }^{\prime }}_j(a_2)-{{\psi }^{\prime }}_j(a_1))\ge 0$$E50

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;

$$x_P=\frac{x_f}{t_f}t-\frac{x_f}{2\pi }\sin (\frac{2\pi }{t_f}t)$$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(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 MATLAB^TM software. Parameters of piezoelectric materials, currently manufactured, are listed in Table 3. The placement position of the PZT actuator is adjusted to a₁=0.66, a₂=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₁=0.4, a₂=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.

Platform	side length mass	0. 1 m 0.2 kg
Slider	m ass	0.2 kg
Linear guide (Ball-screw)	s troke incline angle	0. 4 m 1 5 0 ^o , 270 ^o , 30 ^o
L ink	l ength density modulus cross-section	0. 2 m 2770 kg/m ^373 GPa 0.0 25 m(W) * 0.0 15 m(H)

Table 1.

Specification of the planar parallel manipulator.

k p	10,000 N/m
k d	500 N-sec/m
k I	4,000 V-sec/m for PVDF 1,500 V-sec/m for PZT

Table 2.

Feedback control gains.

	PVDF	PZT
modulus length thickness width density d 31	2 GPa 0.2 m 0.28 mm 0.0 25 m 1800 kg/m 3 22 * 10 -12 m/V	63 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.