Handbook of Manufacturing Engineering and Technology

DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

Manipulator Dynamics
Shaoping Bai*, Lelai Zhou and Guanglei Wu
Department of Mechanical and Manufacturing Engineering, Aalborg University, Aalborg, Denmark

Abstract
This book chapter is about fundamentals of manipulator dynamics and their applications. Two
approaches of manipulator dynamics, namely, recursive Newton-Euler approach and the Lagrange
equations, are introduced and discussed. Examples are included to demonstrate their application in
manipulator dynamics simulations and analysis. This book chapter can provide basic understanding
on manipulator dynamics, which is applicable to manipulators, including serial and parallel
manipulators.

Introduction
Manipulator dynamics concerns with the actuation forces required to generate motion of manipu-
lators. The manipulator dynamics is expressed in terms of equations of motion, which describe the
relationships between the actuation, reaction forces, and the accelerations and motion trajectories of
manipulators. In robotics, dynamics is a fundamental for the manipulator design, analysis, and
control.
Two types of problems are involved in the dynamics study, namely, the inverse and forward
dynamics. The inverse dynamics calculates the required actuating forces at joints for a prescribed
trajectories or known kinematic information (position, velocity, and acceleration). The forward
dynamics, on the other hand, determines the joints’ and end-effector accelerations with applied
actuating forces at joints.
The equations of motion can be established with several approaches. The most used approaches
include recursive Newton-Euler approach and the Lagrange equations, which can be found in
literatures (Murray et al. 1994; Angeles 1997; Featherstone and Orin 2000). In this chapter, the
dynamics equation formulation with the above two approaches is introduced. The formulations are
confined within the rigid-body dynamics. Readers with interests in flexible body manipulators are
directed to literatures (Shabana 2013; Low 1987; Dwivedy and Eberhard 2006). Two examples are
included to demonstrate their applications, one applicable to a serial manipulator and the other to
a parallel robot.

*Email: [email protected]

Page 1 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

Recursive Newton-Euler Method

The motion of a rigid body can be described by the Newton and Euler’s equations, the former for the
translation and the latter for the rotation. Basically, the equations of motion of a manipulator can be
established for each body of the manipulator on the basis of free-body diagrams. Alternatively, the
equations of motion can be formulated through the recursive Newton-Euler method for serial
manipulators, with which the free-body diagrams can be omitted. The formulation presented is
based on the method introduced by Luh et al. (1980). This formulation is well suited for inverse
dynamics problems, i.e., to find the driving torques with known kinematics information.
The recursive Newton-Euler method consists of two steps. The first step is an outward recursion,
from the base to the end link, to calculate velocity and acceleration. The second step is an inward
recursion, from the end link to the base, to calculate the forces and torques at joints. In the following
formulation, the recursive method is applied to a manipulator with n rigid links. Coordinate systems
are established by following Denavit-Hartenberg’s convention (D-H convention) (Denavit and
Hartenberg 1955).

Outward Recursion to Calculate Velocities and Accelerations

In the inverse dynamics problems, the positions, velocities, and accelerations (yi , y_ i , €
yi for a revolute
_ €
joint and di, d i , d i for a prismatic joint) of the joints are known. The outward recursion propa gates
velocities, and accelerations of each link from the base to the end link.

Revolute joints
Referring to Fig. 1, the angular velocity of a revolute joint, propagated from link i
1 to link i, is
given by
:
vi ¼ vi
1 þ yi zi (1)

or in local components
:
Ri v0i ¼ Ri
1 v0i þ y Ri z0i (2)

where Ri stands for the orientation matrix of link i and v'i is the angular velocity of link i with respect
to the ith frame (local coordinate system). Premultiplying both sides of Eq. 2 with RTi and
simplifying yield
:
v0i ¼ Rii
1 v0i þ yi z0i (3)

where Rii
1 is the rotation matrix relating link i
1 and link i. z ' i is the unit vector [0,0,1]T parallel to
the zi axis. Similarly, the propagation of linear velocity from link i
1 to link i can be derived as


v0i ¼ Rii
1 v0i
1 þ v0i
1  pi
1
i (4)

where pi
1
i denotes the position vector of the frame of link i in the (i
1)th frame.
The angular acceleration is calculated through

Page 2 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

Fig. 1 The illustration of the symbols

_ 0i
1 þ €
_ 0i ¼ Rii
1 v
v yi z0i þ Rii
1 v0i
1  y_ i z0i (5)

The linear acceleration with respect to the link-attached coordinate system is given by


0 
_ 0i
1  pi
1
v_ 0i ¼ Rii
1 v_ 0i
1 þ v i þ v 0
i
1  v i
1  p i
1
i (6)

where v_ 0i is the linear acceleration of link i in the ith frame.

The linear acceleration of the center of mass of link i is further calculated from the velocity and
acceleration of the link


v_ 0Ci ¼ v_ 0i þ v
_ 0i  p0Ci þ v0i  v0i  p0Ci (7)

where p0C i denotes the position vector of the center of mass of link i in its own frame.

Prismatic Joints
If the joints in question are prismatic, the angular velocities are identical for link i and link i
1,
which means vi ¼ vi
1. Expressing with local components yields

Ri v0i ¼ Ri
1 v0i
1 (8)

that is,

v0i ¼ Rii
1 v0i
1 (9)

where matrix Rii
1 is independent to rotations.

Likewise, the local angular acceleration of link i is

_ 0i ¼ Rii
1 v
v _ 0i
1 (10)

Furthermore, the local linear velocity and acceleration become

Page 3 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014


v0i ¼ Rii
1 v0i
1 þ v0i
1  pi
1
i þ d_ i z0i (11)



v_ 0i ¼ Rii
1 v_ 0i
1 þ v
_ 0i
1  pi
1
i þ v0i
1  v0i
1  pi
1
i þ 2v0i  d_ i z0i þ d€i z0i (12)

where di is the displacement from axis xi
1 to xi along the zi axis.

Inward Recursion to Calculate Forces and Torques

With the velocities and accelerations found, the inertial force and torque applied to link i at its center
of mass can be calculated by applying the Newton-Euler equations:

f 0C i ¼ mi v_ 0Ci (13)

_ 0i þ v0i  Ii v0i
n0Ci ¼ Ii v (14)

where mi is the mass of the link and Ii denotes the inertia tensor of link i about its center of mass,
calculated in a frame identical to the link’s coordinate system except the origin. In the equations
above, all items are expressed in the local coordinate systems.
The forces and torques acting on the links are calculated through the outward recursion. The joint
torques that generate the forces and torques will be calculated through the inward recursion.
Referring to Fig. 1, the inward recursion propagates the joint forces and torques from the end link
to the base as

f 0i ¼ f 0C i þ Riiþ1 f 0iþ1 (15)

n0i ¼ n0C i þ Riiþ1 n0iþ1 þ p0Ci  f 0C i þ piiþ1  Riiþ1 f 0iþ1 (16)

where f'i and n'i are the force and torque exerted on link i with respect to the ith frame.
The actuating torque on a revolute joint can be calculated by
0
ti ¼ niT z0i (17)

For prismatic joints, the driving force is

0
f i ¼ f iT z0i (18)

Forward Dynamics Problem

In the forward dynamics problem, the joint accelerations are to be calculated with known joint
torques. One method called articulated-body algorithm (ABA) (Featherstone 1983) is applicable to
this type of problem with recursive approach. More dynamics modeling methods were reviewed and
discussed in (Featherstone and Orin 2000).

Lagrange’s Equation
The equations of motion for a manipulator can also be formulated with Lagrangian dynamics
method. Two types of Lagrange’s equations can be formulated, one with dependent coordinates

Page 4 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

(1st kind) and the other with independent coordinates (2nd kind). This section starts with the
Lagrange’s equation of the 1st kind. Hereafter, the vector q represents a set of n unknown dependent
coordinates, m is the total number of independent kinematically constraint equations, and f ¼ n

m is the number of dynamic degrees of freedom (de Jalón and Bayo 1994). The constraint conditions
for the entire system are written in the following general form:

FðqÞ ¼ 0 (19)

Differentiating Eq. 19 with respect to time leads to

dFðqÞ @F dq
¼ ¼ Fq q_ ¼ 0 (20)
dt @q dt

where matrix Fq is the Jacobian matrix of the constraint Eq. 19.

Let L ¼ T
V be the Lagrangian of the system, where T ¼ T(q, q) _ and V ¼ V(q) are the kinetic
and potential energies, respectively, and Qex is the vector of generalized external forces acting along
the dependent coordinates q. Considering a manipulator as a constrained multibody system, the
Lagrange’s equation can be derived by generalizing the Hamilton’s principle by means of Lagrange
multiplier technique(de Jalón and Bayo 1994). Skipping details of derivation, the equations of
motion become
 
d @L @L

þ FTq l ¼ Qex (21)
dt @ q_ @q

which is the Lagrange’s equation of the 1st kind. In the equation, l is a m-dimensional vector of
Lagrange multipliers. Item FTql represents the generalized forces due to reaction forces at joints.
With independent coordinates, the constraint Jacobian vanishes, and the term with Lagrange
multipliers is dropped out. The Lagrange’s equation of the 2nd kind is thus obtained:
 
d @L @L

¼ Qex (22)
dt @ q_ @q

The Lagrange’s equation of the 1st kind stands for a system of n equations for n dependent
co-ordinates and m unknown Lagrange multipliers. To solve the equations, additional acceleration
equations from the kinematic constraints are required, as presented below.
The kinetic energy of a robotic system can be written as follows:

1
T ¼ q_ T MðqÞq_ (23)
2
where M(q) ¼ M is the mass matrix of the system. The first and second terms in Eq. 21 can be
expressed as below:
 
d @L _ q; @L
¼ M€ qþM _ ¼ Lq ¼ T q
V q (24)
dt @ q_ @q

For the general case in which the kinetic energy depends on q, Eq. 21 becomes

Page 5 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

M€ _ q_
q þ FTq l ¼ Qex þ Lq
M (25)

On the other hand, differentiating the constraint Eq. 19 twice with respect to time yields

Fq q _ q q_ ¼ 0
€þF (26)

_ q q_ and n ¼ Qex + Lq
M

Let c ¼
F _ q.
_ Combining Eqs. 25 and 26 leads to:
    
M FTq €
q n
¼ (27)
Fq 0 l c

Equation 27 can solve simultaneously the accelerations and the Lagrange multipliers.

Examples
Two examples are included to apply the equations of motion to robotic manipulators. In the first
example, a lightweight robotic arm consisting of five revolute joints developed at Aalborg Univer-
sity (AAU), Denmark, is considered. In the second example, the dynamics of a spherical parallel
manipulator is presented.

Example I: A 5-dof Lightweight Robotic Arm

The lightweight robotic arm is demonstrated in Fig. 2a. Following the D-H convention, Cartesian
coordinate systems are attached to each link of the manipulator, as shown in Fig. 2b. The D-H
parameters of the manipulator are listed in Table 1.

Dynamics Modeling
The kinematics of the robotic arm can be found in (Zhou et al. 2011). This chapter includes here only
the dynamics formulation.
Jacobian matrix The joint angular velocity can be calculated with the Jacobian matrix

u_ ¼ J
1 vef (28)

h
_ y_ 1 , y_ 2 . . . , y_ n T denotes an n-dimensional (n denotes the number of dof) vector of the joint
where u¼
angular velocities, J is the Jacobian of the robotic arm, and vef the velocity of the end-effector.
For a revolute joint, the Jacobian matrix can be calculated by (Tsai 1999)
 
zi
1  pi
1
J ¼ ½j1 ; j2 ; . . . ; jn , ji ¼ (29)
zi
1

where zi
1 and pi
1 are given by


zi
1 ¼ Rii
1 0 0 1 T , pi
1 ¼ Rii
1 qi
1 þ pi (30)

where qi
1 ¼ [ai cos yi, ai sin yi, di]T. The local coordinates of the end-effector are defined as

Page 6 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

Fig. 2 AAU 5-dof robotic arm: (a) a prototype and (b) coordinate systems

Table 1 D-H parameters of the 5-dof robotic arm

Joint i ai ai di yi
1 p/2 0 h1 y1
2 0 l1 0 y2
3 p/2 0 0 y3
4
p/2 0 l2 y4
5 p/2 0 d1 y5

pn ¼ [0, 0, 0]T. When the desired end-effector velocity vef is given, the joint angular velocity can be
solved by Eq. 28.
€ can be calculated through time derivative of the velocity u.
The joint angular acceleration u _

€ ¼ J_
1 ve f þ J
1 v_ e f
u (31)

Simulation Results
The dynamics of the 5-dof robotic arm was simulated with a program implemented in MATLAB.
The trajectory of the end-effector in the global coordinate system is defined as xef (t) ¼ 5 + 400
(1
cos(t)), yef (t) ¼
990 + 800(1
cos(t/2)), and zef (t) ¼ 280 + 250(cos(t/2)
1), all with
unit of mm, which leads to

Page 7 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

2 3
400 sin ðt Þ
x_ e f ¼ 4 400 sin ðt=2Þ 5

125 sin ðt=2Þ

The Euler angles for the end-effector are given as c ¼ [0, cos(t/180), 0]T, following ZXZ
convention. The rotation implies the end-effector remains horizontal during the prescribed motion.
For the prescribed motion, the rotation matrix R and angular velocity v can be readily found. The
h iT
velocities of the end-effector are vef ¼ x_ Tef ; vT , from which the joint angular velocities and
accelerations are calculated.
To implement the Newton-Euler method to solve the inverse dynamics problem, the equations in
section “Recursive Newton-Euler Method” are summarized in a recursive manner:
Outward recursion: i : 1 ! 5

v0i ¼ Rii
1 v0i
1 þ y_ i z0i (32a)

v0i ¼ Rii
1 v0i
1 þ €

yi z0i þ Rii
1 v0i
1  y_ i z0i (32b)



v_ 0i ¼ Rii
1 v_ ii
1 þ v
_ 0i
1  pi
1
i þ v0i
1  v0i
1  pi
1
i Þ (32c)

v_ 0Ci ¼ v_ 0i
1 þ v
_ 0i
1  p0C i þ v0i
1  v0i
1  p0Ci (32d)

f 0Ci ¼ mi v_ 0C i (32e)

_ 0i þ v0i  Ii v0i
n0Ci ¼ Ii v (32f )

Inward recursion: i : 5 ! 1

f 0i ¼ f 0C i þ Riiþ1 f 0iþ1 (32g)

n0i ¼ n0C i þ Riiþ1 n0iþ1 þ p0Ci  f 0C i þ p0iþ1  Riiþ1 f 0iþ1 (32h)

0
ti ¼ niT z0i (32i)

The joint angular velocities and accelerations are solved through Eq. 28, with results displayed in
Fig. 3. Let the payload at the end-effector of the robot be 5 kg, which implies f6 ¼ 5gN, n6 ¼ 0,
where g is the vector of gravity acceleration. Through the modeling of inverse dynamics of
Eqs. 32g–32i, the joint torques are solved, as shown in Fig. 4.

Example II: A Spherical Parallel Manipulator

In this example, the Lagrange’s equation is illustrated with a spherical parallel manipulator (SPM)
(Gosselin and Angeles 1989). The SPMs are of closed-kinematic chain. The Lagrange’s equation of
the 1st kind is adopted for the dynamics modeling of this manipulator.
A general SPM is shown in Fig. 5a with the parameterized ith leg in Fig. 5b. The ith limb consists
of three revolute joints, whose axes are parallel to the unit vectors ui, vi, and wi. All three limbs have

Page 8 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

Fig. 3 Angular velocities and accelerations of the robotic arm joints

Fig. 4 Joint torques of the robotic arm

identical architectures, defined by angles a1 and a2. Moreover, b and g define the geometry of two
triangular pyramids on the base and the mobile platforms, respectively. The origin O of the base
coordinate system xyz is located at point O. The z axis is normal to the bottom surface of the base
pyramid and points upward, while the y axis is located in the plane made by the z axis and u1.

Kinematic Modeling
Under the prescribed coordinate system, unit vector ui is derived as

ui ¼ ½
sin i sin g cos i sin g
cos g T (33)

where i ¼ 2(i
1)p/3, i ¼ 1, 2, 3.

The unit vector vi of the axis of the intermediate revolute joint in the ith leg is expressed as:

Page 9 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

Fig. 5 A general SPM and the coordinate system, where only one leg is shown for clarity
2 3

si sgca1 þ ðci syi
si cgcyi Þsa1
vi ¼ 4 ci sgca1 þ ðsi syi þ ci cgcyi Þsa1 5 (34)

cgca1 þ sgcyi sa1

where s() ¼ sin() and c() ¼ cos().

The unit vector wi of the top revolute joint in the ith leg is a function of the orientation of the
mobile platform (MP), namely,

wi ¼ wix wiy wiz T ¼ Qwi (35)

where Q is the rotation matrix that carries the MP from its reference orientation to the current one
and w*i is the unit vector of the axis of the top revolute joint in the ith leg when the mobile platform is
in its reference orientation, which is given as

wi ¼
sin i sin b cos i sin b cos b T (36)

Kinematic Jacobian Matrix

Let v ¼ [ox, oy, oz]T denote the angular velocity of the mobile platform; the velocity equation via
the ith leg can be stated as

v ¼ y_ i ui þ c_ i vi þ x_ i wi (37)

To eliminate c_ i and x_ i , dot-multiplying Eq. 37 on both sides with vi  wi yields

ðvi  wi Þ  v ¼ ðui  vi Þ  wi y_ i (38)

from which the SPM velocity equation is obtained as

Page 10 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

Av ¼ Bu_ (39)

with

A ¼ v1  w1 v2  w2 v3  w3  T (40a)

B ¼ diag ½ ðu1  v1 Þ  w1 ðu2  v2 Þ  w2 ðu3  v3 Þ  w3  (40b)

h iT
_ y_ 1 ; y_ 2 ; y_ 3 . Matrices A and B are the forward and inverse Jacobian matrices of the
where u¼
manipulator, respectively. The kinematic Jacobian matrix J of the manipulator can be expressed as
follows as long as matrix A is not singular:
vi  wi
J ¼ B
1 A ¼ ½ j1 j2 j 3 T ; j i ¼ (41)
ui  vi  wi

Equation 38 is thus rewritten for the single joint velocity as

y_ i ¼ jTi v (42)

Velocity and Acceleration Analysis

h iT
The motions of the link and mobile platform are shown in Fig. 6. The angle rates f_ ¼ f;
_ y;
_ s_ and
the angular velocity v are linearly dependent, namely, v¼C f. _ Differentiating the equation with
respect to time yields

v € þ c_ f_
_ ¼ cf (43)

Matrix C is dependent on the rotations. For example, a rotation with the Euler convention ZY Z of
Q ¼ Rz(f)Ry(y)Rz(s
f), matrix C is

Fig. 6 The movement of the mobile platform and links

Page 11 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

2 3

sycf
sf sycf
c ¼ 4
sysf cf sysf 5 (44)
1
cy 0 cf

The velocity c_ of the intermediate joint of ith leg is found by making use of Eq. 37 to eliminate y_ i
and x_ i . Dot-multiplying Eq. 37 on both sides with ui  wi leads to

ðui  wi ÞT
ðui  wi Þ  v ¼ c_ ðui  wi Þ  vi or c_ ¼ v  jTci v (45)
ðui  wi Þ  vi

The angular velocity of the distal link in the ith leg in the reference frame xyz is found as
ri ¼ y_ i ui þ c_ i vi . Let rli denote the corresponding angular velocity in the local frame xiyizi, which
means

ri ¼ ½ eix eiy eiz rli or ri ¼ Ei rli (46)

with

vi þ wi vi þ wi vi þ wi
eix ¼ , eiy ¼ , eiz ¼ (47)
jjvi þ wi jj jjvi þ wi jj jjvi þ wi jj

which yields

rli ¼ ETi ri (48)

Dynamic Modeling
The Lagrange’s equation for the SPM is expressed as
 
d @L @L

þ CTq l ¼ Qex (49)
dt @ q_ @q

where q ¼ [y1, y2, y3, f, y, s]T, Qex ¼ [tT, 0]T 2 ℝ6 is the vector of external forces, and vector
t ¼ [t1, t2, t3]T characterizes the actuator torques. Moreover, l ¼ [l1, l2, l3]T is a vector of
Lagrange multipliers. Matrix Cq is the system’s constraint Jacobian, which can be found from the
velocity Eq. 39, namely,
 
Bu_
Av ¼ B _ T T ¼ 0

Ac  u_ T f (50)

that is,

Cq q_ ¼ 0 (51)

with Cq ¼ [B
AC].

Page 12 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

Lagrangian of the Mobile Platform

The local frame xpypzp of the MP is established with the origin located at point P, i.e., the center of
mass of the MP of position vector p. Henceforth, the Lagrangian of the mobile platform is obtained
as

1 0
Lp ¼ T p
V p ¼ v T Ip v0
mp gT p (52)
2
where Ip denotes the local inertia tensor of the mobile platform and v0 its local angular velocity.
Moreover, g ¼ [0, 0, 9.81]T.

Lagrangian of a Single Leg

The Lagrangian of the ith leg is derived as

1 1
Li ¼ T i
V i ¼ I l1 y_ 2i þ rTli Il2 rli
ml1 gT hi
ml2 gT ki i ¼ 1, 2, 3 (53)
2 2
where Il1 is the proximal link’s mass moment of inertia about ui and Il2 is the distal link’s local inertia
tensor about point O. Moreover, hi and ki indicate the position vectors for the centers of the mass of
the proximal and distal links, respectively.
Substituting the Lagrangian Lp and Li, i ¼ 1, 2, 3, into Eq. 49, the inertia matrix can be derived as
 
M1 03
MðqÞ ¼ (54)
03 M2

with

M1 ¼ diag ½ I l1 I l1 I l1  (55a)
h X3

i
M2 ¼ CT QIp QT þ i¼1
j u
i i
T
þ j v
ci i
T
EI l2 ET
u j
i i
T
þ v j T
i ci C (55b)

and the other terms also can be derived. Moreover, differentiating Eq. 51 with respect to time yields

€ ¼
C_ q q_
Cq q (56)

with
h 
€ ¼ Cq u
Cq q €T € T T ;
f C_ q q_ ¼ B
_
_ AC
A _ q_
C (57)

where

_
A¼ jA1 jA2 jA3 T (58a)

_
B¼diag ½ jB1 jB2 jB3  (58b)

and

Page 13 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

jAi ¼ y_ i ðui  vi Þ  wi þ vi  ðv  wi Þ (59a)

jBi ¼ y_ i ½ui  ðui  vi Þ  wi þ ðui  vi Þ  ðv  wi Þ (59b)

Combining Eqs. 49 and 56 leads to

    
MðqÞ CTq €
q _ ðqÞq_ þ T q
V q
τ
M
¼ (60)
Cq 0 l
C_ q q_

With an external moment vector m, the actuator torques are expressed as

τ a ¼ τ
J
T m (61)

Simulation Results
The previously developed dynamic model is applied to a co-axial SPM (Bai 2010), whose properties
are given in Tables 2 and 3, where inertia tensors Ip and Il2 contain only diagonal entries. Simulations
were conducted for motion of constant angular accelerations at actuating joints, their initial
conditions being specified in Table 2. The results of the actuating torques are shown in Fig. 7.
The results agree well with the simulation results obtained from CAE software, namely, the MSC
Adams.

Summary
This book chapter introduces the fundamental of manipulator dynamics, with focus on the recursive
Newton-Euler approach and the Lagrange equations, which are approaches mostly common used in

Table 2 Parameters of the SPM and the initial simulation condition

a1 [deg] a2 [deg] b [deg] g [deg] [f, y, s] [rad] [y_ 1 , y_ 2 , y_ 3 ] [rad/s] m [Nm]
60 80 75 0 [0, p/6, 0] [
8,
7,
9] [0.1, 0.1, 0.1]

Table 3 Mass and inertia properties of the SPM

Mobile platform Proximal link Distal link
mp [kg] Ip [10
4 kg m2] ml1 [kg] Il1 [10
3 kg m2] ml2 [kg] Il2 [10
4 kg m2]
0.332 [3.855 3.855 7.688] 0.123 1.169 0.114 [1.852 0.089 1.921]

Page 14 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

Fig. 7 The simulation of inverse dynamics: actuator torques

robot dynamics analysis and control. Both inverse and forward dynamic problems are discussed.
Applications to manipulators are demonstrated.
In solving the dynamics equation, numerics problems may arise. One is computational efficiency.
The recursive approach presented in this chapter is quite efficient, while the Lagrange’s equations, as
a closed-form equation, implies high computational cost in finding solutions, which uses numerical
methods, for example, Newton–Raphson method. Other problems include the singular Jacobian and
the instability of the solution, etc., among others. The Jacobian matrix can become singular when a
singular configuration is reached or redundant constraints are included, which may lead to a crash of
simulation or a large error. Instability takes place in numerical integration of the equation of motion,
when the roundoff errors increase with time, which can lead to the constraint equations not satisfied.
A method to improve the stability is the Baumgarte stabilization method, which is introduced in
(de Jalón and Bayo 1994).
The approaches presented in this chapter are closely relevant to multibody dynamics, which was
extensively studied in the last century. With the advances of robotics in the twenty-first century,
manipulator dynamics has significantly been expanded, not limited to multibody dynamics, but get
involved with other subsets of the dynamics. A new development among others is the continuum
manipulators, where continuum mechanics is the major theory applicable (Walker 2013). The new
development of manipulator dynamics, which is not discussed in this chapter, presents new
opportunity and challenges for robotics.

References
Angeles J (1997) Fundamentals of robotic mechanical systems: theory, methods, and algorithms.
Springer, New York
Bai S (2010) Optimum design of spherical parallel manipulators for a prescribed workspace. Mech
Mach Theory 45(2):200–211
de Jalón JG, Bayo E (1994) Kinematic and dynamic simulation of multibody systems: the real-time
challenge. Springer, New York

Page 15 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

Denavit J, Hartenberg RS (1955) A kinematic notation for lower pair mechanisms based on
matrices. ASME J Appl Mech 77:215–221
Dwivedy SK, Eberhard P (2006) Dynamic analysis of flexible manipulators, a literature review.
Mech Mach Theory 41(7):749–777
Featherstone R (1983) The calculation of robot dynamics using articulated-body inertias. Int J Robot
Res 2(1):13–30
Featherstone R, Orin D (2000) Robot dynamics: equations and algorithms. In: Proceedings of the
2000 IEEE International Conference on Robotics & Automation, San Francisco, pp 826–834, Apr
2000
Gosselin CM, Angeles J (1989) The optimum kinematic design of a spherical three-degree-of-
freedom parallel manipulator. ASME J Mech Des 111(2):202–207
Low KH (1987) A systematic formulation of dynamic equations for robot manipulators with elastic
links. J Robot Syst 4(3):435–456
Luh JYS, Walker MW, Paul RPC (1980) On-line computational scheme for mechanical manipula-
tors. J Dyn Sys Meas Control 102(2):69–76
Murray RM, Li Z, Sastry SS (1994) A mathematical introduction to robotic manipulation. CRC
Press, Boca Raton
Shabana AA (2013) Dynamics of multibody systems. Cambridge University Press, New York, USA
Tsai LW (1999) Robot analysis: the mechanics of serial and parallel manipulators. Wiley, New York
Walker ID (2013) Continuous backbone continuum robot manipulators. ISRN Robot Article ID:
726506, 19pp
Zhou L, Bai S, Hansen MR (2011) Design optimization on the drive train of a light-weight robotic
arm. Mechatronics 21(3):560–569

Page 16 of 17
Handbook of Manufacturing Engineering and Technology
DOI 10.1007/978-1-4471-4976-7_91-1
# Springer-Verlag London 2014

Index Terms:
Industrial robots 6
Lagrange's equation 4
Light-weight robot 6
Newton-Euler equations 4, 8
Recursive Newton-Euler method 2
Spherical parallel manipulator 8

Page 17 of 17