Next Article in Journal
Trends in the Control of Hexapod Robots: A Survey
Next Article in Special Issue
Bio-Inspired Conceptual Mechanical Design and Control of a New Human Upper Limb Exoskeleton
Previous Article in Journal
Self-Disclosure to a Robot: Only for Those Who Suffer the Most
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Kinematic Synthesis and Analysis of the RoboMech Class Parallel Manipulator with Two Grippers

by
Zhumadil Baigunchekov
1,2,
Med Amine Laribi
3,*,
Azamat Mustafa
2,* and
Abzal Kassinov
1
1
Department of Mechanics, Al–Farabi Kazakh National University, 71 Al–Farabi Ave., Almaty 050040, Kazakhstan
2
Department of Mechanical Engineering and Modeling, Satbayev University, 22a Satpaev Str., Almaty 050013, Kazakhstan
3
Department GMSC, Institut PPRIME, CNRS, Université de Poitiers, ENSMA, UPR 3346, 86962 Poitiers, France
*
Authors to whom correspondence should be addressed.
Robotics 2021, 10(3), 99; https://doi.org/10.3390/robotics10030099
Submission received: 12 June 2021 / Revised: 30 July 2021 / Accepted: 30 July 2021 / Published: 3 August 2021
(This article belongs to the Special Issue Kinematics and Robot Design IV, KaRD2021)

Abstract

:
In this paper, methods of kinematic synthesis and analysis of the RoboMech class parallel manipulator (PM) with two grippers (end effectors) are presented. This PM is formed by connecting two output objects (grippers) with a base using two passive and one negative closing kinematic chains (CKCs). A PM with two end effectors can be used for reloading operations of stamped products between two adjacent main technologies in a cold stamping line. Passive CKCs represent two serial manipulators with two degrees of freedom, and negative CKC is a three-joined link with three negative degrees of freedom. A negative CKC imposes three geometric constraints on the movements of the two output objects. Geometric parameters of the negative CKC are determined on the basis of the problems of the Chebyshev and least-square approximations. Problems of positions and analogues of velocities and accelerations of the PM with two end effectors have been solved.

1. Introduction

There are technological processes in industry where it is necessary to perform several operations simultaneously or sequentially, for example, in stamping production, in loading and unloading operations. For the simultaneous or sequential execution of several operations, it is advisable to use manipulation robots with many end effectors.
In this paper, a PM with two end effectors is synthesized that can be used to perform reloading operations from one technological equipment to another. This PM with two end effectors replaces two industrial serial robots in the existing production line of cold stamping and it belongs to the RoboMech class PM. The PM, simultaneously setting the laws of motions of the end effectors and actuators, is called the RoboMech class PM [1]. Setting the laws of motion of the actuators monotonously and uniformly but not defining by solving the inverse kinematics problem simplifies the control system and improves dynamics. Replacing two industrial robots with one RoboMech class PM with two end effectors simplifies the control system and increases the productivity and reliability of the technological line.
Since in the RoboMech class PMs simultaneously set the laws of motion of the end effectors and actuators, they work with certain structural schemes and geometric parameters of their links. The existing methods of kinematic analysis and synthesis of mechanisms and manipulators are based on the derivation of loop-closure equations and their study: in kinematic analysis, using known constant geometric parameters of links and variable generalized coordinates, variable parameters characterizing the relative movements of elements of kinematic pairs are determined, and in kinematic synthesis (dimensional or parametric synthesis) for the given positions of the input and output links, constant geometric parameters of the links are determined. Loop-closure equations are derived on the base of vector and matrix methods [2,3,4,5,6,7,8,9,10], and the theory of screws [11,12,13], which are leads to polynomials of higher degrees. Then, examining the resulting polynomials using computers, depending on the assigned tasks, the kinematic analysis or synthesis is performed. McCarthy in his papers [14,15] shows the close relationship between the kinematics, synthesis, polynomials, and computations in the 21st century. In the considered approach of kinematic analysis and synthesis of mechanisms and manipulators, it is rather difficult to obtain the polynomials; moreover, with the complication of the structures of mechanisms and manipulators, the formation of polynomials becomes more complicated and their degree increases. Performance analysis and applications of the PMs and robots are also presented in [16,17,18,19,20,21].
In this paper, kinematic synthesis of the PM with two end effectors is carried out on the basis of a modular approach [22,23], according to which PMs, regardless of their complexity, are formed by connecting the output objects (end effectors) with a base using closing kinematic chains (CKCs), which are structural modules. CKCs can be active, passive, and negative, which have positive, zero, and negative DOFs, respectively. The active and negative CKCs impose geometric constraints on the motions of the output objects, and passive CKCs do not impose geometric constraints. The representation of PMs from separate structural modules simplifies the methods of their investigation.

2. Kinematic Synthesis of the PM with Two Grippers

A PM with two end effectors can be used in a cold stamping technological line for reloading operations between two hydraulic presses [24].
Figure 1 shows a structural scheme of the PM with two end effectors in two positions.
In the first position (Figure 1a), the first gripper P 1 in position P 1 , 1 takes the workpiece after processing in the first hydraulic press for delivery to the store. At this time, the second gripper P 2 in position P 2 , 1 takes the previous workpiece processed in the first hydraulic press for delivery to the second hydraulic press for further processing.
In the second position (Figure 1b), the first gripper P 1 in position P 1 , N delivers the workpiece to the store and the second gripper P 2 in position P 2 , N delivers the previous workpiece to the second hydraulic press. The cycle is then repeated.
The considered positioning PM with two end effectors is formed by connecting two output objects (grippers P1 and P2) with a base using two passive and one negative CKC in the following sequence. First, the grippers P1 and P2 are connected to the base using passive CKCs ABC and DEF with revolute kinematic pairs, respectively, which have two degrees of freedom. Since passive CKCs A B C and D E F have two degrees of freedom, they can reproduce the given laws of motion of the output points P 1 and P 2 . Then, to form a single movable PM with two end effectors, we connect the links B C and E F of the passive CKCs A B C and D E F with the base using a negative CKC GHI with three negative degrees of freedom. Figure 2 shows a block structure of the formed PM with two end effectors.
According to the block structure (Figure 2), the parametric synthesis of the PM with two end effectors (Figure 3) consists of the parametric synthesis of two passive CKCs— A B C and D E F —and one negative CKC—GHI. The parameters of the synthesis of two passive CKCs— A B C and D E F —are X A , Z A , l A B , l B C and X D , Z D , l D E , l E F , respectively, where X A , Z A and X D , Z D are the coordinates of the fixed joints A and D in the absolute coordinate system OXYZ; l A B , l B C , l D E , l E F are the length of the links A B , B C , D E , E F . Let denote these parameters by the vectors p 1   and   p 2 , where   p 1 = [ X A , Z A , l A B , l B C ] T and p 2 = [ X D , Z D , l D E , l E F ] T .
Since the passive CKCs do not impose geometric constraints on the movements of the output points C   and   F , the vectors of the synthesis parameters p 1   and   p 2 are varied by the generator of L P τ sequence [25] to satisfy the constraints of the negative CKC IDH. In this case, the following conditions should be fulfilled:
l A B l B C ρ 1 min ,       l A B + l B C ρ 1 max ,
l D E l E F ρ 2 min ,     l D E + l E F ρ 2 max ,
where:
ρ 1 min = min X P 1 , Y P 1 l A P 1 i ,     ρ 1 max = max X P 1 , Y P 1 l A P 1 i ,
ρ 2 min = min X P 2 , Y P 2 l D P 2 i ,     ρ 2 max = max X P 2 , Y P 2 l D P 2 i .
The variable distances l A P 1 i and l C P 2 i in the Expressions (1)–(4) are determined by the equations:
l A P 1 i = X P 1 i X A 2 + Y P 1 i Y A 2 1 2 , ( i = 1 , 2 , , N ) ,
l D P 2 i = X P 2 i X D 2 + Y P 2 i Y D 2 1 2 .
Let us consider the parametric synthesis of the negative CKC GHI with three negative degrees of freedom, determined by the Chebyshev formula [26]:
W = 3 n 2 p 5 = 3 1 2 3 = 3 ,
where n is number of links, p 5 is the kinematic pairs of the fifth class.
To do this, we preliminarily determine the positions of links 2 and 4 of the passive CKCs A B C and D E F by the equations:
φ 2 i = tg 1 Z P 1 i Z B i X P 1 i X B i ,
φ 4 i = tg 1 Z P 2 i Z E i X P 2 i X E i ,
where:
X B i Z B i = X A Z A + l A B cos φ 1 i sin φ 1 i ,
X E i Z E i = X D Z D + l D E cos φ 3 i sin φ 3 i ,
φ 1 i = φ A P 1 i cos 1 l A B 2 + l A P 1 i 2 l B C 2 2 l A B l A P 1 i ,
φ 3 i = φ D P 2 i + cos 1 l D E 2 + l D P 2 i 2 l E F 2 2 l D E l D P 2 i ,
φ A P 1 i = tg 1 Z P 1 i Z A X P 1 i X A ,
φ D P 2 i = tg 1 Z P 2 i Z D X P 2 i X D .
Let us attach the coordinate systems B x 2 z 2 and   E x 4 z 4 with the links B C and E F of the passive CKCs A B C and D E F , where the axes B x 2 and E x 4 are directed along the links B C and E F , respectively (Figure 3). Then, the synthesis parameters of the negative CKC GHI are x G ( 2 ) , z G ( 2 ) , x H ( 4 ) , z H ( 4 ) ,   l G H , X I , Z I , l G I , l H I where x G ( 2 ) , z G ( 2 ) , x H ( 4 ) , z H ( 4 ) ,   X I , Z I are the coordinates of the joints G, H, I in the coordinate systems B x 2 z 2 , D x 4 z 4 , OXYZ, respectively; l G H , l G I , l H I are the lengths of the sides GH, GI, HI of the link GHI. Let us denote these synthesis parameters through the vector p 3 = [ x G ( 2 ) , z G ( 2 ) , x H ( 4 ) , z H ( 4 ) ,   l G H , X I , Z I , l G I , l H I ] .
Since the three-joined link GHI imposes three geometric constraints on the movements of the links of two passive CKCs A B C and D E F , we derive three functions of weighted differences:
Δ q 1 i = ( x H i ( 2 ) x G ( 2 ) ) 2 + ( z H i ( 2 ) z G ( 2 ) ) 2 l H G 2 ,
Δ q 2 i = ( X G i X I ) 2 + ( Z G i Z I ) 2 l G I 2 ,
Δ q 3 i = ( X H i X I ) 2 + ( Z H i Z I ) 2 l H I 2 ,
where x H i ( 2 ) and z H i ( 2 ) are the coordinates of the joint H in the local coordinate system B x 2 z 2 ; Z G i , X G i and X H i , Z H i are the coordinates of the joints G and H in the absolute coordinate system OXYZ, which are determined by the equations:
x H i ( 2 ) z H i ( 2 ) = cos φ 2 i sin φ 2 i sin φ 2 i cos φ 2 i X H i X B i Z H i Z B i ,
X G i Z G i = X B i Z B i + cos φ 2 i sin φ 2 i sin φ 2 i cos φ 2 i x G ( 2 ) z G ( 2 ) ,
X H i Z H i = X E i Z E i + cos φ 4 i sin φ 4 i sin φ 4 i cos φ 4 i x H ( 4 ) z H ( 4 ) ,
where:
X B i Z B i = X A Z A + l A B cos φ 1 i sin φ 1 i ,
X E i Z E i = X D Z D + l D E cos φ 3 i sin φ 3 i .
The geometric meanings of Functions (16)–(18) are the deviations of the coordinates of the joints H and G from circles with radiuses l H G , l G I , l H I in the relative motion of the plane E   x 4 z 4 and in the absolute motion of link 5.
After replacing the synthesis parameters of the form:
p 1 p 2 = x G ( 2 ) y G ( 2 ) , p 4 p 5 = x H ( 4 ) y H ( 4 ) ,
p 3 = 1 2 ( x G ( 2 ) 2 + z G ( 2 ) 2 + x H ( 4 ) 2 + z H ( 4 ) 2 l G H ( 2 ) ) ,
p 6 p 7 = X I Z I ,   p 8 = 1 2 ( x G ( 2 ) 2 + z G ( 2 ) 2 + X I 2 + Z I 2 l G I 2 ) ,
p 9 = 1 2 ( x H ( 4 ) 2 + z H ( 4 ) 2 + X I 2 + Z I 2 l H I 2 )
Functions (16)–(18) are expressed linearly in the following vectors of synthesis parameters p 3 ( 1 ) = [ p 1 ,   p 2 ,   p 3 ] T , p 3 ( 2 ) = [ p 4 ,   p 5 ,   p 3 ] T , p 3 ( 3 ) = [ p 6 ,   p 7 ,   p 8 ] T , p 3 ( 4 ) = [ p 1 ,   p 2 ,   p 8 ] T , p 3 ( 5 ) = [ p 6 ,   p 7 ,   p 9 ] T , p 3 ( 6 ) = [ p 4 ,   p 5 ,   p 9 ] T as:
Δ q 1 i ( k ) = 2 ( g 1 i ( k ) T p 3 ( k ) g 0 i ( k ) ) , k = 1 , 2 ,
Δ q 2 i ( k ) = 2 ( g 2 i ( k ) T p 3 ( k ) g 0 i ( k ) ) , k = 3 , 4 ,
Δ q 3 i ( k ) = 2 ( g 3 i ( k ) T p 3 ( k ) g 0 i ( k ) ) , k = 5 , 6 ,
where:
g 1 i ( 1 ) = [ Γ 1 ( φ 2 i ) 0 0 0          0 1 ] [ X E i X B i Z E i Z B i 1 ] [ Γ ( φ 4 i φ 2 i ) 0 0 0          0 1 ] [ p 4 p 5 0 ] ,
g 1 i ( 2 ) = [ Γ 1 ( φ 4 i ) 0 0 0          0 1 ] [ X E i X B i Z E i Z B i 1 ] [ Γ 1 ( φ 4 i φ 2 i ) 0 0 0          0 0 ] [ p 1 p 2 0 ] ,
g 0 i ( 1 ) = 1 2 [ ( X E i X B i ) 2 + ( Z E i Z B i ) 2 ] + [ X E i X B i , Z E i Z B i ] Γ ( φ 4 i ) p 4 p 5 ,
g 0 i ( 2 ) = 1 2 [ ( X E i X B i ) 2 + ( Z E i Z B i ) 2 ] [ X E i X B i , Z E i Z B i ] Γ ( φ 2 i ) p 1 p 2 ,
g 2 i ( 3 ) = [ X B i Z B i 1 ] [ Γ ( φ 2 i ) 0 0 0          0 1 ] [ p 1 p 2 0 ] ,
g 2 i ( 4 ) = [ Γ 1 ( φ 2 i ) 0 0 0          0 1 ] [ X B i Z B i 1 ] [ Γ ( φ 2 i ) 0 0 0          0 1 ] [ p 6 p 7 0 ] ,
g 0 i ( 3 ) = 1 2 [ X B i 2 + Z B i 2 ] + [ X B i , Z B i ] Γ ( φ 2 i ) p 1 p 2 ,
g 0 i ( 4 ) = 1 2 [ X B i 2 + Z B i 2 ] [ X B i , Z B i ] p 6 p 7 ,
g 3 i ( 5 ) = [ X E i Z E i 1 ] [ Γ ( φ 4 i ) 0 0 0          0 1 ] [ p 4 p 5 0 ] ,
g 3 i ( 6 ) = [ Γ 1 ( φ 4 i ) 0 0 0          0 1 ] [ X E i Z E i 1 ] + [ Γ ( φ 4 i ) 0 0 0          0 1 ] [ p 6 p 7 0 ] .
g 0 i ( 5 ) = 1 2 [ X E i 2 + Z E i 2 ] + [ X E i , Z E i ] Γ ( φ 4 i ) p 4 p 5 ,
g 0 i ( 6 ) = 1 2 [ X E i 2 + Z E i 2 ] [ X E i , Z E i ] p 6 p 7 .
Furthermore, the synthesis parameters of the negative CKC GHI are determined on the basis of the problems of Chebyshev and least-square approximations [16,17].

3. Kinematic Analysis of the PM with Two Grippers

In the kinematic analysis of the PM with two end effectors (Figure 4) for the given geometric parameters of the links and the input angle φ 1 i , it is necessary to determine the positions and analogues of the velocities and accelerations of the links, including the output points C and F.
The considered PM with two end effectors has the structural formula:
I ( 1 ) II ( 2 , 5 ) , II ( 3 , 4 ) ,
i.e., it contains two dyads II ( 2 , 5 ) and II ( 3 , 4 ) .
According to the structural Formula (40), first, a kinematic analysis of the dyad II ( 2 , 5 ) is carried out, and then of the dyad II ( 3 , 4 ) .

3.1. Kinematic Analysis of the PM with Two Grippers

Let us derive a vector BGI loop-closure equation:
l B G e 2 i l I G e 5 i + l ( I B ) i e ( I B ) i = 0 ,
where:
l ( I B ) i = [ ( X B i X I ) 2 + ( Z B i Z I ) 2 ] 1 2 ,
φ ( I B ) i = tg 1 Z B i Z I X B i X I ,
X B i Z B i = X A Z A + l A B cos φ 1 i sin φ 1 i .
Transfer l B G e 2 i to the right side of Equation (41) and square both sides. As a result, we obtain:
φ 5 i = φ ( I B ) i + cos 1 l ( I B ) i 2 + l I G 2 l B G 2 2 l ( I B ) i l I G .
Next, we define:
X G i Z G i = X I Z I + l I G cos φ 5 i sin φ 5 i ,
φ 2 i = tg 1 Z G i Z B i X G i X B i .
To solve the problem of the positions of the dyad II (3.4), we derive a vector DEH loop-closure equation:
l D E e 3 i + l E H e 4 i l ( D H ) i e ( D H ) i = 0 ,
where:
l ( D H ) i = X H i X D 2 + Z H i Z D 2 1 2 ,
φ ( D H ) i = tg 1 Z H i Z D X H i X D ,
X H i Z H i = X I Z I + cos φ 5 i sin φ 5 i sin φ 5 i cos φ 5 i x G ( 5 ) y G ( 5 ) .
Transfer l E H e 4 i to the right side of Equation (48) and square both sides. As a result, we obtain:
φ 3 i = φ ( D H ) i + cos 1 l D E 2 + l ( D H ) i 2 l E H 2 2 l D E l ( D H ) i .
Next, we define:
X E i Z E i = X D Z D + l D E cos φ 3 i sin φ 3 i ,
φ 4 i = tg 1 Z H i Z E i X H i X E i .
Coordinates of the output points C and F in the absolute coordinate system OXYZ are determined by the equations:
X C i Z C i = X B i Z B i + cos φ 2 i sin φ 2 i sin φ 2 i cos φ 2 i x C ( 2 ) y C ( 2 ) ,
X F i Z F i = X E i Z E i + cos φ 4 i sin φ 4 i sin φ 4 i cos φ 4 i x F ( 4 ) y F ( 4 ) .

3.2. Analogues of Velocities and Accelerations

To determine the analogues of the angular velocities of the PM with two end effectors, we derive the vector ABGI and IHED loop-closure equations:
l A B e 1 i + l B G e 2 i l I G e 5 i l A I e A I = 0
and
l I H e ( I H ) i l E H e 4 i l D E e 3 i l I D e I D = 0
and project them on the axes OX and OZ of the absolute coordinate system OXYZ
l A B cos φ 1 i + l B G cos φ 2 i l I G cos φ 5 i l A I cos φ A I = 0   l A B sin φ 1 i + l B G sin φ 2 i l I G sin φ 5 i l A I sin φ A I = 0  
and
l I H cos ( φ 5 i + α 5 ) l E H cos φ 4 i l D E cos φ 3 i l I D cos φ I D = 0   l I H sin ( φ 5 i + α 5 ) l E H sin φ 4 i l D E sin φ 3 i l I D sin φ I D = 0  
differentiate the systems of Equations (59) and (60) with respect to the generalized coordinate φ 1 i
l A B sin φ 1 i l B G sin φ 2 i φ 2 i + l I G sin φ 5 i φ 5 i = 0   l A B cos φ 1 i + l B G cos φ 2 i φ 2 i l I G cos φ 5 i φ 5 i = 0  
and
l I H sin ( φ 5 i + α 5 ) φ 5 i + l E H sin φ 4 i φ 4 i + l D E sin φ 3 i φ 3 i = 0   l I H cos ( φ 5 i + α 5 ) φ 5 i l E H cos φ 4 i φ 4 i l D E cos φ 3 i φ 3 i = 0  
From the system of Equation (61), we determine the analogues of the angular velocities φ 2 i and φ 5 i
A 1 1 u 1 = b 1
where:
A 1 = Z B i Z G i Z G i Z I X G i X B i X I X G i ,
u 1 = φ 2 i φ 5 i ,             b 1 = Z B i Z A X A X B i .
Substituting the obtained values of the angular velocity analogue φ 5 i into the system of Equation (62), from this system, we determine the angular velocities analogues φ 3 i and φ 4 i
A 2 1 u 2 = b 2 ,
where:
A 2 = Z E i Z D Z H i Z E i X D X E i X E i X H i     ,
u 2 = φ 3 i φ 4 i ,                     b 2 = ( Z H i Z I ) φ 5 i ( X I X H i ) φ 5 i .
Projections of the linear velocities analogues of the output points C and F on the axis of the absolute coordinate system OXYZ are determined by differentiating Equations (55) and (56) with respect to the generalized coordinate φ 1 i
u C i X u C i Z = u B i X u B i Z + sin φ 2 i cos φ 2 i cos φ 2 i sin φ 2 i x C ( 2 ) y C ( 2 ) φ 2 i ,
u F i X u F i Z = u E i X u E i Z + sin φ 4 i cos φ 4 i cos φ 4 i sin φ 4 i x F ( 4 ) y F ( 4 ) φ 4 i ,
where the projections of the linear velocity analogues of the joints B and E are determined by differentiating Equations (44) and (53) with respect to the generalized coordinate φ 1 i
u B i X u B i Z = l A B sin φ 1 i cos φ 1 i ,
u E i X u E i Z = l D E sin φ 3 i cos φ 3 i φ 3 i .
To determine the angular acceleration analogues of the links, we differentiate the systems of Equations (61) and (62) with respect to the generalized coordinate φ 1 i
l A B cos φ 1 i l B G cos φ 2 i φ 2 i 2 l B G sin φ 2 i φ 2 i + l I G cos φ 5 i φ 5 i 2 + l I G sin φ 5 i φ 5 i = 0 l A B sin φ 1 i l B G sin φ 2 i φ 2 i 2 + l B G cos φ 2 i φ 2 i + l I G sin φ 5 i φ 5 i 2 l I G cos φ 5 i φ 5 i = 0
and
l I H cos ( φ 5 i + α 5 ) · φ 5 i 2 l I H sin ( φ 5 i + α 5 ) · φ 5 i + l E H cos φ 4 i · φ 4 i 2 + l E H sin φ 4 i · φ 4 i + l D E cos φ 3 i · φ 3 i 2 + l D E sin φ 3 i · φ 3 i = 0 l I H sin ( φ 5 i + α 5 ) · φ 5 i 2 + l I H cos ( φ 5 i + α 5 ) · φ 5 i + l E H sin φ 4 i · φ 4 i 2 l E H cos φ 4 i · φ 4 i + l D E sin φ 3 i · φ 3 i 2 l D E cos φ 3 i · φ 3 i = 0
From the systems of Equation (69), we determine the angular accelerations analogues φ 2 i and φ 5 i
A 1 1 w 1 = b 3 ,
where:
w 1 = φ 2 i φ 5 i ,           b 3 = X B i X A + X G i X B i φ 2 i 2 + X I X G i φ 5 i 2 Z B i Z A + Z G i Z A φ 2 i 2 + Z I Z G i φ 5 i 2 .
Substituting the obtained values of the angular velocity analogues φ 2 i and φ 5 i into the system of Equation (70), from this system, we determine the angular velocities analogues φ 3 i and φ 4 i
A 2 1 w 2 = b 4 ,
where:
w 2 = φ 3 i φ 4 i , b 4 = X H i X I φ 5 i 2 + X E i X H i φ 4 i 2 + X D X E i φ 3 i 2 Z H i Z I φ 5 i 2 + Z E i Z H i φ 4 i 2 + Z D Z E i φ 3 i 2 .
Projections of the linear velocities analogues of the output points C and F on the axis of the absolute coordinate system OXYZ are determined by differentiating Equations (65) and (66) with respect to the generalized coordinate φ 1 i
w C i X w C i Z = w B i X w B i Z + cos φ 2 i sin φ 2 i sin φ 2 i cos φ 2 i x C ( 2 ) z C ( 2 ) φ 2 i 2 + sin φ 2 i cos φ 2 i cos φ 2 i sin φ 2 i x C ( 2 ) z C ( 2 ) φ 2 i ,
w F i X w F i Z = w E i X w E i Z + cos φ 4 i sin φ 4 i sin φ 4 i cos φ 4 i x F ( 4 ) z F ( 4 ) φ 4 i 2 + sin φ 4 i cos φ 4 i cos φ 4 i sin φ 4 i x F ( 2 ) z F ( 2 ) φ 4 i ,
where the projections of the linear velocity analogues of the joints B and E are determined by differentiating Equations (67) and (68) with respect to the generalized coordinate φ 1 i
w B i X w B i Z = l A B cos φ 1 i sin φ 1 i ,
w E i X w E i Z = l D E cos φ 3 i sin φ 3 i φ 3 i 2 + sin φ 3 i cos φ 3 i φ 3 i .

4. Numerical Results

Table 1 shows N = 11 positions of the grippers P 1 and P 2 of the PM with two end effectors.
Table 2, Table 3 and Table 4 show the obtained values of the synthesis parameters of the two passive CKCs ABC, DEF and negative CKC GHI, respectively.
3D CAD model of the synthesized PM with two grippers is shown in Figure 5.
Table 5 shows the obtained values of the coordinates X P 1 , i , Y P 1 , i , X P 2 , i , Y P 2 , i and projections of analogues of linear velocities u P 1 , i X , u P 1 , i Z , u P 2 , i X , u P 2 , i Z and linear accelerations w P 1 , i X , w P 1 , i Z , w P 2 , i X , w P 2 , i Z of the grippers of the PM with two end effectors.
Positions and modules of the velocities and acceleration analogues of the synthesized PM grippers P 1 and P 2 are also presented with the graphical plots in Figure 6, Figure 7 and Figure 8.

5. Conclusions

Kinematic synthesis and analysis of the PM with two end effectors have been carried out. In the kinematic synthesis according to the given laws of motions (or positions) of two end effectors, the structural scheme, and geometric parameters of links of the synthesized PM are determined. The structural scheme of this PM is formed by connecting two output objects (end effectors) and a base using three CKCs: two passive and one negative CKC. Passive and negative CKCs are structural modules from which the PM is formed. Passive CKCs are two movable serial manipulators, and the negative CKCs is a three-jointed link. Serial manipulators (passive CKCs) do not impose geometric constraints on the movement of the output objects, and the three-jointed link (negative CKC) imposes three geometric constraints. Therefore, the geometric parameters of the links of the negative CKCs are determined, and the geometric parameters of the links of the passive CKCs are varied depending on the imposed geometric constraints of the negative CKC. Kinematic synthesis of the negative CKC was carried out on the basis of the Chebyshev and least-square approximations. Since the structure of the synthesized PM consists of two dyads, the position analysis is solved analytically. Analogues of angular velocities and accelerations are determined from two systems of linear equations obtained by differentiating the loop-closure equations with respect to the generalized coordinate.

Author Contributions

Z.B. established the methods of structural and kinematic synthesis of the RoboMech class parallel manipulator (PM) with two grippers. M.A.L. established the methods of kinematic analysis of this PM. A.M. created the programs in the MatLab for kinematic synthesis of the PM and obtained numerical results. A.K. created the programs in the MatLab for kinematic analysis of the PM and obtained numerical results. All authors contributed equally to the evaluation of the results and writing of this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research is founded by the Science Committee of the Ministry of Education and Science of the Republic of Kazakhstan (Grant No AP08857522).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baigunchekov, Z.; Kalimoldayev, M.; Ibrayev, S.; Izmambetov, M.; Baigunchekov, T.; Naurushev, B.; Aisa, N. Parallel Manipulator of a Class RoboMech. Lect. Notes Electr. Eng. 2017, 408, 547–557. [Google Scholar]
  2. Paul, B. Kinematics and Dynamics of Planar Machinery; Prentice-Hall: Hoboken, NJ, USA, 1979. [Google Scholar]
  3. Hartenberg, R.S.; Denavit, J. Kinematic Synthesis of Linkages; McGraw-Hill: New York, NY, USA, 1964. [Google Scholar]
  4. Rodhavan, M.; Roth, B. Solving Polynomial Systems for the Kinematic Analysis and Synthesis of Mechanisms and Robot Manipulators. J. Mech. Des. 1995, 117, 71–79. [Google Scholar]
  5. Erdman, A.G.; Sandor, G.N.; Kota, S. Mechanism Design: Analysis and Synthesis Lecture Notes; McGill University: Montreal, QC, Canada, 2001. [Google Scholar]
  6. Uicker, J.J.; Pennack, G.R.; Shigley, J.E. Theory of Machines and Mechanisms, 4th ed.; Oxford University Press: New York, NY, USA, 2001. [Google Scholar]
  7. Hernandez, A.; Petya, V. Position Analysis of Planar Mechanisms with R-pairs Using Geometrical-Iterative Method. Mech. Mach. Theory 2004, 39, 50–61. [Google Scholar] [CrossRef]
  8. Pennock, G.R.; Israr, A. Kinematic Analysis and Synthesis of an Adjustable Six-Bar Linkage. Mech. Mach. Theory 2009, 44, 306–323. [Google Scholar] [CrossRef]
  9. McCarthy, J.M.; Soh, G.M. Geometric Design of Linkages, 2nd ed.; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  10. Angeles, J.; Bai, S.H. Kinematic Synthesis, Lecture Notes; McGill University: Montreal, QC, Canada, 2016. [Google Scholar]
  11. Duffy, J. Analysis of Mechanisms and Robot Manipulator; Edward Arnold: London, UK, 1980. [Google Scholar]
  12. Tsai, L.W. Robot Analysis: The Mechanics of Serial and Parallel Manipulators; John Wiley and Sons, Inc.: New York, NY, USA, 1999. [Google Scholar]
  13. Glazunov, V. Design of Decoupled Parallel Manipulators by Means of the Theory of Screws. Mech. Mach. Theory 2010, 45, 209–230. [Google Scholar] [CrossRef]
  14. McCarthy, J.M. Kinematics, Polynomials, and Computers-A Brief History. J. Mech. Robot. 2011, 3, 010201-1. [Google Scholar] [CrossRef] [Green Version]
  15. McCarthy, J.M. 21st Kinematics Synthesis, Compliance, and Tensegrity. J. Mech. Robot. 2011, 3, 020202-1. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, D.; Wang, L.; Wu, J. Physics-Based Mechatronics Modeling and Application of an Industrial-Grade Parallel Tool Head. Mech. Syst. Signal Process. 2021, 148, 107158. [Google Scholar] [CrossRef]
  17. Bourbonnais, F.; Bigras, P.; Bonev, I.A. Minimum-Time Trajectory Planning and Control of a Pick-and-Place Five-Bar Parallel Robot. IEEE/ASME Trans. Mechatron. 2014, 20, 740–749. [Google Scholar] [CrossRef]
  18. Zhang, Z.; Wang, L.; Shao, Z. Improving the Kinematic Performance of a Planar 3-RRR Parallel Manipulator Through Actuation Mode Conversion. Mech. Mach. Theory 2018, 130, 86–108. [Google Scholar] [CrossRef]
  19. Carabin, G.; Scalera, L.; Wongratanaphisan, T.; Vidoni, R. An Energy-Efficient Approach for 3D Printing with a Linear Delta Robot Equipped with Optimal Springs. Robot. Comput. Integr. Manuf. 2021, 67, 102045. [Google Scholar] [CrossRef]
  20. Scalera, L.; Carabin, G.; Vidoni, R.; Wongratanaphisan, T. Energy Efficiency in a 4-DOF Parallel Robot Featuring Compliant Elements . Int. J. Mech. Control 2019, 20, 49–57. [Google Scholar]
  21. Hu, B.; Shi, D.; Xie, T.; Hu, B.; Ye, N. Kinematically Identical Manipulators Derivation for the 2-RPU+ UPR Parallel Manipulator and Their Constraint Performance Comparison. J. Mech. Robot. 2020, 6, 061011. [Google Scholar] [CrossRef]
  22. Baigunchekov, Z.; Ibrayev, S.; Izmambetov, M.; Mustafa, A.; Sagyntay, M. Synthesis of Reconfigurable Positioning Parallel Manipulator of a Class RoboMech. In Proceedings of the 2018 International Conference on Reconfigurable Mechanisms and Robots, ReMAR, Delft, The Netherlands, 20–22 June 2018. [Google Scholar]
  23. Baigunchekov, Z.; Izmambetov, M.; Zhumasheva, Z.; Baigunchekov, T.; Mustafa, A. Parallel Manipulator of a Class RoboMech for Generation of Horizontal Trajectories Family. Mech. Mach. Sci. 2019, 73, 1395–1402. [Google Scholar]
  24. Grabchenko, I.P.; Klepnikov, V.B.; Dobroskok, V.L.; Kryzhni, G.K. Introduction to Mechatronics; HPI: Harkov, Ukraine, 2017. [Google Scholar]
  25. Sobol, I.M.; Statnikov, R.B. Choice of Optimal Parameters in Problems with Many Criteria; "Nauka": Moscow, Russian, 1981. [Google Scholar]
  26. Artobolevskiy, I.I. Theory of Mechanisms and Machines; "Nauka": Moscow, Russian, 1988. [Google Scholar]
Figure 1. PM with two end effectors in two positions: (a) the first position (b) the second position.
Figure 1. PM with two end effectors in two positions: (a) the first position (b) the second position.
Robotics 10 00099 g001
Figure 2. Block structure of the PM with two end effectors.
Figure 2. Block structure of the PM with two end effectors.
Robotics 10 00099 g002
Figure 3. PM with two end effectors in the first position.
Figure 3. PM with two end effectors in the first position.
Robotics 10 00099 g003
Figure 4. PM with two end effectors.
Figure 4. PM with two end effectors.
Robotics 10 00099 g004
Figure 5. 3D CAD model of the PM with two grippers.
Figure 5. 3D CAD model of the PM with two grippers.
Robotics 10 00099 g005
Figure 6. Graphics of the grippers P 1 and P 2 positions.
Figure 6. Graphics of the grippers P 1 and P 2 positions.
Robotics 10 00099 g006
Figure 7. Graphics of the grippers P 1 and P 2 linear velocities analogues.
Figure 7. Graphics of the grippers P 1 and P 2 linear velocities analogues.
Robotics 10 00099 g007
Figure 8. Graphics of the grippers P 1 and P 2 linear accelerations analogues.
Figure 8. Graphics of the grippers P 1 and P 2 linear accelerations analogues.
Robotics 10 00099 g008
Table 1. Positions of the grippers P 1 and P2.
Table 1. Positions of the grippers P 1 and P2.
i1234567891011
X P 1 , i , mm 04.69.213.818.4023.027.6032.2036.8041.4046.3882
Y P 1 , i , mm 0−6.2512−11.0478−14.8398−17.8092−20.1433−21.8922−23.1166−23.8528−24.1065−23.8384
X P 2 , i ,   mm 46.388253.593859.835965.336770.434475.022379.189682.945386.473889.647692.7766
Y P 2 , i ,   mm −23.8384−24.0439−23.2688−21.8398−19.8480−17.4431−14.6785−11.6242−8.1590−4.42280
Table 2. Synthesis parameters of the CKC ABC.
Table 2. Synthesis parameters of the CKC ABC.
X A ,   mm Z A ,   mm l A B ,   mm l B C ,   mm
14.04−57.1146.530943.0
Table 3. Synthesis parameters of the CKC DEF.
Table 3. Synthesis parameters of the CKC DEF.
X D ,   mm Z D ,   mm l D E ,   mm l E F ,   mm
78.74−57.1146.530943.0
Table 4. Synthesis parameters of the CKC GHI.
Table 4. Synthesis parameters of the CKC GHI.
x G 2 ,   mm y G 2 ,   mm x H 4 ,   mm y H 4 ,   mm l G H ,   mm X I ,   mm Z I ,   mm l G I ,   mm l H I ,   mm
6.60729.9607−6.6072−9.960745.046.3957.1150.050.0
Table 5. Positions and analogues of linear velocities and accelerations of the grippers.
Table 5. Positions and analogues of linear velocities and accelerations of the grippers.
i1234567891011
X P 1 , i ,   mm 04.69.213.818.4023.027.6032.2036.8041.4046.3882
Z P 1 , i ,   mm 0−6.2512−11.0478−14.8398−17.8092−20.1433−21.8922−23.1166−23.8528−24.1065−23.8384
X P 2 , i ,   mm 46.388253.593859.835965.336770.434475.022379.189682.945386.473889.647692.7766
Z P 2 , i ,   mm −23.8384−24.0439−23.2688−21.8398−19.8480−17.4431−14.6785−11.6242−8.1590−4.42280
u P 1 , i X ,   mm 0.84400.52040.2432−0.026−0.276−0.495−0.693−0.8729−1.036−1.184−1.325
u P 1 , i Z ,   mm 0.98300.96270.94890.93620.92360.91100.89770.88290.86620.84650.8211
u P 2 , i X ,   mm 1.07331.03801.01730.99980.98450.97100.95830.94590.93350.92120.9094
u P 2 , i Z ,   mm −0.1071−0.0733−0.0650−0.0515−0.0399−0.0294−0.0189−0.00750.00620.02400.0515
w P 1 , i X ,   mm −0.1096−0.1066−0.0868−0.1056−0.1964−0.3570−0.5768−0.8396−1.1265−1.4157−1.6956
w P 1 , i Z ,   mm −0.8688−0.2504−0.02830.0221−0.0380−0.1620−0.3271−0.5223−0.7421−0.9843−1.2642
w P 2 , i X ,   mm −0.1116−0.0885−0.0968−0.1112−0.1353−0.1676−0.2100−0.2661−0.3420−0.4486−0.6182
w P 2 , i Z ,   mm 0.64010.47100.38030.30690.24160.17710.10380.01319−0.1059−0.2713−0.5320
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Baigunchekov, Z.; Laribi, M.A.; Mustafa, A.; Kassinov, A. Kinematic Synthesis and Analysis of the RoboMech Class Parallel Manipulator with Two Grippers. Robotics 2021, 10, 99. https://doi.org/10.3390/robotics10030099

AMA Style

Baigunchekov Z, Laribi MA, Mustafa A, Kassinov A. Kinematic Synthesis and Analysis of the RoboMech Class Parallel Manipulator with Two Grippers. Robotics. 2021; 10(3):99. https://doi.org/10.3390/robotics10030099

Chicago/Turabian Style

Baigunchekov, Zhumadil, Med Amine Laribi, Azamat Mustafa, and Abzal Kassinov. 2021. "Kinematic Synthesis and Analysis of the RoboMech Class Parallel Manipulator with Two Grippers" Robotics 10, no. 3: 99. https://doi.org/10.3390/robotics10030099

APA Style

Baigunchekov, Z., Laribi, M. A., Mustafa, A., & Kassinov, A. (2021). Kinematic Synthesis and Analysis of the RoboMech Class Parallel Manipulator with Two Grippers. Robotics, 10(3), 99. https://doi.org/10.3390/robotics10030099

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop