Abstract This paper presents a joint–coordinate adjoint method for optimal control of multi-rigid-body systems. Initially formulated as set differential-algebraic equations, the system is brought into minimal form by projecting original expressions joint’s motion and constraint force subspaces. Consequently, cumbersome partial derivatives corresponding to joint-space equations are avoided, appr...