Next Article in Journal
Communication-Efficient Wireless Traffic Prediction with Federated Learning
Previous Article in Journal
Study of Flexibility Transformation in Thermal Power Enterprises under Multi-Factor Drivers: Application of Complex-Network Evolutionary Game Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Symbolic Computation of the Lie Algebra se(3) of the Euclidean Group SE(3): An Application to the Infinitesimal Kinematics of Robot Manipulators

by
Jaime Gallardo-Alvarado
1,
Mario A. Garcia-Murillo
2,*,
Juan Manuel Tabares-Martinez
1 and
X. Yamile Sandoval-Castro
3
1
Department of Mechanical Engineering, National Technological Institute of Mexico, Celaya Campus, Celaya 38010, Mexico
2
Department of Mechanical Engineering, DICIS, University of Guanajuato, Salamanca 36885, Mexico
3
School of Engineering and Sciences, Tecnologico de Monterrey, Monterrey 64849, Mexico
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(16), 2538; https://doi.org/10.3390/math12162538
Submission received: 4 June 2024 / Revised: 12 August 2024 / Accepted: 14 August 2024 / Published: 16 August 2024
(This article belongs to the Special Issue Application of Mathematical Analysis and Theory to Robotics)

Abstract

:
This paper reports an application of the Lie algebra s e ( 3 ) of the Euclidean group S E ( 3 ) , which is isomorphic to the theory of screws in the velocity and acceleration analyses of serial manipulators. The symbolic computation of the infinitesimal kinematics allows one to obtain algebraic expressions related to the kinematic characteristics of the end effector of the serial manipulator, while in the case of complex manipulators, numerical computations are preferred owing to the emergence of long terms. The algorithm presented enables the symbolic computation of the velocity and acceleration characteristics of the end effector in serial manipulators in order to allow the compact and efficient modeling of velocity and acceleration analyses of both parallel and serial robotic manipulators. Unlike other algebras, these procedures can be extended without significant effort to higher-order analyses such as the jerk and jounce.

1. Introduction

Over the past few decades, interest in both academia and industry for developing robotic manipulators has increased considerably due to their potential applications. The repetitive and precise applications of robotic manipulators have been complemented by activities in which robots must perform unplanned tasks, but which need to be executed in real time. Thus, nowadays robotic manipulators require highly reliable mathematical resources for their performance.
The study of rigid bodies in motion is one of the most influential topics in modern robotics; e.g., industrial robot arms are a classical reference. A group can be defined as the abstraction of the properties of symmetry operations. In kinematics, the group of rigid transformations is called the Special Euclidean group S E ( 3 ) , which is investigated in depth by means of the Lie algebra s e ( 3 ) . Perhaps the most comprehensive contributions to the rigid body motion related to the Lie algebra s e ( 3 ) can be found in the theory of screws developed by Ball [1] at the beginning of the century. After the pioneering work of Ball, it took several decades for the screw theory to be revalorized thanks to the contributions of Hunt [2] and Phillips [3]. Afterward, the isomorphism between the screw theory and the Lie algebra s e ( 3 ) of the Special Euclidean group S E ( 3 ) led to remarkable productivity within modern differential geometry and Lie group theory; see, for instance, [4,5]. Wu and Carricato [6] reported a method to identify Lie triple screw systems of the Lie algebra s e ( 3 ) based on algebraic Lie group theory and descriptive screw theory and their relationship with the exponential motion manifolds of the Lie triple screw systems in dual quaternion representation. It is known that a curve in S E ( 3 ) whose tangent vector the left or right invariant representation is a curve in s e ( 3 ) can be used to describe the general continuous motion of a rigid body. Zhang et al. [7] developed a matrix method based on the geometric properties of a pair of conjugate axodes equivalent to a curve in s e ( 3 ) in algebras to achieve this result. Sanchez-Garcia et al. [8] investigated the mobility of parallel manipulators by means of screw theory in which the intersection of subalgebras of Lie algebra, s e ( 3 ) , of the special Euclidean group, S E ( 3 ) , associated with the limbs of the parallel manipulator plays a central role. Chen et al. [9] introduced a multisymplectic Lie algebra variational integrator for simulating the dynamics of flexible multibody systems.
Mathematically, a Euclidean group is the characteristic group of isometries of a Euclidean space E n where n denotes the dimension of the space and its transformations preserve the Euclidean distance between any pair of points of the physical space. The configuration of the group depends on the dimension n of the space and is commonly notated as E ( n ) . Note that the Euclidean group E ( n ) consists of all translations, rotations, and reflections of E n , as well as arbitrary finite combinations of these transformations. The Euclidean group can be considered the symmetry group of the space itself and contains the group of symmetries of any subset of that space. A Euclidean isometry may be direct or indirect, depending on whether it preserves the parity of the subsets. In that regard, direct Euclidean isometries can be grouped into the special Euclidean group S E ( n ) , a subgroup of E ( n ) whose elements are called rigid motions or Euclidean motions. It is worth emphasizing that they comprise arbitrary combinations of translations and rotations, but reflections are excluded. The special group S E ( n ) is among the oldest, at least in the cases of dimensions 2 and 3, being implicitly studied long before the concept of the group was envisioned. The contribution is focused on the special Euclidean group S E ( 3 ) . Euclidean groups in addition to being topological groups are also catalogued as Lie groups, after Marius Sophus Lie, and therefore, the concepts of infinitesimal calculus can be immediately adapted to this configuration.
Screw theory is perhaps just another algebra for the kinematic analysis of robotic manipulators; however, its effectiveness is well proven in geometric interpretation. For example, the quaternion algebra works very well in velocity analysis, but its extension to acceleration analysis encounters difficulties that its followers have not solved satisfactorily, so it is often necessary to resort to devices that escape mathematical rigor, as is the case with exponential methods. Screw theory, as a number of investigations have shown, has largely outperformed higher-order analysis. Therefore, the contribution of this possible contribution will certainly meet the expectations of users of three-dimensional kinematics of robotic manipulators.
In mathematics and computer science, symbolic computation, also known as algebraic computation, is a scientific area that refers to the study and development of algorithms and software for manipulating mathematical expressions derived from diverse disciplines of knowledge and other mathematical objects in general. Symbolic computation is a highly reliable strategy in the validation of mathematical procedures that attempt to model the behavior of natural phenomena in which the use of mathematical formulas plays a preponderant role. In addition, and beyond modeling physical phenomena, in this case, the kinematics of robotic systems, algebraic computation is widely used in areas such as control engineering, in model-based strategies that require solving Coupled Nonlinear Systems of Variable-Order Fractional PDEs, especially if higher-order kinematic variables need to be controlled [10,11].
Within kinematic analysis, one of the fields where computational methods are usually employed is in position analysis, especially in parallel manipulators. Nearly forty years ago, Freudenstein referred to the challenge of solving polynomial expressions from the inverse kinematics problem of the 7R spatial linkage as the “Mount Everest of kinematics”. Since then, there has been a continuous increase in the publication of specialized papers introducing new kinematic analysis strategies or refining existing methods, leading to a significant re-evaluation of kinematics. Despite the considerable progress made, challenges persist in addressing these analyses. Research in this area can be categorized into three main approaches: (i) iterative numerical solutions, (ii) the implementation of auxiliary sensors, and (iii) polynomial solutions. Iterative methods are straightforward but face limitations such as incomplete solutions, convergence issues, and the need for good initial estimates. The second approach requires a physical prototype. Polynomial methods, on the other hand, mitigate these challenges and provide a reliable way to determine the feasible poses of a parallel manipulators for a given set of prescribed generalized coordinates.
Moreover, with the advent of symbolic software, it is possible to generate the symbolic equations of practically any engineering problem. In this paper, we present symbolic procedures generated with the software Maple ver. 2022 associated with the operations of the Lie algebra s e ( 3 ) of the Euclidean group S E ( 3 ) . These procedures are general and can be applied without major difficulties in the kinematic analysis of robotic manipulators.
The aim of the paper is not to discuss or provide new insights into the Lie algebra s e ( 3 ) of the Euclidean Group S E ( 3 ) . Rather, the focus and contribution of the paper are as follows: (a) to provide symbolic procedures that allow the compact and efficient modeling of velocity and acceleration analyses of both parallel and serial robotic manipulators; (b) to show that the Lie algebra structure allows us to easily provide computational algorithm operations and analyses compared to other algebras, such as the quaternion algebra, which lead to unintuitive methods that are not always generalizable or that are difficult to implement for complex robotic structures; and (c) to show that these results can easily be extended to higher-order kinematic analyses (third and fourth time derivatives) without significant effort.
The rest of this contribution is organized as follows: in Section 2, the fundamentals of the Lie algebra of the Euclidean group are presented and related to fundamental concepts of the kinematics of multibody systems, such as the velocity state and the reduced acceleration state of a rigid body and open kinematic chains. In Section 3, some operations of Lie Algebra implemented in Maple used in kinematic analysis are reported. Then, in Section 4, applications and numerical examples of this contribution to the kinematic analysis of representative open kinematic chains are presented. Finally, Section 5 concludes this paper.

2. Theory

2.1. Plücker Coordinates

The homogeneous coordinates of a segment line in space are known as Plücker coordinates and were introduced in the 19th century by Julius Plücker [12,13]. The Plücker coordinates are also known as Grassmann coordinates and provide a logical beginning to introduce the basics of the theory of screws.
In physics, a rigid body or object is a solid body in which deformation is zero or so small that it can be disregarded. Thus, the distance between two arbitrary points of a rigid body remains constant even though the body undergoes external forces or moments. The pose of a rigid body refers to the information required to determine the position of the center of mass of the rigid body, as well as its orientation or attitude. In the general case, the pose of a rigid body is defined by six parameters: three for the position and three for the orientation. The number of parameters considered in the pose of the rigid body is called the degrees of freedom of the rigid body. Considering that, in a line segment, the rotation around an axis along the line itself does not affect the pose of the line, then four parameters are sufficient to determine the pose of a line segment.
Figure 1 shows a segment line l passing through point p. According to the reference frame O _ X Y Z , the line l can be characterized entirely by introducing the unit vector l ^ and the position vector p of point p. Then, it is possible to define a vector L R 6 formed with the concatenation of two vectors l ^ and m as follows
L = l ^ m
where m = l ^ × q is the moment vector of the unit vector l ^ about the origin O of the reference frame O _ X Y Z . Note that vector q is normal to vector l ^ . Naturally, any point of the line l may be used to compute the moment vector m . For example, let us consider that q = p + λ l ^ , where λ R . Then, it follows that
m = q × l ^ = ( p + λ l ^ ) × l ^ = p × l ^ + λ ( l ^ × l ^ ) = p × l ^
The parameters involved in Equation (2) are called the Plücker coordinates of the line segment. It is straightforward to demonstrate that Equation (2) has two redundancies: (i) l ^ is a unit vector, which implies that l ^ = 1 ; (ii) it must satisfy that l ^ · m = 0 . Then, a line in Plücker coordinates can be characterized by only four parameters.
Figure 2 shows a line l passing through point p. Let us consider an arbitrary point o. The moment vector m o produced by the line l around the point o may be computed as
m o = ( p o ) × l ^ = p × l ^ o × l ^ = m o × l ^
Hence
m o = m o × l ^
Furthermore, according to point q, the vector m o may be computed as
m o = ( q o ) × l ^
Premultiplying both sides of Equation (5) by l ^ , it follows that
l ^ × m o = l ^ × [ ( q o ) × l ^ ]
However, the expulsion property of the triple cross product leads to
l ^ × [ ( q o ) × l ^ ] = l ^ · l ^ ( q o ) [ l ^ · ( q o ) ] l ^ = q o
Hence,
q = o + l ^ × m o
The reciprocal product between two lines is one of the highlights of Plücker coordinates. Let l 1 and l 2 be two lines whose Plücker coordinates are established as follows:
L 1 = l ^ 1 m 1 , L 2 = l ^ 2 m 2
Moreover, according to Equation (4), it follows that
m 1 = m 2 p 1 × l ^ 2
The application of the inner product between the unit vector l ^ 1 and the right side of Equation (10) yields
l ^ 1 · ( m 2 p 1 × l ^ 2 ) = l ^ 1 · m 2 + l ^ 2 · m 1
where the term l ^ 1 · ( m 2 p 1 × l ^ 2 ) is known as the reciprocal product between two lines, and it is notated as L 1 L 2 . That is to say,
L 1 L 2 = l ^ 1 · m 2 + l ^ 2 · m 1
The geometry of lines found in the screw theory is a formidable ally for its mathematical formulation. For example, a screw, which is notated as $, can be represented in Plücker coordinates as follows:
$ = s ^ s O
where the unit normalized vector s ^ is known as the primal part of the screw, and it occurs along the axis of the screw, which is called the Mozzi’s axis. Meanwhile, s O is the moment produced by s ^ around the reference pole O. The vector s O is called the dual part of the screw, and it is computed as
s O = h s ^ + r Q / O × s ^
where h is the pitch of the screw, while r Q / O is a vector beginning at the reference pole O and ending in arbitrary point Q of Mozzi’s axis. Mozzi’s axis is the Instantaneous Screw Axis (ISA). In the rotational displacement, the pitch h of the screw vanishes, and therefore, in this case, we have that
$ = s ^ r Q / O × s ^
On the other hand, in the translational displacement, the pitch h of the screw goes to infinity, and therefore, it follows that
$ = lim h s ^ h s ^ + r Q / O × s ^ = lim h s ^ lim h ( h s ^ + r Q / O × s ^ ) = 0 s ^

2.2. Chasles’ Theorem

Theorem 1
(Euler’s Theorem). “By arbitrarily rotating a sphere around its center, it is always possible to find a diameter whose position is equal to the initial position in spite of the rotation.” [14]
One element of the Euclidean group that has attracted the most attention is that corresponding to rotational motion. Based on Euler’s pioneering work, it is possible to establish that any spatial displacement of a rigid body that maintains a constant or fixed point produces an axis with constant orientation. Moreover, Euler elegantly concluded that a composition of arbitrary rotations about the same point is equivalent to a single rotation about a new axis, which is designated as Euler’s pole. Motions satisfying these conditions have an algebraic structure which is designated as the special orthogonal group S O ( 3 ) of the Euclidean group S E ( 3 ) and whose dimension is 3. Mozzi extends Euler’s theorem, concluding that the displacement of a rigid body is equivalent to a rotation around an axis passing through the center of mass of the body followed by a translation in an arbitrary direction. As a prelude to the celebrated theorem of Chasles, Mozzi was able to glimpse that in the displacement of the rigid body, there are points of the rigid body that translate with the same intensity along axes parallel to the axis of rotation. In 1904, Whittaker was able to formally state that “A rotation about any axis is equivalent to a rotation by the same angle about any axis parallel to it, together with a simple translation in a direction parallel to the axis itself”.
Theorem 2
(Chasles’ Theorem). “The general motion of a rigid body in space can be expressed as a translational motion along a helical axis, called the Mozzi axis, preceded by a rotational motion around an axis parallel to the Mozzi axis.” [15,16]
Figure 3 shows a rigid body m in general motion with respect to another body j or reference frame O _ X Y Z . Assuming that u ^ is a unit vector along the ISA then the angular velocity vector   j ω m of body m with respect to body j may be expressed as   j ω m =   j ω m u ^ . On the other hand, let us consider two points P and Q attached to the ISA, where, in general, P Q . Then, the velocity vectors of points P and Q may be related as
v P = v Q +   j ω m × r P / Q
where r P / Q is the position vector of point P with respect to point Q. Taking into account that r P / Q = r P / Q u ^ , it follows that
v P = v Q +   j ω m × r P / Q = v Q + (   j ω m u ^ ) × ( r P / Q u ^ ) = v Q
By mathematical induction, it is concluded that all the points on the ISA have the same velocity vector. Other interesting properties of Chasles’ theorem are the following:
  • The velocity component along the Mozzi axis is the same for all points on the body m.
  • The velocity component of any point of the body m along an axis perpendicular to the ISA is null.
  • The velocity vector of any point of the body m lies in a plane defined by two vectors; the first one is parallel to the ISA, while the second one is normal to the ISA and to the position vector of the point of interest and an arbitrary point on the ISA.

2.3. The Lie Algebra s e ( 3 ) of the Euclidean Group S E ( 3 )

The special Euclidean group S E ( 3 ) is a Lie group composed of a set of rigid body motions and has an associated Lie algebra s e ( 3 ) . The real space S E ( 3 ) is a six-dimensional manifold where the dimension is the number of degrees of freedom of a rigid body freely moving in space. In physics, Lie groups emerge as symmetry groups of physical systems, and their Lie algebras may be considered as infinitesimal symmetry motions. In this section, the operations and properties of the Lie algebra s e ( 3 ) of the Euclidean group S E ( 3 ) are briefly explained from the point of view of the kinematic state of a rigid body. To enumerate the operations of the Lie s e ( 3 ) of the Euclidean group S E ( 3 ) it must be considered that it is isomorphic to the theory of screws [1] and the motor algebra [17].
Let us consider that a rigid body m is in general motion with respect to another body j or reference frame O _ X Y Z . Assuming that   j ω m is the angular velocity vector of body m as observed from body j, while v P is the velocity vector of an arbitrary point P of body m, then it is possible to define the velocity state   j V P m of body m with respect to body j as the concatenation of vectors   j ω m and v P as follows:
  j V P m =   j ω m v P
The angular velocity vector   j ω m = P (   j V P m ) is referred to as the primal component of the velocity state   j V P m , while the velocity vector v P = D (   j V P m ) is referred to as the dual part of   j V P m . It is interesting to note that while the primal component of the velocity state is the same in any part of body m, the dual part depends on the point, in turn, of body m. Hence, the velocity state of a rigid body does not keep a unique representation. The correctness of the velocity state representation of a rigid body would be validated by resorting to equivalence relationship concepts, as well as basic concepts, in this case, of kinematics. For example, let us consider that Q is a point of body m where Q P . Then, the velocity state of body m using point Q as a reference to the dual part may be written as follows:
  j V Q m =   j ω m v Q
Expressions (19) and (20) are two representations of the velocity state of body m. The primal parts are the same, which indicates that the angular velocity vector is a property of a rigid body. On the other hand, the dual parts may be related through the properties of helicoidal vector fields [18]. That is to say,
v P = D (   j V Q m ) + P (   j V Q m ) × r P / Q = v Q +   j ω m × r P / Q
where r P / Q is the vector beginning at point Q an ending at point P.
Let
V 1 = ω 1 v 1 , V 2 = ω 2 v 2 , V 3 = ω 3 v 3
be three elements of the Lie algebra s e ( 3 ) of the Euclidean group S E ( 3 ) . For the sake of completeness, the fundamental operations and their corresponding properties of this algebra are enumerated below:
  • Addition (+)
    V 1 + V 2 = ω 1 v 1 + ω 2 v 2 = ω 1 + ω 2 v 1 + v 2
    The operation of addition (+) has the following properties:
    (a)
    Commutative V 1 + V 2 = V 2 + V 1 ;
    (b)
    Associative V 1 + ( V 2 + V 3 ) = ( V 1 + V 2 ) + V 3 ;
    (c)
    Identical additive 0 , V 1 + 0 = V 1 ;
    (d)
    scalar multiplication, λ V 1 = λ ω 1 λ v 1 . λ R .
  • Lie bracket ( )
    V 1 V 2 = ω 1 × ω 2 ω 1 × v 2 ω 2 × v 1
    The Lie bracket has the following operations
    (a)
    Anti-symmetry (anticommutativity), V 1 V 2 = V 2 V 1
    (b)
    Alternativity (nilpotent), V 1 V 1 = 0 ;
    (c)
    Bilinearity, λ V 1 + μ V 2 V 3 = λ V 1 V 3 + μ V 2 V 3 ;
    (d)
    Jacobi identity, [ V 1 V 2 ] V 3 + [ V 3 V 1 ] V 2 + [ V 2 V 3 ] V 1 = 0 ;
    (e)
    product rule, V 1 f V 2 = V 1 ( f ) V 2 + f V 1 V 2 .
  • Bilinearity
    (a)
    Klein form ( { ; } ) , { V 1 ; V 2 } = ω 1 · v 2 + ω 2 · v 1 ;
    (b)
    Killing form ( ; ) , ( V 1 ; V 2 ) = ω 1 · ω 2 .
The Klein form is an operation that has proven to be very convenient in the kinematic analysis of parallel platforms of complex architecture, since it is directly related to the concept of reciprocity of two elements of the screw algebra [19]. Moreover, it is a powerful tool in the kinetic analysis of these robots, when elegantly combining screw theory with the principle of virtual power.
On the other hand, the Killing form has been shown to be a tool for, for example, the expansion of the Aronhold–Kennedy theorem, as detailed in [20].

2.4. Open Kinematic Chains

Topologically, a robot manipulator is a set of rigid bodies called links assembled in such a way that adjacent links are physically connected by means of geometric forms. The Lie algebras are a trusted mathematical resource to approach diverse issues of robot manipulators, for example, kinematic analysis. In this subsection, the velocity, acceleration, and jerk states of rigid body in screw form are developed for robot manipulators. For details relating this subsection, the reader is referred to [19,21].

2.4.1. Velocity

Figure 4 shows two bodies labeled m and j connected by means of a helical pair   j $ m . Assuming that Q is a point on the Mozzi axis, then its velocity vector v Q may be computed as
v Q =   j ω m h   j s ^ m
where h is the pitch of the helical pair or screw and   j s ^ m is a unitary vector in the same direction such that   j ω m . Note that according to Chasles’ theorem, all points on the Mozzi axis have the same velocity vector. Furthermore, the velocity vector of point P of body m, which is located out of the Mozzi axis, can be computed by resorting to elementary kinematics,
v P = v Q +   j ω m × r P / Q
where point P would be considered as the reference pole in turn of body m. Substituting Equation (23) into Equation (24) and taking into account that   j ω m =   j ω m   j s ^ m , it follows that
v P =   j ω m h   j s ^ m +   j ω m   j s ^ m × r P / Q =   j ω m ( h   j s ^ m +   j s ^ m × r P / Q )
Hence,
v P =   j ω m   j s P m
where   j s P m = h   j s ^ m +   j s ^ m × r P / Q is the moment vector produced by   j s ^ m about point P. Thereafter, based on Equation (19) we have
  j V P m =   j ω m   j v P m =   j ω m   j s ^ m   j ω m   j s P m =   j ω m   j s ^ m   j s P m
or
  j V P m =   j ω m   j $ m
where   j $ m =   j s ^ m   j s P m .
Expression (28) proves that the velocity state of a rigid body has a direct connection with the concept of an infinitesimal screw. Certainly, a twist can be represented as a normalized screw axis s ^ denoting the direction of the motion, multiplied by a joint rate ω along the screw axis. The extension of this concept to robot manipulators is the next step.
Theorem 3.
Figure 5 shows an open kinematic chain where the bodies labeled m and j are uninary links. For example, j and m would be considered as the fixed link and the end effector of a serial manipulator, respectively. Between bodies m and j, there are intermediate links forming the spinal column of the serial manipulator. The velocity state of body m as measured from body j considering point O as the reference pole is given by
  j V O m =   j V O j + 1 +   j + 1 V O j + 2 + +   m 2 V O m 1 +   m 1 V O m
Proof. 
Let us consider three bodies, j, k, and m, that are serially connected by helical pairs, as shown in Figure 6. Furthermore, let us consider that β is an arbitrary vector embedded in body m. Then, the time derivative of β may be obtained considering the reference frames j and k as follows:
  j d β d t =   j ω m × β ,   k d β d t =   k ω m × β
On the other hand, to change the reference frame of the time derivative, consider that
  j d β d t =   j ω k × β +   k d β d t
Thus,
  j d β d t =   j ω k × β +   k ω m × β = (   j ω k +   k ω m ) × β
which leads to
  j ω m × β = (   j ω k +   k ω m ) × β
Therefore,
  j ω m =   j ω k +   k ω m
and assuming that j = m 2 and k = m 1 , one obtains
  m 2 ω m =   m 2 ω m 1 +   m 1 ω m
By the induction hypothesis supposes that
  j + 1 ω m =   j + 1 ω j + 2 + +   m 2 ω m 1 +   m 1 ω m
and substituting k = j + 1 into Equation (34), we have
  j ω m =   j ω j + 1 +   j + 1 ω m
Hence, by resorting to Equation (36), one obtains
  j ω m =   j ω j + 1 +   j + 1 ω j + 2 + +   m 2 ω m 1 +   m 1 ω m
To obtain   j v O m , consider Figure 7. In that concern, without loss of generality the reference pole O would be precisely the origin of the reference frame O _ X Y Z .
As shown in the second subfigure of Figure 7, the velocity vector   j v O k of point O, fixed in body k as observed from body j, is given by
  j v O k =   j h k   j ω k +   j ω k × r O / Q
Similarly, according to the last subfigure of Figure 7, it follows that the velocity vector   k v O m of point O, fixed in body m as observed from body k, is given by
  k v O m =   k h m   k ω m +   k ω m × r O / P
Hence, the velocity vector   j v O m of point O, fixed in body body m as observed from body j, results in
  j v O m =   j h k   j ω k +   j ω k × r O / Q +   k h m   k ω m +   k ω m × r O / P
or
  j v O m =   j v O k +   k v O m
Following the trend of the angular velocity vector   j ω m (see Equation (38), for the open kinematic chain of Figure 5), one obtains
  j v O m =   j v O j + 1 +   j + 1 v O j + 2 + +   m 2 v O m 1 +   m 1 v O m
Furthermore, taking into account that the velocity state   j V O m is given by
  j V O m =   j ω m   j v O m
from Equations (38) and (43), we have that
  j V O m =   j ω m   j v O m =   j ω j + 1 +   j + 1 ω j + 2 + +   m 2 ω m 1 +   m 1 ω m   j v O j + 1 +   j + 1 v O j + 2 + +   m 2 v O m 1 +   m 1 v O m
or
  j V O m =   j ω j + 1   j v O j + 1 +   j + 1 ω j + 2   j + 1 v O j + 2 + +   m 2 ω m 1   m 2 v O m 1 +   m 1 ω m   m 1 v O m
which finally leads to Equation (29) completing the proof of Theorem 3.
It is interesting to note that Theorem 3 allows us to formulate the velocity state   j V O m in terms of the infinitesimal screws of the serial kinematic chain by resorting recursively to Equation (28) as follows:
  j V O m =   j ω j + 1   j $ j + 1 +   j + 1 ω j + 2   j + 1 $ j + 2 + +   m 2 ω m 1   m 2 $ m 1 +   m 1 ω m   m 1 $ m
Expression (47) is one of the highlights of the theory of screws. Its introduction in the 1980s [22] marked the successful beginning of screw theory in the kinematic analysis of robotic manipulators.

2.4.2. Acceleration

The acceleration state A of rigid body refers to the set of kinematic elements required to characterize the acceleration of a rigid body. Let m be a rigid body in general motion as observed from body j. At a glance the acceleration state of rigid body m would be considered as the time derivative of the corresponding velocity state. However, this is an erroneous assumption because the helicoidal vectorial field conditions are not satisfied with this hypothesis [18]. Brand [23] defined the acceleration state of body m as measured from body j, the six-dimensional vector   j A P m , as follows:
  j A P m =   j α m a P   j ω m × v P
where   j α m =   j d d t   j ω m is the angular acceleration vector, while a P is the acceleration vector of point P.
  j α m =   j d d t   j ω m =   j d d t   j ω m   j s ^ m =   j α m   j s ^ m
Hence
  j α m =   j α m   j s ^ m
On the other hand, with reference to Figure 4 and Equation (23), the acceleration of point Q on the Mozzi axis is given by
a Q =   j d d t v Q =   j d d t (   j ω m h   j s ^ m )
Thus,
a Q =   j α m h   j s ^ m
Furthermore, by resorting to elementary kinematics, the acceleration of point P can be calculated as
a P = a Q +   j α m × r P / Q +   j ω m × (   j ω m × r P / Q )
However, v P =   j ω m × r P / Q . Therefore,
a P   j ω m × v P = a Q +   j α m × r P / Q
The substitution of Equations (50) and (52) into Equation (55) leads to
a P   j ω m × v P =   j α m h   j s ^ m +   j α m   j s ^ m × r P / Q
or
a P   j ω m × v P =   j α m ( h   j s ^ m +   j s ^ m × r P / Q )
Therefore,
a P   j ω m × v P =   j α m   j s P m
Finally, the substitution of Equations (50) and (57) into Equation (48) yields [21]
  j A P m =   j α m   j $ m
Expression (58) is a confirmation of the correctness of the acceleration state of the rigid body introduced by Rico and Duffy [21], which, by the way, guarantees that the acceleration state, like the velocity state, is a twist about a screw.
Theorem 4.
Figure 5 shows an open kinematic chain where the bodies labeled m and j are unitary links. For example, j and m would be considered as the fixed link and the end effector of a serial manipulator, respectively. Between bodies m and j, there are intermediate links forming the spinal column of the serial manipulator. The acceleration state of body m as measured from body j is given by
  j A O m =   j A O j + 1 +   j + 1 A O j + 2 + +   m 2 A O m 1 +   m 1 A O m +   j A O m +   j A O m
where   j A O m is the Coriolis-Lie vector of acceleration, which is given by
  j A O m =   j V O j + 1   j + 1 V O j + 2 +   j + 2 V O j + 3 + +   m 2 V O m 1 +   m 1 V O m +   j + 1 V O j + 2   j + 2 V O j + 3 +   j + 3 V O j + 4 + +   m 2 V O m 1 +   m 1 V O m + +   m 2 V O m 1   m 1 V O m
Proof. 
Let us consider three bodies, j, k, and m, in relative motion; see Figure 6. The angular acceleration vector   j α m may be obtained upon Equation (34) as follows:
  j α m =   j d d t   j ω m =   j d d t (   j ω k +   k ω m ) =   j d d t   j ω k +   j ω k ×   k ω m +   k d d t   k ω m
Hence,
  j α m =   j α k +   k α m +   j ω k ×   k ω m
Considering that j = m 2 and k = m 1 , then, one obtains that
  m 2 α m =   m 2 α m 1 +   m 1 α m +   m 2 ω m 1 ×   m 1 ω m
Afterwards, based on the induction hypothesis, assume that
  j + 1 α m =   j + 1 α j + 2 + +   m 2 α m 1 +   m 1 α m +   j + 1 ω j + 2 × (   j + 2 ω j + 3 + +   m 1 ω m ) +   j + 2 ω j + 3 × (   j + 3 ω j + 4 + +   m 1 ω m ) + +   m 2 ω m 1 ×   m 1 ω m
Substitution k = j + 1 into Equation (62) yields
  j α m =   j α j + 1 +   j + 1 α m +   j ω j + 1 ×   j + 1 ω m
With the substitution of Equation (64) into Equation (65), we have the result
  j α m =   j α j + 1 + +   m 2 α m 1 +   m 1 α m +   j ω j + 1 × (   j + 1 ω j + 2 + +   m 1 ω m ) +   j + 1 ω j + 2 × (   j + 2 ω j + 3 + +   m 1 ω m ) + +   m 2 ω m 1 ×   m 1 ω m
On the other hand, referring to Figure 6 and Figure 7, consider that the acceleration vector of point O is given by
  j a O m =   j a O k +   k a O m + 2   j ω k ×   k v O m
Thereafter, based on Equation (48) it follows that
  j A O m =   j α m   j a O m   j ω m ×   j v O m =   j α k +   k α m +   j ω k ×   k ω m   j a O k +   k a O m + 2   j ω k ×   k v O m (   j ω k +   k ω m ) × (   j v O k +   k v O m ) =   j α k   j a O k   j ω k ×   j v O k +   k α m   k a O m   k ω m ×   k v O m +   j ω k ×   k ω m   j ω k ×   k v O m   k ω m ×   j v O k
or
  j A O m =   j A O k +   k A O m +   j V O k   k V O m
Finally, by the induction hypothesis, one obtains Equation (59), completing the proof of Theorem 4. It is interesting to note that by resorting to Equation (58), the acceleration state   j A O m may be written in screw form as follows:
  j A O m =   j α j + 1   j $ j + 1 +   j + 1 α j + 2   j + 1 $ j + 2 + +   m 2 α m 1   m 2 $ m 1 +   m 1 α m   m 1 $ m +   j $ a m
where   j $ a m is given by
  j $ a m =   j ω j + 1   j $ j + 1   j + 1 ω j + 2   j + 1 $ j + 2 + +   m 1 ω m   m 1 $ m +   j + 1 ω j + 2   j + 1 $ j + 2   j + 2 ω j + 3   j + 2 $ j + 3 + +   m 1 ω m   m 1 $ m + +   m 2 ω m 1   m 2 $ m 1   m 1 ω m   m 1 $ m

3. Symbolic Computation: Maple Sheets

In this section, Maple sheets for computing some relevant operations of the Lie algebra s e ( 3 ) of the Euclidean group S E ( 3 ) , as well as the velocity and acceleration states of rigid bodies in open kinematic chains, are reported.
Lie product between two six-dimensional vectors:
Mathematics 12 02538 i001
Klein form between two six-dimensional vectors:
Mathematics 12 02538 i002
Killing form between two six-dimensional vectors:
Mathematics 12 02538 i003
Velocity and acceleration states of rigid body:
Mathematics 12 02538 i004

4. Examples

4.1. Planar RPR Serial Manipulator

Figure 8 shows an open kinematic chain playing the role of a three-degree-of-freedom planar serial manipulator.
The base link 0 is connected to the end effector m by means of a prismatic joint P assembled between two revolute joints R. Two rotary actuators q 1 and q 3 associated with the revolute joints endowed with a linear actuator q 2 associated with the lonely prismatic joint are employed to control the position and orientation of the end effector m as observed from the base link 0. More specifically, q 1 denotes the orientation of the skirt piston with respect to the base link, q 2 denotes the variable length of the connecting rod including the corresponding offsets, and q 3 denotes the orientation of the end effector with respect to the piston. The example consists of symbolically computing the velocity and acceleration states of the end effector m as observed from the base link 0, assuming that Q is the reference pole.
Referring to Figure 8, the infinitesimal screws of the RPR serial manipulator are given by
  1 $ 2 =   1 s ^ 2 r O / Q ×   1 s ^ 2 ,   2 $ 3 = 0   2 s ^ 3 ,   3 $ 4 =   3 s ^ 4 r P / Q ×   3 s ^ 4
where the primal parts of the screws representing the revolute joints are given by
  1 s ^ 2 =   3 s ^ 4 = 0 0 1 ,   2 s ^ 3 = cos q 1 sin q 1 0
The application of the algorithm yields the following velocity and acceleration states of body 4 as measured from body 1, taking point Q as the reference pole.
  1 V Q 4 = 0 0   1 ω 2 +   3 ω 4 [ q 2 sin q 1 + a sin ( q 1 + q 3 ) ]   1 ω 2 + q 2 cos q 1   2 ω 3 a sin ( q 1 + q 3 )   3 ω 4 [ q 2 cos q 1 + a cos ( q 1 + q 3 ) ]   1 ω 2 + q 2 sin q 1   2 ω 3 a cos ( q 1 + q 3 )   3 ω 4 0   1 A Q 4 = 0 0   1 α 2 +   3 α 4 ( q 2 sin q 1 a sin q 1 cos q 3 a cos q 1 sin q 3 )   1 α 2 + q 2 cos q 1   2 α 3 a ( sin q 1 cos q 3 + cos q 1 sin q 3 )   3 α 4 q 2 sin q 1   1 ω 2   2 ω 3 + q 2 cos q 1   1 ω 2   3 ω 4 + q 2 sin q 1   2 ω 3   3 ω 4 ( q 2 cos q 1 + a cos q 1 cos q 3 a sin q 1 sin q 3 )   1 α 2 + q 2 sin q 1   2 α 3 + a ( cos q 1 cos q 3 sin q 1 sin q 3 )   3 α 4 + q 2 cos q 1   1 ω 2   2 ω 3 + q 2 sin q 1   1 ω 2   3 ω 4 q 2 cos q 1   2 ω 3   3 ω 4 0
where   1 ω 2 = q ˙ 1 ,   2 ω 3 = q ˙ 2 ,   3 ω 4 = q ˙ 3 ,   1 α 2 = q ¨ 1 ,   2 α 3 = q ¨ 2 , and   3 α 4 = q ¨ 3 .

4.2. PUMA Robot

PUMA stands for Programmable Universal Assembly Machine or Programmable Universal Manipulator Arm. This non-redundant spatial serial manipulator was developed by General Motors based on the research work of Victor Scheinman, which was supported by MIT and Stanford University. Later on, UNIMATION provided the final impulse for the inclusion of the PUMA robot in the industry. Without a doubt, the PUMA is an emblematic robot with which modern robotics manifested its intention to replicate the movement of the human body. By that time, the potential applications of the PUMA were immediately recognized in the industry, raising great expectations, and it quickly went from being a subject academic research to a practical application. In this subsection, the velocity and acceleration analyses of the PUMA robot are approached by means of the theory of screws, which is isomorphic to the Lie algebra s e ( 3 ) of the Euclidean group S E ( 3 ) and the motor algebra of von Mises [17]. For the sake of completeness, the position analysis of the robot is included.

4.2.1. Position Analysis

The PUMA robot is a serial manipulator with six spatial degrees of freedom composed of links articulated by means of revolute joints; see Figure 9.
The generalized coordinates of the PUMA robot are denoted by q i ( i = 1 , 2 , 3 , , 6 ) , while the parameters, known as the Denavit–Hartenberg parameters, are notated as a and d . The handling of several reference frames attached to different links of the robot is a natural way to solve the position analysis. With this in mind, the homogeneous transformation matrices using the Denavit–Hartenberg convention are formulated as follows
  0 T 1 = cos q 1 0 sin q 1 0 sin q 1 0 cos q 1 0 0 1 0 0 0 0 0 1 ,   1 T 2 = cos q 2 sin q 2 0 a 2 cos q 2 sin q 2 cos q 2 0 a 2 sin q 2 0 0 1 d 2 0 0 0 1 ,   2 T 3 = cos q 3 0 sin q 3 a 3 cos q 3 sin q 3 0 cos q 3 a 3 sin q 3 0 1 0 0 0 0 0 1 ,   3 T 4 = cos q 4 0 sin q 4 0 sin q 4 0 cos q 4 0 0 1 0 d 4 0 0 0 1 ,   4 T 5 = cos q 5 0 sin q 5 0 sin q 5 0 cos q 5 0 0 1 0 0 0 0 0 1 ,   5 T 6 = cos q 6 sin q 6 0 0 sin q 6 cos q 6 0 0 0 0 1 d 6 0 0 0 1
Thereafter, the homogeneous matrix of the end-effector m as measured from the base link 0, the matrix   0 T 6 , is computed as follows:
  0 T 6 =   0 T 1   1 T 2   2 T 3   3 T 4   4 T 5   5 T 6
Next, let us consider an arbitrary point Q of the end effector m. The coordinates ( Q X , Q Y . Q Z ) of point Q expressed in the fixed reference frame O _ X Y Z are obtained with the coordinates ( Q X 6 , Q Y 6 , Q Z 6 ) of the same point Q but expressed in the reference frame O 6 X 6 Y 6 Z 6 as follows:
X Q Y Q Z Q 1 =   0 T 6 x Q y Q z Q 1

4.2.2. Velocity Analysis

To approach the velocity analysis, Figure 10 shows the infinitesimal screws associated with the revolute joints of the PUMA robot.
Let ω = ω X i ^ + ω Y j ^ + ω Z k ^ be the angular velocity vector of the end effector, link m = 6 , as observed from the base link 0. Moreover, let v O 6 = v X i ^ + v Y j ^ + v Z k ^ be the velocity vector of point O 6 of body 6 expressed in the global reference frame O _ X Y Z , where the point O 6 is the origin of the reference frame attached to the end effector m. The velocity state   0 V O 6 6 of the end effector as measured from the base link may be written in matrix form as follows:
  0 V O 6 6 = ω X ω Z ω Z v X v Y v Z
where point O 6 is the reference pole. Moreover, the velocity state   0 V O 6 6 may be written in screw form as follows:
q ˙ 1   0 $ 1 + q ˙ 2   1 $ 2 + q ˙ 3   2 $ 3 + q ˙ 4   3 $ 4 + q ˙ 5   4 $ 5 + q ˙ 6   5 $ 6 =   0 V O 6 6
where the infinitesimal screws are given by
  0 $ 1 = k ^ k ^ × r O 6 / O ,   1 $ 2 = r ^ B / O r ^ B / O × r O 6 / O ,   2 $ 3 = r ^ B / O r ^ B / O × r O 6 / O 2 ,   3 $ 4 = r ^ O 4 / O 3 r ^ O 4 / O 3 × r O 6 / O 4 ,   4 $ 5 = ( r O 6 / O 4 × r O 4 / O 3 ) ^ ( r O 6 / O 4 × r O 4 / O 3 ) ^ × r O 6 / O 4 ,   5 $ 6 = r ^ O 6 / O 4 0
and r ^ O m / O n is a unitary vector parallel to r O m / O n .
Expression (76) may be written in a matrix vector form as follows
J Q v = V O 6
where J is the Jacobian matrix of the manipulator, which is given by
J =   0 $ 1   1 $ 2   2 $ 3   3 $ 4   4 $ 5   5 $ 6
On the other hand, Q v = q ˙ 1 q 2 ˙ q ˙ 3 q ˙ 4 q ˙ 5 q ˙ 6 T is the first-order active matrix of the PUMA robot.
Expression (78) is the input–output equation of velocity of the PUMA robot and is valid for both the inverse and the direct velocity analyses of the robot manipulator.

4.2.3. Acceleration Analysis

Let α = α X i ^ + α Y j ^ + α Z k ^ be the angular acceleration vector of the end effector m = 6 as observed from the fixed link 0. Moreover, let a O 6 = a X i ^ + a Y j ^ + a Z k ^ be the acceleration vector of point O 6 . The acceleration state of the end effector, as observed from the base link, taking point O 6 as the reference pole as given by
A O 6 = α a O 6 ω × v O 6
Moreover, the acceleration equation in screw form of the PUMA robot may be written as follows:
q ¨ 1   0 $ 1 + q ¨ 2   1 $ 2 + q ¨ 3   2 $ 3 + q ¨ 4   3 $ 4 + q ¨ 5   4 $ 5 + q ¨ 6   5 $ 6 + $ a =   0 A O 6 6
where
$ a = q ˙ 1   0 $ 1 q ˙ 2   1 $ 2 + q ˙ 3   2 $ 3 + q ˙ 4   3 $ 4 + q ˙ 5   4 $ 5 + q ˙ 6   5 $ 6 + q ˙ 2   1 $ 2 q ˙ 3   2 $ 3 + q ˙ 4   3 $ 4 + q ˙ 5   4 $ 5 + q ˙ 6   5 $ 6 + q ˙ 3   2 $ 3 q ˙ 4   3 $ 4 + q ˙ 5   4 $ 5 + q ˙ 6   5 $ 6 + q ˙ 4   3 $ 4 q ˙ 5   4 $ 5 + q ˙ 6   5 $ 6 + q ˙ 5   4 $ 5 q ˙ 6   5 $ 6
Expression (80) may be written in matrix vector form by introducing the Jacobian matrix J of the velocity analysis as follows:
J Q a + $ a = A O 6
where Q a = q ¨ 1 q ¨ 2 q ¨ 3 q ¨ 4 q ¨ 5 q ¨ 6 T is the active acceleration matrix of the manipulator, which is formed by the generalized accelerations of the robot.

4.2.4. Numerical Application

For this numerical case, let us consider the Denavit–Hartenberg parameters for a standard PUMA manipulator (Figure 9) shown in the Table 1. On the other hand, let us consider that in the reference configuration, all the generalized coordinates of the robot are given by 10 . Furthermore, from its reference configuration, the generalized coordinates of the PUMA manipulator follow periodic functions given by q 1 = 0.2 sin ( t ) [rad], q 2 = 0.15 sin ( t ) [rad], q 3 = 0.25 sin ( t ) [rad], q 4 = 0.25 sin ( t ) [rad], q 5 = 0.15 sin ( t ) [rad] and q 6 = 0.3 sin ( t ) [rad], where 0 t 2 π [s]. With these data, we must compute the time history of the direct kinematics of the PUMA robot. The relevant resulting plots of the numerical example are shown in Figure 11.

5. Conclusions

This work presents an application of the Lie algebra s e ( 3 ) of the Euclidean group S E ( 3 ) , which is isomorphic to the theory of screws and the motor algebra, to the infinitesimal kinematics of serial manipulators. The Lie algebra, as is shown in this contribution, has been consolidated as a highly efficient mathematical tool in the analysis and synthesis of robotic manipulators. The algorithm presented allows us to symbolically compute the velocity and acceleration characteristics of the end effector of serial manipulators. In the case of complex serial manipulators, numerical computations are preferred, owing to the generation of expressions with long terms. To exemplify the presented methodology, two serial manipulators are analysed: (1) a three-dof planar serial manipulator, (2) the six-dof PUMA robot.
These procedures can be extended without significant effort to higher-order analyses such as the jerk and jounce. The fact that the screw algebra, which is isomorphic to the Lie algebra s e ( 3 ) , has been reported and successfully tested in the analysis of robotic manipulators in previous works proves its undoubted usefulness and correctness. Unlike other algebras, the screw algebra has demonstrated its applicability and versatility in higher-order analyses, a real obstacle for other algebras such as the quaternion algebra, among others. Now, the task is to search for areas of opportunity dealing with practical applications. These applications may also require the use of hybrid methods for complex kinematic architectures, for example, in [24], where the positional analysis of a complex parallel platform is investigated by combining the well-known Newton–Raphson algorithm with methods based on homotopic continuation. In that regard, the Maple sheets provided in this work are sufficiently readable so that a non-expert user can apply them without substantial problems in specific robot manipulators. However, especially in problems where it is not possible to generate closed-form solutions, it becomes necessariy to resort to numerical methods.

Author Contributions

Methodology, formal analysis, writing—original draft preparation, J.G.-A. and M.A.G.-M.; software, validation, writing—review and editing, J.M.T.-M. and X.Y.S.-C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ball, R.S. A Treatise on the Theory of Screws; Cambridge University Press: Cambridge, UK, 1900; reprinted 1998. [Google Scholar]
  2. Hunt, K.H. Kinematic Geometry of Mechanisms; Oxford University Press: New York, NY, USA, 1978. [Google Scholar]
  3. Phillips, J. Freedom in Machinery: Introducing Screw Theory; Cambridge University Press: Cambridge, UK, 1984; Volume 1. [Google Scholar]
  4. Karger, A. Singularity analysis of serial robot-manipulators. ASME Mech. Des. 1996, 118, 520–525. [Google Scholar] [CrossRef]
  5. Müller, A. Higher derivatives of the kinematic mapping and some applications. Mech. Mach. Theory 2014, 76, 70–85. [Google Scholar] [CrossRef]
  6. Wu, Y.; Carricato, M. Identification and geometric characterization of Lie triple screw systems and their exponential images. Mech. Mach. Theory 2017, 107, 305–323. [Google Scholar] [CrossRef]
  7. Zhang, L.; Zhao, Y.; Zhao, T. Fundamental equation of mechanism kinematic geometry: Mapping curve in se(3) to counterpart in SE(3). Mech. Mach. Theory 2020, 146, 103732. [Google Scholar] [CrossRef]
  8. Sanchez-Garcia, A.J.; Rico, J.M.; Cervantes-Sanchez, J.J.; Lopez-Custodio, P.C. A mobility determination method for parallel platforms based on the Lie algebra of SE(3) and Its Subspaces. ASME J. Mech. Robot. 2021, 13, 031015. [Google Scholar] [CrossRef]
  9. Chen, J.; Huang, Z.; Tian, Q. A multisymplectic Lie algebra variational integrator for flexible multibody dynamics on the special Euclidean group SE(3). Mech. Mach. Theory 2022, 174, 104918. [Google Scholar] [CrossRef]
  10. Hassani, H.; Avazzadeh, Z.; Agarwal, P.; Ebadi, M.J.; Eshkaftaki, A.B. Generalized Bernoulli—Laguerre Polynomials: Applications in Coupled Nonlinear System of Variable-Order Fractional PDEs. J. Optim. Theory Appl. 2024, 200, 371–393. [Google Scholar] [CrossRef]
  11. Avazzadeh, Z.; Hassani, H.; Ebadi, M.J.; Bayati Eshkaftaki, A.; Hendy, A.S. An optimization method for solving a general class of the inverse system of nonlinear fractional order PDEs. Int. J. Comput. Math. 2024, 101, 138–153. [Google Scholar] [CrossRef]
  12. Hodge, W.V.D.; Pedoe, D. Methods of Algebraic Geometry; Cambridge University Press: Cambridge, UK, 1994; Volume 1. [Google Scholar]
  13. Mason, M.T. Mechanics of Robotic Manipulation; MIT Press: Cambridge, MA, USA, 2001. [Google Scholar]
  14. Palais, B.; Palais, R.; Rodi, S. A Disorienting Look at Euler’s Theorem on the Axis of a Rotation. Am. Math. Mon. 2009, 116, 892–909. [Google Scholar] [CrossRef]
  15. Chasles, M. Note sur les propriétés générales du système de deux corps semblables entr’eux. Bull. Sci. Math. Astron. Phys. Chim. 1830, 4, 321–326. [Google Scholar]
  16. Mozzi, G. Discorso Matematico Sopra il Rotamento Momentaneo Dei Corpi; Nella Stamperia di Donato Campo: Campo, Italy, 1763. [Google Scholar]
  17. von Mises, R. Motorrechnung; ein neues Hilfsmittel der Mechanik. ZAMM 1924, 4, 155–181. [Google Scholar] [CrossRef]
  18. Gallardo, J.; Rico, J.M. Screw theory and helicoidal fields. In Proceedings of the ASME 25th Biennial Mechanisms Conference, Atlanta, GA, USA, 13–16 September 1998. Paper DETC98/MECH-5893. [Google Scholar]
  19. Gallardo-Alvarado, J. Análisis Cinemáticos de Orden Superior de Cadenas Espaciales, Mediante el Álgebra de Tornillos, y sus Aplicaciones. Doctoral Dissertation, Instituto Tecnológico de la Laguna, Torreón, México, April 1999. (In Spanish). [Google Scholar]
  20. Valderrama-Rodríguez, J.I.; Rico, J.M.; Cervantes-Sánchez, J.J. A general method for the determination of the instantaneous screw axes of one-degree-of-freedom spatial mechanisms. Mech. Sci. 2020, 11, 91–99. [Google Scholar] [CrossRef]
  21. Rico-Martinez, J.M.; Duffy, J. An application of screw algebra to the acceleration analysis of serial chains. Mech. Mach. Theory 1996, 31, 445–457. [Google Scholar] [CrossRef]
  22. Sugimoto, K.; Duffy, J. Application of linear algebra to screw systems. Mech. Mach. Theory 1982, 17, 73–83. [Google Scholar] [CrossRef]
  23. Brand, L. Vector and Tensor Analysis; John Wiley & Sons: Hoboken, NJ, USA, 1947. [Google Scholar]
  24. Gallardo-Alvarado, J.; Garcia-Murillo, M.A.; Aguilera-Camacho, L.D.; Alcaraz-Caracheo, L.A.; Sandoval-Castro, X.Y. An approach to solving the forward kinematics of the 5-RPUR (3T2R) parallel manipulator. J. Mech. Sci. Technol. 2023, 37, 1443–1453. [Google Scholar] [CrossRef]
Figure 1. A line l in the physical space passing through point p.
Figure 1. A line l in the physical space passing through point p.
Mathematics 12 02538 g001
Figure 2. A line l in physical space.
Figure 2. A line l in physical space.
Mathematics 12 02538 g002
Figure 3. A rigid body m in general motion.
Figure 3. A rigid body m in general motion.
Mathematics 12 02538 g003
Figure 4. Two bodies m and j connected by a helical pair   j $ m .
Figure 4. Two bodies m and j connected by a helical pair   j $ m .
Mathematics 12 02538 g004
Figure 5. Open kinematic chain whose adjacent links are connected by means of helical pairs.
Figure 5. Open kinematic chain whose adjacent links are connected by means of helical pairs.
Mathematics 12 02538 g005
Figure 6. Three bodies j, k, and m in relative motion.
Figure 6. Three bodies j, k, and m in relative motion.
Mathematics 12 02538 g006
Figure 7. Equivalent motion between bodies j, k, and m.
Figure 7. Equivalent motion between bodies j, k, and m.
Mathematics 12 02538 g007
Figure 8. RPR planar serial manipulator.
Figure 8. RPR planar serial manipulator.
Mathematics 12 02538 g008
Figure 9. Reference frames and parameters of the PUMA manipulator.
Figure 9. Reference frames and parameters of the PUMA manipulator.
Mathematics 12 02538 g009
Figure 10. Screws of the PUMA robot.
Figure 10. Screws of the PUMA robot.
Mathematics 12 02538 g010
Figure 11. Time history of the infinitesimal kinematics of the end effector of the PUMA robot.
Figure 11. Time history of the infinitesimal kinematics of the end effector of the PUMA robot.
Mathematics 12 02538 g011
Table 1. Denavit–Hartenberg parameters of the PUMA robot.
Table 1. Denavit–Hartenberg parameters of the PUMA robot.
Link Offset [mm]Joint AngleLink Length [mm]Twist Angle [rad]
d 1 = 0 θ 1 = q 1 a 1 = 0 α 1 = π 2
d 2 = 144.09 θ 2 = q 2 a 2 = 431.8 α 2 = 0
d 3 = 0 θ 3 = q 3 a 3 = 20.32 α 3 = π 2
d 4 = 433.07 θ 4 = q 4 a 4 = 0 α 4 = π 2
d 5 = 0 θ 5 = q 5 a 5 = 0 α 4 = π 2
d 6 = 56.25 θ 6 = q 6 a 6 = 0 α 6 = 0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gallardo-Alvarado, J.; Garcia-Murillo, M.A.; Tabares-Martinez, J.M.; Sandoval-Castro, X.Y. Symbolic Computation of the Lie Algebra se(3) of the Euclidean Group SE(3): An Application to the Infinitesimal Kinematics of Robot Manipulators. Mathematics 2024, 12, 2538. https://doi.org/10.3390/math12162538

AMA Style

Gallardo-Alvarado J, Garcia-Murillo MA, Tabares-Martinez JM, Sandoval-Castro XY. Symbolic Computation of the Lie Algebra se(3) of the Euclidean Group SE(3): An Application to the Infinitesimal Kinematics of Robot Manipulators. Mathematics. 2024; 12(16):2538. https://doi.org/10.3390/math12162538

Chicago/Turabian Style

Gallardo-Alvarado, Jaime, Mario A. Garcia-Murillo, Juan Manuel Tabares-Martinez, and X. Yamile Sandoval-Castro. 2024. "Symbolic Computation of the Lie Algebra se(3) of the Euclidean Group SE(3): An Application to the Infinitesimal Kinematics of Robot Manipulators" Mathematics 12, no. 16: 2538. https://doi.org/10.3390/math12162538

APA Style

Gallardo-Alvarado, J., Garcia-Murillo, M. A., Tabares-Martinez, J. M., & Sandoval-Castro, X. Y. (2024). Symbolic Computation of the Lie Algebra se(3) of the Euclidean Group SE(3): An Application to the Infinitesimal Kinematics of Robot Manipulators. Mathematics, 12(16), 2538. https://doi.org/10.3390/math12162538

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