Dynamics Of Nonholonomic Systems (2025)

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.

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.

Welcome to the world of , where the rules of classical mechanics get a subtle, often counterintuitive, twist. dynamics of nonholonomic systems

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. This non-integrable velocity constraint is the hallmark of

The resulting equations of motion are:

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 holonomic systems, we can reduce the problem:

[ \frac{d}{dt} \left( \frac{\partial L}{\partial \dot{q}^j} \right) - \frac{\partial L}{\partial q^j} = \lambda_i a^i_j(q) ]

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.

[ \dot{x} \sin \theta - \dot{y} \cos \theta = 0 ]

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.