The Journal of the Astronautical Sciences, vol. 40, pp. 27-50, Jan-March 1992
Spatial Op erator Algebra for Multib o dy
1
System Dynamics
2 2 3
G. Ro driguez , A. Jain and K. Kreutz{Delgado
Abstract
This pap er describ es a new spatial operator algebra for the dynamics of general{top ology
rigid multib o dy systems. Spatial op erators allow a concise and systematic formulation of the
dynamical equations of motion of multib o dy systems and the development of ecient computational
algorithms. Equations of motion are develop ed for progressively more complex systems: serial
chains, top ological trees, and closed{lo op systems. New op erator factorizations and expressions
for the mass matrix and its inverse are derived and used to obtain ecient, spatially recursive
computational algorithms. The algorithms can be easily recon gured in resp onse to changes in
the constraints and the top ology of constituent b o dies. Thus, they are particularly suited for
time{varying multib o dy systems. References are provided for extensions to exible multib o dy
systems. Spatially recursive algorithms, based on the sequential ltering and smo othing metho ds
encountered in optimal estimation theory, provide the computational infrastructure to mechanize
the spatial op erators.
Intro duction
There are two ma jor current challenges in multib o dy system dynamics: achieving com-
putational eciency for increasingly complex systems; and retaining algorithmic eciency while
accounting for p ossible event-dep endentchanges to the constraints and top ology of the constituent
b o dies. These challenges o ccur for example in highly complex and interactive spacecraft and in
space rob otic systems which require that the dynamics algorithms be quickly recon gured in re-
sp onse to con guration changes without sacri cing computational eciency.
In recent years, there has b een signi cant progress in the formulation of the equations of
motion and on the development of ecient dynamics algorithms for multib o dy systems. A large
1 rd
Presented at the 3 Annual Conference on Aerospace Computational Control, August 28{30, 1989, Oxnard, CA
2
Jet Propulsion Lab oratory/California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109
3
AMES Department, University of California, San Diego, La Jolla, CA 92093 1
part of this work, [1 ]{[6 ], has fo cused on the development of recursive computational algorithms
for the dynamics of serial rigid multib o dy systems. References [5 , 7] describ e recursive forward
dynamics algorithms for tree top ology rigid multib o dy systems. References [8 , 9] contain recursive
algorithms sp eci c to closed top ology rigid multib o dy systems consisting of multiple rob ot manipu-
lators grasping a rigid ob ject. References [10 , 11 ] describ e recursive algorithms for general top ology
rigid multib o dy systems.
The spatial operator algebra [8, 9] arose from a recognition of the close parallels b etween
the structure of the equations of motion for serial chain dynamics and those encountered in the area
of optimal estimation theory [12 ]. These parallels naturally led to the intro duction of spatial op er-
ators to obtain a concise formulation of the equations of motion of multib o dy systems. The spatial
op erators have b een used to obtain new op erator factorizations and expressions for the mass matrix
and its inverse. These op erator expressions make it p ossible to recognize high-level mathematical
patterns asso ciated with the mass matrix which the detailed algorithms do not reveal. Therefore,
the number of symb ols that the analyst has to see is reduced signi cantly. In addition, the spa-
tial op erators are mechanized byvery ecient spatially recursive computational algorithms whose
complexity dep ends only linearly on the number of degrees of freedom dofs. These algorithms
closely resemble the algorithms used for recursive ltering and smo othing in the Kalman ltering
and estimation problems. There is the reassuring presence of such familiar concepts as Riccati
equations and Kalman gains. This makes it easier to mechanize the dynamics algorithms and to
monitor their numerical stability and robustness. Much of the exp erience gained over manyyears
of research in estimation theory can now b e used to solvemultib o dy system dynamics problems.
In this pap er, the spatial op erator algebra is advanced as a new systematic pro cedure for
concisely formulating the equations of motion and deriving ecient spatially recursive dynamics
algorithms for general top ology rigid multib o dy systems. Spatial op erator formulations of the
dynamics of serial chain, tree and closed chain top ology multib o dy systems are develop ed in a
progressive sequence. It is seen that the spatial op erator formulation of the dynamics for serial
and tree top ology systems are identical in form. Consequently, such results as the op erator factor-
ization and inversion of the mass matrix for serial chains are directly applicable to tree top ology
systems. Moreover, the recursive computational algorithms for serial chains extend with only minor
mo di cations to tree top ology systems.
Closed top ology systems are mo deled as consisting of primary and secondary tree-top ology
subsystems along with additional kinematical closure constraints within and between them. The
primary system represents the ma jor part of the system whose top ology and constituent b o dies do
not change with time. The secondary system is much smaller in most applications, and represents 2
the part of the system whose top ology and constituent b o dies do change with time.
A recursive computational algorithm is then develop ed for closed{chain systems. This algo-
rithm contains separate steps for the dynamics of the primary and secondary systems using recursive
algorithms for tree{top ology systems, plus additional steps to handle the closure constraints. A
signi cant feature of these algorithms is that only relatively small and lo calized changes are needed
for recon guration in resp onse to changes in either the top ology, the constraints or the constituent
b o dies of the system. The overall computational complexity of the algorithm stays linear in the
number of degrees of freedom in the system. In addition, there is a linear dep endency in the
absence of kinematical singularities on the numb er of closed lo ops in the system.
Extensions of the spatial op erator algebra and the computational algorithms to exible
multib o dy system dynamics are discussed in reference [13 ].
6-Dimensional Spatial Notation
Co ordinate{free 6-dimensional spatial notation used throughout this pap er. Given the linear
3
and angular velo cities v and ! , the linear force F , and moment N in R at a p oint on a body, the
6
inertial spatial velocity V , inertial spatial acceleration and the spatial force f in R are:
0 1 0 1
! N
4 4 4
_
@ A @ A
V = ; = V; f =
v F
Here \ _ " implies the time derivative in an inertial frame. The rigid body transformation operator
66
x; y 2 R for two p oints x and y is:
0 1
~
I l
4 x;y
@ A
l =x; y =
x;y
0 I
3
~
where l 2 R is the vector joining the two p oints x and y . l is the cross{pro duct matrix
x;y x;y
asso ciated with l which acts on a vector to pro duce the cross{pro duct of l with the vector.
x:y x;y
x; y and x; y transform spatial forces and spatial velo cities resp ectively between the two
66
b o dy{ xed p oints x and y on a rigid b o dy. The spatial inertia M 2 R of a rigid b o dy at at a
O
p oint O on the b o dy is
0 1
J O mp~
4
@ A
M O = =