Multibody Simulation of REMORA CubeSat Docking to and Pushing a Spent Rocket Booster

This paper details the multibody simulation of three phases of critical importance to the feasibility of the REMORA CubeSat space debris mitigation concept: the final approach of the CubeSat to a spent rocket booster; the grappling of the spent rocket booster using a robotic arm; and the pushing of the spent rocket booster to divert its course from another on-orbit asset. The extension of a robotic mobility and manipulation modeling toolkit (M3TK) from multibody dynamics simulation of manipulators and ground vehicles to simulation of orbital robotics is outlined. This includes the identification of the appropriate parameters required to concisely and generically describe thruster loads, thruster mixing, spacecraft control, and spacecraft navigation for the purpose of on-orbit robotics simulation. A high-level spacecraft navigator commands maneuvers to target spacecraft states. A PID spacecraft controller takes the target states and calculates desired forces and torques. A thruster mixer solves a quadratic program to determine the optimal thruster firing times for the propulsion system. Pulse width modulated actuation of eight canted cold gas thrusters is used in the simulated approach to a rocket nozzle from a distance of 200 m. A five degree of freedom robotic arm is controlled to position a pair of pincers to grasp the rocket nozzle. Contact dynamics are used to accurately simulate the grasping of the rocket nozzle by the pincers. A similar, but separate simulation is performed to assess the ability of the REMORA CubeSat to push the large spent rocket booster. This diversion maneuver makes use of an additional, larger thruster, and pushes the rocket booster in excess of 400 m. Appropriate motor control gains on the robotic arm are found to be higher during the pushing phase than those which are appropriate during free motion; this increase promotes rigidity of the arm and allows it to properly direct the pushing force. Challenges encountered in time step selection for numerical stability of the simulation are also discussed.