In this chapter, using the ring properties of dual number algebra, vector and tensor calculus, a computing method for the higher-order acceleration vector field properties in general rigid body motion is proposed. The higher-order acceleration field of a rigid body in a general motion is uniquely determined by higher-order time derivative of a dual twist. For the relative kinematics of rigid body motion, equations that allow the determination of the higher-order acceleration vector field are given, using an exponential Brockett-like formula in the dual Lie algebra. In particular cases, the properties for velocity, acceleration, jerk, and jounce fields are given. This approach uses the isomorphism between the Lie algebra of the rigid displacements se(3), of the Special Euclidean group, SE3,and the Lie algebra of dual vectors. The results are coordinate free and in a closed form.
- higher-order kinematics
- dual algebra
- lie group
The kinematic analysis of multibody systems has been traditionally considered as the determination of the positions, velocities, accelerations, jerks and jounces of their constitutive members. This is an old field with a long history, which has attracted the attention of mathematicians and engineers. Michel Chasles discovered (1834) that any rigid body displacement is equivalent to a screw displacement . Screw theory is an efficient mathematical tool for the study of spatial kinematics. The pioneering work of Ball , the treatises of Hunt , and Phillips  and the multitude of contributions appearing in the literature are evidence of this. The isomorphism between screw theory and the Lie algebra, se(3), of the Special Euclidean group, , provide with a wealth of results and techniques from modern differential geometry and Lie group theory [5, 6, 7, 8, 9].
A kinematic mapping relates the motion of a rigid body to the joint motions of a kinematic chain. Its time derivatives yield the twist, acceleration, jerk and jounce etc. of the body. Time derivatives of the twists of members in a kinematic chain and derivatives of screws are essential operations in kinematics. Recognizing the Lie group nature of rigid body motions, and correspondingly the Lie algebra nature of screws, Karger , Rico et al. , Lerbet  and Müller [8, 9] derived closed form expressions of higher-order time derivatives of twist.
In this chapter, using the tensor calculus and the dual numbers algebra, a new computing method for studying the higher-order accelerations field properties is proposed in the case of the general rigid body motion. For the spatial kinematic chains, equations that allow the determination of the order accelerations field are given, using a Brockett-like formula. The crucial observation is that the order time derivative of twist of the terminal body in a kinematic chain can be determined by propagating the order time derivative of twists of the bodies in the chain, for . The results are coordinate-free and in a closed form.
2. Theoretical consideration on rigid body motion
The general framework of this chapter is a rigid body that moves with respect to a fixed reference frame . Consider another reference frame originated in a point that moves together with the rigid body. Let denote the position vector of point with respect to frame , its absolute velocity and its absolute acceleration.
Then the vector parametric equation of motion is:
The results of this section succinctly present the velocity and acceleration vector field in rigid body motion. These results lead to the generalization presented in the next section.
With the denotations that were introduced, the vector fields of velocities and accelerations are described by:
do not depend on the choice of point of the rigid body. They are called the
2.1 The velocity field in rigid body motion
It is described by:
The instantaneous angular velocity of the rigid body may be determined as . The major property that may be highlighted from Eq. (4) is that the velocity of a given point of the rigid may be computed when knowing the velocity tensor and the velocity invariant :
2.2 The acceleration field in rigid body motion
It is described by
The absolute acceleration of a given point of the rigid body may be computed when knowing the acceleration tensor and the acceleration invariant :
The instantaneous angular acceleration of the rigid body may be determined as:
The determinant of tensor is (see ): . It follows that if , then tensor
is invertible and its inverse is (see ):
It follows that if tensor is non-singular, then for an arbitrary given acceleration we may find a point of the rigid that has this acceleration. Its absolute position is given by (see also Eq. (8)):
Particularly, if is non-singular, then there exists a point of zero acceleration, named the
3. The vector field of the order accelerations
This section extends some of the previous considerations to the case of the order accelerations. We define the order acceleration of a point as:
For , it represents the velocity, and for , the acceleration. By derivation with respect to time successively in Eq. (2), it follows that:
does not depend on the choice of the point of the rigid body for which the acceleration is computed. Vector is named the
The next Theorem gives the fundamental properties of the vector field of the order accelerations.
for any point of the rigid body with the absolute position defined by vector .
3.1 Homogenous matrix approach to the field of order accelerations
The set of affine maps, , where is an orthogonal proper tensor and a vector in is a group under composition and it is called
One may remark that the following relations hold true:
We may extend now to , the set of the functions with the domain ℝ and the range . The parametric vector equation of the rigid body motion (1) may be rewritten with the help of a homogenous matrix function in like it follows:
From Eq. (25), it follows that:
and by making the computations and taking into account Eqs. (3) and (4) it follows that:
By using the previous considerations, it follows that Eq. (25) may be extended like:
Eq. (28) represents a unified form of describing the vector field of the order accelerations in rigid body motion. The matrix:
contains both the order acceleration tensor and the vector invariant . Eqs. (20) may be put in a compact form:
If follows that may be written as:
4. Symbolic calculus of higher-order kinematics invariants
We will present a method for the symbolic calculation of higher-order kinematics invariants for rigid motion.
Let be and , vector invariant, respectively, tensor invariant for the order accelerations fields. We denote by
and we have the following relationship of recurrence:
The pair of vectors is also known as
Let be the matrix ring
and the set of polynomials with coefficients in the non-commutative ring . A generic element of has the form
Since it follows the next outcome.
Thus, it follows:
the velocity field invariants
the acceleration field invariants
jerk field invariants
hyper-jerk (jounce) field invariants
4.1 Higher-order acceleration center and vector invariants of rigid body motion
Equation (16) may be written as
This shows us that the vector function
has the same value in every point of the rigid body under the general spatial motion, at a given moment of time
The invariant value of vector is obtained for and it is the n-th order acceleration of the point of the rigid body that passes the origin of the fixed reference frame at a given moment of time: . Eq. (46) becomes:
Let be be the adjugate tensor of uniquely defined by:
From Eq. (46), results another invariant
The value of this invariant is
In the specific case when tensor is non-singular (), from (47) results the position vector having an imposed n-th order acceleration :
In a particular case of the
Assuming that the tensor is non-singular, the previous relations lead to a new vector invariant that characterize the accelerations of n-th and m-th order ():
The value of this invariant is .
The problem of the determination the adjugate tensor of the n-th acceleration tensor and the conditions in which these tensors are inversable is, as the author knows, still an open problem in theoretical kinematics field. We will propose a method based on the tensors algebra that will give a closed form, coordinate- free solution, dependent to the time derivative of spatial twist.
The vector field of the higher-order acceleration is a non-stationary vector field. Differential operator
Let a tensor and we note and . The below theorem takes place.
Let the n-th order acceleration tensor,
The vectors and the symmetric tensors can be obtained with the below recurrence relation:
It follows that:
Velocity field: , ,
is singular for any . In this case,
is nonsingular if and only if . In this case
is nonsingular if and only if . In this case
In Eqs. (65) and (66), the following notation has been used:
then is inversible and
In the hypothesis (68), there is jounce center, determined by
5. Dual algebra in rigid body kinematics
In this section, we will present some algebraic properties for dual numbers, dual vectors and dual tensors. More details can be found in [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25].
5.1 Dual numbers
Let the set of real dual numbers to be denoted by
where is the real part of and the dual part. The sum and product between dual numbers generate a ring with zero divisors structure for .
Any differentiable function can be completely defined on such that:
Based on the previous property, two of the most important functions have the following expressions: ;
5.2 Dual vectors
In the Euclidean space, the linear space of free vectors with dimension 3 will be denoted by . The ensemble of dual vectors is defined as:
where is the real part of and the dual part. For any three dual vectors , the following notations will be used for the basic products:
The magnitude of , denoted by , is the dual number computed from
where is the Euclidean norm. For any dual vector , if then is called unit dual vector.
5.3 Dual tensors
An -linear application of into is called an Euclidean dual tensor:
Let be the set of dual tensors, then any dual tensor can be decomposed as , where are real tensors. Also, the dual transposed tensor, denoted by , is defined by
while , the determinant is
Orthogonal dual tensor maps are a powerful instrument in the study of the rigid motion with respect to an inertial and non-inertial reference frames.
Let the orthogonal dual tensor set be denoted by:
where is the set of special orthogonal dual tensors and is the unit orthogonal dual tensor.
Taking into account the Lie group structure of and the result presented in previous theorem, it can be concluded that any orthogonal dual tensor can be used globally parameterize displacements of rigid bodies.
The parameters and are called the
The Lie algebra of the Lie group is the skew-symmetric dual tensor set denoted by , where the internal mapping is .
The link between the Lie algebra , the Lie group , and the exponential map is given by the following.
is well defined and onto.
and is the inverse of Eq. (81).
Based on Theorems 6 and 7, for any orthogonal dual tensor , a dual vector can be computed and represents the
The form of implies that . The types of rigid displacements that can be parameterized by are:
roto-translation if ;
pure translation if ;
pure rotation if .
6. Higher-order kinematics in dual Lie algebra
Being the rigid body motion given by the following parametric equation in a given reference frame:
with is time variable.
The dual orthogonal tensor that describes the rigid body motion is [13, 24]:
In relation (87), the skew symmetric tensor associated to the vector is denoted by .
It can be easily demonstrated [14, 15], that:
The dual angular velocity for the rigid body motion (86) is given by (87):
It can be demonstrated that:
is the instantaneous angular velocity of the rigid body and
is the linear velocity of a point of the rigid body that coincides with the origin of the reference frame at that given moment.
The dual angular velocity completely characterizes the distribution of the velocity field of the rigid body. The pair () is called “the twist of the rigid body motion” [13, 14].
the dual velocity for a point localized in the reference frame by the position vector .
In (93), is the instantaneous angular velocity of the rigid body and is the linear velocity of the point. Using the next equation,
from (90), (92)–(94), results:
Consequently, is an invariant having the same value for any
Writing this invariant in two different points of the rigid body (noted with P and Q ), results that:
From (97), results:
Relation (97) is true for any P and Q.
Analogue with Eq. (95), the following invariants take place:
where we denoted
with the reduced acceleration, reduced jerk, respectively the reduced hyper-jerk (jounce), in a point given by the position vector :
In (100), , and are, respectively, the acceleration, the jerk, and the hyper-jerk (jounce), in a point given by the position vector
Analogue with Eq. (97) the following equations take place:
The lines corresponding to the dual vectors represent the loci, where the vectors have the minimum module value. Supplementary,
Interesting is the fact that for the plane motion because
All properties are extended for higher-order accelerations. The vector describes completely the helicoidally field of the order reduced accelerations, for :
In Eq. (103) denote the order of the dual reduced acceleration in a given point by the position vector .
It follows that the dual part of the order differentiation of
is the order reduced acceleration of that point of the rigid body that at the given time pass by the origin of the reference frame.
it follows that
with the following notations
for the order acceleration of the point given by the position vector and
for the order reduced acceleration of the same point the equation:
which proves the character of the helicoidally field of the order reduced accelerations field.
For , the relations between the order reduced acceleration and the n order acceleration from point O, the origin of the reference frame, are written
The invert of previous equation is written:
where is the polynomial with the coefficients in the ring of the second order Euclidean tensors and the polynomials follow the recurrence equation:
it follows successively
If we denote and by , , for the case of the velocities, accelerations, jerks and jounces, on obtain (Figure 1):
7. Higher-order kinematics of spatial chain using dual Lie algebra
Consider a spatial kinematic chain of the bodies where the relative motion of the rigid body with respect to is given by the proper orthogonal tensor . The relative motion properties of the body with respect to are described by the orthogonal dual tensor (Figure 2):
Instantaneous dual angular velocity (dual twist) of the rigid body in relation to the reference frame it will be given by the equation
It follows from (110) and (111) that:
Using the denotation
Eq. (118) will be written
where is the dual twist of the relative motion of the body in relation to the body observed from the body .
To determine the field of the order accelerations of a rigid body we have to determine the , .
We denote the order derivative of the relative dual twist , resolved in the body frame of .
In order to determine the order accelerations field of a rigid body we have to determine the , .
To compute we will use the following
Other equivalent forms of Eq. (127) are the following recursive formulas (Figures 3 and 4):
The previous equations are valid in the most general situation where there are no kinematic links between the rigid bodies .
The following identity can be proved:
where is the multinomial coefficient.
From Eq. (131), on obtain the closed form non-recursive coordinate-free formula:
8. Higher-order kinematics for general 2C manipulator
We’ll apply the general results obtained in the previous chapter for the particular case of four degrees of freedom 2C general manipulator. In this case the relative motions of three bodies are given, the spatial motion of the terminal body been described by dual orthogonal tensor as it follows:
In Eqs. (138) and (139), the dual angles and are four times differentiable functions, and unit dual vectors and being constant. To simplify the writing, we will denote:
According to the observations from Section 6, the vector field of the velocity, the acceleration, the jerk, the jounce is uniquely determined by the dual vectors . Taking into account Eq. (133), we will have:
Similarly, the results for six degrees of freedom general 3 C manipulator can be obtained, the calculus being a little longer.
The higher-order kinematics properties of rigid body in general motion had been deeply studied. Using the isomorphism between the Lie group of the rigid displacements and the Lie group of the orthogonal dual tensors , a general method for the study of the field of arbitrary higher-order accelerations is described. It is proved that all information regarding the properties of the distribution of high-order accelerations are contained in the n-th order derivatives of the dual twist of the rigid body. These derivatives belong to the Lie algebra associated to the Lie group
For the case of the spatial relative kinematics, equations that allow the determination of the n-th order field accelerations are given, using a Brockett-like formulas specific to the dual algebra. In particular cases the properties for velocity, acceleration, jerk, hyper-jerk (jounce) fields are given.
The obtained results interest the theoretical kinematics, jerk and jounce analysis in the case of parallel manipulations, control theory and multibody kinematics.
Everett JD. On the kinematics of a rigid body. The Quarterly Journal of Pure and Applied Mathematics. 1875; 13:33-66
Ball RS. Theory of Screws. London: Cambridge University Press Warehouse; 1900
Hunt KH. Kinematic Geometry of Mechanisms. New York: Oxford University Press; 1978
Phillips J. Freedom in Machinery: Introducing Screw Theory. Vol. 1. Cambridge: Cambridge University Press; 1984
Karger A. Singularity analysis of serial robot-manipulators. ASME Journal of Mechanical Design. 1996; 118(4):520-525
Rico JM, Gallardo J, Duffy J. Screw theory and higher order kinematic analysis of open serial and closed chains. Mechanism and Machine Theory. 1999; 34(4):559-586
Lerbet J. Analytic geometry and singularities of mechanisms. Zeitschrift für Angewandte Mathematik und Mechanik. 1999; 78(10):687-694
Müller A. Higher derivatives of the kinematic mapping and some applications. Mechanism and Machine Theory. 2014; 76:70-85
Müller A. An overview of formulae for the higher-order kinematics of lower-pair chains with applications in robotics and mechanism theory. Mechanism and Machine Theory. 2019; 142
Condurache D, Matcovschi M. Computation of angular velocity and acceleration tensors by direct measurements. Acta Mechanica. 2002; 153(3–4):147-167
Angeles J. The application of dual algebra to kinematic analysis, computational methods. Mechanical Systems. 1998; 161:3-31
Bokelberg EH, Hunt KH, Ridley PR. Spatial motion-I: Points of inflection and the differential geometry of screws Raumliche bewegung-I wendepunkte und die differentiale schraubengeometrie. Mechanism and Machine Theory. 1992; 27(1):1-15
Condurache D, Burlacu A. Orthogonal dual tensor method for solving the AX=XB sensor calibration problem. Mechanism and Machine Theory. 2016; 104:382-404
Condurache D, Burlacu A. Dual tensors based solutions for rigid body motion parameterization. Mechanism and Machine Theory. 2014; 74:390-412
Condurache D, Burlacu A. Dual lie algebra representations of the rigid body motion. In: AIAA/AAS Astrodynamics Specialist Conference; San Diego, CA; 2014
Condurache D. A Davenport dual angles approach for minimal parameterization of the rigid body displacement and motion. Mechanism and Machine Theory. 2019; 140:104-122
Condurache D. General rigid body motion parameterization using modified Cayley transform for dual tensors and dual quaternions. In: The 4th Joint International Conference on Multibody System Dynamics; Montreal; 2016
Condurache D. Higher-order accelerations on rigid bodies motions. A tensors and dual lie algebra approach. Acta Technica Napocensis – Series: Applied Mathematics, Mechanics and Engineering. 2018; 61(1):29-38
Condurache D. Higher-order kinematics of rigid bodies. A tensors algebra approach. In: Kecskeméthy A, Flores FG, Carrera E, Elias DA, editors. Interdisciplinary Applications of Kinematics. Mechanisms and Machine Science. Vol. 71. Cham: Springer; 2019. pp. 215-225
Condurache D. Higher-order relative kinematics of rigid body motions: A dual lie algebra approach. In: Lenarcic J, Parenti-Castelli V, editors. Advances in Robot Kinematics. Vol. 8. 2018. pp. 83-91
Condurache D. Higher-order acceleration centers and kinematic invariants of rigid body. In: The 5th Joint International Conference on Multibody System Dynamics; Lisbon; 2018
Condurache D. Higher-order Rodrigues dual vectors. Kinematic equations and tangent operator. In: The 5th Joint International Conference on Multibody System Dynamics; Lisbon; 2018
Pennestrı E, Valentini PP. Linear dual algebra algorithms and their application to kinematics. In: Multibody Dynamics, Computational Methods and Applications. Dordrecht: Springer; 2009. pp. 207-229
Samuel A, McAree P, Hunt K. Unifying screw geometry and matrix transformations. The International Journal of Robotics Research. 1991; 10(5):454-472
Veldkamp GR. Canonical systems and instantaneous invariants in spatial kinematics. Journal of Mechanisms. 1967; 2(3):329-388