Automated optimization of residual reduction algorithm parameters in OpenSim

Abstract

The residual reduction algorithm (RRA) in OpenSim is designed to improve dynamic consistency of kinematics and ground reaction forces in movement simulations of musculoskeletal models. RRA requires the user to select numerous tracking weights for the joint kinematics to reduce residual errors. Selection is often performed manually, which can be time-consuming and is unlikely to yield optimal tracking weights. A multi-heuristic optimization algorithm was used to expedite tracking weight decision making to reduce residual errors. This method produced more rigorous results than manual iterations and although the total computation time was not significantly reduced, this method does not require the user to monitor the algorithm's progress to find a solution, thereby reducing manual tuning. Supporting documentation and code to implement this optimization is freely provided to assist the community with developing movement simulations.

Keywords: Biomechanics; Dynamic Consistency; MATLAB; Movement simulation; OpenSim; Python; RRA.