## Affine Connection Approach to the Realization of Nonholonomic Constraints by Strong Friction Forces

Vaughn Gzenda and Robin Chhabra*

Department of Mechanical and Aerospace Engineering, Carleton University, 1125 Colonel By Dr, Ottawa, K1S 5B6, Ontario, Canada.

*Corresponding author(s). E-mail(s): robin.chhabra@carleton.ca; Contributing authors: vaughngzenda@cmail.carleton.ca;

## Abstract

In this paper, we study an affine connection approach to realizing nonholonomic mechanical systems mediated by viscous friction forces with large coefficients, viewed as a singular perturbation of the nonholonomic system. We show that the associated slow manifold is represented coordinate-free as the image of a section over the nonholonomic distribution. We propose a novel invariance condition based on covariant derivatives and prove that this condition is equivalent to the classical invariance condition based on time derivatives. Accordingly, we propose a novel recursive procedure to approximate the slow manifold based on the covariant derivatives of a formal power series expansion of the section. Using this recurrence relation, we derive, up to second order, approximations of the slip velocities residing in the slow manifold, as well as the associated approximated dynamics up to first order. Lastly, we illustrate our approach with a case study of a vertical rolling disk.

Keywords: Nonholonomic mechanics, Perturbation methods, Reduced-order modelling, Riemannian Geometry

Submitted to the Journal of Nonlinear Dynamics

## 1 Introduction

Rolling without slipping is a constraint that appears in many mechanical systems, such as vehicle systems, and is commonly viewed as an example of an ideal nonholonomic constraint [1]. Nonholonomic constraints normally refer to point-wise linear constraints on the velocity space at each configuration, which cannot be described by a derivative of a set of constraints on the configuration space [1]. Hence, such constraints are then viewed as a nonintegrable linear distribution of the tangent bundle to the configuration manifold [1]. However, the nonholonomic constraint forces are only valid for velocities in the distribution. In practice, nonholonomic constraints in vehicle systems are often violated, resulting in an accumulation of errors for standard nonholonomic model-based control laws. To model the velocities violating the nonholonomic constraints, the ideal nonholonomic constraint forces may be modelled as dissipative friction, whose limit at large coefficients approximates the ideal nonholonomic system [2]. The large friction coefficients introduce a steep linear relationship between the friction force and the motion of the contact point that is analogous to the linear regime of Bakker-Pacejka magic formula [3]. Mathematically, the system with large coefficients may be treated as a slow/fast singular perturbation where the dissipative friction forces act on a much faster time scale than the nominal variables resulting in an exponentially attractive invariant manifold. This manifold - often called the slow manifold - physically represents small violations of the nonholonomic constraints.

Nonholonomic mechanical systems form a class of constrained systems whose velocity constraints have a natural differential geometric structure [1, 4-6]. Classically, the equations of motion are derived using Lagrange-d'Alembert's (L-d'A) variational principle -which is an adaptation of Hamilton's variational principle -whereby the variations are constrained to the nonholonomic distribution and enforced with Lagrange multipliers that do no work on the distribution [1]. Some examples of nonholonomic constraints include rolling without slipping, no sliding, and conserved angular momenta, among others [1, 7]. The kinetic energy of the mechanical system may be interpreted geometrically as a Riemannian metric tensor [8]. Equivalent to the Lagrange-d'Alembert principle, Lewis and Bullo [9, 10] use Riemannian geometry to model the nonholonomic dynamics with affine connections. An affine connection is a geometric object that enables a notion of differentiation of vector fields on a smooth manifold. A unique affine connection is derived from the metric tensor called the Levi-Civita affine connection that encodes Coriolis, centrifugal, and centripetal accelerations. In [7], the nonholonomic constraints may be viewed as the kernel of a projection map onto the complement of the nonholonomic distribution. In the affine connection framework, [7], the nonholonomic constraints are then enforced using Lagrange multipliers that are solved for using a covariant derivative of the complementary projection map. As such, Lewis introduces an affine connection called the nonholonomic connection, based on the Levi-Civita connection and the covariant of the orthogonal projection map, modelling the necessary accelerations required to maintain the nonholonomic constraints. The geometric interplay of affine connections and nonintegrable distributions is studied in detail in [9]. Other equivalent approaches to nonholonomic mechanics include symplectic geometry [4, 11], Poisson geometry [12], Gauss's Principle of Least Constraint [1, 7] and Kane's equations [13].

Singularly perturbed systems (SPS) manifest in various domains of science and engineering, emerging when a dynamical system relies on scalable parasitic parameters like small time constants, masses, friction, etc . Accordingly, SPS have natural applications in electric circuits [14], chemistry [15], neuroscience [16], among others. Slow/fast systems are a particular case of a SPS that involve two distinct time scales, delineated by a singular parameter /epsilon1 . Tikhonov established that when the spectra of the linearized fast dynamics lies in the left half-plane, the flow converges to an invariant (slow) manifold as /epsilon1 tends to 0 [17]. The functional relationship between the slow and fast variables is dictated by an invariance condition and a formal power series expansion is a common approach to approximate the slow manifold [18]. Fenichel examined the persistence of the invariant manifold under small perturbations of ε , utilizing the concept of Normally Hyperbolic Invariant Manifold (NHIM) and under the assumption of compactness in the state space [19]. In [20], Fenichel utilizes a geometric approach to investigate asymptotic expansions of the graph representation of the invariant manifold. In the Riemannian geometric context, Eldering extends Fenichel's findings to non-compact state spaces, assuming regularity under bounded geometry conditions [21]. Eldering's work opens up a pathway for the coordinate-independent examination of such systems. If a mechanical system possesses a noncompact configuration space Q and is endowed with a symmetry group G , where the group action leads to a compact shape space Q/G [22], then it follows that Q has bounded geometry [2, 21]. In practical applications, numerous mechanical systems exhibit symmetry groups, as seen in mobile vehicles that yield compact shape spaces [11]. Hence, Eldering's geometric singular perturbation techniques [21] can be employed for the analysis of such systems.

The physical realization of nonholonomic mechanical systems by strong viscous friction forces was treated classically by Brendelev and Karapetian who mainly focused on convergence to the ideal nonholonomic dynamics [23, 24]. More recently, Eldering used the geometric singular perturbation theory of NHIM's [21] to generalize these convergence results for systems with bounded geometry. Common to all these models, the contact forces are modelled as a Rayleigh dissipation function (RDF) defined by a symmetric (0 , 2) -tensor. Rayleigh dissipation functions have been used to model surface contact forces following first physical principles [25]. To realize the nonholonomic constraints, the RDF is constructed such that its kernel is identified with the nonholonomic distribution ensuring that the contact forces do no work on the distribution, mimicking the L-d'A principle, yet the dynamics are exponentially attractive to the distribution. At the limit of large friction coefficients, Eldering proves the existence of a persistent exponentially attractive NHIM and shows that in a nonstandard frame spanning the distribution and its complement, the slow manifold may be represented as the graph over the nonholonomic variables [2]. It is shown that this graph may be represented as a formal power series, which is used to approximate the slow manifold up to first order [2] based on the classical invariance condition.

Antali and Stepan considered the physical realization of the rolling phenomenon of wheels by dry friction forces [26]. Their analysis is based on the crossing dynamics of extended Filippov systems [27]. A Filippov system is described by a piece-wise smooth vector field with discontinuities on a discontinuity manifold [28]. Extended

Filippov (EF) systems are Filippov systems with codimension-2 discontinuity manifolds and describe Coulomb rigid body dynamics [27]. The EF approach has been used to analyze the spinning and sliding behavior of rigid bodies under contact [29]. Recently, bifurcation analysis of EF systems has also been developed to analyze slip and sliding dynamics of a wheel motion [30]. In contrast to the slow-fast realization by viscous friction forces [2], the Coulomb friction in [26] is both nonsmooth and nonlinear in velocities. The wheel rolling phenomenon has also been captured as a limit of piece-wise continuous holonomic constraints [31]. Other historical approaches to the realization of nonholonomic constraints include coupling of an external force to generate nonholonomic dynamics [32].

In this paper, we study realizations of nonholonomic systems as slow/fast systems using an affine connection approach inspired by the work of Bullo and Lewis [7, 9, 10]. In our geometric approach, we represent the slow manifold coordinate-independent as the image of a section over the nonholonomic distribution. We introduce a novel invariance based on a covariant derivative along the integral curves on the slow manifold. We prove that this novel invariance condition is equivalent to the invariance condition based on time derivatives [2]. We propose a recursion to approximate the velocities in the slow manifold based on covariant derivatives of geometric quantities. Using the recursion, we derive approximations of equations of motion. Lastly, we study an example system of a vertical rolling disk. The major contributions of this work are summarized as follows:

- 1. formulating a coordinate-free representation of the slow manifold as the image of a nonlinear section over the nonholonomic distribution,
- 2. formulating a novel invariance condition based on covariant derivatives, and proving its equivalence to the classical invariance condition,
- 3. proposing a novel recurrence relation to approximate the slip velocities described by the slow manifold,
- 4. using the recurrence relation to derive higher-order perturbations of the equations of motion.

The paper is organized as follows. Section 2 provides a self-contained summary of the differential geometric structures used in the remainder of the article. This includes a review of Riemannian geometry, affine connection systems, nonholonomic affine connection systems, and the geometry of the second tangent bundle including the procedure to take covariant derivatives of general (nonlinear) tangent bundle maps. In Section 3, we derive the equations of motion for a slow fast realization of nonholonomic systems and represent the slow manifold as the image of a section. Using this model, in Section 4, we derive the relationship between the reaction forces and the slow manifold and calculate zeroth, first and second order approximations of the slip velocities, and zeroth and first-order dynamical equations of motion. Lastly, in Section 5 we consider our approach with a case study of a vertical rolling disk.

## 2 Preliminaries

This section provides a self-contained review of mathematical notions used throughout the paper. For references on Riemannian geometry, geometric mechanics and related topics see [1, 7, 8, 10, 33, 34] and references therein.

## 2.1 Riemannian Geometry and Affine Connection Systems

Let Q be an n -dimensional smooth manifold representing the configuration space of a simple mechanical system. Let TQ = ⊔ q ∈ Q T Q q and T ∗ Q = ⊔ q ∈ Q T ∗ q Q be the disjoint union of tangent and cotangent spaces called the tangent and cotangent bundle, respectively. Let π TQ : TQ → Q be the bundle projection map for TQ and let X TQ be the set be the set of vector fields on Q . In local coordinates ( U, q = ( q 1 , . . . , q n ) ) , we denote the coordinate induced basis by { ∂ ∂ q 1 , . . . , ∂ ∂ q n } for T Q q and { d q 1 , . . . , d q n } for T ∗ q Q . The natural pairing 〈 | 〉 . . : T ∗ q Q × T Q q → R between vectors and covectors in the tangent and cotangent spaces, respectively, is defined by

$$\langle \omega _ { q } | X _ { q } \rangle = \sum _ { i = 1 } ^ { n } \omega _ { i } X ^ { i }, \quad \omega _ { q } \in T _ { q } ^ { * } Q, \ X _ { q } \in T _ { q } Q. \quad \ \ ( 1 ) \\. \ \wedge \quad \.$$

where X q = ∑ i X i ∂ ∂ q i and ω q = ∑ i ω i d q i in the canonical basis of the tangent and cotangent bundles, having components X ,.. . , X 1 n ∈ R and ω , . . . , ω 1 n ∈ R , respectively.

Let G q : T Q q × T Q q → R be a Riemannian metric on Q defining the kinetic energy of the system [8]. For each q ∈ Q the metric G induces two isomorphisms between the tangent and cotangent bundles given point-wise by

↦

$$& \mathfrak { f } _ { \mathbb { G } } \, \colon T _ { q } ^ { * } Q \to T _ { q } Q \\ & \omega _ { q } \mapsto \mathbb { G } ( \omega _ { q } ^ { \mathfrak { f } _ { \mathbb { G } } }, Y _ { q } ) = \langle \omega _ { q } | Y _ { q } \rangle \,, \quad \forall Y _ { q } \in T _ { q } Q$$

↦

called the sharp and flat isomorphisms, respectively [7, 8]. In the canonical local coordinates for TQ and T ∗ Q , the metric tensor takes the matrix form such that by abuse of notation we denote the matrix by G q . As the result, the sharp and flat isomorphisms become ω /sharp G q = G -1 q ω q and X /flat G q = G q X q , where ω q and X q are expressed in the canonical bases.

$$& \mu _ { q } \cdot \lambda \vdash \mu _ { q }, \nu _ { q / } \cdot \lambda \vdash \mu _ { q / }, \cdot \lambda \vdash \mu _ { q } \lambda \vdash \\ & \flat _ { \mathbb { G } } \, \colon T _ { q } Q \to T _ { q } ^ { * } Q \\ & X _ { q } \mapsto \left \langle X _ { q } ^ { \flat _ { \mathbb { G } } } | Y _ { q } \right \rangle = \mathbb { G } _ { q } ( X _ { q }, Y _ { q } ), \quad \forall Y _ { q } \in T _ { q } Q \\ \text{the sharp and flat isomorphisms. respectively} \, | 7. \, 8 l. \, In the canonical local coor-$$

The Levi-Civita connection ∇ , is the unique affine connection associated with G and is derived from the assumption that the connection is divergentless and torsionfree [7, 8]. The ∇ X Y be the covariant derivative of a vector field Y along a vector field X with respect to ∇ is given by

$$\begin{array} { c c c } \cdot & \cdot & \cdot \\ & \cdot & \cdot \\ & \nabla _ { X } Y = \left ( \frac { \partial Y ^ { i } } { \partial q ^ { j } } X ^ { j } + \Gamma _ { j k } ^ { i } X ^ { j } Y ^ { k } \right ) \frac { \partial } { \partial q ^ { i } }, \end{array}$$

where the Levi-Civita connection coefficients - also called the Chistoffel symbolsΓ i jk in the coordinate-induced basis is expressed as

$$\Gamma _ { j k } ^ { i } = \frac { 1 } { 2 } \mathbb { G } ^ { i l } \left ( \frac { \partial \mathbb { G } _ { j l } } { \partial q ^ { k } } + \frac { \partial \mathbb { G } _ { k l } } { \partial q ^ { j } } - \frac { \partial \mathbb { G } _ { j k } } { \partial q ^ { l } } \right ).$$

The Christoffel symbols have the physical meaning for a simple mechanical system of Coriolis, centripetal and centrifugal accelerations [7]. The action of ∇ on (1 , 1) -tensors A ∈ X ( T ∗ Q ⊗ TQ ) (such as a projection map) is given by

$$( \nabla _ { X } \mathcal { A } ) = \left ( \frac { \partial \mathcal { A } _ { j } ^ { i } } { \partial q ^ { k } } X ^ { k } + \Gamma _ { k l } ^ { i } \mathcal { A } _ { j } ^ { l } X ^ { k } - \Gamma _ { k j } ^ { l } \mathcal { A } _ { i } ^ { i } X ^ { k } \right ) \frac { \partial } { \partial q ^ { i } } \otimes d q ^ { i }.$$

↦

We may consider the connection coefficients as the point-wise bilinear map Γ : T Q q × T Q q → T Q q as the function ( v q , w q ) → Γ( q v , q ) w q .

The Lagrangian L : TQ → R for a simple mechanical system is defined as

$$L = \frac { 1 } { 2 } \mathbb { G } _ { q } ( v _ { q }, v _ { q } ) - V ( q ).$$

The Euler-Lagrange equations of motion for Lagrangian L are equivalent to forced affine connection system ( ∇ , G , V, F )

$$\nabla _ { v ( t ) } v ( t ) = ( - d V + F ) ^ { \ddagger _ { G } } \,, \quad \frac { \text{d} } { \text{d} t } q ( t ) = v ( t ).$$

defined by the connection and the gradient of the potential and external force F ∈ T ∗ Q [7].

## 2.1.1 Nonholonomic Affine Connection System

We review the basic results of nonholonomic affine connection systems treated. For general references, see [7, 9, 10].

Let { a l ⊂ T ∗ Q | l = 1 , . . . , m } be a collection of everywhere linearly independent one-forms over Q such that they specify an ( n -m ) -dimensional nonholonomic distribution D at each point q ∈ Q by

$$\mathcal { D } _ { q } = \left \{ v _ { q } \in T _ { q } Q \ | \ \langle a _ { q } ^ { l } | v _ { q } \rangle = 0, \ \ l = 1, \dots, m \right \}. \quad \text{(9)} \\ \text{onstraint one-forms are locally given by } a _ { \sigma } ^ { l } \ = \ \sum _ { i = 1 } ^ { n } \, \text{A} _ { i } ^ { l } ( q ) d q ^ { i } \text{ for each } l \ =$$

The constraint one-forms are locally given by a l q = ∑ n i =1 A l i ( q d ) q i for each l = 1 , . . . , m . A basis for D q is chosen by selecting m everywhere linearly independent vector fields, { S 1 ( q ) , . . . , S n -m ( q ) } satisfying the constraint equation (9). Using the Riemannian metric we complete the basis for T Q q by defining

$$\mathcal { D } _ { q } ^ { \perp } \coloneqq \text{span} \{ ( a _ { q } ^ { 1 } ) ^ { \dagger _ { G } }, \dots, ( a _ { q } ^ { p } ) ^ { \dagger _ { G } } \}.$$

Therefore, each tangent space can be G -orthogonally decomposed into T Q q ∼ D ⊕D = q ⊥ q . Define the projection maps P q : T Q q →D q and P ⊥ q : T Q q →D ⊥ q at each point q ∈ Q . The nonholonomic distribution D can be described in terms of the projection maps as

$$\mathcal { D } = \{ w \in T Q \, | \, w = \mathcal { P } ( v ), \, v \in T Q \} = \{ w \in T Q \, | \, \mathcal { P } ^ { \perp } ( w ) = 0 \}. \quad \quad ( 1 1 )$$

To represent the evolution of the trajectory of the system q : [ a, b ] → Q , we form the nonholonomic affine connection system:

$$\nabla _ { v ( t ) } v ( t ) = \lambda ( t ) - ( d V ) ^ { \ddagger _ { G } } \,,$$

$$v \in \mathcal { D } \iff \mathcal { P } ^ { \perp } ( v ( t ) ) = 0,$$

$$\frac { \text{d} } { \text{d} t } q ( t ) = v ( t ),$$

where λ is a section in D ⊥ representing the Lagrange multipliers representing the constraint reaction forces given by

$$\lambda ( t ) = - \left ( \nabla _ { v } \mathcal { P } ^ { \perp } \right ) ( v ) + \mathcal { P } ^ { \perp } \left ( d V \right ) ^ { \# _ { G } }. \\ T h o m h a l o m m i o a n n o n t i o n \, I \gamma \, 1 0 \, 9 1 1 \, i o \, a t i n d \, h o \, a t h o t i n t i o n.$$

Remark 1. The nonholonomic affine connection [7, 10, 34] is obtained by substituting the Lagrange multipliers (15) into the (12) and restricting the equations of motion (12) (14) -to the distribution D .

## 2.2 Geometry of the Second-Tangent Bundle

In this section, we review the geometry of the second tangent bundle TTQ . First, we review the bundle structure and decomposition [35, 36] which is required to compute the covariant derivative of a general nonlinear vector bundle map. The covariant derivative of a general nonlinear vector bundle map is a prerequisite for the main results of this article. It is noteworthy that this specific result isn't explicitly stated in the literature. Lastly, we review the second tangent bundle in the context of geometric mechanics [35]. For general references, see [7, 22, 35, 37, 38].

## 2.2.1 Vertical/Horizontal Decomposition of TTQ

Consider the second tangent bundle of Q with bundle structure given by π TTQ : TTQ → TQ . Local coordinates for T ( q v , q ) TQ by ( q v , q , V ( q v , q ) , W ( q v , q ) ) where ( V ( q v , q ) , W ( q v , q ) ) given in the canonical coordinate induced basis by

$$V _ { ( q, v _ { q } ) } = \sum _ { i = 1 } ^ { n } V ^ { i } \frac { \partial } { \partial \tilde { q } ^ { i } } \Big | _ { ( q, v _ { q } ) }, \quad W _ { ( q, v _ { q } ) } = \sum _ { j = 1 } ^ { n } W ^ { j } \frac { \partial } { \partial v ^ { j } } \Big | _ { ( q, v _ { q } ) }. \quad ( 1 6 ) \\ \intertext { here } \tilde { a } ^ { i } \coloneqq \pi _ { \lambda \lambda } ^ { * } a ^ { i } \text{ The nsh forward of the projection man } \pi _ { \nu \lambda } \text{ indices a bundle}$$

where ˜ q i := π ∗ TQ q i . The push-forward of the projection map π TQ induces a bundle structure ( π TQ ) ∗ : TTQ → TQ . In these coordinates, the bundle projection maps

satisfy

↦

$$\pi _ { T T Q } \colon ( q, v, V, W ) \mapsto ( q, v )$$

↦

$$( \pi _ { T Q } ) _ { * } \, \colon ( q, v, V, W ) \mapsto ( q, V ) \,.$$

The vertical bundle is defined by V TQ = ker ( π TQ ) ∗ and given locally by

$$\text{me vericai number is einemo by } \nu _ { 1 } \, \mathcal { Q } = \ker ( \pi _ { T Q } ) _ { * } \, \text{and given locally by} \\ \mathcal { V } T Q = \left \{ ( q, v, V, W ) \in T T Q \, | \, \left ( \pi _ { T Q } \right ) _ { * } ( q, v, V, W ) = 0 \right \} \\ = \text{span} \left \{ \frac { \partial } { \partial v ^ { i } } \Big | _ { ( q, v _ { q } ) } | \, i = 1, \dots, n \right \}. \quad \text{(19)} \\ \text{The parallel transport equation indices a vector-valued connection one form } \, \omega \, \in \\ \mathfrak { I } ^ { 1 } \left ( T Q ; End ( \mathbb { R } ^ { n } ) \right ) \text{given by (see section 9.7 of } [37])$$

The parallel transport equation induces a vector-valued connection one form ω ∈ Ω ( 1 TQ ; End ( R n )) given by (see section 9.7 of [37])

$$\omega \coloneqq d v ^ { k } + \Gamma _ { i j } ^ { k } v ^ { j } d x ^ { i }.$$

Define the horizontal distribution H TQ = ker( ω ) ⊂ TTQ [37] which is given locally by

$$\mathcal { H } T Q = \text{span} _ { i = 1, \dots, n } \left \{ \frac { \partial } { \partial \tilde { q } ^ { i } } \Big | _ { ( q, v _ { q } ) } - \Gamma _ { i j } ^ { k } v ^ { j } \frac { \partial } { \partial v ^ { k } } \Big | _ { ( q, v _ { q } ) } \right \}. \quad \text{(21)} \\ \text{ore, the second tangent bundle is decomposed as $T T Q \cong \mathcal{H} T Q \cong \mathcal{V} T Q. Define} \\ \text{to} \, \text{if} \, \text{$1$} \, \text{out} \, \text{out} \, \text{add} \, \mathbf V \sim T \Omega \, \text{to} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \,\text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in}\, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{\,out} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \,$$

Therefore, the second tangent bundle is decomposed as TTQ ∼ H = TQ ∼ V = TQ . Define the vertical lift [10] of a vector field X ∈ TQ to the vertical space V ( q v , q ) TQ

↦

In coordinates, the vertical lift is given by X ↑ ver = X i ∂ ∂ v i . The vertical lift is a smooth bijection and hence provides an isomorphism between the V ( q v , q ) TQ and TQ given by ↓ ver := ( ↑ ver ) -1

$$\uparrow \text{ver} \colon \, T Q \to \mathcal { V } _ { ( q, v _ { q } ) } T Q \\ X \mapsto \frac { \text{d} } { \text{d} t } \left ( v _ { q } + t X \right ) \Big | _ { t = 0 }. \text{ \quad } ( 2 2 ) \\ \text{vertical lift is given by $X^{\uparrowver}=\, X^{\frac{0}{\partial v^{\uparrow}}$. The vertical lift is a } \\ | \text{hence } \text{provides an isomorphism between the $\mathcal{V}_{\wedge}$} \cdot \, \text{TO and TO}$$

↦

$$& \downarrow _ { \text{ver} } \colon \mathcal { V } _ { ( q, v _ { q } ) } T Q \to T Q \\ & ( q, v _ { q }, 0, W ) \mapsto ( q, W ) \,.$$

Similarly, we define the horizontal lift [38] with respect to the affine connection ∇ with Christoffel Symbols Γ by

↦

$$& \uparrow _ { \text{hor} } \colon T _ { q } Q \times T Q \to \mathcal { H } _ { ( q, v _ { q } ) } T Q \\ & ( q, v _ { q }, X ) \mapsto ( q, v _ { q }, X, - \Gamma ( q, v _ { q } ) X ) \,.$$

and its inverse ↓ hor := ( ↑ hor ) -1 given by the map

$$\downarrow _ { \text{hor} } \colon \, \mathcal { H } _ { ( q, v _ { q } ) } T Q \to T Q$$

↦

$$( q, v _ { q }, X, - \Gamma ( q, v _ { q } ) X ) \mapsto ( q, X ) \,.$$

Lastly, we define projection maps from the tangent of the tangent bundle onto its vertical and horizontal subbundles. The projection onto the vertical bundle is defined intrinsically by Q V Γ = id V TQ ⊕ Γ and given locally by

↦

$$& \mathcal { Q } _ { \Gamma } ^ { \mathcal { V } } \colon T T Q \to \mathcal { V } T Q \\ & \left ( q, v _ { q }, V, W \right ) \mapsto \left ( q, v _ { q }, 0, W + \Gamma ( q, v _ { q } ) V \right ).$$

From the projection property, the projection onto the horizontal space H TQ is defined intrinsically as Q H Γ := id TTQ -Q V Γ and given locally by

$$& \mathcal { Q } _ { \Gamma } ^ { \mathcal { H } } \, \colon T T Q \to \mathcal { H } _ { ( q, v _ { q } ) } T Q \\ & \left ( q, v _ { q }, V, W \right ) \to ( q, v _ { q }, V, - \Gamma ( q, v _ { q } ) V ) \,.$$

From these projection maps, we define projection from TTQ to two copies of TQ , each of which is isomorphic to V TQ and H TQ , respectively. These projections are defined in the following diagram:

$$T T Q \xrightarrow { \mathcal { P } _ { \nu } } \nu T Q \xrightarrow { \mathcal { V } } T Q \\ T T Q \xrightarrow { \mathcal { Q } _ { \Gamma } ^ { \mathcal { H } } } \mathcal { H } T Q \xrightarrow { \mathcal { H } _ { \Gamma } } T Q$$

Therefore, the isomorphism TTQ ∼ = π ∗ TQ TQ ⊕ π ∗ TQ TQ is established by P H ⊕P V . Proposition 1. The projection maps P H : TTQ → TQ and P V : TTQ → TQ satisfy the following properties:

- · P H = ( π TQ ) ∗
- · P V ( X ∗ ) = ∇ v X , for a section X : Q → TQ .

Proof. See appendix A.1

## 2.2.2 Covariant Derivative of a Tangent Bundle Map

In this subsection, we find the covariant derivative of a general (nonlinear) tangent bundle map. This covariant derivative is essential for the results in subsequent sections. The authors have not found this result explicitly stated in the literature. Lastly, we provide a matrix formulation of the equations for practical computations.

Proposition 2. Let h : TQ → TQ be a general (not necessarily linear) vector bundle map over the tangent bundle π TQ : TQ → Q . Let v , w : Q → TQ be a sections. Then,

the covariant derivative of h ◦ w along v is

$$\nabla _ { v } ( h \circ w ) & = \nabla _ { v } ^ { \mathcal { H } } h ( w ) + D ^ { \mathcal { V } } h \left ( \nabla _ { v } w \right ), \\.. &.. &.. &... &.. &.. &..$$

$$D ^ { \mathcal { V } } h ( \nabla _ { v } w ) = \frac { \partial h ^ { m } } { \partial w ^ { l } } \left ( \nabla _ { v } w \right ) ^ { l } \frac { \partial } { \partial q ^ { m } },$$

where the vertical Jacobian D V h = P V D h ∣ ∣ V TQ is given locally by and the horizontal connection H ∇ v h w ( ) = P V D h ∣ ∣ H TQ is given locally by

$$\frac { \mathcal { H } } { \nabla _ { v } } h ( w ) = \left ( \frac { \partial h ^ { m } } { \partial q ^ { k } } v ^ { k } + \Gamma _ { k i } ^ { m } v ^ { k } h ^ { i } - \Gamma _ { k i } ^ { l } w ^ { k } v ^ { i } \frac { \partial h ^ { m } } { \partial w ^ { l } } \right ) \frac { \partial } { \partial q ^ { m } }.$$

$$P r o o f. \, \text{See appendix A.2}$$

Remark 2. Note that when the bundle map h : TQ → TQ linear in velocities, the vertical derivative D V h = h is linear in its argument, and the horizontal covariant derivative H ∇ X h reduces to the usual covariant derivative ∇ X h e.g., see equation (6) . Therefore, this "chain rule" is a generalization of the traditional chain rule for the covariant derivative of a linear vector bundle map. ,

Remark 3. For practical calculations, in the coordinate induced basis { ∂ ∂ q 1 , . . . , ∂ ∂ q n } the covariant derivative ∇ X Y given in (4) may be calculated in matrix form as

$$\nabla _ { X } Y = \left [ \frac { \partial Y } { \partial q } \right ] X + \left [ \Gamma ( q, X ) \right ] Y,$$

where [ ∂ Y ∂ q ] = [ ∂Y i ∂q k ] is the Jacobian of X and the Christoffel symbols are organized in a matrix [ Γ ( q X , )] = [Γ i kj X j ] . Similarly, the covariant derivative of a (1,1)-tensor A given in (6) may be calculated in matrix form as

$$[ \nabla _ { X } \mathcal { A } ] = \left [ \sum _ { k = 1 } ^ { n } \, \frac { \partial \mathcal { A } } { \partial q ^ { k } } X ^ { k } \right ] + [ \Gamma ( q, X ) ] [ \mathcal { A } ] - [ \mathcal { A } ] [ \Gamma ( q, X ) ].$$

Lastly, the covariant derivative H ∇ X h of a tangent bundle map h q w ( , ) as in (30) may be calculated in matrix form by

$$\overset { \mathcal { H } } { \nabla } _ { X } h = \left [ \frac { \partial h } { \partial q } ( q, w ) \right ] X + [ \Gamma ( q, X ) ] h - \left [ D ^ { \mathcal { V } } h \right ] [ \Gamma ( q, X ) ] w.$$

## 2.2.3 Dynamics of Mechanical Systems on TTQ

A vector field X : TQ → TTQ is a second order system on Q if and only if for each integral curve q : [ a, b ] → TQ of X we have d d t ( π TQ ◦ q ) ( t ) = ( π TQ ) ∗ ˙ q = q ( ) t , summarized in the following diagram [35]

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00095v1-figure-0236c829-c001-4f74-beef-a6db665b7d74.jpg)

If X = ( q v V , , , W ) is a second-order vector field, then it follows that the velocity components satisfy v i = V i .

Given, the decomposition TTQ ∼ H = TQ ⊕ V TQ , we construct the vector field associated with the affine connection system (8). Let ( q v , ) : [ a, b ] → TQ be a curve of satisfying (8). The forces acting on the system ( -d V /sharp G + F /sharp G ) in the right-hand side of (8) are vertically lifted to V ( q v , ) TQ by (22). The vector field associated with the parallel transport equation ∇ v ( ) t v ( ) = t 0 on the left-hand side of (8) is given by the horizontal lift (24) of the velocity v to H ( q v , q ) TQ . Therefore, the time variation of the curve d d t ( q v , ) ∈ T ( q v , ) TQ ∼ H = ( q v , ) TQ ⊕V ( q v , ) TQ satisfies the differential equations

$$\frac { d v } { d t } & = - \Gamma _ { i j } ^ { k } v ^ { i } v ^ { j } \frac { \partial } { \partial v ^ { k } } + \left ( - \frac { \partial V } { \partial q ^ { i } } + F _ { i } \right ) \mathbb { G } ^ { i k } \frac { \partial } { \partial v ^ { k } }, & ( 3 4 ) \\ \mathbb { d } q & = \partial &$$

$$\frac { \mathrm d q } { \mathrm d t } = v ^ { i } \frac { \partial } { \partial \tilde { q } ^ { i } }.$$

Remark 4. When there are no external forces in (34) , the vector field Z : TQ → TTQ

$$Z ( q, v ) = v ^ { i } \frac { \partial } { \partial \tilde { q } ^ { i } } - \Gamma ^ { k } _ { i j } v ^ { i } v ^ { j } \frac { \partial } { \partial v ^ { k } },$$

is called the geodesic spray [7, 10]. The integral curves of the geodesic spray Z in TQ project to geodesics with respect to ∇ in Q .

## 3 Affine Connection System as a Slow/Fast System with Strong Friction

Let Q be a smooth manifold and let D ⊂ TQ be a linear subbundle representing a set of m ideal nonholonomic constraints. To realize the nonholonomic constraints by friction forces, consider a Rayleigh dissipation function R : T Q q × T Q q → R , represented as a symmetric (0,2)-tensor on each tangent space [2]. This function induces a

dissipative friction force given by the (negative) fibre derivative - R F : TQ → T ∗ Q of the dissipation function. Assume that this friction force satisfies D = ker( F R ) to ensure that the nonholonomic distribution is attractive.

We construct this dissipation function R as follows. Consider a R m -valued oneform a q : T Q q → R m consisting of the m nonholonomic constraint one-forms a l ∈ T ∗ Q for l = 1 , . . . , m . For each configuration q ∈ Q , define the friction coefficients as the symmetric linear vector bundle map ¯( µ q ) : R m → R m ∗ ∼ = R m and finally, define the adjoint of the map representing the constraint as a ∗ q : R m ∗ → T ∗ q Q . Define the Rayleigh dissipation function R = 1 2 a ∗ q ◦ µ q ( ) ◦ a q and in local coordinates is given by

$$\mathcal { R } = \frac { 1 } { 2 } \underbrace { \mathbb { A } _ { i } ^ { l } ( q ) \bar { \mu } _ { l k } ( q ) \mathbb { A } _ { j } ^ { k } ( q ) } _ { \mathcal { R } _ { i j } ( q ) } d q ^ { j } \otimes d q ^ { i }. \\ \ e \, \text{Rayleigh friction force with large coefficients is defined as follows. Fix}$$

The dissipative Rayleigh friction force with large coefficients is defined as follows. Fix a small positive real number ε &gt; 0 (which is thought to be the time scale on which this force acts) and rescale the friction coefficients ¯ = 1 µ /εµ to find the friction force given by

$$F _ { \varepsilon } = - \frac { 1 } { \varepsilon } \mathbb { F R } ( v ) = - \frac { 1 } { \varepsilon } \mathcal { R } _ { i j } ( q ) v ^ { i } d q ^ { j } \in T ^ { * } Q.$$

Let q : [ a, b ] → Q be a smooth curve that satisfies the equations of motion formulated by the forced affine connection system

$$\nabla _ { v } v & = - \left ( d V \right ) ^ { \ddagger _ { G } } - \frac { 1 } { \varepsilon } \mathbb { F R } ( v ) ^ { \ddagger _ { G } },$$

$$\frac { \mathrm d } { \mathrm d t } q = v.$$

As this is a singularly perturbed system, it follows that (39-40 emits a normally hyperbolic invariant manifold (NHIM) for ε &gt; 0 [2]. That is, for small values of ε , there exists a perturbed manifold D ε that is C 1 -close and diffeomorphic to the nonholonomic constraint distribution D . The existence of this submanifold has been studied by Eldering and summarize the result [2] as follows. Write the system as a vector field (34) on TQ according to the decomposition TTQ ∼ H = TQ ⊕V TQ :

$$\frac { \text{d} } { \text{d} t } \left ( q, v \right ) = Z ( q, v ) + \left ( - d V ^ { \ddagger _ { G } } \right ) ^ { \uparrow \text{ver} } - \left ( \frac { 1 } { \varepsilon } \mathbb { F } \mathcal { R } ( v ) ^ { \ddagger _ { G } } \right ) ^ { \uparrow \text{ver} }$$

where we have vertically lifted the forces to TTQ using (22). Make the reparameterization of the trajectory to fast-time τ = t/ε to find

$$\frac { \text{d} } { \text{d} \tau } \left ( q, v \right ) = \varepsilon Z ( q, v ) + \varepsilon \left ( - d V ^ { \ddagger _ { G } } + U ^ { \ddagger _ { G } } \right ) ^ { \uparrow \text{ver} } - \left ( \mathbb { F R } ( v ) ^ { \ddagger _ { G } } \right ) ^ { \uparrow \text{ver} }. \quad \quad ( 4 2 )$$

When ε = 0 then the fast time system is of the form

$$\frac { d } { d \tau } \left ( q, v \right ) = - \left ( \mathbb { F R } ( v ) ^ { \natural _ { G } } \right ) ^ { \uparrow \prime }. \quad & ( 4 3 ) \\ \intertext { n force acts on an instantaneous time-scale ( \varepsilon = 0 ) } \text{it is the dom} \text{-}$$

When this friction force acts on an instantaneous time-scale ( ε = 0 ) it is the dominant force acting on the system. By construction, when we restrict the dynamics to the distribution D , the friction force vanishes; thus D is an invariant manifold for the fast-time dynamics (43). When the dynamics are restricted to the G -orthogonal complement D ⊥ , we find that the dynamics are attractive. Therefore, by the theory of normally hyperbolic invariant manifolds [21] there exists a perturbed manifold D ε that is C 1 -close and diffeomorphic to D .

## 3.1 Representation of the Slow Manifold

We represent the perturbed manifold D ε as the image of a section ˜ H ε : D → TQ over the bundle P : TQ → D . By definition, the section ˜ H ε satisfies the following commutative diagram

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00095v1-figure-7256c57f-7359-4b19-ae3e-8b5cdc88ad24.jpg)

That is, P ◦ ˜ H ε = id D . By the projection property id TQ = P ⊕ P ⊥ , we find that the section ˜ H ε is given by

$$\tilde { H } _ { \varepsilon } = \text{id} _ { \mathcal { D } } + \underbrace { \mathcal { P } ^ { \perp } \circ \tilde { H } _ { \varepsilon } } _ { = \colon h _ { \varepsilon } } ( \text{44} ) \\ \text{a section over } \pi _ { T Q } \colon T Q \to Q. \text{ Therefore, the slow manifold } \mathcal { D } _ { \varepsilon }$$

where h ε : D → D ⊥ is a section over π TQ : TQ → Q . Therefore, the slow manifold D ε is defined by the set

$$\mathcal { D } _ { \varepsilon } = \left \{ v \in T Q \, | \, v = v ^ { \mathcal { D } } + h _ { \varepsilon } ( v ^ { \mathcal { D } } ), \, v ^ { \mathcal { D } } \in \mathcal { D } \right \} \,, & & ( 4 5 ) \\ \mathfrak z \, & \text{the projection maps, we then define the sections } H _ { \varepsilon } = \tilde { H } _ { \varepsilon } \circ \mathcal { P } \, \colon T Q \to T Q$$

Using the projection maps, we then define the sections H ε = ˜ H ε ◦ P : TQ → TQ as

$$H _ { \varepsilon } ( v ) = \mathcal { P } ( v ) + h _ { \varepsilon } \circ \mathcal { P } ( v ),$$

where h ε ◦ P = P ⊥ ◦ ˜ H ε ◦ P : TQ → D ⊥ . Therefore, the perturbed manifold D ε = im ( H ε ) is given by

$$\mathcal { D } _ { \varepsilon } = \left \{ w \in T Q \, | \, w = H _ { \varepsilon } ( v ) = \mathcal { P } ( v ) + h _ { \varepsilon } \circ \mathcal { P } ( v ), \, v \in T Q \right \}. \quad \quad ( 4 7 )$$

Note that the form of this function h ε depends on the choice of projection map P . Remark 5. Physically speaking, the slow manifold represents the slip velocities in the orthogonal complement D ⊥ and are a function of the variables in the nonholonomic distribution D .

Lemma 1. The map H ε : TQ → TQ for the system (39) (40) -is a nonlinear projection onto D ε .

Proof. Given a section v : Q → TQ we find

$$H _ { \varepsilon } \circ H _ { \varepsilon } ( v ) & = H _ { \varepsilon } \left ( \mathcal { P } ( v ) + h _ { \varepsilon } \circ \mathcal { P } ( v ) \right ) \\ & = \left ( \tilde { H } _ { \varepsilon } \circ \mathcal { P } \right ) \left ( \mathcal { P } ( v ) + h _ { \varepsilon } \circ \mathcal { P } ( v ) \right ) \\ & = \tilde { H } _ { \varepsilon } \left ( \mathcal { P } \circ \mathcal { P } ( v ) + \mathcal { P } \circ h _ { \varepsilon } \circ \mathcal { P } ( v ) \right ) \\ & = \tilde { H } _ { \varepsilon } \left ( \mathcal { P } \circ \mathcal { P } ( v ) + \mathcal { P } \circ \mathcal { P } ^ { \perp } \circ \tilde { H } _ { \varepsilon } \circ \mathcal { P } ( v ) \right ) \\ & = \tilde { H } _ { \varepsilon } \circ \mathcal { P } ( v ) = H _ { \varepsilon } ( v ).$$

Remark 6. Since H ε defines a nonlinear projection map, we may also define its complement H ⊥ ε := id Q -H ε . This defines a projection map onto the complement of the slow manifold abbreviated by TQ/ D ε .

Lemma 2. The map H ⊥ ε : TQ → TQ/ D ε is a nonlinear projection map that is complementary to H ε .

Proof. First, we show that H ε is a projection map. Let v : Q → TQ be a section. Then,

$$H _ { \varepsilon } ^ { \perp } \circ H _ { \varepsilon } ^ { \perp } ( v ) & = ( i d _ { T Q } - H _ { \varepsilon } ) \circ ( v - H _ { \varepsilon } ( v ) ) \\ & = v - H _ { \varepsilon } ( v ) - H _ { \varepsilon } ( v ) + H _ { \varepsilon } \circ H _ { \varepsilon } ( v ) \\ & = v - H _ { \varepsilon } ( v ) = H _ { \varepsilon } ^ { \perp } ( v ).$$

Lastly, we show that this map is complement to H ε .

$$H _ { \varepsilon } \circ H _ { \varepsilon } ^ { \perp } ( v ) & = H _ { \varepsilon } ( v ) - H _ { \varepsilon } \circ H _ { \varepsilon } ( v ) = 0, \\ H _ { \varepsilon } ^ { \perp } \circ H _ { \varepsilon } ( v ) & = H _ { \varepsilon } \circ H _ { \varepsilon } ( v ) - H _ { \varepsilon } ( v ) = 0.$$

Remark 7. As a consequence of the complementary projection map in Proposition 2, we find that v ∈ D ε if and only if

$$H _ { \varepsilon } ^ { \perp } ( v ) & = v - H _ { \varepsilon } ( v ) = v - \mathcal { P } ( v ) - h _ { \varepsilon } \circ \mathcal { P } ( v ) \\ & = \mathcal { P } ^ { \perp } ( v ) - h _ { \varepsilon } \circ \mathcal { P } ( v ) = 0.$$

This is summarized in the following proposition.

Lemma 3. v ∈ D ε if and only if P ⊥ ( v ) = h ε ◦ P ( v ) .

Proof. Suppose that w ∈ D ε . Then, there exists v ∈ TQ such that w = P ( v ) + h ε ( v ) . We calculate

$$\mathcal { P } ^ { \perp } ( w ) & = \mathcal { P } ^ { \perp } \left ( \mathcal { P } ( v ) + \mathcal { P } ^ { \perp } \circ h _ { \varepsilon } ( v ) \right ) \\ & = \mathcal { P } ^ { \perp } \circ \mathcal { P } ^ { \perp } \circ H _ { \varepsilon } ( v ) = h _ { \varepsilon } ( v ) = \mathcal { P } ^ { \perp } \circ \tilde { H } _ { \varepsilon } \circ \mathcal { P } ( v ).$$

Moreover, we find that w and v project to the same vector

$$\mathcal { P } ( w ) = \mathcal { P } \left ( \mathcal { P } ( v ) + h _ { \varepsilon } \circ \mathcal { P } ( v ) \right ) \right ) = \mathcal { P } ( v ).$$

Therefore, P ⊥ ( w ) = P ⊥ ◦ ˜ H ε ◦ P ( v ) = h ε ◦ P ( w ) .

Conversely, suppose v ∈ TQ such that P ⊥ ( v ) = h ε ◦ P ( v ) . Then, since P is a projection map, we have that id TQ = P ⊕ P ⊥ . Therefore, since we assume that P ⊥ ( v ) = h ε ◦ P ( v ) , we find

$$v = \mathcal { P } ( v ) + \mathcal { P } ^ { \perp } ( v ) = \mathcal { P } ( v ) + h _ { \varepsilon } \circ \mathcal { P } ( v ).$$

Therefore, v ∈ D ε .

By Proposition 3, we find that the distribution can be represented as the set

$$\mathcal { D } _ { \varepsilon } = \left \{ \upsilon \in T Q \, | \, \mathcal { P } ^ { \perp } ( \upsilon ) = h _ { \varepsilon } \circ \mathcal { P } ( \upsilon ) \right \}. \quad & ( 4 9 ) \\ \intertext { 3. \, The equations of motion for a curve } \, \mathbf q \colon [ a, b ] \to Q \, \text{lying on the slow}$$

Proposition 3. The equations of motion for a curve q : [ a, b ] → Q lying on the slow manifold D ε are given by

$$\nabla _ { v ( t ) } v ( t ) = - d V ^ { \ddagger _ { G } } - \frac { 1 } { \varepsilon } \mathbb { F } \mathcal { R } \left ( v ( t ) \right ) ^ { \ddagger _ { G } }$$

$$\frac { d } { d t } q ( t ) = v ( t ), \quad \frac { \lambda } { d t } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 2 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac { 1 } { 3 } \frac$$

$$\mathcal { P } ^ { \perp } ( v ( t ) ) = h _ { \varepsilon } \circ \mathcal { P } ( v ( t ) ) \iff v ( t ) \in \mathcal { D } _ { \varepsilon }.$$

Proof. The equations of motion (50)-(51) are described by a forced affine connection system. The existence of the slow manifold follows from the geometric singular perturbation theory of normally hyperbolic manifolds [21]. By Proposition 3 the slow manifold is represented by (52).

## 4 Approximation of the Slow Subsystem

In this section, we find a generating equation for the invariant manifold D ε using the equations of motion (50) - (52). The strategy is to consider velocities restricted to the slow manifold D ε and expand the section h ε as a power series and match terms to the desired order in /epsilon1 .

## 4.1 Generating equation for the slow-manifold

Let us find an expression for the dissipation function ( F R ) /sharp G in terms of the unknown function h ε describing the perturbation of the ideal constraints. Using the projection maps, we project the dynamical equations (50) onto D⊕D ⊥ ∼ = TQ

$$& \mathcal { P } \left ( \nabla _ { v } v \right ) = - \mathcal { P } \left ( d V ^ { \ddagger _ { G } } \right ) \\ & \mathcal { P } ^ { \perp } \left ( \nabla _ { v } v \right ) = - \mathcal { P } ^ { \perp } \left ( d V ^ { \ddagger _ { G } } \right ) - \frac { 1 } { c } \mathbb { F } \mathcal { R } ^ { \ddagger _ { G } } ( v ).$$

$$\varphi ( \nu \nu ) & = - \varphi ( \nu \nu ) \quad \nu \nu \quad \nu \nu \quad \nu \nu \quad \\ \varphi ^ { \perp } \left ( \nabla _ { v } v \right ) & = - \varphi ^ { \perp } \left ( d V ^ { \ddagger _ { G } } \right ) - \frac { 1 } { \varepsilon } \mathbb { F } \mathcal { R } ^ { \ddagger _ { G } } ( v ). & ( 5 4 ) \\ \varphi _ { \text{audio} } \text{on the non-automation} \mathbb { F } \mathbb { D } ^ { \ddagger _ { G } } \cdot T \, \Delta \quad T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad } \cdot \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\nu \quad \text{ } \cdot \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\matrix } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \ \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb { } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \ :: \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\chbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \alpha \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\tau \quad \cdot \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T\, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text$\mathabb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta\quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbin } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb{ } T \, \Delta \quad \text{$\mathbb
 } \cdot$$

Lemma 4. Consider the transformation F R /sharp G q : T Q q → T Q q , then there exists a transformation Q q : D ⊥ q →D ⊥ q satisfying Q q ◦ F R /sharp G q = P ⊥ q .

Proof. The friction force is a linear transformation F R /sharp G q : T Q q → T Q q satisfying D = ker( F R /sharp G q ) and im ( F R /sharp G q ) = D ⊥ . Thus, by the first isomorphism theorem [39], we have the result.

Remark 8. The map Q q may be constructed as follows. Let Φ( q ) be a frame spanning D ⊕D q ⊥ q . Define linear maps in this frame with the subscript Φ , so, [ Q Φ ] = Φ -1 [ Q ]Φ , [ P ⊥ Φ ] = Φ -1 [ P ⊥ ]Φ and [ F R /sharp G Φ ] = Φ -1 [ F R /sharp G ]Φ and the equation [ P ⊥ Φ ] = [ Q Φ ][ F R /sharp G Φ ] is given by the block matrix equation of the form

$$\begin{array} { c } \cdot & \cdot & \cdot & \cdot \\ \underbrace { \left [ \mathbb { O } \ \mathbb { O } \ \mathbb { I } \right ] } _ { [ \mathbb { P } _ { q } ^ { + } ] } = \underbrace { \left [ \mathbb { O } \ \mathbb { O } \ \mathbb { I } \right ] } _ { [ \mathbb { Q } _ { + } ] } \left [ \mathbb { O } \ \mathbb { O } \ \mathbb { I } \right ] } _ { [ \mathbb { P } \mathbb { R } _ { \Phi } ^ { \pm G } ] }. \\ \cdot & \cdot & \cdot & \end{array}$$

Therefore, we set ˜ Q = ˜ F R -1 and perform the inverse transformation [ Q q ] = Φ[ Q Φ ]Φ -1 to obtain the result.

Theorem 1 (Generating Equation for the Slow Manifold) . Let h ε : D → TQ be the section defining the invariant manifold D ε with respect to the dynamics (50) (52) -. Let h ε = ∑ k ≥ 0 ε k h ( k ) be an expansion of the section by a formal power series. Then, the slow manifold may be recursively approximated from the following generating equation

$$\Lambda _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambla } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lamb { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \hat { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \frac { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \matrix { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } ^ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambd } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambta } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lamb. } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lamb� } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lamb! } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \LambDA } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambrum } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _ { \vec { \Lambda } _$$

Proof. By analogy with the ideal nonholonomic case [10], take the covariant derivative of the latter invariant manifold equation (52)

$$\nabla _ { v } ( \mathcal { P } ^ { \perp } ( v ) ) = \nabla _ { v } ( h _ { \varepsilon } \circ \mathcal { P } ( v ) ) & & ( 5 7 ) \\ \text{affine connection system (50)-(51). The left-hand side is given by}$$

and substitute the affine connection system (50)-(51). The left-hand side is given by

$$\dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots\delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \dots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta \cdots \delta
\nabla _ { v } \left ( \mathcal { P } ^ { \perp } ( v ) \right ) & = \left ( \nabla _ { v } \mathcal { P } ^ { \perp } \right ) ( v ) + \mathcal { P } ^ { \perp } ( \nabla _ { v } v ) \right ) \\ & = \left ( \nabla _ { v } \mathcal { P } ^ { \perp } \right ) ( v ) - \frac { 1 } { \varepsilon } \mathbb { F } \mathcal { R } ( v ) ^ { \dagger _ { G } } - \mathcal { P } ^ { \perp } \left ( d V ^ { \dagger _ { G } } \right ) \cdot \\ \intertext { e other hand, according to Proposition 2 and the chain rule, we find the 1 }$$

On the other hand, according to Proposition 2 and the chain rule, we find the righthand side

$$\text{ sure} \\ \nabla _ { v } ( h _ { \varepsilon } \circ \mathcal { P } ( v ) ) = D ^ { \mathcal { V } } h _ { \varepsilon } \left ( \nabla _ { v } ( \mathcal { P } ( v ) ) \right ) + \frac { \mathcal { H } } { \nabla _ { v } } h _ { \varepsilon } \left ( \mathcal { P } ( v ) \right ) \\ = D ^ { \mathcal { V } } h _ { \varepsilon } \left ( \left ( \nabla _ { v } \mathcal { P } \right ) ( v ) + \mathcal { P } \left ( \nabla _ { v } v \right ) \right ) + \frac { \mathcal { H } } { \nabla _ { v } } h _ { \varepsilon } \left ( \mathcal { P } ( v ) \right ) \\ = D ^ { \mathcal { V } } h _ { \varepsilon } \left ( ( \nabla _ { v } \mathcal { P } ) ( v ) \right ) - D ^ { \mathcal { V } } h _ { \varepsilon } \left ( \mathcal { P } \left ( d V ^ { \mathfrak { t } _ { G } } \right ) \right ) + \frac { \mathcal { H } } { \nabla _ { v } } h _ { \varepsilon } ( \mathcal { P } ( v ) ). \\ \text{ing both sides of this calculation together, we find the following identity fo.}$$

Putting both sides of this calculation together, we find the following identity for the dissipative forces

$$& - \frac { 1 } { \varepsilon } \mathbb { F } \mathcal { R } ^ { \ddagger _ { G } } ( v ) = - ( \nabla _ { v } \mathcal { P } ^ { \perp } ) ( v ) + \mathcal { P } ^ { \perp } ( d V ^ { \ddagger _ { G } } ) & & ( 5 8 ) \\ & + D ^ { \mathcal { V } } h _ { \varepsilon } \left ( ( \nabla _ { v } \mathcal { P } ) ( v ) \right ) + \frac { \mathcal { H } } { \nabla _ { v } } h _ { \varepsilon } ( \mathcal { P } ( v ) ) - D ^ { \mathcal { V } } h _ { \varepsilon } \left ( \mathcal { P } ( d V ^ { \ddagger _ { G } } ) \right ) \\ & \text{Lastly, we restrict (58) to the velocities $v \in \mathcal{D}_{\varepsilon}$ to obtain the generating equation}$$

Lastly, we restrict (58) to the velocities v ∈ D ε to obtain the generating equation for h ε

$$\tau h _ { \varepsilon } \\ \mathbb { F R } ^ { \ddagger _ { G } } h _ { \varepsilon } ( v ^ { \mathcal { D } } ) & = \varepsilon \left ( \nabla _ { v ^ { \mathcal { D } } + h _ { \varepsilon } } \mathcal { P } ^ { \perp } \right ) ( v ^ { \mathcal { D } } + h _ { \varepsilon } ) + \varepsilon \mathcal { P } ^ { \perp } ( d V ^ { \ddagger _ { G } } ) \\ & - \varepsilon D ^ { \mathcal { V } } h _ { \varepsilon } ( - d V ^ { \ddagger _ { G } } ) - \varepsilon \left ( \frac { \mathcal { H } } { \nabla } _ { v ^ { \mathcal { D } } + h _ { \varepsilon } } h _ { \varepsilon } \right ) ( v ^ { \mathcal { D } } ) \\ & + \varepsilon D ^ { \mathcal { V } } h _ { \varepsilon } \left ( ( \nabla _ { v ^ { \mathcal { D } } + h _ { \varepsilon } } \mathcal { P } ) ( v ^ { \mathcal { D } } + h _ { \varepsilon } ) \right ). \\ \intertext { I Lemma 4 and linearity we obtain the generating equation (56). } \quad & \Box$$

By Lemma 4 and linearity we obtain the generating equation (56).

Remark 9. Note that there is an ε in all terms on the right-hand side of (56) . Thus, any power series expansion of h ε results in an order higher on the right-hand side than the left-hand side. Therefore, the generating equation (56) provides a recurrence relation for the approximation of h ε .

Theorem 2 (Equivalence of Reaction Forces) . The generating equation for h ε ( v D ) in (56) results in the same reaction forces -1 ε F R /sharp G h ε ( v D ) as those reaction forces resulting from the classical invariance condition

$$\frac { \text{d} } { \text{d} t } \mathcal { P } ^ { \perp } ( v ) \Big | _ { \mathcal { D } _ { \epsilon } } = \frac { \text{d} } { \text{d} t } h _ { \epsilon } \circ \mathcal { P } ( v ) \Big | _ { \mathcal { D } _ { \epsilon } }. & & ( 6 0 ) \\ \text{loss of generality, assume that the potential $V=0$. The equations of}$$

Proof. Without loss of generality, assume that the potential V = 0 . The equations of motion (50)-(52) can be written as

$$\dot { v } = - \Gamma ( q, v ) v - \frac { 1 } { \varepsilon } \mathbb { F R } ^ { \ddagger _ { G } } ( v ) \\. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \. \quad \. \. \quad \. \. \quad \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \ \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \ }$$

$$\dot { q } = v \in \mathcal { D } _ { \varepsilon } \iff \mathcal { P } ^ { \perp } ( v ) = h _ { \varepsilon } \circ \mathcal { P } ( v ).$$

Taking a Lie derivative of (52) along the along these dynamics we find

$$\text{a Lie derivative of (52) along the along these dynamics we find} \\ \left [ \sum _ { k = 1 } ^ { n } \frac { \partial \mathcal { P } ^ { \perp } } { \partial q ^ { k } } v ^ { k } \right ] v - [ \mathcal { P } ^ { \perp } ] [ \Gamma ( q, v ) ] v - \frac { 1 } { \varepsilon } [ \mathcal { P } ^ { \perp } ] [ \mathbb { F } \mathcal { R } ^ { \ddagger _ { G } } ] v = \\ \left [ \frac { \partial h _ { \varepsilon } } { \partial q } ( q, \mathcal { P } ( v ) ) \right ] v + [ D ^ { \nu } h _ { \varepsilon } ] \left [ \sum _ { k = 1 } ^ { n } \frac { \partial \mathcal { P } } { \partial q ^ { k } } v ^ { k } \right ] v \\ - [ D ^ { \nu } h _ { \varepsilon } ] [ \mathcal { P } ] [ \Gamma ( q, v ) ] v - \frac { 1 } { \varepsilon } [ D ^ { \nu } h _ { \varepsilon } ] [ \mathcal { P } ] [ \mathbb { F } \mathcal { R } ^ { \ddagger _ { G } } ] v. \\ \text{he quantities in the square brackets are the cannot matrix notation intr}$$

where the quantities in the square brackets are the compact matrix notation introduced in Remark 3. As D = ker F R /sharp G , the last term vanishes. Let v = v D ε ∈ D ε be a velocity in the slow manifold, then the reaction forces are restricted to D ε are given by

$$- \, \frac { 1 } { \varepsilon } \mathbb { F R } ^ { \ddagger _ { c } } ( v ^ { \mathcal { D } _ { \epsilon } } ) & = - \left [ \sum _ { k = 1 } ^ { n } \frac { \partial \mathcal { P } ^ { \perp } } { \partial q ^ { k } } v ^ { k } _ { \mathcal { D } _ { \epsilon } } \right ] v ^ { \mathcal { D } _ { \epsilon } } & + [ \mathcal { P } ^ { \perp } ] [ \Gamma ( q, v ^ { \mathcal { D } _ { \epsilon } } ) ] v ^ { \mathcal { D } _ { \epsilon } } \\ & \quad + \left [ \frac { \partial h _ { \epsilon } } { \partial q } ( q, v ^ { \mathcal { D } } ) \right ] v ^ { \mathcal { D } _ { \epsilon } } + [ D ^ { \nu } h _ { \epsilon } ] \left [ \sum _ { k = 1 } ^ { n } \frac { \partial \mathcal { P } } { \partial q ^ { k } } v ^ { k } _ { \mathcal { D } _ { \epsilon } } \right ] v ^ { \mathcal { D } _ { \epsilon } } - [ D ^ { \nu } h _ { \epsilon } ] [ \mathcal { P } ] [ \Gamma ( q, v ^ { \mathcal { D } _ { \epsilon } } ) ] v ^ { \mathcal { D } _ { \epsilon } }. \\ \zeta \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \tau \in \, \dots \,..$$

Contrastingly, the reaction forces h ε resulting from the generating equation (56) are calculated from the equation

$$- \frac { 1 } { \varepsilon } \mathbb { F } \mathcal { R } ^ { \ddagger _ { G } } ( v ) = - ( \nabla _ { \dot { q } } \mathcal { P } ^ { \perp } ) ( v ) + D ^ { \mathcal { V } } h _ { \varepsilon } ( \nabla _ { \dot { q } } \mathcal { P } ) ( v ) + \overset { \mathcal { H } } { \nabla _ { \dot { q } } } h _ { \varepsilon } ( \mathcal { P } ( v ) ).$$

as was observed in equation (58). By Remark 3 the covariant derivatives in this equation may be expressed in matrix form by

$$- \frac { 1 } { \varepsilon } \mathbb { F R } ^ { \ddagger _ { G } } ( v ) = - \left [ \sum _ { k = 1 } ^ { n } \frac { \partial \mathcal { P } ^ { \perp } } { \partial q ^ { k } } v ^ { k } \right ] v - [ \Gamma ( q, v ) ] [ \mathcal { P } ^ { \perp } ] v + [ \mathcal { P } ^ { \perp } ] [ \Gamma ( q, v ) ] v$$

$$& + [ D ^ { \Psi } h _ { \varepsilon } ] \left [ \sum _ { k = 1 } ^ { n } \frac { \partial \mathcal { P } } { \partial q ^ { k } } v ^ { k } \right ] v + D ^ { \Psi } h _ { \varepsilon } [ \Gamma ( q, v ) ] [ \mathcal { P } ] v - [ D ^ { \Psi } h _ { \varepsilon } ] [ \mathcal { P } ] [ \Gamma ( q, v ) ] v \\ & + \left [ \frac { \partial h _ { \varepsilon } } { \partial q } ( q, v ^ { \mathcal { D } } ) \right ] v + [ \Gamma ( q, v ) ] h _ { \varepsilon } ( q, v ^ { \mathcal { D } } ) - [ D ^ { \Psi } h _ { \varepsilon } ] [ \Gamma ( q, v ) ] [ \mathcal { P } ] v. \\ - \quad \cdot \quad \cdot \quad \wedge \quad \cdot \quad - \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot$$

Then, restricting to v = v D ε , and using P ⊥ ◦ h ε = h ε terms cancel and therefore, these friction forces (64) restricted to D ε are identical to (63).

Corollary 2.1 (Equivalence of Tangent Conditions) . The dynamics defined by (50) -(52) restricted to D ε are tangent to the manifold D ε .

Proof. In local coordinates (50)-(52) are equivalent to the equations of motion on TTQ given by (41). Following singular perturbation theory, the classical invariance equation (60) ensures that the reaction forces result in dynamics tangent to D ε . By Theorem 2, the reaction friction forces described by (56) are equivalent to the reaction forces resulting from the classical invariance equation. Therefore, the dynamics (50)-(52) along with the generating equation (56) result in dynamics tangent to D ε .

Corollary 2.2 (Equivalence of Invariance Conditions) . The classical invariance condition (60) for the dynamics (50) (52) -is equivalent to the condition

$$\nabla _ { \dot { q } } ( \mathcal { P } ^ { \perp } ( v ) ) \Big | _ { \mathcal { D } _ { \epsilon } } = \nabla _ { \dot { q } } ( h _ { \epsilon } \circ \mathcal { P } ( v ) ) \Big | _ { \mathcal { D } _ { \epsilon } }. & & ( 6 5 ) \\ \text{orem 2 and Corollary 2.1 both conditions result in the same dynamical} \\ \text{tion reduction to the leap manifold } \mathcal { U } & & \Box$$

Proof. By Theorem 2 and Corollary 2.1 both conditions result in the same dynamical equations of motion restricted to the slow manifold D ε .

Remark 10. The later covariant derivative invariance condition (65) may be geometrically interpreted as the parallel transport of the condition P ⊥ ( v ) = h ε ◦ P ( v ) along trajectories ˙ q = v ∈ D ε lying in the slow manifold.

## 4.2 Approximation of Slip Velocities

In this section, we employ the generating equation (56) from Theorem 1 to compute the zeroth, first and second order approximations of the slow manifold h ε . Expand h ε as a formal power series

$$h _ { \varepsilon } = h ^ { ( 0 ) } + \varepsilon h ^ { ( 1 ) } + \varepsilon ^ { 2 } h ^ { ( 2 ) } + O ( \varepsilon ^ { 3 } ).$$

Further expanding the power series allows for higher order approximations of D ε .

## 4.2.1 Zeroth Order Approximation of Slip Velocities

Lets consider the zeroth order term of the slow manifold h ε = h (0) + O ε ( ) . Observe that when ε = 0 we have

$$\mathcal { D } _ { \varepsilon } \Big | _ { \varepsilon = 0 } = \mathcal { D } _ { 0 } = \{ v \in T Q \, | \, \mathcal { P } ^ { \perp } ( v ) = h ^ { ( 0 ) } \}.$$

Since h ε is a perturbation of the nonholonomic distribution, we find that h (0) = 0 . Then, to find the zeroth order velocity, solve the identity P ⊥ ( v ) = 0 for v = v D ∈ D . Remark 11. The zeroth order approximation of the slow fast velocities are equivalent to the nonholonomic velocities.

## 4.2.2 First Order Approximation of Slip Velocities

Let us assume that section h ε is taken to first order, so we expand h ε = ε h (1) + O ε ( 2 ) . Then according to the generating equation (56) to this first order slip velocity, we find

$$h ^ { ( 1 ) } = \mathbb { Q } _ { q } \circ \left [ ( \nabla _ { v ^ { \mathcal { D } } } \mathcal { P } ^ { \perp } ) ( v ^ { \mathcal { D } } ) - \mathcal { P } ^ { \perp } ( d V ^ { \mathfrak { k } _ { G } } ) \right ].$$

Remark 12. The first order approximation of the slow fast velocities are proportional to the nonholonomic Lagrange multipliers (15) . This is consistent with previous results [21].

## 4.2.3 Second Order Approximation of Slip Velocities

Lastly, assume that section h ε is taken to second order, so we expand the formal power series as h ε = ε h (1) + ε 2 h (2) + O ε ( 3 ) . Then according to the generating equation (56) to this first order slip velocity, we find

$$h ^ { ( 2 ) } & = \mathbb { Q } _ { q } \circ \left [ ( \nabla _ { h ^ { ( 1 ) } } \mathcal { P } ^ { \perp } ) ( v ^ { \mathcal { D } } ) + ( \nabla _ { v ^ { \mathcal { D } } } \mathcal { P } ^ { \perp } ) ( h ^ { ( 1 ) } ) \right ] \\ & \quad - \mathbb { Q } _ { q } \circ \left [ D ^ { \mathcal { V } } h ^ { ( 1 ) } \left ( ( \nabla _ { v ^ { \mathcal { D } } } \mathcal { P } ) ( v ^ { \mathcal { D } } ) \right ) - D ^ { \mathcal { V } } h ^ { ( 1 ) } \mathcal { P } ( d V ^ { \mathfrak { T } _ { G } } ) \right ] - \mathbb { Q } _ { q } \circ \left [ \frac { \mathcal { H } } { \nabla _ { v ^ { \mathcal { D } } } } h ^ { ( 1 ) } ( v ^ { \mathcal { D } } ) \right ]. \\ \in & \quad \cdot \, \quad \cdot \, \quad \cdot \, \quad \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \,.$$

Remark 13. Higher order approximations of h ε may be calculated recursively in terms of higher order covariant derivatives ∇ n . However, in practice, first and second order approximations suffice. This is due to the fact that the friction forces -1 ε F R /sharp G ( v ) in (50) restricted to the second order approximation results in a first order force -ε F R /sharp G ( h (2) ) that is a function of zeroth and first-order velocities.

## 4.3 Approximation of Slip Dynamics

In this section, we use the approximations of the slip velocities h ε = ε h (1) + ε h (2) + O ε ( ) to approximate the zeroth and first order dynamics. Due to the recurrence relation defined by the generating equation (56) on h ε we observe that the zeroth order approximation of the dynamics (50)-(52) depends on the first order approximation h (1) . Furthermore, the first order approximation of the dynamics (50)-(52) depends on the second order approximation h (2) , and so on. Therefore, the order of the slip velocities h ε is always one order above the order of the dynamics.

## 4.3.1 Zeroth Order Dynamics

By the equations of motion (50) and the first order approximation of the slow manifold h ε = ε h (1) + O ε ( 2 ) , we find the zeroth order equations of motion are given by

$$\frac { d v } { d t } = - \Gamma ( q, v ^ { \mathcal { D } } ) v ^ { \mathcal { D } } - ( \nabla _ { v ^ { \mathcal { D } } } \mathcal { P } ^ { \perp } ) ( v ^ { \mathcal { D } } ) - \mathcal { P } ( d V ^ { \dagger _ { G } } )$$

$$\frac { \dddot { \vec { d } q } } { \vec { d } t } = v ^ { \mathcal { D } } \in \mathcal { D } _ { 0 } \cong \mathcal { D }.$$

Remark 14. The zeroth order equations of motion (69) are exactly the nonholonomic equations of motion (12) (14) -. This result is consistent with previous realizations of slow/fast nonholonomic systems [2].

## 4.3.2 First Order Dynamics

Using the generating function (56) we obtain the first and second order approximations h ε = ε h (1) + ε 2 h (2) + O ε ( 3 ) to obtain (67) and (68), the former gives the slip velocity and the later gives the reaction force. Restricting the dynamics to the first order approximation of slip, we find the dynamics:

$$\frac { d v } { d t } & = - \Gamma ( q, v ^ { \mathcal { D } } ) v ^ { \mathcal { D } } - ( \nabla _ { v ^ { \mathcal { D } } } \mathcal { P } ^ { \perp } ) ( v ^ { \mathcal { D } } ) - \mathcal { P } ( d V ^ { \mathfrak { g } } ) \\ & \quad - \varepsilon \left [ \Gamma ( q, v ^ { \mathcal { D } } ) h ^ { ( 1 ) } + \Gamma ( q, h ^ { ( 1 ) } ) v ^ { \mathcal { D } } \right ] + \varepsilon \frac { \mathcal { H } } { \nabla _ { v ^ { \mathcal { D } } } } h ^ { ( 1 ) } ( v ^ { \mathcal { D } } ) \\ & \quad - \varepsilon \left [ ( \nabla _ { v ^ { \mathcal { D } } } \mathcal { P } ^ { \perp } ) ( h ^ { ( 1 ) } ) + ( \nabla _ { h ^ { ( 1 ) } } \mathcal { P } ^ { \perp } ) ( v ^ { \mathcal { D } } ) \right ] \\ & \quad + \varepsilon \, D ^ { \mathcal { V } } h ^ { ( 1 ) } \left [ \left ( ( \nabla _ { v ^ { \mathcal { D } } } \mathcal { P } ) ( v ^ { \mathcal { D } } ) \right ) - \mathcal { P } ( d V ^ { \mathfrak { g } } ) \right ] \\ \frac { d q } { d t } & = v ^ { \mathcal { D } } + \varepsilon h ^ { ( 1 ) } ( v ^ { \mathcal { D } } ) \\ \text{re the first order approximation of the slip } h ^ { ( 1 ) } \text{ is given by } ( 6 7 ), \text{ and vertical}$$

where the first order approximation of the slip h (1) is given by (67), and vertical derivative D V h (1) and horizontal H ∇ v D h (1) is calculated according to equation (29) and (30).

## 5 Case Study: Vertical Rolling Disk

To illustrate the theory, we consider a vertical rolling disk for simplicity. However, our approach may be applied to other systems with (linear) nonholonomic rolling constraints such as wheeled mobile robots [40] and vehicle systems [41]. For the traditional derivation of the nonholonomic equations for the vertical rolling disk using the Lagrange-d'Alembert variational principle see [1, Sec. 1.4] and for the affine-connection approach see [7, Sec. 4.5], and [10].

The configuration space for a planar vertical rolling disk is Q = SE (2) × SO (2) with local coordinates q = ( θ, x, y, ϕ ) and local coordinates for velocities v = ( v , v θ x , v y , v ϕ ) .

The Riemannian metric defining the kinetic energy is given by

$$\mathbb { G } _ { q } = \text{diag} \left ( I, m, m, J \right )$$

where m is the mass of the disk, I is the moment of inertia normal to the plane and J is the moment of inertia about the rolling axis. Note that as G q is constant, the corresponding Christoffel symbols Γ ( q v , ) = 0 for all ( q v , ) ∈ TQ i.e. , the metric G q is flat.

The nonholonomic constraint one-forms

$$a _ { q } ^ { 1 } = d x - R \cos \theta d \varphi, \ \ a _ { q } ^ { 2 } = d y - R \sin \theta d \varphi,$$

enforcing rolling without slipping are given by

$$a _ { q } ^ { 1 } ( v _ { \theta }, v _ { x }, v _ { y }, v _ { \varphi } ) = v _ { x } - R \cos \theta v _ { \varphi } = 0,$$

$$a _ { q } ^ { 2 } ( v _ { \theta }, v _ { x }, v _ { y }, v _ { \varphi } ) = v _ { y } - R \sin \theta v _ { \varphi } = 0,$$

where R is the radius of the disk. This constraint defines a nonholonomic distribution D q = ker( a 1 q ) ∩ ker( a 2 q ) which are spanned by the vector fields

$$\frac { \partial } { \partial \theta }, \ \ R \cos \theta \frac { \partial } { \partial x } + R \sin \theta \frac { \partial } { \partial y } + \frac { \partial } { \partial \varphi }.$$

Using the musical isomorphisms, we find that the G -orthogonal complement spanned by the orthogonal complement D ⊥ = ( D ◦ ) /flat G is spanned by the vector fields

$$\frac { 1 } { m } \frac { \partial } { \partial x } - \frac { R } { J } \cos \theta \frac { \partial } { \partial \varphi }, \quad \frac { 1 } { m } \frac { \partial } { \partial y } - \frac { R } { J } \sin \theta \frac { \partial } { \partial \varphi }.$$

The projection map P : TQ →D onto the distribution D is constructed by the vector fields spanned by (76) and the metric (72) and is given by

$$\begin{array} { c c c c } \alpha \, \L y \, ( \theta ) \, \text{and} \, \text{the main} \, ( \tau ) \, \text{and} \, \text{is given} \, \L y \end{array} \\ \mathcal { P } = \mathcal { I } \begin{bmatrix} \frac { 1 } { \mathcal { I } } & 0 & 0 & 0 \\ 0 & \cos ^ { 2 } \theta & \frac { 1 } { 2 } \sin 2 \theta & \frac { J } { m R } \cos \theta \\ 0 & \frac { 1 } { 2 } \sin 2 \theta & \sin ^ { 2 } \theta & \frac { J } { m R } \sin \theta \\ 0 & \frac { \cos \theta } { R } & \frac { \sin \theta } { R } & \frac { J } { m R ^ { 2 } } \end{bmatrix} \\ \frac { m R ^ { 2 } } { m R } \, \text{is a ratio of inertia.} \end{array}$$

where I := mR 2 J + mR 2 is a ratio of inertia.

Similarly, the projection map P ⊥ : TQ → D ⊥ onto the distribution D ⊥ is constructed by the vector fields spanned by (76) and the metric (72) and is given by

$$\mathcal { P } ^ { \perp } & = \mathcal { I } \begin{bmatrix} 0 & 0 & 0 \\ 0 & \frac { J + m R ^ { 2 } \sin ^ { 2 } \theta } { m R ^ { 2 } } & - \frac { 1 } { 2 } \sin 2 \theta & - \frac { J } { m R } \cos \theta \\ 0 & - \frac { 1 } { 2 } \sin 2 \theta & \frac { J + m R ^ { 2 } \cos ^ { 2 } \theta } { m R ^ { 2 } } & - \frac { J } { m R } \sin \theta \\ 0 & - \frac { \cos \theta } { R } & - \frac { \sin \theta } { R } & 1 \end{bmatrix}.$$

The nonholonomic constraints may then be expressed by the equation P ⊥ ( v ) = 0 , from which we find the constrained velocities v D ∈ D is given by

$$v ^ { \mathcal { D } } = \left [ v _ { \theta } \ R \cos \theta v _ { \varphi } \ R \sin \theta v _ { \varphi } \ v _ { \varphi } \right ] ^ { T } \quad \quad ( 8 0 ) \\ \tau _ { - } \quad - \quad - \ \.$$

We construct the Rayleigh dissipation function according to equations (38) using the nonholonomic constraint form (74), from which the coefficient matrix is given by in the coordinate induced frame { ∂ ∂ θ , ∂ ∂ x , ∂ ∂ y , ∂ ∂ ϕ } .

$$\mathbb { F R } ^ { \ddagger _ { G } } = \begin{bmatrix} 0 & 0 & 0 \\ 0 & \frac { \mu } { m } & 0 & - \frac { \mu R } { m } \cos \theta \\ 0 & 0 & \frac { \mu } { m } & - \frac { \mu R } { m } \sin \theta \\ 0 & - \frac { \mu R } { J } \cos \theta & - \frac { \mu R } { J } \sin \theta & \frac { \mu R ^ { 2 } } { J } \end{bmatrix} \quad & ( 8 1 )$$

where the units of the friction coefficients µ are [ kg/s ] , and the corresponding accelerations given by F ε = -1 ε F R /sharp G v . Defining a frame Φ as the direct sum of the vector fields (76) and (77) we find the matrix Q q according to Remark 8 is given by

$$\mathbb { Q } _ { q } = \begin{bmatrix} 0 & 0 & 0 & 0 \\ 0 & J ^ { 2 } + \gamma m R ^ { 2 } \sin ^ { 2 } \theta & - \frac { 1 } { 2 } \gamma m R ^ { 2 } \sin 2 \theta & J ^ { 2 } R \cos \theta \\ 0 & - \frac { 1 } { 2 } \gamma m R ^ { 2 } \sin 2 \theta & J ^ { 2 } + \gamma m R ^ { 2 } \cos ^ { 2 } \theta & - J ^ { 2 } R \sin \theta \\ 0 & - m R J \cos \theta & - m R J \sin \theta & m R ^ { 2 } J \end{bmatrix}$$

where γ = mR 2 +2 J . One can check that Q q satisfies Lemma 4. As the Christoffel symbols are all zero, the equations of motion are simply given by

$$\frac { \text{d} } { \text{d} t } \left [ v _ { \theta } \\ v _ { y } \right ] & = - \frac { 1 } { \varepsilon } \left [ \begin{array} { c } \frac { \text{u} } { \text{u} } v _ { y } - \frac { \text{u} } { \text{u} } \cos \theta v _ { \varphi } \\ \frac { \text{u} } { \text{m} } v _ { y } - \frac { \text{u} } { \text{m} } \sin \theta v _ { \varphi } \\ - \frac { \text{u} } { \text{u} } \cos \theta v _ { x } - \frac { \text{u} } { \text{u} } \sin \theta v _ { y } + \frac { \text{u} } { \text{u} } \right ], \\ \frac { \text{d} } { \text{d} t } \left [ \begin{array} { c } \theta \\ x \\ y \\ \varphi \end{array} \right ] & = \left [ \begin{array} { c } v _ { \theta } \\ v _ { x } \\ v _ { y } \\ \text{u} \varphi \end{array} \right ].$$

Remark 15. Note that the dynamics of the orientation in the plane θ is of the form ˙ v θ = 0 as v θ ∈ D and D = ker ( F R /sharp G ) , by construction. Hence, the orientation in the plane is invariant for the dynamics of the vertical rolling disk.

## 5.1 Approximations of Slip Velocities

The first order approximation of the slip velocities are calculated according to (67) and are found as

$$\mathfrak { e } \text{vum } \mathfrak { a } \mathfrak { e } \\ \mathfrak { h } ^ { ( 1 ) } = \mathbb { Q } _ { q } \circ ( \nabla _ { v \mathcal { D } } \mathcal { P } ^ { \perp } ) ( v ^ { \mathcal { D } } ) = \begin{bmatrix} 0 \\ \frac { m R } { \mu } \sin \theta v _ { \theta } v _ { \varphi } \\ - \frac { m R } { \mu } \cos \theta v _ { \theta } v _ { \varphi } \\ 0 \end{bmatrix} \mathfrak { e } \text{v}, the second order approximation of the slip velocities is calculated according$$

Further, the second order approximation of the slip velocities is calculated according to (68) using the results of Remark 3 and is found to be

$$to ( 68 ) \, \text{using the results of} \, \text{Remark} \, 3 \, \text{ and is found to be} \\ h ^ { ( 2 ) } & = \mathbb { Q } _ { q } \circ \left [ ( \nabla _ { h ^ { ( 1 ) } } \mathcal { P } ^ { \perp } ) ( v ^ { \mathcal { P } } ) + ( \nabla _ { v ^ { \mathcal { D } } } \mathcal { P } ^ { \perp } ) ( h ^ { ( 1 ) } ) \right ] \\ & \quad - \mathbb { Q } _ { q } \circ \left [ D ^ { \mathcal { V } } h ^ { ( 1 ) } \left ( ( \nabla _ { v ^ { \mathcal { D } } } \mathcal { P } ) ( v ^ { \mathcal { D } } ) \right ) \right ] - \mathbb { Q } _ { q } \circ \left [ \frac { \mathcal { H } } { \nabla _ { v ^ { \mathcal { D } } } } h ^ { ( 1 ) } ( v ^ { \mathcal { D } } ) \right ] \\ & = \frac { m ^ { 2 } R J ^ { 2 } } { \mu ^ { 2 } ( J + m R ^ { 2 } ) ^ { 2 } } \left [ \substack { 0 \\ \sin \theta v ^ { 2 } _ { \theta } v _ { \varphi } \\ - \frac { m J } { J } v ^ { 2 } _ { \theta } v _ { \varphi } } \right ] \\ \\ 5. 2 \, \text{Approximations of} \, \text{Slip Dynamics}$$

## 5.2 Approximations of Slip Dynamics

Using the first order approximation of the slip velocities (85), we find the zeroth order approximation of (83)-(84) using equations (69) are given by

$$\text{ of } ( 8 ) \text{-} ( 8 4 ) \, \text{ using equations } ( 6 9 ) \, \text{ are given by} \\ \frac { d } { d t } \begin{bmatrix} v _ { \theta } \\ v _ { y } \\ v _ { \varphi } \end{bmatrix} = \begin{bmatrix} 0 \\ - R \sin \theta v _ { \theta } v _ { \varphi } \\ R \cos \theta v _ { \theta } v _ { \phi } \\ 0 \end{bmatrix} + O ( \varepsilon ), \quad & ( 8 7 ) \\ \frac { d } { d t } \begin{bmatrix} \theta \\ x \\ y \\ \varphi } \end{bmatrix} = \begin{bmatrix} v _ { \theta } \\ R \cos \theta v _ { \varphi } \\ R \sin \theta v _ { \varphi } \\ v _ { \varphi } \end{bmatrix} + O ( \varepsilon ). \quad & ( 8 8 ) \\ \text{ } \\ \text{The zeroth order equations of motion } ( 8 7 ) - ( 8 8 ) \, \text{ are the same equations} \\ \text{lima from Laarmane-d'Alembert's invariational nrincinite } I 1 \text{ See } I 1 \text{ }$$

Remark 16. The zeroth order equations of motion (87) (88) -are the same equations of motion resulting from Lagrange-d'Alembert's variational principle [1, Sec. 1.4].

Moreover, using the second order approximation of the slip velocities (85) and (86), we find the first order equations (71) are given by

$$\frac { \mathrm d } { \mathrm d t } \begin{bmatrix} v _ { \theta } \\ v _ { x } \\ v _ { y } \\ v _ { \varphi } \end{bmatrix} = \begin{bmatrix} 0 \\ - R \sin \theta v _ { \theta } v _ { \varphi } \\ R \cos \theta v _ { \theta } v _ { \varphi } \\ 0 \end{bmatrix} + \varepsilon \frac { \mathrm m R J } { \mu ( J + m R ^ { 2 } ) } \begin{bmatrix} 0 \\ - \cos \theta v _ { \varphi } v _ { \theta } ^ { 2 } \\ - \sin \theta v _ { \varphi } v _ { \theta } ^ { 2 } \\ \frac { \mathrm m R } { J } v _ { \varphi } v _ { \theta } ^ { 2 } \end{bmatrix} + O ( \varepsilon ^ { 2 } ), \quad ( 8 9 ) \\ \end{bmatrix}$$

$$\frac { d } { d t } \begin{bmatrix} \theta \\ x \\ y \\ \varphi \end{bmatrix} = \begin{bmatrix} v _ { \theta } \\ R \cos \theta v _ { \varphi } \\ R \sin \theta v _ { \varphi } \\ v _ { \varphi } \end{bmatrix} + \varepsilon \frac { m R } { \mu } \begin{bmatrix} 0 \\ \sin \theta v _ { \theta } v _ { \varphi } \\ - \cos \theta v _ { \theta } v _ { \varphi } \\ 0 \end{bmatrix} + O ( \varepsilon ^ { 2 } ).$$

## 6 Conclusion

In this paper, we developed an affine connection approach for realizing nonholonomic mechanical systems. We formulated the slow/fast dynamics in a coordinateindependent manner in terms of sections, projection maps and covariant derivatives. We proposed a novel invariance condition for slow/fast realizations of nonholonomic constraints and proved that this novel condition is equivalent to the classical approach. We proposed a recursive procedure to approximate higher order slip velocities and their associated dynamics. We conclude the paper with some remarks regarding the advantages of the proposed formalism.

- 1. The geometric formalism is coordinate-independent and suitable for symbolic computation.
- 2. The proposed slow-fast decomposition permits a recursive procedure to approximate the invariant manifold representing slip velocities. In the context of vehicle systems, the slow manifold may be physically interpreted as a deviation from the ideal no-slip constraints. Approximations of the slow manifold provide estimates of these deviations and have applications in model-based slip estimation and localization in vehicle systems.
- 3. Higher order approximations of slip velocities are based on objects that can be interpreted geometrically and may give insight into the geometric structure of other equivalent approaches.

The future directions of this research include the following: (i) Designing a recursive controller robust to the uncertainties in the estimation of the friction coefficients and external forces; (ii) Use of higher-order approximations of the slow-manifold for in-the-loop motion estimation of vehicle systems; and (iii) symmetry reduction for slow-fast realizations of nonholonomic systems.

Conflict of Interest: The authors declare that they have no conflict of interest.

## Appendix A Proofs of Statements

## A.1 Proof of Proposition 1

Proof. For the horizontal space

$$\mathcal { P } _ { \mathcal { H } } \left ( q, v _ { q }, U, W \right ) & = \left ( \mathcal { Q } _ { \Gamma } ^ { \mathcal { H } } \left ( q, v _ { q }, U, W \right ) \right ) ^ { \downarrow } \text{hor} \\ & = \left ( q, v _ { q }, U, - \Gamma ( q, v _ { q } ) U \right ) ^ { \downarrow } \text{hor} \\ & = \left ( q, U \right ) = \left ( \pi _ { T Q } \right ) _ { * } \left ( q, v _ { q }, U, W \right ).$$

Therefore we have P H = ( π TQ ) ∗ . For the vertical space, let X : Q → TQ be a section:

$$\mathcal { P } _ { \mathcal { V } } ( X _ { * } ) & = \left ( \mathcal { Q } _ { \Gamma } ^ { \mathcal { V } } ( q, X _ { q }, v _ { q }, D X ( q ) \cdot v _ { q } ) \right ) _ { \text{ver} } ^ { \downarrow } \\ & = ( q, X _ { q }, 0, D X ( q ) \cdot v _ { q } + \Gamma ( q, X _ { q } ) \cdot v _ { q } ) _ { \text{ver} } ^ { \downarrow } \\ & = ( q, D X ( q ) \cdot v _ { q } + \Gamma ( q, X _ { q } ) \cdot v _ { q } ) = \nabla _ { v } X.$$

## A.2 Proof of Proposition 2

Proof. Let h : TQ → TQ be a vector bundle map, i.e., π TQ ◦ h = π TQ . Let w v , : Q → TQ be sections, then h ◦ w : Q → TQ is also a section. The push-forward of this section is a map that satisfies the following commutative diagram:

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00095v1-figure-c92b40ca-0e95-4b4c-b0db-f2147b521762.jpg)

By virtue of Proposition (1), we calculate the covariant derivative along v as follows

$$\nabla _ { v } ( h \circ w ) & = \mathcal { P } _ { \mathcal { V } } \left ( D ( h \circ w ) \right ) = \mathcal { P } _ { \mathcal { V } } \left ( D h \cdot D w \right ) \\ & = \mathcal { P } _ { \mathcal { V } } \left ( D h \cdot \left ( \mathcal { P } _ { \mathcal { H } } \oplus \mathcal { P } _ { \mathcal { V } } \right ) ( D w ) \right ) \\ & = \mathcal { P } _ { \mathcal { V } } \left ( D h \cdot \left ( \mathcal { P } _ { \mathcal { H } } ( D w ) + \mathcal { P } _ { \mathcal { V } } ( D w ) \right ) \right ) \\ & = \mathcal { P } _ { \mathcal { V } } \left ( D h \cdot \left ( \mathcal { P } _ { \mathcal { H } } ( D w ) \oplus 0 \right ) \right ) + \mathcal { P } _ { \mathcal { V } } \left ( D h \cdot \left ( 0 \oplus \nabla _ { v } w \right ) \right ), \\ \text{where we used the direct product isomorphism, twice and } D h \text{ is the Jacobian of } h. \text{In}$$

where we used the direct product isomorphism, twice and D h is the Jacobian of h . In the decomposition TTQ ∼ H = TQ ⊕V TQ , consider the Jacobian of the map h :

$$\L _ { \ast } \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdots \\ \ D h = \begin{bmatrix} \mathcal { P } _ { \mathcal { H } } D h \\ \mathcal { P } _ { \mathcal { V } } D h \end{bmatrix} _ { \mathcal { H } } \begin{bmatrix} \mathcal { P } _ { \mathcal { H } } D h \\ \mathcal { P } _ { \mathcal { V } } D h \end{bmatrix} _ { \mathcal { H } } \begin{bmatrix} \mathcal { P } _ { \mathcal { H } } D h \\ \mathcal { P } _ { \mathcal { V } } D h \end{bmatrix} _ { \mathcal { V } } \end{bmatrix}. \quad \quad ( A 1 ) \\ \L _ { \ast }, \, \text{when restrict $Dh$ to the vertical space $\mathcal{V}TQ=ker(\pi_{TQ})_{*}$ we find}$$

Since P H = ( π TQ ) ∗ , when restrict D h to the vertical space V TQ = ker( π TQ ) ∗ we find

$$\pi _ { \mathcal { H } } D h \Big | _ { \mathcal { V } } = ( \pi _ { T Q } ) _ { * } h _ { * } \Big | _ { \mathcal { V } } = ( \pi _ { T Q } ) _ { * } \Big | _ { \mathcal { V } } = 0. & & ( \Lambda 2 ) \\ & & 2 6$$

Let z : [ a, b ] → TQ be a horizontal curve along q ( ) t in Q . Then, ˙ z ∈ T w ( ) t TQ ∼ = H z ( ) t TQ ⊕V z ( ) t TQ satisfying

$$\dot { z } = \dot { q } ^ { i } \left ( \frac { \partial } { \partial \tilde { q } ^ { i } } - \Gamma ^ { k } _ { i j } v ^ { j } \frac { \partial } { \partial v ^ { i } } \right ).$$

Since h z ( ( )) t ∈ T q ( ) t Q is a curve, we calculate derivative

$$\dot { q } ( t ) = \frac { \text{d} } { \text{d} t } \left ( \pi _ { T Q } \circ h \circ z ( t ) \right ) = ( \pi _ { T Q } ) _ { * } \, D h \, \dot { z } = \mathcal { P } _ { \mathcal { H } } ( D h \, \dot { z } ).$$

Therefore, we find that P H D h ∣ ∣ H is the identity map id. Define the vertical component of the Jacobian by D V h := P V ( D h ) | V and the horizontal covariant derivative ∇ H h := P V D h ∣ ∣ H . Therefore, the covariant derivative of the bundle map h is given by

$$\nabla _ { v } ( h \circ w ) = \overset { \mathcal { H } } { \nabla } _ { v } h ( w ) + D ^ { \mathcal { V } } h ( \nabla _ { v } w ).$$

We calculate local coordinate expressions for the covariant derivative as follows:

$$\text{e calculate local coordinate expressions for the covariant derivative as follows:} \\ \nabla _ { v } ( h \circ w ) & = v ^ { i } \nabla _ { \frac { \sigma } { \sigma ^ { i } } } \left ( ( h \circ w ) ^ { j } \frac { \partial } { \partial q ^ { j } } \right ) \\ & = v ^ { i } \left ( \mathcal { L } _ { \frac { \sigma } { \sigma ^ { i } } } ( h ( q, w ) ^ { j } ) \frac { \partial } { \partial q ^ { j } } + h ( q, w ) ^ { j } \nabla _ { \frac { \sigma } { \sigma ^ { i } } } \frac { \partial } { \partial q ^ { j } } \right ) \\ & = v ^ { i } \left ( \frac { \partial h ^ { j } } { \partial q ^ { i } } \frac { \partial } { \partial q ^ { j } } + \frac { \partial h ^ { j } } { \partial w ^ { k } } \frac { \partial w ^ { k } } { \partial q ^ { i } } \frac { \partial } { \partial q ^ { j } } + h ^ { i } \Gamma _ { i j } ^ { k } \frac { \partial } { \partial q ^ { k } } \right ). \\ \text{the other hand we have that}$$

On the other hand we have that

$$\nabla _ { \frac { \sigma } { \partial q ^ { i } } } w = \left ( \frac { \partial w ^ { k } } { \partial q ^ { i } } + \Gamma _ { i j } w ^ { j } \right ) \frac { \partial } { \partial q ^ { k } }.$$

Putting this expression into the above formula we find

$$\nabla _ { v } ( h \circ w ) & = \frac { \partial h ^ { m } } { \partial w ^ { l } } \left ( \nabla _ { v } w \right ) ^ { l } \frac { \partial } { \partial q ^ { m } } + \left ( \frac { \partial h ^ { m } } { \partial q ^ { k } } v ^ { k } + \Gamma ^ { m } _ { k i } v ^ { k } h ^ { i } - \Gamma ^ { l } _ { k i } w ^ { k } v ^ { i } \frac { \partial h ^ { m } } { \partial w ^ { l } } \right ) \frac { \partial } { \partial q ^ { m } } \\ & = D ^ { \mathcal { V } } h ( \nabla _ { v } w ) + \frac { \mathcal { H } } { \nabla } _ { v } h ( w ).$$

## References

- [1] Bloch, A., Brogliato, B.: Nonholonomic mechanics and control. Appl. Mech. Rev. 57 (1), 3-3 (2004)

| [2]   | Eldering, J.: Realizing nonholonomic dynamics as limit of friction forces. Regular and Chaotic Dynamics 21 (4), 390-409 (2016)                                                    |
|-------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [3]   | Pacejka, H.B., Bakker, E.: The magic formula tyre model. Vehicle system dynamics 21 (S1), 1-18 (1992)                                                                             |
| [4]   | Bates, L., Śniatycki, J.: Nonholonomic reduction. Reports on Mathematical Physics 32 (1), 99-115 (1993)                                                                           |
| [5]   | Bloch, A.M., Marsden, J.E., Zenkov, D.V.: Quasivelocities and symmetries in non-holonomic systems. Dynamical systems 24 (2), 187-222 (2009)                                       |
| [6]   | Bloch, A.: Nonholonomic mechanics, dissipation and quantization. Advances in the Theory of Control, Signals and Systems with Physical Modeling, 141-152 (2010)                    |
| [7]   | Bullo, F., Lewis, A.D.: Geometric Control of Mechanical Systems: Modeling, Analysis, and Design for Simple Mechanical Control Systems. Springer (2019)                            |
| [8]   | Petersen, P.: Riemannian geometry. Springer (2006)                                                                                                                                |
| [9]   | Lewis, A.D.: Affine connections and distributions with applications to nonholo- nomic mechanics. Reports on Mathematical Physics 42 (1-2), 135-164 (1998)                         |
| [10]  | Lewis, A.D.: Simple mechanical control systems with constraints. IEEE Transac- tions on Automatic Control 45 (8), 1420-1436 (2000)                                                |
| [11]  | Chhabra, R., Emami, M.: Nonholonomic dynamical reduction of open-chain multi-body systems: a geometric approach. Mechanism and Machine Theory 82 , 231-255 (2014)                 |
| [12]  | Schaft, A., Maschke, B.: On the hamiltonian formulation of nonholonomic mechanical systems. Reports on mathematical physics 34 (2), 225-233 (1994)                                |
| [13]  | Kane, T.R., Levinson, D.A.: Dynamics, theory and applications. McGraw Hill (1985)                                                                                                 |
| [14]  | Sannuti, P.: Singular perturbations in the state space approach of linear electrical networks. International Journal of Circuit Theory and Applications 9 (1), 47-57 (1981)       |
| [15]  | Zagaris, A., Kaper, H.G., Kaper, T.J.: Analysis of the computational singular perturbation reduction method for chemical kinetics. Journal of Nonlinear Science 14 , 59-91 (2004) |
| [16]  | Rubin, J.E., Terman, D.: Geometric singular perturbation analysis of neuronal dynamics. In: Handbook of Dynamical Systems vol. 2, pp. 93-146 (2002)                               |

| [17]   | Tikhonov, A.N.: Systems of differential equations containing small parameters in the derivatives. Matematicheskii sbornik 73 (3), 575-586 (1952)                                |
|--------|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [18]   | Verhulst, F.: Methods and applications of singular perturbations. Springer (2005)                                                                                               |
| [19]   | Fenichel, N., Moser, J.: Persistence and smoothness of invariant manifolds for flows. Indiana University Mathematics Journal 21 (3), 193-226 (1971)                             |
| [20]   | Fenichel, N.: Geometric singular perturbation theory for ordinary differential equations. Journal of differential equations 31 (1), 53-98 (1979)                                |
| [21]   | Eldering, J.: Normally Hyperbolic Invariant Manifolds: the Noncompact Case. Springer (2013)                                                                                     |
| [22]   | Marsden, J., Ratiu, T.: Introduction to mechanics and symmetry (1998)                                                                                                           |
| [23]   | Brendelev, V.: On the realization of constraints in nonholonomic mechanics. Journal of Applied Mathematics and Mechanics 45 (3), 351-355 (1981)                                 |
| [24]   | Karapetian, A.: On realizing nonholonomic constraints by viscous friction forces and celtic stones stability. Journal of Applied Mathematics and Mechanics 45 (1), 30-36 (1981) |
| [25]   | Minguzzi, E.: Rayleigh's dissipation function at work. European Journal of Physics 36 (3), 035014 (2015)                                                                        |
| [26]   | Antali, M., Stepan, G.: Nonsmooth analysis of three-dimensional slipping and rolling in the presence of dry friction. Nonlinear Dynamics 97 , 1799-1817 (2019)                  |
| [27]   | Antali, M., Stepan, G.: Sliding and crossing dynamics in extended filippov systems (2018)                                                                                       |
| [28]   | Filippov, A.F.: Differential Equations with Discontinuous Righthand Sides: Control Systems vol. 18. Springer, ??? (2013)                                                        |
| [29]   | Antali, M., Varkonyi, P.L.: The nonsmooth dynamics of combined slip and spin motion under dry friction. Journal of Nonlinear Science 32 (4), 58 (2022)                          |
| [30]   | Antali, M.: Bifurcations of the limit directions in extended filippov systems. Physica D: Nonlinear Phenomena 445 , 133622 (2023)                                               |
| [31]   | Ruina, A.: Nonholonomic stability aspects of piecewise holonomic systems. Reports on mathematical physics 42 (1-2), 91-100 (1998)                                               |
| [32]   | Bloch, A.M., Rojo, A.G.: Quantization of a nonholonomic system. Physical review letters 101 (3), 030402 (2008)                                                                  |

- [33] Conlon, L.: Differentiable manifolds. Springer (2001)

| [34]   | Lewis, A.D.: Nonholonomic and constrained variational mechanics. Journal of Geometric Mechanics 12 (2), 165 (2020)                                      |
|--------|---------------------------------------------------------------------------------------------------------------------------------------------------------|
| [35]   | Abraham, R., Marsden, J.E.: Foundations of Mechanics. American Mathematical Soc. (2008)                                                                 |
| [36]   | Ryan, J.M.: Geometry of horizontal bundles and connections. PhD thesis, Wichita State University (2014)                                                 |
| [37]   | Frankel, T.: The Geometry of Physics: an Introduction. Cambridge University Press (2011)                                                                |
| [38]   | Yano, K., Ishihara, S.: Differential Geometry of Tangent and Cotangent Bundles. Dekker, New York (1973)                                                 |
| [39]   | Dummit, D.S., Foote, R.M.: Abstract Algebra. Wiley Hoboken (2004)                                                                                       |
| [40]   | Samson, C., Morin, P., Lenain, R.: Modeling and control of wheeled mobile robots. Springer handbook of robotics, 1235-1266 (2016)                       |
| [41]   | Chhabra, R.: Dynamical reduction and output-tracking control of the lunar exploration light rover (lelr). In: IEEE Aerospace Conference, pp. 1-8 (2016) |
Vaughn Gzenda Robin Chhabra
Affine Connection Approach to the Realization of Nonholonomic Constraints by Strong Friction Forces
In this paper, we study an affine connection approach to realizing
nonholonomic mechanical systems mediated by viscous friction forces with large
coefficients, viewed as a singular perturbation of the nonholonomic system. We
show that the associated slow manifold is represented coordinate-free as the
image of a section over the nonholonomic distribution. We propose a novel
invariance condition based on covariant derivatives and prove that this
condition is equivalent to the classical invariance condition based on time
derivatives. Accordingly, we propose a novel recursive procedure to approximate
the slow manifold based on the covariant derivatives of a formal power series
expansion of the section. Using this recurrence relation, we derive, up to
second order, approximations of the slip velocities residing in the slow
manifold, as well as the associated approximated dynamics up to first order.
Lastly, we illustrate our approach with a case study of a vertical rolling
disk.