Weakly Nonholonomic Systems

• Weakly nonholonomic systems are those where ideal nonholonomic constraints are approximated via singular limits, friction, or higher-order Lie brackets.
• They encompass multiple regimes—such as Tatarinov-type, strong-friction, and control-based models—where constraint activation evolves gradually.
• Applications range from geometric mechanics and stochastic perturbations to control accessibility, revealing subtle drift and reduced integrability under perturbations.

Searching arXiv for recent and foundational papers relevant to weakly nonholonomic systems. {"query":"weakly nonholonomic systems arXiv Tatarinov stochastic nonholonomic constraints strong friction weakly degenerate Chaplygin", "max_results": 10} {"query":"weakly nonholonomic systems Tatarinov stochastic nonholonomic constraints strong friction weakly degenerate Chaplygin", "max_results": 10} Weakly nonholonomic systems are mechanical or control systems in which nonholonomic structure appears in a singular, approximate, or low-order sense rather than as a fixed ideal distribution. In the cited literature, the phrase is used in more than one technical meaning: constraints may depend on a small parameter and become integrable at the singular limit; large viscous friction may realize an exponentially attracting slow manifold that is nearly nonholonomic but permits small slip; or driftless control-affine systems may require Lie brackets up to second order to recover full accessibility. Closely related developments treat stochastic deformations of constraints, weakly degenerate Chaplygin reduction, and weak variational formulations for impacts, which together place weak nonholonomy within geometric mechanics, singular perturbation theory, and nonlinear control (Kuleshov et al., 14 Aug 2025, Gzenda et al., 2024, 1811.09120).

1. Terminological scope and principal meanings

A common source of confusion is that “weakly nonholonomic” does not denote a single definition across the literature. The term labels distinct, though related, mathematical regimes.

Sense	Defining feature	Source
Tatarinov-type weak nonholonomy	At $\varepsilon=0$, the constraints become completely integrable	(Kuleshov et al., 14 Aug 2025)
Strong-friction weak nonholonomy	For small $\varepsilon$ and large $\mu=1/\varepsilon$, the system is nearly nonholonomic with $h_\varepsilon=O(\varepsilon)$ slip	(Gzenda et al., 2024)
Second-degree weak nonholonomy	Lie brackets up to order $2$ span the tangent space	(1811.09120)
Weakly degenerate Chaplygin systems	$g$ is degenerate on $TQ$ but $g|_{\Delta_Q}$ is positive-definite and the reduced Lagrangian is non-degenerate	(Leok et al., 2011)

In the Tatarinov framework, one studies linear velocity constraints

$$f_s(x,\dot x,\varepsilon)=\sum_{i=1}^{n+m} a_{si}(x,\varepsilon)\,\dot x_i=0,\qquad s=1,\dots,m,$$

with $\operatorname{rank}(a_{si})=m$, and declares the system weakly nonholonomic when the constraints become completely integrable at $\varepsilon$0. In the strong-friction framework, the terminology refers to the regime in which viscous friction with large coefficients realizes a slow manifold $\varepsilon$1 close to the ideal nonholonomic distribution $\varepsilon$2. In the control-theoretic framework, “weakly nonholonomic of degree 2” refers not to approximate integrability but to controllability generated by vector fields together with first- and second-order Lie brackets.

This multiplicity of meanings implies that the phrase is best interpreted relative to the analytic mechanism under study: singular integrable limit, singular frictional realization, or bracket-generating accessibility. A plausible implication is that the unifying theme is not a single definition but a family of “near-nonholonomic” regimes in which ideal nonholonomic motion is recovered only after reduction, averaging, or bracket closure.

2. Small-parameter weak nonholonomy and the transgression effect

In the sense developed from Tatarinov’s theory, a mechanical system with $\varepsilon$3 generalized coordinates $\varepsilon$4 and Lagrangian $\varepsilon$5 is weakly nonholonomic when, at $\varepsilon$6, its $\varepsilon$7 linear constraints become completely integrable. Precisely, there exist $\varepsilon$8 independent functions $\varepsilon$9 and an invertible matrix $\mu=1/\varepsilon$0 such that

$\mu=1/\varepsilon$1

Hence the singular limit integrates to

$\mu=1/\varepsilon$2

and the original system reduces to an $\mu=1/\varepsilon$3-degree-of-freedom holonomic system depending on the parameters $\mu=1/\varepsilon$4. After introducing adapted coordinates

$\mu=1/\varepsilon$5

the constraints take the form

$\mu=1/\varepsilon$6

or equivalently

$\mu=1/\varepsilon$7

The reduced $\mu=1/\varepsilon$8 dynamics is assumed to be a completely integrable Hamiltonian system with Hamiltonian $\mu=1/\varepsilon$9 and action-angle variables $h_\varepsilon=O(\varepsilon)$0 satisfying

$h_\varepsilon=O(\varepsilon)$1

For $h_\varepsilon=O(\varepsilon)$2, the formerly frozen coordinates $h_\varepsilon=O(\varepsilon)$3 become slow variables, and to first order one obtains

$h_\varepsilon=O(\varepsilon)$4

In the Hamiltonian form,

$h_\varepsilon=O(\varepsilon)$5

and averaging over the fast torus yields the leading-order normalized equations

$h_\varepsilon=O(\varepsilon)$6

The paper identifies this slow $h_\varepsilon=O(\varepsilon)$7 drift of the former integration constants as the transgression effect (Kuleshov et al., 14 Aug 2025).

The worked example with coordinates $h_\varepsilon=O(\varepsilon)$8, quadratic potential $h_\varepsilon=O(\varepsilon)$9, and weak constraint

$2$0

shows the mechanism explicitly. At $2$1, $2$2 is constant and $2$3 execute decoupled harmonic oscillations. For $2$4, the in-plane motion remains

$2$5

while

$2$6

Thus the holonomic approximation is not stationary on times $2$7; instead, the system remains close to an integrable family with slowly drifting parameters. The stated validity regime is “small” nonintegrability $2$8 and time intervals up to $2$9.

3. Strong friction, slow manifolds, and weak realization of constraints

A different usage of weak nonholonomy arises when ideal nonholonomic constraints are realized as the singular limit of strong viscous friction. The geometric setup consists of a configuration manifold $g$0, a Riemannian metric $g$1, a nonholonomic distribution $g$2 defined by linearly independent one-forms $g$3, the orthogonal complement $g$4, and a Rayleigh dissipation function

$g$5

With large friction scaled by $g$6, the forced affine-connection equations are

$g$7

where $g$8. In the limit $g$9, the fast dynamics drive $TQ$0, hence $TQ$1. For $TQ$2, there exists an exponentially attracting slow manifold $TQ$3, represented by a nonlinear section

$TQ$4

with decomposition

$TQ$5

The paper proposes a novel invariance condition based on covariant derivatives and proves that it is equivalent to the classical invariance condition based on time derivatives. The covariant formulation is

$TQ$6

From this one obtains a coordinate-free generating equation for $TQ$7 and a recursive expansion

$TQ$8

with

$TQ$9

and an explicit formula for $g|_{\Delta_Q}$0. The associated slip-velocity approximation is

$g|_{\Delta_Q}$1

while the zeroth-order reduced dynamics reproduces the standard nonholonomic affine-connection equation

$g|_{\Delta_Q}$2

The first-order correction introduces viscous drift terms (Gzenda et al., 2024).

The vertical rolling disk illustrates the construction concretely. With $g|_{\Delta_Q}$3, coordinates $g|_{\Delta_Q}$4, flat metric $g|_{\Delta_Q}$5, and constraints

$g|_{\Delta_Q}$6

the distribution $g|_{\Delta_Q}$7 is spanned by

$g|_{\Delta_Q}$8

The first-order slip term is

$g|_{\Delta_Q}$9

and the $$f_s(x,\dot x,\varepsilon)=\sum_{i=1}^{n+m} a_{si}(x,\varepsilon)\,\dot x_i=0,\qquad s=1,\dots,m,$$0 dynamics reproduces rolling without slip. The paper states that, in the small-$$f_s(x,\dot x,\varepsilon)=\sum_{i=1}^{n+m} a_{si}(x,\varepsilon)\,\dot x_i=0,\qquad s=1,\dots,m,$$1 large-$$f_s(x,\dot x,\varepsilon)=\sum_{i=1}^{n+m} a_{si}(x,\varepsilon)\,\dot x_i=0,\qquad s=1,\dots,m,$$2 regime, the system behaves nearly nonholonomically with $$f_s(x,\dot x,\varepsilon)=\sum_{i=1}^{n+m} a_{si}(x,\varepsilon)\,\dot x_i=0,\qquad s=1,\dots,m,$$3 and explicitly calls this weakly nonholonomic. Its stated validity condition is that $$f_s(x,\dot x,\varepsilon)=\sum_{i=1}^{n+m} a_{si}(x,\varepsilon)\,\dot x_i=0,\qquad s=1,\dots,m,$$4 depending on geometry so that the $$f_s(x,\dot x,\varepsilon)=\sum_{i=1}^{n+m} a_{si}(x,\varepsilon)\,\dot x_i=0,\qquad s=1,\dots,m,$$5 terms remain small over the time-scale of interest.

4. Second-degree weak nonholonomy in nonlinear control

In nonlinear control, weak nonholonomy refers to a bracket-generation property of a driftless control-affine system

$$f_s(x,\dot x,\varepsilon)=\sum_{i=1}^{n+m} a_{si}(x,\varepsilon)\,\dot x_i=0,\qquad s=1,\dots,m,$$6

with $$f_s(x,\dot x,\varepsilon)=\sum_{i=1}^{n+m} a_{si}(x,\varepsilon)\,\dot x_i=0,\qquad s=1,\dots,m,$$7. The system is called weakly nonholonomic of degree $$f_s(x,\dot x,\varepsilon)=\sum_{i=1}^{n+m} a_{si}(x,\varepsilon)\,\dot x_i=0,\qquad s=1,\dots,m,$$8, or second-degree nonholonomic, if Lie brackets of order up to $$f_s(x,\dot x,\varepsilon)=\sum_{i=1}^{n+m} a_{si}(x,\varepsilon)\,\dot x_i=0,\qquad s=1,\dots,m,$$9 span the tangent space at every $\operatorname{rank}(a_{si})=m$0. Equivalently, there exist index sets

$\operatorname{rank}(a_{si})=m$1

with $\operatorname{rank}(a_{si})=m$2 such that

$\operatorname{rank}(a_{si})=m$3

The corresponding rank condition is

$\operatorname{rank}(a_{si})=m$4

where $\operatorname{rank}(a_{si})=m$5 is the $\operatorname{rank}(a_{si})=m$6 matrix formed by the selected vector fields, first brackets, and second brackets. This is the Lie-algebra-rank condition up to second order (1811.09120).

The control objective in the cited work is obstacle avoidance. Obstacles are embedded into the free space

$\operatorname{rank}(a_{si})=m$7

and one chooses a smooth navigation function $\operatorname{rank}(a_{si})=m$8 whose sublevel sets stay away from $\operatorname{rank}(a_{si})=m$9, with a unique minimum at the goal $\varepsilon$00 and growth to “$\varepsilon$01” at obstacle boundaries. The target unconstrained dynamics is

$\varepsilon$02

To approximate this with the underactuated system, the feedback is chosen in explicit time-varying form

$\varepsilon$03

with coefficients

$\varepsilon$04

For small $\varepsilon$05,

$\varepsilon$06

The constant terms steer along $\varepsilon$07, the $\varepsilon$08 oscillations generate first brackets, and the $\varepsilon$09 bi-periodic products generate second brackets. A non-resonance condition on the chosen frequencies ensures that no spurious bracket directions appear.

The main theorem states that, under the second-degree rank condition and mild smoothness and boundedness assumptions, there exists $\varepsilon$10 such that for all $\varepsilon$11, the sampled solution remains in $\varepsilon$12 for all $\varepsilon$13 and converges to

$\varepsilon$14

If $\varepsilon$15 is a navigation function with a single non-degenerate minimum at $\varepsilon$16, then $\varepsilon$17. The proof uses a one-step Volterra expansion yielding

$\varepsilon$18

followed by a discrete-time LaSalle argument. Here weak nonholonomy is therefore an accessibility property generated by iterated brackets, not a small-slip or small-parameter perturbation of a mechanical constraint.

5. Stochastic deformations of nonholonomic constraints

A nearby line of research studies stochastic perturbations of nonholonomic constraints through a stochastic extension of the Lagrange–d’Alembert framework. Two Stratonovich formulations are introduced. In the affine stochastic constraint, one replaces

$\varepsilon$19

by

$\varepsilon$20

The corresponding stochastic Lagrange–d’Alembert equations are

$\varepsilon$21

In the ideal stochastic constraint, one imposes

$\varepsilon$22

and obtains

$\varepsilon$23

In the ideal case, the mechanical energy

$\varepsilon$24

is exactly preserved along Stratonovich trajectories, provided $\varepsilon$25 do not do work. In both cases, elimination of the multipliers yields Stratonovich SDEs on $\varepsilon$26, and the constraint force picks up extra terms from $\varepsilon$27 or from $\varepsilon$28, thereby coupling noise into the momentum evolution (Gay-Balmaz et al., 2017).

The geometric interpretation is given in terms of stochastic connections and curvature. Deterministic nonholonomic constraints define an Ehresmann connection $\varepsilon$29 with horizontal distribution $\varepsilon$30. In the affine stochastic case, the constraint distribution becomes the affine shift

$\varepsilon$31

so $\varepsilon$32 is replaced by $\varepsilon$33 and holonomy along loops acquires a stochastic line integral of $\varepsilon$34. In the ideal stochastic case, the perturbed connection $\varepsilon$35 yields stochastic curvature

$\varepsilon$36

On Lie groups with left-invariant constraints, these deformations give stochastic Euler–Poincaré–Suslov equations on $\varepsilon$37.

The paper develops these constructions for the stochastic Suslov problem on $\varepsilon$38 and $\varepsilon$39. For example, the affine stochastic Suslov constraint is

$\varepsilon$40

while the ideal version is

$\varepsilon$41

The analysis then tracks the fate of classical integrals. For the Kharlamova integral $\varepsilon$42, the affine case is nonconservative in general, whereas in the ideal case one has

$\varepsilon$43

and conservation follows if $\varepsilon$44. For the Clebsch–Tisserand integral $\varepsilon$45, the affine case gives

$\varepsilon$46

while the ideal case gives

$\varepsilon$47

The qualitative discussion distinguishes the two noise models sharply: in the affine case the energy generally drifts and the system may slide off the ideal constraint manifold, whereas in the ideal case the energy is exactly conserved by construction. This suggests that stochastic deformation provides a controlled way to study how weak violations of ideal constraints alter holonomy, momentum drift, and long-time transport.

6. Geometric, Hamilton–Jacobi, and impact formulations

The broader theory of weakly nonholonomic phenomena intersects with two additional frameworks: Dirac reduction for degenerate constrained systems and weak variational formulations for impacts. In the Lagrange–Dirac setting, a regular constraint distribution $\varepsilon$48 and a possibly degenerate Lagrangian $\varepsilon$49 define an induced Dirac structure $\varepsilon$50. A Lagrange–Dirac system is a triple $\varepsilon$51 such that

$\varepsilon$52

which in coordinates yields the implicit Euler–Lagrange equations

$\varepsilon$53

The associated Dirac–Hamilton–Jacobi equation for a section

$\varepsilon$54

is

$\varepsilon$55

with generalized energy $\varepsilon$56. When $\varepsilon$57 is connected and $\varepsilon$58 is bracket-generating, this is equivalent to $\varepsilon$59 (Leok et al., 2011).

Within this framework, the cited paper singles out weakly degenerate Chaplygin systems: Chaplygin systems for which

$\varepsilon$60

has degenerate $\varepsilon$61 on $\varepsilon$62, but $\varepsilon$63 is positive-definite so that the reduced Lagrangian $\varepsilon$64 is non-degenerate. Reduction yields an almost Hamiltonian system on $\varepsilon$65,

$\varepsilon$66

where the non-closed two-form $\varepsilon$67 encodes the curvature of the principal connection and the momentum map. The reduced Dirac–HJ problem requires

$\varepsilon$68

This is not itself a definition of weak nonholonomy, but it supplies an exact integration framework for systems in which degeneracy and nonholonomic reduction coexist.

A different extension is the weak-solution theory for collisions in nonholonomic systems. There, one replaces classical trajectories by curves $\varepsilon$69 satisfying the constraint almost everywhere and the weak integral identity

$\varepsilon$70

for all compactly supported virtual displacements $\varepsilon$71 tangent to the constraint distribution. If a collision occurs on a submanifold $\varepsilon$72 with stronger instantaneous constraint $\varepsilon$73, then the momentum jump satisfies

$\varepsilon$74

equivalently $\varepsilon$75 with $\varepsilon$76. Under energy conservation and reversibility, the unique rebound law is

$\varepsilon$77

where $\varepsilon$78 is the $\varepsilon$79-orthogonal projector onto $\varepsilon$80 (Treschev et al., 2012).

These frameworks broaden the meaning of “weak” in nonholonomic mechanics. In one direction, weak degeneracy allows reduction to almost Hamiltonian form and exact Hamilton–Jacobi integration. In another, weak solutions permit discontinuous velocities and collision laws inside the Lagrange–d’Alembert formalism. Taken together with the small-parameter, strong-friction, second-degree, and stochastic settings, they show that weakly nonholonomic systems are best understood as a family of singular or generalized limits in which ideal nonholonomic structure is either approached, perturbed, or reconstructed rather than imposed once and for all.