August 28, 2006
We present a combination of tools which allows for investigation of the coupled orbital and rotational dynamics of two rigid bodies with nearly arbitrary shape and mass distribution, under the influence of their mutual gravitational potential. Methods for calculating that mutual potential and resulting forces and moments for a polyhedral body representation are simple and efficient. Discrete equations of motion, referred to as the Lie Group Variational Integrator (LGVI), preserve the structure of the configuration space, SE(3), as well as the geometric features represented by the total energy and the total angular momentum. The synthesis of these approaches allows us to simulate the full two body problem accurately and efficiently. Simulation results are given for two octahedral rigid bodies for comparison with other integration methods and to show the qualities of the results thus obtained. A significant improvement is seen over other integration methods while correctly capturing the interesting effects of strong orbit and attitude dynamics coupling, in multiple scenarios.
Similar papers 1
June 22, 2018
We used binary octahedrons to investigate the dynamical behaviors of binary asteroid systems. The mutual potential of the binary polyhedron method is derived from the fourth order to the sixth order. The irregular shapes, relative orbits, attitude angles, as well as the angular velocities of the binary asteroid system are included in the model. We investigated the relative trajectory of the secondary relative to the primary, the total angular momentum and total energy of the ...
July 7, 2021
A new technique that utilizes surface integrals to find the force, torque and potential energy between two non-spherical, rigid bodies is presented. The method is relatively fast, and allows us to solve the full rigid two-body problem for pairs of spheroids and ellipsoids with 12 degrees of freedom. We demonstrate the method with two dimensionless test scenarios, one where tumbling motion develops, and one where the motion of the bodies resemble spinning tops. We also test th...
August 19, 2005
We develop the equations of motion for full body models that describe the dynamics of rigid bodies, acting under their mutual gravity. The equations are derived using a variational approach where variations are defined on the Lie group of rigid body configurations. Both continuous equations of motion and variational integrators are developed in Lagrangian and Hamiltonian forms, and the reduction from the inertial frame to a relative coordinate system is also carried out. The ...
May 6, 2008
This paper formulates optimal control problems for rigid bodies in a geometric manner and it presents computational procedures based on this geometric formulation for numerically solving these optimal control problems. The dynamics of each rigid body is viewed as evolving on a configuration manifold that is a Lie group. Discrete-time dynamics of each rigid body are developed that evolve on the configuration manifold according to a discrete version of Hamilton's principle so t...
December 8, 2016
This work presents an elegant formalism to model the evolution of the full two rigid body problem. The equations of motion, given in a Cartesian coordinate system, are expressed in terms of spherical harmonics and Wigner D-matrices. The algorithm benefits from the numerous recurrence relations satisfied by these functions allowing a fast evaluation of the mutual potential. Moreover, forces and torques are straightforwardly obtained by application of ladder operators taken fro...
January 2, 2019
The two full body problem concerns the dynamics of two spatially extended rigid bodies (e.g. rocky asteroids) subject to mutual gravitational interaction. In this note we deduce the Euler-Poincare and Hamiltonian equations of motion using the geometric mechanics formalism.
October 8, 2010
This paper presents an analytical model and a geometric numerical integrator for a tethered spacecraft model that is composed of two rigid bodies connected by an elastic tether. This model includes important dynamic characteristics of tethered spacecraft in orbit, namely the nonlinear coupling between deformations of tether, rotational dynamics of rigid bodies, a reeling mechanism, and orbital dynamics. A geometric numerical integrator, referred to as a Lie group variational ...
March 23, 2021
Spin-orbit coupling of planetary systems plays an important role in the dynamics and habitability of planets. However, symplectic integrators that can accurately simulate not only how orbit affects spin but also how spin affects orbit have not been constructed for general systems. Thus, we develop symplectic Lie-group integrators to simulate systems consisting gravitationally interacting rigid bodies. A user friendly package (GRIT \url{https://github.com/GRIT-RBSim/GRIT}) is ...
November 2, 2016
We introduce a symplectic dual quaternion variational integrator(DQVI) for simulating single rigid body motion in all six degrees of freedom. Dual quaternion is used to represent rigid body kinematics and one-step Lie group variational integrator is used to conserve the geometric structure, energy and momentum of the system during the simulation. The combination of these two becomes the first Lie group variational integrator for rigid body simulation without decoupling transl...
January 10, 2005
An approach is developed to find approximate solutions to the classical Newtonian problem of N bodies. Sets of N gravitating bodies having spherically symmetric mass distributions, small angular velocities (< 1 rad/s) and bounded position vectors have been taken into consideration. In addition, it is assumed that the masses form an isolated system in free space and perform free gravitating motion. Although the problem is not exactly solvable, a new approach will be developed ...