A variational model for data fitting on manifolds by minimizing the acceleration of a Bézier curve