SUMMARYThis paper discusses optimization algorithms in movement simulations for models of humans, humanoid robots or other mechanisms. Targeted movements between two configurations define a dynamically redundant system, for which there is freedom in the choice of control force time variations. A previously developed formulation for the treatment of targeted dynamics for mechanisms was used as a basis. The paper describes the development of an algorithm related to the method of moving asymptotes for the necessary optimization. The algorithm is specifically adapted to problems which are large and non-linear but sparse, and which include very high numbers of design variables as well as constraints. In particular, non-linear equality constraints from dynamic equilibrium equations are important. The optimization algorithm was developed to include these, but also in order to allow successively increasing penalty factors for constraint violations. The resulting setting was shown to be able to handle the systems established, robustly giving convergence to at least a local minimum also for very distant start iterates. The existence of very closely situated local optima, representing very similar movements, was discovered for the problem formulation, calling for an ad hoc method for finding the best of these local optima.