A physics engine for constrained multi-flexible-body dynamics