Licence et Magist`ere de Physique

Universite Joseph Fourier

Cours de M
ecanique Analytique
Jonathan Ferreira

Annee Universitaire 2002-2003

Laboratoire dAstrOphysique de Grenoble

http://www-laog.obs.ujf-grenoble.fr/ferreira

Contents
1 M
ecanique de Lagrange
1.1 Coordonnees generalisees . . . . . . . . . . . . . .
1.2 Equations de la dynamique . . . . . . . . . . . . .
1.2.1 Principe de dAlembert . . . . . . . . . . .
1.2.2 Equations de Lagrange . . . . . . . . . . . .
1.2.3 Exemple 1: le pendule . . . . . . . . . . . .
1.2.4 Exemple 2: masse sur une tige avec ressort
1.3 Lois de conservation . . . . . . . . . . . . . . . . .
1.3.1 Variables cycliques . . . . . . . . . . . . . .
1.3.2 Lagrangien independant du temps . . . . .
1.3.3 Theor`eme de Noether . . . . . . . . . . . .
1.4 Une application: force centrale entre deux corps . .
1.4.1 Invariance par translation dans le temps . .
1.4.2 Invariance par translation dans lespace . .
1.4.3 Invariance par rotation dans lespace . . . .
1.4.4 Loi horaire . . . . . . . . . . . . . . . . . .
1.4.5 Trajectoire . . . . . . . . . . . . . . . . . .
1.4.6 Force en 1/r2 , Loi de Kepler . . . . . . . .
1.5 Petites oscillations . . . . . . . . . . . . . . . . . .
1.5.1 Syst`emes `a 1 degre de liberte . . . . . . . .
1.5.2 Syst`emes `a n degres de liberte . . . . . . .
1.5.3 Oscillations forcees . . . . . . . . . . . . . .

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

1
1
3
3
4
7
7
7
8
8
10
10
11
11
11
12
12
13
13
13
14
16

2 Principe variationnel
2.1 Le principe de Hamilton . . . . . . . . . . . . . . . . . . . . . . .
2.2 Deduction des equations de Lagrange . . . . . . . . . . . . . . . .
2.3 Exemples simples de calcul variationnel . . . . . . . . . . . . . .
2.3.1 Plus petite distance dans un plan . . . . . . . . . . . . . .
2.3.2 La brachistochrone . . . . . . . . . . . . . . . . . . . . . .
2.4 Generalisation des equations de Lagrange . . . . . . . . . . . . .
2.4.1 Forces non conservatives . . . . . . . . . . . . . . . . . . .
2.4.2 Contraintes non holonomes: multiplicateurs de Lagrange .
2.5 Expressions du lagrangien en fonction de lespace-temps . . . . .
2.5.1 Mecanique non relativiste . . . . . . . . . . . . . . . . . .
2.5.2 Mecanique relativiste . . . . . . . . . . . . . . . . . . . . .
2.5.3 Remarques epistemologiques . . . . . . . . . . . . . . . .

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

19
19
20
21
21
21
22
22
23
25
25
27
29

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

ii

CONTENTS

3 M
ecanique de Hamilton
3.1 Hamiltonien dun syst`eme . . . . . . . . . . . . . . . . . . . . .
3.2 Equations canoniques de Hamilton . . . . . . . . . . . . . . . .
3.3 Principe variationnel . . . . . . . . . . . . . . . . . . . . . . . .
3.4 Etude dun cas simple: pendule 1D . . . . . . . . . . . . . . . .
3.4.1 Ecriture de lhamiltonien . . . . . . . . . . . . . . . . .
3.4.2 Le portrait de phase . . . . . . . . . . . . . . . . . . . .
3.4.3 Etude au voisinage de points particuliers . . . . . . . . .
3.4.4 Remarques dordre general . . . . . . . . . . . . . . . .
3.5 Theorie de Hamilton-Jacobi . . . . . . . . . . . . . . . . . . . .
3.6 Transformations canoniques . . . . . . . . . . . . . . . . . . . .
3.6.1 Fonctions generatrices . . . . . . . . . . . . . . . . . . .
3.6.2 Quelques transformations canoniques remarquables . . .
3.7 Les crochets de Poisson . . . . . . . . . . . . . . . . . . . . . .
3.7.1 Denition . . . . . . . . . . . . . . . . . . . . . . . . . .
3.7.2 Proprietes . . . . . . . . . . . . . . . . . . . . . . . . . .
3.7.3 Invariance canonique . . . . . . . . . . . . . . . . . . . .
3.7.4 Interpretation geometrique . . . . . . . . . . . . . . . .
3.8 Lespace des phases . . . . . . . . . . . . . . . . . . . . . . . . .
3.8.1 Flot hamiltonien . . . . . . . . . . . . . . . . . . . . . .
3.8.2 Incompressibilite du ot . . . . . . . . . . . . . . . . . .
3.8.3 Theor`eme de Liouville: lien avec la physique statistique
3.9 Syst`emes integrables . . . . . . . . . . . . . . . . . . . . . . . .
3.9.1 Theor`eme de Arnold-Liouville . . . . . . . . . . . . . . .
3.9.2 Cartes et atlas symplectiques . . . . . . . . . . . . . . .

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

31
31
32
33
33
33
34
34
35
35
37
37
38
39
39
40
41
42
44
44
45
47
48
48
49

4 Syst`
emes hamiltoniens
4.1 Lequation de Hamilton-Jacobi . . . . . . . . . . . . . .
4.1.1 La fonction principale de Hamilton . . . . . . . .
4.1.2 Laction hamiltonienne . . . . . . . . . . . . . . .
4.1.3 Methode generale de resolution . . . . . . . . . .
4.1.4 Methode de separation des variables . . . . . . .
4.1.5 Applications `a quelques probl`emes simples . . . .
4.1.6 Le principe de Maupertuis . . . . . . . . . . . . .
4.1.7 Mecanique ondulatoire de Louis de Brooglie . . .
4.2 Variables canoniques angles-actions . . . . . . . . . . . .
4.2.1 Syst`emes fermes periodiques . . . . . . . . . . . .
4.2.2 Variables angulaires . . . . . . . . . . . . . . . .
4.2.3 Variables dactions . . . . . . . . . . . . . . . . .
4.2.4 Fonction generatrice des variables angles-actions
4.2.5 Applications . . . . . . . . . . . . . . . . . . . .

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.

51
51
51
52
53
53
55
58
59
63
63
64
65
67
68

5 Description lagrangienne des milieux continus

5.1 Exemple dun passage `a la limite continue . . .
5.1.1 Corde elastique 1D . . . . . . . . . . . .
5.1.2 Retour au lagrangien . . . . . . . . . . .
5.2 Formulation lagrangienne des milieux continus
5.2.1 Conjectures . . . . . . . . . . . . . . . .
5.2.2 Equations de Lagrange du champ . . . .
5.3 Theorie classique des champs . . . . . . . . . .
5.3.1 Cadre general . . . . . . . . . . . . . . .
5.3.2 Exemple: electrodynamique classique .
5.3.3 Tenseur energie-impulsion dun champ .

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

73
73
73
74
75
75
76
78
78
79
81

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

CONTENTS
5.4

Complements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.4.1 Formulation relativiste de la theorie des champs . . . . . . . . . . . . . . . . . . . . .
5.4.2 Densite dhamiltonien . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

iii
82
82
83

iv

CONTENTS

Chapter 1

M
ecanique de Lagrange
1.1

Coordonn
ees g
en
eralis
ees

La mecanique de Newton se base sur trois postulats:

1. Principe dinertie: le mouvement dun corps isole est rectiligne uniforme dans un referentiel galileen.
2. Principe de la dynamique, orant une denition de la force
p
= m
r = F

3. Principe daction et de la reaction.
A laide de ces trois principes, la mecanique de Newton a montre sa puissance de description dans de
nombreux cas. Le mouvement dun syst`eme quelconque de N particules est ainsi obtenu par la resolution
de N equations vectorielles dierentielles du 2eme ordre, mettant en jeu 6N constantes dintegrations, correspondant aux positions et vitesses initiales des N particules.
Par ailleurs, cette description du reel suit un autre principe, commun `a toute la physique et quon peut
appeler principe de relativite. Ce principe stipule que les lois de la physique doivent etre independantes
de lobservateur, ce qui se traduit par une invariance de la forme des equations lors dun changement de
referentiel.
Il y a cependant des circonstances o`
u lapplication de la mecanique de Newton est delicate. Cest lorsquun
syst`eme poss`ede des contraintes internes (dues `a des forces de liaison), limitant le mouvement du syst`eme et
diminuant ainsi ses degres de liberte.
Exemple 1: corps rigide
Dans un corps indeformable, la distance entre deux points doit rester constante, cest `a dire (ri rj )2 = c2ij .
Exemple 2: pendule
Pendule de longueur l, bougeant dans le plan. Ses coordonnees obeissent `a la contrainte x2 + y 2 = l2 : il y a
donc 2 1 = 1 seul degre de liberte du syst`eme, langle . Par ailleurs, ce pendule peut devenir parametrique
si l = l(t) impose par lexterieur (ex: encensoir de Compostelle).
Exemple 3: perle sur un cerceau
Cerceau tournant avec une vitesse angulaire imposee, perle glissant sur le cerceau. La position de la perle
est reperee par les coordonnees
x

R sin cos

y
z

=
=

R sin sin
R cos
1


CHAPTER 1. MECANIQUE
DE LAGRANGE

Le seul degre de liberte de la perle est . Dans cet exemple, on suppose evidemment que la perle ne deforme
pas le cerceau.
Exemple 4: disque vertical roulant sans glisser
Soit un disque de rayon R roulant sans glisser sur un plan horizontal, de vitesse v constante. Ceci ne peut
evidemment se produire que si une force maintient un contact avec le sol, mais sans introduire un frottement.
On rep`ere la position de son centre par les deux coordonnees x et y, le plan du disque par langle que fait
laxe de rotation avec laxe Ox et un point M du disque par un angle : 4 coordonnees (x,y,, ). La vitesse
v du centre du disque verie
x =

R cos

R sin

On obtient ainsi les contraintes

La condition de roulement sans glissement introduit la contrainte v = R.
suivantes
dx Rcosd =

dy Rsind =

que lon ne peut integrer quapr`es avoir resolu le probl`eme complet: on a deux relations dierentielles du
).
type f (x,
,
Les forces de liaison nous sont le plus souvent inconnues et ne nous interessent pas: on voudrait
simplement pouvoir calculer le mouvement de notre syst`eme soumis `a des forces exterieures (appliquees)
et qui, elles, sont connues.
Par ailleurs, sil y a k contraintes, les degres de liberte reels du syst`eme se reduisent `a n = 3N k.
Cela signie que, dans la formulation newtonnienne, on resoud trop dequations (un nombre k dentre
elles se deduisent des autres).
Lidee simple est alors dexprimer les lois de la mecanique en fonction, non pas des coordonnees habituelles
de position
ri avec i=1,. . . ,N, mais des coordonnees dites generalisees independantes qj , j=1, . . . ,n. Les
coordonnees generalisees les plus naturelles correspondent aux n degres de liberte du syst`eme. Il sut a
priori didentier les coordonnees q et de faire ensuite toute la cinematique avec elles,

ri =
ri (q1 , . . . , qn , t)
D
enition: On appelle contraintes holonomes, toutes contraintes obeissant `a une relation du type
f (
r1 , . . . ,
rN , t) = 0
dierentiable en tout point. Si les contraintes sont holonomes, alors on peut exprimer une ou plusieurs
coordonnees en fonction des autres, et ceci doit etre vrai partout.
Les contraintes sont dites scleronomes si elles ne dependent pas explicitement du temps, rheonomes dans
le cas contraire.
Remarques:
(1) Dans les exemples precedents, les 1,2 et 3 sont holonomes tandis que le 4 est non-holonome. Lexemple
3 est rheonome.
(2) Un syst`eme rheonome est un syst`eme ouvert. Un syst`eme ferme (autonome) est necessairement decrit
par des contraintes scleronomes.
(3) Les probl`emes holonomes ont toujours (au moins formellement) une solution. Par contre, il nexiste
pas de methode generale pour traiter les probl`emes non-holonomes.
(4) la physique moderne est essentiellement sub-atomique et la notion de contrainte y est rare. Quand
elle apparait, cest souvent sous la forme dune modelisation holonome.

1.2. EQUATIONS DE LA DYNAMIQUE

1.2
1.2.1

Equations de la dynamique
Principe de dAlembert

En choisissant de faire la cinematique avec les coordonnees generalisees, nous sommes s

urs de travailler avec
n variables independantes. Mais il nous reste maintenant `a voir sil est possible dexprimer les lois de la
dynamique dun syst`eme `a n degres de liberte en fonction uniquement des forces exterieures F
ext . Autrement
dit, comment faire disparaitre les forces de liaison F
l ?
Principe de dAlembert:
Lors dun deplacement virtuel dun syst`eme, les forces de liaison ne travaillent pas.

i
Un deplacement virtuel correspond `a un deplacement de chaque vecteur position
ri dune quantite r


`a un instant t donne. Un deplacement reel dr, met en jeu une translation correspondante dans le temps
(ainsi quun eventuel travail des forces de liaison). Le principe de dAlembert stipule donc que les seuls
deplacements virtuels possibles sont ceux qui sont compatibles avec les forces (internes ou non) de liaison et
donc nengendrent aucun travail.
Si les contraintes ne varient pas au cours du temps (contraintes scleronomes), alors le deplacement virtuel
est equivalent `a un deplacement reel.
Le principe de dAlembert ne se verie que par lexperience. Pour se convaincre malgre tout de sa validite,
examinons quelques cas:
Cas du pendule:
Les seuls deplacements possibles de la masse seectuent selon un angle . Lors dun deplacement virtuel ,
la tension de la tige exerce un travail nul (deplacement perpendiculaire `a la force de liaison).
Cas de la boule:
Une boule roulant sans glisser sur un plan peut se deplacer selon deux directions x et y. La force de
liaison qui lempeche de glisser est dirigee selon laxe z et ne va donc pas engendrer de travail (en fait,
lapproximation sans glissement signie quon neglige toute forme de dissipation par rapport `a lenergie
cinetique de la boule).
Cas dune contrainte mobile:
Prenons le cas dune particule contrainte de se deplacer sur une courbe, elle-meme mobile. La force de

de la particule pendant
contrainte (`a t xe) est normale `a la courbe instantanee, mais le deplacement dr
lintervalle dt nest pas tangent `a la courbe. Consequence: la force de contrainte nest pas normale au
deplacement reel et produit donc un travail.
A partir de maintenant, on utilisera les indices grecs () pour caracteriser une particule parmi les N
constituant le syst`eme, et les indices latins (i, k) pour caracteriser une coordonnee generalisee parmi les n.
En vertu du principe de dAlembert, le travail virtuel des forces totales sur un syst`eme est donc simplement
W =

=
(F
ext + F
l ) r

F
ext, r

Or, le deplacement virtuel verie

=
r

r
k

qk

qk

puisque t = 0 dans un deplacement virtuel. Le travail des forces exterieures sexprime alors en fonction des
coordonnees generalisees
W

N

F
ext,

=1

k

r
k

Qk qk

qk

qk


CHAPTER 1. MECANIQUE
DE LAGRANGE

o`
u nous avons introduit la force generalisee dont la k-i`eme composante secrit
Qk =

N

r
F
ext,
qk

=1

1.2.2

(1.1)

Equations de Lagrange

Dapr`es la relation fondamentale de la dynamique, `a ce travail des forces exterieures lors dun deplacement
virtuel correspond une variation
due `a une variation dimpulsion mesuree dans un referentiel
 denergie

. Il nous reste donc `a calculer le terme de droite en fonction
galileen, cest `a dire W = m
r r
des coordonnees generalisees.
On a

m
v r

m
v

,k

r
qk =
Ak qk
qk
k

o`
u les coecients Ak (parfois appeles accelerations generalisees) sont
Ak

d
dt

m
v

r
qk

r
m
v
qk

d
m
v
dt

r
qk

Or, on peut intervertir les derivees par rappport `a des coordonnees independantes. En eet, la vitesse secrit

v =
r =

r
k

qk

qk +

r
t

ce qui nous fournit la relation utile

v

r
=
qk
qk
Par ailleurs, le deuxi`eme terme se simplie



2
r
d
r
2
r
qi +
=
dt qk
qi qk
tqk
i




r

r
=
qi +
qk
qi
t
i
=
car


r
qk

v
qk

est une fonction des qi et du temps uniquement. En regroupant tout, on obtient

Ak

=
=







v

v
m
v
m
v

qk
k






1

1
d

2
2
m v
m v
dt qk
2
qk
2

d
dt

cest `a dire
Ak =

d T
T

dt qk
qk

1.2. EQUATIONS DE LA DYNAMIQUE

o`
u
T =

1

m v2

est lenergie cinetique du syst`eme. Le principe de dAlembert permet donc de reecrire la relation fondamentale
de la dynamique sous la forme
n

(Ak Qk )qk = 0

k=1

o`
u les n deplacements virtuels qk sont quelconques et independants. Cette independance decoule directement
du fait que les contraintes sont holonomes. Pour des contraintes non-holonomes, on ne pourrait rien dire.
Ainsi, les qk etant independants, lequation ci-dessus ne peut etre satisfaite que si chaque coecient est
lui-meme nul. On obtient donc n equations algebriques independantes
T
d T

= Qk
dt qk
qk
u les forces generalisees sont obtenues, soit par lequation (1.1), soit en calculant le travail W =
 o`
Q
k k qk .
Remarques:
(1) Obtenues `a partir de la RFD, ces equations ne sont valables que dans des referentiels galileens. Dans un
referentiel non galileen R
, la RFD secrit
d

p

= F
+ f
in
dt
o`
u f
in sont les forces dinertie. Les equations de la dynamique dans R
seront alors
Ak = Qk + Qin
k

(1.2)


r

in sont les forces generalisees dinertie. Par exemple, dans un referentiel anime dune
o`
u Qin
k = f
qk
acceleration
a par rapport `a un referentiel galileen, la force generalisee sera



Qin
a
m
r

k =

qk

par rapport `a un referentiel galileen, la force generalisee sera

Si R
est en rotation uniforme
Qin
k =

d U
U
T

+
dt qk
qk
qk

L


est un potentiel generalise, L


=  m
r

v
est le moment cinetique total du syst`eme
o`
u U =



r
)2 .
(
v
est la vitesse relative vue dans R
) et T
= 12 m (

(2) La RFD soccupe des forces: il faut donc faire le bilan de lensemble des forces pour calculer le
comportement dynamique dun syst`eme. Lapproche ci-dessus est energetique: il sut de ne prendre en
compte que les forces qui travaillent.
(3) Ces equations sont algebriques et non vectorielles: cest une simplication appreciable. . .

Forces conservatives: V (q)

Si la force (totale) exterieure qui sexerce sur chaque particule du syst`eme derive dun potentiel V(q), cest
`a dire si F
= V (ou encore F,i = rV,i pour i=1,2,3) alors la force generalisee secrit


CHAPTER 1. MECANIQUE
DE LAGRANGE

Qk

N



V r,i

r
=
F
ext,
qk
r,i qk
=1
=1 i

V
r
=

r qk

N

(Attention cest une notation: on ne divise pas par un vecteur!). On note que ceci est exactement lexpression
de la derivee partielle dune fonction V (
r1 , . . . ,
rN ) par rapport `a qk ,ce qui montre que la force generalisee
sexprime directement sous la forme
V
Qk =
qk
Sans perte de generalite, puisque V(q) uniquement (et pas du temps!), on peut denir la grandeur
L = T V ou lagrangien, et les n equations du mouvement prennent la forme connue sous le nom dequations
de Lagrange
d L
L

=0
dt qk
qk
Cette equation doit evidemment redonner la relation fondamentale de la dynamique.
L
V
= q
= Qk , k-i`eme composante de la force (generalisee),
Pour une particule, cela implique que q
k
k
doit etre egale `a la derivee temporelle de limpulsion (generalisee). On denit ainsi
pk =

L
qk

(1.3)

comme etant limpulsion generalisee ou moment conjugue de qk .

Potentiels g
en
eralis
es: V (q, q)

On remarque que lon peut encore mettre les equations du mouvement sous la forme lagrangienne ci-dessus
meme si le syst`eme nest pas conservatif dans le sens usuel. Il sut que lon puisse denir un potentiel
generalise V (q, q)
tel que la force generalisee secrive
Qk =

d V
V

dt qk
qk

+
v B),

derive du potentiel generalise
Exemple: montrer que la force de Lorentz, qui secrit F
= q(E
suivant

r , t))
V = q(U (
r, t)
v A(

(1.4)

sont les potentiels scalaire et vecteur (E

= U
o`
u U et A


A
t

= A).


et B

Forces dissipatives: fonction de Rayleigh

Si toutes les forces sexercant sur un syst`eme ne derivent pas dun potentiel (meme generalise), on peut
toujours ecrire les equations de Lagrange sous la forme
d L
L

= Qk
dt qk
qk
o`
u les Qk sont les forces generalisees qui ne derivent pas dun potentiel. Un cas particulier important
concerne les forces de frottement qui secrivent sous la forme Fi = ki vi . Les forces de ce type peuvent en
eet sobtenir `a partir dune fonction, appelee fonction de dissipation de Rayleigh denie par
1

2
2
2
(kx vx
+ ky vy
+ kz vz
)
F=
2
On montre sans dicultes que la force generalisee de frottement est alors Qk = F
qk .

1.3. LOIS DE CONSERVATION

1.2.3

Exemple 1: le pendule

Soit un pendule de longueur l avec une masse m place dans un champ de pesanteur
g et astreint `a se deplacer
dans un plan (x,y). Ce syst`eme poss`ede donc 2 dimensions et 1 contrainte x2 + y 2 = l2 , donc 1 seul degre
de liberte. On choisit comme coordonnee generalisee.
u (
uz ). Lenergie cinetique vaut alors T = 1 mv 2 = 1 ml2 2 .

ur = l


=

La vitesse secrit
v = l
ur = l

2
2
On a ensuite deux methodes possibles de resolution.
Methode 1: On ne connait pas lexpression du potentiel. On calcule donc le travail lors dun deplacement

= l
u . La seule force qui travaille est le poids, on a donc
virtuel r

= mgl sin = Q
W = m
g r
ce qui fournit lequation
d T
T
= Q = mgl sin

dt

cest `a dire + 2 sin = 0, avec 2 = g/l.

Methode 2: On sait que le potentiel secrit (`a une constante pr`es)
V = mgl(1 cos )
Le lagrangien est L = T V et on ecrit directement lequation de Lagrange
d L L

=0
dt

qui redonne evidemment le meme resultat.

1.2.4

Exemple 2: masse sur une tige avec ressort

Soit une masse m astreinte `a se deplacer sur une tige indeformable, faisant un angle avec la verticale Oz,

=
uz . La masse est attachee `a un ressort de constante de
en rotation imposee avec un vecteur vitesse
raideur k et de longueur `a vide l0 et glisse sans frottement. Elle est par ailleurs soumise au poids. Ce syst`eme
est `a 1 degre de liberte, on choisit la distance r = OM comme coordonnee generalisee. Le referentiel choisi
est celui du laboratoire, donc galileen.


ur . On obtient une energie cinetique
La vitesse secrit
v = r

ur + r
T =

m 2
(r + r2 2 sin2 )
2

Les equations de Lagrange secrivent

d T
T

= Qr
dt r
r
o`
u Qr est la force generalisee totale associee `a la coordonnee r. Un deplacement virtuel, ie. compatible avec

= r
ur , ce qui nous donne un travail virtuel d
u au poids et au ressort
les forces de liaison, est de la forme r

k(r l0 )
ur r

= (mg cos k(r l0 ))r = Qr r
W = m
g r
Lequation du mouvement de la masse est donc
r = ( 2 2 sin2 )r g cos + 2 l0
en posant 2 = k/m.

1.3

Lois de conservation

Lensemble des lois qui vont suivre ne servent qu`a simplier la resolution des equations de Lagrange.


CHAPTER 1. MECANIQUE
DE LAGRANGE

1.3.1

Variables cycliques

D
enition: Une variable qi est dite cyclique si le lagrangien L ne depend pas explicitement de cette variable.
Th
eor`
eme: Si qi est cyclique, alors son moment conjugue pi est une constante du mouvement ou integrale
premi`ere.
Demonstration: Lequation de Lagrange secrit
L
d L
=
=0
dt qi
qi
Do`
u pi = L
qi est une constante du mouvement.
Linterpretation est aisee: si L ne depend pas de qi , cela signie que le syst`eme mecanique lui-meme ne
depend pas de cette variable. On voit donc apparaitre ici un lien entre les symetries dun syst`eme et ses
invariants.
Si qi = x est une longueur, alors pi = mx est limpulsion associee. Linvariance de L par rapport `a une
translation selon x se traduit donc par la conservation de la quantite de mouvement.
Si qi = est un angle, alors pi = I est le moment cinetique associe. Linvariance de L par rapport `a une
rotation dangle implique la conservation dune composante du moment cinetique.

1.3.2

Lagrangien ind
ependant du temps

Le temps joue un role particulier puisquon fait des derivations par rapport `a lui. Que se passe-t-il si le
lagrangien L = L(q, q)
ne depend pas explicitement du temps ?

Exprimons lenergie cinetique T = 12 m v2 en fonction des coordonnees generalisees:
v2

=
v
v =

=

r
t





r

r

r

qi +
qj +
qi
t
q
t
j
j




r
i

2
+2


r  
r 


r  
r 
qi +

qi qj

t
qi
qi
qj
i
i,j

Dans le cas de contraintes scleronomes t = 0, lenergie cinetique est uniquement une fonction quadratique
de qi , cest `a dire
1

T =
mij (q)qi qj
(1.5)
2 i,j
o`
u la matrice (m) est reelle, symetrique et

mij = mji =


m

r
qi

 


r

qj


r
q1

..
.


r
q1

...
..
.


r
qn

..
.


r
q1


r
q1


r
qn

...


r
qn


r
qn

Dapr`es cette expression, on voit que lenergie cinetique T est une fonction homog`ene du second degre en
qi (T (q, q)
= 2 T (q, q))
1.
1 Th
eor`eme

dEuler: si f (qi , Qi ) = n f (qi , Qi ), o`

u Qi sont toutes les autres variables, alors

Ceci se demontre en posant ui = qi , calculant

df
d

qi

f
= nf
qi

f dui
i ui d

= nn1 f puis en posant = 1.

1.3. LOIS DE CONSERVATION

La variation du lagrangien avec le temps est alors donnee par



 L
dL
L
=
qi +
qi
dt
qi
qi
i


 d  L 
L
qi
=
qi +
dt qi
qi
i





d L
d

=
qi =
pi qi
dt qi
dt
i
i
o`
u lon a utilise les equations de Lagrange dans le passage `a la deuxi`eme ligne. On obtient donc bien une
grandeur invariante, appelee integrale de Jacobi,


dH
=0
(1.6)
pi qi L avec
H=
dt
i
mais quelle est sa signication ?

T
(a) Pour V(q), on a pi = L
j mij (q)qj (le terme 1/2 disparait puisque T est une fonction
qi = qi =
homog`ene du second degre, voir note sur theor`eme dEuler). Do`
u



H=
mij (q)qj qi T + V = 2T T + V = T + V = E
i

est lenergie mecanique totale du syst`eme. On a donc conservation de lenergie dun syst`eme mecanique si
celui-ci est invariant par translation dans le temps (syst`eme autonome ou ferme: contraintes et potentiel ne
dependant pas explicitement du temps). Il faut noter quici V ne contient que le travail des forces externes
ou appliquees (absence des forces de contrainte).
(b) Lorsque V (q, q),
on a
H

L
qi L
qi
T
V
= qi
qi
T +V
qi
qi
V
= T + V qi
qi
= pi qi L =

Dans le cas de la force de Lorentz, par exemple, on obtient ainsi H = T +qU = E, somme de lenergie cinetique
et de lenergie electrique. Dans dautres cas, H est bien conserve (integrale premi`ere), mais H = T + V .
Remarques:
(1) Si les contraintes sont scleronomes mais V depend explicitement du temps (par exemple: particules
placees dans un champ exterieur variable), alors H = E et lenergie varie comme
L
V
dE
=
=
(1.7)
dt
t
t
Un tel cas correspond `a un syst`eme ouvert, recevant ou perdant de lenergie par lintermediaire du champ
impose.
t) dun syst`eme peut se mettre sous la forme L = L1 (q, q)
t),
(2) Si le lagrangien L(q, Q, q,
Q,
+ L2 (Q, Q,
L1
alors H1 = q q L1 est une integrale premi`ere.
(3) Un syst`eme mecanique ferme (=autonome) nest donc possible que pour V ne dependant pas explicitement du temps. Sil poss`ede n degres de liberte, alors il y a au plus 2n-1 integrales premi`eres independantes:
elles correspondent aux 2n conditions initiales moins une, servant `a xer le choix de lorigine des temps:
qi (t) = qi (C1 , . . . , C2n , t) = qi (C1 , . . . , C2n1 , t + t0 )
(4) Toutes les grandeurs conservatives liees aux proprietes de lespace-temps sont additives (energie E,

Cela est d
impulsion p
et moment cinetique L).
u au fait que leur denition ne depend pas de lexistence ou
non dune interaction entre les particules.


CHAPTER 1. MECANIQUE
DE LAGRANGE

10

1.3.3

Th
eor`
eme de Noether

Enonc
e: Soit un jeu de coordonnees generalisees qi (s) dependant contin
ument dun param`etre s et tel que
q , q , t) = L(q, q,
t), alors
qi (0) = qi . Si Le lagrangien L est independant de s, cest `a dire si L(


L d
qk
I(qk , qk ) =
(1.8)
qk ds s=0
k

est une constante du mouvement.

Demonstration: L independant de s secrit
dL
L d
L dqi
qi
=
+
=0
ds
qi ds
qi ds
i
Or,

dqi
ds

qi
d d
dt ds

et, avec lequation de Lagrange

dL
ds

L
qi

d L
dt qi ,

on obtient

d L d
L d d
qi
qi
+

dt qi ds
qi dt ds
i


qi
d
L d
=0
=
i ds
dt

ce qui, pour s = 0, prouve le theor`eme.

Ce theor`eme ore le lien rigoureux entre symetries et integrales premi`eres que nous avons vu precedemment.
Supposons quun syst`eme soit invariant par translation dans une direction x. On peut alors faire un changement de coordonnees tel que qi (s) = qi +s, pour i tel que qi soit associe `a la coordonnee x de chaque particule
du syst`eme et qk (s) = qk pour les autres. Alors, dapr`es le theor`eme de Noether,



L
=
p,x
qk

I=

est un invariant: cest la somme des composantes x des impulsions generalisees.

Si un syst`eme est invariant par translation dans les trois directions, alors on peut repeter ce jeu et on
obtient que limpulsion totale




P
=
p
=
m
r
(1.9)

est une integrale premi`ere. A noter que ceci est valable pour les quantites de mouvement comme pour les
moments cinetiques dun syst`eme. Invariance par rotation autour dun axe Oz implique la conservation de
la composante selon z du moment cinetique total du syst`eme. Un syst`eme possedant une symetrie spherique
a un moment cinetique total




J
=
J
=
m
r
r
(1.10)

conserve.

1.4

Une application: force centrale entre deux corps

Soit un syst`eme mecanique ferme, constitue de deux particules de masses m1 et m2 , situees respectivement
en
r1 et
r2 et interagissant par lintermediaire dun potentiel V (
r1 ,
r2 ). Dans un referentiel galileen, le
lagrangien de ce syst`eme secrit
L=T V =


2
1  2
m1
r1 + m2
r 2 V (
r1 ,
r2 )
2

1.4. UNE APPLICATION: FORCE CENTRALE ENTRE DEUX CORPS

11

Cest un syst`eme `a n = 6 degres de liberte, on a donc 6 equations dierentielles du second ordre, couplees,
et donc 12 constantes dintegration. La resolution de ce probl`eme va etre grandement simpliee en utilisant
les lois de conservation. Nayant pas de contraintes particuli`eres, on choisit comme coordonnees generalisees
qi = ri . Reste `a choisir le syst`eme de coordonnees, par exemple cartesien q1 = r1,x , . . . , q4 = r2,x , . . . , q6 =
r2,z .

1.4.1

Invariance par translation dans le temps

Le syst`eme est ferme, L ne dependant pas explicitement du temps, donc il y a conservation de lenergie, cest
`a dire


pi qi L = T + V = E
H=
i

1.4.2

Invariance par translation dans lespace

Si le potentiel dinteraction est tel que V (
r1 ,
r2 ) = V (
r2
r1 ), alors il est evident que L reste inchange par
translation dans les trois directions. En vertu du theor`eme de Noether, cela implique que limpulsion totale
du syst`eme
P
=
p1 + p
2 = m1
r 1 + m2
r 2
(1.11)
est une constante (en norme et en direction). Cela signie que le centre de masse du syst`eme a un mouvement
rectiligne uniforme. Cela nous sugg`ere donc de faire le changement de variable suivant

r


R

=
r2
r1
m1
r1 + m2
r2
=
M

o`
u M = m1 + m2 est la masse totale du syst`eme. La premi`ere coordonnee porte sur la distance entre les

= P
est une constante. Le mouvement du
deux particules et la seconde est celle du centre de masse: M R
centre de masse est donc sans interet et compl`etement determine par 6 conditions initiales.
Dans ce nouveau jeu de coordonnees, le lagrangien secrit
L=

1
2 1 2
M R +
r V (
r)
2
2

sont cycliques, donc que

o`
u = m1 m2 /M est la masse reduite. On retrouve que les 3 coordonnees de R


M R est une constante. Du coup, on peut simplier le lagrangien qui devient
L=

1 2

r V (
r)
2

(1.12)

cest `a dire celui dune particule (ctive) de masse , situee en
r, soumise `a une force dirigee vers le centre
de masse. La nouvelle energie (integrale premi`ere) est
E=

1 2

r + V (
r)
2

(1.13)

2

, correspondant `a lenergie de translation du centre de masse.
et ne comporte plus le terme 12 M R
Ces resultats sont generaux pour deux particules en interaction, pourvu que la force soit radiale et ne
depende que de la distance entre les particules.

1.4.3

Invariance par rotation dans lespace

On consid`ere maintenant le cas V = V (r) uniquement, cest `a dire le cas dune force centrale (independante
de la direction). Le probl`eme de la particule ctive devient alors un probl`eme `a symetrie spherique, aucune


CHAPTER 1. MECANIQUE
DE LAGRANGE

12

direction nest privilegiee. Cela signie que le syst`eme est invariant par rotation dans les trois directions et
implique donc que le moment cinetique de la particule ctive,
L


=
r m
r =

r

(1.14)

est une constante (en norme et en direction). Cela xe ainsi 3 autres constantes dintegration: il ne reste
plus qu`a xer lenergie E et la position (r) et vitesse (r)
initiales de la particule ctive, et les 12 constantes
dintegration necessaires auront bien ete utilisees.
Soit Oz laxe porte par :


r perpendiculaire `a 

`a tout instant nest possible que si le mouvement
seectue dans le plan xOy. Dans les coordonnees polaires, le lagrangien secrit alors
L=

1
(r 2 + r2 2 ) V (r)
2

On pourrait ecrire les equations de Lagrange et tenter ensuite de resoudre le probl`eme (une fois V (r) specie).
Mais il est plus utile de partir des integrales premi`eres, sachant que le moment conjugue secrit pi = L
qi ,


E

1.4.4

= p = r2
=

1
(r 2 + r2 2 ) + V (r)
2

(1.15)

Loi horaire

Si lon sinteresse `a la loi horaire r(t) alors on peut reecrire lenergie sous la forme
E=

1 2
r + Vef f (r)
2

o`
u

2
Vef f (r) = V (r) +

2r2

Le probl`eme est ainsi ramene au calcul du mouvement dune particule dans un potentiel eectif `a 1 dimension.
On obtient alors
 r
dr

(1.16)
t t0 =
2
r0
(E

V
(r))
ef
f

relation qui fournit t = t(r) et qui, du moins en principe, permet dobtenir r(t) par inversion.

1.4.5

Trajectoire

Si lon sinteresse `a la trajectoire r(), alors il est judicieux de faire le changement de variable u() = 1/r().
A partir de lexpression de lenergie (1.15) on obtient
dV
u + u = 2
du

(1.17)

o`
u u = d2 u/d2 . Cette equation porte le nom dequation de Binet et est valable pour tout potentiel
central. Cependant, toutes les formes de potentiels V (r) ne donnent pas lieu `a des equations integrables
analytiquement. Les cas simples (et les plus etudies) sont des potentiels de la forme
V (r) = rn
Pour n = 2 (oscillateur harmonique), n = 1 (Kepler), n = 2, les solutions sont analytiques. Pour dautres
valeurs de n, on tombe sur des integrales elliptiques.

1.5. PETITES OSCILLATIONS

13

Force en 1/r 2 , Loi de Kepler

1.4.6

Dans le cas de la loi de Kepler, le potentiel secrit V (u) = Ku avec K = Gm1 m2 . Lintegration de
lequation de Binet est alors immediate et fournit
r=

p
1 + e cos( 0 )

(1.18)

2
lequation dune conique, lexcentricite e et 0 etant deux constantes dintegration et p = K le param`etre
de la conique. Pour e = 0 on obtient des cercles, 0 < e < 1 des ellipses, e = 1 des paraboles et e > 1 des
hyperboles. Il est ensuite aise de relier e `a la valeur de lenergie,

2
2E
e = 1+
K 2

1.5

Petites oscillations

Le formalisme de Lagrange se prete particuli`erement bien au traitement des petites oscillations, cest `a dire
de faible amplitude, au voisinage dune position dequilibre. Dans cette section nous allons obtenir des
resultats tr`es generaux.

1.5.1

Syst`
emes `
a 1 degr
e de libert
e

Soit un syst`eme mecanique decrit par une coordonnee generalisee q (pour simplier, elle sera consideree du
type distance), soumis `a un potentiel V(q). Son energie cinetique est T = 12 mq2 et lequation du mouvement
est fournie par lequation de Lagrange
d L
dt q

m
q =

L
q
dV

dq

Il existe une position dequilibre qe si, par denition, le potentiel y est extremal cest `a dire

dV
=0
dq q=qe
Mais cet equilibre est-il stable? Autrement dit, si on donne au syst`eme une vitesse initiale q0 lorsquil
est place en qe , va-t-il seloigner ou revenir vers qe ? Pour repondre, on fait un developpement de Taylor du
potentiel `a lordre 2 (approximation harmonique), en posant x = q qe ,

d2 V x2
+ O(3)
V (x) = V (0) +
dx2 0 2
Lequation de Lagrange linearisee est alors

dV
d2 V
x + O(2)
m
x=
=
dx
dx2 0

(a) Si

d2 V
dx2 0

> 0 lequilibre est stable. En eet, on pose

2 =

1 d2 V
m dx2 0


CHAPTER 1. MECANIQUE
DE LAGRANGE

14

et le mouvement est celui dun oscillateur harmonique, de pulsation bien denie , x(t) = acos(t + ) =
it
Re(beit + ceit
) = Re(Ae ) pour a,b,c reels ou A complexe.
d2 V
(b) Si dx2 < 0 lequilibre est instable. On pose
0

r2 =

1 d2 V
m dx2 0

et la solution secarte exponentiellement de la position dequilibre, avec un temps caracteristique 1/r, x(t) =
aert + bert .
2

(c) Si ddxV2 = 0 lequilibre est indierent. On a x
= 0 et donc tout depend des conditions initiales,
0
x = x0 + v0 t.

1.5.2

Syst`
emes `
a n degr
es de libert
e

La generalisation du traitement precedent `a des syst`emes `a n degres de liberte est simple dans le principe,
mais met en jeu des techniques de calcul matriciel. On rappelle que lon ne sinteresse qu`a des syst`emes
possedant une ou plusieurs positions dequilibre, ie.

V
=0
i = 1, . . . , n
qi qi =qie
Posons que V passe par un extremum en qi = qie et introduisons les petits deplacements xi = qi qie .
Un developpement de Taylor `a lordre 2 au voisinage de lequilibre donne
V = V (
0) +

1

Vij xi xj + O(3)
2 i,j

o`
u la matrice (V), delements Vij =

2V
xi xj 0

T =

est reelle et symetrique. Lenergie cinetique secrit

1

mij (
q)qi qj
2 i,j

A lordre 2, au voisinage de la position dequilibre (variables xi ), lenergie devient

T =

1

Mij x i x j
2 i,j

o`
u la matrice (M) delements Mij = mij (
0) est egalement symetrique et reelle. On a

L
j
j Mij x
x i =

L
j Vij xj
xi =
et les equations de Lagrange linearisees deviennent alors un syst`eme de n equations


(Mij x
j + Vij xj ) = 0
j

La methode generale de resolution consiste `a rechercher des solutions de la forme xj = Xj eit , ce qui
fournit un syst`eme de n equations algebriques



Vij 2 Mij Xj = 0
(1.19)
j

1.5. PETITES OSCILLATIONS

15

La resolution de ce syst`eme est fournie par lequation caracteristique, cest `a dire

det (V M ) = 0

(1.20)

o`
u lon a pose = 2 . Cest une equation de degre n par rapport `a . En denitive, la determination des
frequences propres se ram`ene `a une operation qui ressemble au calcul des valeurs propres de matrices. Les
matrices (V) et (M) etant symetriques et reelles, elles poss`edent n racines reelles (eventuellement multiples).
Dans la base des vecteurs propres associes aux valeurs propres, les matrices (M) et (V) sont toutes deux
diagonales. Cette diagonalisation correspond geometriquement `a un changement de coordonnees lineaire
(matrice (A)), faisant passer dun syst`eme daxes `a un autre. Mais `a ce stade, il nest pas du tout evident
quil existe une matrice (A) diagonalisant `a la fois (V) et (M)!
Quoi quil en soit, on peut toujours resoudre lequation caracteristique (1.20) et obtenir ainsi les n valeurs
propres a . On reporte ensuite chaque valeur propre a dans le syst`eme dequations



Vij a2 Mij Xa,j = 0
(1.21)
j

a correspondants. Il faut repeter cette operation

an dobtenir les n composantes Xa,j du vecteur propre X
pour obtenir lensemble des n vecteurs propres.
u Ca est un nombre complexe
La solution particuli`ere recherchee est donc de la forme xj = Ca Xa,j eia t o`

a eia t) ). La solution particuli`ere compl`ete est donc la partie reelle de la somme des
(ou encore
x = Ca X
solutions precedentes, `a savoir






Ca Xa,j eia t =
Aaj za
xj = Re
a

o`
u



za = Re Ca eia t

(1.22)

constitue une solution oscillante triviale, avec la pulsation propre a puisquelle obeit `a lequation harmonique
za + a2 za = 0

(1.23)

Quest ce que cela signie? La variation de chacune des coordonnees xj avec le temps apparait donc
comme une superposition lineaire des n etats propres independants za . On appelle ces etats les modes
propres, associes aux axes principaux ou normaux du syst`eme. Le changement de coordonnees seectue par
la matrice de passage (A), dont les elements sont Aaj = Xa,j : cest tout simplement la matrice constituee

a
des n vecteurs propres X

Xn,1
X1,1

..
.
Aaj = ..

.
.
.
.

X1,n
Xn,n
Les coordonnees normales za etant independantes, cela signie que le lagrangien du syst`eme, ecrit avec
ces coordonnees, doit pouvoir secrire comme la somme de lagrangiens independants. Chaque mode evolue
sans interagir avec les autres, cest `a dire

ma 

za2 a2 za2
(1.24)
L=
2
a
Autrement dit, la matrice de passage (A) diagonalise bien simultanement (M) et (V).
En resume, la determination des frequences propres se ram`ene toujours au calcul des valeurs propres
dune matrice. Celle-ci etant symetrique, les valeurs propres sont toujours reelles, positives ou negatives.
Les valeurs propres positives correspondent `a des modes oscillants (equilibre stable). Les modes propres
sont des modes collectifs doscillation `a une seule frequence, pouvant etre excites independamment les uns


CHAPTER 1. MECANIQUE
DE LAGRANGE

16

des autres. Dans le cas o`

u toutes les valeurs propres sont positives, on peut considerer le mouvement complet
dun syst`eme comme etant obtenu en excitant les divers oscillateurs harmoniques avec des amplitudes et des
phases dierentes.
Les valeurs propres negatives correspondent `a des solutions secartant exponentiellement de la position
dequilibre (ex: selle). Un syst`eme isole devant conserver son energie, cela signie que notre traitement
echoue pour les grandes amplitudes.
Cette remarque est egalement valable dans le cas doscillations de trop grande amplitude. Dans ce cas,
on excite egalement les harmoniques des frequences fondamentales (traitement des grandes amplitudes par
les series ou transformees de Fourier).
Enn, les valeurs propres nulles forment ce quon appelle des modes mous (en elasticite, ils sont appeles
modes rigides). Ils correspondent `a za = 0 et donc `a une vitesse de translation constante dans la direction de
la coordonnee principale associee. Puisque cette vitesse est constante, cela signie que limpulsion associee
est un invariant. A linverse, un syst`eme invariant par translation dans une direction donnee aura donc un
mode propre mou correspondant. Meme chose pour la rotation.

1.5.3

Oscillations forc
ees

Le traitement des oscillations forcees est particuli`erement simple dans les coordonnees normales. On reprend
tout depuis le debut avec une force generalisee exterieure, ne dependant pas necessairement dun potentiel2 .
Soit Qi la force generalisee exprimee avec les coordonnees xi = qi qie . Les equations de Lagrange linearisees
secrivent alors


(Mij x
j + Vij xj ) = Qi
j

et peuvent se mettre sous la forme matricielle suivante

+VX =Q
MX
o`
u X est le vecteur des coordonnees xi . Puisque A est la matrice de passage des axes principaux vers les
axes usuels, on a X = AZ. On pose alors Q = AF o`
u les Fi sont les coordonnees de la force generalisee
A1 et V = AV A1 les matrices exprimees
exprimees dans la base principale zi . Par ailleurs, soit M = AM
)ij = mi ij et (V )ij = mi 2 ij .
dans la base principale, donc diagonales: on peut donc poser (M
i
Les equations de Lagrange se ram`enent alors
+VX
MX
+ AV A1 X
A X
AM
A1 AZ + AV A1 AZ
AM
Z + V Z)
A(M
1

= Q
= AF
= AF
= AF

cest `a dire `a la forme simple

zi + i2 zi =

Fi
mi

(1.25)

En general, on sinteresse au cas particulier o`

u la force est elle-meme periodique, de pulsation , Fi =
Fi0 cos(t+ i ). La solution generale de lequation de Lagrange est alors la somme dune solution particuli`ere
et de la solution generale de lequation homog`ene. Cette derni`ere est une oscillation (les cas instables nont
pas dinteret car leur comportement est evident) `a la frequence popre i . Une telle solution est souvent un
regime transitoire. Dans ce qui suit, on ne sinteresse qu`a la solution particuli`ere.
Le moteur exterieur imposant son mouvement, on recherche des solutions particuli`eres de la forme zi =
Zi0 cos(t + i ). En les reinjectant dans lequation du mouvement on obtient les amplitudes
Zi0 =
2 On

Fi0 /mi
i2 2

ext
peut cependant ecrire le lagrangien L = T V
avec V
= V + Vext et Qi = V
.
x
i

1.5. PETITES OSCILLATIONS

17

ce qui fournit un mouvement complet (sans la partie oscillations propres) decrit par
xi (t) =

j

Aij zj =

Aij Fj (t)/mj
j

j2 2

(1.26)

Consequences:
(1) Si la force na pas de composante dans la direction de vibration dun mode normal particulier, alors
Fi0 = 0: une force exterieure ne peut exciter un mode propre particulier que si elle agit dans le meme sens
de vibration.
(2) Il y a resonance lorsque = i : lamplitude du mode propre correspondant se met `a croitre. Dans
une situation reelle, il y a toujours de la dissipation qui limite la croissance de lamplitude et la maintient
nie.

18

CHAPTER 1. MECANIQUE
DE LAGRANGE

Chapter 2

Principe variationnel
2.1

Le principe de Hamilton

Jusqu`a present, nous avons decrit letat dun syst`eme mecanique `a n degres de liberte par la donnee, `a
chaque instant t, des n coordonnees generalisees qi (t). A un instant donne, on peut donc representer letat1
de ce syst`eme par un point dans un espace cartesien de dimension n, appele espace des congurations.
A chaque axe de cet espace correspond une coordonnee qi . Un syst`eme mecanique subissant une evolution
entre deux instants t1 et t2 va donc dessiner une courbe dans lespace des congurations entre un point 1 et
un point 2, quon appelle (faute de mieux) trajectoire. Le temps peut alors etre pris comme param`etre de
cette courbe.
La trajectoire reelle est celle qui correspond eectivement `a la dynamique suivie par le syst`eme, elle sera
donc obtenue en resolvant les equations de Lagrange. Mais quavons-nous fait pour les obtenir? En suivant
le principe de dAlembert, nous avons considere le travail lors de deplacements virtuels. Graphiquement,
cela revient `a considerer des chemins2 dierents mais tr`es proches, reliant les points 1 `a 2. Cest cette constatation qui nous am`ene au principe de Hamilton.
Enonc
e: Le mouvement dun syst`eme reel, depuis linstant t1 jusqu`a linstant t2 , est tel que lintegrale


t2

S=

L dt

(2.1)

t1

o`
u L = T V est le lagrangien du syt`eme, est extremale. Par extremale, on entend quelle reste stationnaire,
cest `a dire S = 0, lors dune variation fonctionnelle des chemins. On appelle cette integrale laction du
syst`eme.
Remarques:
(1) Dun point de vue mathematique, calculer laction dun syst`eme revient `a faire une application qui, `a
une fonction

q (t) (la trajectoire), associe un nombre. Ce nombre depend donc de la fonction utilisee: on dit
que S est une fonctionnelle de la trajectoire.
(2) Laction doit etre extremale, cest `a dire posseder un minimum ou un maximum. Un grand nombre
de conditions physiques sont en fait decrites par un minimum, ce qui fait que le principe de Hamilton est
parfois appele principe de moindre action.
(3) Nous allons demontrer que ce principe permet de retrouver les equations de Lagrange. Or, celles-ci
sont elles-memes equivalentes aux equations de Newton. On peut donc remplacer les trois principes de
Newton par celui-ci!
1 L
etat

du systme ne consiste en sa position que si les variables qi sont des variables despace.
va maintenir cette distinction: le terme de trajectoire 
q (t) decrira la solution des equations de la dynamique, tandis que
celui de chemin decrira une fonction 
q (t) virtuelle.
2 On

19

20

CHAPTER 2. PRINCIPE VARIATIONNEL

(4) Cette formulation de la mecanique poss`ede en outre lavantage detre independante du syst`eme de
coordonnees choisi pour exprimer L.

2.2

D
eduction des
equations de Lagrange

Nous allons prouver que le principe de Hamilton permet eectivement de retrouver les equations de Lagrange.
Autrement dit, quil contient les equations de Lagrange. Soit laction


f (y1 , . . . , yn , y1
, . . . , yn
, x)dx

S=
1

Soit
y(x) la trajectoire recherchee, cest `a dire celle qui assure S(y) = 0. On peut construire des chemins
voisins de cette solution en introduisant un param`etre
et une fonction
(x) quelconque telle que

y (x,
) =
y (x) +


(x)
o`
u
(x1 ) =

(x2 ) =
0: tous les chemins passent par les points 1 et 2.
La variation de laction (qui est devenue une fonction de
) est alors
S =

dS

=
d

x2

x1



 f yi
f y

+
i
dx
yi

yi

On int`egre par parties le deuxi`eme terme,



x2

x1

f yi

dx
yi

do`
u

x2

f 2 yi
dx


x1 yi
x


x

 x2
f yi 2
f yi
d
=
dx

yi

x1
yi

x1 dx

x2

S =
x1

yi
i

f
d

yi dx

f
yi


dx

i
o`
u lon peut identier y
etre verie quelque soit le

= yi . Le principe de Hamilton (S = 0) ne peut
chemin (les yi sont independants), que si, pout tout indice i

d
dx

f
yi

f
=0
yi

(2.2)

ce qui redonne les equations de Lagrange lorsquon eectue le changement de variable x t.

Th
eor`
eme: La fonction de Lagrange L dun syst`eme mecanique nest determinee qu`a une fonction pr`es,
secrivant comme la derivee totale par rapport au temps dune fonction quelconque f (q, t) du temps et des
coordonnees.

d
t) = L(q, q,
t) + dt
f (q, t). Laction associee `a ce lagrangien est alors S
= L
dt = S +
Soit L
(q, q,
f (q2 , t2 ) f (q1 , t1 ). Do`
u S
= S = 0.
Consequence: ce theor`eme tr`es utile permet de simplier un lagrangien de tous les termes qui peuvent
d
se mettre sous la forme dt
f (q, t). On peut se ramener ainsi `a resoudre des equations considerablement plus
simples.

2.3. EXEMPLES SIMPLES DE CALCUL VARIATIONNEL

2.3

21

Exemples simples de calcul variationnel

On vient incidemment de demontrer que toute une classe de probl`emes se ramenant au calcul de variations
(cest `a dire lorsquon cherche une fonction de une ou plusieurs variables rendant extremale une quantite)
peut se resoudre par les equations de Lagrange. Ceci nest valable, evidemment, que si la fonctionnelle `a
extremiser secrit bien sous la forme
 x2
f (y, y
, x)dx

S=
x1

cest `a dire ne met en jeu ni des derivees secondes y ou dordre plus eleve, ni des variables supplementaires
dintegration, par exemple,
 
f (y, y
, x1 , x2 )dx1 dx2

S=

Dans ces cas l`a, la methode utilisee ci-dessus fonctionne, mais donne des equations qui ne sont pas celles de
la mecanique de Lagrange.

2.3.1

Plus petite distance dans un plan

Quelle est la courbe y(x) qui minimise la distance

 entre deux points A et B dans un plan ?
Lelement innitesimal de distance est ds = dx2 + dy 2 . La grandeur `a minimiser est alors


xB

ds =

S=
A


1 + y
2 dx

xA

Ici, la variable x joue le role du temps et le lagrangien est


L(y
) = 1 + y
2
La variable y est cyclique donc le moment conjugue py est constant, cest `a dire
py =

L
y

=


y
1 + y
2

On obtient donc que y
= dy/dx est une constante, cest `a dire lequation dune droite.

2.3.2

La brachistochrone

Sous laction de la pesanteur seule, un point materiel de masse m glisse sans frottement dans un plan vertical.
Quelle est lequation de la courbe joignant 2 points, O et A, dans le temps le plus court ?
La grandeur `a minimiser est evidemment
 xA 
 A
 A
ds
1 + y
2

=
S=
dt =
dx
2gy
O
O v
0
o`
u la variable x joue encore le role du temps et le lagrangien est

1 + y
2


L(y, y ) =
2gy
Ici, x est cyclique, ce qui implique la conservation de la grandeur
H = y

1
L
L= 


y
2gy(1 + y
2 )

u a est une constante.

ce qui fournit y(1 + y
2 ) = a o`

22

CHAPTER 2. PRINCIPE VARIATIONNEL

On peut resoudre cette equation en donnant une forme parametrique y(), x(). Pour cela, on pose
u = y
, ce qui fournit

ay
u=
y
Essayons la fonction dessai y = a sin2 , cest `a dire u =
cos 2)d et on obtient apr`es integration
x

cos
sin .

Alors, dx =

dy
u

2a sin cos d
u

= a(1

a
(2 sin 2)
2
a
(1 cos 2)
2

qui est lequation dune cycloide.

2.4
2.4.1

G
en
eralisation des
equations de Lagrange
Forces non conservatives

Nous avons vu que le principe de Hamilton permet de deduire les equations de Lagrange pour des syst`emes
soumis des forces conservatives, cest `a dire pour lesquels toutes les forces en presence derivent dun potentiel
(eventuellement generalise). Que se passe-t-il pour une force F
non conservative?
Si on reprend la demonstration des equations de Lagrange, on voit que celles-ci secrivent


T
d T
= Qi

dt qi
qi

r
o`
u T est lenergie cinetique et Qi = F
q
la force generalisee.
i
On peut egalement obtenir ces equations `a partir dun principe variationnel en ecrivant que laction est
cette fois-ci denie par

t2

(T + W ) dt

S=

(2.3)

t1

o`
u W = F

r.
En eet, le principe variationnel S = 0 implique
 t2

T dt =
t1

W dt

t1

Or, le terme de gauche, on la vu, secrit

 t2

T dt =
t1

t2

t2

t1


qi

d
dt

T
qi


T

dt
qi

tandis que celui de droite vaut



t2

W dt

t2

t1

t1
 t2

=
t1

t2

=
t1

dt
F
r


i

r
qi F

qi
qi Qi

ce qui redonne eectivement les equations de Lagrange dissipatives.

 ERALISATION

2.4. GEN
DES EQUATIONS
DE LAGRANGE

2.4.2

23

Contraintes non holonomes: multiplicateurs de Lagrange

Nous avons vu que la condition de contraintes holonomes (il existe k relations du type fk (
r1 , . . . ,
rN , t) = 0)
etait essentielle pour letablissement des equations de Lagrange, puisque cest elle qui assure lindependance
des coordonnees generalisees qi . Pour les syst`emes non holonomes, les qi ne sont pas independantes les unes
des autres. Comment faire?
M
ethode des multiplicateurs
Si les m equations de contrainte peuvent se mettre sous la forme dierentielle suivante,


alk dqk + alt dt = 0

(2.4)

(l = 1, . . . , m) alors on peut utiliser la methode dite des multiplicateurs de Lagrange.

Lors dun deplacement virtuel dt = 0 et les equations qui doivent etre satisfaites sont seulement


alk qk = 0 .

(2.5)

On introduit alors m constantes indeterminees l (qui peuvent etre des fonctions du temps) et les m relations
suivantes sont evidemment veriees


alk qk = 0
l
k

ainsi que leur version integrale


1

l

alk qk dt = 0

Pour un syst`eme lagrangien, le principe de Hamilton fournit



 2

L
d L
dt

S =
qk = 0
qk
dt qk
1
k

do`
u, la somme de ces deux variations doit egalement etre nulle, `a savoir





 2

n
d L
L
dt

l alk qk = 0
+
qk
dt qk
1
k=1

(2.6)

Dans cette equation, les n qk sont dependants les uns des autres, relies par les m contraintes (2.5). On peut
toujours choisir les n-m premi`eres coordonnees comme etant les coordonn
ees independantes. Les m derni`eres

coordonnees generalisees seront ensuite xees par les m relations k alk qk = 0. Mais les m constantes l
que nous avons introduites sont libres. On peut donc les choisir de telle sorte que



m
d L
L

l alk = 0
(2.7)
+
qk
dt qk
l=1

pour k = n m + 1, . . . , n. Si on introduit cela dans lequation (2.6) on obtient






 2 nm

L
d L
dt

l alk qk = 0
+
qk
dt qk
1
k=1

o`
u les qk mis en jeu sont cette fois-ci tous independants. Du coup, on obtient les n-m equations suivantes
pour k = 1, . . . , n m,



m
L
d L

l alk = 0
(2.8)
+
qk
dt qk
l=1

24

CHAPTER 2. PRINCIPE VARIATIONNEL

Les equations (2.7) et (2.8) se condensent sous la forme des n equations
d
dt

L
qk

L
=
l alk
qk
m

k = 1, . . . , n

(2.9)

l=1

o`
u il y a n+m inconnues: les n coordonnees qk et les m constantes l . Pour resoudre le syst`eme complet il
faut rajouter les m equations (dierentielles) de contraintes
n

alk qk + alt = 0

(2.10)

k=1

Quelle est la signication physique des multiplicateurs de Lagrange l ? Les equations de Lagrange en
presence de forces conservatives et non conservatives secrivent


L
d L
= Qk

dt qk
qk
o`
u les Qk sont les forces generalisees associees aux forces non conservatives. Les l decrivent donc les
forces (inconnues) generalisees de contrainte. Celles-ci sont donc obtenues lors de la resolution compl`ete du
probl`eme (etendu `a n+m variables).
Remarque: Cette methode des multiplicateurs de Lagrange est egalement applicable pour des syst`emes
holonomes. En eet, toute relation du type f (q1 , . . . , qn , t) = 0 peut secrire apr`es dierentiation

f
f
dt = 0
dqk +
qk
t
k

cest `a dire alk =

f
qk

et alt =

f
t .

On peut donc utiliser cette methode lorsque

il nest pas commode de ramener tous les qk `a des coordonnees independantes.

lon souhaite obtenir les forces de liaison internes `a un syst`eme.
Exemple
Soit un cerceau de rayon R et de masse M roulant sans glisser sur un plan incline dangle , sous leet de
son poids. Lenergie cintique du solide est
T =

1
1
M x 2 + I 2
2
2



o`
u I = r2 dm = R2 dm = M R2 est le moment dinertie du cerceau par rapport `a laxe passant par son
centre dinertie G, x est la vitesse de G et la vitesse angulaire de rotation du cerceau sur lui-meme. Ce
probl`eme poss`ede donc 2 coordonnees x et , reliees entre elles par la contrainte de roulement sans glissement
x = R
Cette contrainte peut se mettre sous la forme holonome x = R + C, permettant ainsi de traiter le probl`eme
par les equations de Lagrange usuelles. Mais ce faisant, nous ne serons pas capables de calculer la force de
contrainte qui permet justement au cerceau de ne pas glisser. Si on veut la calculer, il faut introduire un
multiplicateur de Lagrange et conserver les deux variables comme coordonnees generalisees. Lequation de
contrainte fournit
dx Rd = 0
Le potentiel secrit
V = M g sin x

2.5. EXPRESSIONS DU LAGRANGIEN EN FONCTION DE LESPACE-TEMPS

et les equations de Lagrange generalisees deviennent


L
d L

=
dt x
x


d L
L
=

dt

25

a1x =
a1 = R

On obtient ainsi un syst`eme de 3 equations `a 3 inconnues (x, , )

Mx
M g sin =
M R =
x =

cest `a dire
g sin
2
x
=
R
M g sin
=
2
plus petite que celle quil aurait en labsence des frottements

x
=

Le cerceau descend avec une acceleration 2 fois

dus `a la contrainte , dirigee selon x.

2.5
2.5.1

Expressions du lagrangien en fonction de lespace-temps

M
ecanique non relativiste

Nous allons suivre une demarche deductive (`a la Landau) qui, `a partir de principes premiers tr`es simples,
va nous permettre de determiner lensemble des principes de la mecanique de Newton.
Tout dabord, il est necessaire de choisir un syst`eme de reference (un observateur) pour y exprimer les
lois de la physique, puisque tout evenement physique est relatif `a un observateur. Cependant, les lois memes
de la physique doivent etre independantes de ce choix. On peut donc choisir celui o`
u les lois y adoptent la
forme la plus simple. Un referentiel galileen est ainsi un syst`eme de reference privilegie dote des proprietes
suivantes:
(1) lespace est homog`ene et isotrope;
(2) le temps y est uniforme (= le meme partout). On dit egalement que le temps est absolu.
Lensemble de ces premisses constitue le principe de relativite de Galilee.
En vertu du principe de Hamilton, le mouvement dune particule materielle se deplacant librement dans
lespace, est tel que son action
 t2
L dt
S=
t1

est extremale. La fonction de Lagrange L ne peut etre fonction

- ni de
r: espace homog`ene (inv par translation, donc L aussi)
- ni de t: temps uniforme (inv par translation). On doit avoir L(
r
, t
) = L(
r, t). Cela signie que L ne peut
etre quune fonction de la vitesse
v =
r . Mais lespace etant isotrope (toutes ses directions sont equivalentes),
L ne peut dependre de la direction de
v . Elle depend donc de sa valeur absolue, cest `a dire
L = L(v 2 )
Les equations de Lagrange (qui decoulent du principe de Hamilton) fournissent ensuite


d L
L
=0
=
dt qi
qi

(2.11)

26

CHAPTER 2. PRINCIPE VARIATIONNEL

v est une constante (car L

donc L
qi est une constante, ce qui implique que la vitesse

qi est une fonction de
qi et v uniquement). Ainsi, dans un referentiel galileen, le mouvement dune particule libre seectue avec
une vitesse uniforme. Nous venons de demontrer le principe dinertie, postule dans le cadre de la mecanique
newtonienne.
Soit L = L(v 2 ) dans un referentiel galileen R. Reste `a trouver la dependance fonctionnelle de L en v 2 .

par rapport `a
Soit
v
la vitesse de la particule dans un autre referentiel galileen R
, anime dune vitesse V
R. Les formules de changement de referentiel sont donnees par la transformation de Galilee,

r

t
=
r
+ V

= t

(2.12)

Si V
est faible par rapport `a
v
, on a

L(v 2 ) =


+ V 2)
L(v
+ 2v

V
L
2
L(v
) +
2 2v

V

v
2

d
Or, puisque L(q, q,
t) = L
(q
, q
, t) + dt
f (q
, t), cela signie que le deuxi`eme terme du developpement limite
doit secrire comme une derivee temporelle totale dune fonction de
r
et du temps, cest `a dire

d
L dr



2
V f (r

, t)
dt
v
2 dt
L
ependant de
r
et de t, cest `a dire une constante a. Donc L = av
, o`
u a doit
Cela implique que v

2 est ind
etre une caracteristique intrins`eque de la particule. Ce qui est vrai dans R
lest egalement dans R (donc

quelconque. En eet, on a
L = av 2 ) et pour une vitesse V
2

L

=
=

)2 = av 2 2a
v V
+ aV 2
av
2 = a(
v V

d 
2a
r V
aV 2 t
L
dt

Puisque L = av 2 pour une particule libre dans tout referentiel galileen, les equations de lagrange secrivent


d L
L
=0
=
dt vi
ri
d
(2avi ) = 2av i = 0
dt
d
v
=
0
2a
dt
Pour que les equations de Lagrange dune particule libre soient compatibles (et, de ce fait, demontrent)
la relation fondamentale de la dynamique de Newton, il sut de poser
a=

m
2

2
Le lagrangien devient alors L = m
2 v = T pour une particule de masse m.
Pour un syst`eme de particules ninteragissant pas entre elles, les equations du mouvement de lune ne
peut contenir des grandeurs se reportant aux autres: cela implique que le lagrangien du syst`eme est une
somme de lagrangiens independants et

m


v2
L =
L=
2

Lorsque les particules interagissent entre elles (syst`eme ferme), il sut de denir le lagrangien comme
etant
L=T V

2.5. EXPRESSIONS DU LAGRANGIEN EN FONCTION DE LESPACE-TEMPS

27

o`
u le potentiel V ne depend que des positions
r . Les equations de Lagrange secrivent alors
m

d
v
= F

dt

V
o`
u F
= 
r (notation vectorielle) est la force qui sexerce sur la particule et qui est la RFD.
Lhomogeneite de lespace implique quun syst`eme ferme reste invariant par translation densemble. Cela
implique que le lagrangien egalement reste invariant lorsquon applique une translation
r =

identique `a
tous les vecteurs position, cest `a dire

L =
On obtient alors

L

L

r =


=0

r

r

 V 


L
F
=
0
=

=


r


r

Dans le cas particulier dun syst`eme ferme de deux particules, cela fournit F
1 + F
2 =
0, qui est le principe
daction et de reaction.
On voit donc que le principe de Hamilton, associe `a des considerations sur les proprietes de lespacetemps (principe de relativite de Galilee) permettent de redemontrer les trois principes fondamentaux de la
mecanique de Newton:
(1) Principe dinertie;
(2) Principe de la dynamique (RFD);
(3) Principe daction et de reaction.
Par ailleurs, il est important de noter que L = T V est une consequence directe de notre a priori galileen
sur lespace-temps. Dans un referentiel non-galileen, on ne pourra utiliser L = T V , `a moins dintroduire
`a la main le potentiel qui serait associe `a la force dinertie correspondante3. En pratique, il vaut mieux
ecrire le lagrangien dans un referentiel galileen, puis de faire un simple changement de variables.

2.5.2

M
ecanique relativiste

On peut faire de meme dans le cadre de la mecanique relativiste, en remplacant le principe de relativite de
Galilee par celui de Lorentz. Nous allons obtenir de cette facon le lagrangien relativiste en suivant plusieurs
etapes:
(1) Laction dune particule libre doit etre en eet denie de facon independante du referentiel: elle doit
donc etre invariante par transformation de Lorentz. La facon la plus simple est quelle soit lintegrale dun
scalaire, lui-meme etant un invariant de Lorentz.
(2) Cette fonction scalaire doit cependant mettre en jeu des dierentielles du premier ordre: en eet, les
equations de Lagrange qui resultent du principe de Hamilton (et qui mettent en jeu des derivees de cette
fonction), produisent des equations avec des derivees secondes des positions.
Le seul scalaire construit `a partir de
u a > 0 est une constante
dierentielles du premier ordre est ads o`
caracteristique de la particule et ds = c2 dt2 dl2 lintervalle despace-temps de Minkowski.
(3) On pose que laction dune particule libre relativiste secrit alors

S = a

ds

(2.13)

Le signe moins se justiera plus tard. Il indique quici S doit etre minimale (principe de moindre action)
pour un mouvement reel.
3 Dans

un referentiel en rotation, lespace nest plus isotrope puisquil existe une direction privilegiee: laxe de rotation.

28

CHAPTER 2. PRINCIPE VARIATIONNEL

(4) On peut ensuite identier le lagrangien L dune particule libre comme etant

v2
L = ac 1 2
c

(2.14)

(5) Reste `a determiner la constante a. La limite newtonienne est obtenue en faisant tendre la vitesse de
la lumi`ere c vers linni, cest `a dire


v2
a
L  ac 1 2 = ac + v 2
2c
2c
Le premier terme est une constante sans inuence tandis que le second redonne bien lenergie cinetique de
la particule si on pose
a = mc
(incidemment, on voit lutilite du signe moins: la mecanique relativiste obeit `a un principe de moindre
action).
R
esum
e: Laction et la fonction de Lagrange dune particule materielle libre de masse m sont
 2
ds
(2.15)
S = mc
1

v2
L = mc2 1 2
(2.16)
c
Limpulsion generalisee dune particule libre est pi =

L
vi ,

cest `a dire

m
v
p
= 
2
1 vc2

(2.17)

qui tend vers m
v lorsque

 c .
Lenergie est E = i pi vi L, ce qui donne
mc2
E= 
2
1 vc2

(2.18)

Lorsque c , lenergie devient E  mc2 + T est (`a une constante pr`es) coherente avec lexpression
newtonnienne. Cette constante est lenergie que poss`ede la particule meme en labsence de vitesse: cest
lenergie de masse au repos. On remarque donc quici, la valeur de lenergie dune particule est parfaitement
determinee. Dans la mecanique de Newton, elle nest connue qu`a une constante pr`es. Par ailleurs on voit
quaucune particule materielle ne peut aller `a v = c (divergence). Lenergie cinetique de la particule secrit
T = ( 1)mc2
o`
u= 

1
2

1 vc2

est le facteur de Lorentz.

Enn, en combinant les deux relations ci-dessus, on peut un lien direct entre energie et impulsion,

p =

E

v
c2

(2.19)

Cette equation est interessante, car elle est independante de la masse de la particule. Si p
= p
u, alors une
particule de masse nulle, allant `a la vitesse de la lumi`ere et possedant une energie E verie p
= Ec
u. Enn,
utilisant cette relation, on peut egalement construire un autre invariant
E2
= p 2 + m 2 c2
c2
qui nest autre que la quadri-norme de lenergie-impulsion. . .

(2.20)

2.5. EXPRESSIONS DU LAGRANGIEN EN FONCTION DE LESPACE-TEMPS

2.5.3

29

Remarques
epist
emologiques

Que retenir de cette approche? Quil a su dassocier au principe de Hamilton un a priori (principe de
relativite de Galilee ou de Lorentz) sur la structure de lespace-temps dans lequel se produit un evenement
pour nous permettre de reconstruire lensemble des outils necessaires `a la description de la dynamique.
Certaines constantes apparues lors de cette procedure ont ete ensuite identiees `a des grandeurs communes
(ex, la masse m, ou le produit mc), en imposant simplement une continuite de notre description.
Si nous avons des raisons de croire que lespace-temps est encore plus complique (par ex, dote dune
metrique complexe ds2 = g dx dx ou de proprietes geometriques particuli`eres), alors la procedure reste
absolument identique et aboutirait `a de nouvelles equations de la mecanique. Il faudrait ensuite comparer
les predictions de cette theorie avec les experiences.

30

CHAPTER 2. PRINCIPE VARIATIONNEL

Chapter 3

M
ecanique de Hamilton
Ce que nous avons vu du formalisme lagrangien sut amplement `a traiter lensemble des probl`emes de la
mecanique classique (et relativiste). Lapproche de Hamilton, que nous allons developper dans la suite du
cours, napporte rien de nouveau du point de vue du contenu physique. Mais elle ore un cadre theorique puissant, permettant une interpretation geometrique de la mecanique. Cest dans ce cadre que sest developpee
la mecanique quantique et la physique moderne (en particulier la theorie des champs) et cest dans ce cadre
egalement que setudient tous les phenom`enes de chaos.

3.1

Hamiltonien dun syst`

eme

Nous savons modeliser tout syst`eme `a n degres de liberte qk soumis `a des forces conservatives, grace `a un
lagrangien L(qk , qk , t) = T V . Dans cette formulation, qk = dqk /dt et donc depend `a priori de qk , de meme
que limpulsion generalisee pk = L
qk .
Mais, in ne, lexpression qk (t) depend de sa valeur initiale qui est independante des qk . Letat complet
dun syst`eme depend ainsi des positions

q (t = 0) et des vitesses
q(t
= 0) initiales qui sont, elles, totalement
independantes les unes des autres. Cest notre formalisme de la mecanique qui a cree ce lien, les qk et les pk
etant, en denitive, deux jeux de coordonnees independantes.
Imaginons par exemple que nous voudrions connaitre tous les comportements possibles dun syst`eme
dynamique: il sut pour cela de specier de facon independante les positions et les vitesses initiales. Cela
peut se faire aisement dans le formalisme lagrangien, mais naurait-on pas plutot interet `a formuler la
mecanique de telle sorte que qk et pk soient demblee independants?
Travailler sur certaines variables puis en changer pour dautres plus pertinentes (tout en conservant la
notion de dierentielle totale) est une demarche courante en thermodynamique: on appelle cela faire une
transformee de Legendre. On cherche ainsi `a obtenir une fonction g(q, p, t) construite `a partir de L(q, q,
t).
La facon la plus simple est de chercher une fonction h triviale, telle que
g(q, p, t) = L(q, q,
t) + h(q, q,
p, t)
Si on dierencie cette expression, on obtient


 g
g
g
dt
dqk +
dpk +
dg =
qk
pk
t
k



 


 L
L
L h
h
h
h
=
+
+
+
dpk +
dqk +
dqk +
dt
qk
qk
qk
qk
pk
t
t
k

ce qui se traduit par les contraintes suivantes

g
qk

L
h
+
qk
qk
31


CHAPTER 3. MECANIQUE
DE HAMILTON

32

h
pk
L
h
0 =
+
qk
qk
L h
g
=
+
t
t
t

La troisi`eme equation permet dobtenir h(q, q,
p, t) = k pk qk + f (q, p, t). On voit ensuite quil existe une
solution triviale en posant f = 0. On obtient alors


pk qk
g(q, p, t) = L(q, q,
t)
g
pk

On reconnait alors, au signe pr`es, lintegrale premi`ere obtenue dans le formalisme lagrangien et associee `a la
translation dans le temps. On choisit donc pour nouvelle fonction, lexpression


H(q, p, t) =
pk qk L
(3.1)
k

appelee le hamiltonien du syst`eme.

3.2

Equations canoniques de Hamilton

Grace au hamiltonien H, les coordonnees q et les moments conjugues p sont independants. Reste donc `a
reformuler les equations de Lagrange avec cette nouvelle grandeur. Par construction, on a


 H
dH
H
H
=
qk +
p k +
dt
qk
pk
t
k

Par ailleurs, la denition de H fournit

dH
dt




L
L
L
=
qk
qk
p k qk + pk qk
qk
qk
t
k




L
L
=
qk + qk p k

qk
t
k

ce qui, par identication et puisque p k =

L
qk

(equations de Lagrange), donne les equations suivantes

qk

p k

H
pk
H

qk

(3.2)

Ce jeu dequations est appele equations canoniques de Hamilton (canoniques car simples et symetriques).
Elles decoulent directement des equations de Lagrange et de la denition du hamiltonien. Nous avons
remplace n equations dierentielles du second ordre (Lagrange) par 2n equations dierentielles du premier
ordre (Hamilton), ce qui est un gain appreciable en termes de resolution mathematique!
On peut donc choisir de resoudre un probl`eme mecanique en utilisant le formalisme hamiltonien. Pour
cela, il faut
(1) Etablir le lagrangien L(q, q,
t) = T (q, q,
t) V (q, q,
t)
(2) Calculer les moments conjugues pk = L
 qk
(3) Calculer le hamiltonien H(q, p, t) = k pk qk L
(4) Resoudre les 2n equations canoniques de Hamilton.

3.3. PRINCIPE VARIATIONNEL

33

Lensemble des proprietes de symetrie et des lois de conservation veriees par le lagrangien se retrouvent
dans lhamiltonien. Le formalisme hamiltonien apparait cependant particuli`erement adapte au traitement
egrale premi`ere.
des coordonnees cycliques: si qi est cyclique, cest `a dire si H
qi = 0, alors pi est une int
Enn, le hamiltonien obeit `a lequation suivante
dH
H
L
=
=
dt
t
t

(3.3)

et, on la vu, se confond avec lenergie (pertinente) totale du syst`eme H = E = T + V pour des syst`emes
fermes (autonomes) et V(q).

3.3

Principe variationnel

Nous allons montrer ici que les equations de Hamilton se deduisent egalement dun principe variationnel.
Laction est en eet denie par



pk qk H)dt
S = L dt = (
k

Le principe de Hamilton stipule que laction est stationnaire (S = 0) lors dune variation des chemins q et
p entre deux points xes. La variation daction est (pour n=1 pour alleger lecriture)


dq
S = (pq H)dt = (pq + p H)dt
dt

d(q) H
H
=
(qp
+p

q
p)dt
dt
q
p

 
H
H
2
=
(q
)p (p +
)q dt + [pq]1
p
q
o`
u le dernier terme est nul puisque tous les chemins passent par les extremites. Ainsi, une action extremale
pour des variations q et p independantes correspond aux equations suivantes
q

H
p
H

qui sont bien les equations de Hamilton. On remarque accessoirement que les termes coordonnees et
moments sont trompeurs car ils semblent donner plus dimportance aux coordonnees, alors que dans le
formalisme hamiltonien les q et les p jouent un role equivalent.

3.4
3.4.1

Etude dun cas simple: pendule 1D

Ecriture de lhamiltonien

Nous avons vu que le potentiel dun tel pendule secrit (`a une constante sans importance pr`es)
V = mgl cos
ce qui fournit un lagrangien L = 12 ml2 2 + mgl cos . La coordonnee generalisee choisie est un angle q =
En posant I = ml2 et 2 = g/l, on obtient un hamiltonien de la
et le moment conjuge est p = L
= ml2 .

forme
p2
2 I cos q
(3.4)
H(q, p) = T + V =
2I


CHAPTER 3. MECANIQUE
DE HAMILTON

34

La connaissance de la dynamique de ce syst`eme passe a priori par la resolution des equations de Hamilton.
Mais avec ce formalisme, on a maintenant acc`es `a une information plus globale du syst`eme. En eet, ce
syst`eme est conservatif car H ne depend pas explicitement du temps: H = E est une integrale premi`ere. La
nature du mouvement suivi par le pendule va alors dependre de la valeur prise par E, ce qui signie que lon
peut prendre E comme un param`etre. On va voir que, sans resoudre les equations de Hamilton, on a acc`es
`a des informations importantes.

3.4.2

Le portrait de phase

On peut representer la dynamique de ce syst`eme par des courbes diso-energie,



E
p = I 2 cos q + 2
I

(3.5)

dans un espace `a 2 dimensions dont les axes sont p et q, appele lespace des phases. Trois cas se presentent
alors:
(1) 0 < E < 2 I; il ny a de solution que pour cos q E2 I , autrement dit que pour un domaine borne
de q: on doit assister `a des oscillations libres. Cest ce quon appelle un mouvement de libration. Les
points o`
u p = 0 sont appeles points tournants ou points de rebroussement.
(2) E > 2 I; il y a une solution quelle que soit la valeur de q, par contre p reste bornee: on assiste `a un
mouvement de rotation (on dit aussi circulation) complet du pendule. La courbe p(q) est q0 -periodique
(ici q0 = 2).
eparant les 2 regimes precedents (qui sont topologiquement
(3) E = 2 I; correspond au cas limite
s

eparatrice.
dierents). La courbe decrite par p = I 2 1 + cos q sappelle la s

3.4.3

Etude au voisinage de points particuliers

Quelles sont les positions dequilibre de ce syst`eme? Il y en a deux, lune stable `a qe = 0 et lautre instable `a
qe = . Faisons une etude `a proximite de ces equilibres en posant x = q qe pour x petit et developpons
le potentiel `a lordre 2 (approximation harmonique).
En qe = 0, on obtient
2I 2
p2
+
x
E=
2I
2
Ainsi, la trajectoire est une ellipse et le pointqe = 0 est dit point elliptique. Remarquons quen faisant un
changement de coordonnees approprie, Q = Ix, P = L
, le hamiltonien devient plus simple
Q
E=

2
(P + Q2 )
2

et fournit lequation dun cercle. Les equations de Hamilton,

x
p

p
H
=
p
I
H
=
cest `a dire
x
=

Ix
= 2 Ix

fournissent une solution oscillante x(t) = Re(Aeit ) (A nombre complexe).

En qe = , on obtient
2I 2
p2

x = (P 2 Q2 )
E=
2I
2
2
cest `a dire lequation dune hyperbole, do`
u lappellation de point hyperbolique pour qe = . Les axes
de lhyperbole sont les separatrices elles-memes, p = Ix. Les equations de Hamilton fournissent une
solution divergente x(t) = Re(Aewt ).


3.5. THEORIE
DE HAMILTON-JACOBI

3.4.4

35

Remarques dordre g
en
eral

Que faut-il retenir de general dans ce cas particulier?

(1) Dans le cadre hamiltonien, letat complet dun syst`eme `a n degres de liberte est represente par un
point dans lespace des phases, espace `a 2n dimensions qi , pi .
(2) Le trace du portrait de phase permet de voir graphiquement le comportement dun syst`eme, sans
meme resoudre les equations du mouvement. Il ore ainsi une vision globale de sa dynamique. Chaque portrait de phase depend du potentiel V(q). Cette approche de geom`etre associe ainsi une equation dierentielle
`a une courbe (dite courbe integrale) dans lespace des phases.
(3) Au voisinage dun point dequilibre stable ou elliptique (minimum local de V), les trajectoires sont
des courbes fermees, des ellipses. Au voisinage dun point dequilibre instable ou hyperbolique (maximum
local de V), les trajectoires se comportent comme des hyperboles. Les deux types de solutions occupent des
regions distinctes de lespace des phases, separees par des courbes appelees separatrices.
(4) Nous venons de voir un point important: letude locale (au voisinage des points particuliers) a fourni
des cartes locales, valables uniquement dans un domaine restreint. Mais les solutions etant uniques en des
points reguliers (theor`eme de Cauchy), on peut faire un prolongement analytique des solutions `a lexterieur
de leur domaine. En pratique, on peut tracer qualitativement le portrait de phase dun syst`eme en etudiant
seulement son comportement au voisinage des points singuliers.
(5) Poincare a introduit une classication des points singuliers:
Centre ou point elliptique: toutes les solutions forment `a son voisinage des courbes fermees (ellipses).
Col ou selle ou encore point hyperbolique: les solutions sont des hyperboles `a son voisinage.
Noeud: cest un point traverse par une innite de solutions.
Foyer: cest un point vers lequel les solutions convergent, formant des spirales.
(6) Lintroduction dune dissipation aurait comme consequence la diminution de lenergie sur une echelle
de temps . La trajectoire de ce pendule dissipatif peut sapprehender grace au portrait de phase du pendule
conservatif, de periode T . Si  T , la trajectoire va converger tr`es rapidement vers le point centre. Chaque
trajectoire va dependre des conditions initiales: on peut observer soit une courbe pratiquement rectiligne,
soit une forme en spirale. Si T , il faudra beaucoup de periodes pour que le pendule perde de lenergie.
Ainsi, partant par exemple dune solution circulante, on arriverait, de proche en proche `a une solution du
type libration, dont lamplitude decroit au l du temps. Le pendule dissipatif poss`ede ce quon appelle
un attracteur: cest le point elliptique (0,0). Un autre dexemple dattracteur apparait pour un syst`eme
soumis `a laction dun moteur. Le syst`eme, apr`es une phase transitoire liee `a sa dynamique interne, sera
force de suivre un certain comportement impose par le moteur. Ce comportement dessine une trajectoire
dans lespace des phases qui nest rien dautre quun attracteur.

3.5

Th
eorie de Hamilton-Jacobi

Nous avons tout dabord montre (chapitre I) que tout syst`eme mecanique ayant n degres de liberte q, soumis
ou non `a des contraintes holonomes, sur lequel sexercent des forces conservatives peut etre decrit par un
lagrangien L(q, q,
t). Les equations de la dynamique sont constituees de n equations dierentielles du second
ordre, les equations de Lagrange. Letat du syst`eme `a un instant donne t est decrit par un point dans
lespace des congurations (dim n). On vient de voir que lon peut remplacer ces equations par les equations
canoniques de Hamilton, 2n equations dierentielles du premier ordre seulement, portant sur le hamiltonien
H(q, p, t) = pq L. Letat complet du syst`eme `a un instant donne t est alors deni par un point dans
lespace des phases (dim 2n). Ces deux approches sont parfaitement equivalentes, elles decoulent toutes
deux du meme principe variationnel (chapitre II).


CHAPTER 3. MECANIQUE
DE HAMILTON

36

Mais la simplicite meme des equations canoniques de Hamilton

H
qk =
pk
H
p k =
qk
sugg`ere une methode generale de resolution des probl`emes mecaniques pour des syst`emes dont les forces
derivent dun potentiel (meme generalise).
En eet, rechercher levolution temporelle dun syst`eme est equivalent, dans le formalisme hamiltonien,
`a faire un changement de coordonnees (p, q) (P, Q) dans lespace des phases. Envisageons deux cas pour
se convaincre de ceci.
Soit H
(Q, P, t) le nouvel hamiltonien exprime avec les nouvelles variables. Supposons que lon trouve
des variables telles que H
= 0 (alors que H(q, p, t) = 0). Dapr`es les equations canoniques de Hamilton on
aurait alors
H

Q k =
=0
Pk
H

Pk =
=0
Qk
ce qui implique que toutes les Q ainsi que les P sont des invariants. Un syst`eme `a n degres de liberte
aurait ainsi 2n invariants ! Cela peut paraitre etrange mais cest eectivement toujours le cas: ce sont les 2n
conditions initiales. Ainsi, notre changement de variable (p, q) (P, Q) permet, apr`es inversion, dexprimer
qk (t) =

qk (Q, P, t) = qk (q1 (0), . . . , qn (0), p1 (0), . . . , pn (0), t)

pk (t) =

pk (Q, P, t) = pk (q1 (0), . . . , qn (0), p1 (0), . . . , pn (0), t)

Dun point de vue geometrique, cela signie que nous avons fait un changement de coordonnees qui, `a un
point de lespace des phases caracterise par le param`etre t = 0, associe un nouveau point de param`etre t. Le
probl`eme consiste `a trouver ce fameux changement de variables. . .
Dans lexemple ci-dessus, le hamiltonien peut dependre explicitement du temps. Considerons maintenant
le cas plus simple dun syst`eme conservatif H = H(q, p). Si lon trouve de nouvelles variables (Q, P ) telles
que H
= H
(P ), autrement dit telles que les Q soient toutes cycliques, alors les equations de Hamilton
donnent
H

P k =
=0
Qk
H

Q k =
= k (P )
Pk
Cela signie que les n nouveaux moments Pk sont des invariants tandis que levolution temporelle des Qk
est triviale
Qk = k (P )t + Qk,0
Il sut ensuite dexprimer (q, p) en fonction de (Q, P ) pour avoir levolution temporelle des anciennes coordonnees.
Cette approche geometrique est connue sous le nom de theorie de Hamilton-Jacobi. Nous verrons au
chapitre suivant sur les syst`emes hamiltoniens que chercher un changement de variables tel que
(i) H
= 0, nous conduit `a lequation de Hamilton-Jacobi;
(ii) toutes les Q sont cycliques, nous conduit aux variables canoniques dangles et daction.
Mais le raisonnement que nous avons suivi ne tient que si, avec les nouvelles variables, on peut encore
ecrire
H

Q k =
Pk
H

Pk =
Qk

3.6. TRANSFORMATIONS CANONIQUES

37

autrement dit, qu`a la seule condition que la structure formelle des equations de hamilton soit conservee.
Avant donc daborder la mani`ere de rechercher les bons changement de variables, il faut nous munir doutils
permettant de manipuler des grandeurs quelconques dans lespace des phases: transformations canoniques
et crochets de Poisson.

3.6
3.6.1

Transformations canoniques
Fonctions g
en
eratrices

Jusqu`a present, nous navons considere que des transformations de coordonnees dites ponctuelles, cest `a
dire de la forme
Qi = Qi (qk , t)
(par exemple: passage cartesiennes en polaires). Dans le formalisme de Hamilton, les pk sont egalement des
coordonnees, il faut donc elargir le concept de transformation. On appelle transformation de contact toute
transformation de la forme
Qi
Pi

=
=

Qi (qk , pk , t)
Pi (qk , pk , t)

Avec ce nouveau jeu de variables, on doit pouvoir ecrire de nouvelles equations de Hamilton, portant sur un
nouvel hamiltonien H
(P, Q, t). Or, cet hamiltonien decrit le meme syst`eme physique que H(p, q, t) et donc,
obeit au meme principe variationnel,



S = Ldt = (pk qk H)dt = (Pk Q k H
)dt
H et H viennent donc du meme lagrangien, `a une derivee totale par rapport au temps pr`es, cest `a dire
dG
H
= H + Pk Q k pk qk +
dt

(3.6)

D
enition: une transformation canonique est une transformation de contact satisfaisant la condition
(3.6). On appelle G la fonction generatrice de la transformation canonique, denie sur lespace des phases
du syst`eme.
Limportance de la fonction generatrice reside dans le fait que sa connaissance determine compl`etement
la transformation canonique. Ceci nest pas evident et resulte de la contrainte canonique.
A 1 degre de liberte, il y a 4 classes de fonctions generatrices interessantes melant les deux jeux de
coordonnees: G1 (q, Q, t), G2 (q, P, t), G3 (p, Q, t) et G4 (p, P, t). Le choix depend des circonstances. A n
degres de liberte, le nombre de fonctions generatrices melant les anciennes et les nouvelles coordonnees est
evidemment beaucoup plus eleve, mais la methode qui suit reste toujours la meme.
Supposons que nous ayons une relation du type G1 (q, Q, t), cest `a dire que nous sachions exprimer
p = p(q, Q, t) et P = P (q, Q, t). La condition (3.6) signie que si on les remplace dans le membre de gauche
et dans le membre de droite, on arrive `a une egalite fonctionnelle quelles que soient les valeurs des variables
independantes q et Q. Plus precisement, on doit avoir
pdq Hdt = P dQ H
dt +

G1
G1
G1
dq +
dQ +
dt
q
Q
t

cest `a dire
(p

G1
G1
G1
)dq + (H
H
)dt (P +
)dQ = 0
q
t
Q


CHAPTER 3. MECANIQUE
DE HAMILTON

38

pour des variations independantes des variables q et Q et du param`etre t. On a donc les equations suivantes:
p

H

G1
q
G1

Q
G1
H+
t

Par souci de completude, regardons ce quil advient des autres fonctions generatrices. Par exemple pour
G2 (q, P, t), on aurait
G2
G2
G2
(p
)dq + (H
H
)dt P dQ
dP = 0
q
t
P
Or, ici nous sommes supposes savoir exprimer Q en fonction de (q, P, t), on peut donc en theorie eliminer le
terme P dQ. Pour cela, il sut de prendre pour fonction generatrice la fonction G = G2 P Q, ce qui revient
`a prendre la transformee de Legendre de G2 . On obtient alors
p

H

H

H

G2
q
G2
P
H+

G2
t

pour G3 (p, Q, t) + pq,

G3
Q
G3

p
G3
H+
t

et pour G4 (p, P, t) + pq P Q
G4
p
G4
P
G4
H+
t

Remarques:
(1) Si lon desire que H
soit un invariant, on choisira des fonctions generatrices telles que G
t = 0.
(2) Meme si lon ne trouve pas des nouvelles coordonnees o`
u toutes les Qi sont cycliques, il est avantageux
de rechercher celles qui en orent le maximum.

3.6.2

Quelques transformations canoniques remarquables

Identit
e
Soit la transformation canonique suivante
G2 (q, P ) =

i

qi Pi

(3.7)

3.7. LES CROCHETS DE POISSON

39

Dapr`es les relations obtenues plus haut, on a

G2
= Pk
qk
G2
Qk =
= qk
Pk
G2
=H
H
= H +
t
On dit que G2 engendre la transformation identite.
pk

Transformations ponctuelles
Une transformation canonique de la forme
G2 (q, P ) =

fi (q, t)Pi

(3.8)

fournit

fi
G2
=
Pi
qk
qk
i

pk

G2
= fk (q, t)
Pk
G2
H
= H +
t
engendre une transformation ponctuelle puisque les nouvelles coordonnees Qk ne dependent pas des moments.
Les fi etant arbitraires, on en deduit que toute transformation ponctuelle est canonique (ce qui netait pas
evident).
Qk

Echange du r
ole
Soit la transformation
G1 =

qi Qi

(3.9)

fournissant les relations de passage suivantes

pk
Pk
H

G1
= Qk
qk
G1
=
= qk
Qk
= H
=

Cette transformation eectue un echange entre coordonnee generalisee et moment. En mecanique hamiltonienne, il faut perdre lhabitude de considerer q comme une coordonnee spatiale et p comme une impulsion.

3.7
3.7.1

Les crochets de Poisson

D
enition

Une nouvelle op
eration
Soit une fonction f (q, p, t) denie dans lespace des phases. Sa derivee par rapport au temps est


 f
f
df
f
=
qk +
p k +
dt
qk
pk
t
k


CHAPTER 3. MECANIQUE
DE HAMILTON

40

Remplacant qk et p k par les equations de Hamilton, on obtient une nouvelle expression qui peut se mettre
sous la forme compacte
f
df
= {f, H} +
(3.10)
dt
t
o`
u lon a introduit la notation suivante, appelee crochet de Poisson pour H et f,


 f H
f H
{f, H} =

(3.11)
qk pk
pk qk
k

Etablie `a partir des equations de Hamilton, lequation (3.10) porte en elle la meme information sur la
dynamique dun syst`eme. Mais linteret essentiel de lequation (3.10) reside dans le calcul des integrales
erier
premi`eres. En eet, si f est une integrale premi`ere, df
dt = 0 et f doit donc v
{f, H} +

f
=0
t

(3.12)

Si f = f (q, p) ne depend pas explicitement du temps (ce qui est toujours vrai pour des syst`emes autonomes),
alors la condition pour etre une integrale premi`ere se ram`ene `a
{f, H} = 0

(3.13)

Ainsi, une facon de verier si une expression quelconque des variables dynamiques est un invariant est de
calculer son crochet de Poisson avec le hamiltonien. Les crochets de Poisson nous orent donc un test general
pour la recherche et lidentication des constantes du mouvement.
D
enition g
en
erale
Soient deux fonctions f (q, p) et g(q, p) quelconques denies dans lespace des phases, les crochets de Poisson
sont par denition


 f g
f g
{f, g} =

(3.14)
qk pk
pk qk
k

A partir de cette denition generale, on peut construire des crochets dits fondamentaux, en prenant pour
fonctions f et g les variables qk et pk . On obtient alors les relations suivantes


 pi pj
pi pj
{pi , pj } =

=0
qk pk
pk qk
k

{qi , qj } =
{qi , pj } =

0



 qi pj
qi pj

ik jk = ij
=
qk pk
pk qk
k

3.7.2

(3.15)

Propri
et
es

Tout un ensemble de proprietes decoulent de la denition meme des crochets de Poisson:

{f, g} = {g, f }
{f, c} = 0 o`
u c est une constante
{f1 + f2 , g} = {f1 , g} + {f2 , g}
{f1 f2 , g} = f1 {f2 , g} + f2 {f1 , g}
f
g

{f, g} = { , g} + {f,
}
t
t
t
f
{f, qi } =
pi
f
{f, pi } =
qi

(3.16)
(3.17)
(3.18)
(3.19)
(3.20)
(3.21)
(3.22)

3.7. LES CROCHETS DE POISSON

41

Enn, il est aise (bien que penible) de verier que les crochets de Poisson satisfont lidentite de Jacobi
{f, {g, h}} + {h, {f, g}} + {g, {h, f }} = 0

(3.23)

Une application importante resulte directement de cette identite. En eet, soient deux fonctions quelconques f et g denies dans lespace des phases. On a dapr`es (3.10)
d
{f, g} =
dt
=
=
=

{{f, g}, H} +

{f, g}
t

{f, g}
t
df
f

dg g

} {g,

} + {f, g}
{f,
dt
t
dt
t
t
dg
df
{ , g} + {f, }
dt
dt
{f, {g, H}} + {g, {H, f }} +

Th
eor`
eme de Poisson: Le crochet de Poisson de deux invariants est lui-meme un invariant.
dg
d
La demonstration est au-dessus: si f et g sont des invariants alors df
dt = dt = 0 et donc dt {f, g} = 0.
Note: cette methode pour trouver de nouveaux invariants a cependant ses limitations, car on tombe
souvent sur des fonctions triviales sans interet.

3.7.3

Invariance canonique

Th
eor`
eme: Les crochets de Poisson sont independants du syst`eme de coordonnees canoniques dans lequel
ils sont exprimes, autrement dit
(3.24)
{f, g}q,p = {f, g}Q,P
La demonstration compl`ete (tr`es calculatoire) ne sera pas donnee, je me contente de largumentation
donnee par Landau (et qui constitue egalement une demonstration).
Lorsquon fait une transformation canonique (q, p) (Q, P ), le temps nintervient jamais explicitement.
Il ne joue eventuellement que le role dun param`etre. En consequence, si on demontre (3.24) pour des
grandeurs ne dependant pas explicitement du temps, le theor`eme sera egalement vrai dans le cas general
(puisquil sagit dune propriete independante du temps). Par ailleurs, on peut toujours considerer que g est
formellement le hamiltonien dun syst`eme ctif. En vertu de lequation (3.10) on obtient alors (pour g = H)
df
= {f, H}q,p
dt
Or, le taux de variation de f ne peut dependre du syst`eme de coordonnees choisi. Cela implique donc
{f, g}q,p = {f, g}Q,P . On ne mettra donc plus les indices q, p aux crochets.
Regardons ce que cela implique. Soit une transformation canonique q = q(Q, P ) et p = p(Q, P ). Les
crochets de Poisson de deux fonctions quelconques secrivent


 f g
f g
{f, g}q,p =

qk pk
pk qk
k





f
 g Qi
g Pi
g Pi
f
g Qi
=
+
+

qk i Qi pk
Pi pk
pk i Qi qk
Pi qk
k





g
f Qi

g
f Pi
f Qi
f Pi
=

+
Qi
qk pk
pk qk
Pi
qk pk
pk qk
i
i
k
k

g

g
=
{f, Qi }q,p +
{f, Pi }q,p
Qi
Pi
i
i


CHAPTER 3. MECANIQUE
DE HAMILTON

42

Cette derni`ere relation va nous etre tr`es utile. En remplacant f par Qi et g par f , on obtient
{f, Qi }q,p = {Qi , f }q,p =

f

f
{Qi , Qj }q,p
{Qi , Pj }q,p
Qj
Pj
j
j

o`
u lon reconnait des crochets fondamentaux. La condition (3.24) sapplique evidemment aux crochets
fondamentaux, ce qui fournit
{Qi , Qj }q,p = 0 {Qi , Pj }q,p = ij

{Pi , Pj }q,p = 0

(3.25)

On obtient ainsi
{f, Qi }q,p

{f, Pi }q,p

f
Pi
f
Qi

ce qui nous fournit bien

{f, g}q,p =

g f

g f
+
= {f, g}Q,P
Qi Pi
Pi Qi
i
i

Ainsi, on voit quil sut que les crochets fondamentaux de Poisson soient des invariants canoniques (ie.
laisses invariants lors dune transformation canonique) pour que cela reste vrai pour des fonctions f et g
quelconques. On peut donc eriger en theor`eme la proposition suivante.
Th
eor`
eme: Une transformation sera canonique si elle verie (3.24) pour des fonctions f et g quelconques
ou si, de mani`ere equivalente, les relations (3.25) sont satisfaites.

3.7.4

Interpr
etation g
eom
etrique

Nous avons vu dans le formalisme lagrangien quil y avait un lien entre symetrie et invariants. Le theor`eme
de Emmy Noether nous a meme fourni une methode permettant de calculer linvariant lorsque nous savons
quel changement de variable operer.
Nous allons maintenant voir que ce lien entre invariants et proprietes de symetrie est une propriete
geometrique de lespace des phases1 .
G
en
erateurs de transformations canoniques innit
esimales
Toute operation de symetrie se traduit par un changement de variables (nous ne considerons evidemment
que des transformations canoniques). Soit f (q, p, t) une fonction quelconque denie sur lespace des phases.
La modication due `a une transformation innitesimale ne dependant pas explicitement du temps est
f

= f (q + q, p + p, t) f (q, p, t)


 f
f
=
qi +
pi
qi
pi
i

o`
u les qi et pi ne sont pas des variations au sens du calcul variationnel mais des modications elementaires.
Soit la transformation canonique innitesimale
Qi
Pi

=
=

qi + qi
pi + pi

1 Ce faisant, nous allons g

eneraliser le theor`eme de Noether de la meme mani`ere quune transformation de contact generalise
la notion de transformation ponctuelle.

3.7. LES CROCHETS DE POISSON

43

quelle est la fonction

generatrice G associee? Nous avons vu que la transformation identite est engendree

par G2 (q, P ) = i qi Pi . On peut donc ecrire
G(q, P ) =

qi Pi +
F

(3.26)

o`
u
est un petit param`etre et F une fonction encore inconnue. Or, G doit etre une transformation canonique,
= G P Q)
elle verie donc (en posant G

pdq Hdt = P dQ H
dt + dG
quelles que soient les variations des Pi et des qi . On obtient ainsi les relations
pi

Qi

H

G
F
= Pi +

qi
qi
G
F
= qi +

Pi
Pi
H

do`
u
pi

qi

F
qi
F
F
Qi qi =

Pi
pi
Pi pi =

(3.27)

Ces expressions doivent etre toutes deux du premier ordre en

, ce qui etait bien le cas pour pi , alors que
qi mettait en jeu une derivation par rapport `a Pi . On voit donc quil sut de se donner une fonction
F (q, p), dependant des anciennes variables, pour engendrer une transformation canonique innitesimale de
param`etre
.
Forts de ce resultat, on obtient que au premier ordre


 f F
f F

f =

qi pi pi qi
i
=
{f, F }
Loperation {f, F } a donc une signication precise: f est la modication apportee `a la fonction f lors
dun changement innitesimal de coordonnees canoniques, de fonction generatrice F . Cette transformation
innitesimale est caracterisee par un unique param`etre
.
Relation avec les invariants dun syst`
eme
Th
eor`
eme: Les constantes du mouvement dun syst`eme autonome sont les fonctions generatrices des transformations canoniques innitesimales qui laissent H invariant.
Soit G une transformation canonique quelconque (ne dependant pas explicitement du temps). Si on a
{H, G} = 0

(3.28)

alors, cela signie dapr`es lequation (3.10) que G est une integrale premi`ere. Mais en vertu de la relation
precedente, cela signie aussi que H = 0, autrement dit, que G laisse H invariant.
Ce que nous savions dej`a grace `a Lagrange apparait ici comme un cas particulier dun ensemble tr`es vaste
de transformations canoniques possibles.


CHAPTER 3. MECANIQUE
DE HAMILTON

44
Exemple 1: translation dans le temps

Nous avons vu que lhamiltonien H dun syst`eme est invariant si ce dernier est symetrique par translation
dans le temps.
Prenons donc F = H et regardons ce que cela implique. Les relations (3.27) deviennent
pi

qi

H
=
p i
qi
H

=
qi
pi

ce qui indique que

= dt. Ainsi, le hamiltonien dun syst`eme est bien le generateur dune translation
innitesimale dans le temps, faisant passer celui-ci de linstant t `a linstant t + dt.
En consequence, le mouvement dun syst`eme entre deux instants t1 et t2 peut etre decrit par une succession
de transformations canoniques innitesimales dont le generateur est H.
Exemple 2: translation dun qi
Nous avons vu auusi que si une variable qk est cyclique (ex: invariance par translation dans une direction,
rotation autour dun axe, mais cela peut etre plus general), alors le moment conjugue pi est un invariant.
Prenons F = pk , les relations (3.27) deviennent
pi

qi

pk
=0
qi
pk

=
ik
pi

ce qui montre que pk est bien le generateur dune translation dans la direction qk dune quantite
. Limpulsion
est le generateur du mouvement de translation du syst`eme, tandis que le moment cinetique celui de rotation.

3.8

Lespace des phases

Nous avons vu que letat complet dun syst`eme `a un instant donne t est un point x(q1 , . . . , qn , p1 , . . . , pn )
dans un espace `a 2n dimensions, appele lespace des phases. Cet espace na pas la structure dun espace
vectoriel (cest une variete dierentiable, une classe dequivalence datlas). Lespace des phases va nous
permettre dapprehender tout un ensemble de proprietes formelles des syst`emes dynamiques et den tirer des
interpretations geometriques simples.

3.8.1

Flot hamiltonien

Levolution temporelle du syst`eme est regie par les equations de hamilton qui peuvent se mettre sous la
forme condensee suivante
t
(x, t)
(3.29)
x = gH
o`
u le vecteur

t
gH
=

H

p

H

q

represente le champ de vitesses au point x. En resolvant les equations de Hamilton, on dessine une
trajectoire dans lespace des phases, de la meme mani`ere quon le ferait en suivant de proche en proche le
t
le ot hamiltonien, par analogie avec la mecanique
vecteur vitesse dans un espace euclidien. On appelle gH
des uides.
Le ot hamiltonien a une structure de groupe (loi interne, associativite, element neutre, inverse).

3.8. LESPACE DES PHASES

45

= t0 ), la solution
Th
eor`
eme de Cauchy (non demontre): Pour des conditions initiales donnees x 0 = x(t
x(t)

de lequation (3.29) existe et est unique pour un temps t ni.

Ce theor`eme nest valable que pour des positions dites non-singuli`eres de lespace des phases. Un point
hyperbolique est un exemple de point singulier.
Consequences:
(1) Deux trajectoires dans lespace des phases ne peuvent se couper. Attention: un point commun `a deux
trajectoires signie le meme point de lespace des phases au meme instant. Si cetait le cas, alors cela
signierait que deux conditions initiales dierentes conduiraient `a un meme etat commun, ce qui est en
contradiction avec lunicite des solutions. A noter que cela ne concerne pas les trajectoires associees `a des
temps de parcours innis.
(2) Il est dicile de representer un espace `a plus de 2 dimensions (n=1). Lorsquon le fait, on fait une
projection et on peut alors y observer des trajectoires qui se coupent.

3.8.2

Incompressibilit
e du ot

Conservation du volume
t
, cest `a
Th
eor`
eme: un volume V quelconque de lespace des phases est conserve par le ot hamiltonien gH
dire
dV
=0
(3.30)
dt
Soit un point quelconque de V . Lors dune evolution temporelle, ce point va suivre une trajectoire precise
dans lespace ces phases, imposee par les equations canoniques de Hamilton. Il en va de meme pour tous les
autres points appartenant `a V . On va donc assister, entre les instants initial et nal, `a une deformation de V ,
analogue `a celle dun element de uide pris dans un ecoulement. Le theor`eme ci-dessus stipule simplement
que la deformation due au ot hamiltonien conserve le volume (ecoulement incompressible).
Nous avons dej`a vu que la variation de p et q au cours du mouvement peut etre consideree comme une
transformation canonique. Donc, montrer que le volume V reste invariant par le ot hamiltonien, revient `a
montrer que le volume est un invariant canonique.
Lespace des phases etant un espace `a 2n dimensions, le produit de dierentielles

dV = dq1 . . . dqn dp1 . . . dpn



peut etre considere comme un element innitesimal de volume et le volume V secrit V = . . . dq1 . . . dqn dp1 . . . dpn .
Pour n = 1, V est en fait une surface, tandis que pour n 2 V est un hypervolume.
Cest un invariant canonique si




. . . dQ1 . . . dQn dP1 . . . dPn
. . . dq1 . . . dqn dp1 . . . dpn =


=
. . . Jdq1 . . . dqn dp1 . . . dpn
o`
u

(Q1 , . . . , Qn , P1 , . . . , Pn )
J=
=

(q1 , . . . , qn , p1 , . . . , pn )

Q1
q1

...

Qn
q1

P1
q1

...

Pn
q1

..
.

...

..
.

..
.

...

..
.

Q1
pn

...

Qn
pn

P1
pn

...

Pn
pn

est le determinant de Jacobi ou jacobien. La demonstration du theor`eme revient `a demontrer que J = 1

pour toute transformation canonique.
La manipulation des jacobiens jouit de certaines proprietes. On peut en particulier reecrire J sous la
forme
J=

(Q1 ,...,Qn ,P1 ,...,Pn )

(q1 ,...,qn ,P1 ,...,Pn )
(q1 ,...,qn ,p1 ,...,pn )
(q1 ,...,qn ,P1 ,...,Pn )


CHAPTER 3. MECANIQUE
DE HAMILTON

46
Par ailleurs, on peut reduire le rang dun jacobien si

(x,y,a)
(u,v,a)

(x,y)
(u,v)

est veriee. Dans notre cas, on obtient

J=

(Q1 ,...,Qn )
(q1 ,...,qn ) P =Cst

(p1 ,...,pn )
(P1 ,...,Pn ) q=Cst

J1
J2

o`
u J1 = det(A) et J2 = det(B), les matrices A et B etant denies par
Aij =

Qi
qj

et

Bij =

pi
Pj

Si on choisit comme fonction generatrice de la transformation canonique G2 (q, P ), on a les relations suivantes
Qi

pi

G2
Pi
G2
qi

ce qui fournit
Aij =

2 G2
qj Pi

et

Bij =

2 G2
Pj qi

Ainsi, B = t A est la matrice transposee de A et det(B) = det(A), autrement dit J = 1.

Invariants int
egraux de Poincar
e
Soit la grandeur suivante
I1 =

 

dqi dpi

o`
u lintegration seectue sur une surface `a 2 dimensions S quelconque (en fait une variete) de lespace des
phases. On peut demontrer tr`es facilement que I1 est un invariant canonique. Soient deux coordonnees u et
v independantes, permettant de caracteriser tout point appartenant `a la surface S. Pour tous indices i et k
nous sommes capables deectuer les changements de variables suivants
qi

= qi (u, v)

pi
Qk

= pi (u, v)
= Qk (u, v)

Pk

= Pk (u, v)

o`
u les variables (Q, P ) sont reliees aux (q, p) par une transformation canonique. Linvariance canonique de
I1 est assuree d`es lors que
 

 

I1 =
dqi dpi =
Ji dudv
=

i
 

dQk dPk =

  i

autrement dit, si on montre que

i

Ji =

k

J
k

J
k dudv

3.8. LESPACE DES PHASES

47

i ,pi )


avec les jacobiens Ji = (q
(u,v) et Jk =
generatrice G2 (q, P, t), il vient

Ji =

(qi ,pi )
(qi ,Pi )
(u,v)
(qi ,Pi )

(Qk ,Pk )
(u,v) .

1
=
Ai

Utilisant les r`egles de calcul des jacobiens et la fonction

qi
qi
qi
Pi

pi
qi
pi
Pi

2

= 1 pi = 1 G2
Ai Pi
Ai Pi qi

Qk
qk
Qk
Pk

Pk
qk
Pk
Pk

2

= 1 Qk = 1 G2
Ak qk
Ak qk Pk

tandis que lautre jacobien devient

J
k =

(Qk ,Pk )
(qk ,Pk )
(u,v)
(qk ,Pk )

1
=
Ak



Il est donc evident que i Ji = k J
k .
On peut recommencer de la meme mani`ere pour la grandeur suivante
   

I2 =
dqi dpi
i

o`
u lintegration seectue sur une hypersurface S de dimension 4 (ce qui suppose un espace des phases de
dimension 2n avec n 2). Dans ce cas, il faut 4 coordonnees independantes pour caracteriser un point de
S et on reproduit la mme demonstration que ci-dessus.
Dune facon generale, Poincare a montre que toutes les integrales de la forme
 
 

Is =
...
dqi dpi
(3.31)
i

etendues `a des varietes S `a 2s dimensions sont des invariants canoniques, dans un espace des phases `a 2n
dimensions (n s).
Remarques: La conservation du volume nest quun cas particulier, obtenu pour s = n. Ce nest dailleurs
que pour ce cas l`a que lon nest pas obliges de passer par des variables intermediaires (comme u et v), puisque
tous les qi et pi sont mis en jeu (ce qui simplie la demonstration comme on la vu).

3.8.3

Th
eor`
eme de Liouville: lien avec la physique statistique

Th
eor`
eme: Soit D = dN/dV la densite detats dun syst`eme mecanique au voisinage dun point de lespace
des phases. Cette densite detats verie
dD
=0
(3.32)
dt
et est donc une grandeur qui se conserve lors de levolution du syst`eme.
Soit dN le nombre detats contenus dans le volume dV de lespace des phases. Nous avons vu que dV
est un invariant canonique, ce qui signie quil reste constant lors de levolution du syst`eme. Par ailleurs,
le nombre detats dN reste egalement constant, car aucune trajectoire ne peut traverser la surface fronti`ere
denissant le volume dV . En eet, si cetait le cas, alors deux conditions initiales dierentes m`eneraient au
meme etat (celui situe `a la fronti`ere), ce qui est en contradiction avec lunicite des solutions. Le quotient
dN/dV est donc une constante.
Le theor`eme de Liouville a surtout une importance en physique statistique, o`
u le nombre de particules
N est tr`es grand. Levolution temporelle de la densite detats (ie. densite de particules dans un etat donne)
secrit
D
dD
=
{D, H} = {H, D}
t
dt


CHAPTER 3. MECANIQUE
DE HAMILTON

48

A lequilibre statistique, le nombre de particules dans un etat donne doit etre constant. Comme dV reste
a dire
egalement constant, on obtient D
t = 0 cest `
{H, D} = 0

(3.33)

Pour assurer alors lequilibre statistique, il sut de choisir D comme fonction des constantes du mouvement. . .
Exemple: Ainsi lensemble microcanonique est obtenu pour D = D0 E0 , cest `a dire une constante pour
une energie donnee E0 , nulle sinon.

3.9
3.9.1

Syst`
emes int
egrables
Th
eor`
eme de Arnold-Liouville

De nos jours, il est facile de resoudre par ordinateur les equations du mouvement de nombreux syst`emes,
une fois choisies les conditions initiales. Mais tr`es souvent, ce nest pas cela que nous desirons.
Cette methode ne nous permet en eet que de calculer une trajectoire particuli`ere dans lespace des
phases, alors que nous pouvons desirer avoir une image de lensemble des comportements possibles du
syst`eme (comme pour le pendule). Par exemple: savoir si laxe de rotation de la Terre est stable ou chaotique,
connaitre les regimes de fonctionnement non chaotique des reacteurs etc. . .
Par ailleurs, il est des syst`emes (ex: probl`eme `a N corps) pour lesquels lintegration numerique est (encore) hors de portee et leur connaissance passe donc par une etude analytique. Enn, on peut ne sinteresser
qu`a certaines caracteristiques dun syst`eme, comme par exemple vouloir connaitre sa (ou ses) periode(s) et
savoir `a quelle(s) condition(s) celle(s)-ci reste(nt) stable(s).
D
enition: Un syst`eme est dit integrable si lon peut caracteriser qualitativement son comportement
(les trajectoires) dans lespace des phases.
Th
eor`
eme de Arnold-Liouville (1963) (non demontre): un syst`eme mecanique `a n degres de liberte
est integrable sil poss`ede les trois proprietes suivantes:
1- Il existe n integrales premi`eres Ii ;
2- Elles sont independantes;
3- Elles sont en involution.
La premi`ere propriete necessite la recherche au prealable des n integrales premi`eres dun syst`eme. Il faut
donc rechercher lensemble des proprietes dinvariance liees `a lespace-temps, utiliser les crochets de Poisson
et surtout se laisser guider par lintuition physique ou mathematique (symetries cachees). Si le syst`eme est
ferme, alors H = I1 est un invariant. Alors, Ii est une integrale premi`ere si
{Ii , I1 } = 0

i = 2, . . . , n

On a donc n 1 relations veriees par I1 .

La seconde propriete stipule lindependance des invariants: lespace forme par lintersection des surfaces
Ii = Cte doit etre de dimension n. Chaque invariant doit apporter une information supplementaire.
La troisi`eme propriete est assez contraignante. Les integrales premi`eres sont en involution si
{Ii , Ij } = 0

i, j n

(3.34)

autrement dit si Ii reste constante le long du ot gItj . Cela fournit n(n 1)/2 relations en tout, ce qui pour
n > 2 est plus contraignant que les n 1 relations. Pour n = 2, la condition dinvolution est equivalente `a
celle dindependance.
Une consequence immediate de ce theor`eme est que tout syst`eme conservatif `a 1 degre de liberte est
integrable. A plus de dimensions, les syst`emes integrables sont lexception plutot que la r`egle.

3.9. SYSTEMES
INTEGRABLES

3.9.2

49

Cartes et atlas symplectiques

Nous lavons vu `a maintes reprises, le bon choix des coordonnees (q, p) est celui o`
u la dynamique est la
plus simple. Faire une choix de coordonnees, cest se doter dune carte (langage de geographe). Mais cette
carte peut navoir dinteret quau voisinage dun point. Par exemple dans le cas du pendule `a 1D, nous avons
introduit une transformation ponctuelle ramenant le hamiltonien `a lexpression particuli`erement simple
H=

2
(P + Q2 )
2

mais valable uniquement au voisinage du point elliptique.

Pour representer correctement lensemble de lespace des phases dun syst`eme dynamique complexe, il est
en general necessaire davoir un ensemble de cartes. Elles correspondent par exemple aux divers changements
de coordonnees au voisinage des divers extrema (minima et maxima) locaux du potentiel et autres points
singuliers. On dit que ces cartes forment un atlas si elles poss`edent des zones de recouvrement, autrement
dit si elles sont compatibles (les relations de passage Q(q, p, t) et P (q, p, t) de lune `a lautre doivent etre
dierentiables).
Cet atlas sera dit symplectique si chacune de ces cartes correspond `a une transformation canonique,
autrement dit si, sur chacune delles, la structure formelle des equations de Hamilton est preservee.
Dapr`es ce que nous avons vu precedemment, les cartes les plus appropriees sont celles o`
u les coordonnees
sont des invariants du syst`eme. Il nous faudrait donc trouver un procede permettant dobtenir le plus grand
nombre possible dinvariants dun syst`eme: lutilisation des crochets de Poisson mais surtout la recherche
systematique des proprietes de symetrie du syst`eme devient necessaire. Cette demarche, initiee par les
travaux de Hamilton et de Jacobi, est celle qui perdure desormais en physique.

50

CHAPTER 3. MECANIQUE
DE HAMILTON

Chapter 4

Syst`
emes hamiltoniens
4.1

L
equation de Hamilton-Jacobi

4.1.1

La fonction principale de Hamilton

La facon la plus simple de sassurer quen faisant un changement de variables on tombe sur des nouvelles
variables P et Q cycliques, cest naturellement dimposer que le nouvel hamiltonien H
(P, Q) soit nul. Le
mouvement est alors trivial puisque
Q k

Pk

H

=0
Pk
H

=0
Qk

ce qui signie quil y a 2n invariants, fonctions des 2n conditions initiales. La condition de transformation
canonique (3.6) secrit alors


dG
=0
(4.1)
H(q, p, t)
pk qk +
dt
k

Ici, on a le choix de la fonction generatrice. Si lon choisit G = G2 (q, P, t), alors

G2

G2
dG2
G2
Pk +
=
qk +
dt
qk
Pk
t
k

2
Le deuxi`eme terme est nul par construction, puisque tous les Pk sont des invariants, tandis que pk = G
qk .
On obtient alors lequation suivante,
G
G
H(qk ,
=0
(4.2)
, t) +
qk
t

connue sous le nom dequation de Hamilton-Jacobi. La solution de cette equation est justement la fonction
G(q, t; P ), fonction de n + 1 variables (les n qk et le temps: les n Pk agissent comme des param`etres
puisquinvariants)1. Ainsi, la resolution de lequation de Hamilton-Jacobi nous fournit la fonction generatrice
de la transformation canonique souhaitee. La solution est appelee fonction principale de Hamilton.
Remarque: Cette approche reste valable meme si lhamiltonien de depart H(p, q, t) depend explicitement
du temps.
1 La fonction g
eneratrice G depend de n + 1 variables mais seulement de n invariants. Les invariants sont simplement les
constantes dintegration qui apparaissent lors de la resolution de lequation de Hamilton-Jacobi. A n + 1 variables correspondent
donc n + 1 constantes. Mais lune dentre elles est additive. Comme G nintervient que par ses derivees dans lequation de
Hamilton-Jacobi, cette constante nintervient pas et seules n sont pertinentes.

51

`
CHAPTER 4. SYSTEMES
HAMILTONIENS

52

4.1.2

Laction hamiltonienne

Une nouvelle vision de laction

Peut-on attribuer un sens `a la fonction principale de Hamilton? Remarquons que
dG
G
G

=
=
qk +
pk qk H = L
dt
qk
t
k

o`
u L est le lagrangien du syst`eme, lui meme etant la derivee totale de laction S par rapport au temps. La
fonction principale de Hamilton nest autre que laction hamiltonienne, denie par
 t
L dt
(4.3)
S(q, t; P ) =
t1

Laction hamiltonienne nest pas tout `a fait laction S(q) que nous avions introduit avec le principe variationnel, `a savoir
 t2
S(q) =
L dt
t1

Au chapitre II, laction etait vue comme une fonctionnelle des chemins q, rendue extremale par la trajectoire
reelle entre deux points xes q1 (t1 ) et q2 (t2 ). Ici, laction hamiltonienne S = S(q, t; P ) depend aussi du
temps (et des conditions initiales). On peut donc la voir comme une fonctionnelle des chemins q o`
u seul le
point de depart est maintenu xe.
Mais si on laisse le point darrivee libre, on perd une contrainte par rapport au calcul variationnel eectue
au chapitre II. Il faut donc la remplacer par une autre. Dans notre derivation de lequation de HamiltonJacobi, nous avions impose que notre changement de variables conservait les equations de Hamilton. Or
celles-ci sont equivalentes aux equations de Lagrange. Il est donc normal detre obliges de reutiliser cette
information. On va donc exiger que chaque chemin suivi soit une trajectoire possible mais aboutissant `
a un
point darrivee laisse libre.
Principe variationnel
Demontrons cette conjecture, `a savoir que lequation de Hamilton-Jacobi decoule bien dun principe variationnel, mais construit avec laction hamiltonienne.
Soit laction ecrite sous la forme suivante
 t
L(q, q,
t) dt
S(q, t; q1 , t1 ) =
t1

On part dun point xe dans lespace des congurations, mais le point darrivee est laisse libre. Par contre,
on suit la trajectoire reelle (et pas des chemins virtuels), ce qui signie que qi (t) doit obeir aux equations de
Lagrange. La dierentielle totale de laction est par denition

S
S
dt = Ldt
(4.4)
dqi +
dS =
qi
t
i
Le premier terme correspond `a une variation St de laction `a t constant, donc ce dont nous avons dej`a
lhabitude. Evaluons-le en faisant une variation des trajectoires
q =


, telles que
(t1 ) =
0. La variation
St associee est alors

 t

S
qi =
L dt
St =
qi t
t1
i


 t

 t

L
d L
L
qi +
qi dt =
qi dt
=
qi
qi
t1 i
t1 i dt qi



L
qi =
pi qi
=
qi
i
i


4.1. LEQUATION
DE HAMILTON-JACOBI

53

puisque t est xe et que qi (t1 ) = 0. On en deduit donc

S
= pi
qi

(4.5)

Repartant de la dierentielle totale de laction, on obtient

S
dt = Ldt
pi dqi +
dS =
t
i


S
L =
pi qi +
t
i
0

S
S
, t) +
qi
t

H(qi ,

(4.6)

qui est bien lequation de Hamilton-Jacobi.

Remarque: En suivant la demonstration, on voit quil a su dexiger q = 0 en t1 , rien concernant p.
En eet, il nest pas possible dexiger quoi que ce soit dautre, puisquon demande que chaque chemin suivi
corresponde `a une trajectoire reelle aboutissant `a un point darrivee dierent. Il faut donc se laisser libre
la direction de depart, `a savoir les pi `a t = t1 . Ceci se comprendra mieux plus loin (avec lanalogie de la
construction de Huygens).

4.1.3

M
ethode g
en
erale de r
esolution

Tout comme les equations de Lagrange ou les equations canoniques de Hamilton, lequation de HamiltonJacobi est le point de depart dune nouvelle methode generale de resolution des probl`emes mecaniques.
Montrons que le probl`eme est eectivement formellement resolu d`es lors que la fonction generatrice est
connue. Soit S(qk , t; Pk ) cette fonction, obtenue par resolution de lequation de Hamilton-Jacobi (4.2). On a
pk

Qk

S
qk
S
Pk

(4.7)
(4.8)

Les equations (4.7) sont valables `a tout instant. Ces n equations fournissent `a t = 0 un lien entre les
S
conditions initiales et les Pk . On connait donc, en principe, la valeur des Pk et par consequent de P
. Par
k
construction les Qk sont aussi des invariants. Les n equations (4.8) fournissent egalement `a t = 0 un lien
entre les Qk et les conditions initiales, permettant de connaitre ces Qk . On a presque ni: les expressions
pk (t) sont alors fournies par (4.7) tandis que les qk (t) sont obtenues en inversant le syst`eme des n equations
(4.8).
En pratique cependant, on utilise une autre facon de calculer les qk (t), mettant en jeu la resolution de n
ODE du premier ordre (voir les applications).
On a remplace les 2n ODE du premier ordre de Hamilton par la resolution dune PDE (equation aux
derivees partielles) du premier ordre en temps, non lineaire (si S1 et S2 sont solutions, S1 + S2 ne lest pas),
suivie ensuite par la resolution de n ODE du premier ordre et n inversions. Ce nest pas forcement un gain,
car la resolution dune PDE est en general tr`es dicile mathematiquement. Concr`etement, seuls les syst`emes
separables seront interessants pour cette methode.

4.1.4

M
ethode de s
eparation des variables

Notion de s
eparabilit
e
Si lequation de Hamilton-Jacobi peut se mettre sous la forme (apr`es avoir remplace lexpression de lhamiltonien
S
, t))
H(qi , q
i



S
S
S
, F1 q1 ,
n qi ,
, t,
=0
(4.9)
qi
t
q1

`
CHAPTER 4. SYSTEMES
HAMILTONIENS

54
alors on peut chercher une solution de la forme
S = S
(qi , t) + S1 (q1 )

pour i = 1

(4.10)

En eet, lequation (4.9) se met sous la forme




S

S

dS1
, F1 q1 ,
n qi ,
, t,
=0
qi
t
dq1
qui doit etre veriee quel que soit q1 . Cela implique que les termes ne dependant que de q1 sont constants,
cest `a dire


dS1
= 1
F1 q1 ,
dq1


S

n1 qi ,
, t; 1
= 0
qi
La premi`ere equation se resoud sans probl`eme (ODE du premier ordre), tandis que nous avons diminue de
un le nombre de variables intervenant dans la seconde. Il faut esperer pouvoir continuer cette tactique le
plus longtemps possible.
Dans le cas particulier o`
u le syst`eme est conservatif, le temps nintervient pas explicitement dans le
hamiltonien, H = E et lequation de Hamilton-Jacobi peut se lire S
egre facilement et
t = E, qui sint`
fournit
(4.11)
S(q, t; P ) = S0 (q; P ) Et
eristique de Hamilton. Elle est determinee
La nouvelle fonction S0 (q; P ) est parfois appelee fonction caract
par
S0
H(qi ,
)=E
qi
Un probl`eme sera enti`erement separable si la fonction caracteristique de Hamilton peut secrire


S0 =
Si (qi ; P1 , . . . , Pn )

(4.12)

Dans ce cas, la resolution de lequation de Hamilton-Jacobi se ram`ene `a la simple resolution de n ODE du

premier ordre de la forme


dSi
Fi qi ,
= i
dqi
o`
u les i sont des constantes.
Variables cycliques
Si une variable q1 est cyclique, elle napparait jamais dans le hamiltonien de telle sorte que lon doit resoudre


dS1
F1
= 1
dq1
S
La solution est evidemment S1 =
1 q1 . Or, on a egalement p1 = q
=
1 , ce qui determine la constante
1
(rien de nouveau: si q1 est cyclique alors p1 est un invariant). Laction hamiltonienne secrit alors

S(q, t; P ) = S
(qi , t) + p1 q1
o`
u lun des n invariants Pi vient detre identie: P1 = p1 .
generalisation du cas conservatif.

(4.13)

Noter egalement que ceci constitue une


4.1. LEQUATION
DE HAMILTON-JACOBI

55

Commentaires
Il nexiste pas de crit`ere simple pour savoir si lequation de Hamilton-Jacobi sera separable ou pas. La
separabilite depend du syst`eme de coordonnees choisi. La principale diculte consiste donc `a introduire des
variables commodes, mais nous ne sommes guides par aucune r`egle. . .
Ainsi, le probl`eme de la force centrale est separable en coordonnees polaires, mais pas en cartesiennes.
Le potentiel V = r F z decrivant un champ coulombien superpose `a un champ homog`ene est separable
en coordonnees paraboliques. Le potentiel V = r11 + r22 , decrivant les champs coulombiens crees par deux
centres immobiles, est separable en coordonnees elliptiques (Landau, section 48). Le probl`eme `a 3 corps est
non separable, mais la plupart des probl`emes en physique atomique le sont.

4.1.5

Applications `
a quelques probl`
emes simples

La particule libre 1D
Soit une particule libre de masse m avec une seul degre de liberte q. Son hamiltonien secrit
H(q, p) =

p2
2m

Lequation de Hamilton-Jacobi (HJ) est alors

1
2m

S
q

2
+

S
=0
t

Lhamiltonien ne dependant pas explicitement du temps, on a H = E conserve et la fonction generatrice

secrit S(q; P ) = S0 (q) Et. Lequation (HJ) fournit

dS0
= 2mE
dq

cest `a dire S(q; P ) = 2mEq Et. La fonction principale engendre une transformation canonique o`
u
toutes les variables sont cycliques, cest `a dire P = 0, Q = 0, donc des invariants. Linvariant evident ici est
E. Si on choisit P = E, alors la nouvelle coordonnee est

S
m
Q=
=
qt
E
2E
ce qui, par inversion, fournit


q=

2E
(t Q)
m

Le signe est impose par les conditions initiales.

Une deuxi`eme methode, plus rapide, consiste `a utiliser le lagrangien. En eet,

S
dS0
L
= mq =
=
= 2mE
q
q
dq

ce qui fournit immediatement q(t) = 2E
ethode est preferable, car elle evite dintegrer
m t + q0 . Cette m
p=

dS0
dq

S
pour obtenir lexpression de S0 , puis de calculer E
pour avoir Q(q) et enn de linverser. Mais elle
necessite la connaissance du lagrangien, ce qui peut ne pas etre le cas.

La chute libre 1D
Soit une particule de masse m soumise `a un potentiel gravitationnel constant. Lhamiltonien secrit
H(q, p) =

p2
+ mgq
2m

`
CHAPTER 4. SYSTEMES
HAMILTONIENS

56
et fournit donc lequation de Hamilton-Jacobi (HJ)
1
2m

S
q

2
+ mgq +

S
=0
t

Lhamiltonien ne dependant pas explicitement du temps, on a H = E conserve et la fonction principale

secrit S(q, t) = S0 (q) Et. Lequation (HJ) devient

dS0
= 2m(E mgq)
dq
On a deux methodes possibles de resolution pour obtenir q(t).
(a) Dans la premi`ere, on int`egre cette equation et on obtient
2
S0 (q) =
3g

2
(E mgq)3/2
m

La fonction generatrice S etant maintenant connue, on choisi P = E (par ex), et la nouvelle coordonnee est
obtenue en calculant

1 2
S
=
(E mgq)1/2 t
Q=
E
g m
ce qui, apr`es inversion, donne
q=

g
E
(t + Q)2
mg
2

(b) Dans la deuxi`eme methode, on utilise le Lagrangien

p

=
=

L
= mq
q

S
dS0
=
= 2m(E mgq)
q
dq

ce qui nous fournit une ODE du premier ordre qui sint`egre facilement et redonne le meme resultat.
Loscillateur harmonique 3D
Soit lhamiltonien suivant


3  2


pi
x2i
H(q, p) =
+ ki
2m
2
i=1

correspondant `a un oscillateur harmonique de masse m et de constante de raideur ki dans la direction xi .

Ce syst`eme etant ferme, lenergie se conserve et la fonction generatrice secrit
S(x1 , x2 , x3 , t) = S0 (x1 , x2 , x3 ) Et
Lequation de Hamilton-Jacobi est separable, en eet elle secrit


2

3


S
x2i
1
=E
+ ki
2m xi
2
i=1
On cherche donc une solution de la forme S0 =

3
i=1

Si (xi ). On obtient alors

F1 (x1 ) + F2 (x2 ) + F3 (x3 ) = E


4.1. LEQUATION
DE HAMILTON-JACOBI

57

qui doit etre satisfaite quels que soient les xi . Cela nest possible que si Fi (xi ) = i avec
i etant des invariants xes par les conditions initiales. On obtient alors



x2
dSi
= 2m i ki i
dxi
2

3
i=1

i = E, les

A ce stade, les deux methodes sorent `a nous. Soit on int`egre cette equation et on obtient ainsi la fonction
S0 , soit on repart du Lagrangien.

(a) La premi`ere donne en posant u =

ki
2i xi

Si =

= sin ,
i
(2 + sin 2)
2i

o`
u i2 = ki /m. En choisissant ensuite pour valeurs des nouveaux moments Pi = i , on obtient
Qi

=
=

S
Si
E
=

t
i
i
i

1
ki
arcsin
xi t
i
2i

cest `a dire, apr`es inversion,


xi =

2i
sin i (t + Qi )
ki

(b) La deuxi`eme methode est beaucoup plus directe. On obtient en eet

pi

L
= mx i
x i

S
dSi
=
=
xi
dxi



x2i
2m i ki
2

fournissant

du

= i dt
1 u2
cest `a dire = i t + i , ie. le meme resultat.
La force centrale
Nous avons vu que le mouvement dune particule soumise `a une force centrale seectue dans un plan et est
decrite par le lagrangien (en coordonnees spheriques)
L=

m 2
(r + r2 2 ) V (r)
2

correspondant `a un hamiltonien
H(r, , pr , p ) =

1 2 p 2
(p + 2 ) + V (r)
2m r
r

independant du temps et H = E. La fonction generatrice de Hamilton-Jacobi secrit alors

S(r, , t) = S0 (r, ) Et
et on cherche des solutions de la forme S0 (r, ) = Sr (r) + S (). La variable etant cyclique, on a
immediatement
S () = p

`
CHAPTER 4. SYSTEMES
HAMILTONIENS

58

et S(r, , t) = Sr (r) + p Et. Il ne reste donc plus qu`a resoudre lequation de Hamilton-Jacobi, mise sous
la forme

2
p2
1
dWr
+
+ V (r) = E
2m
dr
2mr2
cest `a dire
Sr =

 
p 2
2m(E V ) 2 dr
r

(a) En choisissant P = p et Pr = E, on obtient

Q

=
=

Qr

=
=

S
Sr
=
+
p
p

p dr


r2 2m(E V )
Sr
S
=
t
E
E
mdr

2m(E V )

p2
r2

p2
r2

qui, par inversion, nous fourniront la trajectoire r() et la loi horaire r(t) (Q et Qr sont des invariants).
(b) La methode du lagrangien fournit

p 2
S
dSr
L
= mr =
=
= 2m(E V ) 2
pr =
r
r
dr
r
qui sint`egre par rapport au temps et donnera t(r) puis r(t) apr`es inversion. De meme, on a
p =
et donc un lien (r):

4.1.6

d
dr

r .

L
= mr2

On obtient ainsi une integrale identique `a celle pour Q .

Le principe de Maupertuis

Nous avons vu que laction hamiltonienne, denie par

 t

S(q, t; P ) =
(
pi qi H)dt
t1

permettait de rendre compte de lequation de Hamilton-Jacobi par un principe variationnel. Cest egalement
la fonction generatrice de la transformation canonique qui produit des variables toutes cycliques. Pour un
syst`eme ferme (H ne dependant pas explicitement du temps), H est un invariant, de telle sorte que
 t

pi dqi H(t t1 )
S(q, t; P ) =
t1

Par ailleurs, laction secrit aussi dans le cas conservatif

S(q, t; P ) = S0 (q; P ) Et
A une constante sans importance pr`es, ces deux expressions sont egales et on obtient comme expression
generale de la fonction caracteristique de Hamilton
 q

S0 (q; P ) =
pi dqi
(4.14)
q
1


4.1. LEQUATION
DE HAMILTON-JACOBI

59

appelee plus couramment action reduite. Le param`etre temps a evidemment disparu, lintegration se
faisant le long dune trajectoire denie par les points
q1 et
q dans lespace des congurations.
Principe de Maupertuis: La trajectoire dun syst`eme conservatif est determinee par lextremisation
de laction reduite S0 .
En pratique, on peut calculer la trajectoire suivie par un syst`eme dynamique conservatif en utilisant
laction reduite. Il sut en eet dexprimer les impulsions en fonction des coordonnees (pi = L
qi ) puis
dutiliser les techniques de calcul variationnel pour lexpression S0 . Prenons un syst`eme soumis `a un potentiel
V (q). Dans ce cas, H = T + V = E et on a



L





pi dqi =
dqi =
mij qj dqi =
mij qj qi dt
qi
i
i
i
j
i,j
=

2(E V )dt

Exprimer lintegrale de S0 en fonction du temps nest pas habile car E est constant tandis que V = V (q). Il
vaut mieux exprimer dt en fonction des dqi . Utilisant la conservation de lenergie, on obtient

mij dqi dqj
dt2 =
2(E V )
i,j
do`
u


q

2(E V )

S0 (

q) =
q1



mij dqi dqj =

q


2(E V )ds

(4.15)

q1


i,j

o`
u le temps a disparu, remplace par lenergie E du syst`eme et o`
u lon a pose


ds2 =
mij dqi dqj
i,j

et dont linterpretation est le carre de lelement de distance innitesimale entre deux points de lespace des
u dl est
congurations (variete dierentiable). Ainsi, pour un syst`eme `a 1 degre de liberte ds2 = mdl2 o`
lelement innitesimal de longueur le long de la trajectoire.
Remarque: Il faut comprendre que lintegrale est denie sur lespace des congurations, cest `a dire
q (t1 ) vers un point

q (t). Laction reduite dans un cas `a n degres de liberte est
dun point

q1 =

 t

 t
 t
 q
 q
S0 (

q) =
pi qi dt =
p1 q1 dt + . . . +
pn qn dt =
p1 dq1 + . . . +
pn dqn
t1

4.1.7

t1

t1

q1


q
1

M
ecanique ondulatoire de Louis de Brooglie

Les exemples precedents montrent que, bien quorant en elle-meme une methode de resolution des probl`emes
mecaniques, lutilisation de lequation de Hamilton-Jacobi nest pas des plus pratiques. En general, il vaut
mieux utiliser les equations de Lagrange ou les equations canoniques de Hamilton.
Concr`etement, on peut neanmoins utiliser cette methode pour construire une transformation canonique
de fonction generatrice G(q, t; P ). En eet, il sut de la choisir de telle facon que
G
t
G
qk

= H(q, p, t)
= pk

G
avec Pk et P
constants.
k
En fait, limportance de cette equation reside surtout dans la comprehension profonde quelle a apporte
de la mecanique et dans le role quelle a joue dans lelaboration de la mecanique quantique.

`
CHAPTER 4. SYSTEMES
HAMILTONIENS

60
Analogie avec la construction de Huygens

Nous avons vu que pour un syst`eme conservatif, laction secrit

S(q, t) = S0 (q) Et
o`
u laction reduite S0 (q) est independante du temps. En consequence, les surfaces daction reduite S0 (q)
constante ont des positions xes dans lespace des congurations. Inversement, les surfaces degale action S
obeissent `a lequation suivante
S0 (q) = S(q, t) + Et
Ainsi, prenons la surface S(q, t) = a o`
u a est une constante. A t = 0 cette surface coincide avec la surface
S0 (q) = a, denissant ainsi un lieu dans lespace des congurations. Un instant dt plus tard, la surface
S(q, t + dt) = a coincide avec une autre surface, celle denie par S0 (q) = a + Edt. Autrement dit, les surfaces
degale action se deplacent dans lespace des congurations. Ces surfaces peuvent etre considerees comme
des fronts dondes se propageant dans lespace des congurations.
Par ailleurs, nous avons dej`a rencontre une propriete curieuse de lequation de Hamilton-Jacobi que nous
navions pas explique. En eet, elle provient bien dun principe variationnel mais avec laction hamiltonienne
S(q, t): le point darrivee est laisse libre. En fait, ce point darrivee depend des conditions initiales, en
S
. Ceci peut secrire vectoriellement
particulier des valeurs des pi = q
i

p =

S
= S0

q

(4.16)

Or, rien na ete impose aux pi au niveau du point de depart. Tout se passerait donc comme si le point
darrivee (autrement dit la trajectoire reelle suivie par le syst`eme) dependait de la direction initiale choisie.
Et cette direction est orthogonale `a la surface S(q, t) = a.
Ces deux elements, propagation de surfaces dondes et trajectoires orthogonales, ressemblent `a sy
meprendre `a la construction de Huygens pour loptique geometrique. Quelles seraient alors les proprietes
des rayons lumineux que lon pourrait rapprocher des proprietes des syst`emes mecaniques?
Les surfaces dondes lumineuses sont decrites par une phase constante (
r, t) =
k
r t,
k etant le
vecteur donde (k = 2/) et sa pulsation.
Dans un milieu dindice de refraction n(
r) inhomog`ene, la vitesse de phase u = nc des ondes lumineuses
depend de la position, c etant la vitesse de la lumi`ere dans le vide. En consequence les surfaces dondes
sont deformees lorsquelles se propagent dans un milieu inhomog`ene.
En optique geometrique, la trajectoire des rayons lumineux est orthogonale aux surfaces donde. Elle
obeit au principe de Fermat


cdt = nds = 0
(4.17)
etant une distance innitesimale dans la direction
o`
u u = ds
dt est la vitesse de phase de londe, ds
normale aux surfaces dondes.
Pour suivre lanalogie jusquau bout, regardons ce qui se passe pour une particule de masse m dans
un potentiel V (q). Nous avons vu que les surfaces S(q, t) = Constante seraient lequivalent des surfaces
donde. Cela signie une analogie entre action et phase. On y reviendra plus loin, pour linstant abordons
des questions simples:
A quelle vitesse de phase u se d
eplace la surface S = Cste?
La variation dS0 due `a un deplacement ds dans la direction orthogonale `a cette surface est donnee par
son gradient
dS0 = S0 M
M
= |S0 |ds


4.1. LEQUATION
DE HAMILTON-JACOBI

61

Or, le gradient est fourni par lequation de Hamilton-Jacobi

1
(S0 )2 + V = E
2m

(4.18)

tandis que dS0 = Edt. On obtient donc une vitesse de deplacement de la surface
u=

ds
E
E
=
=
dt
|S0 |
2m(E V )

(4.19)

Puisque cette vitesse depend de V (q) donc de la position, cela signie que la surface va se deformer. En
mecanique, lespace des congurations (milieu de propagation de londe) est inhomog`ene. . .

Cette
vitesse
ne
correspond
pas
`
a
la
vitesse
v
de
la
particule.
En
eet,
celle-ci
serait
v
=
2T /m =

2(E V )/m, ce qui donne une relation
E
(4.20)
uv =
m
analogue `a celle reliant vitesse de phase u et vitesse de groupe v pour une onde. Si on pousse lanalogie `a
terme, cest `a dire si on exige que la mecanique se comporte comme lelectromagnetisme, on obtient
E = mc2

(4.21)

Quel est lindice de r

efraction du milieu m
ecanique?
Nous venons de voir que la vitesse de phase u de la surface donde S = Cste depend de la position.
Cela implique un indice de refraction inhomog`ene. En optique geometrique, une onde se propageant dans
un milieu dindice de refraction n est caracterisee par une longueur de chemin optique L veriant lequation
iconale (de icone, du grec image)
(4.22)
(L)2 = n2
Ici, lequation maitresse est evidemment lequation de Hamilton-Jacobi, qui peut se mettre sous la forme
(S0 )2 = 2m(E V )

(4.23)

Puisquelle poss`ede la meme structure formelle, on peut associer

L
n

S0

2m(E V )

(4.24)
(4.25)

Du coup, on comprend enn lanalogie entre le principe de Fermat et le principe de Hamilton. Ils sont
en eet equivalents dans le cas qui nous interesse ici puisque ce dernier se ram`ene au principe de Maupertuis
 

2m(E V )ds = 0
(4.26)
nds = 0 S0 =
Quelle est la fr
equence de londe associ
ee `
a la particule?
Pour que lanalogie soit compl`ete il nous faut trouver lexpression de la frequence de londe associee `a
une particule. Les fronts donde etant denis par S = S0 Et = Cste correspondent `a des phases constantes.
On peut donc faire lanalogie
=
k
r t = 2(

L(r)
t) S = S0 (q) Et
0

o`
u 0 est la longueur donde dans le vide. Cela signie, puisque S0 joue le role du chemin optique `a une
constante multiplicative pr`es, que lenergie mecanique E doit etre proportionnelle `a la frequence de londe
associee. On peut donc poser
E = h
(4.27)

`
CHAPTER 4. SYSTEMES
HAMILTONIENS

62

o`
u la constante h est appelee maintenant la constante de Planck2. La longueur donde etant = u/ =
(E/mv)(h/E), on obtient
h
p = = hk
(4.28)

puisque p = mv pour une particule. Cette derni`ere relation est connue sous le nom de relation de De
Brooglie.
Vers une m
ecanique ondulatoire ?
Nous avons vu, grace `a lequation de Hamilton-Jacobi, lanalogie profonde de la mecanique avec loptique
geometrique. Mais il se trouve que loptique geometrique nest quune approximation des ondes electromagnetiques,
decrites par les equations de Maxwell3. En particulier, tous les phenom`enes ondulatoires, tels que la diraction ou les interferences, ne peuvent etre decrits dans le cadre de loptique geometrique. Cela impliquerait-il
alors que la mecanique classique ne soit quune approximation (du type optique geometrique) dune realite
plus vaste (du type ondulatoire)?
Lequation donde generale, issue des equations de Maxwell, secrit


n2 d2
2 2 2 = 0
c dt

(4.29)

et devient pour des ondes planes de la forme = 0 ei(krt)



2
2 + 2 = 0
u
Si lon suppose (comme nous lavons dej`a fait) que lequation ondulatoire (donc fondamentale) de la mecanique
est equivalente `a celle de lelectromagnetisme, alors on obtient quune fonction en mecanique devrait verier


2
2 + 2 = 0
u
p2
2 + 2 = 0
h
h2 2
 + (E V ) = 0
2m


h2 2
 + V = E

2m
H = E
2 En

consequence, laction normalisee S/

h o`
u
h = h/2 a la dimension dune phase.
les milieux inhomog`enes, une onde plane nest pas une solution car le front donde se deforme. En posant k0 = 2/0
o`
u 0 est la longueur donde dans le vide, on recherche des solutions de la forme
3 Dans

=
0 eA(r) eik0 (L(r)ct)

pour des ondes se propageant dans une direction donnee, A(r) etant une mesure de lamplitude de londe (constante pour une
onde plane) et L(r) du chemin optique (egal `
a nr pour un indice n constant). Lequation donde se ram`ene alors `
a un syst`eme
de deux equations

2 A + (A)2 + k02 n2 (L)2

L + 2A L
2

Soit l lechelle spatiale de variation de lamplitude A, donc celle de lindice n (puisque il en est la cause). La deuxi`eme equation
montre que le chemin optique varie alors sur la meme echelle. Lapproximation de loptique geometrique consiste `
a supposer
que l  0 . Cela signie alors que le 3`eme terme de la premi`ere equation est dominant, ce qui implique
(L)2 = n2
et qui est lequation iconale. Lanalogue en mecanique serait une longueur donde beaucoup plus petite que lechelle de variation
du potentiel.

4.2. VARIABLES CANONIQUES ANGLES-ACTIONS

63

o`
u lon reconnait la version stationnaire de lequation de Schrodinger de la mecanique quantique. . .
Lequivalence de lequation de Hamilton-Jacobi et de lequation iconale fut etablie par Hamilton luimeme en 1834. Mais il ny avait `a cette epoque aucune preuve experimentale de lexistence de phenom`enes
ondulatoires en rapport avec la mati`ere. Lequation donde correspondante fut decouverte par De Broglie et
Schrodinger presque un si`ecle plus tard, en 1926.
Le fait que lon ne trouve, par ces analogies, que la forme stationnaire de la mecanique quantique est assez
comprehensible. Tout se passe en eet comme si nous ne traitions, en mecanique classique, que de grandeurs
moyennees sur le temps et donc stationnaires du point de vue de la mecanique quantique. Il sav`ere que la
mecanique classique est eectivement une approximation de la mecanique quantique, obtenue dans la limite
h 0 (h = (E/c)). Remarquer cependant que, dans notre approche classique, il ny a aucune necessite
dintroduire la moindre quantication. Autant lenergie que laction peuvent prendre des valeurs continues.
Tout ce qui a ete dit ici, cest que la mati`ere pourrait en fait obeir `a une equation ondulatoire et donc
posseder des proprites en rapport (interferences, diraction). La notion de quantication est une notion
supplementaire qui sera abordee en mecanique quantique. Il sut ici de savoir que cette quantication se
traduit (entre autres) par une discretisation des valeurs prises par laction S et lenergie E dune particule.

4.2

Variables canoniques angles-actions

Dans de nombreux cas de syst`emes integrables, le mouvement sera periodique ou quasi-periodique. Nous
allons voir quil existe dans ce cas l`a un nouveau jeu de variables particuli`erement adapte, les variables
canoniques. Ce sont egalement les variables privilegiees pour letude des mouvements perturbes.

4.2.1

Syst`
emes ferm
es p
eriodiques

Pour un syst`eme `a 1 degre de liberte, il existe deux types de mouvements periodiques possibles, dej`a rencontres avec le pendule:
la libration: q(t) et p(t) sont des fonctions periodiques de meme frequence. La trajectoire dans lespace
des phases decrit une courbe fermee. On rencontre ce type de mouvement autour dun minimum de
potentiel par exemple.
la rotation (ou circulation): q(t) nest pas periodique mais q + q0 laisse le syst`eme invariant, tandis que
p(t) reste bornee. La trajectoire dans lespace des phases est ouverte mais q0 -periodique.
Pour un syst`eme `a n degres de liberte, le point representatif du syst`eme va decrire une trajectoire compliquee dans lespace des phases, notre representation se limitant `a une projection dans un plan (qi , pi ) par
exemple.
D
enition 1: Le mouvement dun syst`eme `a n degres de liberte sera dit periodique si la projection de
la trajectoire sur chaque plan (qi , pi ) est periodique dans le sens deni pour un mouvement `a 1 degre de liberte.
D
enition 2: Si la periode Ti associee `a un degre de liberte qi est la meme pour tous les qi , le mouvement du syst`eme sera dit simplement periodique. Si ce nest pas le cas, le mouvement sera dit multiplement
periodique ou quasi-periodique.
La consequence essentielle ici, cest que de tels syst`emes ont un espace des phases borne. Par ailleurs, le
mouvement etant periodique, cela signie que le syst`eme est invariant par translation dans le temps et donc
que lhamiltonien H est une integrale premi`ere.
Dans la suite, on ne consid`ere que des syst`emes de ce type.

`
CHAPTER 4. SYSTEMES
HAMILTONIENS

64

4.2.2

Variables angulaires

Motivation
La methode de Hamilton-Jacobi (`a ne pas confondre avec lequation du meme nom), consiste `a rechercher
une bonne transformation canonique. Dapr`es le theor`eme de Arnold-Liouville, un syst`eme integrable `a n
degres de liberte poss`ede n invariants Ik independants. Il nous sut alors de choisir pour nouveaux moments
Pi = Pi (Ik ), cest `a dire des fonctions quelconques des invariants, pour avoir
P i
Q i
H

H

=0
Qi
H

=
= i (Pi )
Pi
= H(Pi )
=

o`
u i sont egalement des invariants, puisque toutes les Qi seront alors cycliques et que H ne depend pas
explicitement du temps. La dependance temporelle des Qi est alors excessivement simple
Qi (t) = i t + Qi0
La question qui se pose est alors la suivante: Comment choisir les nouveaux moments Pi (Ik ) (on suppose
evidemment connues les n integrales Ik ) ? Dit autrement, existe-t-il un choix plus naturel que dautres ?
Un espace des phases born
e
Puisque les integrales premi`eres restent constantes tout au long du mouvement, celui-ci a lieu dans un sousespace de lespace des phases, de dimension n, note Tn , et qui est lintersection des surfaces Ik (qi , pi ) = Cste.
Nous avons vu que les Qi varient lineairement avec le temps. Si ces nouvelles coordonnees sont du type
distance, alors elles peuvent prendre des valeurs qui croissent indeniment jusqu`a linni. Cest ce qui
se passerait si Tn etait un hyperplan. Mais nous ne traitons ici que le cas de syst`emes (quasi-)periodiques,
cest `a dire nexplorant quune partie nie (bornee) de lespace des phases. Cela implique donc que les Qi
sont du type angle, variant de 0 `a 2.
Topologiquement, cela implique que Tn est un hypertore de dimension n, plonge dans un espace de
dimension 2n. Un hypertore de dimension n est le produit tensoriel de n cercles. Soit un point M Tn ,
alors M est caracterise par n angles (1 , . . . , n ), o`
u i est langle sur le i-i`eme cercle4 .
Sans faire de la topologie, peut-on comprendre intuitivement pourquoi est-ce un tore et non une sph`ere?
Dans un probl`eme `a 1 dimension comme celui du pendule par exemple, H = E est un invariant et un
mouvement de libration du pendule fournit une ellipse. Nous avons vu quil etait possible de se ramener `a un
cercle de rayon R, directement relie `a la valeur de E. On peut alors changer de coordonnees (q, p) (, R)
et se placer en coordonnees polaires: le rayon R(E) est un invariant et la deuxi`eme coordonnee est un angle
.
Si le probl`eme est `a 2 dimensions, nous aurons 2 invariants, chacun assimilable `a un rayon, les 2
autres coordonnees assimilables `a des angles. La position dun point M (1 , 2 ) de lespace des phases
peut sapprehender de la mani`ere suivante. Sur le premier cercle de rayon R1 se trouve, situee en un angle
1 , lintersection du deuxi`eme cercle de rayon R2 . Le point M se trouve alors en un angle 2 sur ce deuxi`eme
cercle. La trajectoire decrite par M est ainsi une helice qui senroule autour du tore T2 .
La position dun point M sur une sph`ere serait egalement decrite par deux angles (latitude et longitude).
Mais il ny a quun seul invariant associe dans ce cas l`a,le rayon de la sph`ere, et pas deux.
4 Ce cercle correspond `
a un lacet irreductible sur le tore, cest `
a dire `
a un contour ferme non assimilable par deformation
continue `
a un point.

4.2. VARIABLES CANONIQUES ANGLES-ACTIONS

4.2.3

65

Variables dactions

Nous venons donc de montrer que pour les syst`emes periodiques les nouvelles variables Qi doivent etre des
angles i , variant entre 0 et 2. Cela contraint alors automatiquement le choix des nouvelles variables Pi
puisque
H
i (Pi ) =
Pi
doivent etre des pulsations.
Les nouvelles Pi doivent donc avoir la dimension dune energie divisee par une pulsation, cest `a dire la
dimension dun moment cinetique. Or, cest precisement la dimension de laction. Dans la suite, on notera
Ji ces variables dactions.
Cas dun syst`
eme `
a 1 degr
e de libert
e
Lespace des phases est de dimension 2 et lon cherche la transformation canonique telle que les nouvelles
variables verient
(t)

= t + 0

J(t)

= J

cest `a dire, soient respectivement un angle variant lineairement dans le temps et une action J constante.
Lhamiltonien dun tel syst`eme secrit toujours
H=

p2
+ V (q)
2m

Puisque nous ne nous interessons qu`a des syst`emes periodiques, cela signie que lespace des phases est
borne et donc, que le portrait de phase engendre par lensemble des contours

p(q) = 2m(E V (q))
dessine une aire A(E) pour une energie E donnee du syst`eme. Or, cette aire est precisement un invariant
canonique de Poincare, pour lequel on a
 
 
dqdp = =
dQdP
 A
 A
dqdp =
pdq
A

o`
u est un contour ferme de lespace des phases, oriente dans le sens gauche ou indirect (important)5 .
Utilisant ces proprietes et le fait que est un angle variant entre 0 et 2 tandis que J doit etre un invariant,
5 La

deuxi`eme propriete peut se demontrer comme cas particulier de la formule de Stokes:

 



 ( A) dS =


 dl
A

 = dxdy
est une surface delimitee par un coutour ferme oriente dans le sens direct. Prenons le cas 2D o`
u dS
k,
o`
u
 = dxi + dyj et A
 = (y; x; 0)/2. Alors A
 = (0; 0; 1) et on obtient
dl

 

 



 ( A) dS

 dxdy

 =
 dl
A

xdy ydx
2

Par ailleurs, sur un contour ferme, cest `

a dire partant dun point M pour y revenir, on a

d(xy)

[xy]M
M =0

`
CHAPTER 4. SYSTEMES
HAMILTONIENS

66
on obtient

A(E) =

pdq =

Jd = J

d = J2
0

cest `a dire lexpression suivante pour la variable action


1
J=
pdq
2

(4.30)

Cas dun syst`

eme s
eparable `
a n degr
es de libert
e
Nous avons vu que laction hamiltonienne dun syst`eme conservatif `a n degres de liberte peut secrire
S(q, t; P ) = S0 (q) Et
o`
u S0 (q) est laction reduite ou fonction caracteristique de Hamilton. Etant parfaitement separable, elle
verie




pi dqi =
Si (qi ; Ik )
S0 (q; P ) =
i

les n Ik etant nos n integrales premi`eres independantes. Or, chaque moment conjuge pi est obtenu par
pi =

Si
= pi (qi ; Ik )
qi

et donc ne depend que de la variable conjuguee qi (et des n invariants). Cela signie que chaque paire de
variables canoniquement conjuguees (qi , pi ) est independante des autres, laction reduite etant la somme de
n termes independants.
Dans chaque plan (qi , pi ), le mouvement est borne (syst`eme periodique) et laire Ai (Ik ) formee par la
projection du mouvement sur ce plan est conservee. Il est donc naturel dintroduire la grandeur denie par

1
Ji =
pi dqi
(4.31)
2 i
o`
u lintegration seectue sur un contour ferme de ce plan, parcouru dans le sens gauche. On aura ainsi
Ji = Ji (Ik ) uniquement, puisque la variable qi nintervient que comme variable dintegration. Cela signie
donc que Ji est un invariant du mouvement. Et comme chaque paire (qi , pi ) est independante des autres, les
n Ji forment egalement n invariants independants.
La transformation canonique compl`ete est alors
H
=0
i
H
= i (Ji )
Ji

Ji

ce qui est la transformation recherchee, les i (Ji ) etant bien de la dimension dune pulsation.

=
do`
u

xdy =

(xdy + ydx)

ydx et on obtient nalement

 


 dxdy =


ydx =

xdy

En general, on associe x `
a la variable q et y `
a p, le hamiltonien permettant dexprimer facilement p = p(q). Cela signie quil
faut utiliser lexpression du milieu. Le signe - disparait si lon decide dorienter le contour dans le sens gauche (cest `
a dire
indirect, contraire `
a celui donne par la r`egle du bonhomme dAmp`ere, dans le sens des aiguilles dune montre).

4.2. VARIABLES CANONIQUES ANGLES-ACTIONS

67

Cas g
en
eral (non d
emontr
e)
Le mouvement du syst`eme a lieu dans un hypertore Tn , sous-espace de lespace des phases de dimension n.
Sur ce tore, on admettra que lon peut y trouver n contours irreductibles (non reduisibles `a un point) i
dierents, cest `a dire ne pouvant etre transformes lun en lautre par une deformation continue.
Les variables actions, denies par


1
Ji =
pk dqk
(4.32)
2 i
k

forment ainsi n integrales premi`eres du mouvement (puisque fonctions uniquement des n Ik ), independantes.
Le fait quelles soient independantes nest pas evident, cest lie `a lindependance des n contours i . Cest
assure par la condition dinvolution invoquee dans le theor`eme de Arnold-Liouville.

4.2.4

Fonction g
en
eratrice des variables angles-actions

Le choix particulier des variables canoniques angles-actions correspond `a prendre pour fonction generatrice
laction hamiltonienne reduite


S0 (q; Pk ) =
pi (qk ; Pk )dqi
i

En eet, cest une fonction generatrice du type G2 (q, P, t) donnant les relations generales suivantes
pi

Qi

H

S0
qi
S0
Pi
H+

S0
t

Dans notre cas le syst`eme est conservatif, Qi i et Pi Ji . On a donc bien H
= H et les n Ji etant n
invariants independants, on peut exprimer pi = pi (qk , Jk ). Avec ces notations on obtient
pi

S0
qi
S0
Ji

Si on montre que, les i etant des angles, les Ji sont alors des variables daction, alors on aura montre que
S0 est bien la bonne fonction generatrice.
u Ti =
On suppose donc que i est un angle recevant un accroissement i = 2 entre t et t + Ti o`
2/i (Jk ) est lune des periodes du syst`eme. On a alors
i

S0
S0
S0
= i (t + Ti ) i (t) =
(t + Ti )
(t) =
= 2
Ji
Ji
Ji






 t

 Ti

t+Ti

pk dqk
pk dqk =
pk dqk =
pk dqk
=
Ji
Ji
Ji
0
0
0
i
k

(2Ji )
Ji

ce qui demontre que S0 est bien la fonction generatrice de notre transformation canonique.
Si lon veut exprimer les variables angles-actions en fonction des anciennes variables (qi , pi ), il sut de
remplacer les invariants par leur expression en fonction des (qi , pi ).

`
CHAPTER 4. SYSTEMES
HAMILTONIENS

68

Regardons le cas `a 1 degre de liberte. On remplace E par H(q, p) ce qui fournit

J

J(E) = J(H(q, p)) = J(q, p)

(J)t + 0 = (q, p)t + 0

Exprimer (q, p) en fonction des (, J) est plus laborieux et plusieurs methodes sont disponibles. On peut
tout dabord essayer dinverser les relations ci-dessus. Si ce nest pas facile, il sut de rechercher q = q(, J)
puisque nous avons dej`a exprime p = p(q, E) `a partir de lhamiltonien. Deux methodes sorent `a nous:
0
(1) Si lon connait S0 (q; E), on lexprime en fonction de S0 (q; J) et on peut ainsi calculer S
J . Cela nous
fournit une relation = (q; J) qui, apr`es inversion, nous donnera q = q(; J).
(2) Sans calculer la fonction generatrice, on peut utiliser les equations canoniques, cela donne
 t
 q
dq
0 = t =
dt =
0
q0 q
 q
dq
= (J)
H
q0

qui, par inversion, nous fournira q = q(; J).

4.2.5

Applications

M
ethode
La methode est assez simple, une fois que lon connait n invariants Ij satisfaisant les conditions du theor`eme
de Arnold-Liouville (independance, involution):
(1) On exprime les pi en fonction des qi et des invariants;
(2) On choisit n contours irreductibles i , puis on calcule les variables actions Ji ;
(3) On exprime H = H(Ji );
(4) On calcule ensuite les variables angles par
i (t) = i t + i0

avec

i =

H
Ji

(4.33)

On a donc ainsi acc`es aux n pulsations propres i dun syst`eme periodique `a n degres de liberte sans
meme avoir eu recours `a la resolution dequations dierentielles. Nous avons fait des inversions de formules
(etapes 1 et 3), des integrations (etape 2) et des derivations (etape 4). On appelle cela resoudre les equations
du mouvement par quadratures.
Exemple 1: Chute libre 1D
Soit une particule de masse m en chute libre dans un champ de gravitation g, rebondissant sur le sol de
facon elastique. Le hamiltonien secrit
p2
+ mgq
H=
2m
et le portrait de phase est decrit par les courbes

p(q; E) = 2m(E mgq)
La solution negative correspond `a la chute depuis une altitude q0 = E/mg, la solution positive au rebond
depuis q = 0 jusqu`a q0 , formant ainsi un contour ferme parcouru dans le sens gauche. La variable action
est alors

 q0 
1
2
J =
p(q)dq =
2m(E mgq)dq
2
2 0

2
2 3
E
=
3g m

4.2. VARIABLES CANONIQUES ANGLES-ACTIONS

69

u la pulsation est donnee par =

La variable angle est (t) = t + 0 , o`
H(J). Par inversion de lexpression de J, on obtient

H(J) =

9 2 mg 2
8

do`
u une pulsation
H
=
= g
J

H
J .

Il faut donc exprimer dabord

1/3
J 2/3

m
2E

Si lon veut (q, p) et J(q, p) il sut de remplacer E par H(q, p) dans leurs expressions. Inversement,
pour revenir aux variables (q, p), il sut de connaitre q(, J) puisquon a dej`a p(q; E) et E(J). Lequation
eecrire
canonique de Hamilton q = H
p , peut se r

0 =

dt =
0

q=0

dq
=
q

dq

H
q=0 p

du

= ( 1 u 1)
1u

avec u = mgq/E = q/q0 . On obtient ainsi la relation recherchee



2 
0
E(J)
1 1+
q(, J) =
mg

Plus usuellement, cette expression est equivalente `a q(t) = 12 gt2 + at + b, a et b devant etre determinees
par les conditions initiales (ici on a arbitrairement choisi q(t = 0) = 0).
Exemple 2: Oscillateur harmonique 1D
Le portrait de phase est decrit par les courbes dequation

p(q; E) =

2m(E

k 2
q )
2

cest `a dire une ellipse, de valeur maximale a = 2E

2mE sur p. Cette ellipse est parcourue
k sur q et b =
dans le sens gauche. La variable action est


ab
m
A(E)
1
=
=E
p(q)dq =
J=
2
2
2
k

k
On en deduit immediatement H(J) = J m
et une pulsation
H
=
=
J

k
m

Lexpression J(q, p) est immediate puisquil sut de remplacer E par H(q, p). On peut obtenir (q, p)
en procedant comme suit:





dq
dq
dq
du
m

=
m
=
arcsin(q
=
=
=
)
H
2
q
p
2E
1

u
p
q
= arcsin 
2
p
2 2
m2 + q

On a ainsi `a la fois (q, p) mais aussi q(; E) =

2E
k

sin(t + 0 ).

`
CHAPTER 4. SYSTEMES
HAMILTONIENS

70
Syst`
emes quasi-p
eriodiques

Prenons lexemple simple de loscillateur harmonique 3D. Dans ce type de probl`eme, on consid`ere generalement
que les trois directions sont decouplees et quon a donc aaire `a un syst`eme compl`etement separable. La
projection du mouvement dans chaque plan (qi , pi ) de lespace des phases est donnee par
p2i
ki
+ qi2
2m
2

Ii =

o`
u ki est la constante de raideur dans cette direction et Ii une constante du mouvement telle que E =
La variable action est alors

m
Ji = Ii
ki
tandis que la variable angulaire secrit

i (t) = i t + i0 = arcsin(qi

3

i=1 Ii .

ki
)
2Ii


avec une pulsation i = kmi . Chaque coordonnee qi est donc une fonction periodique de i et, dans ce cas
simple, egalement une fonction periodique du temps, de periode Ti = 2/i .
Considerons maintenant une fonction F quelconque de letat dynamique du syst`eme. Cette fonction
sexprimera en fonction des coordonnees generalisees, dune facon qui peut etre non triviale. Ce sera une
fonction periodique des i , puisqu`a chaque periode Ti la coordonnee qi reprend sa valeur initiale (dans le cas
dune libration, sinon modulo qi0 dans le cas dune rotation). Mais il nest pas dit que ce soit une fonction
periodique du temps. Peut-on en eet encore trouver un temps T tel que F (t+ T ) = F (t)? Si loscillateur est
isotrope, les trois periodes Ti sont identiques et on assiste `a des trajectoires fermees. Par exemple `a 2D, les
trajectoires dessineraient des courbes de Lissajou fermees. Mais si loscillateur est anisotrope les trajectoires
ne sont plus des courbes fermees (`a moins peut-etre dattendre tr`es longtemps. . . ).
Reprenons le cas dun syst`eme `a n degres de liberte pour lequel on a trouve n variables angulaires i .
Soit F (q, p) une fonction quelconque detat du syst`eme. Cest une fonction periodique des i de periode 2
pour chacune delles (i = 2). On peut donc faire un developpement multiple en series de Fourier6




F (1 , . . . , n ) =
...
Ck1 ...kn ei(k1 1 +...+kn n )
k1

kn

o`
u les Ck1 ...kn sont des nombres complexes. Exprimee en fonction du temps, cette expression devient



F (t) =
...
Ak1 ...kn eit(k1 1 +...+kn n )
k1
6 Une

kn

fonction f (t) T -periodique peut se decomposer en serie de Fourier de la forme

a0
(ak cos kt + bk sin kt)
+
2

f (t)

k=1

ck eikt

k=

o`
u les coecients verient
ak

bk

ck

2
T
2
T
1
T

f (t) cos ktdt

0

f (t) sin ktdt

0

f (t)eikt dt
0

4.2. VARIABLES CANONIQUES ANGLES-ACTIONS

71

Le terme en argument dans lexponentielle peut secrire symboliquement

k o`

u
k est un vecteur dont
les n composantes sont des nombres entiers relatifs independants. On ne verra apparaitre de la periodicite
que si toutes ou partie des pulsations i sont commensurables entre elles.
D
enition 1: Deux pulsations sont commensurables si leur rapport est egal `a un nombre rationnel, cest
n
1
`a dire si
u n et m sont des nombres entiers7 .
2 = m o`
D
enition 2: Un syst`eme `a n degres de liberte est dit partiellement degenere si 2 ou plus de ses pulsations propres sont commensurables. Il est dit compl`etement degenere si elles sont toutes commensurables.
D
enition 3: Un syst`eme quasi-periodique est un syst`eme partiellement degenere.
Dapr`es ces denitions, un syst`eme compl`etement degenere est un syst`eme pour lequel on a n1 relations
du type


k = 0
puisquon pourra alors exprimer les n 1 i en fonction dune seule pulsation. Si lon na que m < n 1
relations de ce type, alors on a aaire `a un syst`eme partiellement degenere ou m fois degenere.
Prenons le cas de deux pulsations commensurables 1 et 2 . On a alors legalite suivante
t =

2
2
=
k1 1
k2 2

cest `a dire k1 k2 t = k2 T1 = k1 T2

La vraie periode T commune `a ces deux mouvements est obtenue en divisant k1 k2 t par le PPCM de k1 et
k2 .

7 Cette

condition est equivalente `

a k1 1 + k2 2 = 0, o`
u k1 et k2 sont des nombres entiers relatifs.

72

`
CHAPTER 4. SYSTEMES
HAMILTONIENS

Chapter 5

Description lagrangienne des milieux

continus
5.1
5.1.1

Exemple dun passage `

a la limite continue
Corde
elastique 1D

Considerons une corde continue de longueur l tendue `a ses deux extremites. Nous connaissons sa forme `a
lequilibre mais que devient-elle lorsquon tire leg`erement sur la corde ? Pour traiter ce probl`eme, utilisons
les methodes vues precedemment.
On va ainsi modeliser la corde comme une succession de N masses m reliees entre elles par des ressorts de
raideur k. Cette modelisation est legitime puisque nous nous placons au voisinage de la position dequilibre:
tout potentiel dinteraction est en eet assimilable `a un potentiel harmonique (voir chapitre I). Soit i lindice
decrivant une masse particuli`ere, i = 0 et i = N + 1 correspondant aux deux extremites. La coordonnee
generalisee, notee i (t), decrit le deplacement de la masse i par rapport `a sa position dequilibre. Enn,
notons que, `a lequilibre, la distance entre deux masses est a = l/(N + 1).
Ceci etant pose, le lagrangien de la corde secrit

N 


1
1
m i2 k(i+1 i )2
L=
2
2
i=1
On obtient ainsi N equations de Lagrange
d L
dt k

= m
k
=

L
= k(k k1 k+1 + k )
k
k
[(k+1 k ) (k k1 )]
m

Maintenant, il est evident que si lon pouvait seviter de resoudre N equations, N etant grand, on le
ferait. Or justement, N est amene `a tendre vers linni dans notre mod`ele dun syst`eme mecanique continu.
Que se passe-t-il dailleurs lorsque N ? On a a 0 et m/a , etant la masse lineique de la corde.
Tout se passe donc comme si lon seloignait de la corde et que lon ne pouvait plus separer deux masses
lune de lautre. Continuer `a utiliser lindice k na donc plus de sens: il faut le remplacer par la position x de
la masse m, x etant maintenant une variable continue. Ainsi, on a k (t) = (xk , t) (x, t) lorsque a 0.
Par ailleurs,

k+1 k
(xk+1 , t) (xk , t)

lim
= lim
=
a0
a0
a
xk+1 xk
x k+1/2
73

74

CHAPTER 5. DESCRIPTION LAGRANGIENNE DES MILIEUX CONTINUS

k k1
lim
a0
a

(xk , t) (xk1 , t)
lim
=
a0
xk xk1
x k1/2

En reprenant lequation de Lagrange pour la masse k et en faisant a 0, on obtient





a2 k
2 (x)

lim k =
= lim

a0
a0 ma
t2
x k+1/2
x k1/2
ak 2 (xk )
a0
x2

lim

Comment varie le terme ak lorsque a 0? La description du potentiel cree par un ressort provient de
la loi experimentale de Hooke. Celle-ci stipule que lallongement dun ressort est proportionnel `
a la force
appliquee par unite de longueur. Mathematiquement, si est un allongement (en m`etres) et F la force
appliquee (en Newtons), cela secrit
F
= E
lim
a0 a
o`
u E, appele le module dYoung, caracterise la durete du materiau. Si lon revient `a la corde, on obtient
lima0 ak = E, ce qui fournit lequation suivante pour la corde 1D,

 2

E 2

(x, t) = 0
(5.1)
t2
x2
Cette equation montre que la corde 1D va subir des deplacements longitudinaux
oscillants, cest `a dire

quelle va etre parcourue par une onde longitudinale de vitesse de phase c = E/: plus la corde est rigide
(E eleve) et/ou plus elle est leg`ere ( faible) et plus c sera elevee.
Nous voyons apparaitre sur cet exemple un trait caracteristique du passage `a la limite continue. Nous
sommes en eet passes dun nombre N eleve (en fait N ) dequations aux derivees ordinaires (EDO)
`a une seule equation aux derivees partielles (EDP). Mais celle-ci contient la meme information que linnite
dEDO ! Cest d
u au fait que la solution de (5.1) est le deplacement (x, t), deni en tout x, variable continue.

5.1.2

Retour au lagrangien

Ce que lon voudrait avoir, cest une nouvelle facon decrire le lagrangien dun syst`eme continu de telle sorte
que lon puisse trouver une equation analogue `a (5.1) sans avoir `a refaire, `a chaque fois, explicitement le
passage `a la limite. Reprenons le cas de la corde pour laquelle le lagrangien secrit
L=

N 


1
i=1

1
m 2 k(i+1 i )2
2 i
2

Au vu de la methode employee plus haut, on est tentes de mettre le lagrangien du syst`eme continu sous la
forme suivante


2 


i+1 i
1m 2 1
ak
a
L = lim
a0
2 a i
2
a
i


 2
 2


E
=
dx

(5.2)
2 t
2 x

cest `a dire L = dx L(t , x )1 . Si cette conjecture est correcte, alors nous devrions pouvoir retrouver
lequation (5.1) directement, `a partir des equations dEuler-Lagrange. Mais ces derni`eres ne sappliquent
qu`a des syst`emes discrets. Il nous faut donc trouver leur equivalent pour des syst`emes continus.
1 Pour

des raisons de commodite decriture, on note t =

similairement, xx =

,
x2

xy =

.
xy

et x =

.
x

Des derivees dordre superieur seraient notees

5.2. FORMULATION LAGRANGIENNE DES MILIEUX CONTINUS

75

Dans lexpression du lagrangien L de la corde interviennent les derivees de (x, t). Il faut bien realiser
quici x nest pas une coordonnee associee `a la corde elle-meme: x decrit une position sur la corde relative
`a notre mode de mesure (cest un indice dans le cas discret), au meme titre que le temps t est relatif `a
notre horloge. Ainsi, x et t sont deux param`etres continus relatifs `a notre description de lespace-temps
et non des variables decrivant la corde. La variable decrivant la corde, autrement dit la vraie coordonnee
generalisee, reste (x, t). Par denition, on appelle champ toute grandeur denie sur lespace-temps. Ainsi,
(x, t) est le champ (scalaire) de deplacement de la corde 1D. Pour cette raison, les equations de Lagrange
continues sont appelees les equations du champ.
Si lon avait permis des deplacements de la corde dans les trois directions despace, on aurait tout
simplement un champ (x, y, z, t) ainsi quune generalisation de lequation (5.1), cest `a dire


3
E
2
2

(5.3)
(xi , t) = 0
t2
i=1 x2i
o`
u xi vaut pour x, y, z.

5.2
5.2.1

Formulation lagrangienne des milieux continus

Conjectures

On suppose que tout syst`eme continu, cest `a dire ayant une innite de degres de liberte, peut etre decrit
par un lagrangien se mettant sous la forme
 x2  y 2  z2
L=
L dxdydz
(5.4)
x1

y1

z1

o`
u L est appelee la densite de lagrangien. Le volume dintegration est arbitraire: il doit contenir le syst`eme
etudie mais peut etre etendu `a linni (il sut de faire tendre L vers zero en dehors du syst`eme).
Il faut (de nouveau) bien realiser que les variables x, y et z ne sont que des param`etres permettant
didentier une region de lespace et non des variables decrivant le syst`eme lui-meme. Celui-ci est donc vu
comme etant plonge dans un continuum (lespace-temps).
Le syst`eme physique est suppose enti`erement decrit par un nombre ni n de champs i (
r, t) (i = 1, . . . , n).
On suppose que la densite de lagrangien peut secrire comme



,
,
,
, x, y, z, t
(5.5)
L = L ,
x y z t
et non en fonction des derivees dordre superieur de . Cette hypoth`ese est tr`es importante car elle a de
profondes consequences. Concr`etement, elle impose directement la forme des equations du champ (voir cidessous), philosophiquement elle signie que les champs interagissent seulement localement: en eet, pour
calculer une derivee dordre deux, par exemple, il faut aller chercher deux points de part et dautre du point
o`
u lon veut la calculer alors quil nen faut quun seul pour une derivee dordre un. Cette hypoth`ese ne peut
se justier par elle-meme: il faut verier experimentalement que les equations ainsi obtenues fournissent bien
une description correcte des phenom`enes.
Le fait que les coordonnees x, y, z et t ne jouent que le role dindexation de lespace-temps et non celui
de description du syst`eme lui-meme introduit une subtilite. Soit x une coordonnee quelconque de lespacetemps, si par exemple L = L(, x ) alors on a
L
L
dL
=
+
dx
x
x

(5.6)

Autrement dit, il faut faire attention `a la notion de derivee droite ou de derivee partielle. La variation de
L entre x et x + dx met en eet en jeu la variation liee uniquement au changement de position (terme

L
egalement celle introduite par la dependance de (terme x
).
x ) mais

76

5.2.2

CHAPTER 5. DESCRIPTION LAGRANGIENNE DES MILIEUX CONTINUS

Equations de Lagrange du champ

On part du principe de Hamilton, principe fondamental de la physique. Il stipule que les champs realises
dans la nature sont ceux qui rendent extremale laction
 t2  x2  y2  z2
 t2
L dt =
L dxdydzdt
(5.7)
S=
t1

t1

x1

y1

z1

cest `a dire tels que S = 0. Ici, S est la variation fonctionnelle de laction naectant pas les bornes (ni
le temps, ni le volume) et obtenue en faisant une variation i de tous les champs independants i telle
que i = 0 aux bornes dintegration2. An dalleger lecriture, on fait ici la demonstration pour i = 1, la
generalisation `a n champs independants etant ensuite immediate.
La variation fonctionnelle de S secrit donc
 t2  x2  y2  z2
 
S =
L dxdydzdt =
L dtd3 v
t1

x1

y1

z1



 
o`
u lon a dt d3 v L = d3 v dt L puisque les coordonnees x, y, z et t sont independantes. Par
ailleurs, on a
3


L
L
L
+
t +
x
L =

(t )
(x )
=1
avec x1 = x, x2 = y, x3 = z. Regardons ce que donne lintegration de chacun de ces termes sur le
(quadri)volume despace-temps. Ainsi, le deuxi`eme terme secrit







L
L
L ()
3
3
3
t =

= d v dt
d v dt
dt d v
(t )
(t ) t
(t ) t
t2 


 !

L
d L
3

dt
=
d v

(t )
dt (t )
t1
o`
u lon a integre par parties `a la deuxi`eme ligne. Le premier terme de droite est nul car = 0 aux bornes.
Remarquer lutilisation de la derivee droite d/dt: lors de lintegration par parties, on prend en compte toutes
les variations (explicites et implicites) lors dune modication t t + dt. Faisant le meme travail sur la
variable x, on obtient






L
L ()
x =
dt dy dz dx
dt d3 v
(x )
(tx ) x
x 2 


 !



d L
L

dx
=
dt dy dz

(x )
dx (x )
x1
o`
u le premier terme du membre de droite est nul et on utilise la derivee droite d/dx pour les memes raisons
que precedemment. Le calcul est identique pour les variables y et z. Rassemblant ces resultats on obtient
!


3


d L
d
L
L
3

= 0
S =
dt d v

dt (t ) =1 dx (x )
Ceci devant etre verie quelle que soit la variation , cela implique

d
d L
L
L
+
=
dt (t ) =1 dx (x )

(5.8)

2 Dans le cas discret, on dit que tous les chemins possibles doivent partir de t et arriver en t , autrement dit q = 0 en ces
1
2
i
deux bornes. Dans le cas continu, on dit que tous les champs possibles doivent se comporter de la meme mani`ere aux bornes.

5.2. FORMULATION LAGRANGIENNE DES MILIEUX CONTINUS

77

Cette equation est lequation de Lagrange pour un syst`eme decrit par un champ continu . Il est remarquable que le role du temps apparaisse exactement sous la meme forme que celle pour une variable spatiale
quelconque. La raison profonde reside tout betement dans le fait que les coordonnees de temps et despace
ne font que decrire le continuum espace-temps sur lequel evolue le syst`eme physique. On peut donc ecrire
cette equation sous forme plus compacte en etendant le nombre de coordonnees `a 4 (espace + temps). Pour
un nombre n de champs independants i (
r, t), on obtient alors les n equations du champ
3


L
d
L
=
dx
(

i
i
=0

(5.9)

et x
: cest la
o`
u lon a utilise cette fois-ci la notation = x . Noter quil ny a pas dambiguite entre dx
meme chose puique = (x, y, z, t) et que les 4 coordonnees despace-temps sont independantes.

Exemple: Reprenons la corde elastique 1D et supposons quelle soit decrite par la densite de lagrangien
suivante
 2
 2
E

L=
2 t
2 x
Les diverses derivees secrivent alors
L

d L
dt (t )
d L
dx (x )

= 0
=
=

d
(t ) = tt
dt
d
(Ex ) = Exx
dx

ce qui fournit lequation tt Exx = 0 pour le champ qui est bien identique `a lequation (5.1). Cette
equivalence (ainsi que pour bien dautres exemples) justie, dune part, nos conjectures initiales et, dautre
part, notre conance dans lutilisation de la densite de lagrangien choisie pour decrire la corde.
Densit
es de lagrangiens
equivalentes: Nous avons vu que dans le cas discret deux lagrangiens L et
L
decrivent en fait le meme syst`eme sils ne di`erent que par la derivee totale par rapport au temps dune
fonction f dependant des coordonnees generalisees qi et du temps. Dans le cas continu, deux densites de
lagrangiens L
et L correspondent au meme syst`eme si
3

L = L+
f (i ,
r, t)
x

=0

(5.10)

o`
u f pour = 0, . . . , 3 sont quatre fonctions scalaires quelconques, denies sur lespace-temps et dependant
des champs (mais pas de leurs derivees). Le terme supplementaire est en fait la 4-divergence dun quadrivecteur,
que lon peut ecrire de facon parfaitement equivalente
L
= L +

f0
+ div F

t

(5.11)

o`
u lon a construit un vecteur F
tel que chaque composante verie Fk = fk (i ,
r, t) pour k = 1, 2, 3.
La demonstration sappuie sur le principe de Hamilton puisque lon doit avoir S
= S = 0. Ainsi,
   


S =
L
dxdydzdt




f0
+ dt d3 v div F

= S + d3 v dt
t



t
3
2
= 0 +
d v [f0 ]t1 +
dt d3 v div F

78

CHAPTER 5. DESCRIPTION LAGRANGIENNE DES MILIEUX CONTINUS

Le dernier terme de droite met en jeu une integrale sur le volume de la divergence de F
. Dapr`es le theor`eme
de Gauss-Ostrogradski, cette integrale est equivalente au ux de F
`a travers la surface A englobant le volume
dintegration (qui peut etre inni). On a ainsi
S


=

d v

 n

f0
i=1

t2
i


+

t1




F
dA

dt
A

Le premier terme est manifestement nul de meme que le second (moyennant la precaution supplementaire
dexiger que les fk tendent susamment rapidement vers zero en A). Remarquer que si les f dependaient
aussi des derivees des i , ceci ne serait plus vrai en toute generalite comme cest le cas ici.

5.3
5.3.1

Th
eorie classique des champs
Cadre g
en
eral

Jusqu`a present, notre demarche a ete dutiliser le formalisme lagrangien pour decrire un syst`eme mecanique
continu (ie. avec un nombre inni de degres de liberte). Mais on peut tout aussi bien utiliser la meme
formulation pour decrire nimporte quel objet physique pourvu quil soit decrit en termes dun ou plusieurs
champs. Par champ, on entend toute grandeur i (
r, t) denie sur lespace-temps. Ceci nest pas une surprise
puisque meme notre exemple de la corde 1D est en fait ramenee `a la connaissance du champ de deplacement
(x, t).
Pour un syst`eme mecanique bien deni, nous avons vu quen mecanique non-relativiste on avait L = T V .
Dans le cas dun syst`eme continu on aurait donc


3
L = L d v = (T V) d3 v
o`
u T et V sont respectivement la densite denergie cinetique et potentielle. Cest eectivement ce que lon a
dans le cas de la corde 1D
E

L = t2 x2
2
2
Mais pour un syst`eme physique non associe `a un syst`eme mecanique, L = T V et lon peut donc utiliser toute
expression conduisant aux equations recherchees. Ainsi, on peut trouver des densites de lagrangiens L telles
que les equations du champ (equations de Lagrange correspondantes) fournissent les equations de Maxwell
(electrodynamique), lequation dEuler (hydrodynamique), lequation dEinstein (relativite generale) ou
meme lequation de Schrodinger (mecanique quantique).
En mecanique non-relativiste, il y a parfaite equivalence entre force et potentiel (qui nest rien dautre
quun champ). Lutilisation dun potentiel est dans ce cas une simple commodite de description mais ce
nest plus le cas en relativite. En eet, la valeur nie de la vitesse de la lumi`ere (et de toute autre vitesse)
introduit un delai de propagation entre toute cause (ex: presence dune particule chargee en un point) et son
eet (champ electromagnetique vu en un autre point). Ce delai est celui pris par la propagation de cette
information, propagation qui se traduit par la deformation du champ dinteraction. Celui-ci acquiert donc
une realite intrins`eque dans le contexte de la relativite et cest pourquoi la theorie classique des champs a
une si grande importance dans la description des interactions fondamentales.
Une autre consequence importante de la relativite est que toute particule fondamentale est necessairement
ponctuelle. En eet, considerons une particule spatialement etendue et faisons-la se mouvoir sous leet
dune force appliquee en lun de ses points. A linstant meme o`
u elle recoit cette impulsion et y reagit en se
deplacant, le cote oppose de cette particule na pas encore pris connaissance de cette action et ne peut donc
se deplacer. Cela signie que cette particule se deforme dabord avant detre capable de se deplacer de facon
rigide. Si elle se deforme, cest quelle-meme poss`ede une constitution: une particule etendue ne peut donc
etre consideree comme fondamentale.
Parmi les quatre interactions fondamentales reconnues actuellement, seules deux peuvent etre decrites par
la theorie classique des champs: lelectrodynamique et la gravitation (relativite generale). Les deux autres


5.3. THEORIE
CLASSIQUE DES CHAMPS

79

(interactions faible et forte) necessitent de prendre en compte des aspects quantiques. La theorie quantique
des champs est le cadre theorique `a linterieur duquel les physiciens cherchent `a unier les 4 interactions. Ils
ont ainsi reussi `a produire une description quantique des champs electromagnetique, nucleaires fort et faible.
A lheure actuelle, seule la gravitation se derobe `a la quantication. . .
En resume, la formulation lagrangienne a plusieurs avantages:
La possibilite dinventer de nouveaux types de champs L d`es lors quils sont des invariants de Lorentz.
Toutes les analyses concernant les invariants, y compris le theor`eme de Noether, peuvent etre etendus
dans le cadre de la theorie des champs. Cela ore souvent des resultats generaux importants.
La formulation lagrangienne est le passage oblige vers la quantication.

5.3.2

Exemple:
electrodynamique classique

Proposee en 1903 par Schwarzschild, la formulation lagrangienne de lelectrodynamique permet (entre autres)
de ramener toute linformation contenue dans lensemble des quatre equations de Maxwell `a la seule donnee
de la densite de lagrangien L.
Les quatre equations de Maxwell secrivent (dans le vide):


div B


rot E


div E


rot B

= 0

(5.12)

B
=
t

=
0

(5.13)
(5.14)

= 0
j + 0 0

E
t

(5.15)

Les deux premi`eres equations sont automatiquement satisfaites si lon pose

B

rot A

E

(5.16)


A
t

(5.17)

sont appeles respectivement les potentiels scalaire et vecteur. On voit donc quil sut de se donner
o`
u et A


4 grandeurs (champs) independants A (tels que A0 soit relie `a , les 3 autres etant les composantes de A)




pour decrire compl`etement E et B. Les deux autres equations de Maxwell (5.14) et (5.15) relient les champs
aux sources (densite locale de charges et courants
j). Ce sont donc ces deux derni`eres equations quil faut
retrouver comme equations du champ.
La densite de lagrangien permettant de le faire secrit
L=

0 2


(E c2 B 2 ) +
j A
2

(5.18)

Les equations (5.16) et (5.17) secrivent tensoriellement

Bi
Ei

=
ijk j Ak
= i t Ai

o`
u lon a introduit, par commodite decriture, le tenseur dordre 3 de Levi-Civita3
ijk . Lequation de
Lagrange du champ est obtenue en calculant les derivees suivantes:
L

3 Ce

tenseur est tel que 123 = 1, ijj = 0 i, j, ijk = jik = jki i, j, k.

80

CHAPTER 5. DESCRIPTION LAGRANGIENNE DES MILIEUX CONTINUS

L
(t )
L
(i )

= 0
0 E 2
Ej
= 0 Ej
= 0 Ei
2 (i )
(i )

Lequation de Lagrange pour le champ est alors

3


d
(0 Ei ) =
dxi
i=1

et nest autre que lequation (5.14). Calculons maintenant les derivees associees au champ Ai :
L
Ai
L
(t Ai )
L
(j Ai )

ji

0 Ej

Ej
= 0 Ei
(t Ai )



m
A
x
0 c2 B 2
B
l
k


= 0 c2 Bk
= 0 c2 Bk
klm 
2 (j Ai )
(j Ai )
Ai
xj

0 c2 Bk
klm im lj = 0 c2
kji Bk =

ijk Bk
0

Lequation de Lagrange pour le champ Ai est alors

3


=0

(0 Ei ) +
t

3


j=1

L
Ai

d
L
dx ( Ai )

ijk Bk
xj 0

= ji

ijk j Bk

= 0 ji + 0 0 t Ei

qui nest autre que la i-`eme composante de lequation (5.15). On vient ainsi de prouver que lexpression
(5.18) de L proposee porte bien en elle toute linformation contenue dans les equations de Maxwell. On
voit dailleurs que cette expression comporte deux termes. Le premier met en jeu uniquement les champs
en labsence de toute source. Le deuxi`eme met en jeu et
j et decrit donc linteraction du champ avec la
mati`ere.
En presence dun grand nombre de particules (ponctuelles) de densites de charge et vitesses v
= r

dierentes, on a naturellement

j

v

Le lagrangien du champ secrit alors



L d3 v =
Lchamp =

!



0 2
2 2


(E c B )
v
A
d3 v
2

Cette expression est `a rapprocher de celle, dej`a connue, dun syst`eme de particules ponctuelles de charges
q soumises `a un champ electromagnetique exterieur

#

"1
2


Lmat = T V =
m v q v
A
2


5.3. THEORIE
CLASSIQUE DES CHAMPS

81

Dans cette expression, le deuxi`eme terme decrit linteraction du champ sur la mati`ere. En vertu du principe
daction et de reaction, ce doit etre le meme terme que dans lexpression de la densite de lagrangien pour le
champ. Pour concilier les deux expressions il sut de realiser que, pour une particule ponctuelle situee en
r
, la densite de charge associee est simplement
= q (
r r
)

(5.19)

dans le cas continu du symbole de Kronecker. La

o`
u (
r r
) est la distribution de Dirac, generalisation

condition de normalisation est tout simplement d3 v = q .
Ainsi, moyennant cette precaution liee `a laspect ponctuel des charges, on peut construire un seul lagrangien decrivant `a la fois le comportement de particules chargees et celui du champ electromagnetique
total (incluant celui cree par les charges):
Ltot =



1
m v2 +
2

0 2
(E c2 B 2 ) d3 v
2






d3 v
q (
r r
) v
A

(5.20)

Le premier terme est le lagrangien de particules isolees, le second decrit un champ electromagnetique en
labsence de sources et le troisi`eme decrit linteraction mati`ere-rayonnement. Ce lagrangien est une fonction
L = L(, Ai , , Ai , r
, v
, x, y, z, t), la dependance eventuelle dans les coordonnees x, y, z, t intervenant
ici uniquement via le volume dintegration. Noter egalement que si lon desire une expression relativiste pour
les particules (ce qui est souhaitable pour une description correcte), Lmat = T V mais plutot
Lmat =


m c

5.3.3

v2
V
c2

Tenseur
energie-impulsion dun champ

Nous avons donc la possibilite de decrire tout syst`eme physique continu en utilisant une densite de lagrangien
L (si on la connait). Connaitre la dynamique dun tel syst`eme consiste alors `a obtenir les n champs i en
tout point de lespace-temps et ceci, grace `a la resolution des n equations du champ associees
3


L
d
L
=
dx ( i )
i
=0

(5.21)

Mais, comme on la vu dans le cas des syst`emes discrets, certaines proprietes de symetrie peuvent considerablement simplier la resolution, voire meme sure `a faire emerger des comportements remarquables.
On va donc suivre la meme demarche que dans le chapitre I.
Tout dabord, si k est un champ cyclique, cest `a dire si L ne depend pas explicitement de k , alors
3
3




d L
d
L
d
L
=
+
=0
dx
(

)
dt
(

)
dx
(

k
k

k )
=0
=1

(5.22)

qui peut sinterpreter comme la nullite dune quadri-divergence dans lespace-temps, donc la conservation
dun quadrivecteur (evidemment associe au champ k ).
On peut egalement rechercher lequivalent de lintegrale premi`ere de Jacobi H associee `a un syst`eme
continu. Pour un syst`eme discret de lagrangien L on avait


d L
L
dL
=
q +
dt
dt q
t


L
d L
q L
=
(5.23)
dt q
t

82

CHAPTER 5. DESCRIPTION LAGRANGIENNE DES MILIEUX CONTINUS

Pour un syst`eme continu L, lespace joue le meme role que le temps et il faut rechercher une grandeur
invariante selon toute variation selon une coordonnee x . Il nous faut donc calculer
dL
dx

L d
L d
L
+
+
dx
dx
x


d
L d2
L
L
d
+
+
dx dx
dx dx
x


L d
d
L
+
dx dx
x

=
=
=

qui peut aussi secrire

d
dx

L d
L
dx


=

L
x

(5.24)

(5.25)

o`
u est le symbole de Kroneker. On voit donc que si L ne depend pas explicitement de la coordonnee x ,
alors la quantite
L d
T =
L
(5.26)
dx
est bien conservee, cest `a dire
dT
=0
dx

(5.27)

Quelle est la signication de T ? Cest un tenseur dordre deux deni dans lespace-temps, donc `a 16
composantes. Dans le cas general, cest un tenseur non symetrique T = T . On lappelle tenseur energieimpulsion du syst`eme car ses composantes decrivent la densite denergie et dimpulsion transportees par
celui-ci (ex: champ electromagnetique). A ce stade, sans en faire ici la demonstrations, il nous sut de
savoir que ce tenseur peut toujours se decomposer de la facon suivante:
T00

5.4.1

L=H

densite denergie

Ti0

composantes du vecteur densite dimpulsion

T0j

composantes du vecteur courant denergie

Tij

5.4

Tji

tenseur des contraintes

Compl
ements
Formulation relativiste de la th
eorie des champs

La theorie classique des champs peut tr`es facilement sexprimer de facon covariante, cest `a dire satisfaisant
au principe de relativite dEinstein. En fait, une telle description est meme presque naturelle, comme latteste
lexpression de laction elle-meme
   
S=
L dxdydzdt
(5.28)
Dans cette expression, lelement innitesimal de volume de lespace-temps (dxdydzdt) est lui-meme un invariant de Lorentz. Il sut donc davoir une densite de lagrangien invariante par transformee de Lorentz
pour que toute equation du champ, obtenue via le principe de Hamilton, soit elle-meme invariante, donc
covariante. Le gain est appreciable: avoir L scalaire invariant de Lorentz garantit une description satisfaisant
la relativite. Il sut cependant de faire attention aux bornes dintegration. . .


5.4. COMPLEMENTS

5.4.2

83

Densit
e dhamiltonien

Ce qui a ete fait pour le lagrangien peut evidemment etre fait pour lhamitonien. Une telle construction
na cependant de reel interet que pour la quantication des champs. Nous avons vu quil y avait une
equivalence formelle totale entre loptique geometrique et la mecanique. Cest de cette equivalence et de
lensemble des experiences montrant que la mati`ere avait eectivement un caract`ere ondulatoire et quantie
quest nee la mecanique quantique. Or, cette equivalence ne se voit que dans le cadre de la mecanique
hamiltonienne. La quantication de la mati`ere elle-meme (mecanique quantique, premi`ere quantication)
a opere par correspondance entre des objets propres `a cette nouvelle physique et des concepts developpes
au sein de la mecanique hamiltonienne. Cest certainement la raison pour laquelle nous ne savons pas
faire de quantication en dehors de la description hamiltonienne. La quantication des champs ou seconde
quantication nechappe pas `a cette r`egle. . .
Quoi quil en soit, et independamment de toute utilite pratique, on peut se demander ce que devient le
concept de hamiltonien dans le cas dun syst`eme ayant une innite de degres de liberte. Dans le cas discret,
nous avions introduit le moment conjugue `a partir du lagrangien. Dans le cas de la corde 1D par exemple,
le moment conjugue de la coordonnee generalisee i est
pi =
puisquon peut ecrire L =

N

L
Li
=a
i
i

aLi . Lhamiltonien discret est alors





 Li
H =
pi i L =
i aLi
a
i
i
i


 Li
a
i Li
=
i
i

i=1

N

aHi

i=1

Notre experience du passage `a la limite continue, cest `a dire lorsque a 0 (donc N ), pour le lagrangien
nous indique que lon obtient lhamiltonien





L
aHi = dx
t L
H = lim
a0
t
i

=
dx H
On denit ainsi H comme etant la densite dhamiltonien et
=

L
t

(5.29)

la densite de moment. Noter que la densite de moment conjuguee du champ est une grandeur consideree
comme ne dependant que de lespace-temps, = (x ) (independante de et de ses derivees).
Dans le cas general dun syst`eme continu decrit par n champs i (
r, t) independants, on passe donc de la
densite de lagrangien L `a la densite dhamiltonien H en faisant la transformation


H=
i t i L
(5.30)
i

les i (
r, t) etant les densites de moments conjugues des champs i (
r, t). Le hamiltonien complet du syst`eme
est calcule de la meme mani`ere, `a savoir
 x2  y 2  z2
H=
H dxdydz
(5.31)
x1

y1

z1

84

CHAPTER 5. DESCRIPTION LAGRANGIENNE DES MILIEUX CONTINUS

o`
u le volume dintegration est le meme que pour L (ie. arbitraire). La transformation (5.30) correspond `a
une transformation de Legendre qui, `a partir dune fonction L = L(, x , y , z , t , x, y, z, t) construit une
fonction
(5.32)
H = H(, x , y , z , , x, y, z, t)
Pour sen convaincre, il sut de calculer la dierentielle totale de H et de voir queectivement elle ne met
plus en jeu t . Pour un syst`eme decrit par un seul champ (commodite) on a ainsi
d(t L) =

t d + dt
t d

3
3




L
L
L
L
d
dx
dt
dx

x
x
t

=1
=0

3
3




L
L
L
d
dx
dx

x
x

=1
=0

o`
u les x correspondent aux coordonnees x, y, z, t. Or, dapr`es la denition de la densite dhamiltonien, cette
dierentielle est egale `a
dH =

3
3




H
H
H
H
d +
d +
dx +
dx

x
x

=1
=0

Ces deux expressions ne peuvent etre egales que si lon a

t
H

H
x

=
=
=

= 0, . . . , 3

La premi`ere de ces equations est lanalogue de lune des equations canoniques de Hamilton (q = H
p ). La
derni`ere indique que la dependance explicite de H en fonction des coordonnees de lespace-temps est lopposee
L
equation ci-dessus va nous
de celle de L (generalisation du cas discret o`
u lon avait H
t = t ). La seconde
es
fournir une generalisation de la seconde equation canonique de Hamilton (p = H
q ). En eet, dapr`
lequation de Lagrange des champs (5.8), on a

d
L
d L
L
=
+

dt (t ) =1 dx (x )
3

do`
u lon obtient

H
d
L
=

dx
(

x )
=1
3

En resume, les equations canoniques pour les champs secrivent

(5.33)

H
d L

dxi (xi )
i=1
3

(5.34)

Bien quelles aient perdu leur symetrie qui les rendaient si attrayantes dans le cas discret, ces equations
conservent les memes proprietes. En particulier, on peut les redemontrer directement `a partir du principe
de Hamilton, comme pour le cas discret.