Next Article in Journal
Magnetohydrodynamics Stagnation-Point Flow of a Nanofluid Past a Stretching/Shrinking Sheet with Induced Magnetic Field: A Revised Model
Next Article in Special Issue
Analysis of a Poro-Thermo-Viscoelastic Model of Type III
Previous Article in Journal
New Characterizations of the Clifford Torus and the Great Sphere
Previous Article in Special Issue
Hysteretically Symmetrical Evolution of Elastomers-Based Vibration Isolators within α-Fractional Nonlinear Computational Dynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Synthesis on the Acceleration Energies in the Advanced Mechanics of the Multibody Systems

by
Iuliu Negrean
* and
Adina-Veronica Crișan
*
Faculty of Machine Building, Department of Mechanical Systems Engineering, Technical University of Cluj-Napoca, 400641 Cluj-Napoca, Romania
*
Authors to whom correspondence should be addressed.
Submission received: 16 July 2019 / Revised: 12 August 2019 / Accepted: 22 August 2019 / Published: 27 August 2019
(This article belongs to the Special Issue Symmetry in Applied Continuous Mechanics)

Abstract

:
The present paper’s objective is to highlight some new developments of the main author in the field of advanced dynamics of systems and higher order dynamic equations. These equations have been developed on the basis of the matrix exponentials which prove to have undeniable advantages in the matrix study of any complex mechanical system. The present paper proposes some new approaches, based on differential principles from analytical mechanics, by using some important dynamics notions, regarding the acceleration energies of the first, second and third order. This study extended the equations of the higher order, which provide the possibility of applying the initial motion conditions in the positions, velocities and accelerations of the first and second order. In order to determine the time variation laws for the generalized variables, the driving forces and acceleration energies of the higher order are applied by the time polynomial functions of the fifth order. According to inverse kinematics also named control kinematics of the robots, the applications of polynomial functions lead to the kinematic control functions of mechanical motions, especially the transitory motions. They influence the dynamic behavior of multibody systems, in which robot structures are included.

1. Introduction

Until the present, the dynamic behavior of rigid body and of the multibody systems can be studied by considering the fundamental theorems from dynamics or the differential principles of analytical mechanics [1,2,3,4]. Both are based on angular momentum, the mechanical work and kinetic energy, as well as on their fundamental theorems. In mechanics, the defining expressions for the advanced notions are based on fundamental input parameters that are represented by kinematic parameters along with their differential transformations, as well as by the mass properties highlighted by the inertia tensor and corresponding generalized variation laws.
The advanced concepts include higher-order accelerating energies known in the scientific literature [5,6,7] as kinetic energies of higher-order accelerations. The kinetic energy is used, in Newtonian mechanics, as a central function of Lagrange Euler equations. The existence of higher-order accelerations has been proven by the fact that they occur every time the bodies/multi-body systems are subjected to fast and sudden motions or, when particularly, this type of excitation produces vibration modes with more frequencies of resonance. The acceleration itself cannot produce vibration being just a static load. In fact, the sudden motions are generated by higher order accelerations. For example, during the operation of a cam and cam follower mechanism at high speeds and accelerations, the jumping off of the cam follower from the camshaft causing accidents may occur. Another example is the case of manufacturing processes that require rapid changes in the acceleration of a cutting tool. This leads to a quick tool wear and results in poor quality of the manufactured surfaces [8,9]. The two situations presented above are generated by the effect of higher order accelerations. These effects have to be also considered each time when a mechanical system is in a transient motion phase (take-off/landing, acceleration/deceleration, etc.). The novelty of the paper consists in presenting a new approach regarding the dynamic modeling of rigid bodies and multibody mechanical systems. Based on scientific literature [10,11,12,13,14,15,16,17,18,19,20], the input kinematic and dynamic parameters have been formulated. These parameters are compulsorily included in the dynamic equations of the higher order, corresponding to fast and sudden motions. This type of motion characterizes the multibody systems, a category where the robots are also included. The important formulations regarding the acceleration energies of the first, second, third and higher orders are also presented in explicit and matrix forms. The expressions for the acceleration energies of the higher order were determined, among others, based on matrix exponential functions [14,19]. Although at first sight they seem difficult to put into practice, they proved to have undisputable advantages in terms of their use in the kinematic and dynamic study of complex mechanical systems.
These expressions are implemented into the differential equations of the higher order, thus resulting the time variations of the generalized forces, which influences the dynamic behavior of rigid body or of multibody systems. To illustrate the validity of the proposed mathematical model, an experimental study, in which the higher order accelerations are highlighted, is performed.

2. Input Parameters in Advanced Dynamic Modeling

2.1. Position and Orientation Parameters

In Newtonian mechanics, the dynamics of the rigid bodies are studied based on the assumption that they do not deform under the action of the external forces.
To study the dynamical behavior of a rigid body, its linear and angular position in the Cartesian space must be known at each moment during the motion. To define the linear and angular position, [13,20], the last one being also known as orientation, the simplest mechanical model, which is the material point, is analyzed (Figure 1).
The following notations are considered:
χ = { u ; v ; w } ;    χ 0 = { u 0 ; v 0 ; w 0 }      w h e r e ​​​   u = { x ; y ; z } ; v = { y ; z ; x } u ; w = { z ; x ; y } v
χ ¯ = { u ¯ ;    v ¯ ;    w ¯ } ;    χ ¯ 0 = { u ¯ 0 ;    v ¯ 0 ;    w ¯ 0 }      w h e r e      u ¯ = { i ¯ ; j ¯ ; k ¯ } ;      v ¯ = { j ¯ ; k ¯ ; i ¯ } u ¯ ;      w ¯ = { k ¯ ; i ¯ ; ​  j ¯ } v ¯
δ χ = { α χ ; β χ ; γ χ } ;   cos δ χ = c δ χ ; sin δ χ = s δ χ
O 0 x 0 y 0 z 0 { 0 } ;   O 0 x 0 y 0 z 0 { 0 } ;   O x y z { S }
In the Expression (1), the coordinates or axes of the Cartesian frame, symbols from (2) representing the unit vectors, and (3) the angles and the direction cosines are defined.
The dynamic study of the rigid solid is carried out, according to Figure 1, by considering two reference systems: O 0 x 0 y 0 z 0 = { 0 } , which is considered fixed and O x y z = { S } which is a mobile reference system, invariably linked to the body, having the origin in an arbitrary point O belonging to the rigid body. The system O x 0 y 0 z 0 = { 0 } , represented in the same figure, is a system which originated in point O , and whose orientation is maintained constant all throughout the movement and is identical with the orientation of the fixed system, { 0 } , that is { 0 } O R { 0 } O R .
The geometrical state of any material point (for example a certain point O ) characterizes the position, according to Figure 1a. The position is usually defined by means of the position vector:
r ¯ 0 = [ x 0 y 0 z 0 ] T ,    r e l a t i v e t o { 0 } f r a m e
The three linear coordinates from (5) are independent in the case of a free material point, and represents the degrees of freedom (d.o.f.). Further, the study is extended to a vector or an axis belonging to the Cartesian frame (see Figure 1b). This geometrical state is called orientation (angular position). Orientation (angular position) is defined through unit vectors. For any given unit vector χ ¯ { S } , in relation to one of the frames { 0 } / { 0 } , the orientation is defined by the direction cosines:
χ ¯ = χ ¯ T ( i ¯ 0 j ¯ 0 k ¯ 0 ) = ( c α χ c β χ c γ χ ) ( c α c β c γ ) χ ,   where   χ ¯ T χ ¯ = c 2 α χ + c 2 β χ + c 2 γ χ = 1
According to linear algebra, in expression (6), the symbol χ ¯ T defines the transposed matrix.
Further, considering the second expression from (6), the orientation of any vector or axis is characterized by two independent angles. The above geometrical aspects are extended on an orthogonal and right oriented frame (see Figure 1b, O u v w O x y z { S } ) relative to { 0 } . In this case, the geometrical state is defined by position and orientation. The position is defined by means of position vector (5), while for orientation the rotation matrix is considered [10]:
[ R ] S 0 = [ i ¯ j ¯ k ¯ ] = [ ( c α c β c γ ) x ( c α c β c γ ) y ( c α c β c γ ) z ] .
The resultant orientation matrix defined with (7), contains the unit vectors of { S } frame relative to { 0 } . The rotation matrix or the matrix of direction cosines describes the orientation of each axis of the mobile frame attached to the rigid body, relative to a fixed frame. Between the nine direction cosines contained by the rotation matrix, six mathematical relations can be established. Thus, the result is that the orientation of a mobile frame with respect to a fixed frame is defined by up to three independent parameters, which are defining the orienting angles. Therefore, the resultant orientation of a frame { S } with respect to another frame, { 0 } / { 0 } is characterized by three independent orientation angles (three d.o.f), [13]:
ψ ¯ ¯ ( t ) ( 3 x 1 ) = [ α u ( t ) β v ( t ) γ w ( t ) ] T
The angles from (8) are the components of the column matrix of orientation ψ ¯ ¯ ( t ) , and they describe, from a geometrical point of view, dihedral angles between two geometrical planes:
χ 0 = { u 0 ; v 0 ; w 0 } = c s t f i x e d p l a n e { 0 } / { 0 }   and   χ = { u ; v ; w } = 0 m o b i l e p l a n e { S }
Physically, the three angles defined with (8) expresses a simple rotation around one of the axes of the Cartesian reference system: χ = { u ; v ; w } . Based on research from [10,11,12,13,14], when combining the three simple rotations, there is a result of twelve sets of orientation angles (8). Taking χ = { x ; y ; z } , the expressions of definition are further developed for the three simple rotation matrices symbolized as:
R ( χ ¯ ; δ χ ) = { R ( x ¯ ; α x ) ; R ( y ¯ ; β y ) ;   ​​ R ( z ¯ ; γ z ) } .
The paper proposes the following form of mathematical representation for the generalized matrix:
R ( χ ¯ ; δ χ ) = { R ( x ¯ ; α x ) ; R ( y ¯ ; β y ) ;   ​​ R ( z ¯ ; γ z ) } = [ c ( δ χ Δ y z ) s ( δ χ Δ z ) s ( δ χ Δ y ) s ( δ χ Δ z ) c ( δ χ Δ z x ) s ( δ χ Δ x ) s ( δ χ Δ y ) s ( δ χ Δ x ) c ( δ χ Δ x y ) ] ,
w h e r e      Δ u v { χ = { u ; v } } = { Δ y z ; Δ z x ; Δ x y } = { { 1 0 }    i f    δ χ = { { ( β y ; γ z ) α x } ; { ( γ z ; α x ) β y } ; { ( α x ; β y ) γ z } } } ,
a n d             Δ u { χ = u } = { Δ x ;    Δ y ;    Δ z } = 1 Δ v w . { χ = { v ; w } }
By substituting (11) and (12) in the generalized expression (10), the simple rotation matrices (9) are obtained. The generalized matrix can be written in a new formulation as follows:
R ( χ ¯ ; δ χ ) = D i a g [ Δ u v ¯ ¯ ] + [ Δ u ¯ ¯ × ] ,
where   Δ u v ¯ ¯ ( 3 × 1 ) = [ c ( δ χ Δ y z )    c ( δ χ Δ z x )    c ( δ χ Δ x y ) ] T ,
and   Δ u ¯ ¯ ( 3 × 1 ) = [ s ( δ χ Δ x )    s ( δ χ Δ y )      s ( δ χ Δ z ) ] T ,
In the expression (13), the symbol [ Δ u ¯ ¯ × ] defines the skew-symmetric matrix associated to (15), while D i a g [ Δ u v ¯ ¯ ] represents the diagonal matrix, written below:
D i a g [ Δ u v ¯ ¯ ] ( 3 x 3 ) = [ c ( δ χ Δ y z ) 0 0 0 c ( δ χ Δ z x ) 0 0 0 c ( δ χ Δ x y ) ] .
By applying (14) and (15), an expression identical with the classical Rodriguez formula is obtained [17]:
R ( χ ¯ ; δ χ ) = χ ¯ χ ¯ T ( 1 c δ χ ) + I 3 c δ χ + ( χ ¯ × s δ χ ) .
According to research from [3,4,5,6,7,8,9,10], the three simple rotations from (8) are performed either around the moving axes or fixed axes belonging to { S } or { 0 } / { 0 } . The resultant rotation matrix, symbolized as [ R ] S 0 , defines the orientation of the system { S } ( attached to the rigid body), with respect to a fixed system { 0 } , and χ = { u ; v ; w } is represented in the following matrix form:
[ R ] S 0 = R ( u ¯ ; α u ) R ( v ¯ ; β v ) ​​  R ( w ¯ ; γ w ) ,
[ R ] S 0 = { χ ¯ ; δ χ } R ( χ ¯ ; δ χ ) = { χ = { u ; v ; w } } exp [ χ ¯ × δ χ ] = { u ; v } { I 3 Δ u v ¯ ¯ + [ Δ u ¯ ¯ × ] } .
Using the research in the field of matrix exponentials [10,14], the following form for the resultant rotation matrix (18) is written as:
[ R ] S 0 = { χ ¯ ; δ χ } R ( χ ¯ ; δ χ ) = { χ = { u ; v ; w } } exp [ χ ¯ × δ χ ] = exp [ u ¯ × α u ] exp [ v ¯ × β v ] exp [ w ¯ × γ w ]
The above mathematical Expressions (5) and (8), which define the position and orientation in the case of a rigid body, are written in a general form. Based on the definitions from the first paragraph of this section, a rigid body consists of an infinite number of material points and also, an infinite number of geometrical axes, parallel and perpendicular one to another, characterized by a continuous distribution in the entire volume of the rigid ( S ) [14,15]. The rigid body is also made up of an infinite number of assemblies consisting of three orthogonal geometrical plans, continuously distributed in the entire volume of the rigid body.
In terms of geometry, in order to define a right oriented reference frame having the origin in an arbitrary point O of the rigid body, a single ensemble composed of three orthogonal geometrical plans is chosen. According to (4) and Figure 2, this reference frame is symbolized as O x y z { S } . The reference frame { S } is linked to the body due to its rigid nature. Expressions (5) and (8) define the position and orientation of this frame.
Whereas the present paper proposes a study on advanced dynamics, it considers two material points belonging to a rigid body, so that the following conditions are met: M O and C { O ; M } . In this case, the expressions defining the position can be written as time functions, according to [12]:
r ¯ M ( t ) = r ¯ 0 ( t ) + ρ ¯ M ( t ) = r ¯ 0 ( t ) + [ R ] S 0 ( t ) S ρ ¯ M ,
r ¯ C ( t ) = r ¯ 0 ( t ) + ρ ¯ C ( t ) = r ¯ 0 ( t ) + [ R ] S 0 ( t ) S ρ ¯ C ,   where   ρ ¯ M ( t ) ρ ¯ C ( t )   and   r ¯ M ( t ) r ¯ C ( t )
Based on (18) and (19), the position equation is written, by considering the matrix exponentials:
r ¯ M ( t ) = r ¯ 0 ( t ) + { χ = { u ; v ; w } } exp [ χ ¯ × δ χ ] S ρ ¯ M = r ¯ 0 ( t ) + { exp [ u ¯ × α u ] exp [ v ¯ × β v ] exp [ w ¯ × γ w ] } S ρ ¯ M .
The position equation for any material point from the rigid body can be determined only when the position r ¯ 0 ( t ) and orientation [ R ] S 0 ( t ) of the moving frame O x y z { S } are well known. By analyzing the Expressions (21) and (22), the result is that the orientation is invariant for all the points of the rigid solid. Thus, by considering the geometrical aspects, the body is substituted by the moving frame O x y z { S } . This is defined from a geometrically point of view by means of the six independent parameters (six d.o.f.), which are included in the following symbol:
X ¯ ¯ ( 6 × 1 ) ( t ) = [ r ¯ 0 ( t ) ψ ¯ ¯ ( t ) ] = [ [ x 0 ( t ) y 0 ( t ) z 0 ( t ) ] T [ α u ( t ) β v ( t ) γ w ( t ) ] T ] = [ q i ( t ) δ i ,      i = 1 6 w h e r e    δ i = { ( 1 Δ i ) , f o r    q i l i n e a r , Δ i , f o r    q i a n g u l a r } ]
where q j ( t ) is the generalized coordinate and Δ i = { ( 1 , q i a n g u l a r ) , ( 0 , q i l i n e a r ) } is an operator which highlights the type of driving joint. In this paper, the symbol X ¯ ¯ ( 6 × 1 ) ( t ) , is substituted as:
X ¯ ¯ ( 6 × 1 ) ( t ) = θ ¯ ( t ) = [ q i ( t ) , i = 1 6 ] T , w h e r e θ ¯ θ ¯ ( 0 )
Considering Expression (24), the following notations are implemented:
{ θ ¯ ( t ) ; θ ¯ ˙ ( t ) ;    θ ¯ ¨ ( t ) ; ; θ ¯ ( m ) ( t ) } = { q i ( t ) ; q ˙ i ( t ) ; q ¨ i ( t ) ; ; q i ( m ) ( t ) i = 1 6 ,    m 1 } .
where m 1 represents the order of the time derivative.
In advanced mechanics, instead of (8), in the case of the angular vector of orientation, the following expression of definition is used:
ψ ¯ ( t ) = J 0 ψ [ α u ( t ) β v ( t ) γ w ( t ) ] ψ ¯ ¯ ( t ) = ψ ¯ [ q j ( t ) Δ j ;      j = 1 k = 6 ,      t ] ,
J 0 ψ [ α u ( t ) β v ( t ) γ w ( t ) ] = [ u ¯ R ( u ¯ ; α u ) v ¯ R ( u ¯ ; α u ) R ( v ¯ ; β v ) w ¯ ] .
where J 0 ψ represents the angular transfer matrix which is a function of the orienting angles.
In mechanics, the position and orientation of a moving frame O x y z { S } , with respect to another frame, for example a fixed frame { 0 } , is represented in a matrix form, according to [5,8] by means of homogeneous transformations which are developed using matrix exponentials:
[ T ] S 0 ( 4 × 4 ) ( t ) = [ [ R ] S 0 ( t ) p ¯ ( t ) 0 0 0 1 ] = [ { χ = { u ; v ; w } } exp [ χ ¯ × δ χ ] p ¯ ( t ) 0           0           0 1 ] ,
p ¯ ( t ) = { χ = { u ; v ; w } } { { χ = { u ; v ; w } } exp [ χ ¯ × δ χ ] } b ¯ χ + { { χ = { u ; v ; w } } exp [ χ ¯ × δ χ ] } p ¯ ( 0 ) Δ p ,   where   Δ p = { { 0 ; p ¯ = r ¯ 0 } ; { 1 ; p ¯ = r ¯ M } } ,
b ¯ χ = [ χ ¯ ( 0 ) χ ¯ ( 0 ) T ( δ χ s δ χ ) + I 3 s δ χ + ( χ ¯ ( 0 ) × ) ( 1 c δ χ ) ] [ p ¯ ( 0 ) × χ ¯ ( 0 ) Δ χ ] .
Expression (30) defines the position vector of the mobile system { S } relative to fixed system { 0 } , while (31) represents a position vector, which according to [11,20], is defined by means of homogenous coordinates as a function of screw parameters. The conclusions and expressions of definition, synthetically disseminated in this introductory section, are compulsory applied in the advanced kinematics and dynamics of mechanical systems. For example, a kinetic element from this mechanical structure of a robot is considered (Figure 3).
In case of both rigid bodies and multibody systems, the dynamical study of current and sudden motions is performed. The study is based on the differential principles from analytical dynamics of systems, as well as on the advanced notions from the dynamics of solid bodies: Momentum and angular momentum, kinetic energy, acceleration energies of different orders and their absolute time derivatives of the high order.

2.2. The Kinematic Parameters of Higher Order

The parametric equations of motion are defined according to (24). The absolute position Equation (21) shows the variable distribution from a point to another material point of the body. By applying the first order time derivative on Expression (21), the result is:
r ¯ ˙ M ( t ) = r ¯ ˙ 0 ( t ) + ρ ¯ ˙ M ( t ) = r ¯ ˙ 0 ( t ) + [ R ˙ ] S 0 ( t ) ρ ¯ S M = r ¯ ˙ 0 ( t ) + [ R ˙ ] S 0 ( t ) [ R ] S 0 T ( t ) [ R ] S 0 ( t ) ρ ¯ S M ,
The following property is mandatory for the motion equations defined with (32):
[ R ˙ ] S 0 ( t ) [ R ] S 0 T ( t ) = ( ω ¯ × )
where considering [1,10], ( ω ¯ × ) is the skew symmetric matrix associated to the angular velocity vector.
The time derivatives of the higher order, applied on the property (32), are defined as follows:
[ R ] S 0 ( k ) [ R ] S 0 T = j = 1 k 3 d j d t j { [ R ] S 0 [ k ( j + 1 ) ] [ R ] S 0 T [ ω ¯ × ] } δ k 1 + d k 2 d t k 2 { [ ω ¯ × ] [ ω ¯ × ] } δ k + [ ω ¯ × ] ( k 1 ) + [ R ] S 0 ( k 1 ) [ R ] S 0 T [ ω ¯ × ] ,   w h e r e      k 2 ,    a n d δ k = { ( 0 ,      w h e n    k 2 )    a n d    ( 1 ,      w h e n    k > 2 ) }
In advanced kinematics and dynamics, the time derivatives of the higher order for rotation matrices and position vectors must be used according to [10,11,12,13] as follows:
[ R ] S 0 [ q j ( t ) Δ j ] ( k ) = j = 1 6 { q j ( m ) { [ R ] S 0 } ( m ) Δ j q j ( k ) } + j = 1 6 r = 1 k 1 { p = 1 r ( k p ) p ! { p ! m ! ( m + p ) ! q j ( m ) { [ R ] S 0 ( m + p ) } Δ j q j ( k p ) } } ,
v ¯ M ( k 1 ) = d k r ¯ M ( t ) d t k = j = 1 6 [ r ¯ M ( m ) q j ( m ) q j ( k ) ] + j = 1 6 r = 1 k 1 { p = 1 r ( k p ) p ! [ p ! m ! ( m + p ) ! r ¯ M ( m + p ) q j ( m ) q j ( k p ) ] }
where   k 1 ;      k = { 1 ; 2 ; 3 ; 4 ; 5 ; }    a n d    m ( k + 1 ) ;    m = { 2 ; 3 ; 4 ; 5 ; }
The symbols: ( k ) and ( m ) represent the order of the time derivatives applied on expressions (35) and (36).
The expressions of definition for angular velocities, and then for the angular accelerations are:
ω ¯ ˙ [ α u ( t ) β v ( t ) γ w ( t ) ] = d d t { ω ¯ [ α u ( t ) β v ( t ) γ w ( t ) ] } = d d t { α ˙ u ( t ) { exp [ 0 ] } u ¯ ( 0 ) + + β ˙ v ( t ) { exp [ u ¯ ( t ) × α u ( t ) ] } v ¯ ( 0 ) + γ ˙ w ( t ) { exp [ u ¯ ( t ) × α u ( t ) ] exp [ v ¯ ( t ) × β v ( t ) ] } w ¯ ( 0 ) } .
By analyzing (38), it is noticed that the expressions of the kinematic parameters that define the resultant rotation motion are written based on the matrix exponential functions.
Both velocities and accelerations of the higher order are established based on the following vectors:
r ¯ 0 ( t ) + exp [ { χ = { u ; v ; w } } exp [ χ ¯ × δ χ ] ] ρ ¯ S C = r ¯ C [ q j ( t ) ;      j = 1 k = 6 ,      t ] = r ¯ C ( t ) ,
ψ ¯ ( t ) = ψ ¯ [ q j ( t ) Δ j ; j = 1 k = 6 ,    t ] = α u ( t ) { exp [ 0 ] } u ¯ ( 0 ) + β v ( t ) { exp [ u ¯ ( t ) × α u ( t ) ] } v ¯ ( 0 ) + + γ w ( t ) { exp [ u ¯ ( t ) × α u ( t ) ] exp [ v ¯ ( t ) × β v ( t ) ] } w ¯ ( 0 )
where the orientation vector is defined according to the following expressions: (8), (19), (27) and (28).
An essential component (28) included in (40), known as the angular transfer matrix, can be defined as a function of a set of the orientation angles. Considering (39) and (40), it is observed that they are functions of generalized variables (24) and (25). Actually, the six generalized variables are the independent parameters of the position and orientation from (24). Using the research from [10,11,12,13,14], the on-time functions associated to position (39) and orientation (40) vectors and the differential properties which are compulsory and applied in advanced kinematics and dynamics, have been developed:
d k 1 d t k 1 ( r ¯ C q j ) = ( k 1 ) ! m ! ( m + k 1 ) ! r ¯ C ( m + k 1 ) q j ( m ) ,      k 1 ;      m ( k + 1 ) ;    m = { 2 ; 3 ; 4 ; 5 ; } ,
d k 1 d t k 1 ( ψ ¯ i q j Δ j ) = ( k 1 ) ! m ! ( m + k 1 ) ! ε ¯ ( m + k 3 ) q j ( m ) Δ j = ( k 1 ) ! m ! ( m + k 1 ) ! ψ ¯ ( m + k 1 ) q j ( m ) Δ j ,  
ε ¯ ( k 1 ) ( t ) = ω ¯ ( k ) ( t ) = j = 1 k = 6 d k 1 d t k 1 [ ψ ¯ ( t ) ( m ) q j ( m ) Δ j q ¨ j ( t ) ] + 1 m + 1 j = 1 k = n d k 1 d t k 1 [ ψ ¯ ( m + 1 ) q j ( m ) Δ j q ˙ j ( t ) ] = ψ ¯ ( k + 1 ) ( t ) ,
a ¯ C ( k 1 ) ( t ) = v ¯ C ( k ) ( t ) = j = 1 k = 6 d k 1 d t k 1 [ r ¯ C ( t ) ( m ) q j ( m ) q ¨ j ( t ) ] + j = 1 k = 6 d k 1 d t k 1 [ 1 m + 1 r ¯ C ( t ) ( m + 1 ) q j ( m ) q ˙ j ( t ) ] = r ¯ C ( k + 1 ) ( t ) ,
The symbols k and m from (41) and (42) define the deriving orders with respect to time.
The expression (44) defines the angular and linear accelerations of the higher order for a rigid solid in general motion. The input expressions and parameters of the higher order are compulsory and applied in the definition of the dynamic notions of the higher order, such as the momentum, angular momentum, kinetic energy and acceleration energy of the higher order.
They have been included in the dynamics theorems which applies to current and sudden motions.
The next step consists of establishing the time variation law for the generalized coordinates and generalized variables of the higher order. Considering [9], the result is that all parameters from advanced kinematics are functions of generalized variables as well as of their time derivatives and they can be developed by using polynomial interpolation functions [13,20]. The following general form for the polynomial interpolation functions of the higher order is used:
q j i ( τ ) ( m p ) = ( 1 ) p ( τ i τ ) p + 1 t i ( p + 1 ) ! q j i 1 ( m ) + ( τ τ i 1 ) p + 1 t i ( p + 1 ) ! q j i ( m ) + δ p k = 1 p τ p k ( p k ) ! a j i k ,
{ w h e r e      p = 0 m ,      m d e r i v i n g o r d e r , m 2 ,    m = 2 , 3 , 4 , 5 , ,      δ p = { ( 0 ,      p = 0 ) ; ( 1 ;      p 1 ) } j = 1 n ,      ( degrees of freedom - ( d . o . f . ) ) ;      i = 1 s ,    ( intervals of motion    trajectories ) τ actual time variable ,          t i = τ i τ i 1 ( time to each trajectory interval ) }
For every trajectory interval ( i = 1 s ) , the number of unknowns is ( m + 1 ) , and their significance is:
{ ( a j i k )    f o r    k = 1 m ; ​        a n d ( q j i 1 ( m ) )    f o r    i = 2 s w h e r e ( a j i k ) integration constants ,      and      ( q j i 1 ( m ) ) generalized accelerations of ( m ) order } ,
To determine the unknowns, (50) is required to apply geometrical and kinematical constraints:
{ ( τ 0 ) q j 0 ( m p ) ,      p = 0 m ;          ( τ s ) { q j s ( m ) ,      q j s } ( τ i ) { q j i ( 2 ) generalized accelerations    q j i ( τ + ) ( m p ) = q j i + 1 ( τ ) ( m p ) ,           p = 0 m ,      continuity conditions all conditions are applied to each ( τ i ) ,      where      i = 1 s 1 ,
The kinematical constrains that can be applied in the case of polynomials of the higher order are:
( τ 0 ) { h 0 k = q j 0 k ;    v 0 k = q ˙ j 0 k ;    a 0 k = q ¨ j 0 k ; a ¨ 0 k = q j 0 k } ,          ( τ n ) { h n k = q j n k ;    a ¨ n k = q j n k } ,
( τ i ) { a i k = q ¨ j i k ;    i = 1 n 1 ;      i k ;    k = 1 m { h i k ( t + ) = h i + 1 k ( t ) ;      v i k ( t + ) = v i + 1 k ( t ) ;    a ˙ i k ( t + ) = a ˙ i + 1 k ( t ) } w h e r e      k = 1 m      i s t h e s e q u e n c e o f t h e    w o r k i n g p r o c e s s
Finally, Expression (45) is substituted in the advanced notions of kinematics and dynamics.

2.3. The Mass Properties

In order to find an exact geometrical solution, it is considered that the rigid body consists of an infinite number of material particles for which the distance between two points is maintained constant, regardless of the forces that act upon it. The elementary particles, characterized by elementary and infinitesimal mass, have a continuous distribution in the geometrical shape of the solid body.
When the density is constant inside the rigid structure, the solid rigid is homogeneous. When the integration limits applied on the geometrical outline are well defined, it results in a homogeneous body with a simple or regular shape. In this last case, the geometrical and mass integrals are applied.
An essential aspect in advanced dynamics is represented by mass properties [13] The mass and the position of the mass center is determined in relation with O 0 x 0 y 0 z 0 { 0 } frame, as presented:
M = d m ;    w h e r e d m = { ρ σ d σ } ; σ = { V ; A ; L } ,
ρ ¯ C ( t ) = ρ ¯ M ( t ) d m d m = ρ ¯ M ( t ) d m M ,
In the expression (51), ρ σ defines the density of the material, where the symbol σ represents, by case, the volume (V), the area (A) or the length (L) of the rigid body subjected to study
The position of the mass center, with respect to the O 0 x 0 y 0 z 0 { 0 } reference frame, is expressed at first in a classical form, and then by using matrix exponentials as:
r ¯ C ( t ) = r ¯ M ( t ) d m d m = r ¯ M ( t ) d m M = r ¯ 0 ( t ) + ρ ¯ C ( t ) = r ¯ 0 ( t ) + [ R ] S 0 ( t ) ρ ¯ S C = = r ¯ 0 ( t ) + exp [ { χ = { u ; v ; w } } exp [ χ ¯ × δ χ ] ] ρ ¯ S C = r ¯ 0 ( t ) + { exp [ u ¯ × α u ] exp [ v ¯ × β v ] exp [ w ¯ × γ w ] } ρ ¯ S C .
By applying the time derivative on (53), the linear velocity and acceleration of the mass center are determined:
v ¯ C ( t ) = r ¯ ˙ C ( t ) = r ¯ ˙ 0 ( t ) + ρ ¯ ˙ C ( t ) = r ¯ ˙ 0 ( t ) + [ R ˙ ] S 0 ( t ) [ R ] S 0 T ( t ) [ R ] S 0 ( t ) ρ ¯ S C = = v ¯ 0 ( t ) + ω ¯ ( t ) × ρ ¯ C ( t ) = d d t { r ¯ 0 ( t ) + { exp [ { χ = { u ; v ; w } } exp [ χ ¯ × δ χ ] ] } ρ ¯ S C } = = d d t { r ¯ 0 ( t ) + { exp [ u ¯ × α u ] exp [ v ¯ × β v ] exp [ w ¯ × γ w ] } ρ ¯ S C } .
Using classical Formulations (52) and (53), the result is the expressions for the acceleration of mass center:
a ¯ C = v ¯ ˙ C ( t ) = v ¯ ˙ 0 ( t ) + d d t [ ω ¯ ( t ) × ρ ¯ C ( t ) ] = a ¯ 0 + ε ¯ × ρ ¯ C + ω ¯ × ω ¯ × ρ ¯ C
In dynamic modeling, the following expression is necessary:
d d t ( ρ ¯ M ( t ) ) = d d t ( ρ ¯ C ( t ) + r ¯ ( t ) ) = d d t ( [ R ] S 0 ( t ) ( ρ ¯ S C + r ¯ S ) ) = = ω ¯ ( t ) × ρ ¯ M ( t ) = ω ¯ ( t ) × ρ ¯ C ( t ) + ω ¯ ( t ) × r ¯ ( t ) .
In the case of rigid bodies involved in rotation motion, the inertia property is highlighted by the mechanical moments of inertia [12,13,14]. According to Figure 2, the position of the elementary mass ( d m ) relative to the mass center is defined by means of the vector: r ¯ = r ¯ ( t ) .
Considering (52), the following mass property is obtained:
r ¯ ( t ) d m = [ R ] S 0 ( t ) r ¯ S d m = 0 .
Based on research [12,13], the inertial tensor and its variation law relative to the concurrent frames applied in the mass center: { S } a n d { 0 } O R { 0 } O R is established, according to:
I S = ( r ¯ × ) ( r ¯ × ) T d m = [ R ] S 0 I S S [ R ] S 0 T ,
where I S S is inertial tensor axial and centrifugal of the body ( S ) in with relation with { S } frame, applied in the mass center ( C ) , having the property: { S } O R { S } O R :
I S S = ( r ¯ S × ) ( r ¯ S × ) T d m ,
Since the dynamic study refers to the absolute motion, the use of the following expressions of the inertial tensor are required:
I S = [ R ] S 0 I S S [ R ] S 0 T = I S C + I S , ( ρ ¯ S M × ) ( ρ ¯ S M × ) T d m = I S S   and   I S C = M ( ρ ¯ C × ) ( ρ ¯ C × ) T .
The matrix expression (60) represents the generalized variation law of the inertial tensor axial and centrifugal with respect to the frame { 0 } . Further, in the expression (60), I S C is the inertia matrix axial and centrifugal of the mass center relative to { 0 } . The inertial tensor axial and centrifugal is defined with respect to the absolute frame { 0 } O R = { 0 } O R , as presented below:
I S = ( r ¯ M × ) ( r ¯ M × ) T d m = M ( r ¯ 0 × ) ( r ¯ 0 × ) T + I S = I S O + I S = I S O + I S C + I S ,
In the Expression (61), I S O represents the inertia matrix containing axial and centrifugal inertia moments, which is determined with respect to the fixed system { 0 } .

3. Higher Order Acceleration Energies

The term, advanced notions, from analytical dynamics presented in this paper refers to the motion energies whose central functions are the higher order accelerations. They have been developing in any sudden and transitory motion of mechanical systems. By having Appell’s function as the starting point [1,2], a function which was highlighted in 1899 and that is also known as the kinetic energy of accelerations, new mathematical formulations on the expressions for acceleration energies of the first, second, and third order have been developed [10,14]. In this section, they are presented only in an explicit form. From a physical point of view, the acceleration energy is expressed as J / s 2 and represents the mechanical power developed by a mechanical system per time unit. By applying certain constraints, the acceleration energy can also be defined as the second order time derivative applied on the kinetic energy corresponding to the whole mechanical system. In the following, based on [10,14], the acceleration energies of the higher order are defined. They characterize the complex mechanical systems and include among others, the robots. The author established the acceleration energy corresponding to a rigid body involved in general motion, in a generalized form [10,14], also called the first order acceleration energy and defined with:
E A ( 1 ) = 1 2 a M 2 d m = 1 2 v ¯ ˙ M T v ¯ ˙ M d m = 1 2 T r a c e [ v ¯ ˙ M v ¯ ˙ M T ] d m ,      a ¯ M = a ¯ 0 + ε × ρ ¯ M + ω ¯ × ω ¯ × ρ ¯ M
By developing Expression (71), the first three terms are obtained:
1 2 a ¯ 0 T a ¯ 0 d m = 1 2 M a ¯ 0 T a ¯ 0 = 1 2 M a 0 2 ,  
1 2 a ¯ 0 T ( ε ¯ × ρ ¯ M ) d m = 1 2 ( ε ¯ × ρ ¯ M ) T a ¯ 0 d m = 1 2 M a ¯ 0 ( ε ¯ × ρ ¯ C ) ,
1 2 a ¯ 0 T ( ω ¯ × ω ¯ × ρ ¯ M ) d m = 1 2 ( ω ¯ × ω ¯ × ρ ¯ M ) T a ¯ 0 d m = 1 2 M a ¯ 0 ( ω ¯ × ω ¯ × ρ ¯ C ) .
The following three components define the resultant rotation. Among these, the first two contain the angular velocity and acceleration. The defining expressions are shown below:
1 2 ( ε ¯ × ρ ¯ M ) T ( ε ¯ × ρ ¯ M ) d m = 1 2 ε ¯ T { ( ρ ¯ M × ) ( ρ ¯ M × ) T d m } ε ¯ = 1 2 ε ¯ T I S ε ¯
Performing some transformations and developing the terms from (76), the result is:
E A ( 1 ω ε ) = 1 2 ( ε ¯ × ρ ¯ M ) T ( ω ¯ × ω ¯ × ρ ¯ M ) d m = 1 2 ε ¯ T ( ω ¯ × ) [ ( ρ ¯ M × ) ( ρ ¯ M × ) T d m ] ω ¯ = 1 2 ε ¯ T ( ω ¯ × I S ω ¯ ) .
The last component, corresponding to rotation motion, exclusively contains the angular velocity:
E A ( 1 ω 4 ) = 1 2 ω ¯ T [ ω ¯ T I S ω ¯ ] ω ¯ = 1 2 ω ¯ T [ ω ¯ T T r a c e ( I p S ) ω ¯ ω ¯ T I p S ω ¯ ] ω ¯ .  
The acceleration energy of the first order can be written in the following final form:
E A ( 1 ) = 1 2 M a ¯ 0 T a ¯ 0 + M a ¯ 0 T ( ε ¯ × ρ ¯ C ) + M a ¯ 0 T ( ω ¯ × ω ¯ × ρ ¯ C ) + 1 2 ε ¯ T I S ε ¯ + + ε ¯ T ( ω ¯ × I S ω ¯ ) + 1 2 ω ¯ T [ ω ¯ T I S ω ¯ ] ω ¯ .
If the following conditions are met: O C ,      ρ ¯ C = 0 ,    a n d      I S I S , then the mobile reference frame is applied in the mass center of the rigid body and the position of the elementary particle ( d m ) with respect to this frame is known, as well as the inertial tensor axial and centrifugal I S I S with respect to the mass center. Thus, (69) defines the first order acceleration energy in a generalized form:
E A ( 1 ) = 1 2 M a ¯ C T a ¯ C + 1 2 ε ¯ T I S ε ¯ + ε ¯ T ( ω ¯ × I S ω ¯ ) + 1 2 ω ¯ T [ ω ¯ T I S ω ¯ ] ω ¯ ,  
The expression (70) represents a reinterpretation of König’s theorem for the acceleration energy of the first order. In the case of multibody systems (which also includes robots), the expression of acceleration energy of the first order (69) is determined according to:
E A ( 1 ) [ θ ¯ ( t ) ; θ ¯ ˙ ( t ) ; θ ¯ ¨ ( t ) ] = = ( 1 ) Δ M 1 Δ M 1 + 3 Δ M i = 1 n { { 1 2 M i v ¯ ˙ ( i ) C i T v ¯ ˙ ( i ) C i } + Δ M 2 { 1 2 ω ¯ ˙ ( i ) T i [ I ( i ) i ω ¯ ˙ ( i ) i + ( ω ¯ ( i ) i × I ( i ) i ω ¯ ( i ) i ) ] } } + + Δ M 2 i = 1 n { 1 2 ω ¯ ˙ ( i ) i T ( ω ¯ ( i ) i × I ( i ) i ω ¯ ( i ) i ) } + Δ M 2 i = 1 n { 1 2 ω ¯ ( i ) i T [ ω ¯ ( i ) i T I ( i ) i ω ¯ ( i ) i ] ω ¯ ( i ) i } where   Δ M = { ( 1 , General motion ) ;    ( 0 , Translation motion ) ;    ( 1 , Rotation motion ) } , and   i = 1 n   is   the   number   of   kinetic   assemblies   that   compose   the   mechanical   system .
Considering the advanced kinematics notions, the components of acceleration energy of the first order are:
E A ( 1 ) T R = 1 2 i = 1 n M i a C i 2 = = 1 2 i = 1 n M i { j = 1 n p = 1 n [ 2 r ¯ C i ( m ) q j ( m ) q p ( m ) q ¨ j q ¨ p + 1 ( m + 1 ) 2 2 r ¯ C i ( m + 1 ) q j ( m ) q p ( m ) q ˙ j q ˙ p + 1 m + 1 r ¯ C i ( m ) q j ( m ) r ¯ C i ( m + 1 ) q p ( m ) q ¨ j q ˙ p ] }
E A ( 1 ) R O T ε = 1 2 i = 1 n ε ¯ i T I i ε ¯ i = { 1 2 i = 1 n j = 1 n [ ψ ¯ i ( m ) q j ( m ) Δ j q ¨ j + 1 m + 1 ψ ¯ i ( m + 1 ) q j ( m ) Δ j q ˙ j ] T I i ε ¯ i } ,
I i ε ¯ i = I i j = 1 n [ ψ ¯ i ( m ) q j ( m ) Δ j q ¨ j + 1 m + 1 ψ ¯ i ( m + 1 ) q j ( m ) Δ j q ˙ j ] ,  
E A ( 1 ) R O T ω ε = i = 1 n ε ¯ i T ( ω ¯ i × I i ω ¯ i ) = i = 1 n j = 1 n [ ψ ¯ i ( m ) q j ( m ) Δ j q ¨ j + 1 m + 1 ψ ¯ i ( m + 1 ) q j ( m ) Δ j q ˙ j ] T E A ( 1 ) R O T ω ω ,
E A ( 1 ) R O T ω ω = ( ω ¯ i × I i ω ¯ i ) = i = 1 n j = 1 n p = 1 n [ ψ ¯ i ( m ) q j ( m ) ] × [ I i ψ ¯ i ( m ) q p ( m ) ] Δ j Δ p q ˙ j q ˙ p .
According research [10,14], the sudden motion of MBS (multibody system), the transient motion phases, as well as the mechanical systems subjected to the action of a system of external forces, with a time variation law [19,20], are characterized by linear and angular accelerations of the higher order. Thus, the acceleration energy of the second order was also developed. According to the same research, the explicit form for the acceleration energy of the second order is presented below:
E A ( 2 ) [ θ ¯ ( t ) ; θ ¯ ˙ ( t ) ; θ ¯ ¨ ( t ) ; θ ¯ ( t ) ] = = ( 1 ) Δ M 1 Δ M 1 + 3 Δ M i = 1 n { 1 2 M i v ¯ ¨ i C i T v ¯ ¨ i C i } + E A ( 2 ) [ θ ¯ ( t ) ; θ ¯ ˙ ( t ) ; θ ¯ ¨ 2 ( t ) ] + E A ( 2 ) [ θ ¯ ( t ) ; θ ¯ ˙ 6 ( t ) ; θ ¯ ¨ ( t ) ] + + Δ M 2 i = 1 n { 1 2 ω ¯ ¨ i i T I i i ω ¯ ¨ i i + 2 ω ¯ i i T ( ω ¯ ¨ i i × I i p i ω ¯ ˙ i i ) + ω ¯ ˙ i i T ( ω ¯ ¨ i i × I i p i ω ¯ i i ) ω ¯ i i T ( ω ¯ ¨ i i T I i i ω ¯ i i ) ω ¯ i i } ,
w h e r e          E A ( 2 ) [ θ ¯ ( t ) ; θ ¯ ˙ ( t ) ; θ ¯ ¨ 2 ( t ) ] = = Δ M 2 i = 1 n { 2 ω ¯ ˙ i i T ( ω ¯ i i T I i i ω ¯ ˙ i i ) ω ¯ i i + 2 ω ¯ i i T [ ω ¯ ˙ i i T I i p i ω ¯ ˙ i i ] ω ¯ i i 5 ( ω ¯ ˙ i i T I i p i ) ( ω ¯ ˙ i i T ω ¯ i i ) ω ¯ i i + + 5 2 ( ω ¯ ˙ i i T ω ¯ i i ) T r a c e ( I i p i ) ( ω ¯ ˙ i i T ω ¯ i i ) + 1 2 ω ¯ ˙ i i T [ ω ¯ i i T I i p i ω ¯ i i ] ω ¯ ˙ i i } ,
a n d      E A ( 2 ) [ θ ¯ ( t ) ; θ ¯ ˙ 6 ( t ) ; θ ¯ ¨ ( t ) ] = = Δ M 2 { i = 1 n ω ¯ i i T [ ω ¯ i i T ( ω ¯ ˙ i i × I i p i ω ¯ i i ) ] ω ¯ i i + 1 2 ω ¯ i i T [ ω ¯ i i T ( ω ¯ i i T I i i ω ¯ i i ) ω ¯ i i ] ω ¯ i i }
The study of advanced dynamics is extended to the acceleration energy of the third order. According to [10,14], an explicit equation for the third order acceleration energy is developed:
E A ( 3 ) [ θ ¯ ( t ) ; θ ¯ ˙ ( t ) ; θ ¯ ¨ ( t ) ; θ ¯ ( t ) ; θ ¯ ( t ) ] = ( 1 ) Δ M 1 Δ M 1 + 3 Δ M i = 1 n { 1 2 M i v ¯ i C i T v ¯ i C i } + + Δ M 2 i = 1 n { 1 2 ω ¯ i i T I i i ω ¯ i i + 3 ω ¯ i T ( ω ¯ i × I p i ω ¯ ¨ i ) + 3 ω ¯ ˙ i T ( ω ¯ i × I p i ω ¯ ˙ i ) + 3 ω ¯ ¨ i T ( ω ¯ i × I p i ω ¯ i ) + + 2 ( ω ¯ i × ω ¯ i ) T I p i ( ω ¯ ˙ i × ω ¯ i ) 5 ω ¯ ˙ i T [ ω ¯ i T I i ω ¯ i ] ω ¯ i ω ¯ i T [ ω ¯ ˙ i T I i ω ¯ i ] ω ¯ i + + ω ¯ i T [ ω ¯ i T I p i ( ω ¯ i × ω ¯ i ) ] ω ¯ i }
An important aspect that has to be mentioned is that the component E A ( 3 ) [ θ ¯ ( t ) ; θ ¯ ˙ ( t ) ; θ ¯ ¨ ( t ) ; θ ¯ 2 ( t ) ] of the acceleration energy of the third order is not included in the definition Equation (80). In the case of the multibody systems, the acceleration energies of the higher order are determined based on the following generalized expression:
E A ( p ) [ θ ¯ ( t ) ; θ ¯ ˙ ( t ) ; ; θ ¯ ( t ) ( p + 1 ) ] = 1 2 i = 1 n { v ¯ i T ( p ) v ¯ i ( p ) d m } = 1 2 i = 1 n { T r a c e [ r ¯ i ( p + 1 ) r ¯ i T ( p + 1 ) ] d m } = = 1 2 i = 1 n T r a c e { [ R ] i 0 ( p + 1 ) [ r ¯ i i r ¯ i i T d m + r ¯ i C i r ¯ i C i T d m ] [ R ] i 0 T ( p + 1 ) } + 1 2 i = 1 n T r a c e [ p ¯ i ( p + 1 ) p ¯ i T ( p + 1 ) ] d m = = 1 2 i = 1 n T r a c e { [ R ] i 0 ( p + 1 ) [ I i p i + M i r ¯ i C i r ¯ i C i T ] [ R ] i 0 T ( p + 1 ) } + 1 2 i = 1 n T r a c e [ p ¯ i ( p + 1 ) p ¯ i T ( p + 1 ) ] M i ,      w h e r e      p = 1 , 2 , 3 , 4 ,
w h e r e p 1 , k 1 , { p ; k } = { 1 ; 2 ; 3 ; } , ( p ) o r d e r o f a c c e l e r a t i o n , ( k ) d e r i v a t i o n o r d e r a n d    E A ( p ) ( 0 ) [ θ ¯ ( t ) ; θ ¯ ˙ ( t ) ; ; θ ¯ ( t ) ( p + k ) ] = E A ( p ) [ θ ¯ ( t ) ; θ ¯ ˙ ( t ) ; ; θ ¯ ( t ) ( p + 1 ) ]
Expression (81) includes the inertia tensor planar and centrifugal, I i p i = r ¯ i i r ¯ i i T d m .
The expressions for the acceleration energies of the higher order, according to [13] are further included in the equations of advanced dynamics corresponding to the mechanical systems which are subjected to fast and sudden movements.

4. Discussion on the Results

This section was introduced exclusively for highlighting the existence of the accelerations and acceleration energies of the higher order. For this purpose, according to Figure 4, the rotation motion of the arm of the serial robot Fanuc LR Mate 100 iB, on the angular interval ( 0 , π ) was considered. The total time of the rotation motion was of 1.7 s.
For highlighting the time variation law of the acceleration energies of the higher order, the polynomial interpolating functions of the fifth order have been applied in accordance with (48)–(51). This method requires, according to kinematical constraints (51), the establishment through measurements of the angular acceleration corresponding to robot arm, q ¨ 3 . As a result, in the first step, the time variation law corresponding to the tangential component of the acceleration for the characteristic point belonging to the robot arm has been determined at every 0.033 s. The experimental data was collected by using a mono-axial accelerometer. The robot arm (more precisely the third kinetic ensemble of the robot) was engaged in rotation motion at high speeds and by using the data collected by the accelerometer. The time variation law for the tangential component of the acceleration of the characteristic point was graphically represented. Considering the rotation motion, in the next step the time variation law of the angular acceleration q ¨ 3 ( τ ) was also represented. Both graphical representations are illustrated in Figure 5.
According to Figure 5, the maximum value for the linear acceleration exceeds 15 m / s 2 , which defines a sudden motion. The authors have considered, for this application, that the motion can be described by using the polynomial interpolating functions of the fifth order:
q j i ( τ )    =    τ i τ t i q j i ( τ i 1 ) + τ τ i 1 t i q j i ( τ i ) ,
q j i ( τ ) =    ( τ i τ ) 2 2 t i q j i 1 + ( τ τ i 1 ) 2 2 t i q j i +    a j i 1 ,
q ¨ j i ( τ )    =    ( τ i τ ) 3 6 t i q j i 1 + ( τ τ i 1 ) 3 6 t i q j i +    a j i 1 τ +    a j i 2 ,
q ˙ j i ( τ )    =    ( τ i τ ) 4 24 t i q j i 1 + ( τ τ i 1 ) 4 24 t i q j i + a j i 1 τ 2 2 +    a j i 2 τ + a j i 3 ,
q j i ( τ ) = ( τ i τ ) 5 120 t i q j i 1 + ( τ τ i 1 ) 5 120 t i q j i + a j i 1 τ 3 6 + a j i 2 τ 2 2 + a j i 3 τ + a j i 4 .
where a j i p , p = 1 4 , are the integration constants which are determined from the geometrical and kinematical constraints (50), with an important role in ensuring the continuity of the rotation motion on the angular interval ( 0 , π ) , characterized by 51 interpolation segments. Every interpolation segment contains four unknowns, represented by the integration constants. It is obvious that the total number of the unknowns is 204. In order to determine the unknowns, according to (50), the numerical values for the angular acceleration (Figure 5b) were substituted in the functions (83)–(87) presented above. Considering the large number of results, more exactly 255 polynomial functions, only a few sequences are illustrated in Table 1 and Table 2.
As a result, the numerical values for the angular acceleration (see Figure 5) were applied, finally obtaining all of the integration constants. In order to illustrate this step, some sequences containing the time variation laws of the functions described by Equations (83)–(87) are presented. These results are included in Table 1 and Table 2. Similarly, the polynomial functions are obtained for all the 51 interpolation segments. As a result, the time variation law for the angular acceleration of the second and third order are also graphically represented according to Figure 7:
Thus, by customizing Expressions (71), (77)–(80) for the considered robot structure, the following expressions characterizing the acceleration energies of the first, second and third order are:
E A i k ( 1 ) ( τ ) = 1 2 ( M 3 x C 3 2 + M 3 z C 3 2 + I 3 y ) ( q ¨ 3 i k 2 ( τ ) + q ˙ 3 i k 4 ( τ ) ) ,
E A i k ( 2 ) ( τ ) = 1 2 ( M 3 x C 3 2 + M 3 z C 3 2 + I 3 y ) [ q 3 i k 2 ( τ ) 2 q ˙ 3 i k 3 ( τ ) q 3 i k ( τ ) + 9 q ˙ 3 i k 2 ( τ ) q ¨ 3 i k 2 ( τ ) + q ˙ 3 i k 6 ( τ ) ]
E A i k ( 3 ) ( τ ) = 1 2 { ( M 3 x C 3 2 + M 3 z C 3 2 + I 3 y ) [ q ˙ 3 i k 8 ( τ ) 8 q ˙ 3 i k 5 ( τ ) q 3 i k ( τ ) + 30 q ˙ 3 i k 4 ( τ ) q ¨ 3 i k 2 ( τ ) 12 q ˙ 3 i k 2 ( τ ) q ¨ 3 i k ( τ ) q 3 i k ( τ ) + 30 q ˙ 3 i k 4 ( τ ) q ¨ 3 i k 2 ( τ ) 12 q ˙ 3 i k 2 ( τ ) q ¨ 3 i k ( τ ) q 3 i k ( τ ) + + 16 q ˙ 3 i k 2 ( τ ) q 3 i k 2 ( τ ) + 24 q ˙ 3 i k ( τ ) q ¨ 3 i k 2 ( τ ) q 3 i k ( τ ) + 9 q ¨ 3 i k 4 ( τ ) + q 3 i k 2 ( τ ) ] }
The obtained results (see example from Table 1 and Table 2) regarding the higher order polynomial functions have been included in the above expressions of the acceleration energies that characterize the fast rotation motion of the robot arm. The results are represented in Figure 8 and Figure 9:
In the above graphical representations (Figure 8 and Figure 9), the time variation laws for the kinetic energy are illustrated, as well as the acceleration energies of the first, second and third order for an application in which the serial Fanuc robot, presented in the Figure 4, is implemented.

5. Conclusions

The main objective of this paper was to highlight some new developments in the field of advanced dynamics of mechanical systems and higher order dynamic equations. These equations were developed by using matrix exponentials which proved to have indisputable advantages in the matrix study of any complex mechanical system. In this paper, new approaches have been proposed based on differential principles from analytical mechanics by using some important dynamics notions, regarding the acceleration energies of the first, second and third order. The study was extended to the equations of the higher order, which provided the possibility of applying the initial motion conditions in the positions, velocities and accelerations of the first and second order. In order to determine the time variation laws for the generalized variables and acceleration energies of the higher order, the time polynomial functions of the fifth order were applied. In Newtonian dynamics, the establishment of dynamic equations as second order differential equations are performed by considering the kinetic energy as the main function. In the case of fast and sudden movements, the presence of higher order accelerations leads to higher order differential equations which define the variation laws in relation to time in the case of the forces and the driving moments that impose the motion of the mechanical systems. As a result, the application of acceleration energies highlights, in an explicit and direct form, the presence of higher-order accelerations, unlike the case where the study was performed considering kinetic energy.
In accordance to inverse kinematics, also known as the control kinematics of the robots, the use of polynomial functions leads to the control of the mechanical motion, especially of transitory motions. These aspects were exemplified in the previous section. Therefore, it was noticed that on the time interval of 1.7 s, a number of 51 measurements of the tangential acceleration were performed at every 0.033 s. Regarding the polynomial functions, this aspect led to a very good approximation of the mechanical motion. At the same time, the polynomial functions were substituted in the expressions of the acceleration energies of the higher order. Therefore, their time variation laws were obtained as shown in graphical representations from the previous section.

6. Contributions of the Authors

The purpose of this paper was to present some original defining expressions for acceleration energies of the higher order. These expressions also include the accelerations of the higher order which define the fast and sudden motions of mechanical systems. The novelty of this paper was in determining in an explicit form of the expressions for the acceleration energies based on advanced kinematic notions and using polynomial interpolation functions of the higher order.
Therefore, the contributions of the main author are highlighted by a great number of expressions included in the following sections of this paper:
  • 2.1. Position and Orientation Parameters. The new expressions are: (10)–(16), (19) and (20), as well as (24), (27)–(30).
  • 2.2. The Kinematic Parameters of Higher Order. The new equations are: (34)–(38) and (42)–(50);
  • 3. Higher Order Acceleration Energies. This section contains the new expressions of definition for the acceleration energy of the first, second and third order in an explicit form: (69)–(82).

Author Contributions

Conceptualization, I.N.; methodology, I.N.; investigation and software A.-V.C; validation, I.N.; formal analysis, I.N; resources, I.N., A.-V.C.; data curation, I.N., A.-V.C.; writing—original draft preparation, A.-V.C; writing—review and editing, I.N., A.-V.C.; visualization, A.-V.C; supervision, I.N.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank Technical University of Cluj Napoca, Romania for providing the technical support in carrying out the experimental analysis.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Appell, P. Sur Une Forme Générale des Equations de la Dynamique, 1st ed.; Gauthier-Villars: Paris, France, 1899. [Google Scholar]
  2. Appell, P. Traité de Mécanique Rationnelle, 1st ed.; Garnier frères: Paris, France, 1903. [Google Scholar]
  3. Craig, J.J. Introduction to Robotics: Mechanics and Control, 3rd ed.; Pearson Prentice Hall: Upper Sadelliffe, NJ, USA, 2005. [Google Scholar]
  4. Yuan, B.S.; Book, W.J.; Huggins, J.D. Dynamics of Flexible Manipulator Arms: Alternative derivation, Verification, and Characteristics for Control. ASME J. Dyn. Syst. Meas. Control 1993, 115, 394–404. [Google Scholar] [CrossRef]
  5. Bhatti, M.; Lu, D. Analytical Study of the Head-On Collision Process between Hydroelastic Solitary Waves in the Presence of a Uniform Current. Symmetry 2019, 11, 333. [Google Scholar] [CrossRef]
  6. Hussain, F.; Ellahi, R.; Zeeshan, A. Mathematical Models of Electro-Magnetohydrodynamic Multiphase Flows Synthesis with Nano-Sized Hafnium Particles. Appl. Sci. 2018, 8, 275. [Google Scholar] [CrossRef]
  7. Marin, M.; Vlase, S.; Ellahi, R.; Bhatti, M. On the Partition of Energies for the Backward in Time Problem of Thermoelastic Materials with a Dipolar Structure. Symmetry 2019, 11, 863. [Google Scholar] [CrossRef]
  8. Thompson, P. Snap, Crackle, and Pop, (AIAA Info.); Systems Technology: Hawthorne, CA, USA, 2011. [Google Scholar]
  9. Eager, D.; Pendrill, A.M.; Reinstad, N. Beyond velocity and acceleration: jerk, snap and higher derivatives. Eur. J. Phys. 2016, 37, 065008. [Google Scholar] [CrossRef]
  10. Negrean, I.; Negrean, D.C. The Acceleration Energy to Robot Dynamics. In Proceedings of the A&QT-R International Conference on Automation, Quality and Testing, Robotics, Cluj-Napoca, Romania, 23–25 May 2002; pp. 59–64. [Google Scholar]
  11. Negrean, I.; Negrean, D.C. Matrix Exponentials to Robot Kinematics. In Proceedings of the 17th International Conference on CAD/CAM, Robotics and Factories of the Future, CARS&FOF, Durban, South Africa, 10–12 July 2001; pp. 1250–1257. [Google Scholar]
  12. Negrean, I. New Formulations on Motion Equations in Analytical Dynamics. Appl. Mech. Mater. 2016, 823, 49–54. [Google Scholar] [CrossRef]
  13. Negrean, I. Advanced Notions in Analytical Dynamics of Systems. Acta Tech. Napoc. Ser. Appl. Math. Mech. Eng. 2017, 60, 491–502. [Google Scholar]
  14. Negrean, I. Energies of Acceleration in Advanced Robotics Dynamics. Appl. Mech. Mater. 2014, 762, 67–73. [Google Scholar] [CrossRef]
  15. Vlase, S.; Teodorescu, P.P. ElastoDynamics of a Solid with A General “Rigid” Motion Using Fem Model Part I. Theoretical Approach. Rom. J. Phys. 2013, 58, 872–881. [Google Scholar]
  16. Pars, L.A. A Treatise on Analytical Dynamics; Heinemann: London, UK, 2007; Volume 1, pp. 1–122. [Google Scholar]
  17. Jazar, R.N. Theory of Applied Robotics: Kinematics, Dynamics, and Control, 2nd ed.; Springer Nature Switzerland AG: Basel, Switzerland, 2010. [Google Scholar]
  18. Ardema, M.D. Analytical Dynamics. Theory and Applications; Springer Nature Switzerland AG: Basel, Switzerland, 2005. [Google Scholar]
  19. Park, F.C. Computational Aspects of the Product-of-Exponentials Formula for Robot Kinematics. IEEE Trans. Autom. Control 1994, 39, 643–647. [Google Scholar] [CrossRef]
  20. Brener, P. Technical Concepts. Orientation, Rotation, Velocity and Acceleration, and the SRM. User’s manual, version 2.0; SEDRIS. 2008. [Google Scholar]
Figure 1. (a) The position and orientation for a point; (b) The position and orientation of a system.
Figure 1. (a) The position and orientation for a point; (b) The position and orientation of a system.
Symmetry 11 01077 g001
Figure 2. The representation of a free rigid body in the Cartesian frame.
Figure 2. The representation of a free rigid body in the Cartesian frame.
Symmetry 11 01077 g002
Figure 3. The representation of a kinetical element from robot mechanical structure.
Figure 3. The representation of a kinetical element from robot mechanical structure.
Symmetry 11 01077 g003
Figure 4. Fanuc LR Mate 100 iB robot. The location of the transducer on the robot arm.
Figure 4. Fanuc LR Mate 100 iB robot. The location of the transducer on the robot arm.
Symmetry 11 01077 g004
Figure 5. (a) Tangential accelerations variation law. (b) Generalized accelerations variation law. The above presented parameters were also determined on the same angular interval ( 0 , π ) , by applying the sudden stop. The graphical representation is illustrated in Figure 6.
Figure 5. (a) Tangential accelerations variation law. (b) Generalized accelerations variation law. The above presented parameters were also determined on the same angular interval ( 0 , π ) , by applying the sudden stop. The graphical representation is illustrated in Figure 6.
Symmetry 11 01077 g005
Figure 6. (a) The tangential first order accelerations variation law; (b) The generalized accelerations of first order variation law.
Figure 6. (a) The tangential first order accelerations variation law; (b) The generalized accelerations of first order variation law.
Symmetry 11 01077 g006
Figure 7. (a) The variation law of the angular accelerations of second order. (b) The variation law for the angular accelerations of third order.
Figure 7. (a) The variation law of the angular accelerations of second order. (b) The variation law for the angular accelerations of third order.
Symmetry 11 01077 g007
Figure 8. (a) The variation law for the kinetic energy. (b) The time variation law of the acceleration energy of the first order.
Figure 8. (a) The variation law for the kinetic energy. (b) The time variation law of the acceleration energy of the first order.
Symmetry 11 01077 g008
Figure 9. (a) The variation law for the acceleration energy of second order; (b) The time variation law of the acceleration energy of third order.
Figure 9. (a) The variation law for the acceleration energy of second order; (b) The time variation law of the acceleration energy of third order.
Symmetry 11 01077 g009
Table 1. The polynomial time functions for the generalized coordinates for the third rotation joint.
Table 1. The polynomial time functions for the generalized coordinates for the third rotation joint.
Time   Interval τ  ​ s Polynomial Time Functions for Generalized Coordinates
{ k = 1 17 , i = 1 51 }
q 3 i k r a d q ˙ 3 i k s 1
[ 0 , 0.0195 ] 2478.916 τ 5 2.736 τ 3 + 0.047 τ 2 12394.582 τ 4 8.21 τ 2 + 0.094 τ
[ 0.585 , 0.605 ] 3792.503 τ 5 11239.525 τ 4 + + 13301.344 τ 3 7864.32 τ 2 + + 2331.436 τ 277.457 18962.519 τ 4 44958.1 τ 3 + + 39904.034 τ 2 15728.643 τ + + 2331.436
[ 1.191 , 1.21 ] 235.591 τ 5 1429.192 τ 4 + + 3470.091 τ 3 4216.837 τ 2 + + 2565.994 τ 622.767 1177.959 τ 4 5716.771 τ 3 + + 10410.273 τ 2 8433.675 τ + + 2565.994
[ 1.698 , 1.718 ] 62.678 τ 5 + 538.206 τ 4 1849.467 τ 3 + 3177.645 τ 2 2727.52 τ + 939.274 313.393 τ 4 + 2152.825 τ 3 5548.403 τ 2 + 6355.291 τ 2727.52
Table 2. The Polynomial Time Functions for Generalized Variables in Case of the Third Rotation Joint.
Table 2. The Polynomial Time Functions for Generalized Variables in Case of the Third Rotation Joint.
Time   Interval τ  ​ s Polynomial Time Functions for Generalized Variables
{ k = 1 17 , i = 1 51 }
q ¨ 3 i k s 2 q 3 i k s 3 q 3 i k s 4
[ 0 , 0.0195 ] 49578.331 τ 3 16.421 τ + 0.094 148734.995 τ 2 16.421 297469.99 τ
[ 0.585 , 0.605 ] 75850.076 τ 3 134874.302 τ 2 + + 79808.069 τ 15728.643 227550.23 τ 2 269748.605 τ + 79808.069 455100.461 τ 269748.605
[ 1.191 , 1.21 ] 4711.839 τ 3 17150.314 τ 2 + + 20820.546 τ 8433.675 14135.518 τ 2 34300.628 τ + 20820.546 28271.037 τ 34300.628
[ 1.698 , 1.718 ] 1253.573 τ 3 + + 6458.476 τ 2 11096.807 τ + 6355.291 3760.72 τ 2 + + 12916.952 τ 11096.807 7521.4407 τ + + 12916.9529

Share and Cite

MDPI and ACS Style

Negrean, I.; Crișan, A.-V. Synthesis on the Acceleration Energies in the Advanced Mechanics of the Multibody Systems. Symmetry 2019, 11, 1077. https://0-doi-org.brum.beds.ac.uk/10.3390/sym11091077

AMA Style

Negrean I, Crișan A-V. Synthesis on the Acceleration Energies in the Advanced Mechanics of the Multibody Systems. Symmetry. 2019; 11(9):1077. https://0-doi-org.brum.beds.ac.uk/10.3390/sym11091077

Chicago/Turabian Style

Negrean, Iuliu, and Adina-Veronica Crișan. 2019. "Synthesis on the Acceleration Energies in the Advanced Mechanics of the Multibody Systems" Symmetry 11, no. 9: 1077. https://0-doi-org.brum.beds.ac.uk/10.3390/sym11091077

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