Dynamics of A Novel Robotic Leg Based On The Peauc
Dynamics of A Novel Robotic Leg Based On The Peauc
Abstract
This article deals with the kinematics and dynamics of a novel leg based on the Peaucellier–Lipkin mechanism, which is
better known as the straight path tracer. The basic Peaucellier–Lipkin linkage with 1 degree of freedom was transformed
into a more skillful mechanism, through the addition of 4 more degrees of freedom. The resulting 5-degree-of-freedom
leg enables the walking machine to move along paths that are straight lines and/or concave or convex curves. Three
degrees of freedom transform the leg in relation to a reachable center of rotation that the machine walks around. Once
the leg is transformed, the remaining 2 degrees of freedom position the foot at a desirable Cartesian point during the
transfer or support phase. We analyzed the direct and inverse kinematics developed for the leg when the foot describes
a straight line and found some interesting relationships among the motion parameters. The dynamic model equations of
motion for the leg were derived from the Lagrangian dynamic formulation to calculate the required torques during a par-
ticular transfer phase.
Keywords
Peaucellier–Lipkin mechanism, transfer phase, straight line path, Lagrangian dynamics, robotic leg, kinematics, dynamics
Creative Commons CC-BY: This article is distributed under the terms of the Creative Commons Attribution 3.0 License
(http://www.creativecommons.org/licenses/by/3.0/) which permits any use, reproduction and distribution of the work without
further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/
open-access-at-sage).
2 Advances in Mechanical Engineering
Shkolnik11 1 Yes No
Garrec12 1 Yes No
NPU 5 Yes Yes
LAB 0.15
LBC 0.15
LAD = LAE 0.42
LCD = LCE = LDF = LEF 0.20
0 1 0 1
N xF LAF cos a
N
F = @ yF A = @ LAF sin a A ð14Þ
zF 0
where
yC 1
a = tan ð15Þ
xC
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
LAC = x2C + y2C ð16Þ
1 L2AD + L2AC L2CD
b = cos ð17Þ
2LAD LAC
g=a+b ð18Þ
d=a b ð19Þ
Figure 4. Description of the joints for the robotic leg. LAF = 2LAD cos b LAC ð20Þ
The angular accelerations for the links are calculated describes the straight lines LAB and LBC , they are con-
as derivatives of the angular velocities with respect to stant and equal. It is possible to calculate b and LAC by
€ 2 is calculated as equation (30)
time (t). For example, v means of equations (41) and (42)
d d2 1 L2AD + LAF (t)2 L2DF
€ 2 = v_ 2 = 2 v2
v ð30Þ b(t) = cos ð41Þ
dt dt 2LAD LAF (t)
and a similar process is applied to obtain the other LAC = 2LAD cos b LAF (t) ð42Þ
accelerations.
To calculate the linear velocities of the links in the Finally, using equation (43), the value of u(t) is deter-
centers of mass, it is necessary to describe the center of mined, which in turn determines the position, velocity,
mass positions from CM1 to CM8 with respect to the and acceleration for each link in the direct kinematics
frame fN g, for example, ( N rCM1 1 )T =
N T
TAE ( AErCM1 1 ) . A similar process is used to deter- LAC (t) cos a(t)
u(t) = tan1 ð43Þ
mine the absolute description of the mass centers LAC (t) sin a(t) LAB
N
rCM2 , N rCM3 , N rCM4 , N rCM5 , N rCM6 , N rCM7 , and N rCM8 .
The linear velocities in the mass centers of the links are
Trajectory generation
described in equations (31)–(38), where d_ = d^ _ zN ,
_ = u^
u _ zN , s_ = s^ _ zN , m _ = m^ _ zN , 2_ = 2^
_ zN , v _ 2 = v_ 2^zN , The main purpose of the foot trajectory is to determine
_ 3 = v_ 3^zN , v
v _ 5 = v_ 5^zN , v_ 6 = v_ 6^zN , and u_ = u^_ zN the required joint variable u(t) that would be used by
the actuator to allow the leg to follow the required
v1 = d_ N rCM1 ð31Þ path. During the transfer phase describing straight
lines, the leg uses one DOF. To determine the straight
v2 = d_ N E
+v
_ 2 N rCM2 N E ð32Þ line trajectory of the foot, a parametric equation (44) is
used, where the initial and final point positions are
v3 = d_ N E
+v
_ 3 N rCM3 N E ð33Þ o = ð 0 0:35 0 ÞT
P and Pf = ð 0 0:35 0 ÞT ,
respectively, and the function p(t) 2 ½ 0, 1 represents
d_ N E
v4 = + (d_ + s N E)
_ ) (N F ð34Þ
the trajectory profile. For this particular case, p(t) is
v5 = g +v
_ N D _ 5 N rCM5 N D ð35Þ the polynomial parametric profile described in equation
(45)19
v6 = g +v
_ N D _ 6 N rCM6 N D ð36Þ
o + P
T (t) = P f P
o p(t) ð44Þ
v7 = g _ N rCM7 ð37Þ
!
v8 = u_ N N rCM8 N B
ð38Þ p(t) =
70 3 315 4 546 420 120
t 4 t + 5 t5 6 t6 + 7 t7
tf3 tf tf tf tf
Due to the symmetry of the mechanism, it is !
expected that the angular displacement utilized to 32 160 288 224 64
+ vmax 2 t3 + 3 t4 4 t5 + 5 t6 6 t7
describe links AD, AE, CD, CE, DF, and EF would be tf tf tf tf tf
dependent on position, velocity, and angular
ð45Þ
acceleration.
where
Inverse kinematics c
vmax = ð46Þ
This section describes how u is determined as a func- tf
tion of time, when the foot moves along a straight line
on a known path. Once the position The index c 2 ½ 1:458334, 2:1875 produces natural
N
F = ( xF (t), yF (t), zF (t) ) is known, the variables movements in the mechanism. For this particular case,
LAF and a are determined in equations (39) and (40), the index used is c = 1:875 and the resulting polynomial
respectively parametric profile becomes a fifth degree profile. The
total time for the step is tf = 5 s.
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
LAF (t) = xF (t)2 + yF (t)2 ð39Þ
Dynamics of the propulsion unit
yF (t)
a(t) = tan1 ð40Þ Due to the simplicity of the mechanism, the leg motion
xF (t)
equations are derived by applying the typical
The lengths of links LAD , LAE , LCD , LCE , LDF , and LEF Lagrangian dynamic formulation for mechanical sys-
are constant in time. In this study only, when joint F tems,20 but the results are similar to those using
6 Advances in Mechanical Engineering
1 1
Ti = mi v2cmi + Icmi v2i ð50Þ
formulations employed by Hollerbach21 or Kane22 that 2 2
are normally used in robotic mechanical systems. In the
where mi is the mass of the ith link; vcmi is the linear
dynamic analysis, the mass of each link is considered to
velocity of the center of mass of the ith link; Icmi is the
be concentrated at the point called the center of mass,
mass moment of inertia of the ith link, measured in the
described in Table 3 and shown in Figure 5. The mass
center of mass; and vi is the angular velocity of the ith
moment of inertia is calculated at the center of mass of
link, measured in the inertial base.
each link, with respect to its local frame. It is calculated
Finally, to determine the total torque required to
using a CAD model and the values are shown in
move the leg by joint B, it is necessary to add all the
Table 4.
torques according to fT = f1 + f2 + + fi , i = 8.
To obtain the dynamic model, the following consid-
erations are assumed:
Results and discussion
1. Friction in the joints is negligible due to the use
of bearings. The physical parameters employed in the simulation
2. The links of the mechanism move on the plane are described in Tables 2–4. The positions of the links
formed by ^xN ^yN . and their centers of mass when the foot follows a linear
3. The densities of the links are constant. path using a polynomial profile can be seen in Figure 6.
The time used to follow this trajectory (equation (44))
is tf = 5 s with k = 40 discretizations. Progress in rela-
tion to time is calculated with equation (51)
Development of the dynamic equations
tf
The dynamic analysis of the leg is developed for the tn = n , n = 0, 1, . . . , k ð51Þ
k
transfer phase, in which the forces of the ground reac-
tion to the foot are equal to zero. This is because there Figure 6 shows the two positions of the leg when
is no terrain interaction and the positions of the mass n = 8 and n = 32, the points on the straight line, are
Núñez-Altamirano et al. 7
Figure 10. Angular acceleration of the links. Figure 11. Linear velocity of the links.
s_ = m;
_ f_ = 2_ ð55Þ
v_ 5 = v_ 2 ; v_ 6 = v_ 3 ð56Þ
d_ n = g_ kn ð57Þ
Figure 12. Torque required by the links.
Figure 10 shows the angular acceleration of the leg
links calculated according to the process described by
equation (30). Because the acceleration is the variation linear velocity of the links decreases when they are near
with the time of the velocity, consequently it is possible the frame fN g.
to generate interesting relationships between the angu- The previously described parameters are necessary
lar accelerations. The relationships of the relative and for the dynamic model so that the torque can be calcu-
absolute accelerations are presented in equations (58) lated and the leg moved. According to these para-
and (59), respectively, and are confirmed if equations meters, it is possible to estimate the link that will
(55) and (56) are derived with respect to time. At least a require more torque, assuming the mass as the more
double reflection relationship in time between the angu- important parameter, but also considering position,
lar acceleration €
d and g
€ is represented in equation (60) velocity, and acceleration. It is postulated that the foot
requires more torque because it has the biggest mass.
€ = €
s m; € = €n
u ð58Þ Another important link to consider is the DF link, due
to its mass, velocity, and acceleration.
€5 = v
v € 2; €6 = v
v €3 ð59Þ
Torque is calculated independently for each link
€dn = g
€kn ð60Þ using equation (35). The values of the required torque
for the links, with the exception of the foot, are illu-
The linear velocity of each link was calculated at its strated in Figure 12. It shows that link DF( f5 ) is the
center of mass, shown in Figure 11. The equations used one that is expected to require the highest torque.
to calculate the velocity for the foot are described in Figure 13 shows the torque that is required for mov-
equations (31)–(38). ing the foot f4 and the total torque fT that is necessary
According to the leg configuration, the center of for moving the leg. One can see how the torque required
mass located in the foot is the furthest from the inertial for the foot is more than half of the total torque needed
frame fN g, and consequently its linear velocity is the by the leg, as postulated above, due to the foot’s large
highest in the leg, as shown in Figure 11, where the mass.
Núñez-Altamirano et al. 9
References
1. Peaucellier C. Note sur une question de géométrie de
compas. Nouv Ann Math 1873; 2: 71–78.
2. Dijksman E. On the history of focal mechanisms and
their derivatives. In: Proceedings of the international
symposium on history of machines and mechanisms, Cas-
sino, 12–15 May 2004, pp.303–314. The Netherlands:
Springer, Kluwer Academic Publishers. doi: 10.1007/1-
4020-2204-2_24
3. Luquin F. Máquinas de vapor, mecanismos y P.L. Che-
byshev. In: Un paseo por la geometria - Matematicalia:
revista digital de divulgación matemática (ed Raúl Ibáñez
Torres y Marta Macho Stadler), 2005. ISSN: 1699-7700,
pp.3–19. http://www.matematicalia.net/index.php?option=
Figure 13. Foot torque and total torque in the leg. com_content&task=view&id=20&Itemid=5
4. Gao X-S, Zhu C-C, Chou S-C, et al. Automated genera-
Conclusion tion of Kempe linkages for algebraic curves and surfaces.
Mech Mach Theory 2001; 36: 1019–1033.
Our article describes the kinematics and dynamics for a 5. Shigley JE and Uicker JJ. Theory of machines and
novel robotic leg based on the PL mechanism when the mechanisms. Theory of machines and mechanisms. In:
foot follows a linear path on a transfer phase. International edition (ed. Julienne V. Brown and James
Interesting relationships in the positions, velocities, W. Bradley). USA: McGrawHill, 1981, pp.23–25.
6. Kempe AB. How to draw a straight line. London: Nature
and angular accelerations (complementary, supplemen-
Series, 1877.
tary, reflective, and double reflective) between the
7. Dijksman EA. A fourfold generalization of Peaucellier’s
mechanism links were found that are useful for decreas- inversion cell. Meccanica 1996; 31: 407–420.
ing the number of kinematic equations and thus reduc- 8. Buckley J. An investigation into straight-line mechan-
ing the quantity of calculations. isms: optimizing the Peaucellier mechanism (faculty advi-
The use of the Lagrangian equations of motion for sor. Dr. Ming Huang). In: Proceedings of the National
the mass of each link of the mechanism enabled the Conference on Undergraduate Research (NCUR) 2011,
quantification of the required torque for each link, and Ithaca College, New York, 31 March–2 April 2011.
the proposal of improvements in mechanism design, 9. Chalfin NL. Phonograph pickup arm. Patent 2977126 A,
such as the reduction of mass in the foot. USA, 1961.
10. Kunzler PA, Schmitt P, Poblano RDV, et al. Wheel
Most robotic legs are based on open or serial kine-
embedded suspension. Patent 7938210 B2, USA, 2011.
matic chains. Due to the fact that the PL robotic leg is 11. Shkolnik N. Walking apparatus. Patent 4462476 A,
based on a parallel linkage with eight bars, it distributes USA, 1984.
forces and torques in a manner superior to that of the 12. Garrec P. Device for transmitting movement between a
open or serial kinematic chains. Given that the robotic solid and a member, in particular for a robot able to be
leg is based on the PL mechanism, which traces exact moved on legs. Patent 5219410 A, USA, 1993.
paths as linear or circular curves, it requires the employ- 13. Roy SS and Pratihar DK. Effects of turning gait para-
ment of only two actuators to position the foot. These meters on energy consumption and stability of a six-
are the main advantages to using the apparently com- legged walking robot. Robot Auton Syst 2012; 60: 72–82.
plex PL robotic leg, compared with other mechanical 14. Roy SS and Pratihar DK. Soft computing-based
approaches to predict energy consumption and stability
leg architectures.
margin of six-legged robots moving on gradient terrains.
Appl Intell 2012; 37: 31–46.
Declaration of conflicting interests 15. Roy SS and Pratihar DK. Kinematics, dynamics and
power consumption analyses for turning motion of a six-
The author(s) declared no potential conflicts of interest with
legged robot. J Intell Robot Syst 2013; 74: 663–688.
respect to the research, authorship, and/or publication of this
16. Roy SS and Pratihar DK. Dynamic modeling, stability
article.
and energy consumption analysis of a realistic six-
legged walking robot. Robot CIM: Int Manuf 2013; 29:
Funding
400–416.
The author(s) disclosed receipt of the following financial sup- 17. Wang Z, Ding X, Rovetta A, et al. Mobility analysis of
port for the research, authorship, and/or publication of this the typical gait of a radial symmetrical six-legged robot.
article: The authors acknowledge to Consejo Nacional Mechatronics 2011; 21: 1133–1146.
de Ciencia y Tecnologı́a (CONACyT) for the scholarship 18. Juárez Campos I, Núñez Altamirano DA, Márquez Pérez
support. L, et al. Unidad de propulsión reconfigurable para
10 Advances in Mechanical Engineering
máquinas caminantes. Patent application MX2014008184- 20. Wells DA. Theory and problems of Lagrangian dynamics.
A. IMPI (MEX), México, 2014. http://lp.espacenet.com/ New York: McGraw-Hill, 1967.
searchResults?locale=es_LP&AP=MX20140008184&D 21. Hollerbach JM. A recursive Lagrangian formulation of
B=lp.espacenet.com&ST=advanced maniputator dynamics and a comparative study of
19. Júarez Campos I. Sobre las Regiones Geome´tricas Articu- dynamics formulation complexity. IEEE T Syst Man
lares y su Aplicación en el Mejoramiento del Desempeño Cyb 1980; 10: 730–736.
de un Manipulador Mediante Estrategias Evolutivas. PhD 22. Kane TR. Lagrange’s equations for a rigid body. J Appl
Thesis, National Autonomous University of Mexico Mech 1973; 40: 1118–1119.
(UNAM), Mexico, 2002. pp. 160–167. 23. Thomas GB. Thomas’ calculus. 12th ed. Boston, MA:
Pearson, 2009, p. 118.