The growth of model-based control strategies for robotics platforms has led to the need for additional rigid-body-dynamics algorithms to support their operation. Toward addressing this need, this article summarizes efficient numerical methods to compute the Coriolis matrix and underlying Christoffel Symbols (of the first kind) for tree-structure rigid-body systems. The resulting algorithms can be executed purely numerically, without requiring any partial derivatives that would be required in symbolic techniques that do not scale. Properties of the presented algorithms share recursive structure in common with classical methods such as the Composite-Rigid-Body Algorithm . The algorithms presented are of the lowest possible order: O(Nd) for the Coriolis Matrix and O(Nd^2) for the Christoffel symbols, where N is the number of bodies and d is the depth of the kinematic tree. A method of order O(Nd) is also provided to compute the time derivative of the mass matrix. A numerical implementation of these algorithms in C/C++ is benchmarked showing computation times on the order of 10-20 \mus for the computation of the Coriolis matrix and 40-120 \mus for the computation of the Christoffel symbols for systems with 20 degrees of freedom. These results demonstrate feasibility for the adoption of these numerical methods within control loops that need to operate at 1kHz rates or higher, as is commonly required for model-based control applications.
This article presents methods to efficiently compute the Coriolis matrix and underlying Christoffel symbols (of the first kind) for tree-structure rigid-body systems. The algorithms can be executed purely numerically, without requiring partial derivatives as in unscalable symbolic techniques. The computations share a recursive structure in common with classical methods such as the Composite-Rigid-Body Algorithm and are of the lowest possible order: O(Nd) for the Coriolis matrix and O(Nd^2) for the Christoffel symbols, where N is the number of bodies and d is the depth of the kinematic tree. Implementation in C/C++ shows computation times on the order of 10-20 \mus for the Coriolis matrix and 40-120 \mus for the Christoffel symbols on systems with 20 degrees of freedom. The results demonstrate feasibility for the adoption of these algorithms within high-rate (> 1kHz ) loops for model-based control applications.
R
The growth of model-based control strategies for robotics platforms has led to the need for additional rigid-body-dynamics algorithms to support their operation. Toward addressing this need, this article summarizes efficient numerical methods to
This article presents methods to efficiently
R
Symbols
symbols
D
resulting
R
any partial derivatives that would be required in symbolic techniques that do not scale. Properties of the presented algorithms share
partial derivatives as in unscalable symbolic techniques. The computations share a
R
. The algorithms presented
and
R
Matrix
matrix
R
A method of order O(Nd) is also provided to compute the time derivative of the mass matrix. A numerical implementation of these algorithms
Implementation
R
is benchmarked showing
shows
D
computation of the
R
computation of the Christoffel symbols for
Christoffel symbols on
R
These
The
R
numerical methods within control loops that need to operate at
algorithms within high-rate (>
R
rates or higher, as is commonly required
) loops