On the nonlinear modeling of systems with higher-order nonholonomic constraints