Nonholonomic Systems - Dynamics Of
This non-integrable velocity constraint is the hallmark of a nonholonomic system. The skateboard can access all possible $(x, y, \theta)$ configurations—no positional restriction—but it cannot move arbitrarily between them. Its velocity is constrained at every instant. In holonomic systems, we can reduce the problem: express velocities in terms of a smaller set of generalized coordinates and their derivatives. Lagrange’s equations then apply directly.
Welcome to the world of , where the rules of classical mechanics get a subtle, often counterintuitive, twist.
The resulting equations of motion are:
[ \frac{d}{dt} \left( \frac{\partial L}{\partial \dot{q}^j} \right) - \frac{\partial L}{\partial q^j} = \lambda_i a^i_j(q) ] dynamics of nonholonomic systems
Most introductory physics courses teach constraints through the lens of a bead on a wire or a pendulum. These are holonomic constraints: they reduce the number of independent coordinates (degrees of freedom) needed to describe the system. A bead on a fixed wire has 1 degree of freedom instead of 3. Simple.
This leads to the , which differs from the standard Euler-Lagrange equations in a crucial way: the constraint forces do no work under virtual displacements, but real displacements (which must satisfy the constraints) may still lead to energy-conserving but non-integrable motion.
This is a differential equation. Can you integrate it to find a relationship between $x, y,$ and $\theta$ alone? No. Because you can change the skateboard’s orientation without changing its position (spin in place), and you can move it along a closed loop and return to the same orientation but a different position (think parallel parking). This non-integrable velocity constraint is the hallmark of
The Lie brackets of constraint vector fields generate directions not initially allowed. That’s why you can parallel park: the bracket of “move forward” and “turn” gives “sideways slide” at the Lie algebra level, and through a sequence of motions, you achieve net motion in the forbidden direction.
where $a^i_j$ are coefficients of the velocity constraints $\sum_j a^i_j(q) \dot{q}^j = 0$, and $\lambda_i$ are Lagrange multipliers.
And yet, at the fundamental level, they remind us that constraints in physics are not merely simplifications—they are active shapers of possibility. The wheel that refuses to slip, the blade that refuses to slide, the ice skater’s edge—all carve out a geometry of motion richer than any set of fixed coordinates can capture. In holonomic systems, we can reduce the problem:
In nonholonomic dynamics, the map is not the territory. The path is not reducible to positions. And the dance is, quite literally, in the derivatives. If you’d like to go further: look into the “Chaplygin sleigh,” “rolling penny,” or the “nonholonomic integrator” in geometric numerical integration. The rabbit hole is deep, and the wheels never slip.
In nonholonomic systems, we cannot. The constraints are linear in velocities, so we can use Lagrange multipliers to enforce them. But here’s the deep part: (in the ideal case). That means D’Alembert’s principle still holds—but only for virtual displacements consistent with the constraints.
Imagine trying to push a shopping cart sideways. No matter how hard you push, it stubbornly resists, rolling only forward or backward. Or consider a car on an icy road: you can turn the wheels, but the car might continue sliding straight. Contrast this with a helicopter’s swashplate or a cat falling upright. These are not just different problems in mechanics—they represent a fundamental split in how constraints shape motion.
But nonholonomic constraints are different. They restrict the velocities of a system, not its positions, in a way that cannot be integrated into a positional constraint. The classic example? A rolling wheel without slipping. Take a skateboard. Its position in the plane is given by $(x, y)$ and its orientation by $\theta$. That’s 3 degrees of freedom. Now impose the “no lateral slip” condition: the wheel’s velocity perpendicular to its orientation must be zero.