Nonlinear characteristic analysis of high-speed spatial parallel mechanism

Chen, Xiulong; Li, Yunfeng; Deng, Yu; Li, Wenbin; Liang, Xiaoxia; Wang, Qing

Journal of Vibroengineering

Browse Journal

Submit article

Published: 15 November 2015

Check for updates

Nonlinear characteristic analysis of high-speed spatial parallel mechanism

Xiulong Chen¹

Yunfeng Li²

Yu Deng³

Wenbin Li⁴

Xiaoxia Liang⁵

Qing Wang⁶

^{1, 2, 3, 4, 5}College of Mechanical and Electronic Engineering, Shandong University of Science and Technology, Qingdao 266590, China

⁶College of Civil Engineering and Architecture, Shandong University of Science and Technology, Qingdao 266590, China

Corresponding Author:

Xiulong Chen

Cite the article Download PDF

Downloads 1217

Abstract

In order to grasp the nonlinear characteristic of high-speed spatial parallel mechanism, the analysis of nonlinear characteristics for spatial parallel mechanism is investigated. The nonlinear elastic dynamic equation of 4-UPS-UPU high-speed spatial parallel mechanism is derived by kineto-elastodynamics theory, the dynamic equation is solved by numerical method, the nonlinear characteristic of the parallel mechanism is analyzed by phase diagram. Numerical results show that 4-UPS-UPU high-speed spatial parallel mechanism exhibits typical nonlinear characteristic during exercise, the factors, such as the motion trajectory of parallel mechanism, the material of driving limbs, the diameter of driving limbs and the mass of moving platform, are also have effect on nonlinear characteristics of parallel mechanism. Therefore the reasonable choice of the above factors can weaken the chaos motion. This researches provide important theoretical base of the chaos suppression for spatial parallel mechanism.

1. Introduction

The spatial parallel mechanism, which has advantages of high speed, high precision, good dynamic performance, high stiffness and tight construction, has been paid more attention [1, 2]. The stiffness of all parts in spatial parallel mechanism is different, for example, the stiffness of telescopic rod in driving limbs is far below the stiffness of moving platform. In high speed and high acceleration conditions, the spatial parallel mechanism is a rigid-flexible dynamic coupled system due to the geometry nonlinearity and coupling effect between rigid body motion and elastic deformation, and exhibits strong nonlinear characteristics. Therefore, in order to reveal and improve the dynamic behaviors, meet the needs of the design, the nonlinear characteristics of spatial parallel mechanisms must be investigated.

Up to now, the overwhelming majority of the analysis of elastodynamic behaviors for parallel mechanism are implemented within the linear theory framework [3-8], so the results obtained, which have bigger error [9-11], can not reflect the influencing factors of dynamic behaviors and actual dynamic behaviors, and then affect the dynamic optimal design. Besides, the existing studies are mainly focused on planar parallel mechanism, while less on spatial parallel mechanism. And that the study of the nonlinear characteristics of high-speed spatial parallel mechanism is quite few.

The elastic dynamics equation of 4-UPS-UPU 5-DOF high-speed spatial parallel mechanism (see Fig. 1) is derived by kineto-elastodynamics theory in this paper, the nonlinear characteristics are analyzed by Newmark method and phase diagram, the relationship between nonlinear characteristics and the factors, such as the motion trajectory of parallel mechanism, the material of driving limbs, the diameter of driving limbs and the mass of moving platform, are discussed.

2. Elastic dynamics model of spatial parallel mechanism

When defining the telescopic rod in driving limbs as elastomers, the nonlinear dynamic model of the 4-UPS-UPU high-speed spatial parallel mechanism is established as follows.

Fig. 1Mechanism diagram of 4-UPS-UPU

Fig. 2Model of beam element

2.1. Model of rectangle beam element

The rectangle beam element is adopted to model the spatial parallel mechanism, as shown in Fig. 2. The elastic displacement and elastic corner of arbitrary point on the axis of element are given by:

1

W_{x} (x, t) = N_{A}^{T} δ,

2

W_{y} (x, t) = N_{B}^{T} δ,

3

W_{z} (x, t) = N_{C}^{T} δ,

4

ψ_{x} (x, t) = N_{D}^{T} δ,

5

ψ_{y} (x, t) = {(\frac{\partial N_{C}}{\partial x})}^{T} δ = {\dot{N}}_{C}^{T} δ,

6

ψ_{z} (x, t) = {(\frac{\partial N_{B}}{\partial x})}^{T} δ = {\dot{N}}_{B}^{T} δ,

where $W_{x} (x, t)$ , $W_{y} (x, t)$ , $W_{z} (x, t)$ are elastic displacement along $x$ , $y$ and $z$ axis, $ψ_{x} (x, t),$ $ψ_{y} (x, t)$ , $ψ_{z} (x, t)$ are elastic angle around $x$ , $y$ and $z$ axis. $δ$ is generalized coordinate vector for unit. $N_{A}$ , $N_{B}$ , $N_{C}$ and $N_{D}$ are displacement functions for unit.

The speed of arbitrary point on the axis of element is given by:

7

{\dot{W}}_{a x} (x, t) = {\dot{W}}_{r x} + {\dot{W}}_{x} (x, t),

8

{\dot{W}}_{a y} (x, t) = {\dot{W}}_{r y} (x, t) + {\dot{W}}_{y} (x, t),

9

{\dot{W}}_{a z} (x, t) = {\dot{W}}_{r z} (x, t) + {\dot{W}}_{z} (x, t),

10

{\dot{ψ}}_{a x} (x, t) = {\dot{ψ}}_{r x} (x, t) + {\dot{ψ}}_{x} (x, t),

where ${\dot{W}}_{a x} (x, t)$ , ${\dot{W}}_{a y} (x, t)$ , ${\dot{W}}_{a z} (x, t)$ are absolute velocity along $x$ , $y$ , $z$ direction, ${\dot{W}}_{r x} (x, t)$ , ${\dot{W}}_{r y} (x, t)$ , ${\dot{W}}_{r z} (x, t)$ are rigid body velocity along $x$ , $y$ , $z$ direction; ${\dot{W}}_{x} (x, t)$ , ${\dot{W}}_{y} (x, t)$ , ${\dot{W}}_{z} (x, t)$ are elastic velocity along $x$ , $y$ , $z$ direction, ${\dot{ψ}}_{a x} (x, t)$ , ${\dot{ψ}}_{r x} (x, t)$ , ${\dot{ψ}}_{x} (x, t)$ are absolute angular velocity, rigid body angular velocity and elastic velocity around $x$ direction.

2.2. Dynamic model of rectangle beam element

2.2.1. Kinetic energy of beam element

The kinetic energy of beam element without considering the rotational energy can be expressed as:

11

T = \frac{1}{2} \int_{0}^{L} m (x) [{(\frac{d W_{a x} (x, t)}{d t})}^{2} + {(\frac{d W_{a y} (x, t)}{d t})}^{2} + {(\frac{d W_{a z} (x, t)}{d t})}^{2}] d x

+ \frac{1}{2} {\int_{0}^{L} ρ I_{p} (\frac{d ψ_{a x} (x, t)}{d t})}^{2},

where $L$ is the length of the beam element, $ρ$ is the mass density of the beam element, $A$ is the cross sectional area, $m (x)$ is mass of the beam element, and $m (x) = ρ A$ , $I_{p}$ is polar moment of inertia of cross sectional area to $x$ axis.

2.2.2. Deformation energy of the beam element

The deformation energy of the beam element considering geometric nonlinearity is written as:

12

V = \frac{1}{2} E \int_{0}^{L} [I_{z} {(\frac{\partial^{2} W_{y} (x, t)}{\partial x^{2}})}^{2} + I_{y} {(\frac{\partial^{2} W_{z} (x, t)}{\partial x^{2}})}^{2}] d x + \frac{1}{2} {\int_{0}^{L} G I_{p} (\frac{\partial ψ_{a x} (x, t)}{\partial x})}^{2} d x

+ \frac{1}{2} E A {\int_{0}^{L} [\frac{\partial W_{x} (x, t)}{\partial x} + \frac{1}{2} {(\frac{\partial W_{x} (x, t)}{\partial x})}^{2} + \frac{1}{2} {(\frac{\partial W_{y} (x, t)}{\partial x})}^{2} + \frac{1}{2} {(\frac{\partial W_{z} (x, t)}{\partial x})}^{2}]}^{2} d x,

where $E$ is elastic modulus of tension and compression; $G$ is shearing modulus of elasticity; $I_{p}$ , $I_{y}$ , and $I_{z}$ are principal moment of inertia of unit cross section to $x$ , $y$ , $z$ axis.

2.2.3. Dynamic model of the beam element

Putting Eq. (11) and Eq. (12) into Lagrange equation, we can get:

13

\frac{d}{d t} (\frac{\partial T}{\partial \dot{δ}}) - \frac{\partial T}{\partial δ} + \frac{\partial V}{\partial δ} = F,

where $\dot{δ}$ is the first-order derivative of unit’s generalized coordinate vector $δ$ ; $F$ is the generalized force for unit.

From Eq. (13), the elastic dynamic model of beam element is given by:

14

M_{e} \ddot{δ} + K_{e}^{*} δ = F_{e} + P_{e} + Q_{e},

where $F_{e}$ is generalized force array of element caused by external applied load, $P_{e}$ is force array of research element caused by connecting beam element, $Q_{e} = - M \ddot{δ_{r}}$ is the rigid inertial force array, ${\ddot{δ}}_{r}$ is unit rigid body acceleration described by generalized coordinate vector, $\ddot{δ}$ is the second-order derivative of unit’s generalized coordinate vector $δ$ ; $M_{e}^{*}$ is the unit mass matrix, $K_{e}^{*}$ is the unit stiffness matrix with considering nonlinear factors.

2.3. Kinematic constraint equations

The kinematic constraint equations of spatial parallel mechanism are written as:

15

q_{S i} = [\begin{array}{l} 1 & 0 & 0 & 0 & {}^{A}{Z_{S i}} & - {}^{A}{Y_{S i}} \\ 0 & 1 & 0 & - {}^{A}{Z_{S i}} & 0 & {}^{A}{X_{S i}} \\ 0 & 0 & 1 & {}^{A}{Y_{S i}} & - {}^{A}{X_{S i}} & 0 \end{array}] q_{0},

where ${}^{A}{X_{S i}}$ , ${}^{A}{Y_{S i}}$ , ${}^{A}{Z_{S i}}$ is the coordinate of $S_{i}$ , $q_{S i} = {[\begin{array}{l} {\bar{δ}}_{i m 10} & {\bar{δ}}_{i m 11} & {\bar{δ}}_{i m 12} \end{array}]}^{T}$ is the elastic displacement vector of $S_{i}$ on $U_{i} P_{i} S_{i}$ driving limb, $q_{0}$ is the displacement of moving platform caused by elastic deformation, ${\bar{δ}}_{i m 10}$ , ${\bar{δ}}_{i m 11}$ and ${\bar{δ}}_{i m 12}$ are unit (the unit $m$ on driving limb $i$ ) generalized coordinate vector in coordinate system $\{A\}$ .

2.4. Dynamic constraint equations

The dynamic constraint equations of spatial parallel mechanism are written as:

16

[\begin{matrix} m_{0} & 0 & 0 & 0 & 0 & 0 \\ 0 & m_{0} & 0 & 0 & 0 & 0 \\ 0 & 0 & m_{0} & 0 & 0 & 0 \\ 0 & 0 & 0 & I_{x x} & I_{x y} & I_{x z} \\ 0 & 0 & 0 & I_{y x} & I_{y y} & I_{y z} \\ 0 & 0 & 0 & I_{z x} & I_{z y} & I_{z z} \end{matrix}] [\begin{array}{l} {\ddot{q}}_{1} \\ {\ddot{q}}_{2} \\ {\ddot{q}}_{3} \\ {\ddot{q}}_{4} \\ {\ddot{q}}_{5} \\ {\ddot{q}}_{6} \end{array}] = [\begin{matrix} \sum F_{i x} \\ \sum F_{i y} \\ \sum F_{i z} \\ \sum M_{i x} \\ \sum M_{i y} \\ \sum M_{i z} \end{matrix}] + [\begin{matrix} \sum F_{o x} \\ \sum F_{o y} \\ \sum F_{o z} \\ \sum M_{o x} \\ \sum M_{o y} \\ \sum M_{o z} \end{matrix}]

- [\begin{matrix} m_{0} & 0 & 0 & 0 & 0 & 0 \\ 0 & m_{0} & 0 & 0 & 0 & 0 \\ 0 & 0 & m_{0} & 0 & 0 & 0 \\ 0 & 0 & 0 & I_{x x} & I_{x y} & I_{x z} \\ 0 & 0 & 0 & I_{y x} & I_{y y} & I_{y z} \\ 0 & 0 & 0 & I_{z x} & I_{z y} & I_{z z} \end{matrix}] [\begin{matrix} {\ddot{X}}_{B} \\ {\ddot{Y}}_{B} \\ {\ddot{Z}}_{B} \\ \ddot{γ} \\ \ddot{β} \\ \ddot{α} \end{matrix}],

where $\sum F_{o x}$ , $\sum F_{o y}$ , $\sum F_{o z}$ , $\sum M_{o x}$ , $\sum M_{o y}$ , $\sum M_{o z}$ are the force and moment of moving platform caused by external applied load, $\sum F_{i x}$ , $\sum F_{i y}$ , $\sum F_{i z}$ , $\sum M_{i x}$ , $\sum M_{i y}$ , $\sum M_{i z}$ are the force and moment of moving platform caused by driving limbs, $I_{x x}$ ,…, $I_{z z}$ are the rotational inertia of moving platform, $m_{0}$ is the mass of moving platform, ${\ddot{q}}_{1}$ , ${\ddot{q}}_{2}$ ,…, ${\ddot{q}}_{6}$ are the acceleration of a moving platform, which contain 3 linear acceleration and 3 angular acceleration, ${\ddot{X}}_{B}$ , ${\ddot{Y}}_{B}$ , ${\ddot{Z}}_{B}$ are the second-order derivative of position coordinates; $\ddot{α}$ , $\ddot{β}$ and $\ddot{γ}$ are the second-order derivative of Euler-angle coordinates.

2.5. Dynamic model of the spatial parallel mechanism

When assembling the dynamic model of each element, we can get the dynamic equation of each driving limb. The dynamic model of 4-UPS-UPU High-speed spatial parallel mechanism, which is obtained by combining the dynamic equation of each driving limb, kinematic constraint equations and dynamic constraint equations, is given by:

17

M \ddot{q} + \hat{C} \dot{q} + K q = Q,

where $M = \sum_{i = 1}^{5} R_{i}^{T} M_{i} R_{i} + R_{0}^{T} M_{0} R_{0}$ is the mass matrix, $\hat{C} = \sum_{i = 1}^{5} R_{i}^{T} C_{i} R_{i}$ is the gross damping matrix, $K = \sum_{i = 1}^{5} R_{i}^{T} K_{i} R_{i}$ is the stiffness matrix, $Q = \sum_{i = 1}^{5} R_{i}^{T} P_{i} + R_{0}^{T} Q_{0}$ is the generalized force matrix, $q$ is the generalized coordinate of system, $R_{i}$ is the transfer matrix from driving limb’s generalized coordinates to system’s generalized coordinates, $R_{0}$ is the relation matrix between displacement vector of moving platform and generalized coordinates of the system.

3. The analysis method of nonlinear dynamics characteristic

The elastic dynamic equation of 4-UPS-UPU high-speed spatial parallel mechanism is a variable coefficient differential equation, the equation is commonly solved by time-discrete method. In other words, the motion time of mechanism ( $T$ ) is divided into a number of time units ( $Δ t = T / n$ ), and the elastic dynamic equation, which is regard as a second-order differential equation with constant coefficients, is solved in each time unit. Newmark method is an effective method for solving this type of variable coefficient differential equations. The basic idea of the Newmark method is described as follows.

1) The method does not require any time to meet the motion equations, just require the time in discrete points to meet the motion equations.

2) The relationship between displacement, velocity, acceleration and time are assumed in the scope of $Δ t$ .

The phase diagram, which describes the relationship between the vibration displacement and vibration speeds, is one of the commonly used method to study the nonlinear characteristics of the nonlinear dynamic system. In this paper, the nonlinear characteristics analysis of the 4-UPS-UPU high-speed parallel mechanism are achieved by combining the Newmark method with the phase diagram. That is applying Newmark method to solve nonlinear equations, getting vibration displacement and vibration velocity, and then applying phase diagram to analyze nonlinear characteristics.

When knowing the solution of dynamic differential equation in $t$ moment, the solutions of dynamic differential equation in $t + Δ t$ moment obtained by Newmark method are written as:

18

{\dot{q}}_{t + Δ t} = {\dot{q}}_{t} + [(1 - λ_{1}) {\ddot{q}}_{t} + λ_{1} {\ddot{q}}_{t + Δ t}] Δ t,

19

q_{t + Δ t} = q_{t} + {\dot{q}}_{t} Δ t + [(\frac{1}{2} - λ_{2}) {\ddot{q}}_{t} + λ_{2} {\ddot{q}}_{t + Δ t}] Δ t^{2},

where $0 \leq λ_{1} \leq 1$ , $0 \leq λ_{2} \leq 1 / 2$ , $λ_{1}$ and $λ_{2}$ are parameters related to the accuracy and stability.

In $t + Δ t$ moment, the unknown quantity $q_{t + Δ t}$ , ${\dot{q}}_{t + Δ t}$ and ${\ddot{q}}_{t + Δ t}$ meet the system elastic dynamic equation, so we can get:

20

M {\ddot{q}}_{t + Δ t} + C {\dot{q}}_{t + Δ t} + K q_{t + Δ t} = Q_{t + Δ t} .

The Eqs. (18), (19) and (20) are the basic equations of the Newmark method. The equations of $q_{t + Δ t}$ , ${\dot{q}}_{t + Δ t}$ and ${\ddot{q}}_{t + Δ t}$ can be derived by the Eqs. (18), (19) and (20). The elastic displacement vector $q_{t + Δ t}$ in $t + Δ t$ moment can be obtained by Eq. (20), and the ${\dot{q}}_{t + Δ t}$ , ${\ddot{q}}_{t + Δ t}$ in $t + Δ t$ moment can be obtained by Eqs. (18) and (19), respectively.

The main steps of nonlinear dynamics analysis are as follows:

Step 1. Calculate the total system mass matrix $M$ , stiffness matrix $K$ and damping matrix $C$ .

Step 2. Give the initial value of $q_{t = t_{0}}$ , ${\dot{q}}_{t = t_{0}}$ , ${\ddot{q}}_{t = t_{0}}$ .

Step 3. Select the time step $Δ t$ , parameters $λ_{1}$ and $λ_{2}$ , and calculate the constants of integration:

a_{0} = \frac{1}{λ_{2} Δ t^{2}}, a_{1} = \frac{λ_{1}}{λ_{2} Δ t}, a_{2} = \frac{1}{λ_{2} Δ t}, a_{3} = \frac{1}{2 λ_{2}} - 1, a_{4} = \frac{λ_{1}}{λ_{2}} - 1,

a_{5} = Δ t (\frac{λ_{1}}{2 λ_{2}} - 1), a_{6} = Δ t (1 - λ_{1}), a_{7} = λ_{1} Δ t .

Step 4. Form the effective stiffness matrix of system:

K^{*} = K + a_{0} M + a_{1} C .

Step 5. Calculate the payload of $t + Δ t$ time:

Q_{t + Δ t}^{*} = Q_{t + Δ t} + M (a_{0} q_{t} + a_{2} {\dot{q}}_{t} + a_{3} {\ddot{q}}_{t}) + C (a_{1} q_{t} + a_{4} {\dot{q}}_{t} + a_{5} {\ddot{q}}_{t}) .

Step 6. Calculate displacement of $t + Δ t$ time:

K^{*} q_{t + Δ t} = Q_{t + Δ t}^{*} .

Step 7. Calculate acceleration and velocity of $t + Δ t$ time:

{\ddot{q}}_{t + Δ t} = a_{0} (q_{t + Δ t} - q_{t}) - a_{2} {\dot{q}}_{t} - a_{3} {\ddot{q}}_{t}, {\dot{q}}_{t + Δ t} = {\dot{q}}_{t} + a_{6} {\ddot{q}}_{t} + a_{7} {\ddot{q}}_{t + Δ t} .

Step 8. Draw the curve of $q_{t + Δ t} / {\dot{q}}_{t + Δ t}$ changes with time, get the phase diagram.

4. Calculation and analysis of nonlinear characteristics for the parallel mechanism

4.1. The parameters of 4-UPS-UPU spatial parallel mechanism

In the 4-UPS-UPU spatial parallel mechanism, the distance between the first Hooke joint of the fixed platform and the center of the fixed platform is 780 mm, the other Hooke joints of the fixed platform distribute around the circle with a radius of 720 mm. The hinge points of moving platform distribute around the circle with a radius of 200 mm. The driving limbs are made of steel, mass density is 7.801×10³ kg/m³, elastic modulus $E$ = 2.1×10¹¹ Pa.

Shear modulus $G =$ 8.0×10¹⁰ Pa, Poisson’s ratio is 0.29, the mass of moving platform is 36.64 kg. The length of oscillating rod in driving limbs is 0.76 m, the length of expansion link in driving limbs is 0.88 m. The moment of inertia of the moving platform are $I_{x x} =$ 0.948 kg∙m², $I_{x y} =$ 0.81×10^-8 kg∙m², $I_{x z} =$ 0.56×10^-2 kg⋅m², $I_{y y} =$ 0.66424 kg⋅m², $I_{y z} =$ 1.166×10^-7 kg⋅m², $I_{z x} =$ 0.56×10^-2 kg⋅m², $I_{z y} =$ 1.166×10^-7 kg⋅m² and $I_{z z} =$ 0.64 kg⋅m², respectively.

4.2. Numerical example 1

The law of motion for 4-UPS-UPU spatial parallel mechanism is defined as (unit: s, m):

21

\{\begin{array}{l} X = 1.20, \\ Y = - 0.15 c o s (2 t), \\ Z = 0.15 s i n (2 t), \\ α = - \frac{π}{36} + \frac{π}{72} t, \\ β = - \frac{π}{36} + \frac{π}{72} t, (0 \leq t \leq 1) . \end{array}

According to the given motion Eq. (21) of 4-UPS-UPU high-speed spatial parallel mechanism, when the material of driving limbs are steel and aluminum, the phase diagrams are shown in Fig. 3. From Fig. 3, the motion in the $x$ -axis direction is irregular and bounded, so its movement is in the form of chaotic motion, the motion in the $y$ -axis direction is regular and cyclical, but it has diffusion in the local scope, so its movement is in the form of quasi-periodic motion, the motion in the $z$ -axis direction is also in the form of quasi-periodic motion. By and large, the trajectory of 4-UPS-UPU high-speed spatial parallel mechanism, which is very complex and has a long-term unpredictability, is a typical chaotic motion. And as shown in Fig. 3, the different material of driving limbs have a certain effect on the phase diagram.

According to the given motion Eq. (21) of 4-UPS-UPU high-speed spatial parallel mechanism, when the diameter of driving limbs are 10 mm and 40 mm, the phase diagrams are shown in Fig. 4. From Fig. 4, the diameter of driving limbs have a certain effect on the phase diagram.

According to the given motion Eq. (21) of 4-UPS-UPU high-speed spatial parallel mechanism, when the mass of the moving platform are 36 kg and 105 kg, the phase diagram are shown in Fig. 5. From Fig. 5, the motion in the $x$ -axis direction is chaotic motion, the motion in the $y$ -axis and $z$ -axis direction are periodic motion. The motion of 4-UPS-UPU high-speed spatial parallel mechanism is chaotic motion, and the different mass of the moving platform have a very obvious effect on the phase diagram.

Fig. 3The phase diagram comparison with different material of driving limbs

Fig. 4The phase diagram comparison with different diameter of driving limbs

Fig. 5The phase diagram comparison with different mass of the moving platform

4.3. Numerical example 2

The law of motion for 4-UPS-UPU spatial parallel mechanism is defined as (unit: s, m):

22

\{\begin{array}{l} X = 0.95 + 0.09 t, \\ Y = - 0.15 c o s (π t), \\ Z = 0.15 s i n (π t), \\ α = 0, \\ β = 0, (0 \leq t \leq 1) . \end{array}

According to the given motion Eq. (21) of 4-UPS-UPU high-speed spatial parallel mechanism, when the material of driving limbs are steel and aluminum, the phase diagram are shown in Fig. 6. From Fig. 6, the motion in the $x$ -axis direction is irregular and bounded, so its movement is in the form of chaotic motion, the motion in the $y$ -axis direction has the excessive trend from periodic motion to chaotic motion, so its movement is in the form of quasi-periodic motion, when the material of driving limbs is steel, the motion in the $z$ -axis direction is regular, when the material of driving limbs is aluminum, the motion in the $z$ -axis direction is chaotic motion. And as shown in Fig. 6, the different material of driving limbs have an obvious influence on the phase diagram.

According to the given motion Eq. (21) of 4-UPS-UPU high-speed spatial parallel mechanism, when the diameter of driving limbs are 10 mm and 40 mm, the phase diagrams are shown in Fig. 7. From Fig. 7, the different diameter of driving limbs have an obvious effect on the phase diagram.

According to the given motion Eq. (22) of 4-UPS-UPU high-speed spatial parallel mechanism, when the mass of the moving platform are 36 kg and 105 kg, the phase diagram are shown in Fig. 8. From Fig. 8, the motion in the $x$ -axis direction is chaotic motion, the motion in the $y$ -axis direction is quasi-periodic motion, when the mass of the moving platform is 105 kg, the motion in the $z$ -axis direction is regular, when the mass of the moving platform is 36 kg, the motion in the $z$ -axis direction is chaotic motion.

Fig. 6The phase diagram comparison with different material of driving limbs

4.4. Numerical example 3

The law of motion for 4-UPS-UPU spatial parallel mechanism is defined as (unit: s, m):

23

\{\begin{array}{l} X = 1.25, \\ Y = 0.16 t, \\ Z = 0.16 t^{2}, \\ α = 0, \\ β = 0, (0 \leq t \leq 1) . \end{array}

According to the given motion Eq. (23) of 4-UPS-UPU high-speed spatial parallel mechanism, when the material of driving limbs are steel and aluminum, the phase diagrams are shown in Fig. 9. From Fig. 9, the motion in the $x$ -axis direction is irregular and bounded, so its movement is in the form of chaotic motion, the motion in the $y$ -axis direction is regular and cyclical, but it has diffusion in the local scope, so its movement is in the form of quasi-periodic motion, the motion in the $z$ -axis direction is also in the form of quasi-periodic motion. By and large, the trajectory of 4-UPS-UPU high-speed spatial parallel mechanism, which is very complex and has a long-term unpredictability, is a typical chaotic motion. And as shown in Fig. 9, the different material of driving limbs have a certain effect on the phase diagram.

Fig. 7The phase diagram comparison with different diameter of driving limbs

Fig. 8The phase diagram comparison with different mass of the moving platform

Fig. 9The phase diagram comparison with different material of driving limbs

Fig. 10The phase diagram comparison with different diameter of driving limbs

Fig. 11The phase diagram comparison with different mass of the moving platform

According to the given motion Eq. (23) of 4-UPS-UPU high-speed spatial parallel mechanism, when the diameter of driving limbs are 10 mm and 40 mm, the phase diagrams are shown in Fig. 10. From Fig. 10, the diameter of driving limbs have a certain effect on the phase diagram.

According to the given motion Eq. (23) of 4-UPS-UPU high-speed spatial parallel mechanism, when the mass of the moving platform are 36 kg and 105 kg, the phase diagram are shown in Fig. 11. From Fig. 11, the motion in the $x$ -axis direction is chaotic motion, the motion in the $y$ -axis and $z$ -axis direction are periodic motion. The motion of 4-UPS-UPU high-speed spatial parallel mechanism is chaotic motion, and the different mass of the moving platform have a very obvious effect on the phase diagram.

By comparing the phase diagrams of the given motion Eqs. (1)-(3) of 4-UPS-UPU high-speed spatial parallel mechanism, the chaotic motion occurs in the movement, the parallel mechanism shows a typical nonlinear characteristics, and the different trajectory has an obvious effect on the phase diagram of parallel mechanism.

5. Conclusions

1) The nonlinear elastic dynamic equation of 4-UPS-UPU spatial parallel mechanism is established.

2) The dynamic equation is solved by Newmark method, the nonlinear characteristics of the parallel mechanism is analyzed by phase diagram.

3) 4-UPS-UPU high-speed parallel mechanism shows chaotic motion in the movement. The trajectory of the mechanism, the material of driving limb, the diameter of driving limb and the mass of the moving platform have a certain impact on the nonlinear characteristics of mechanism, respectively.

References

Santosha K. D., Peter E. Dynamic analysis of flexible manipulators, a literature review. Mechanism and Machine Theory, Vol. 41, 2006, p. 749-777.

Search CrossRef
Haihong I., Yang Zhiyong, Huang Tian Dynamics and elasto-dynamics optimization of a 2-DOF planar parallel pick-and-place robot with flexible links. Structural and Multidisciplinary Optimization, Vol. 38, Issue 2, 2009, p. 195-204.

Search CrossRef
Zhang X. P., Mills J. K., Cleghorn W. L. Dynamic modeling and experimental validation of a 3-PRR parallel manipulator with flexible intermediate links. Journal of Intelligent and Robotic Systems, Vol. 50, 2007, p. 323-340.

Search CrossRef
Zhang Minghui, Huang Tian Research on elastic dynamic analysis and test for Diamond mechanism. Journal of Machine Design, Vol. 21, Issue 11, 2004, p. 6-8, (in Chinese).

Search CrossRef
Zhao Yongjie, Gao Feng, Dong Xingjian, et al. Elastodynamic characteristics comparison of the 8-PSS redundant parallel manipulator and its non-redundant counterpart – the 6-PSS parallel manipulator. Mechanism and Machine Theory, Vol. 45, Issue 2, 2010, p. 291-303.

Search CrossRef
Yao Jianxin, Chen Yong Kineto-elastodynamic analysis of parallel industrial robotic manipulators. Robot, Vol. 18, Issue 6, 1996, p. 328-331, (in Chinese).

Search CrossRef
Nidal Farhat, Vicente Mata, Álvaro Page, et al. Identification of dynamic parameters of a 3-DOF RPS parallel manipulator. Mechanism and Machine Theory, Vol. 43, Issue 4, 2008, p. 1-17.

Search CrossRef
Liu Shanzeng, Yu Yueqing, Su Liying, et al. Dynamics modeling and frequency analysis of a 3-RRS flexible parallel manipulator. China Mechanical Engineering, Vol. 19, Issue 10, 2008, p. 1219-1223.

Search CrossRef
Chen Xiulong, Wei Deyong, et al. Nonlinear elastodynamic behaviour analysis of high-speed spatial parallel coordinate measuring machines. International Journal of Advanced Robotic Systems, Vol. 9, Issue 140, 2012, p. 1-6.

Search CrossRef
Shiau T. N., Tsai Y. J., Tsai M.-S. Nonlinear dynamic analysis of a parallel mechanism with consideration of joint effects. Mechanism and Machine Theory, Vol. 43, Issue 4, 2008, p. 491-505.

Search CrossRef
Li Yangmin, Xu Qingsong Modeling and performance evaluation of a flexure-based XY parallel micromanipulator. Mechanism and Machine Theory, Vol. 44, Issue 12, 2009, p. 2127-2152.

Search CrossRef

About this article

Received

21 April 2015

Accepted

28 July 2015

Published

15 November 2015

SUBJECTS

Chaos, nonlinear dynamics and applications

Keywords

spatial parallel mechanism

kineto-elastodynamics

nonlinear characteristic

chaos

Acknowledgements

This research is supported by the National Natural Science Foundation of China (Grant Nos. 51005138, 11272190, 11272167), Shandong Young Scientists Award Fund (Grant No. BS2012ZZ008), Taishan Scholarship Project of Shandong Province (No. tshw20130956), the Science Foundation of SUST (Grant No. 2011KYJQ102), the Project of Jiangsu Key Laboratory of Digital Manufacturing Technology (Grant No. HGDML-1104).

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.