Dynamics analysis of a novel two-translation parallel mechanism

Hongzhou Wang1 , Mengmeng Song2 , Mengwei Li3 , Ru Yang4 , Zhen Jiang5

1, 4School of Mechatronic Engineering and Automation, Shanghai University, Shanghai, People’s Republic of China

2, 3College of Information, Mechanical and Electrical Engineering, Ningde Normal University, Ningde, People’s Republic of China

5College of Mechanical and Electrical Engineering, Fujian Agriculture and Forestry University, Fuzhou, People’s Republic of China

2, 5Corresponding authors

Vibroengineering PROCEDIA, Vol. 33, 2020, p. 193-198. https://doi.org/10.21595/vp.2020.21686
Received 7 September 2020; received in revised form 14 November 2020; accepted 22 September 2020; published 19 October 2020

Copyright © 2020 Hongzhou Wang, et al. 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
Table of Contents Download PDF Acknowledgements References
Cite this article
Views 12
Reads 5
Downloads 190
CrossRef Citations 0
Abstract.

In this paper, we take a novel two-translation parallel mechanism based on variable length parallelogram structure as the research object, and carry out the research on the dynamic characteristics of the mechanism. Based on the principle of virtual work, the dynamic model of the new two-translation parallel mechanism is established. By numerically solving the dynamic model, the dynamic behavior of the new mechanism with impact characteristics is found, which provides a theoretical basis for optimizing this new mechanism.

Dynamics analysis of a novel two-translation parallel mechanism

Highlights
  • A dynamic model of the new two-translation parallel mechanism is established.
  • Dynamic behavior of the new two-translation parallel mechanism with impact characteristics is revealed.
  • A prototype of the new mechanism was manufactured based on the dynamic characteristics.

Keywords: parallel mechanism, dynamics analysis, dynamic modeling.

1. Introduction

In industry, in order to complete the handling and picking and placing of workpieces, only two translational degrees of freedom are required to meet the work requirements [1-3]. Liu proposed a planar two-degree-of-freedom parallel mechanism with a constant direction [1], in which each branch chain contains a planar four-bar mechanism. Kim also proposed a new type of two-translation parallel mechanism composed of two groups of PPa structures (P: moving pair, Pa: parallelogram closed-loop sub-chain) [2]. Professor Huang Tian of Tianjin University used the properties of the parallelogram structure to propose a fully articulated parallel mechanism called Diamond with two translational degrees of freedom. It is a two-dimensional form of the Delta mechanism. The overall performance of the robot is almost comparable to that of Delta. The mechanism is comparable to the single-axis feed mechanism and can be connected in series to form a hybrid mechanism with three translational degrees of freedom [4]. However, the two translational parallel mechanism proposed above needs to connect the two branches in a symmetrical embracing manner, so it has the disadvantages of a large lateral area and small vertical movement space. Based on this, this paper proposes a new two-translation parallel mechanism based on a variable-length parallelogram structure, as shown in Fig. 1(a), and takes the new structure as the research object to carry out the dynamic analysis of the mechanism, and then demonstrate the feasibility of the mechanism.

2. Dynamic modeling of novel two-translation parallel mechanism

Since the robot system is a non-linear and complex dynamic system, its dynamic characteristics must be analyzed in the robot dynamic real-time control system. Rigid body dynamics analysis is mainly to carry out inverse dynamics analysis of rigid body, that is, given the geometric and inertial parameters and motion laws of the mechanical system, the driving force (torque) required to realize this kind of motion is solved. For this reason, this section mainly uses the principle of virtual work to construct the inverse dynamics model of the new two-translation parallel mechanism.

Fig. 1. Schematic diagram of kinetic analysis of novel two-translation parallel mechanism

 Schematic diagram of kinetic analysis of novel two-translation parallel mechanism

a) The overall schematic diagram

 Schematic diagram of kinetic analysis of novel two-translation parallel mechanism

b) Schematic diagram of branch chain

In order to facilitate analysis, as shown in Fig. 1, a coordinate system o-xyz is established at the center of the fixed platform A1A2, which is marked as A, where the y axis is parallel to A1A2 and the x axis is parallel to the axis of the rotating pair A1A2, and the z axis is determined by the right-hand rule. The coordinate system p-uvw established on the moving platform B1B2 is recorded as B, where the v axis is parallel to B1B2, the u axis and the rotating pair B1B2 are in the same direction, and the w axis follows the right-hand rule.

First, suppose the external inertial generalized forces acting on the moving platform and each branch chain are F^P and F^i (i=1,2). Since each branch chain in the new two-translation parallel mechanism is composed of sliders and rails, the external inertial force acting on the branch chain can be divided into two parts F^i=F^i1+F^i2, where the subscript i1 represents the variables related to the sliders, and the subscript i2 represents variables related to retractable rails. Then the Jacobian matrix of the branch is determined. Finally, the driving force Γ=[F1,τ1]T is considered, therefore, the principle of virtual work applied to the new two-translation parallel mechanism can be expressed as:

(1)
J T Γ + F ^ P + i = 1 2 J i 1 T F ^ i 1 + J i 2 T F ^ i 2 = 0 .

As shown in Fig. 1, l1=l2=l, Q=[l,β]T, which is the rate of change of the input joint, χ=[  AyP,  AzP]T, which represents the speed of the output point P of the moving platform. Therefore, the closed loop equation of each branch can be written as:

(2)
    A P +     A R B     B B i =     A l s ^ i +     A A i .

Derive Eqs. (2) with respect to time, the outcome is:

(3)
  A v P +   A R ˙ B B B i +   A R B B B ˙ i =   A l ˙ s ^ i +   A l s ^ ˙ i +   A A ˙ i .

Because   BBi=0,   AAi=0,   ARB=0, there is:

(4)
s ^ ˙ i = β ˙ k ^ × s ^ i ,

where s^i is the unit direction vector of the branch chain i, and k^ is the unit vector of the direction of the axis of the rotating pair. In order to eliminate β, both sides of Eq. (3) are multiplied by s^i at the same time to obtain the following branch chain linear velocity expression:

(5)
l ˙ = v P s ^ i .

In order to eliminate l˙, the two sides of Eq. (3) are cross-multiplied by s^i at the same time to obtain the angular velocity of the branch chain:

(6)
β ˙ k ^ = 1 l s ^ i × v P .

So we can write Q˙=Jχ˙.

Derive Eq. (3) with respect to time, the outcome is:

(7)
a P = l ¨ s ^ i + 2 l ˙ β ˙ k ^ × s ^ i + l β ¨ k ^ × s ^ i - l β ˙ 2 s ^ i .

In order to eliminate β¨ to obtain the expression of l¨, both sides of Eq. (7) are multiplied by s^i at the same time, the outcome is:

(8)
l ¨ = a P s ^ i + l β ˙ 2 .

In order to eliminate l¨ to obtain the expression of β¨, both sides of the Eq. (7) are cross-multiplied by s^i at the same time, and after finishing, Eq. (8) can be expressed as:

(9)
l ¨ = a P s ^ i + l β ˙ 2 .

Since the branch chain can be decomposed into two independent parts, the slider and the guide rail, their masses are represented by mi1 and mi2 respectively, as shown in Fig. 1(b). The position vectors of the two centroids can be determined by the following formula:

(10)
p i 1 = A i + c i 1 s ^ i ,
(11)
p i 2 = A i + l - c i 2 s ^ i .

After deriving the above two formulas, the linear velocity expressions at the centroids of the two parts of the branch chain are obtained:

(12)
v ^ i 1 = c i 1 β ˙ k ^ × s ^ i ,
(13)
v ^ i 2 = l - c i 2 β ˙ k ^ × s ^ i + l ˙ s ^ i .

Derive Eq. (12) and Eq. (13), we can get the expression of linear acceleration velocity at the centroid of each branch:

(14)
a ^ i 1 = c ˙ i 1 β ˙ k ^ × s ^ i + c i 1 β ¨ k ^ × s ^ i + c i 1 β ˙ k ^ ˙ × s ^ i + c i 1 β ˙ k ^ × s ^ ˙ i ,
(15)
a ^ i 2 = l ˙ - c ˙ i 2 β ˙ k ^ × s ^ i + l - c i 2 β ¨ k ^ × s ^ i + l - c i 2 β ˙ k ^ ˙ × s ^ i
            + l - c i 2 β ˙ k ^ × s ^ ˙ i + l ¨ s ^ i + l ˙ s ^ ˙ i .

Since the end moving platform only has the degree of freedom to move along the y axis and the z axis, the external inertial torsion force acting on the moving platform can be expressed as:

(16)
F ^ P = f d y + M g y ^ - M y ¨ P f d z + M g z ^ - M z ¨ P ,

where M is the mass of the moving platform, fdy and fdz are the perturbation torsion force acting on the moving platform along the y axis and the z axis. M(gy^) and M(gz^) are the gravity along the y axis and the z axis respectively.

External inertial torsion force acting on the branch chain is:

(17)
F ^ i 1 = f ^ i 1 n ^ i 1 = m i 1 g ^ - a ^ i 1 - A I c i 1 β ¨ k ^ - β ˙ k ^ × A I c i 1 β ˙ k ^ ,           F ^ i 2 = f ^ i 2 n ^ i 2 = m i 2 g ^ - a ^ i 2 - A I c i 2 β ¨ k ^ - β ˙ k ^ × A I c i 2 β ˙ k ^ .

Among them, g^=[0,0,-g]T. mi1 and mi2 are the masses of the slider and the guide rail respectively. Ici1 represents the inertia matrix of the slider whose coordinate origin is on the point ci1 and parallel to A. And also Ici2 represents the inertia matrix of the guide rail whose coordinate origin is on the point ci2 and parallel to A. Because the coordinate system of the point ci1 and the point ci2 is moving relative to the coordinate system fixed on the fixed platform, the inertia matrix can be obtained by the following transformation:

(18)
  A I c i 1 =   A R A i A i I c i 1 A R A i T ,             A I c i 2 =   A R A i A i I c i 2 A R A i T .

It is assumed here that the branches are completely symmetric,   AiIci1 and   AiIci2 are also diagonal matrixs.

From the linear velocity and angular velocity of the slider and the guide rail, the Jacobian matrix Ji1 and Ji2 of the branch chain can be obtained:

(19)
χ i 1 = J i 1 χ ˙ ,             χ i 2 = J i 2 χ ˙ .

Among them, χi1=[v^i1,β˙k^]T, χi2=[v^i2,β˙]T. In the same coordinate system [s^i,n^i,k^]T, the Jacobian matrix Ji1 and Ji2 of the branch can be derived.

Since the speed and acceleration of the branch chain 1 and the branch chain 2 in the two-translation parallel mechanism are the same, the dynamic equation can be derived after the result obtained by the above calculation is put into Eq. (1), and finally the driving force and driving torque of the two-translation parallel mechanism can be obtained.

3. Dynamic simulation of a novel two-translation parallel mechanism

The structural parameters of the novel two-translation parallel mechanism are shown in Table 1.

Table 1. Structural parameters of the novel two-translation parallel mechanism

The mass of the moving platform M (kg)
The mass of the slider m1 (kg)
The mass of the rail m2 (kg)
The distance of the slider centroid and the hinge point of fixed platform ci1 (mm)
The distance of the rail centroid and the hinge point of moving platform ci2 (mm)
0.947
2.042
1.059
59.700
77.520

The inertia matrix in the novel two-translation parallel mechanism (kg⋅mm2):

  A i I c i 1 = 35872.354 0 0 0 35229.282 0 0 0 1032.678 ,
  A i I c i 2 = 1.082 × 10 7 0 0 0 1.787 × 10 5 0 0 0 1.071 × 10 7 .

Assume that the motion law of the moving platform of the novel two-translation parallel mechanism is:

(20)
y P = 150 c o s 4 π t z P = 100 s i n 4 π t + 50 .

Fig. 2. Driving force and driving torque

 Driving force and driving torque

a) Driving force

 Driving force and driving torque

b) Driving torque

Fig. 3. The prototype of 2T parallel mechanism

 The prototype of 2T parallel mechanism

a) Frontal partial view

 The prototype of 2T parallel mechanism

b) Back partial view

 The prototype of 2T parallel mechanism

c) Overall view

In addition, the velocity and acceleration of the moving platform can be obtained by calculating the first derivative and the second derivative of the time with the trajectory function. By writing the above dynamic model into MATLAB, the input force and input torque curves of the driver can be obtained, respectively as shown in Fig. 2(a) and 2(b). It can be seen that under this trajectory, an impact will occur in the movement of the novel two-translation parallel mechanism, which will affect the normal operation of the mechanism. For this reason, the structure and motion trajectory of the mechanism need to be further optimized to improve the dynamic characteristics of the mechanism. The prototype of the optimized novel two-translation parallel mechanism is shown in Fig. 3.

4. Conclusions

In this work, a new type of two translational parallel mechanism based on variable length parallelogram structure is taken as the research object, and the dynamic model of the new mechanism is established by using the principle of virtual work. Through numerical solution, the dynamic characteristics such as driving force and driving torque of the new two translational parallel mechanism are obtained, which lays the foundation for optimizing the structure of the mechanism. Based on the dynamic characteristics, a prototype of a new mechanism was manufactured.

Acknowledgements

This paper was supported by the following research projects: by the Special Project of Ningde Normal University in 2018 (Grant No. 2018ZX401, Grant No. 2019ZX401, Grant No. 2018Q102) and Research project for Yong, Middle-aged Teacher in Fujian Province (Grant No. JT180597). These supports are gratefully acknowledged.

References

  1. Liu X. J., Wang Q. M., Wang J. Kinematics, dynamics and dimensional synthesis of a novel 2-DoF translational manipulator. Journal of Intelligent and Robotic Systems, Vol. 41, Issue 4, 2005, p. 205-224. [Publisher]
  2. Pham V. B. N., Kim H. S. Dynamics analysis of a 2-DOF planar translational parallel manipulator. Journal of the Korean Society for Production and Manufacturing Systems, Vol. 22, Issue 2, 2013, p. 185-191. [Search CrossRef]
  3. Yang Y., Peng Y., Pu H., et al. Design of 2-degrees-of-freedom (DOF) planar translational mechanisms with parallel linear motion elements for an automatic docking device. Mechanism and Machine Theory, Vol. 121, 2018, p. 398-424. [Publisher]
  4. Huang T., Liu S., Mei J., et al. Optimal design of a 2-DOF pick-and-place parallel robot using dynamic performance indices and angular constraints. Mechanism and Machine Theory, Vol. 70, Issue 6, 2013, p. 246-253. [Publisher]