In order to model the behaviour of muscles that work around the joints in musculoskeletal simulations, information about the muscle paths are needed. Typically, musculotendon paths, their lengths, and their force directions are modeled as the shortest connection wrapping around obstacles representing bones and adjacent tissue. In this work, the shortest path problem is described via constrained variational dynamics and is applied to a biomechanical example with muscle paths on non-uniform rational B-spline (NURBS) surfaces representing the wrapping obstacles. The results are compared with muscle paths from a G1-continuous combination of geodesics [2,4,5].