Patrick Ballard

The dynamics of discrete mechanical systems with perfect unilateral constraints

The dynamics of discrete mechanical systems with perfect unilateral constraints is formulated in a very general setting. The well-posedness of the resulting evolution problem is studied. It is proved that existence and uniqueness of a maximal solution is ensured provided strong assumptions on the regularity of the data: they are supposed to be analytic. Simple examples show that this regularity assumption may not be relaxed. Sufficient conditions to ensure that the maximal solution is defined for all time are supplied. The continuous dependence of the solution on initial conditions is also studied and the numerical computation of the solution is evocated.

Introduction

The aim of the Dynamics of Discrete Mechanical Systems (sometimes called Rational Mechanics or, after Lagrange, Analytical Mechanics) is the prediction of the motion of collections of bodies supposed to be perfectly indeformable. The theory classically distinguishes two types of interactions between the bodies themselves and between the bodies and the rest of the universe: the efforts and the constraints. The constraints are kinematical specifications of the motion with which some efforts are associated. A constraint is said perfect or ideal if the associated efforts do not dissipate energy. A constraint is said bilateral (respectively unilateral) if the kinematical specification gives rise to equalities (respectively inequalities). A typical occurrence of unilateral constraints is the handling of non-penetration conditions.

When all the constraints are bilateral and perfect, the motion is classically governed by a second-order ordinary differential equation on a finite dimensionnal riemannian manifold. When the data are smooth enough, Cauchy-Lipschitz theorem guarantees that a unique motion is associated with any given initial state of the system.

When dealing with the dynamics of discrete mechanical systems with unilateral constraints, there is no such a theorem, although many steps towards this direction have been made during the past twenty years. To my knowledge, the first investigation of this question using modern mathematical tools (i.e. introducing motions whose acceleration is a measure with respect to time) is that of Schatzman [START_REF] Schatzman | A Class of Nonlinear Differential Equations of Second Order in Time, Nonlinear Analysis[END_REF]. She studied the particular case where the configuration space is R d equipped with its canonical euclidean structure and the admissible configuration set is convex. Her setting was also limited to the elastic impact constitutive equation. Using Yosida type regularization and compactness arguments, she was able to prove the existence of solutions under very weak regularity assumptions. She also discussed uniqueness but proved it only in a very specific case. Further investigation on uniqueness was performed by Percivale in [START_REF] Percivale | Uniqueness in the Elastic Bounce Problem, I[END_REF] and [START_REF] Percivale | Uniqueness in the Elastic Bounce Problem, II[END_REF]. He is the first to introduce analyticity hypothesis in this respect. But, his results apply also only to very specific cases. The formulation of the problem with completely inelastic impacts has been extensively studied by Moreau [START_REF] Moreau | Standard inelastic shocks and the dynamics of unilateral constraints, in Unilateral problems in structural analysis[END_REF]. An existence result was proved by Monteiro-Marques [START_REF] Monteiro Marques | Differential Inclusions in Nonsmooth Mechanical Problems[END_REF] in the particular case in which the configuration space is euclidean R d and the unilateral constraints are described by a single smooth function. Very recently, Schatzman [START_REF] Schatzman | Uniqueness and continuous dependence on data for one dimensional impact problems[END_REF] studied the general one degree-offreedom problem with arbitrary impact constitutive law. In this case, she proved uniqueness under analyticity assumption on the data.

None of these results has the generality required by Mechanics. The existence and uniqueness results are proved under assumptions which are obviously not fulfilled in most of discrete mechanical systems which may generally be encountered, except the last result of Schatzman, but it is limited to the one degree-of-freedom problem.

In this paper, the dynamics of discrete mechanical systems with perfect unilateral constraints is formulated in a very general setting. To reach full generality, the configuration space is supposed to be an arbitrary riemannian manifold instead of an euclidean space. However, only the most elementary level of differential geometry is needed. The resulting general evolution problem is studied. The existence and uniqueness of a solution associated with given initial condition is proved provided the data are analytic.

In section 2, we give a precise mathematical definition of what we call discrete mechanical system and system of bilateral constraints. We also recall some basic results connected to these definitions that we shall use subsequently.

In section 3, a formulation of the equations of the dynamics of discrete mechanical systems with perfect unilateral constraints is presented. The content of this section follows very closely the work of Moreau [START_REF] Moreau | Standard inelastic shocks and the dynamics of unilateral constraints, in Unilateral problems in structural analysis[END_REF]. It is included since Moreau restricts himself to completely inelastic impacts. More generality, including the case of elastic impacts, is obtained here with no supplementary difficulty.

In section 4, we prove a local existence and uniqueness result concerning the general problem of the dynamics of discrete mechanical systems with perfect unilateral constraints, under the single assumption that the data are analytic. Existence and uniqueness of a maximal solution follows immediately. A sufficient condition to ensure that this maximal solution is defined for all time is also presented.

In section 5, three examples are discussed. One is due to Moreau and another one to Schatzman. They are included for sake of completeness. The aim of these examples is to show that the regularity assumptions made in the previous section are, in some sense, minimal.

In section 6, we illustrate the generality of the theorems of section 3 in applying them to simple examples issuing from Mechanics.

In section 7, the continuous dependence of the solution on initial conditions is discussed. Dependence on initial conditions is seen to be not continuous in general. However, a restrictive case where continuity holds is exhibited.

In section 8, the numerical computation of the solution is discussed. Problems arise in connection with non-continuous dependence on initial conditions. However, we recall an algorithm, which was first described by Moreau, and prove its convergence in some restrictive cases.

The main results in this paper were announced in Ballard [START_REF] Ballard | Dynamique des systèmes mécaniques discrets avec liaisons unilatérales parfaites[END_REF].

Discrete mechanical systems and perfect bilateral constraints

The aim of this section is to give a precise definition of what we call a discrete mechanical system, to introduce notations and to recall some basic results that we shall use later on. For a comprehensive presentation, the reader is referred to Arnold [START_REF] Arnold | Mathematical Methods of Classical Mechanics[END_REF] and Abraham & Marsden [START_REF] Abraham | Foundations of Mechanics[END_REF].

Discrete mechanical systems

Definition 1 A discrete mechanical system is:

-A Hausdorff, smooth (of class C p with 2 ≤ p ≤ ∞) connected manifold Q of dimension d whose topology has a countable basis. Q is called the configuration space of the discrete mechanical system. d is its number of degrees of freedom. The tangent bundle T Q of Q is called the phase space or the state space. A point q of Q is a configuration of the system and a point of T Q a state of the system. T * Q denotes the cotangent bundle,

Π Q : T Q → Q and Π * Q : T * Q → Q the natural projections.
T q Q is, as usual, the tangent space at q and, to designate an element v of T Q, we shall often use the redundant notation (q, v) where q = Π Q (v) and v ∈ T q Q. A curve on Q (i.e. a continuous mapping from a real interval I to Q) is also called a motion of the system. If a motion q : I → Q admits a tangent vector at t, it will be denoted by (q(t), q(t)). This notation is an abuse consecrated by tradition. The dot will also be used in general to denote a derivative with respect to time. A local chart on Q is also called a local parametrization of the system.

-A riemannian metric on Q denoted by (•, •) q . The mapping

K T Q → R + (q, v) → 1 2 (v, v) q = 1 2 v 2 q (1)
is the kinetic energy of the system.

-A real interval I and a smooth (of class C p with 1 ≤ p ≤ p) mapping f : T Q × I → T * Q such that:

∀(q, v) ∈ T Q, ∀t ∈ I, Π * Q (f (q, v; t)) = Π Q (q, v) = q.
The mapping f is called virtual power of internal, external and inertial efforts acting on the system or, in short, the efforts mapping. We will denote by •, • q the local duality product on T * q Q × T q Q and (and = -1 its inverse) the isomorphism of vector bundles from T Q onto T * Q canonically associated with the riemannian metric on Q.

The Fundamental Principle of Dynamics asserts that any motion of the system is of class C 2 and has to satisfy: ∀t ∈ I, D q(t) dt = f (q(t), q(t); t) [START_REF] Arnold | Mathematical Methods of Classical Mechanics[END_REF] where D dt denotes the operator of covariant derivation along q(t) canonically associated with the riemannian metric of Q.

In the sequel, for (U, ψ) a local chart on Q, (e 1 (q), e 2 (q), • • • , e d (q)) and e 1 (q), e 2 (q), • • • , e d (q) will denote the dual basis of T q Q and T * q Q naturally associated with the considered chart. ψ(q), that we shall abusively continue to denote by q, is an element q 1 , q 2 , • • • , q d of R d . If q(t) is a smooth motion on Q, q1 (t), q2 (t), • • • , qd (t) will be the components of its tangent vector (also said velocity) in the local basis: q(t) = qi (t)e i (q(t)), where Einstein's summation convention applies. It will always apply unless explicitly stated. No confusion induced by this notation should be expected since:

∀i ∈ {1, 2, • • • , d} , qi (t) = d dt q i (t).
In general, we shall use the same notation to denote a function and its representative in a chart. As usual, g ij (q) will denote the covariant components of the metric in the considered chart and g ij (q) its contravariant components. Γ i jk (q) will be the associated Christoffel symbols:

Γ i jk (q) = 1 2
g ih (q) ∂g hk ∂q j (q) + ∂g jh ∂q k (q) -∂g jk ∂q h (q) .

(

) 3
Proposition 2 (Lagrange) Let (U, ψ) be a local chart and q(t) a C 2 motion on Q. One has:

D q(t) dt = d dt ∂ ∂ qi K(q(t
), q(t)) -∂ ∂q i K(q(t), q(t)) e i (q(t)).

Proof. It is straightforward since: Coming back to the equation of motion [START_REF] Arnold | Mathematical Methods of Classical Mechanics[END_REF], suppose we are given in supplement an element t 0 of I, said initial instant, and an element (q 0 , v 0 ) of T Q, said initial state. Then, we obtain the following Cauchy problem C on Q:

D q dt =
C ⎧ ⎨ ⎩ D q dt
= f (q(t), q(t); t) (q(t 0 ), q(t 0 )) = (q 0 , v 0 ) .

Cauchy-Lipschitz theorem guarantees existence and uniqueness of a maximal C 2 solution (J m , q m ) where J m is an open subinterval of I including t 0 and q m a C 2 motion defined on J m . This expresses the fact that any other solution (J, q) of C is necessarily a restriction of q m : J ⊂ J m and q m|J = q

This result allows us to associate with any discrete mechanical system a dynamical system, that is a two real parameters collection F s,t of mappings from T Q into T Q such that:

F t 3 ,t 2 • F t 2 ,t 1 = F t 3 ,t 1 and F t,t = Id.
To illustrate these basic definitions and results, we give a simple example that we shall reuse later on in a slightly different context. Consider a plane system of two homogeneous rigid bars 1 and 2. The bar 1, of length l 1 and mass m 1 is connected to a fixed support by means of a perfect ball-and-socket joint equipped with a spiral spring of stiffness k 1 . The bar 2, of length l 2 and mass m 2 is connected to the free extremity of the bar 1 by means of another ball-and-socket joint also equipped with a spiral spring of stiffness k 2 . A force acts on the free extremity of the bar 2. This force remains parallel to the direction of the bar 2 and is of constant magnitude λ > 0 (see figure 2.1). With this system is associated the following discrete mechanical system: -The configuration space is R 2 equipped with its canonical structure of C ∞ manifold (it is not the 2-torus since the spiral springs impose to be able to count the 'number of turns'). This manifold may be represented by a single chart; in other terms, there exists a global parametrization of the system. In the sequel, we shall only use the chart (q 1 , q 2 ) defined by the angular measures associated with each of the joints. -The kinetic energy is:

l , m 1 1 1 q 1 k k 2 l , m 2 2 2 q λ
K = 1 2 l 1 0 m 1 l 1 s 2 q1 2 ds + 1 2 l 2 0 m 2 l 2 l 2 1 q1 2 + s 2 q2 2 + 2l 1 s cos q 1 -q 2 q1 q2 ds = 1 2 1 3 m 1 l 2 1 q1 2 + m 2 l 2 1 q1 2 + 1 3 m 2 l 2 2 q2 2 + m 2 l 1 l 2 cos q 1 -q 2 q1 q2 .
This kinetic energy defines a riemannian structure on the configuration space. The expression of the metric tensor in the considered chart is:

g 11 q 1 , q 2 = m 1 3 + m 2 l 2 1 , g 12 q 1 , q 2 = 1 2 m 2 l 1 l 2 cos q 1 -q 2 = g 21 q 1 , q 2 , g 22 q 1 , q 2 = 1 3 m 2 l 2 2 .
-The efforts mapping has for expression in the considered chart:

f (q, q; t) = λl 1 sin q 1 -q 2 -(k 1 + k 2 ) q 1 + k 2 q 2 e 1 (q) + k 2 q 1 -k 2 q 2 e 2 (q).
Proposition 2 allows us to form easily the equation of motion in the considered chart:

⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ m 1 3 + m 2 l 2 1 q1 + 1 2 m 2 l 1 l 2 cos q 1 -q 2 q2 + 1 2 m 2 l 1 l 2 sin q 1 -q 2 q2 2 = λl 1 sin q 1 -q 2 -(k 1 + k 2 ) q 1 + k 2 q 2 1 2 m 2 l 1 l 2 cos q 1 -q 2 q1 + 1 3 m 2 l 2 2 q2 -1 2 m 2 l 1 l 2 sin q 1 -q 2 q1 2 = k 2 q 1 -q 2
(4) The deterministic conclusion of Cauchy-Lipshitz theorem on the dynamic evolution of the system is illusive. Indeed, if we add to the differential system (4) the initial condition:

q 1 (0) = q 2 (0) = q1 (0) = q2 (0) = 0,
it is easily seen that the maximal solution is the identically vanishing function on the real line. But, Poincaré-Lyapunov theory shows that this solution is unstable for some value of λ and the real motion will differ in this case from this trivial solution.

The correct analysis of the motion should in this case refer to some investigation of topological nature on the dynamical system generated by the equation of motion.

In any case, one has to abandon the objective of predicting exactly the motion of the system. One has to accept to obtain only partial information on this motion: this is a consequence of the over-idealization made during the modelling process. However, Cauchy-Lipschitz theorem is at the basis of any further analysis which has to be performed on the equation of motion. This fact will be discussed with more details in section 7 in the context of the dynamics of discrete mechanical systems with perfect unilateral constraints.

Bilateral constraints

One may introduce on discrete mechanical system another type of efforts, not taken into account by the efforts mapping f . Indeed, one may specify some efforts by their kinematical effects: one speaks of constraint. A constraint induces a restriction on the admissible motions of the system which is expressed by means of a finite number n of smooth real functions defined on Q:

∀i ∈ {1, 2, • • • , n} , ϕ i (q) = 0 (5)
The word constraint in the singular will be used indifferently to speak either of a constraint specifically associated with a single function ϕ i or of the constraint associated with all the functions ϕ i . In this terminology, a set of constraints is still a constraint. In formula [START_REF] Brezis | Opérateurs Maximaux Monotones et Semi-groupes de Contrations dans les Espaces de Hilbert[END_REF], the constraint is said holonomic (because it applies on the configuration and not on the state), scleronomic (because it does not depend explicitly on time) and bilateral (because it is expressed only by equalities and not inequalities). We denote by S the following subset of Q:

S = {q ∈ Q ; ∀i ∈ {1, 2, • • • , n} , ϕ i (q) = 0} ,
and we add the assumption that the functions ϕ i are functionally independent: for all q in S, the dϕ i (q

) (i ∈ {1, 2, • • • , n}) are linearly independent in T * Q. As a consequence, S is a submanifold of Q of dimension d -n.
The realization of kinematical specifications (5) necessarily involves a virtual power of reaction efforts mapping R taking values in T * Q. It is a priori unknown. Now, consider an initial instant t 0 in I and an initial state (q 0 , v 0 ) compatible with the constraint (i.e. (q 0 , v 0 ) ∈ T S ⊂ T Q). The evolution problem associated with the discrete mechanical system with bilateral constraint is: find

T > t 0 , q ∈ C 2 ([t 0 , T [; Q) and R ∈ C 0 ([t 0 , T [; T * Q) such that: ⎧ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎩ ∀t ∈ [t 0 , T [, D q(t) dt = f (q(t), q(t); t) + R(t), ∀t ∈ [t 0 , T [, q(t) ∈ S,
(q(t 0 ), q(t 0 )) = (q 0 , v 0 ).

These equations fail to determine the motion of the system: one has to precise the mapping R by means of a phenomenological assumption on the way the constraint acts. A constraint will be said perfect if the associated reaction efforts do not produce work in any virtual velocity compatible with the constraint:

∀v ∈ {v ∈ T q M ; ∀i ∈ {1, 2, • • • , n} , dϕ i (q), v q = 0} T S, R, v q = 0.
As a result:

∃ (λ i ) i=1,2,•••,n ∈ R n R = n i=1 λ i dϕ i (q).
Therefore, if the bilateral constraint is perfect, the evolution problem may be written as: find

T > t 0 , q ∈ C 2 ([t 0 , T [; Q) and (λ i ) i=1,2,•••,n ∈ C 0 ([t 0 , T [; T * Q) n such that: E Q ⎧ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎩ ∀t ∈ [t 0 , T [, D Q q(t) dt = f (q(t), q(t); t) + n i=1 λ i (t)dϕ i (q(t)), ∀t ∈ [t 0 , T [, q(t) ∈ S, (q(t 0 ), q(t 0 )) = (q 0 , v 0 ),
where

D Q
dt is the operator of covariant derivation on Q. Let q be a point of Q, v a vector of T q Q, and E a subspace of T q Q. The orthogonal projection of v on E for the scalar product of T Q induced by the riemannian structure of Q will be denoted by Proj q [v; E]. Similarly, Proj * q [v * ; E * ] will denote the orthogonal projection of the 1-form v * on the subspace E * of T * q Q. Then, consider the evolution problem E S : find T > t 0 and q ∈ C 2 ([t 0 , T [; S) such that:

E S ⎧ ⎨ ⎩ ∀t ∈ [t 0 , T [, D S q(t) dt = Proj * q(t)
f (q(t), q(t); t); T * q(t) S , (q(t 0 ), q(t 0 )) = (q 0 , v 0 ), where T * q S is considered as a subspace of T * q Q and D S dt is the operator of covariant derivation on S equipped with the riemannian structure inherited from that of Q. We have: Proof. First, let us identify T q S and T * q S to subspaces of T q Q and T * q Q. One has T * q S = (T q S). T * q S and n i=1 R dϕ i (q) are complementary orthogonal subspaces of T * q Q and (Chavel [START_REF] Chavel | Riemannian Geometry: a Modern Introduction[END_REF], p. 54):

D S q dt = Proj q D Q q dt ; T q S .
Now, let q be a solution of E Q :

Proj * q D Q q dt ; T * q S = Proj * q f (q, q; t) + n i=1 λ i dϕ i (q); T * q S .
But, Proj * q f (q, q; t) + n i=1 λ i dϕ i (q); T * q S = Proj * q f (q, q; t); T * q S , and,

Proj * q D Q q dt ; T * q S = Proj q D Q q dt ; T q S = D S q dt ,
which show that q is a solution of E S .

Reciprocally, let q be a solution of E S . From:

D S q dt = D Q q dt + n i=1 α i dϕ i (q), Proj * q f (q, q; t); T * q S = f (q, q; t) + n i=1 β i dϕ i (q),
we deduce the existence of n functions λ i : [t 0 , T [→ R such that:

D Q q dt = f (q, q; t) + n i=1 λ i dϕ i (q)
It follows that:

⎛ ⎜ ⎜ ⎝ . . . λ i . . . ⎞ ⎟ ⎟ ⎠ = ⎛ ⎜ ⎜ ⎝ . . . • • • (dϕ i (q), dϕ j (q)) q • • • . . . ⎞ ⎟ ⎟ ⎠ -1 ⎛ ⎜ ⎜ ⎜ ⎝ . . . D Q q dt -f (q, q; t), dϕ i (q) q . . . ⎞ ⎟ ⎟ ⎟ ⎠
where the Gram matrix is invertible because of the assumption on the functions ϕ i . This shows that the function λ i are uniquely determined and that they are continuous. Therefore, q generates a solution of E Q .

The second part of Proposition 3 follows from standart results on ordinary differential equations (see for example Coddington & Levinson [8]).

The moral of Proposition 3 is that adding a perfect bilateral constraint to a discrete mechanical system generates another discrete mechanical system with smaller number of degrees of freedom.

Discrete mechanical system with perfect unilateral constraints

This section deals with the formulation of the equation of motion of a discrete mechanical system when some perfect unilateral constraints are added. All the basic ideas of this section are due to Moreau [START_REF] Moreau | Standard inelastic shocks and the dynamics of unilateral constraints, in Unilateral problems in structural analysis[END_REF]. It is included since Moreau restricts himself to the special case of completely inelastic impacts and also because Moreau does not consider the general case of an arbitrary configuration manifold equipped with an arbitrary riemannian structure.

Kinematical setting

Consider a discrete mechanical system according to section 2.1 and suppose that a finite number n of unilateral constraints are taken into account:

∀i ∈ {1, 2, • • • , n} , ϕ i (q) ≤ 0, (6)
where the ϕ i : Q → R are C 1 functions. The closed subset A of Q defined by:

A = {q ∈ Q ; ∀i ∈ {1, 2, • • • , n} , ϕ i (q) ≤ 0}
is called the admissible configuration set. We define the mapping J by:

J Q → P ({1, 2, • • • , n}) q → J(q) = {i ∈ {1, 2, • • • , n} ; ϕ i (q) ≥ 0}
where P ({1, 2, • • • , n}) denotes the set of all subsets of {1, 2, • • • , n}. The set J(q) is called the set of all active constraints in the configuration q. As in the case of bilateral constraints, a functionally independence assumption is made on the functions ϕ i :

∀q ∈ A, (dϕ i (q)) i∈J(q) is linear independent in T * q Q (7)
As an easy consequence of the regularity assumptions made on the functions ϕ i , the boundary ∂A and the interior

• A of A in Q are such that: ∂A ⊂ n i=1 ϕ -1 i ({0}) , ( 8
) • A = J -1 ({∅}) (9)
Consider a motion in A (i.e. a continuous mapping from a real interval I to A) and assume that a right velocity q+ (t) ∈ T q(t) Q exists for all instant t of I. We necessarily have:

∀i ∈ {1, 2, • • • , n} , ∀t ∈ I, ϕ i (q(t)) = 0 =⇒ dϕ i (q(t)), q+ (t) q(t) ≤ 0 or, equivalently, ∀i ∈ {1, 2, • • • , n} , ∀t ∈ I, ϕ i (q(t)) = 0 =⇒ ∇ϕ i (q(t)), q+ (t) q(t) ≤ 0
where ∇ϕ i (q) is the gradient of ϕ i at q defined by:

∇ϕ i (q) = (dϕ i (q)) .
Thus, if the system has configuration q, then the right velocity q+ is necessarily in the closed convex cone V (q) of T q Q defined by:

V (q) = {v ∈ T q Q ; ∀i ∈ J(q), dϕ i (q), v q ≤ 0} .
V (q) is called the cone of admissible right velocities at the configuration q. In particular,

q ∈ • A (i.e. J(q) = ∅) =⇒ V (q) = T q Q.
Similarly, if a left velocity q-∈ T q Q exists, then, q-∈ -V (q)

Equation of motion

As for bilateral constraints, the realization of the constraints induces some reaction effort R. The following hypothesis are made:

-H1: the unilateral constraints are of type contact without adhesion:

∀v ∈ V (q), R, v q ≥ 0
-H2: the unilateral constraints are perfect: ∀v ∈ {v ∈ T q M ; ∀i ∈ J(q), dϕ i (q), v q = 0} , R, v q = 0.

There results from hypothesis H1 and H2 and Farkas' lemma (see e.g. Rockafellar [START_REF] Rockafellar | Convex Analysis[END_REF], p. 200) that:

∃ (λ i ) i=1,2,•••,n ∈ R n R = n i=1 λ i dϕ i (q), i ∈ J(q) ⇒ λ i ≤ 0, i ∈ J(q) ⇒ λ i = 0.
Thus, the reaction effort R ∈ T * Q must be such that:

-R ∈ N * (q) def = n i=1
λ i dϕ i (q) ; ∀i ∈ J(q), λ i ≥ 0, ∀i ∈ J(q), λ i = 0 .

(10) N * (q) is a closed convex cone of T * q Q and it is the polar cone of V (q) in the duality T q Q, T * q Q . We will also have to consider the polar cone N (q) of V (q) for the euclidean structure of T q Q: N (q) = n i=1 λ i ∇ϕ i (q) ; ∀i ∈ J(q), λ i ≥ 0, ∀i ∈ J(q), λ i = 0 . Now, consider a motion q(t) starting at q 0 ∈ • A at time t 0 with velocity v 0 . Assumed to be continuous, q(t) remains in • A on a right neighbourhood of t 0 . By formula [START_REF] Monteiro Marques | Differential Inclusions in Nonsmooth Mechanical Problems[END_REF], the reaction effort R vanishes as long as q(t) is in • A and the motion is governed by the ordinary differential equation:

⎧ ⎨ ⎩ D q dt = f (q, q; t), (q(t 0 ), q(t 0 )) = (q 0 , v 0 ).
Suppose that the solution of this Cauchy problem meets ∂A at some instant greater than t 0 . Denote by T the smallest of these instants. The motion admits a left velocity vector v - T at time T . Of course, there may happen: v - T ∈ V (q(T )). In this case, no differentiable prolongation of the motion can exist in A for t greater than T . The requirement of differentiability has to be dropped. An instant such T is called an instant of impact. However, we are going to still require the existence of a right velocity vector q+ (t) ∈ V (q(t)) at every instant t. The right velocity need not to be a continuous function of time and the equation of motion

D q+ dt = f (q, q+ ; t) + R,
should be understood in sense of Schwartz's distribution. Actually, we require R to be a vector valued measure rather than a general distribution. We denote by M M A(I; Q) (motions with measure acceleration) the set of all absolutely continuous motions q(t) from a real interval I to Q admitting a right velocity q+ (t) at every instant t of I and such that the function q+ (t) has locally bounded variation over I. Naturally, bounded variation is classically defined only for functions taking values in a normed vector space. However, for any absolutely continuous curve q(t) on a riemannian manifold, parallel translation along q(t) classically provides intrinsic identification of the tangent spaces at different points of the curve and so, the definitions can easily be carried over to this case. The precise mathematical setting is postponed to the appendix. The reader will notice from the appendix that with any motion q ∈ M M A(I; Q) is intrinsically associated the covariant Stieljes measure D q+ of its right velocity q+ . The equation of motion takes the form:

D q+ = f (q, q+ ; t) dt + R,
where dt denotes the Lebesgue measure. We have to give a precise meaning to condition [START_REF] Monteiro Marques | Differential Inclusions in Nonsmooth Mechanical Problems[END_REF] with R being a vector valued measure. By convention, we shall write: R ∈ -N * (q(t))

to mean: if θ ∈ L 1 loc (I, q, |R| ; T * Q) is the density of measure R with respect to its modulus measure |R| defined by proposition 25 of the appendix, then, one has:

θ(t) ∈ -N * (q(t))
for |R| -a.a.

t ∈ I. ( 11
)
This requirement is easily seen to be equivalent to the requirement of the existence of n nonpositive real measures λ i such that:

R = n i=1 λ i dϕ i (q(t)), ∀i ∈ {1, 2, • • • , n} , Supp λ i ⊂ {t; ϕ i (q(t)) = 0} . ( 12
)
Using this convention, the final form of the equation of motion is:

R = D q+ -f (q(t), q+ (t); t) dt ∈ -N * (q(t)) (13)

The impact constitutive equation

We begin this section by an example. Consider the one degree-of-freedom mechanical system whose configuration space is R equipped with its canonical euclidean structure. The efforts mapping f vanishes identically and the unilateral constraint is represented by the single function ϕ 1 (q) = q so that the admissible configuration set A is R -. At initial time t 0 = 0, we consider an initial state (q 0 , v 0 ) such that q 0 < 0 and v 0 > 0. It is readily seen from the equation of motion (13) that an impact necessarily occurs at time t = -q 0 /v 0 . At this time, the left velocity is v 0 . But, the right velocity can take any negative value and whatever it is, it is compatible with the equation of motion.

The reason for this indetermination lies in the phenomenological nature of the interaction of the system with the obstacle. Thus, we are led to make the general following hypothesis:

-H3: the interaction of the system with the obstacle at time t is completely determined by the present configuration q(t) and the present left velocity q-(t).

In other terms, we postulate the existence of a mapping F : T Q → T Q describing the interaction of the system with the obstacle during an impact:

∀t, q+ (t) = F q(t), q-(t) (14)
To ensure compatibility with the equation of motion [START_REF] Moreau | Bounded variation in time[END_REF], the mapping F should satisfy:

∀q ∈ A, ∀v -∈ -V (q), F (q, v -) ∈ V (q), F (q, v -) -v -∈ -N (q). (15)
First, consider the particular case of a motion with no more than one active constraint at any time (∀t, CardJ(q(t)) ≤ 1). The normal cone N (q(t)) is either {0} or a half-line and hypothesis H3 is equivalent to postulate the existence of an impact function φ : T Q → R such that: ∀t, q+ (t) = q-(t) -1 + φ q(t), q-(t) Proj q(t) q-(t); N (q(t)) . ( 16)

Equation ( 16) admits the equivalent form:

q+ (t) = Proj q(t)
q-(t); V (q(t))φ q(t), q-(t) Proj q(t) q-(t); N (q(t)) .

(17) For the general case where more than one constraint may be active at a time, we recall the following (Moreau [START_REF] Moreau | Décomposition orthogonale d'un espace hilbertien selon deux cones mutuellement polaires[END_REF]):

Lemma 4 (Moreau) Let V and N be two closed convex polar cones of a real Hilbert space H. Then,

∀x ∈ H, x = Proj[x; V ] + Proj[x; N ] and (Proj[x; V ], Proj[x; N ]) H = 0
As a consequence, the 'impact constitutive equations' ( 16) and ( 17) still make sense and are still equivalent when more than one constraint may be active at a time. Therefore, it is natural to retain only the particular forms ( 16) and ( 17) of the general impact constitutive equation [START_REF] Percivale | Uniqueness in the Elastic Bounce Problem, I[END_REF]. As a result of this further hypothesis, the phenomenology of the interation of the system with the obstacle during an impact is described by the single impact function φ : T Q → R. The impact function is also often called 'restitution coefficient'. Naturally, the impact function φ cannot be arbitrary and has to satisfy some consistency conditions. For example, the normality condition in [START_REF] Percivale | Uniqueness in the Elastic Bounce Problem, II[END_REF] requires:

∀q, q-, φ q, q-≥ -1.

But, this is not enough, we have to impose supplementary conditions on φ in order to ensure:

q-∈ -V (q) =⇒ q+ ∈ V (q). (18)
With respect to this, we have:

Proposition 5 Let V and N be two closed convex polar cones of a real Hilbert space H. Consider v -∈ -V such that Proj[v -; N ] = 0 and φ ∈ R. Then, v + = v --(1 + φ)Proj[v -; N ] ∈ V ⇐⇒ [φ ≥ 0]
Proof. For the 'if' part, suppose φ ≥ 0. By lemma 4, one gets:

Proj[v -; N ] = v --Proj[v -; V ] ∈ -V. But, v + = Proj[v -; V ] + φ -Proj[v -; N ] ,
and therefore, v + ∈ V , since V is a convex cone.

For the 'only if' part, we have by hypothesis,

Proj[v -; V ] -φProj[v -; N ] ∈ V.
Evaluating the scalar product with Proj[v -; N ] and using lemma 4, one gets:

-φ Proj[v -; N ]
2 H ≤ 0, and therefore the desired conclusion φ ≥ 0.

There results from proposition 5, that we have to require that the impact function φ should be nonnegative. This consistency assumption ensures that conditions [START_REF] Percivale | Uniqueness in the Elastic Bounce Problem, II[END_REF] and [START_REF] Schatzman | A Class of Nonlinear Differential Equations of Second Order in Time, Nonlinear Analysis[END_REF] will automatically be fulfilled.

At this stage, it should be underlined that hypothesis H3 implies the general forms [START_REF] Rockafellar | Convex Analysis[END_REF] or [START_REF] Rudin | Real and complex analysis[END_REF] for the impact constitutive equation only in the restrictive case where only at most one constraint is active at a time. In case of multiple impacts, the choice we made is only motivated by aesthetic considerations and also to fix ideas, since the concept of restitution coefficient is so firmly anchored in minds. We shall discuss more completely the relevance of that choice in section 6.4. Now, let us examine another example. Consider the one degree of freedom discrete mechanical system whose configuration space is R equipped with its canonical structure of riemannian manifold. The efforts mapping is supposed to be constant: f (q, q; t) ≡ 2. To this discrete mechanical system, we add the unilateral constraint described by the single function ϕ 1 (q) = q. Thus, A = R -. The impact constitutive equation is given by formula [START_REF] Rockafellar | Convex Analysis[END_REF] where the impact function is supposed to be the constant 1/2: φ ≡ 1/2. This mechanical system is a formal description of the physical occurence of a single particle subjected to gravity and bouncing on the floor. Consider the initial instant t 0 = 0 and the initial state (q 0 , v 0 ) = (-1, 0). It is readily seen that the function q : R + → R -defined by:

∀t ∈ [0, 1], q (t) = t 2 -1, ∀t ∈ [1, 2], q (t) = t 2 -3t + 2, ∀t ∈ 3 -1 2 n-1 , 3 -1 2 n , q(t) = t 2 + -6 + 3 2 n t + 3 -1 2 n-1 3 -1 2 n , ∀t ∈ [3, +∞[, q (t) = 0 (n ∈ N) belongs to M M A(R + ; R -)
, satisfies the equation of motion [START_REF] Moreau | Bounded variation in time[END_REF] and also the impact constitutive equation [START_REF] Rockafellar | Convex Analysis[END_REF]. Note, by the way, that this motion exhibits an infinite number of impacts on a compact time subinterval. It could easily be proved that no motion, defined on [0, ∞[, with finite number of impact on every compact interval can exist. Now, we are going to analyse what happens when the flow of time is reversed. Let define q by: q

[0, 4] → R - t → q(4 -t)
Considering the initial state (q 0 , v 0 ) = (0, 0) at t 0 = 0, it is easily seen that q satisfies both the equation of motion and the impact constitutive equation as soon as the impact function is replaced by φ ≡ 2. But, q ≡ 0 is also seen to satisfy the same initial condition, the equation of motion and the impact constitutive equation.

To eliminate this pathology of non-uniqueness, we are led to add the following hypothesis:

-H4: the kinetic energy of the system can not increase during an impact:

∀t, 1 2 q+ (t) 2 q(t) ≤ 1 2 q-(t) 2 q(t) . ( 19
)
Taking into account the impact constitutive equation ( 16), condition [START_REF] Schatzman | Uniqueness and continuous dependence on data for one dimensional impact problems[END_REF] can be rewritten as:

Proj q q-; V 2 + φ 2 Proj q q-; N 2 ≤ Proj q q-; V 2 + Proj q q-; N 2 ,
which implies φ ≤ 1 as soon as Proj q [ q-; N ] = 0. The final form of the impact constitutive equation is therefore:

∀t, q+ (t) = q-(t) -1 + φ q(t), q-(t) Proj q(t) q-(t); N (q(t)) ,
where the impact function φ is an arbitrary function from T Q to [0, 1]. The two extreme cases φ ≡ 0 and φ ≡ 1 are called respectively the completely inelastic and the elastic impact function.

Formulation of the evolution problem

In this subsection, the results of the previous subsections are brought together in order to formulate the resulting evolution problem which will be studied in the subsequent sections. We add an assumption on the regularity of the data: they are supposed to be real-analytic. This assumption will be motivated by the counterexamples of section 5. The precise mathematical setting is:

-Q is an analytic riemannian manifold of dimension d. -ϕ i (i = 1, 2, • • • , n) are n real analytic functions defined on Q. One defines: J(q) = {i ∈ {1, 2, • • • , n} ; ϕ i (q) ≥ 0} A = {q ∈ Q ; ∀i ∈ {1, 2, • • • , n}, ϕ i (q) ≤ 0} V (q) = {v ∈ T q Q ; ∀i ∈ J(q), dϕ i (q), v q ≤ 0} T A + = {(q, v) ∈ T Q ; q ∈ A and v ∈ V (q)} T A -= {(q, v) ∈ T Q ; q ∈ A and v ∈ -V (q)} N * (q) = n i=1 λ i dϕ i (q) ; ∀i ∈ J(q), λ i ≥ 0, ∀i ∈ J(q), λ i = 0 N (q) = n i=1 λ i ∇ϕ i (q) ; ∀i ∈ J(q), λ i ≥ 0, ∀i ∈ J(q), λ i = 0
The functions ϕ i are assumed to be functionally independent in the sense that:

∀q ∈ A, (dϕ i (q)) i∈J(q) is linearly independent in T * q Q (20)
-The impact function φ is an arbitrary function from T A -into [0, 1]. No regularity assumption is made on φ. -I is a real interval and O an open neighbourhood of T A + in T Q and the efforts mapping is supposed to be an analytic mapping from O×I into T * Q such that:

∀(q, v) ∈ O, ∀t ∈ I, Π * Q (f (q, v; t)) = Π Q (q, v) = q.
-We are given an initial time t 0 in I such that I contains a right neighboorhood of t 0 and an initial state (q 0 , v 0 ) in T A + .

According to the previous subsections, the evolution problem associated with the dynamics of discrete mechanical systems with perfect unilateral constraints can be formulated as:

Problem P: find T ∈ I ∪ {+∞} , T > t 0 and q ∈ M M A([t 0 , T [; Q) such that: • q(t 0 ), q+ (t 0 ) = (q 0 , v 0 ) (21) • ∀t ∈ [t 0 , T [ q(t), q+ (t) ∈ T A + (22) • R = D q+ -f q, q+ ; t dt ∈ -N * (q) for |R| -a.a. t ∈ [t 0 , T [ (23) • ∀t ∈]t 0 , T [, q+ = q--1 + φ(q, q-) Proj q q-; N (q) (24)
where equation ( 23) is to be understood in the sense of convention [START_REF] Moreau | Décomposition orthogonale d'un espace hilbertien selon deux cones mutuellement polaires[END_REF].

The existence and uniqueness of solutions for problem P will be studied in section 4. Before studying this question, let us state two almost obvious results.

Proposition 6 Any solution (T, q) of problem P satisfies:

-Supp R ⊂ {t ∈ [t 0 , T [; q(t) ∈ ∂A} -For all open subinterval J of [t 0 , T [ such that q(J) ⊂ • A , q |J is analytic and: D q(t) dt = f (q(t), q(t); t) , ∀t ∈ J.
Proof.

-

Let J be an open subinterval of [t 0 , T [ such that q(J) ⊂ • A . By equality (9), one has: ∀t ∈ J, N * (q(t)) = {0}.
As a consequence of relation ( 23) and convention [START_REF] Moreau | Décomposition orthogonale d'un espace hilbertien selon deux cones mutuellement polaires[END_REF], we get:

∀ϕ ∈ C 0 c J, q |J ; T Q , J ϕ(t), dR q(t) = 0, which is R |J = 0 or Supp R ⊂ [t 0 , T [\J.
The first item of proposition 6 follows.

-We have:

D q+ |J = f (q, q+ ; t) dt, which is, D q+ |J = • f (q, q+ ; t) dt.

There results from proposition 28 that q+

|J is locally absolutely continuous, and, therefore, ∀t ∈ J, q+ (t) = q-(t) = q(t), by proposition 32. We get

D q dt = D q+ dt = f (q, q; t), for dt-a.a. t ∈ J,
again by proposition 28. The conclusion follows by use of classical results on ordinary differential equations.

Proposition 7 (Energy inequality) Any solution (T, q) of problem P satisfies:

∀t 1 , t 2 ∈ [t 0 , T [, t 1 ≤ t 2 , K q(t 2 ), q+ (t 2 ) -K q(t 1 ), q+ (t 1 ) = 1 2 q+ (t 2 ) 2 q(t 2 ) - 1 2 q+ (t 1 ) 2 q(t 1 ) ≤ t 2 t 1 f q(s), q+ (s); s , q+ (s) q(s) ds
Proof. We have the following equality between real measures:

q+ (t) + q-(t) 2 , D q+ q(t) = q+ (t) + q-(t) 2 , f q(t), q+ (t); t q(t) dt + q+ (t) + q-(t) 2 , R q(t) .
Integrating over ]t 1 , t 2 ] and using propositions 30 and 32, one gets:

1 2 q+ (t 2 ) 2 q(t 2 ) - 1 2 q+ (t 1 ) 2 q(t 1 ) = ( 25
) ]t 1 ,t 2 ] q+ (t) + q-(t) 2 , f q(t), q+ (t); t q(t) dt + ]t 1 ,t 2 ] q+ + q- 2 , dR q . Consider D = t ∈]t 1 , t 2 ] ; q+ (t) + q-(t) 2 = q+ (t) .
D is (at most) countable and therefore Lebesgue-negligible. There results:

- D q+ (t) + q-(t) 2 , f q(t), q+ (t); t q(t) dt = 0 -Similarly: ]t 1 ,t 2 ]\D q+ (t) + q-(t) 2 , f q(t), q+ (t); t q(t) dt = t 2 t 1 q+ (t), f q(t), q+ (t); t q(t) dt
The dynamics of discrete mechanical systems with perfect unilateral constraints 19

-Let us denote by θ R the density of measure R with respect to its modulus measure |R| provided by proposition 26. Since:

∀t ∈]t 1 , t 2 ] \ D, q+ (t) + q-(t) 2 = q+ (t) = q-(t),
we get:

]t 1 ,t 2 ]\D q+ (t) + q-(t) 2 , dR q(t) = ]t 1 ,t 2 ]\D q+ (t), θ R (t) q(t) d |R| = ]t 1 ,t 2 ]\D q-(t), θ R (t) q(t) d |R| (26) But: θ R (t) ∈ -N * (q(t)) for |R| -a.a. t ∈]t 1 , t 2 ] \ D,
and therefore the second integral in equation ( 26) is nonnegative whereas the third is nonpositive since V (q(t) and N * (q(t)) are polar cones. As a consequence:

]t 1 ,t 2 ]\D q+ (t) + q-(t) 2 , dR q(t) = 0
-The following integral:

D q+ (t) + q-(t) 2 , dR q(t) = D q+ (t) + q-(t) 2 , D q+ q(t) = 1 2 t∈D q+ (t) 2 q(t) -q-(t) 2 q(t)
is nonpositive by virtue of hypothesis H4.

The proposition results from equation (25) and from the estimation of these four integrals.

Existence and uniqueness of solutions for problem P

This section is devoted to prove existence and uniqueness of a maximal solution for problem P. Sufficient conditions to ensure that this maximal solution is defined for all time are also given. More precisely, we are going to prove the following results.

Theorem 8 There is local existence and uniqueness of solution of problem P in the sense that:

there exists a solution (T, q) of problem P. Actually, there exists T > t 0 and an analytic function q

: [t 0 , T [→ Q which is a solution of problem P. -if (T 1 , q 1
) and (T 2 , q 2 ) are two solutions of problem P, then:

∃T, t 0 < T ≤ min{T 1 , T 2 }, q 1|[t 0 ,T [ = q 2|[t 0 ,T [
Then, a standart argument yields:

Corollary 9 Problem P admits a unique maximal solution (T m , q m ) (t 0 < T m ≤ +∞) in the sense that if (T, q) denotes an arbitrary solution of problem P, then:

T ≤ T m and q = q m|[t 0 ,T [ . Moreover, for each t ∈ [t 0 , T m [, there exists a right neighbourhood [t, t + η[ of t such that the restriction of q m to [t, t + η[ is analytic.
We shall say that the maximal solution of problem P is global if it is defined on

I ∩ [t 0 , +∞[.
Theorem 10 Assume that the configuration space Q is a complete riemannian manifold and that the efforts mapping f admits the estimate:

∀(q, v) ∈ T A + , for dt-a.a. t ∈ I ∩ [t 0 , +∞[, f (q, v; t) q ≤ l(t) 1 + d(q, q 0 ) + v q ,
where l(t) is a (necessarily nonnegative) function of L 1 loc (R; R). Then, the maximal solution of problem P is global.

Let us say a word of how the proof of these results is going to be structured. First, we construct T a > t 0 and an analytic function q a : [t 0 , T a [→ Q such that (T a , q a ) is a solution of problem P: this is the object of section 4.1. In section 4.2, we prove that if q ∈ M M A([t 0 , T [; Q) is any other solution, then q and q a coincide identically on a right neighbourhood of t 0 . This is the most difficult part to prove but it is also the crucial one. For the proof of theorem 10, we first notice that for q ∈ M M A([t 0 , T [; Q) (T finite) satisfying the equation of motion (23), boundedness of q+ implies finiteness of Var ( q+ ; [t 0 , T [): this is the object of proposition 18 of section 4.3. Note that the impact constitutive equation (24) plays no role in this property. Then, theorem 10 is deduced from the energy inequality (proposition 7) and Gronwall-Bellman lemma.

In the proof of these results, we shall use the following notations. If J is any subset of {1, 2, • • • , n}, Gram(J) will be the Gram matrix:

Gram(J) = ⎛ ⎜ ⎜ ⎝ . . . • • • (∇ϕ i (q 0 ), ∇ϕ j (q 0 )) q 0 • • • . . . ⎞ ⎟ ⎟ ⎠ i,j∈J
.

If x is an arbitrary element of R J whose components are x i with i ∈ J, then (x i ) i∈J will denote the columm matrix:

(x i ) i∈J = ⎛ ⎜ ⎜ ⎝ . . . x i . . . ⎞ ⎟ ⎟ ⎠ i∈J ,
and T (x i ) i∈J the associated row matrix:

T (x i ) i∈J = • • • x i • • • i∈J ,
The dynamics of discrete mechanical systems with perfect unilateral constraints 21

Proof of local existence

Local existence is rather easy to prove in the setting of analytic data. The proof is a little bit lengthy but involves no specific difficulty. We begin by technical lemmas.

Let X(t) be a C ∞ vector field along a C ∞ curve q(t) on Q. The covariant derivative D dt X(t) of X along q defines a C ∞ vector field along q. So, one may consider its covariant derivative along q which will be denoted by D 2 dt 2 X(t). By induction, we get the definition of D i dt i X(t) (i ∈ N * ). One has: Lemma 11 Let X be a C ∞ vector field on Q and q I , q II two C ∞ curves on Q. m being a nonnegative integer, one assumes:

q I (t 0 ) = q II (t 0 ), qI (t 0 ) = qII (t 0 ),
and:

∀i ∈ {1, 2, • • • , m} , D i dt i qI (t 0 ) = D i dt i qII (t 0 ) Then, ∀i ∈ {1, 2, • • • , m + 1} , D i dt i X(q I (t 0 )) = D i dt i X(q II (t 0 ))
Proof. Consider a local chart at q I (t 0 ) = q II (t 0 ). If q(t) is either q I (t) or q II (t): q(t) = qi (t)e i (q(t)), X(q(t)) = X i (q(t))e i (q(t)), D dt X(q(t)) = ∇X i (q(t)), q(t) q(t) + Γ i jk (q(t))X j (q(t)) qk (t) e i (q(t)).

Then,

D 2 dt 2 X(q(t)) = D dt ∇X i (q(t)), q(t) q(t) + ∇X i (q(t)), D dt q(t) q(t) + ∇Γ i jk (q(t)), q(t) q(t) X j (q(t)) qk (t) + Γ i jk (q(t)) ∇X j (q(t)), q(t) q(t) qk (t) + Γ i jk (q(t))X j (q(t)) D q(t) dt k -Γ k lm (q(t)) ql (t) qm (t) + Γ i jk (q(t)) DX(q(t)) dt j qk (t) e i (q(t)),
which gives the desired conclusion for the case m = 1. For arbitrary m, an easy induction based on the same type of computation in a local chart shows the existence of functions h i : (T Q) i-1 → T Q independent on the considered curve q(t) and such that:

D i X(q(t)) dt i = h i q(t), q(t), D q(t) dt , • • • , D i-1 q(t) dt i-1 .
Exactly the same technique applies to prove

Lemma 12 Let X : T Q × I → T Q a C ∞ mapping such that: Π Q (X(q, v; t)) = Π Q (q, v) = q
, where I denotes a real interval containing t 0 . Let m be an arbitrary nonnegative integer and q I , q II two C ∞ curves on Q such that:

q I (t 0 ) = q II (t 0 ), qI (t 0 ) = qII (t 0 ), and:

∀i ∈ {1, 2, • • • , m} , D i dt i qI (t 0 ) = D i dt i qII (t 0 ) Then, ∀i ∈ {1, 2, • • • , m} , D i dt i X q I (t 0 ), qI (t 0 ); t 0 = D i dt i X q II (t 0 ), qII (t 0 ); t 0
Lemma 13 Consider (q 0 , v 0 ) ∈ T A + and J ⊂ J(q 0 ) an arbitrary subset of: {i ∈ J(q 0 ); dϕ i (q 0 ), v 0 q 0 = 0} .

We denote by q u and q c some local solutions of problems:

E u ⎧ ⎨ ⎩ D qu dt = f (q u , qu ; t) (q u (t 0 ), qu (t 0 )) = (q 0 , v 0 ) E c ⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ D qc dt = f (q c , qc ; t) + i∈J(q 0 ) λ i (t)dϕ i (q c ) ∀i ∈ J ϕ i (q) ≡ 0 ∀i ∈ J(q 0 ) \ J λ i (t) ≡ 0 (q c (t 0 ), qc (t 0 )) = (q 0 , v 0 )
furnished respectively by Cauchy-Lipschitz theorem and proposition 3. Then:

Gram(J(q 0 )) (λ i (t 0 )) i∈J(q 0 ) = d 2 dt 2 ϕ i (q c (t 0 )) - d 2 dt 2 ϕ i (q u (t 0 )) i∈J(q 0 )
Morover, if:

∃m ∈ N * , ∀i = 0, 1, • • • , m -1, ∀j ∈ J(q 0 ), d i dt i λ j (t 0 ) = 0,
then:

Gram(J(q 0 )) d m dt m λ i (t 0 ) i∈J(q 0 ) = d m+2 dt m+2 ϕ i (q c (t 0 )) - d m+2 dt m+2 ϕ i (q u (t 0 )) i∈J(q 0 )
The dynamics of discrete mechanical systems with perfect unilateral constraints 23

Proof. First, from:

(q u (t 0 ), qu (t 0 )) = (q c (t 0 ), qc (t 0 )) = (q 0 , v 0 ), it follows: ∀i ∈ J(q 0 ), D dt ∇ϕ i (q u (t 0 )) = D dt ∇ϕ i (q c (t 0 )),
on one hand, and:

D dt qu (t 0 ) - D dt qc (t 0 ) = - i∈J(q 0 ) λ i (t 0 )∇ϕ i (q 0 ),
on the other hand. Therefore:

∀i ∈ J(q 0 ), d 2 dt 2 ϕ i (q c (t 0 )) - d 2 dt 2 ϕ i (q u (t 0 )) = D dt ∇ϕ i (q c (t 0 )), v 0 q 0 + ∇ϕ i (q c (t 0 )), D dt qc (t 0 ) q 0 - D dt ∇ϕ i (q u (t 0 )), v 0 q 0 -∇ϕ i (q u (t 0 )), D dt qu (t 0 ) q 0 , = j∈J(q 0 ) λ j (t 0 ) (∇ϕ i (q 0 ), ∇ϕ j (q 0 )) q 0 ,
which is the announced result. Second, assume:

∀j ∈ J(q 0 ), ∀i = 0, 1, • • • , m -1, d i dt i λ j (t 0 ) = 0.
An easy induction based on lemmas 11 and 12 gives:

∀i = 1, 2, • • • , m, D i dt i qu (t 0 ) = D i dt i qc (t 0 ), D m+1 dt m+1 qu (t 0 ) = D m+1 dt m+1 qc (t 0 ) - j∈J(q 0 ) d m dt m λ j (t 0 )∇ϕ j (q 0 ),
and,

∀j ∈ J(q 0 ), ∀i = 1, 2, • • • , m + 1, D i dt i ∇ϕ j (q u (t 0 )) = D i dt i ∇ϕ j (q c (t 0 )).
Therefore:

∀i ∈ J(q 0 ), d m+2 dt m+2 ϕ i (q c (t 0 )) - d m+2 dt m+2 ϕ i (q u (t 0 )) = j∈J(q 0 ) d m dt m λ j (t 0 ) (∇ϕ i (q 0 ), ∇ϕ j (q 0 )) q 0 .
Proposition 14 Considering the data of problem P, we denote by P the following evolution problem. Problem P : find T ∈ I (T > t 0 ), an analytic curve q : [t 0 , T [→ Q and n analytic functions

λ i : [t 0 , T [→ R such that: ⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ ∀t ∈ [t 0 , T [, D q(t) dt = f (q(t), q(t); t) + n i=1 λ i (t)dϕ i (q(t)), ∀t ∈ [t 0 , T [, ∀i = 1, 2, • • • , n, λ i (t) ≤ 0, ϕ i (q(t)) ≤ 0, λ i (t)ϕ i (q(t)) = 0 (q(t 0 ), q(t 0 )) = (q 0 , v 0 )
Then, problem P admits a solution (T, q, λ 1 , • • • , λ n ) unique in the sense that any other solution is either a restriction or an analytic extension of (T, q, λ 1 ,

• • • , λ n ).
Proof. First, let us precise once for all that the meaning of an analytic function on a non necessarily open set S is that there is an analytic extension to an open set O containing S.

Step 1. Construction of some functions q and λ i . Define:

J 0 = {i ∈ {1, 2, • • • , n} ; ϕ i (q 0 ) = 0 and dϕ i (q 0 ), v 0 q 0 = 0} ,
and I 0 = K 0 = ∅. We denote by q (1) a solution of the cauchy problem:

C (1)
⎧ ⎨ ⎩ D q(t) dt = f (q(t), q(t); t), (q(t 0 ), q(t 0 )) = (q 0 , v 0 ).

Define:

C (1) = (λ * i ) ∈ R J 0 ; ∀i ∈ J 0 , λ * i ≤ 0 and ∀i ∈ K 0 , λ * i = 0 = (R -) J 0 C (1) = (μ * i ) ∈ R J 0 ; ∀i ∈ I 0 , μ * i = 0 and ∀i ∈ J 0 , μ * i ≤ 0 = (R -) J 0 Let λ (1) i i∈J 0
∈ C (1) be the solution of the variational inequality:

∀ (λ * i ) i∈J 0 ∈ C (1) , T λ (1) i i∈J 0 Gram(J 0 ) λ * i -λ (1) i i∈J 0 ≥ T - d 2 dt 2 ϕ i (q (1) (t 0 )) i∈J 0 λ * i -λ (1) i i∈J 0
furnished by Lions-Stampacchia theorem (see [START_REF] Lions | Variational inequalities[END_REF]). Let μ

(1) i i∈J 0 ∈ C (1) be defined by:

μ (1) i i∈J 0 = Gram(J 0 ) λ (1) i i∈J 0 + d 2 dt 2 ϕ i (q (1) (t 0 )) i∈J 0 ( 27
)
The dynamics of discrete mechanical systems with perfect unilateral constraints 25

and I 1 , J 1 , K 1 by:

I 1 = I 0 ∪ i ∈ J 0 ; λ (1) i < 0 and μ (1) i = 0 , J 1 = i ∈ J 0 ; λ (1) i = 0 and μ (1) i = 0 , K 1 = K 0 ∪ i ∈ J 0 ; λ (1) i = 0 and μ (1) i < 0 , Now suppose q (n) , λ (n) i , μ (n) i
, I n , J n and K n are constructed. Then, q (n+1) is defined to be a local solution of the Cauchy problem:

C (n+1) ⎧ ⎪ ⎨ ⎪ ⎩ D q(t) dt = f (q(t), q(t); t) + j∈J 0 n i=1 λ (i) j (t -t 0 ) j-1 (j -1)! dϕ j (q(t)),
(q(t 0 ), q(t 0 )) = (q 0 , v 0 ).

C (n+1) = (λ * i ) ∈ R J 0 ; ∀i ∈ J n , λ * i ≤ 0, and ∀i ∈ K n , λ * i = 0 , C (n+1) = (μ * i ) ∈ R J 0 ; ∀i ∈ I n , μ * i = 0, and ∀i ∈ J n , μ * i ≤ 0 . λ (n+1) i i∈J 0 ∈ C (n+1
) is defined to be the solution of the variational inequality:

∀ (λ * i ) i∈J 0 ∈ C (n+1) , T λ (n+1) j i∈J 0 Gram(J 0 ) λ * i -λ (n+1) i i∈J 0 ≥ T - d n+2 dt n+2 ϕ i (q (n+1) (t 0 )) i∈J 0 λ * i -λ (n+1) i i∈J 0 μ (n+1) i i∈J 0 ∈ C (n+1
) is defined by:

μ (n+1) i i∈J 0 = Gram(J 0 ) λ (n+1) i i∈J 0 + d n+2 dt n+2 ϕ i (q (n+1) (t 0 )) i∈J 0
and I n+1 , J n+1 , K n+1 by:

I n+1 = I n ∪ i ∈ J n ; λ (n+1) i < 0 and μ (n+1) i = 0 , J n+1 = i ∈ J n ; λ (n+1) i = 0 and μ (n+1) i = 0 , K n+1 = K n ∪ i ∈ J n ; λ (n+1) i = 0 and μ (n+1) i < 0 , Thus, the sequences q (n) , λ (n) i i∈J 0 , μ (n) i i∈J 0
, I n , J n and K n are defined by induction for n ∈ N * and for all n in N * , I n , J n , K n is a partition of J 0 . Moreover, one has:

∀n ∈ N, I n ⊂ I n+1 , J n+1 ⊂ J n , K n ⊂ K n+1 .
Define:

I = ∞ n=0 I n , J = ∞ n=0 J n , K = ∞ n=0 K n .
It is readily seen that I, J, K form a partition of J 0 . We denote by (q, (λ i ) i∈I ) a local solution of the evolution problem:

C ⎧ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎩ D q(t) dt = f (q(t), q(t); t) + i∈I λ i (t)dϕ i (q(t)), ∀i ∈ I, ϕ i (q(t)) ≡ 0, (q(t 0 ), q(t 0 )) = (q 0 , v 0 ),
furnished by proposition 3. The functions q and λ i are analytic. For any i in {1, 2, • • • , n} \ I, the functions λ i are defined to be the identically vanishing function:

∀i ∈ {1, 2, • • • , n} \ I, λ i ≡ 0.
Step 2. We have:

∀j ∈ J 0 , ∀i ∈ N, d i dt i λ j (t 0 ) = λ (i+1) j , ∀j ∈ J 0 , ∀i ∈ N, d i+2 dt i+2 ϕ j (q(t 0 )) = μ (i+1) j .
Indeed, applying lemma 13 to Cauchy problems C (1) and C yields, thanks to equation (27),

μ (1) j - d 2 dt 2 ϕ j (q(t 0 )) j∈J 0 = Gram(J 0 ) λ (1) j -λ j (t 0 ) j∈J 0
But, by definition of I,

I 1 ⊂ I ⊂ J 0 \ K 1 ,
and, so, ∀j ∈ I, μ

(1)

j = d 2 dt 2 ϕ j (q(t 0 )) = 0, ∀j ∈ J 0 \ I, λ (1)
j = λ j (t 0 ) = 0. Therefore, T λ (1) j -λ j (t 0 ) j∈J 0 Gram(J 0 ) λ (1) j -λ j (t 0 ) j∈J 0 = 0,
and the conclusion follows for i = 0, since the Gram matrix is positive definite. For i ≥ 1, we only have to apply successively lemma 13 to Cauchy problems C (i+1) and C.

Step 3. The functions q and λ i define a solution of problem P . Indeed, by construction of the real numbers λ (j) i and μ (j) i and by step 2, we have:

∀i ∈ I, ∃n i ∈ N, d n i dt n i λ i (t 0 ) < 0 and ∀n < n i , d n dt n λ i (t 0 ) = 0,
and,

∀i ∈ K, ∃n i ≥ 2, d n i dt n i ϕ i (q(t 0 )) < 0 and ∀n < n i , d n dt n ϕ i (q(t 0 )) = 0, ∀i ∈ J 0 \ K, ∀n ∈ N, d n dt n ϕ i (q(t 0 )) = 0.
Each function λ i (t) and ϕ i (q(t)) being real-analytic, there results:

∃α > 0, ∀t ∈ [t 0 , t 0 + α[, ∀i ∈ J 0 , λ i (t) ≤ 0, and ϕ i (q(t)) ≤ 0.
Actually, α > 0 is assumed to be sufficiently small to ensure:

∀i ∈ {1, 2, • • • , n} \ J 0 , ∀t ∈]t 0 , t 0 + α[, ϕ i (q(t)) < 0,
which is possible simply by continuity. Now, it is easily seen that t 0 + α, q, (λ i ) i∈{1,2,•••,n} defines a solution of problem P .

Step 4. Uniqueness part of the proposition. By Cauchy-Lipshitz theorem, q is uniquely determined by the functions λ j (j = 1, 2, • • • , n). Being analytic, these functions λ j are uniquely determined by the collection of real numbers

d i λ j (t 0 ) dt i , (i ∈ N, j ∈ {1, 2, • • • , n}).
Therefore, to prove uniqueness, one has only to show that these real numbers are determined by the data of the evolution problem.

Consider an arbitrary analytic solution (T, q, λ 1 , • • • , λ n ) of problem P . A repeated use of lemma 13, similar to the one of step 2 yields:

∀j ∈ J 0 , ∀i ∈ N, d i dt i λ j (t 0 ) = λ (i+1) j . Moreover, ∀j ∈ {1, 2, • • • , n} \ J 0 , ∀i ∈ N, d i dt i λ j (t 0 ) = 0,
and the conclusion follows.

Proof of the local existence part of theorem 8.

Let (T a , q a , λ 1 a , • • • , λ n a )
be an analytic solution of problem P . It is readily seen that (T a , q a ) is a local solution of problem P.

Proof of local uniqueness

Local uniqueness is the most difficult part of theorem 8. First, we recall a standart result:

Lemma 15 (Gronwall-Bellman) Consider two functions m 1 ∈ BV ([0, T ]; R) and m 2 ∈ L 1 (0, T ; R) such that: for a.a. t ∈]0, T [, m 2 (t) ≥ 0. Let φ ∈ BV ([0, T ]; R) such that: ∀t ∈ [0, T ], φ(t) ≤ m 1 (t) + t 0 m 2 (s)φ(s) ds Then, ∀t ∈ [0, T ], φ(t) ≤ m 1 (t) + t 0 m 1 (s)m 2 (s)e t s m 2 (σ) dσ ds
We have the following corollary of Gronwall-Bellman lemma :

Lemma 16 Let m be a nonnegative integer, and ψ : [0, T ] → R an integrable function. If φ : [0, T ] → R is any absolutely continuous function such that φ(t) = o(t m+1 ) when t tends towards 0 and such that there exists a nonnegative real constant C such that:

for dt-a.a. t ∈]0, T [, t d dt φ(t) ≤ (1 + m + Ct) φ(t) + t m+2 ψ(t), then, ∀t ∈ [0, T ], φ(t) ≤ t m+1 e Ct t 0 ψ(s)e -Cs ds.
Proof. This is almost obvious. Dividing each member of the inequality by t m+2 , we obtain:

for dt-a.a. t ∈]0, T [, d dt φ(t) t m+1 ≤ C φ(t) t m+1 + ψ(t).
After integration, Gronwall-Bellman lemma yields:

∀t ∈]0, T ], φ(t) t m+1 ≤ t 0 ψ(s) ds + t 0 Ce C(t-s) s 0 ψ(σ) dσ ds.
Then, an integration by part gives the desired conclusion.

Proof of the local uniqueness part of theorem 8. Consider, on one hand, the analytic solution (T a , q a , λ 1 a , • • • , λ n a ) of problem P supplied by proposition 14, and on the other hand, an arbitrary solution (T, q) of problem P. We have to prove that q and q a identically coincide on a right neighbourhood of t 0 .

Step 1. Parametrization of the problem and notations.

Consider a local chart ψ : U ⊂ Q → R d on Q centered at q 0 such that the cardJ(q 0 ) first components of ψ(q) are (ϕ i (q)) i∈J(q 0 ) . Recall that such a chart exists since (dϕ i (q 0 )) i∈J(q 0 ) is linearly independent in T * q 0 Q. We choose α > 0, sufficiently small to have:

• ∀t ∈ [t 0 , t 0 + α], q a (t) ∈ U, q(t) ∈ U, • ∀i ∈ J(q 0 ), ∀t ∈ [t 0 , t 0 + α], d dt ϕ i (q a (t)) = dϕ i (q a (t)), qa (t) q a (t) ≤ 0 • ∀i ∈ {1, 2, • • • , n} \ J(q 0 ), ∀t ∈ [t 0 , t 0 + α], ϕ i (q a (t)) < 0, ϕ i (q(t)) < 0.
(28) Such a choice for α is possible because:

the functions q a (t) and ϕ i (q a (t)) are real analytic, the functions q(t) and ϕ i (q(t)) are continuous.

We denote by f i the components of f in the natural basis (e i ) associated with the chart under consideration. Since q a is an analytic local solution of problem P, we have:

∀i ∈ {1, 2, • • • , d} , ∀s ∈ [t 0 , t 0 + α], g ij (q a ) qj a + Γ j kl (q a ) qk a ql a -f i (q a , qa ; s) = λ i a (s), (29)
after appropriate renumbering of the functions λ i a . In the sequel, d 0 will stand for cardJ(q 0 ). There results from these choices that: ∀i > d 0 , λ i a ≡ 0. We denote by |.| the standart euclidean norm on R d . Confusing (abusively) q and ψ(q), we shall write:

|q| 2 = d i=1 q i 2 ,
and:

q+ 2 = d i=1 q+i 2 .
Step 2. There exists some positive real constants C 1 and C 2 such that the following estimate:

∀t ∈ [t 0 , t 0 + α], t t 0 |q -q a | 2 (s) + q+ -qa 2 (s) ds ≤ - 1 C 1 t t 0 e C 2 (t-s) s t 0 d 0 i=1 λ i a (σ) q+i (σ) dσ ds. ( 30
)
holds.

To prove this assertion, we first write the equation of motion (23) in the chart under consideration using proposition 29:

∀i ∈ {1, 2, • • • , d}, ∀t ∈ [t 0 , t 0 + α], g ij (q) d q+j + Γ j kl (q) q+k q+l dt = f i (q, q+ ; t) dt + d 0 j=1 δ ij μ j ,
where the μ j are nonpositive real measures. But, by propositions 29 and 30,

d 1 2 q+i -qi a g ij (q) q+j -qj a = q-i + q+i 2 -qi a g ij (q) d q+j -qj a dt + Γ j kl (q) q+k q+l -ql a dt .
Therefore,

d 1 2 q+i -qi a g ij (q) q+j -qj a = q+i -qi a f i (q, q+ ; t)dt -q+i -qi a g ij (q) qj a + Γ j kl (q) q+k ql a dt + d 0 j=1 q-j + q+j 2 -qj a μ j . But, ∀j ∈ {1, 2, • • • , d 0 } , ∃i ∈ J(q 0 ), ∀t ∈ [t 0 , t 0 + α], qj a (t) = d dt ϕ i (q a (t)) ≤ 0,
by formulae (28), and,

d 0 j=1 q-j + q+j 2 μ j = q-+ q+ 2 , R q ,
which is a nonpositive real measure by proposition 7. Therefore,

d 1 2 q+i -qi a g ij (q) q+j -qj a ≤ q+i -qi a f i (q, q+ ; t) -q+i -qi a g ij (q) qj a + Γ j kl (q) q+k ql a dt,
in the sense of ordering of real measures. Integrating over ]t 0 , t] (t ∈ [t 0 , t 0 + α]), we get:

1 2 q+i -qi a g ij (q) q+j -qj a ≤ t t 0 q+i -qi a f i (q, q+ ; s) -q+i -qi a g ij (q) qj a + Γ j kl (q) q+k ql a ds.
The term within the integral sign is an analytic function of the three variables q, q+ and s. Therefore, it is also an analytic function of the three variables qq a , q+ -qa and s. It is written under the form:

q+i -qi a F i (qq a , q+ -qa ; s).

The dynamics of discrete mechanical systems with perfect unilateral constraints 31 But, each function F i can be decomposed under the form:

F i (q -q a , q+ -qa ; s) = F i (0, 0; s) + G i (q -q a , q+ -qa ; s),
where the G i are analytic and G i (0, 0; s) ≡ 0. Hence, there exists d positive constants M i such that:

∀t ∈ [t 0 , t 0 + α], G i (q(s) -q a (s), q+ (s) -qa (s); s) ≤ M i |q(s) -q a (s)| 2 + | q+ (s) -qa (s)| 2
Defining M to be the maximum of the constants M i , we have proved:

∀t ∈ [t 0 , t 0 + α], 1 2 q+i -qi a g ij (q) q+j -qj a ≤ t t 0 q+i -qi a f i (q a , qa ; s) -g ij (q a ) qj a + Γ j kl (q a ) qk a ql a +M d q+ -qa |q -q a | 2 + | q+ -qa | 2 ds.
Moreover, by a compactness argument:

∃C 1 > 0, ∀t ∈ [t 0 , t 0 + α], 1 2 q+i -qi a g ij (q) q+j -qj a ≥ C 1 q+ -qa 2 ,
and therefore:

∀t ∈ [t 0 , t 0 + α], q+ (t) -qa (t) 2 ≤ 1 C 1 t t 0 q+i -qi a f i (q a , qa ; s) -g ij (q a ) qj a + Γ j kl (q a ) qk a ql a ds + M d C 1 t t 0 q+ -qa |q -q a | 2 + | q+ -qa | 2 ds.
Moreover, by use of Cauchy-Schwartz inequality:

∀t ∈ [t 0 , t 0 + α], |q(t) -q a (t)| 2 ≤ α t t 0 q+ (s) -qa (s) 2 ds.
We obtain:

∀t ∈ [t 0 , t 0 + α], |q -q a | 2 (t) + q+ -qa 2 (t) ≤ M d C 1 + α t t 0 |q -q a | 2 (s) + q+ -qa 2 (s) ds - 1 C 1 t t 0 d 0 i=1 λ i a (s) q+i -qi a ds, ( 31
)
where formulae (29) have been used. We define:

C 2 = M d C 1 + α.
Notice that, actually:

∀i ∈ {1, 2, • • • , d 0 } , λ i a q i a ≡ 0,
and, so, by the analyticity of functions q i a and λ i a :

∀i ∈ {1, 2, • • • , d 0 } , λ i a qi a ≡ 0.
Multiplying both terms of inequality ( 31) by e -C 2 t and integrating, we get:

∀t ∈ [t 0 , t 0 + α], t t 0 |q -q a | 2 (s) + q+ -qa 2 (s) ds ≤ - 1 C 1 t t 0 e C 2 (t-s) s t 0 d 0 i=1 λ i a (σ) q+i (σ) dσ ds,
which is nothing but estimate (30).

Step 3. Estimate (30) implies that the function t → d 0 i=1 λ i a (t) q+i (t) vanishes identically on a right neighbourhood of t 0 Indeed, by estimate (30):

∀t ∈ [t 0 , t 0 + α], t t 0 e -C 2 s s t 0 d 0 i=1 λ i a (σ) q+i (σ) dσ ds ≤ 0,
which is, after integration by part:

∀t ∈ [t 0 , t 0 + α], t t 0 e -C 2 s d i=1 λ i a (s)q i (s) ds ≤ t t 0 e -C 2 s s t 0 d 0 i=1 q i (σ) λi a (σ) dσ ds. (32) But, since, ∀i ∈ {1, 2, • • • , d 0 } , ∀s ∈ [t 0 , t 0 + α], λ i a (s) ≤ 0 and q i (s) ≤ 0,
the two members of inequality (32) are nonnegative and, therefore, the inequality is preserved when taking the absolute value of each member. We get:

∀t ∈ [t 0 , t 0 + α], t t 0 e -C 2 s d 0 i=1 λ i a (s)q i (s) ds ≤ t t 0 e -C 2 s s t 0 d 0 i=1 q i (σ) λi a (σ) dσ ds, ≤ t t 0 s t 0 e -C 2 σ d 0 i=1 q i (σ) λi a (σ) dσ ds.
We define:

Q i (s) = -e -C 2 (s+t 0 ) q i (s + t 0 ), L i (s) = -λ i a (s + t 0 )
. With these notations, we obtain:

∀t ∈ [0, α], t 0 d 0 i=1 L i (s)Q i (s) ds ≤ t 0 s 0 d 0 i=1 Li (s) Q i (s) dσ ds, ( 33
)
where the L i are nonnegative real-analytic functions and the Q i are nonnegative continuous functions which all vanish at t = 0 and which are differentiable at the origin. We are going to prove that inequality (33) implies that:

∃β ∈]0, α], ∀t ∈ [0, α], ∀i ∈ {1, 2, • • • , d 0 } , L i (t)Q i (t) = 0.
The functions L i being nonnegative real-analytic, there exists nonnegative integers

n 1 < n 2 < • • • < n m , a partition I 1 , I 2 , • • • , I m of {1, 2, • • • , d 0 }
, and nonnegative real-analytic functions G i such that:

∀k ∈ {1, 2, • • • , m}, ∀i ∈ I k , L i (s) = s n k G i (s),
with either G i (0) > 0 or G i ≡ 0. Inequality (33) may be rewritten as:

∀t ∈ [0, α], t 0 m k=1 i∈I k σ n k G i (σ)Q i (σ) dσ ≤ t 0 s 0 m k=1 i∈I k n k σ n k -1 G i (σ)Q i (σ) dσ ds + t 0 s 0 m k=1 i∈I k σ n k Ġi (σ) Q i (σ) dσ ds.
But, by the analyticity of the functions G i :

∃β > 0, ∃N > 0, ∀i ∈ J(q 0 ), ∀σ ∈ [0, β], Ġi (σ) ≤ N G i (σ).
Therefore,

∀t ∈ [0, β], t 0 m k=1 i∈I k σ n k G i (σ)Q i (σ) dσ ≤ t 0 s 0 m k=1 i∈I k n k σ n k -1 G i (σ)Q i (σ) dσ ds + N t t 0 s 0 m k=1 i∈I k σ n k -1 G i (σ)Q i (σ) dσ ds.
Integrating by part the left member of the inequality, we obtain:

∀t ∈ [0, β], t t 0 m k=1 i∈I k σ n k -1 G i (σ)Q i (σ) dσ ≤ t 0 s 0 m k=1 i∈I k (n k + 1)σ n k -1 G i (σ)Q i (σ) dσ ds + N t t 0 s 0 m k=1 i∈I k σ n k -1 G i (σ)Q i (σ) dσ ds. ( 34
)
Since each function

G i (σ)Q i (σ)/σ is bounded over [0, β],
there exists a nonnegative real constant H such that:

∀k ∈ {1, 2, • • • , m}, ∀t ∈ [0, β], t 0 s 0 i∈I k σ n k -1 G i (σ)Q i (σ) dσ ds ≤ Ht n k +2 .
Inequality (34) gives:

∀t ∈ [0, β], t t 0 i∈I 1 σ n 1 -1 G i (σ)Q i (σ) dσ ≤ (1 + n 1 + N t) t 0 s 0 i∈I 1 σ n 1 -1 G i (σ)Q i (σ) dσ ds + H 1 t n 2 +2 ,
where H 1 is a non negative real constant. As a consequence of lemma 16, we obtain:

t 0 s 0 i∈I 1 σ n 1 -1 G i (σ)Q i (σ) dσ ds = O(t n 2 +2
).

Coming back to inequality (34), we get:

∀t ∈ [0, β], t t 0 2 k=1 i∈I k σ n k -1 G i (σ)Q i (σ) dσ ≤ (1 + n 2 + N t) t 0 s 0 2 k=1 i∈I k σ n k -1 G i (σ)Q i (σ) dσ ds + H 2 t n 3 +2 .
Applying once more lemma 16, we obtain:

t 0 s 0 2 k=1 i∈I k σ n k -1 G i (σ)Q i (σ) dσ ds = O(t n 3 +2 ).
Proceeding inductively, we obtain:

t 0 s 0 m-1 k=1 i∈I k σ n k -1 G i (σ)Q i (σ) dσ ds = O(t n m +2
).

The dynamics of discrete mechanical systems with perfect unilateral constraints 35

But, by inequality (34),

∀t ∈ [0, β], t t 0 m k=1 i∈I k σ n k -1 G i (σ)Q i (σ) dσ ≤ (1 + n m + N t) t 0 s 0 m k=1 i∈I k σ n k -1 G i (σ)Q i (σ) dσ ds.
Using a last time lemma 16, we get:

∀t ∈ [0, β], t 0 s 0 m k=1 i∈I k σ n k -1 G i (σ)Q i (σ) dσ ds = 0, which implies: ∀i ∈ {1, 2, • • • , d 0 } , ∀t ∈ [0, β], G i (σ)Q i (σ) = 0,
which is nothing but:

∀i ∈ {1, 2, • • • , d 0 } , ∀t ∈ [t 0 , t 0 + β], λ i a (σ)q i (σ) = 0. But, the analyticity of the functions λ i a implies: ∀i ∈ {1, 2, • • • , d 0 } , ∀t ∈ [t 0 , t 0 + β],
λ i a (σ) q+i (σ) = 0, and the assertion of step 3 is proved.

Step 4.

Conclusion of the proof of local uniqueness.

Bringing together the results of steps 2 and 3, we get:

∀t ∈ [t 0 , t 0 + β], t t 0 |q -q a | 2 (s) + q+ -qa 2 (s) ds ≤ 0,
which gives the desired conclusion:

∀t ∈ [t 0 , t 0 + β], q(t) = q a (t).

Global solutions: proof of theorem 10

First, we recall a classical lemma whose proof may be found for example in [START_REF] Brezis | Opérateurs Maximaux Monotones et Semi-groupes de Contrations dans les Espaces de Hilbert[END_REF], p. 157.

Lemma 17 Let m be in L 1 (0, T ; R) such that m(t) ≥ 0 for almost all t in ]0, T [ and a be a real nonnegative constant. Consider φ ∈ BV ([0, T ]; R) such that:

∀t ∈ [0, T ], 1 2 φ 2 (t) ≤ 1 2 a 2 + t 0 m(s)φ(s) ds, then: ∀t ∈ [0, T ], |φ(t)| ≤ a + t 0 m(s) ds.

Proposition 18

The riemannian manifold Q is assumed to be complete. Let (T, q) be a solution of problem P such that:

-T ∈

• I (and ,in particular, T = +∞), -q+ (t) q(t) is bounded:

∃V m , ∀t ∈ [t 0 , T [, q+ (t) q(t) ≤ V m ,
then q+ has bounded variation over [t 0 , T [:

Var q+ ; [t 0 , T [ < ∞.
Proof. We denote by d the distance function associated with the metric space Q. Since,

-∀s 1 , s 2 ∈ [t 0 , T [, s 1 ≤ s 2 , d(q(s 1 ), q(s 2 )) ≤ s 2 s 1 q+ (σ) q(σ) dσ, -∀σ ∈ [t 0 , T [, q+ (σ) q(σ) ≤ V m , -Q is complete,
we deduce that lim t→T -q(t) exists in Q. It is denoted by:

q T = lim t→T -q(t).
Let (U, ψ) a local chart at q T on Q such that the cardJ(q T ) first components of ψ(q) in R d are (ϕ i (q)) i∈J(q T ) . Consider a compact neighbourhood K of q T in Q such that:

-K ⊂ U, -∀q ∈ K, J(q) ⊂ J(q T ).

One defines:

t 0 = min {t ∈ [t 0 , T [ ; ∀s ∈ [t, T [, q(s) ∈ K} .
Since [t 0 , t 0 ] is compact, one has:

Var q+ ; [t 0 , t 0 ] < ∞,
therefore, it remains only to prove:

Var q+ ; ]t 0 , T [ < ∞.
λ max (resp. λ min ) will denote the maximum (resp. the minimum) of the greatest (resp. least) eigenvalue of the matrix (g ij (q)) i,j=1,2,•••,d when q wanders in K.

With these notations, one obtains immediately:

∀i ∈ {1, 2, • • • , d} , ∀t ∈ [t 0 , T [, g ij (q(t)) q+j (t) ≤ √ λ max V m , q+i (t) ≤ V m √ λ min . ( 35
)
We denote by B q (0, V m ) the closed ball of T q Q with radius V m and centered at the origin. Considering the following compact subset K of T Q:

K = q∈K B q (0, V m ),
we define the following nonnegative real constant:

F = max i∈{1,2,•••,d} max (q,v;t)∈K ×[t 0 ,T ] |f i (q, v; t)| ,
and:

G = max i,j,k∈{1,2,•••,d} max q∈K ∂g ij (q) ∂q k .
Writing inclusion (23) in the local chart (U, ψ), we obtain:

∀i ∈ {1, 2, • • • , d} , g ij (q) d q+j + Γ j kl (q) q+k q+l dt = f i (q, q+ ; t) dt + λ i ,
where the λ i are d nonpositive real measures on ]t 0 ; T [. Expressing the Christoffel symbols in terms of the metric, one has:

∀i ∈ {1, 2, • • • , d} , (36)
g ij (q)d q+j + ∂g ij (q) ∂q k q+j q+k dt - 1 2 ∂g kl (q) ∂q i q+k q+l dt = f i (q, q+ ; t) dt + λ i , or, equivalently, ∀i ∈ {1, 2, • • • , d} , ( 37
)
d g ij (q) q+j = 1 2
∂g kl (q) ∂q i q+k q+l dt + f i (q, q+ ; t) dt + λ i .

One deduces:

∀i ∈ {1, 2, • • • , d} , ∀s 1 , s 2 ∈ [t 0 , T [, s 1 < s 2 , ]s 1 ,s 2 ] (-λ i ) = g ij (q(s 1 )) q+j (s 1 ) -g ij (q(s 2 )) q+j (s 2 ) + s 2 s 1 f i (q, q+ ; t) + 1 2 ∂g kl (q) ∂q i q+k q+l dt ≤ 2 √ λ max V m + F + d 2 GV 2 m 2λ min (s 2 -s 1 ). ( 38
)
There results that the λ i are d bounded measures on ]t 0 , T [. Thanks to equation (36), it is readily seen that the measures d q+i are also bounded measures on ]t 0 , T [. Therefore, the d functions q+i :]t 0 , T [→ R have bounded variation over the interval ]t 0 , T [. By proposition 29, there results that q+ has also bounded variation over ]t 0 , T [.

Proof of theorem 10. We assume that the maximal solution q of problem P is defined on [t 0 , T [ with T in • I and try to obtain contradiction. By proposition 7, this maximal solution satisfies:

∀t ∈ [t 0 , T [, 1 2 q+ (t) 2 q(t) - 1 2 v 0 2 q 0 ≤ t t 0
f (q(s), q+ (s); s), q+ (s) q(s) ds.

Thus, ∀t ∈ [t 0 , T [, 1 2 q+ (t) 2 q(t) ≤ 1 2 v 0 2 q 0 + t t 0
f (q(s), q+ (s); s) q(s) q+ (s) q(s) ds.

By lemma 17, we obtain:

∀t ∈ [t 0 , T [, q+ (t) q(t) ≤ v 0 q 0 + t t 0
f (q(s), q+ (s); s) q(s) ds, which gives, using the hypothesis of the theorem:

∀t ∈ [t 0 , T [, q+ (t) q(t) ≤ v 0 q 0 + t t 0 l(s) 1 + d(q(s), q 0 ) + q+ (s) q(s) ds. But, ∀t ∈ [t 0 , T [, d(q(t), q 0 ) ≤ t t 0 q+ (s) q(s) ds, therefore, ∀t ∈ [t 0 , T [, d(q(t), q 0 ) + q+ (t) q(t) ≤ v 0 q 0 + t 0 l(s) ds + t t 0
(1 + l(s)) d(q(s), q 0 ) + q+ (s) q(s) ds.

By Gronwall-Bellman lemma (lemma 15), one gets:

∀t ∈ [t 0 , T [, d(q(t), q 0 ) + q+ (t) q(t) ≤ v 0 q 0 + t 0 l(s) ds e t t 0 (1+l(s)) ds ,
which shows that the function t → q+ (t) q(t) is bounded over [t 0 , T [. By the completeness of Q, we deduce, on one hand that

q T = lim t→T - q(t)
exists in Q and, on the other hand, that

Var q+ ; [t 0 , T [ < ∞,
thanks to proposition 18. Thus,

(q T , v - T ) = lim t→T -(q(t), q+ (t)) exists in T Q. Define: v T = v - T -1 + φ(q T , v - T ) Proj q T v - T ; N (q T
) . Then, theorem 8 furnishes T ∈ I with T > T and a prolongation of q on [T, T [ such that q ∈ M M A ([t 0 , T [; Q) is a solution of problem P. But, this contradicts the definition of T .

Three counterexamples

The existence and uniqueness of solution for problem P has been proved under the assumption of functional independence for the constraint and of analyticity for the data. The three examples which are developed in this section aim at showing that these assumptions cannot be weakened very much. In example 1, we show that, in case that the functional independence of the constraints does not hold, the existence of solution may be lost. For the question of the regularity assumptions on the data, the existence of solution can be proved with much weaker assumptions. However, the uniqueness of solutions is generally lost in such a case as seen in examples 2 and 3. In these examples, the data are supposed to have only regularity C ∞ and two different solutions can be exhibited.

Example 1 is extracted from Moreau [START_REF] Moreau | Standard inelastic shocks and the dynamics of unilateral constraints, in Unilateral problems in structural analysis[END_REF] and example 2 is due to Schatzman [START_REF] Schatzman | A Class of Nonlinear Differential Equations of Second Order in Time, Nonlinear Analysis[END_REF], but an earlier counterexample in the same spirit is also to be found in Bressan [START_REF] Bressan | Incompatibilità dei Teoremi di Esistenza e di Unicità del Moto per un Tipo molto Comune e Regolare di Sistemi Meccanici[END_REF].

Example 1

Consider a discrete mechanical system whose configuration space is euclidean R 3 . The unilateral contraints are kinematically described by the three following functions (n = 3): ϕ 1 (q) = -q 1 , ϕ 2 (q) = q 1q 2 .q 3 , ϕ 3 (q) = -q 2q 3 , where q = (q 1 , q 2 , q 3 ) ∈ R 3 . The initial instant is supposed to be t 0 = 0 and the initial state is given by q 0 = (0, 0, 0) and v 0 = (0, 2, -1). It follows that:

J(q 0 ) = {1, 2, 3} , V (q 0 ) = v = (v 1 , v 2 , v 3 ) ∈ R 3 ; v 1 = 0 and v 2 + v 3 ≥ 0 .
It is readily seen that v 0 belongs to V (q 0 ).

Let now α > 0 be an arbitrary positive real number. Any motion q(t) in M M A([0, α[; R 3 ) compatible with this initial data may be written as:

q 1 (t) = o(t) q 2 (t) = 2t + o(t) q 3 (t) = -t + o(t) Therefore: ϕ 1 (q(t)) + ϕ 2 (q(t)) = 2t 2 + o(t 2 ),
which can not be compatible with:

∀t ∈ [0, α[ ϕ 1 (q(t)) + ϕ 2 (q(t)) ≤ 0.
We deduce that no motion in M M A([0, α[; R 3 ) can be compatible with this initial data whatever α > 0 is. Note that in this particular case: dϕ 1 (q 0 ) = -dϕ 2 (q 0 ) and the unilateral constraints are not functionally independent.

Example 2

Consider a discrete mechanical system whose configuration space is R equipped with its canonical structure of riemannian manifold. This is the configuration space of a particle with unit mass constrained to move along a line. A fixed obstacle at the origin is taken into consideration. It gives rise to a unilateral constraint kinematically described by the single function (n = 1):

ϕ 1 (q) = q
Therefore, the admissible configuration set is A = R -. It is assumed that the impact constitutive equation is the elastic one: φ (q, q-) ≡ 1 and that the efforts mapping f does not depend on the state but only on time. It will be denoted by f (t). The initial instant is t 0 = 0 and the initial state is (q 0 , v 0 ) = (0, 0). Denoting by RCLBV (I; R) the space of right continuous functions with locally bounded variation from a real interval I to R, problem P admits here the equivalent formulation: find T > 0 and v ∈ RCLBV ([0, T [; R) such that:

• v(0) = 0, • q(t) = t 0 v(s) ds ∈ R -, ∀t ∈ [0, T [, • R = dv -f (t) dt is a nonpositive real measure such that: Supp R ⊂ {t ∈ [0, T [ ; q(t) = 0} • ∀t ∈]0, T [, q(t) = 0 ⇒ v(t) = v -(t) q(t) = 0 ⇒ v(t) = -v -(t)
We investigate uniqueness under the assumption that f is of class C ∞ . Suppose, in addition, that f is nonnegative:

∀t ∈ R + , f(t) ≥ 0.
It is readily seen that the null function v ≡ 0 on R + is a solution of problem P whatever is the nonnegative C ∞ function f . Now, we are going to construct an explicit example of such a function f in such a way that the associated problem P admits another solution, different from the identically vanishing one. First, let define a function ρ by:

ρ ⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ R → R x → 0 if x ∈] -∞, 0] ∪ [1, +∞[ e 1 x(x-1) 1 0 e 1 x(x-1) dx if x ∈]0, 1[ We have: ρ ∈ C ∞ (R; R + ) Suppρ = [0, 1] ∀n ∈ N, d n dx n ρ(0) = d n dx n ρ(1) = 0 2 1 0 (1 -s)ρ(s) ds = 1 (39)
The last assertion comes from the fact that:

1 0 sρ(s) ds = 1 0 (1 -s)ρ(s) ds, so, 1 0 sρ(s) ds = 1 2 1 0 ρ(s) ds = 1 2
Consider also the real convergent serie:

(n + 5) 2 (n + 1)(n + 2)(n + 3)(n + 4) n∈N
We define:

T = ∞ n=0 (n + 5) 2 (n + 1)(n + 2)(n + 3)(n + 4) > 0 a n = ∞ i=n (i + 5) 2 (i + 1)(i + 2)(i + 3)(i + 4)
Clearly, a 0 = T and the sequence (a n ) n∈N decreases strictly and converges towards 0 when n tends toward infinity. Actually,

a n ∼ 1 n when n → +∞ (40)
by a very classical and elementary argument. We denote by (δ n ) n∈N , (f n ) n∈N , (v n ) n∈N the real sequences defined by:

δ n = n + 5 (n + 1)(n + 2)(n + 4) (i.e. δ n = n + 3 n + 5 (a n -a n+1 ) < a n -a n+1 ), f n = 1 n! , v n = - 1 (n + 3)! ,
and by f (t), v(t) the functions from [0, T [ to R defined by:

f (0) = 0 f (t) = 0 if t ∈ [a n+1 , a n+1 + δ n [ f n 2 ρ t-a n+1 -δ n a n -a n+1 -δ n if t ∈ [a n+1 + δ n , a n [ (41)
and:

v(0) = 0 v(t) = v n+1 if t ∈ [a n+1 , a n+1 + δ n [ v n+1 + f n 2 t a n+1 +δ n ρ s-a n+1 -δ n a n -a n+1 -δ n ds if t ∈ [a n+1 + δ n , a n [ (42) First, we claim that the function f belongs to C ∞ ([0, T [; R). Proof. The only thing which is not obvious is that f is C ∞ at 0. Since: ∀t ∈ [a n+1 , a n ] |f (t)| ≤ f n 2 max s∈[0,1] |ρ(s)| ,
then, lim t→0 + f (t) = 0 and f is continuous at 0. Now, we are going to prove:

∀r ∈ N lim t→0 + 1 t d r dt r f (t) = 0 (43)
which will imply by an easy induction that f ∈ C ∞ ([0, T [; R) and:

∀r ∈ N, d r dt r f (0) = 0.
Let fix an arbitrary r in N, we have:

∀t ∈ [a n+1 , a n ] , 1 t d r dt r f (t) ≤ f n 2a n+1 (n + 5) r 2 r (a n -a n+1 ) r max s∈[0,1] d r ρ(s) dt r (t) .
Therefore, to prove (43), it suffices to verify:

lim n→∞ f n (n + 5) r a n+1 (a n -a n+1 ) r = 0,
but, by estimate (40), one has:

f n (n + 5) r a n+1 (a n -a n+1 ) r ∼ n 3r+1 n! .
Second, we claim that:

-v ∈ RCLBV([0, T [;R) -dv -f (t) dt is a real nonpositive measure on [0, T [ whose support is {0} ∪ {a n ; n ∈ N * } -v is continuous on [0, T [\ {a n ; n ∈ N * }and: ∀n ∈ N * v(a n ) = -v -(a n ) Proof. It is clear that v is continuous on each interval ]a n+1 , a n [ and right contin- uous on [0, T [. Moreover: v -(a n ) = v n+1 + f n 2 a n a n+1 +δ n ρ s -a n+1 -δ n a n -a n+1 -δ n ds = v n+1 + f n 2 (a n -a n+1 -δ n ) = - 1 (n + 4)! + 1 n! n + 5 (n + 1)(n + 2)(n + 3)(n + 4) = 1 (n + 3)! = -v(a n ) Since v is nondecresing on each interval [a n+1 , a n [: Var(v; [0, T [) = ∞ n=0 v(a n+1 ) -v -(a n+1 ) + v(a n+1 ) -v -(a n ) = ∞ n=0 (-3v n+1 -v n ) = 3 ∞ n=0 1 (n + 4)! + ∞ n=0 1 (n + 3)! < +∞
Denoting δ t the dirac measure located at t, one has:

dv -f (t) dt = -2 ∞ n=1 δ a n (n + 3)! , which is a (bounded) nonpositive measure whose support is {0} ∪ {a n ; n ∈ N * }.
Third, we claim that: if q is defined by:

∀t ∈ [0, T [ q(t) = t 0 v(s) ds, then: ∀t ∈ [0, T [ q(t) ≤ 0, {t ∈ [0, T [; q(t) = 0} = {0} ∪ {a n ; n ∈ N * } .
Proof. An easy calculation using last assertion of formulae (39) shows:

a n a n+1 v(s) ds = 0 ∀t ∈]a n+1 , a n [ t a n+1 v(s) ds < 0
We deduce, that making the choice described by relations (41) for the function f , then the function v defined by relations ( 42) is a solution of the corresponding problem P whereas the identically vanishing function is also a solution. Therefore, the uniqueness of solution does not hold in general if f and the functions ϕ i are supposed to be of class C ∞ only.

Example 3

In example 2, we considered a particle at rest at the initial instant and in contact with the obstacle. Then, a force acts on the particle, constantly pushing it against the obstacle (f ≥ 0). For the particular choice of the function f we made, immobility is a possible motion whereas a bouncing motion is also possible. It is intuitively clear that the assumed elastic impact constitutive equation plays a central role in such a phenomenon. The question arises to know whether such a pathology is possible with the completely inelastic impact constitutive equation φ(q, q-) ≡ 0.

Sticking to notations of example 2, the evolution problem takes in this case the form: find T > 0 and v ∈ RCLBV([0, T [; R) such that:

• v(0) = 0 • q(t) = t 0 v(s) ds ∈ R - ∀t ∈ [0, T [ • R = dv -f (t) dt is a nonpositive real measure such that: Supp R ⊂ {t ∈ [0, T [ ; q(t) = 0} • ∀t ∈]0, T [, q(t) = 0 ⇒ v(t) = v -(t) q(t) = 0 ⇒ v(t) = 0
where δ 1,n has to be the smallest root of this second degree equation:

δ 1,n = v n+1 -v 2 n+1 -2f 1,n q n+1 f 1,n (45)
We have also to express that, under the action of f on [ 1 n+1 , 1 n ], a particle at rest with position q = 0 at time t = 1 n+1 should have position q = -q n and velocity v = v n at time t = 1 n . That is:

v n = -f 1,n δ 1,n + f 2,n δ 2,n , -q n = -1 2 f 1,n δ 2 1,n -f 1,n δ 1,n 1 n(n+1) -δ 1,n + 1 2 f 2,n δ 2 2,n ,
which is:

v n = -f 1,n δ 1,n + f 2,n δ 2,n , -q n = 1 2 f 1,n δ 2 1,n -f 1,n δ 1,n 1
n(n+1) + 1 2 f 1,n δ 1,n δ 2,n + 1 2 v n δ 2,n . ( 46
)
Now, let us try to make the following choice:

∀n ∈ N * , q n = 1 n 4 2 n , v n = 1 2 n , f 1,n = n 3 2 n . ( 47
)
There results from formula (45) that, for sufficiently great n: 4 (48)

δ 1,n = 1 2n 3 1 -1 - 4n 3 (n + 1)
which gives the estimate:

δ 1,n ∼ 1 n 4 when n → ∞ (49)
Equations ( 46) allow us to determine δ 2,n and f 2,n :

δ 2,n = 2n 2 n+1 δ 1,n -n 3 δ 2 1,n -2 n 4 1 + n 3 δ 1,n , f 2,n = f 1,n δ 1,n δ 2,n + v n δ 2,n ,
which provide the estimates:

δ 2,n ∼ 2 n 3 f 2,n ∼ n 3 2 n+1 when n → ∞ (50)
From estimates (49) and ( 50), one gets:

∃n 0 , n ≥ 0 ⇒ 0 < δ 1,n < 1 2n(n + 1) 0 < δ 2,n < 1 2n(n + 1)
We define T = 1 n 0 . In exactly the same way as for example 2, it is readily seen from estimate (50) that: f ∈ C ∞ ([0, T [; R). Define: u I (0) = 0, u II (0) = 0, and for n ≥ n 0 :

u I (t) = v n+1 -f 1,n t 1 n+1 ρ s-1 n+1 δ 1,n ds if t ∈ [ 1 n+1 , 1 n+1 + δ 1,n [ 0 if t ∈ [ 1 n+1 + δ 1,n , 1 n [ u II (t) = -f 1,n t 1 n+1 ρ s-1 n+1 δ 1,n ds t ∈ [ 1 n+1 , 1 n+1 + δ 1,n [ -f 1,n δ 1,n t ∈ [ 1 n+1 + δ 1,n , 1 n -δ 2,n [ -f 1,n δ 1,n + f 2,n t 1 n -δ 2,n ρ s-1 n +δ 2,n δ 2,n ds t ∈ [ 1 n -δ 2,n , 1 n [ and: v I (0) = 0, v II (0) = 0, v I (t) = u I (t) v II (t) = u II (t) if t ∈ [ 1 2p+1 , 1 2p [ (2p ≥ n 0 ) v I (t) = u II (t) v II (t) = u I (t) if t ∈ [ 1 2p , 1 2p-1 [ (2p -1 ≥ n 0 )
Proceeding as in example 2, it is readily seen that the two functions v I and v II belong to RCLBV ([0, T [;R) and furnish two distinct solutions of the problem P associated with the C ∞ function f defined by equations (44).

Illustrative examples and comments

Punctual particle subjected to gravity and bouncing on the floor. Accumulation of impacts.

Let us come back to the example of section 3.3. The configuration space is R equipped with its canonical structure of riemannian manifold, the unilateral constraint is described by the single function ϕ 1 (q) = q (which gives A = R -). The efforts mapping is supposed to be constant: f (q, q; t) ≡ 2 and the impact function (restitution coefficient) is the constant 1/2: φ ≡ 1/2. Considering the initial instant t 0 = 0 and the initial state (q 0 , v 0 ) = (-1, 0), we have seen in section 3.3 that the function q : R + → R -defined by:

∀t ∈ [0, 1], q (t) = t 2 -1, ∀t ∈ [1, 2], q (t) = t 2 -3t + 2, ∀t ∈ 3 -1 2 n-1 , 3 -1 2 n , q(t) = t 2 + -6 + 3 2 n t + 3 -1 2 n-1 3 -1 2 n , ∀t ∈ [3, +∞[, q (t) = 0 (∀n ∈ N) belongs to M M A(R + ; R -)
and is readily seen to be the maximal solution, according to corollary 9, of the corresponding problem P. The solution q(t) is represented on figure 6.1. It is seen that infinitely many impacts accumulate in any left neighbourhood of instant t = 3. However, as predicted by corollary 9, for each instant t ∈ R + , there exists a right neighbourhood [t, t + η[ of t, such that the restriction of q to [t, t + η[ is analytic. A straightforward and general consequence of this is the following.

Proposition 19 Let q be the maximal solution of problem P furnished by corollary 9. Although infinitely many impacts can accumulate at the left of a given instant, this phenomenon can never occur at the right of any instant. Morover, in the particular case where the impact constitutive equation is the elastic one (φ ≡ 1), the instants of impact are isolated and therefore in finite number in any compact interval of time.

Proof. Since for each instant t ∈ [t 0 , T [, there exists a right neighbourhood [t, t+η[ of t, such that the restriction of q to [t, t + η[ is analytic, we get the first part of the proposition. For the second part, let τ be an arbitrary instant in ]t 0 , T [ and consider the problem P associated with the initial condition (q(τ ), -q-(τ )), the elastic constitutive impact equation and the effort mapping g(q, v; t) defined by:

g(q, v; t) = f (q, -v; τ -t)
which is analytic. By theorem 8, there exists an analytic function q a : [0, T a [→ Q which is a solution of this problem P. Another solution of problem P coincides with q a on a right neighbourhood of t = 0. Actually, as seen in the proof of local uniqueness (section 4.2), a little bit more is proved: any function q ∈ M M A([0, T [; Q) satisfying the initial condition (21), the unilateral constraint ( 22), the equation of motion (23) and the energy inequality (proposition 7) has to coincide with q a on a right neighbourhood of t = 0. But, it is readily seen that the function q (t) = q(τt), t ∈ [0, τt 0 [ fulfill these requirements. Thus, q can not have right accumulation of impacts at t = τ and, therefore, q can not have left accumulation of impacts at t = τ and the instants of impact are isolated. Of course, if q is the maximal solution defined on [t 0 , T [, impacts can still accumulate at the left of T , as seen on simple examples. QED.

The fact that infinitely many impacts can accumulate at the left of a given instant but not at the right is a specific feature of the analytical setting that is lost in the C ∞ setting as seen in counterexamples 2 and 3. Actually, these counterexamples show that pathologies of nonuniqueness in the C ∞ setting are intimately connected to the possibility of right accumulations of impacts. The fact that the analytical setting prevents from such right accumulations is the thorough reason why we could prove uniqueness in this case.

The double pendulum

In this section, we come back to the double pendulum described in section 2.1 but we add to the system a rigid obstacle on the vertical coordinate axis as represented on figure 6.2. This obstacle may be represented by two analytic functions whose expressions in the global chart of Q described in section 2.1 is:

ϕ 1 (q 1 , q 2 ) = -l 1 sin q 1 ≤ 0, ϕ 2 (q 1 , q 2 ) = -l 1 sin q 1 -l 2 sin q 2 ≤ 0.
It is readily seen that, except in the particular case where l 1 = l 2 , these constraints are functionally independent: ∀q ∈ A, (dϕ i (q)) i∈J(q) is linear independent in T * q Q These unilateral constraints are assumed to be perfect and we consider an impact function φ supposed to be constant on T A -:

∀(q, v -) ∈ T A -, φ(q, v -) ≡ φ ∈ [0, 1].
The constant φ is often called restitution coefficient (of normal velocities). We recall that the particular cases φ = 0 and φ = 1 describe the completely inelastic and the elastic impact constitutive equations. An initial state (q 0 , v 0 ) ∈ T A + is given at time t 0 = 0. This initial state is represented in the considered chart by four real numbers (q 1 0 , q 2 0 ; v 1 0 , v 2 0 ). According to section 3, the motion of the system is governed by the evolution problem:

l , m 1 1 1 q 1 k k 2 l , m 2 2 2
q λ Fig. 3. Double pendulum with obstacle.

find T ∈]0, +∞] and q ∈ M M A([0, T [; Q) such that:

• (q(0), q+ (0)) = (q 0 , v 0 ),

• ∀t ∈ [0, T [, (q(t), q+ (t)) ∈ T A + , • R = D q+ -f (q(t), q+ (t); t) dt, ∈ -N * (q(t)) for |R| -a.a. t ∈ [0, T [, • ∀t ∈]0, T [, q+ (t) = q-(t) -(1 + φ) Proj q(t) [ q-(t); N (q(t))] ,
where the riemannian structure on Q and the mapping f are those described in section 2.1. Corollary 9 ensures existence and uniqueness of a maximal solution. Now, we are going to check that assuptions of theorem 10 are satisfied so that the maximal solution is defined all over R + . First, Q is a complete riemannian manifold since the quotient topology on the torus T 2 derives from a riemannian structure and T 2 is compact and therefore complete. Second, we have the estimate:

∀(q, v) ∈ T Q, v q ≥ α |(v 1 , v 2 )| , (51)
where,

α = 1 9 m 1 m 2 l 2 1 l 2 2 + 1 12 m 2 2 l 2 1 l 2 2 1 3 m 2 l 2 2 + m 1 3 + m 2 l 2 1 Indeed, v 2 q ≥ λ min (q) (v 1 , v 2 )
where λ min (q) is the least eigenvalue of the matrix (g ij (q)) i,j=1,2 . But:

λ min (q) = 1 2 1 3 m 2 l 2 2 + m 1 3 + m 2 l 2 1 × ⎡ ⎣ 1 -1 - 4 1 9 m 1 m 2 l 2 1 l 2 2 + 1 3 m 2 2 l 2 1 l 2 2 -1 4 m 2 2 l 2 1 l 2 2 cos(q 1 -q 2 ) 1 3 m 2 l 2 2 + m 1 3 + m 2 l 2 1 2 ⎤ ⎦ Using ∀x ∈ [0, 1], 1 - √ 1 -≥ x 2 ,
we get:

λ min (q) ≥ 1 4 4 1 9 m 1 m 2 l 2 1 l 2 2 + 1 3 m 2 2 l 2 1 l 2 2 -1 4 m 2 2 l 2 1 l 2 2 cos(q 1 -q 2 ) 1 3 m 2 l 2 2 + m 1 3 + m 2 l 2 1 ≥ α 2 ,
which achieves the proof of estimate (51). Now, let q I , q II be two points of Q represented by their components in the considered chart (q 1 I , q 2 I ) and (q 1 II , q 2 II ). Q being complete, there is a geodesic g : [s 1 , s 2 ] → Q of minimal length between them. One has:

d(q I , q II ) = s 2 s 1 ġ(s) g(s) ds ≥ s 2 s 1 α | ġ(s)| ds ≥ α (q 1 I -q 1 II ) 2 + (q 2 I -q 2 II ) 2
Moreover, recalling:

f 1 (q 1 , q2) = λl 1 sin(q 1 -q 2 ) -(k 1 + k 2 )q 1 + k 2 q 2 , f 2 (q 1 , q 2 ) = k 2 q 1 -k 2 q 2 ,
one has:

f (q) 2 q ≤ 1 λ min (q) |(f 1 , f 2 )| 2 .
Therefore,

f (q) q ≤ 1 α |(f 1 , f 2 )| , ≤ 1 α λl 1 + (k 1 + k 2 ) q 1 + k 2 q 1 + 2k 2 q 2 , ≤ 1 α λl 1 + 4(k 1 + k 2 ) (q 1 0 , q 2 0 ) + 4(k 1 + k 2 ) (q 1 -q 1 0 , q 2 -q 2 0 ) , ≤ 1 α λl 1 + 4(k 1 + k 2 ) (q 1 0 , q 2 0 ) + 4(k 1 + k 2 ) α 2 d(q, q 0 ), ∀q ∈ Q
By virtue of theorem 10, the motion of the system is defined for all t ∈ R + .

Boltzmann's gas

Consider a collection of N rigid homogeneous balls of mass m and radius R in a rigid rectangular box. The balls cannot interpenetrate. The balls are free of internal or external forces except for the reaction efforts induced by the unilateral constraints. The impact constitutive equation is supposed to be the elastic one. Such a system was introduced by Boltzmann to model the interactions between molecules in a gas in order to perform a statistical analysis to connect the microscopical and macroscopical point of view.

Let us describe the discrete mechanical system associated with this situation. The configuration space is R 3N . Indeed, any configuration is described by the coordinates of the center of the balls in the three-dimensional ambiant space equipped with an origin. Strictly speaking, the configuration space should be R 3N × (SO3) N to incorporate the possible rotations of the balls. But, in this case, it would be readily seen that the rotation velocity of any ball in any motion of the system keeps its value at the initial instant. Therefore, rotations play no role in the motion of the system and we may consider only the restricted configuration space R 3N equipped with its canonical riemannian structure. The forces mapping vanishes identically f (q, q+ ; t) ≡ 0. There are N (N + 11)/2 functions ϕ i , since N (N -1)/2 of them are necessary to express the noninterpenetration constraints:

∀i, j ∈ {1, 2, • • • N }, i = j, x i -x j 2 + y i -y j 2 + z i -z j 2 ≥ R 2 ,
and 6N of them are necessary to express that the balls remains inside the box:

∀i, j ∈ {1, 2, • • • N }, -a + R ≤ x i ≤ a -R, -b + R ≤ y i ≤ b -R, -c + R ≤ z i ≤ c -R,
where 2a, 2b and 2c are the lengths of the sides of the box. The functions ϕ i are defined by arbitrary numbering. They are easily seen to be analytic and functionally independent. Adding the elastic impact constitutive equation φ(q, q-) ≡ 1, and an initial condition at time t 0 = 0, the corresponding evolution problem turns out to belong to the class of problem P formulated at the beginning of section 4. Then, corollary 9 and theorem 10 state that, to any initial condition compatible with the constraints, there corresponds a unique maximal motion and it is defined all over R + . By proposition 19, we may also state that there are at most finitely many impacts on any bounded time interval. As a conclusion, the results developed in this paper allow us to associate a dynamical system with Boltzmann's gas. Related to this question, let us mention Boltzmann's famous ergodic hypothesis. Roughly speaking, Boltzmann postulated that in any motion of the system, time averages can be replaced by space averages. The modern mathematical transcript is: for almost every initial condition in an energy level set of the phase space, the associated phase curve spends an amount of time in every measurable piece of the level set proportional to the measure of that piece. The question to know whether Boltzmann's gas is ergodic or not is still an open question. However, a positive answer was given in 1970 by Sinai ([20]) for a two balls gas in a plane rectangular box. Let us underline that this question makes sense only when we are able to associate a dynamical system with Boltzmann's gas.

Newton's balls and the impact constitutive equation

In section 3.3, it has been derived from the two phenomenological assumptions H3 and H4, that the general constitutive impact equation

q+ = F q, q- (52)
should satisfy:

∀q ∈ A, ∀v -∈ -V (q), F (q, v -) ∈ V (q) F (q, v -) -v -∈ -N (q) F (q, v -) q ≤ v - q (53)
In the particular case of a motion no more than one active constraint at any time (∀t, CardJ(q(t)) ≤ 1), it has been seen in section 3.3 that the general impact constitutive equation ( 52) necessarily takes the form:

q+ = Proj q q-; V (q) -φ q, q-Proj q q-; N (q) , (54)
with φ an arbitrary function taking values in the interval [0, 1]. Actually, equation (54) makes sense even in case of multiple impacts and it is a simple example of an impact constitutive equation satisfying requirements (53). For sake of simplicity, we have adopted this particular form of the impact constitutive equation even in the case where multiple impacts occur. However, the reader should keep in mind the arbitraryness of this choice and we shall show that it could be irrelevant in some cases. A simple occurrence of multiple impact is Newton's balls experiment. The principle of Newton's balls experiment is well-known. It is sketched on figure 6.4a. As a result of this multiple impact experiment, we have the familiar picture drawn of figure 6.4b. But, testing the simple impact constitutive equation (54) (with φ ≡ 1) to predict the issue of the experiment, we get the situation drawn on figure 6.4c.

The question arises to known whether the results of section 4 remain true if we abandon the simple impact constitutive equation (54) and adopt the general impact constitutive equation (52) defined by an arbitrary function F fulfilling requirements (53). Actually, a careful examination of the proofs of section 4 shows that the impact constitutive equation is only used through the energy inequality (proposition 7). Moreover, it is readily seen that proposition 7 still holds when the simple impact constitutive equation ( 24) is replaced by a general one (equation (52)) provided requirements (53) hold true. As a result, all the results of section 4, and in particular, theorem 8, corollary 9 and theorem 10 remain true if we adopt an arbitrary impact constitutive equation instead of equation 24 in the definition of problem P.

A general impact constitutive equation will be said elastic if the last requirement in (53) is replaced by:

∀q ∈ A, ∀v -∈ -V (q), F q, v - q = v - q .
It is readily seen that proposition 19 still holds with an arbitrary impact constitutive equation. In particular, for a solution of problem P with an arbitrary elastic impact constitutive equation, the impacts are isolated.

Continuous dependence on initial conditions

The theory developed in the previous sections allows us to replace the analysis of the motion of a collection of rigid bodies subjected to perfect constraints either bilateral or unilateral by the analysis of the motion of a point in a piece of a d-dimensional manifold bounded by analytic hypersurfaces which intersect transversally. With appropriate regularity assumptions on the data, the motion is completely determined by the initial condition.

The picture seems to be fairly good and the generalization of the dynamics of discrete systems with bilateral constraints to the case of unilateral constraints seems to be achieved. However, there remains a big difference between unilateral and bilateral dynamics of discrete systems that we want to underline in this section.

A pleasant feature of a dynamical system generated by the flow of an ordinary differential equation is that it is smooth. More precisely, if F t,t 0 is the mapping which associates the state of the system at time t with an arbitrary initial condition at time t 0 , then the mapping F t,t 0 is a local diffeomorphism. In particular, the state of the system at a given instant t depends in a differentiable way of the state at time t 0 . Of course, this smooth dependence may be stiff. In such a case, a small uncertainty on the initial state will produce a big one on the actual state and the motion of the system may turn out to be quantitatively unpredictable on both physical and numerical point of view for large time. In certain circumstances, the theory of smooth dynamical systems allow us to get some qualitative information on the motion for large time.

As we shall see, the picture is strongly different in the case of the dynamics of discrete systems with perfect unilateral constraint. The theorems of section 4 allow us to define a mapping F t,t 0 similar to the flow generated by an ordinary differential equation. But, the mapping F t,t 0 is not smooth any more, it is not even continuous in general. In other terms, the generated dynamical system does not belong to the large class of topological dynamical systems.

Let us check this assertion on a simple example. Consider as a configuration space R 2 supplied with its canonical structure of riemannian manifold. A configuration is denoted by a pair (x, y). No forces act on the system: f ≡ 0. Consider a unilateral constraint associated with the two functions:

ϕ 1 (x, y) = y ≤ 0 ϕ 2 (x, y) = x + y ≤ 0
and the elastic impact constitutive equation φ ≡ 1. At time t 0 = 0, we consider the following set of initial conditions:

{(-1 + ε, -1; 1, 1); ε ∈] -1, 1[}
A straightworward calculation gives the state of the system for all instant in R + .

In particular, for t greater than 1, one gets:

F t,0 (-1 + ε, -1; 1, 1) = (-1 + ε + t, 1 -t, 1, -1) if ε ∈] -1, 0[, F t,0 (-1 + ε, -1; 1, 1) = (1 -t, 1 -ε -t, -1, -1) if ε ∈ [0, 1[.
It is readily seen on this example that, if t is greater than 1, the mapping F t,0 is not continuous at initial condition (-1, -1, 1, 1) (see figure 7). Coming back to x y Fig. 5. The generated dynamical system is not continuous in general.

the two examples of section 6, such a situation occurs if, during the motion of the double pendulum, the two bars hit the obstacle at the same time. In the case of Boltzmann's gas, the pathology occurs when three balls hit at the same time. Let us underline that if we consider an initial condition such as the one in the above example, the solution of the associated problem P has no physical meaning. In such a case, one has to abandon the will of predicting the motion of the system: this is a consequence of the over-idealization made in the indeformability assumption.

However, in the particular case of the one degree-of-freedom problem, where no multiple impacts are possible, Schatzman [START_REF] Schatzman | Uniqueness and continuous dependence on data for one dimensional impact problems[END_REF] proved that continuous dependence on initial conditions hold. In the general case, her result admits the following generalization which is proved along the same lines: Theorem 20 Consider the problem P described in section 3.4. Assume furthermore that the impact function φ is constant. Considering the initial condition (q 0 , v 0 ) ∈ T A + at initial instant t 0 , we denote by (T, q) the corresponding maximal solution of problem P. We make the following hypothesis:

∀t ∈ [t 0 , T [, (dϕ i (q(t))) i∈J(q(t)) is orthogonal in T * q(t) Q,
(with the convention that the empty set is orthogonal). Let us onsider a sequence (q 0n , v 0n ) of elements of T A + converging towards (q 0 , v 0 ). For all n, we denote by (T n , q n ) the maximal solution of the problem P associated with the initial condition (q 0n , v 0n ) at instant t 0 . Then,

1. lim inf n→+∞ T n ≥ T ,
2. q n converges towards q uniformly on every compact subset of [t 0 , T [:

∀τ ∈ [t 0 , T [, lim n→+∞ sup t∈[t 0 ,τ ] d(q n (t), q(t)) = 0, 3. (q n (t), q+ n (t)) converges towards (q(t), q+ (t)) in T Q for almost all t in [t 0 , T [.
Proof. The proof of theorem 20 is divided into five steps. Before describing these steps, let us introduce a few notations. We fix, once for all, an arbitrary τ in [t 0 , T [ and a compact neighbourhood K of the compact subset q([t 0 , τ ]) of Q. We define:

V = 1 + sup t∈[t 0 ,τ ] q+ (t) q(t) , and, K = (q, v) ∈ T Q; q ∈ K and v q ≤ 4V .
The subset K of T Q is compact in T Q. We define also:

F = 1 + max (q,v;t)∈K×[t 0 ,τ ] f (q, v; t) q and, d 0 = min (q ,t)∈∂K×[t 0 ,τ ] d(q , q(t)),
and,

δ = min V F , d 0 6V .
Notice that we have δ > 0.

Step 1. Consider t 1 ∈ [t 0 , τ [. We denote q(t 1 ) by q 1 and q+ (t 1 ) by v 1 . Consider an element (q 1 , v 1 ) of T A + such that:

d(q 1 , q 1 ) ≤ d 0 4 and v 1 q 1 ≤ 2V.
Then, the maximal solution q of the problem P associated with the initial condition (q 1 , v 1 ) at initial instant t 1 is defined on an interval containing [t 1 , min(τ, t 1 + δ)] and is such that:

∀t ∈ [t 1 , min(τ, t 1 + δ)], q (t), q + (t) ∈ K.

Proof of step 1. Let us denote by [t 1 , T 1 [ the maximal definition interval of q . Define:

t 1 = sup t ∈ [t 1 , T 1 [; ∀s ∈ [t 1 , t], q (s), q + (s) ∈ K .
We have to prove:

t 1 ≥ min (τ, t 1 + δ) .
Assume the contrary is true:

t 1 < min (τ, t 1 + δ) .
By proposition 7 and lemma 17, we have:

∀t ∈ [t 1 , t 1 [, q + (t) q (t) ≤ v 1 q 1 + t t 1 F ds ≤ 2V + F (t 1 -t 1 ) ≤ 3V
We deduce:

t 1 < T 1 ,
by proposition 18, and:

q + (t 1 ) q (t 1 ) ≤ q -(t 1 ) q (t 1 ) = lim t→t - 1 q + (t) q (t) ≤ 3V,
by proposition 32. Moreover, d (q (t 1 ), q 1 ) ≤ d (q (t 1 ), q 1 ) + d (q 1 , q 1 )

≤ 3V (t 1 -t 1 ) + d 0 4 ≤ 3 4 d 0
By the continuity of the function t → d(q (t), q 1 ) and the right-continuity of the function t → q + (t) q (t) , there results:

∃α > 0, ∀t ∈ [t 1 , t 1 + α], q (t), q + (t) ∈ K.
But, this contradicts the definition of t 1 and achieves the proof of step 1.

Step 2. For n large enough, q n is defined on (an interval containing) the interval [t 0 , min(τ, t 0 + δ)]. Moreover, there exists a subsequence of (q n ), also denoted by (q n ), such that:

q n converges uniformly on [t 0 , min(τ, t 0 +δ)] towards a function q lim belonging to M M A([t 0 , min(τ, t 0 + δ)]; Q), -(q n (t), q+ n (t)) converges towards (q lim (t), q+

lim (t)) in T Q for almost all t in [t 0 , min(τ, t 0 + δ)].
Proof of step 2. For all q in K ∩ A, there exists a compact neighbourhood K q of q which is included in the domain U q of a local chart (U q , ψ q ) at q such that:

-∀q ∈ U q , J(q ) ⊂ J(q), -∀q ∈ U q , the cardJ(q) first components of ψ q (q ) are the ϕ i (q ) (i ∈ J(q)).

Being compact, K ∩ A can be covered by a finite number, say L, of • K q l . We denote by λ max (resp. λ min ) the maximum (resp. the minimum) of the greatest (resp. least) eigenvalue of the matrix (g ij (q)) i,j=1,2,•••,d when q wanders in K q l and l in {1, 2, • • • , L}. We define:

G = max i,j,k∈{1,2,•••,d} l∈{1,2,•••,L} max q∈K q l ∂g ij (q) ∂q k .
We pick an integer N 0 , large enough to ensure:

∀n ≥ N 0 , d(q 0 , q 0n ) ≤ d 0 4 v 0n q 0n ≤ 2V .
By step 1, there results:

∀n ≥ N 0 , T n ≥ min (τ, t 0 + δ) , ∀t ∈ [t 0 , min (τ, t 0 + δ)] , (q n (t), q+ n (t)) ∈ K.
By a compactness argument, we have:

∃α > 0, ∀l ∈ {1, 2, • • • , L} , ∀q ∈ ∂K q l , ∃l , B(q, α) ⊂ K q l .
As a consequence, for n larger than N 0 , the interval [t 0 , min(τ, t 0 + δ)] is the disjoint union of a finite number, say N n , of intervals I ni such that:

∀i ∈ {1, 2, • • • , N n } , ∃l ∈ {1, 2, • • • , L} , q n (I ni ) ⊂ K q l .
Moreover, the intervals I ni can be constructed in such a way that:

∀n ≥ N 0 , N n ≤ 1 + 4V δ α .
Furthermore, recalling:

∀n, D q+ n = f q n , q+ n ; t dt + n i=1 λ ni dϕ i (q n (t)), (55)
where the λ ni are nonpositive real measures on [t 0 , min(τ, t 0 +δ)], and performing the same job as in the proof of proposition 18 (estimate 38), one obtains:

∀n ≥ N 0 , ∀i ∈ {1, 2, • • • N } , ∀j ∈ {1, 2, • • • N n } , I nj (-λ ni ) ≤ 2 √ λ max (4V ) + F + d 2 G(4V ) 2 2λ min δ.
There results:

∀n ≥ N 0 , ∀i ∈ {1, 2, • • • N } , ( 56
) [t 0 ,min(τ,t 0 +δ)] (-λ ni ) ≤ 1 + 4V δ α 2 √ λ max (4V ) + F + d 2 G(4V ) 2 2λ min δ
The measures λ ni are uniformly bounded with respect to n. By the equation of motion (55), we obtain that the real numbers Var ( q+ n ; [t 0 , min(τ, t 0 + δ)]) are uniformly bounded with respect to n, for n larger than N 0 . The assertion of step 2 is now a direct consequence of proposition 34.

Step 3. The function q lim constructed in step 2 satisfies the equation of motion: R lim = D q+ limf q lim , q+ lim ; t dt ∈ -N * (q lim ).

Moreover, the real measure R lim , q+ lim + qlim q lim is a nonpositive measure on the interval [t 0 , min(τ, t 0 + δ)[. Proof of step 3. We denote by M b ([a, b], R) the Banach space of all bounded real measures on an interval [a, b]. By estimate 56, we can find N bounded real measures λ ilim such that:

lim n→+∞ λ in = λ ilim in M b ([t 0 , min(τ, t 0 + δ)], R)
D q+ lim = f q lim , q+ lim ; t dt + n i=1 λ ilim dϕ i (q lim ),
the weak* topology being Hausdorff. Now, by formulae [START_REF] Moreau | Standard inelastic shocks and the dynamics of unilateral constraints, in Unilateral problems in structural analysis[END_REF] we have to prove:

Suppλ ilim ⊂ {t ∈ [t 0 , min(τ, t 0 + δ)]; ϕ i (q lim (t)) = 0} (57)
Consider ]a, b[⊂ [t 0 , min(τ, t 0 + δ)] such that:

∀s ∈]a, b[, ϕ i (q lim (s)) < 0.

The interval ]a, b[ is the union of the compact intervals K j = [a + 1/j, b -1/j] (j ∈ N * ). Fix j ∈ N * . For n large enough:

∀s ∈ K j , ϕ i (q n (s)) < 0, so λ in|K j = 0. We deduce:

∀g ∈ C 0 c (K j ; R), K j g dλ ilim = 0
Therefore, λ ilim|]a,b[ = 0 and this achieves the proof of inclusion (57) and therefore of the first assertion of step 3.

For the second assertion of step 3, we are going to prove actually:

∀t 1 , t 2 ∈ [t 0 , min(τ, t 0 + δ)[, t 1 < t 2 , ]t 1 ,t 2 ]
R lim , q+ lim + qlim q lim ≤ 0.

(58) Fix such t 1 , t 2 and arbitrary ε > 0. We have:

]t 1 ,t 2 ]
R lim , q+ lim + qlim q lim = q+ lim (t 2 )

2 q lim (t 2 ) -q+ lim (t 1 )

2 q lim (t 1 ) -2 t 2 t 1
f q lim (t), q+ lim (t); t , q+ lim (t) dt.

By the right-continuity of the function t → q+ lim (t 2 ) q lim (t) and the results of step 2, we can find t 1 , t 2 ∈ [t 0 , min(τ, t 0 + δ)[ (t 1 < t 2 ) and an integer N 0 such that:

t i ≤ t i ≤ t i + ε 8V F
, and

∀n ≥ N 0 , q+ lim (t i ) 2 q lim (t i ) -q+ n (t i ) 2 q n (t i ) ≤ ε 8 (i = 1, 2).
Moreover, by Lebesgue's Dominated Convergence Theorem, N 0 may be assumed large enough to ensure:

∀n ≥ N 0 , t 2 t 1 f q lim (t), q+ lim (t); t , q+ lim (t) -f q n (t), q+ n (t); t , q+ n (t) dt ≤ ε 8 .
It is easily deduced that:

∀n ≥ N 0 , ]t 1 ,t 2 ] R lim , q+ lim + q- lim q lim - ]t 1 ,t 2 ] R n , q+ n + q- n q n ≤ ε.
Since ε is arbitrary and ]t 1 ,t 2 ] R n , ( q+ n + qn ) q n is nonpositive (proposition 7), the conclusion (assertion (58)) follows.

Step 4. Consider an arbitrary instant t g ∈]t 0 , min(τ, t 0 + δ)[ such that:

(dϕ i (q lim (t g ))) i∈J(q lim (t g )) is orthogonal in T * q lim (t g ) Q.
Then, q lim satisfies the impact constitutive equation at instant t g :

q+ lim (t g ) = q- lim (t g ) -(1 + φ)Proj q lim (t g ) q- lim (t g ); N (q lim (t g )) .
Proof of step 4. Consider a local chart (U, ψ) centered at q lim (t g ) such that:

the cardJ(q lim (t g )) first components of ψ(q) are α i ϕ i (q) (i ∈ J(q lim (t g ))), where the α i are some fixed positive real constants. -∀q ∈ U, J(q) ⊂ J(q lim (t g )), the matrix (g ij ( lim (t g ))) is the identity matrix:

(g ij ( lim (t g ))) = δ ij
We have to prove:

q+i lim (t g ) = -φ q-i lim (t g ), 1 ≤ i ≤ CardJ(q lim (t g )) q+i lim (t g ) = q-i lim (t g ), CardJ(q lim (t g )) + 1 ≤ i ≤ d (59)
First, we are going to prove:

∀ε > 0, ∃N 0 , η > 0, ∀n ≥ N 0 , ∀t 1 , t 2 ∈ [t g -η, t g + η], t 1 < t 2 , q+i n (t 2 ) ≤ q+i n (t 1 ) + ε, (60) and ∀ε > 0, ∃N 0 , η > 0, ∀n ≥ N 0 , ∀t 1 , t 2 ∈ [t g -η, t g + η], t 1 < t 2 , ∀t ∈ [t 1 , t 2 ], q i n (t) < 0 =⇒ q+i n (t 2 ) -q+i n (t 1 ) ≤ ε .
(61) Fix ε > 0 arbitrary, and pick a positive real number η small enough and an integer N 0 large enough to ensure: ∀t ∈ [t gη, t g + η], ∀n ≥ N 0 , q lim (t) ∈ U and q n (t) ∈ U.

Let V be a positive real constant, large enough to majorize all the quantities q +i n (t) and Var q +i n ; [t gη, t g + η] ,

when i, t and n wander respectively in the sets {1, 2, • • • , d}, [t gη, t g + η] and {n ∈ N ; n ≥ N 0 }. We may assume that η is small enough and N 0 large enough to ensure:

∀t ∈ [t g -η, t g +η], ∀n ≥ N 0 , |g ij (q n (t)) -δ ij | ≤ min ε 4dV , ε 2 8dV 2 .
Multiplying the equation of motion (36) by ( q+i n + q+i n )/2 and integrating over ]t 1 , t 2 ], we obtain easily:

∀n ≥ N 0 , ∀t 1 , t 2 ∈ [t g -η, t g + η], t 1 < t 2 , 1 2 q+i n (t 2 ) 2 ≤ 1 2 q+i n (t 1 ) 2 + 1 2 ε 2 2 + t 2 t 1 F + 3 2 d 2 GV 2 q+i n (s) ds, which gives, ∀n ≥ N 0 , ∀t 1 , t 2 ∈ [t g -η, t g + η], t 1 < t 2 , q+i n (t 2 ) ≤ q+i n (t 1 ) + ε 2 + 2η F + 3 2 d 2 GV 2 ,
by lemma 17 and the desired conclusion (60) for sufficiently small η. For the second assertion (61), suppose we have in addition:

∀t ∈ [t 1 , t 2 ], q i n (t) < 0.
There results that λ in vanishes on [t 1 , t 2 ] and integration of the equation of motion (37) gives:

q+i n (t 2 ) -q+i n (t 1 ) ≤ ε 2 + 2η F + 3 2 d 2 GV 2 ,
and, therefore the desired conclusion (61) for sufficiently small η. Now, let us come back to the proof of assertions (59). Fix i ∈ {1, 2, • • • , d}. Only the following four cases are possible:

case 1: CardJ(q lim (t g )) + 1 ≤ i ≤ d case 2: 1 ≤ i ≤ CardJ(q lim (t g )) and q-i lim (t g ) = 0 case 3: 1 ≤ i ≤ CardJ(q lim (t g )), q-i lim (t g ) > 0 and φ = 0 case 4: 1 ≤ i ≤ CardJ(q lim (t g )), q-i lim (t g ) > 0 and φ > 0 We examine them successively.

-Case 1: CardJ(q lim (t g )) + 1 ≤ i ≤ d.

Fix ε > 0 arbitrary. By assertion (61), we can pick a positive real number η small enough and an integer N 0 large enough to ensure:

∀n ≥ N 0 , ∀t 1 , t 2 ∈ [t g -η, t g + η], t 1 < t 2 , q+i n (t 2 ) -q+i n (t 1 ) ≤ ε, since ∀t ∈ [t g -η, t g + η], ∀n ≥ N 0 , q i n (t) < 0,
q+i n (t k ) = q+i lim (t k ) (k = 1, 2)
and, therefore, N 0 can be assumed large enough to ensure:

∀n ≥ N 0 , q+i n (t k ) -q+i lim (t k ) ≤ ε (k = 1, 2)
. Then, we have:

q+i lim (t g ) -q-i lim (t g ) ≤ q+i lim (t g ) -q-i lim (t 2 ) + q+i lim (t 2 ) -q-i n (t 2 ) + q+i n (t 2 ) -q-i n (t 1 ) + q+i n (t 1 ) -q-i lim (t 1 ) + q+i lim (t 1 ) -q-i lim (t g ) ≤ 5ε.
Since ε is arbitrary, we get the desired conclusion: q+i lim (t g ) = q-i lim (t g ) -Case 2: 1 ≤ i ≤ CardJ(q lim (t g )) and q-i lim (t g ) = 0. Fix ε > 0 arbitrary. By assertion (60), we can pick a positive real number η small enough and an integer N 0 large enough to ensure:

∀n ≥ N 0 , ∀t 1 , t 2 ∈ [t g -η, t g +η], t 1 < t 2 ,
q+i n (t 2 ) ≤ q+i n (t 1 ) +ε. Exactly as in case 1, η is assumed sufficiently small to ensure:

∀t ∈ [t g -η, t g [, q+i lim (t) ≤ ε, ∀t ∈]t g , t g + η],
q+i lim (t) -q+i lim (t g ) ≤ ε, and N 0 large enough to have:

∀n ≥ N 0 , q+i n (t k ) -q+i lim (t k ) ≤ ε (k = 1, 2), for some t 1 ∈ [t g -η, t g [ and some t 2 ∈]t g , t g + η]. We get: q+i lim (t g ) ≤ q+i lim (t g ) -q-i lim (t 2 ) + q+i lim (t 2 ) -q+i n (t 2 ) + q+i n (t 2 ) , ≤ q+i n (t 1 ) + 3ε, ≤ 5ε,
which gives the desired conclusion q+i lim (t g ) = 0, since ε is arbitrary.

-Case 3: 1 ≤ i ≤ CardJ(q lim (t g )), q-i lim (t g ) > 0 and φ = 0. Fix ε arbitrary in ]0, q-i lim (t g )/16]. We pick η and N 0 such that both assertions (60) and (61) hold. Actually, η is assumed small enough to ensure:

∀t ∈ [t g -η, t g [, q i lim (t) t -t g -q-i lim (t g ) ≤ ε, ∀t ∈ [t g -η, t g [, q+i lim (t) -q-i lim (t g ) ≤ ε, ∀t ∈]t g , t g + η], q+i lim (t) -q+i lim (t g ) ≤ ε,
and, by step 2, N 0 is assumed large enough to get:

∀n ≥ N 0 , ∀t ∈ [t g -η, t g + η], q i n (t) -q i lim (t) ≤ ηε, ∀n ≥ N 0 , q+i n (t 1 ) -q+i lim (t 1 ) ≤ ε, ∀n ≥ N 0 , q+i n (t 2 ) -q+i lim (t 2 ) ≤ ε, for some fixed t 1 ∈ [t g -η/2, t g -η/4] and t 2 ∈ [t g + 3η/4, t g + η].
From these inequalities, it is easily deduced:

- 17 16
η 2 q-i lim (t g ) ≤ q i lim (t 1 ) ≤ - 15 16
η 4 q-i lim (t g ),
and therefore, ∀n ≥ N 0 , -10 16 η q-i lim (t g ) ≤ q i n (t 1 ) ≤ -

2 16 η q-i lim (t g ). (62)
Furthermore,

q+i n (t 1 ) ≥ q+i lim (t 1 ) -2ε ≥ 14 16 q-i lim (t g ). (63)
Then, by estimates (62) and (63) and assertion (61), it is readily seen that:

∀n ≥ N 0 , ∃t n ∈]t 1 , t 1 + η[, q i n (t n ) = 0.
But, since φ = 0, we have:

∀n ≥ N 0 , q+i n (t n ) = 0,
and, therefore,

∀n ≥ N 0 , q+i n (t 2 ) ≤ ε,
by assertion (60). We deduce:

q+i lim (t g ) ≤ 3ε,
and the desired conclusion q+i lim (t g ) = 0, since arbitrarily small ε can be chosen.

-Case 4: 1 ≤ i ≤ CardJ(q lim (t g )), q-i lim (t g ) > 0 and φ > 0. Fix ε arbitrary in ]0, φ q-i lim (t g )/16]. We pick η, N 0 , t 1 and t 2 exactly in the same way as for case 3. As in step 3, we have:

∀n ≥ N 0 , ∃t n ∈]t 1 , t 1 + η[, q i n (t n ) = 0,
but, here, it is readily seen that t n is the unique instant in [t 1 , t g + η] such that q i n (t n ) = 0. Now, we obtain:

q+i lim (t g ) + φ q-i lim (t g ) ≤ q+i lim (t g ) -q+i n (t 2 ) + q+i n (t 2 ) -q+i n (t n ) + φ q-i n (t n ) -q+i n (t 1 ) + φ q+i n (t 1 ) -q-i lim (t g ) ≤ 6ε
, by use of assertion (60). Since ε can be arbitrarily small, we have the desired conclusion: q+i lim (t g ) = -φ q-i lim (t g ). This achieves the proof of step 4.

Step 5. Conclusion of the proof of theorem 20

First, we are going to prove: ∀t ∈ [t 0 , min(τ, t 0 + δ)], q lim (t) = q(t).

(64) Define:

t 1 = sup {t ∈ [t 0 , min(τ, t 0 + δ)]; ∀s ∈ [t 0 , t], q lim (s) = q(s)} .
Notice that the set in the above definition is non empty, since it contains t 0 . By continuity, we have: ∀t ∈ [t 0 , t 1 ], q lim (t) = q(t). Now , assume:

t 1 < min(τ, t 0 + δ).
By the assumption made on q in the theorem and by step 4, the function q lim is readily seen to satisfy the impact constitutive equation at instant t 1 . Therefore, q lim (t 1 ), q+ lim (t 1 ) = q(t 1 ), q+ (t 1 ) .

Furthermore, we have seen in step 3 that q lim satisfies the equation of motion and that R lim , q+ lim + qlim q lim is a nonpositive real measure. But, the proof of local uniqueness (theorem 8) uses nothing more than that. We deduce that there exists a right-neighbourhood of t 1 on which the functions q lim and q coincide identically. But, this contradicts the definition of t 1 and achieves the proof of assertion (64). As a result, the function q lim is uniquely determined and the conclusions of step 2 are valid not only for a subsequence but for the whole sequence (q n ). Now, if t 0 + δ < τ, we pick t 0 ∈ [t 0 + δ/2, t 0 + δ[ such that:

lim n→+∞ q n (t 0 ), q+ n (t 0 ) = q(t 0 ), q+ (t 0 ) .
Performing the same job for instant t 0 instead of t 0 , we extend the conclusion to interval [t 0 , min(τ, t 0 + 3δ/2)]. Processing so inductively a large enough number of times, we obtain the desired conclusion.

Remark. A straightforward modification of the proof of step 4 shows that the conclusions of theorem 20 hold if we only assume that φ is continuous and constant on each fiber:

∀q, ∀v 1 , v 2 ∈ T q Q, φ(q, v 1 ) = φ(q, v 2 ).
The conclusions of theorem 20 also hold if φ is only assumed continuous and if, moreover, we have:

∀t ∈ [t 0 , T [, CardJ(q(t)) ≤ 1.

Indications on the numerical computation of the solutions

Consider the problem P described in section 3.4. Assume furthermore, for sake of simplicity, that the impact function φ is constant. The maximal solution associated with the initial condition (q 0 , v 0 ) at time t 0 = 0, is denoted by (T m , q). We consider a local chart (U, ψ) at q 0 and a positive real number T such that:

∀t ∈ [0, T ], q(t) ∈ U.
By assumption 20, we may assume: ∀q ∈ U, (dϕ i (q)) i∈J(q) is linear independent in T * q Q, taking a smaller U if necessary. We consider a sequence of approximants, defining for every n ≥ 1:

• h n = 2 -n T, • t n,k = kh n = k2 -n T (k = 0, 1, 2, • • • , 2 n ), • (q n,0 , v n,0 ) = (q 0 , v 0 ), • q n,k = q n,k-1 + h n v n,k-1 (k = 1, 2, • • • , 2 n ), • v α n,k = v α n,k-1 + g αβ (q n,k )f β (q n,k , v n,k-1 ; t n,k ) -Γ α βγ (q n,k )v β n,k-1 v γ n,k-1 h n (k = 1, 2, • • • , 2 n , α = 1, 2, • • • , d), • v n,k = v n,k -(1 + φ) Proj q n,k [v n,k , N(q n,k )] (k = 1, 2, • • • , 2 n ), • v n (t) = v n,k , if t ∈ [t n,k , t n,k+1 [ with k = 0, 1, • • • 2 n -1, v n,2 n , if t = T = t n,2 n , • q n (t) = q 0 + t 0 v n (s) ds.
Actually, there may happen that the function q n cannot be defined on [0, T ] if there exists an integer k n such that q n,k n +1 ∈ U . In such a case, the function q n is defined only on [0, t n,k n ].

This type of algorithm was introduced by Moreau and used without further justifications. It should be stressed that one cannot hope that the sequence of approximants (q n ) converges in general towards the solution q, since continuous dependence on initial condition does not hold in general. Actually, it is easy to build an explicit example, in the spirit of the example of section 7, where the sequence (q n ) does not converge pointwisely towards any function at all. However, in the special case where all the multiple impacts are orthogonal, things behave nicely and we have:

Theorem 21 Suppose that the solution q is such that all multiple impacts are orthogonal:

∀t ∈ [0, T ], (dϕ i (q(t))) i∈J(q(t)) is orthogonal in T * q(t) Q,
(with the convention that the empty set is orthogonal). Then, there exists an integer N 0 such that the function q n is well-defined on [0, T ] for n ≥ N 0 . Moreover, the sequence (q n ) converges uniformly on [0, T ] towards q (or more precisely towards ψ(q)).

Theorem 21 can be proved along the same steps as those of the proof of theorem 20. The necessary adaptation of the details is left to the reader.

Appendix: the class of motion M M A(I; Q)

In this section, we are going to define the concept of vector field with bounded variation along a locally absolutely continuous curve on a riemannian manifold. The definition and basic properties of absolutely continuous functions and functions with bounded variation from a real interval to a finite-dimensional normed vector space are supposed to be known. The reader is refered to Rudin [START_REF] Rudin | Real and complex analysis[END_REF] and Moreau [START_REF] Moreau | Bounded variation in time[END_REF]. These concepts are intimately connected with measure theory. Two expositions of measure theory compete: the set-theoretic approach (see for example Rudin [START_REF] Rudin | Real and complex analysis[END_REF]) and the duality approach (see for example Bourbaki [START_REF] Bourbaki | Intégration[END_REF]). Both approach are connected by Riesz's representation theorem. In this paper, we stick to the duality approach. If I is a real interval and E a real finite-dimensional normed vector space, C 0 c (I; E) will denote the space of continuous functions from I to E with compact support. A measure on I with values in E (resp. E * ) will be any linear form μ on When the constant M a,b can be found independent on a and b, the measure μ is said bounded. For all what concerns measure theory, the reader is refered to Bourbaki [START_REF] Bourbaki | Intégration[END_REF] where he will note the definition of the support Supp μ of a measure μ (Bourbaki [START_REF] Bourbaki | Intégration[END_REF], p. 64).

The following list of definitions and propositions aims at carrying these concepts over riemannian manifolds. The cornerstone is, of course, the identification of tangent spaces at different points of a curve by means of parallel translation.

This appendix is also an occasion to state precisely the classical theorems which are used in this paper.

Definition 22 Let I be a real interval and q : I → Q a curve on Q. q is said locally absolutely continuous if, for all t in I, there exists a compact neighbourhood J of t in I and a chart (U, ψ) such that:

q(J) ⊂ U , ψ • q : J → R d is absolutely continuous.

Since Q can be covered by a countable collection of chart domains, there results from Lebesgue's theorem that q(t) admits a tangent vector q(t) ∈ T q(t) Q for dt-almost all t in I where dt denotes the Lebesgue measure on the real line (and also its restriction on I). The riemannian structure on Q and the Cauchy-Lipschitz-Caratheodory theorem allow us to define classically a parallel translation operator along q, τ t,s : T q(s) Q → T q(t) Q (see, for example, Chavel [START_REF] Chavel | Riemannian Geometry: a Modern Introduction[END_REF], p. 7). τ t,s is defined for all (s, t) ∈ I 2 .

Definition 23 Let q be a locally absolutely continuous curve from I to Q. One calls vector (resp. 1-form) field on q(t) any mapping X (resp. X * ) from I to T Q (resp. T * Q) such that: ∀t ∈ I, Π Q (X(t)) = q(t) (resp. Π * Q (X * (t)) = q(t)).

A vector field X on q(t) (resp. a 1-form field X * on q(t)) will be said locally absolutely continuous (resp. absolutely continuous, resp. locally with bounded variation, resp. with bounded variation) if there exists t 0 in I such that the mapping:

θ t 0 I → T q(t 0 ) Q s → τ t 0 ,s (X(s)) resp. θ * t 0 I → T * q(t 0 ) Q s → • τ t 0 ,s ( • X * (s))
, is locally absolutely continuous (resp. absolutely continuous, resp. locally with bounded variation, resp. with bounded variation). If X has bounded variation on I, its variation over I is by definition:

Var(X(s); I) = Var(τ t 0 ,s (X(s)); I). (65)

From the identity: ∀s 1 , s 2 , t 1 , t 2 ∈ I, τ t 1 ,s 1 (X(s 1 ))τ t 1 ,s 2 (X(s 2 )) q(t 1 ) = τ t 2 ,s 1 (X(s 1 ))τ t 2 ,s 2 (X(s 2 )) q(t 2 ) , it is easily deduced that the above definition is independent on a particular choice of t 0 and so is the real number Var(X(s); I).

The covariant derivative of a locally absolutely continuous vector field X along q can be defined for dt-almost all t in I by: DX(t) dt = d ds (τ t,s (X(s)))| s=t for dt-a.a. t ∈ I.

Definition 24 Let (I, q) be a continuous curve on Q. We denote by C 0 c (I, q; T Q) (resp. C 0 c (I, q; T * Q)) the space of continuous functions ϕ from I to T Q (resp. T * Q) with compact support and such that: ∀t ∈ I, Π Q (ϕ(t)) = q(t) (resp. Π * Q (ϕ(t)) = q(t)).

One calls measure on the curve (I, q) taking values in T Q (resp. T * Q) any linear form μ on C 0 c (I, q; T * Q) (resp. C 0 c (I, q; T Q)) enjoying the following continuity property: The real number μ(ϕ) will also be denoted by I ϕ(t), dμ q(t) .

Proposition 25 Let (I, q) a continuous curve on Q and μ a measure on q taking values in T * Q. For any nonnegative function f of C 0 c (I; R), one defines: |μ| (f ) = sup g∈C 0 c (I,q;T Q) g(t) q(t) ≤f (t) I g(t) dμ , where the supremum is finite thanks to the continuity properties included in the definition of measures. For arbitrary f in C 0 c (I; R), one defines:

|μ| (f ) = |μ| f + -|μ| f -,
where x ± = max{±x, 0} are the classical positive and negative parts. Then, the functional |μ| is a real measure called the modulus measure of μ.

The proof is omitted since it is completely identical to the proof of the similar statement for complex measures (see Bourbaki [START_REF] Bourbaki | Intégration[END_REF], p. 54).

The support Supp μ of a measure μ on q(t) taking values in T * Q is, by definition, the support Supp |μ| of its modulus measure.

We define L 1 loc (I, q, |μ| ; T * Q) by the space of functions θ defined for |μ|almost all t in I, taking values in T * Q and such that:

-Π * Q (θ(t)) = q(t) for |μ|-almost all t ∈ I, -∀ϕ ∈ C 0 c (I, q; T Q), t → ϕ(t), θ(t) q(t) ∈ L 1 (I, |μ| ; R).

Proposition 26 Let μ be a measure on q(t) taking values in T * Q. Then, there exists a unique (class of) function l μ ∈ L 1 loc (I, q, d |μ| ; T * Q) such that: -Π * Q (l μ (t)) = q(t) for d |μ|-almost all t ∈ I,

-∀ϕ ∈ C 0 c (I, q; T Q), I ϕ(t), dμ q(t) = I ϕ(t), l μ (t) q(t) d |μ|. This fact will be denoted by: dμ = l μ d |μ|. We shall say that l μ is the density of measure μ with respect to measure |μ|.

Proof. For measure taking values in a finite-dimensional vector space, the above statement is a classical direct consequence of Lebesgue-Radon-Nikodym theorem (see Rudin [START_REF] Rudin | Real and complex analysis[END_REF]). It is readily carried over manifolds by means of a locally finite partition of unity modelled on chart domains.

Definition 27 Let X be a vector field with locally bounded variation on an absolutely continuous curve (I, q) and t 0 an arbitrary element of I. We denote by d t 0 X the Stieljes measure (see Moreau [START_REF] Moreau | Bounded variation in time[END_REF]) associated to the mapping with locally bounded variation:

θ t 0 I → T q(t 0 ) Q s → τ t 0 ,s (X(s)) .
For Y ∈ C 0 c (I, q; T Q) and Y * ∈ C 0 c (I, q; T * Q), the linear forms:

Y → I (τ t 0 ,s (Y (s)), d t 0 X) q(t 0 ) and Y * → I (τ t 0 ,s ( • Y * (s)), d t 0 X) q(t 0 )
turns out to be independent on a particular choice of t 0 and define measures on q taking respectively values in T * Q and T Q. They are denoted by DX and DX and called the covariant and contravariant representative of the covariant Stieljes measure associated with X.

Proposition 28 If X is a locally absolutely continuous vector field on a locally absolutely continuous curve from I to Q then:

DX = DX dt dt and DX = DX dt dt (66)
Reciprocally, if X is locally with bounded variation and such that its covariant Stieljes measure DX admits a density with respect to the Lebesgue measure, then X is locally absolutely continuous and relations (66) hold.

Proof. This is an immediate consequence of definition 27 and of the similar statement for functions taking values in a finite-dimensional normed vector space.

Proposition 28 ensures the consistency of our notations. Let us now turn to practical calculations in charts.

Proposition 29 Let (U, ψ) be a chart on Q, (I, q) an absolutely continuous curve on Q such that q(I) ⊂ U and X a vector field on (I, q). The components (X i ) (i = 1, 2, • • • d) of X in the natural chart of T Q associated with ψ are real functions defined on I. The vector field X is locally absolutely continuous (resp. absolutely continuous, resp. locally with bounded variation, resp. with bounded variation) if and only if every function X i is locally absolutely continuous (resp. absolutely Proposition 33 Let q ∈ M M A(I; Q) with q(I) ⊂ U domain of a chart. Then, D q+ = d ∂K(q(t), q+ (t)) ∂ q+i -∂K(q(t), q+ (t)) ∂q i dt e i (q(t)).

Proof. Reproduce the proof of proposition 2 with the help of proposition 29.

Proposition 34 Let (q n ) n∈N a sequence of elements of M M A([0, T ]; Q) such that:

there exists a compact subset K of T Q such that:

∀n ∈ N, ∀t ∈ [0, T ], (q n (t), q+ n (t)) ∈ K, -∃M > 0, ∀n ∈ N, Var ( q+ n ; [0, T ]) ≤ M. Then, there exists a subsequence of (q n ) n∈N , also denoted by (q n ) n∈N , such that:

-(q n ) n∈N converges uniformly on [0, T ] for the riemannian metric towards a function q lim belonging to M M A([0, T ]; Q), -The sequence (q n (t), q+ n (t)) converges towards q lim (t), q+ lim (t) in T Q for almost all t in ]0, T [. Proof. This is a generalization of Helly's theorem to the case of a riemannian manifold. The set K = Π Q (K) being compact, there exists ε > 0 such that (cf Chavel [START_REF] Chavel | Riemannian Geometry: a Modern Introduction[END_REF], p. 23):

for all q in K , B(q, ε) (= {q ∈ K ; d(q, q ) < ε}) is the domain of a chart ψ q , for all q in K , the distance defined by |ψ q (q 1 )ψ q (q 2 )| and the riemannian distance d are equivalent on B(q, ε).

First, we extract a subsequence of (q n ), also denoted by (q n ), such that: lim n→+∞ (q n (0), q+ n (0)) = (q 0 , v 0 ) in T Q, and there exists N 0 ∈ N large enough to have:

∀n ≥ N 0 , d(q 0 , q n (0)) < ε 2 ,
Now, by: ∀t ∈ [0, T [, ∀n ∈ N, q+ n (t) q n (t) ≤ q+ n (0) q n (0) + Var q+ n ; [0, T ] , (67) there exists α 0 (0 < α 0 ≤ T ) small enough to have: ∀t 0 ∈ [0, T ], ∀t ∈ [t 0 , min(T, t 0 + α 0 )], ∀n ∈ N, d(q n (t), q n (t 0 )) < ε 2 .

Then, it is easily checked that the functions ψ q 0 (q n (t)) |[0,α 0 ] (n ≥ N 0 ) satisfy the hypothesis of Helly's theorem and therefore the conclusion of the proposition holds on [0, α 0 ].

The dynamics of discrete mechanical systems with perfect unilateral constraints 73

Now, choose t 1 ∈ [α 0 /2, α 0 ] such that: lim n→+∞ q n (t 1 ), q+ n (t 1 ) = q lim (t 1 ), q+ lim (t 1 ) in T Q,

and N 1 large enough to have: ∀n ≥ N 1 , d(q lim (t 1 ), q n (t 1 )) < ε 2 .

Performing the same job as above on the chart of domain B(q lim (t 1 ), ε), one obtains that the conclusion of the proposition holds on [0, min(T, 3α 0 /2)]. Processing so inductively a large enough number of times, we obtain the desired conclusion.

Remark. If the riemannian manifold Q is assumed to be complete, the first hypothesis in proposition 34 can be weakened and replaced by:

there exists a compact subset K 0 of T Q such that: ∀n ∈ N, (q n (0), q+ n (0)) ∈ K 0 , Indeed, this hypothesis allows us to extract a subsequence of (q n ) such that: lim n→+∞ (q n (0), q+ n (0)) = (q 0 , v 0 ) in T Q.

By estimate 67, we get:

∃D > 0, ∀t ∈ [0, T ], ∀n ∈ N, q+ n (t) q n (t) ≤ D, and d (q 0 , q n (t)) ≤ D.

The riemannian manifold Q being complete, by Hopf-Rinow theorem (cf Chavel [START_REF] Chavel | Riemannian Geometry: a Modern Introduction[END_REF], p. 26), the functions (q n , q+ n (t)) take values in a compact subset K of T Q.

Fig. 1 .

 1 Fig. 1. Geometry of the double pendulum.

Fig. 2 .

 2 Fig.2. Motion of a punctual particle subjected to gravity and bouncing on the floor.

Fig. 4 .

 4 Fig. 4. Newton's balls experiment.

  C 0 c (I; E * ) (resp. C 0 c (I; E)) satisfying the following continuity property: ∀a, b ∈ I, a < b ∃M a,b ≥ 0, ∀ϕ with Supp ϕ ⊂ [a, b], |μ(ϕ)| ≤ M a,b max t∈I ϕ(t) .

  ∀a, b ∈ I, a < b ∃M a,b ≥ 0, ∀ϕ with Supp ϕ ⊂ [a, b], |μ(ϕ)| ≤ M a,b max t∈I ϕ(t) q(t) .

  Proposition 3 Problems E Q and E S are equivalent: any solution of E Q generates a solution of E S and reciprocally. Moreover, if Q and the functions ϕ i are of class C p (p ≥ 2), and f of class C p-1 then the unique maximal solution of E Q and E S is of class C p . If Q, f and the ϕ i are analytic functions then so is the maximal solution of E Q and E S .

  by the choice of the chart we made. Actually, η can be assumed small enough to ensure:∀t ∈ [t gη, t g [, q+i lim (t) -q-i lim (t g ) ≤ ε, ∀t ∈]t g , t g + η],q+i lim (t) -q+i lim (t g ) ≤ ε, by proposition 32. By step 2, we can find t 1 ∈ [t gη, t g [ and t 2 ∈]t g , t g + η] such that: lim

	n→+∞

Acknowledgements. I would like to thank particularly Professor Michelle Schatzman for her always relevant remarks. The final form of this article owes very much to her circumspect and disinterested advices.

If we still assume in this case that f is nonnegative, then it is easy to see that the only possible motion is immobility. Indeed, if, ∃t 2 , q(t 2 ) < 0, define t 1 = inf {t ∈ R + ; ∀s ∈]t, t 2 ] q(s) < 0}. Then, by continuity of q: t 1 < t 2 and q(t 1 ) = 0. By the completely inelastic impact constitutive equation, one gets: v(t 1 ) = 0, and, so:

Nevertheless, we are going to construct an example similar to example 2, which shows that even in case of the completely inelastic impact constitutive equation and f of class C ∞ , we can obtain multiple solutions for the corresponding problem P. Of course, f should not be of constant sign.

The function f is searched under the form:

where

) n∈N * are positive real sequences which are to be determined. We demand:

These sequences are to be determined in such a way that the corresponding problem P admits two distinct solutions v I and v II . We demand that v I , v II and the corresponding functions q I , q II are such that :

where (q n ) n∈N * and (v n ) n∈N * are positive real sequences which are to be determined.

Consider the time interval

, the position of a particle which is at q = -q n+1 with velocity v = v n+1 at time t = 1 n+1 should increase from -q n+1 to 0. This is written as:

continuous, resp. locally with bounded variation, resp. with bounded variation). Moreover, in such a case, we have:

Proof. This is an immediate consequence of definition 27.

Proposition 30 Let X be a vector field with locally bounded variation of an absolutely continuous curve (I, q). Then, for any t 0 in I, the two limits lim t→t - 0 X(t) and lim t→t + 0 X(t) exist in T Q and are such that:

These limits are denoted by X -(t 0 ) and X + (t 0 ) and can be different only on an at most countable subset of I. The mapping:

has locally bounded variation and:

Proof. It is a direct consequence of the similar statement for functions taking values in euclidean R d (see Moreau [START_REF] Moreau | Bounded variation in time[END_REF]) and of definition 27.

Definition 31

We denote by M M A(I; Q) (motions with measure acceleration) the set of all locally absolutely continuous motions q(t) from I to Q such that the right velocity q+ (t) exists for all t in I and defines a vector field with locally bounded variation on q(t).

Proposition 32 Let q be in M M A(I; Q). Then, q+ : I → T Q is right continuous: ∀t ∈ I, q+ (t) + = q+ (t).

Moreover, q(t) admits a left velocity vector at each instant of I and:

∀t ∈ I, q-(t) = q+ (t) -.

Proof. Use the Mean Value Inequality in a local chart.