Nonlinear characteristic analysis of high-speed spatial parallel mechanism

Xiulong Chen1 , Yunfeng Li2 , Yu Deng3 , Wenbin Li4 , Xiaoxia Liang5 , Qing Wang6

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

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

1Corresponding author

Journal of Vibroengineering, Vol. 17, Issue 7, 2015, p. 3974-3987.
Received 21 April 2015; received in revised form 14 July 2015; accepted 28 July 2015; published 15 November 2015

Copyright © 2015 JVE International Ltd. 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.
Creative Commons License
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.

Keywords: spatial parallel mechanism, kineto-elastodynamics, nonlinear characteristic, chaos.

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. 1. Mechanism diagram of 4-UPS-UPU

 Mechanism diagram of 4-UPS-UPU

Fig. 2. Model of beam element

 Model 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 = N C x T δ = N ˙ C T δ ,
(6)
ψ z x , t = N B x T δ = N ˙ B T δ ,

where Wx(x,t), Wy(x,t), Wz(x,t) are elastic displacement along x, y and z axis, ψxx,t,ψy(x,t), ψz(x,t) are elastic angle around x, y and z axis. δ is generalized coordinate vector for unit. NA, NB, NC and ND are displacement functions for unit.

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

(7)
W ˙ a x x , t = W ˙ r x + W ˙ x x , t ,
(8)
W ˙ a y x , t = W ˙ r y x , t + W ˙ y x , t ,
(9)
W ˙ a z x , t = W ˙ r z x , t + W ˙ z x , t ,
(10)
ψ ˙ a x x , t = ψ ˙ r x x , t + ψ ˙ x x , t ,

where W˙axx,t, W˙ayx,t, W˙azx,t are absolute velocity along x, y, z direction, W˙rxx,t, W˙ryx,t, W˙rzx,t are rigid body velocity along x, y, z direction; W˙xx,t, W˙yx,t, W˙zx,t are elastic velocity along x, y, z direction, ψ˙axx,t, ψ˙rxx,t, ψ˙xx,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 = 1 2 0 L m x d W a x x , t d t 2 + d W a y x , t d t 2 + d W a z x , t d t 2 d x
          + 1 2 0 L ρ I p 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, Ip 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 = 1 2 E 0 L I z 2 W y x , t x 2 2 + I y 2 W z x , t x 2 2 d x + 1 2 0 L G I p ψ a x x , t x 2 d x
          + 1 2 E A 0 L W x x , t x + 1 2 W x x , t x 2 + 1 2 W y x , t x 2 + 1 2 W z x , t x 2 2 d x ,

where E is elastic modulus of tension and compression; G is shearing modulus of elasticity; Ip, Iy, and Iz 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)
d d t T δ ˙ - T δ + V δ = F ,

where δ˙ 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 δ ¨ + K e * δ = F e + P e + Q e ,

where Fe is generalized force array of element caused by external applied load, Pe is force array of research element caused by connecting beam element, Qe=-Mδr¨ is the rigid inertial force array, δ¨r is unit rigid body acceleration described by generalized coordinate vector, δ¨ is the second-order derivative of unit’s generalized coordinate vector δ; Me* is the unit mass matrix, Ke* 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 = 1 0 0 0 Z S i     A - Y S i     A 0 1 0 - Z S i     A 0 X S i     A 0 0 1 Y S i     A - X S i     A 0 q 0 ,

where XSi A, YSi A, ZSi A is the coordinate of Si, qSi=δ-im10δ-im11δ-im12T is the elastic displacement vector of Si on UiPiSi driving limb, q0 is the displacement of moving platform caused by elastic deformation, δ-im10, δ-im11 and δ-im12 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)
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 q ¨ 1 q ¨ 2 q ¨ 3 q ¨ 4 q ¨ 5 q ¨ 6 = F i x F i y F i z M i x M i y M i z + F o x F o y F o z M o x M o y M o z
            - 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 X ¨ B Y ¨ B Z ¨ B γ ¨ β ¨ α ¨ ,

where Fox, Foy, Foz, Mox, Moy, Moz are the force and moment of moving platform caused by external applied load, Fix, Fiy, Fiz, Mix, Miy, Miz are the force and moment of moving platform caused by driving limbs, Ixx,…, Izz are the rotational inertia of moving platform, m0 is the mass of moving platform, q¨1, q¨2,…, q¨6 are the acceleration of a moving platform, which contain 3 linear acceleration and 3 angular acceleration, X¨B, Y¨B, Z¨B are the second-order derivative of position coordinates; α¨, β¨ and γ¨ 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 q ¨ + C ^ q ˙ + K q = Q ,

where M=i=15RiTMiRi+R0TM0R0 is the mass matrix, C^=i=15RiTCiRi is the gross damping matrix, K=i=15RiTKiRi is the stiffness matrix, Q=i=15RiTPi+R0TQ0 is the generalized force matrix, q is the generalized coordinate of system, Ri is the transfer matrix from driving limb’s generalized coordinates to system’s generalized coordinates, R0 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)
q ˙ t + Δ t = q ˙ t + ( 1 - λ 1 ) q ¨ t + λ 1 q ¨ t + Δ t Δ t ,
(19)
q t + Δ t = q t + q ˙ t Δ t + 1 2 - λ 2 q ¨ t + λ 2 q ¨ t + Δ t Δ t 2 ,

where 0λ11, 0λ21/2, λ1 and λ2 are parameters related to the accuracy and stability.

In t+Δt moment, the unknown quantity qt+Δt, q˙t+Δt and q¨t+Δt meet the system elastic dynamic equation, so we can get:

(20)
M q ¨ t + Δ t + C 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 qt+Δt, q˙t+Δt and q¨t+Δt can be derived by the Eqs. (18), (19) and (20). The elastic displacement vector qt+Δt in t+Δt moment can be obtained by Eq. (20), and the q˙t+Δt, 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 qt=t0, q˙t=t0, q¨t=t0.

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

a 0 = 1 λ 2 Δ t 2 ,       a 1 = λ 1 λ 2 Δ t ,       a 2 = 1 λ 2 Δ t ,       a 3 = 1 2 λ 2 - 1 ,       a 4 = λ 1 λ 2 - 1 ,
a 5 = Δ t λ 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 q ˙ t + a 3 q ¨ t + C a 1 q t + a 4 q ˙ t + a 5 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:

q ¨ t + Δ t = a 0 q t + Δ t - q t - a 2 q ˙ t - a 3 q ¨ t ,       q ˙ t + Δ t = q ˙ t + a 6 q ¨ t + a 7 q ¨ t + Δ t .

Step 8. Draw the curve of qt+Δt/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×103 kg/m3, elastic modulus E = 2.1×1011 Pa.

Shear modulus G= 8.0×1010 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 Ixx= 0.948 kg∙m2, Ixy=0.81×10-8 kg∙m2, Ixz= 0.56×10-2 kg⋅m2, Iyy=0.66424 kg⋅m2, Iyz=1.166×10-7 kg⋅m2, Izx= 0.56×10-2 kg⋅m2, Izy=1.166×10-7 kg⋅m2 and Izz= 0.64 kg⋅m2, respectively.

4.2. Numerical example 1

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

(21)
X = 1.20 , Y = - 0.15 c o s 2 t , Z = 0.15 s i n 2 t , α = - π 36 + π 72 t , β = - π 36 + π 72 t ,       0 t 1 .  

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. 3. The phase diagram comparison with different material of driving limbs

 The phase diagram comparison with different material of driving limbs
 The phase diagram comparison with different material of driving limbs
 The phase diagram comparison with different material of driving limbs
 The phase diagram comparison with different material of driving limbs

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

 The phase diagram comparison with different diameter of driving limbs
 The phase diagram comparison with different diameter of driving limbs
 The phase diagram comparison with different diameter of driving limbs
 The phase diagram comparison with different diameter of driving limbs

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

 The phase diagram comparison with different mass of the moving platform
 The phase diagram comparison with different mass of the moving platform
 The phase diagram comparison with different mass of the moving platform
 The 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)
X = 0.95 + 0.09 t , Y = - 0.15 c o s π t , Z = 0.15 s i n π t , α = 0 , β = 0 ,       0 t 1 .        

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. 6. The phase diagram comparison with different material of driving limbs

The phase diagram comparison with different material of driving limbs
The phase diagram comparison with different material of driving limbs
The phase diagram comparison with different material of driving limbs
The 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)
X = 1.25 , Y = 0.16 t , Z = 0.16 t 2 , α = 0 , β = 0 ,       0 t 1 .  

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. 7. The phase diagram comparison with different diameter of driving limbs

 The phase diagram comparison with different diameter of driving limbs
 The phase diagram comparison with different diameter of driving limbs
 The phase diagram comparison with different diameter of driving limbs
 The phase diagram comparison with different diameter of driving limbs

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

 The phase diagram comparison with different mass of the moving platform
 The phase diagram comparison with different mass of the moving platform
 The phase diagram comparison with different mass of the moving platform
 The phase diagram comparison with different mass of the moving platform

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

The phase diagram comparison with different material of driving limbs
The phase diagram comparison with different material of driving limbs
The phase diagram comparison with different material of driving limbs
The phase diagram comparison with different material of driving limbs

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

 The phase diagram comparison with different diameter of driving limbs
 The phase diagram comparison with different diameter of driving limbs
 The phase diagram comparison with different diameter of driving limbs
 The phase diagram comparison with different diameter of driving limbs

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

 The phase diagram comparison with different mass of the moving platform
 The phase diagram comparison with different mass of the moving platform
 The phase diagram comparison with different mass of the moving platform
 The 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.

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

References

  1. 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]
  2. 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]
  3. 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]
  4. 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]
  5. 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]
  6. 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]
  7. 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]
  8. 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]
  9. 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]
  10. 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]
  11. 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]