Introduction au principe variationnel et

à la mécanique analytique

Notes de Cours de Licence L3-Phytem

Nicolas Sator
Laboratoire de Physique Théorique de la Matière Condensée
Université Pierre et Marie Curie Paris 6

LPTMC - Septembre 2009

 Depuis le XVII eme siècle, l’approche variationnelle permet de décrire des
phénomènes physiques à l’aide d’un principe d’économie, appelé en optique
principe de Fermat et en mécanique principe de moindre action : le comporte-
ment observé d’un système correspond à la minimisation (ou à la maximisation)
d’une certaine grandeur.
Ce cadre abstrait et général permet une reformulation simple et élégante de
nombreux problèmes en physique (optique géométrique, mécanique classique,
relativiste, électromagnétisme). Mais au-delà de considérations esthétiques ou
techniques, l’approche variationnelle forme la pierre angulaire de la mécanique
quantique, de la physique statistique, de la théorie du chaos et de la théorie
des champs. Dans ces notes de cours, nous aborderons quelques notions de
calcul variationnel (et en particulier les équations d’Euler-Lagrange) que nous
appliquerons par la suite à la mécanique analytique.

Bibliographie

– Mécanique quantique, 1. Fondements et premières applications, C. Aslan-

gul (de Boeck, Bruxelles, 2007), chapitre 7
– Analytical Mechanics, R.R. Fowles and G. L. Cassiday (Brooks/Cole Thom-
son learning, 1999), chapter 10
– Le cours de physique de Feynman, tome 1, R. Feynman (Dunod, 1999),
chapitre 19
– Mécanique, L. Landau et E. Lifchitz (Mir, Moscou, 1982)
– Mécanique, De la formulation Lagrangienne au chaos Hamiltonien, C. Gig-
noux et B. Silvestre-Brac (EDP Sciences, 2002)
– Mécanique analytique, R. Dandoloff (Publibook, 2005)

Table des matières

1 Introduction : le principe de Fermat 4

2 Approche variationnelle 5
2.1 Équation d’Euler-Lagrange . . . . . . . . . . . . . . . . . . . . . 6
2.1.1 Notion de fonctionnelle . . . . . . . . . . . . . . . . . . . 6
2.1.2 Démonstration de l’équation d’Euler-Lagrange . . . . . . 7
2.1.3 Formule de Beltrami . . . . . . . . . . . . . . . . . . . . . 8
2.1.4 Exemple : La courbe brachistochrone . . . . . . . . . . . . 9
2.2 Approche variationnelle avec des contraintes . . . . . . . . . . . . 11
2.2.1 Méthode des multiplicateurs de Lagrange . . . . . . . . . 11
2.2.2 Contraintes de type Holonôme . . . . . . . . . . . . . . . 15
2.2.3 Contraintes de forme intégrale . . . . . . . . . . . . . . . 17
2.2.4 Exemple : La chainette . . . . . . . . . . . . . . . . . . . 18

2
3 Rudiments de mécanique analytique 20
3.1 De Newton à Lagrange : une reformulation de la mécanique . . . 20
3.2 Le formalisme Lagrangien . . . . . . . . . . . . . . . . . . . . . . 20
3.2.1 Le principe de moindre action . . . . . . . . . . . . . . . . 20
3.2.2 Le Lagrangien . . . . . . . . . . . . . . . . . . . . . . . . 22
3.3 Le formalisme hamiltonien . . . . . . . . . . . . . . . . . . . . . . 24
3.3.1 Les équations de Hamilton . . . . . . . . . . . . . . . . . 24
3.3.2 Les crochets de Poisson . . . . . . . . . . . . . . . . . . . 27
3.4 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.5 Exercices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3.5.1 L’oscillateur harmonique . . . . . . . . . . . . . . . . . . . 29
3.5.2 Le cerceau à vitesse angulaire constante . . . . . . . . . . 30
3.5.3 L’atome hydrogénoı̈de . . . . . . . . . . . . . . . . . . . . 33
3.5.4 La particule plongée dans un champ électromagnétique . 34

3
1 Introduction : le principe de Fermat
Les lois de l’optique géométrique – propagation rectiligne de la lumière
dans un milieu homogène, principe de retour inverse de la lumière et lois de
la réflexion et de la réfraction de Snell (1621) et Descartes (1637)1 – ont une
origine phénoménologique.
En 1657, Pierre de Fermat propose une autre approche, plus abstraite, basée
sur l’idée que “La nature agit toujours par les voies les plus courtes et les
plus simples”2 . Le principe de Fermat s’énonce comme un principe de moindre
temps3 : La lumière se propage d’un point à un autre sur des trajectoires telles
que la durée du parcours est stationnaire. Il s’agit bien d’un principe variationnel
car la durée du parcours doit être extrémale, en générale minimale, par rapport
à une petite variation du trajet.

A l’aide de ce principe, on démontre immédiatement la loi de propagation

rectiligne de la lumière dans un milieu homogène. Les lois de Snell-Descartes se
reformulent à la Feynman par une analogie cinématique.

Le problème du maı̂tre nageur

Telle la lumière qui se propage moins vite dans l’eau que dans l’air, un maı̂tre
nageur court plus vite qu’il ne nage. Il se trouve au point A lorsqu’il aperçoit
une jolie fille qui se noie en B. Comment arriver en B le plus vite possible ?
Sachant que le maı̂tre nageur court en ligne droite à la vitesse v1 et nage
à la vitesse v2 < v1 , quel est le point I où il doit plonger ? Il faut trouver un
compromis entre la ligne droite et le parcours qui rend minimal le trajet dans
l’eau. La réponse mathématique a été donnée par Maupertuis en 1744.
Comme on le voit sur la figure 1, la trajectoire du maı̂tre nageur est con-
stituée de deux droites AI et IB, où I(x, 0) est le point où le maı̂tre nageur
plonge. A priori la distance AI sera plus grande que la distance IB car il court
plus vite qu’il ne nage.
Le temps T (x) mis par le maı̂tre nageur pour aller de A en B est :
AI IB
T = +
v1 v2
soit, √ p
x2 + a2 (d − x)2 + b2
T (x) = + .
v1 v2
1 En fait le mathématicien arabe Ibn Sahl aurait découvert la loi de la réfraction dès 984.
2 Lettre de Clerselier (cartésien, éditeur de Descartes) à Fermat (1662) : “Le principe que
vous prenez pour fondement de votre démonstration, à savoir que la nature agit toujours par
les voies les plus courtes et les plus simples, n’est qu’un principe moral et non point physique,
qui n’est point et qui ne peut être la cause d’aucun effet de la nature”.
3 Une formulation plus moderne repose sur la notion de chemin optique (L =
R
C ndl, où
n = c/v est l’indice de réfraction et c la vitesse de la lumière dans le vide) : Le trajet suivi
par la lumière entre deux points A et B est celui qui correspond à une valeur stationnaire du
chemin optique.

4
 Fig. 1 – Le problème du maı̂tre nageur en A et de la fille qui se noie en B.

Le temps est minimal (il peut pas s’agir d’un maximum !) quand

dT 1 x 1 d−x
= √ − p = 0,
dx v1 x2 + a2 v2 (d − x)2 + b2

soit
1 1
sin(i1 ) = sin(i2 ).
v1 v2
La vitesse de la lumière dans un milieu d’indice de réfraction n étant v = c/n, où
c est la vitesse dans le vide, on retrouve la loi de la réfraction de Snell-Descartes.
On montrerait de même la loi de la réflexion.

2 Approche variationnelle
Dans le cas de la loi de la réfraction, la trajectoire du rayon lumineux était
déterminée par la valeur de x qui minimise T (x). En général, c’est toute une
fonction (par exemple une trajectoire) que l’on recherche.

Par exemple, quel est le plus court chemin entre deux points A et B d’un plan ?

En notant s l’abscisse curviligne le long d’un chemin, on a

p p
ds = dx2 + dy 2 = 1 + y ′2 dx,

5
 y

A B

ds
dy
dx

x
Fig. 2 – Le chemin le plus court entre deux points A et B.

dy
où y ′ (x) = dx . La longueur l du chemin est donnée par
Z B Z xB p
l[y] = ds = 1 + y ′2 dx.
A xA

On cherche donc la fonction y(x) qui rend l’intégrale l[y] minimale telle que
y(xA ) = y(xB ) = cte. La réponse est bien sûr y(x) = cte, mais comment le
démontrer ?

Comme nous le verrons plus loin, un autre exemple plus intéressant est celui
de la courbe brachistochrone qui est la trajectoire qui minimise le temps de
parcours entre deux points dans le champ de pesanteur.

2.1 Équation d’Euler-Lagrange

2.1.1 Notion de fonctionnelle
La fonctionnelle est une généralisation de la fonction. Plutôt que de dépendre
d’une (ou plusieurs) variable, une fonctionnelle dépend d’une (ou plusieurs)
fonction, elle même dépendant d’une (ou plusieurs) variable. C’est donc une
“fonction de fonction”.
Considérons une fonction y(x) dépendant d’une
variable x et une fonction f
dy
à trois variables non indépendantes f y(x), y ′ (x), x , où y ′ (x) = dx . La fonction
y(x) étant connue, la fonction f prend donc une valeur déterminée pour une
valeur de x donnée. Nous verrons que la dépendance en la dérivée y ′ (x) intervient
fréquemment, en particulier en mécanique.
Remarquons que la fonction f dépend de x explicitement, mais aussi im-
plicitement à travers y(x) et y ′ (x). La dérivée totale de f par rapport à x est

6
donc :
df ∂f ∂f dy ∂f dy ′ ∂f ∂f ′ ∂f ′′
= + + ′ = + y + ′y .
dx ∂x ∂y dx ∂y dx ∂x ∂y ∂y
Définissons la fonctionnelle S[y] par l’intégrale :
Z x2


S[y] = f y(x), y ′ (x), x dx.
x1

La fonctionnelle S[y] est donc un nombre qui dépend de la fonction y(x). L’ap-
proche variationnelle va nous permettre de déterminer la fonction y(x) telle que
S[y] soit stationnaire (extrémale), sachant que y(x1 ) et y(x2 ) sont donnés.

Nous allons montrer que la fonction y(x) qui remplit cette condition doit
vérifier l’équation d’Euler-Lagrange :

∂f d ∂f
− = 0.
∂y dx ∂y ′

2.1.2 Démonstration de l’équation d’Euler-Lagrange

Supposons que l’on connaisse la fonction y0 (x), qui rende S extrémale.
Puisque S[y0 ] est stationnaire, une petite variation η(x) de la fonction y(x)
implique une variation δS = 0 au premier ordre en η(x). Posons explicitement :

y(x) = y0 (x) + η(x),

où ∀x, η(x) << y0 (x) et calculons la variation induite de la fonctionnelle pour
une valeur de x fixée :
Z x2  
′ ′

′


δS = f y0 (x) + η(x), y0 (x) + η (x), x − f y0 (x), y0 (x), x dx.
x1

Au premier ordre en η(x) et en η ′ (x), on a

∂f ∂f
f y0 (x) + η(x), y0′ (x) + η ′ (x), x ≃ f y0 (x), y0′ (x), x + η(x) + ′ η ′ (x).
∂y ∂y
Donc x2  
∂f ∂f
Z
δS = η(x) + η ′ (x) ′ dx.
x1 ∂y ∂y
En intégrant par parties la seconde intégrale, on obtient :
Z x2  
∂f d ∂f  ∂f x2
δS = η(x) − η(x) dx + η(x) ′ x1 .
x1 ∂y dx ∂y ′ ∂y

Puisque y(x1 ) et y(x2 ) sont fixés, η(x1 ) = η(x2 ) = 0 et le dernier terme de

l’équation s’annule. Il reste :
Z x2  
∂f d ∂f
δS = η(x) − dx.
x1 ∂y dx ∂y ′

7
Comme δS = 0, quelle que soit η(x) on doit avoir :

∂f d ∂f
− = 0.
∂y dx ∂y ′
C’est l’équation d’Euler-Lagrange, qui s’écrit plus explicitement :

∂f ∂2f ∂2f ∂2f

− ′ − ′ y ′ − ′2 y ′′ = 0.
∂y ∂y ∂x ∂y ∂y ∂y

Reprenons l’exemple du chemin le plus court :

Z xB
l[y] = f (y ′ )dx,
xA
p
où f (y ′ ) = 1 + y ′2 avec y(xA ) = y(xB ) = cte. Puisque f ne dépend pas de y,
on a d’après l’équation d’Euler :

∂f y′
=p = c,
∂y ′
1 + y ′2

où c est une constante. On trouve donc y ′ = cte, autrement dit, l’équation
d’une droite qui passe par les points A et B. Dans cet exemple la réponse était
évidente, mais l’approche variationnelle est très utile lorsque l’espace-temps est
doté d’une métrique plus complexe.

2.1.3 Formule de Beltrami

Si la fonction f ne dépend pas explicitement de la variable x, ( ∂f
∂x = 0), on
obtient la formule de Beltrami :
∂f
f − y′ = C, (1)
∂y ′
où C est une constante.
Démonstration :
Calculons :
d  ∂f  ∂f ′ ∂f ′′ ∂f  ∂ 2 f ′ ∂ 2 f ′′ 
f − y′ ′ = y + ′ y − y ′′ ′ − y ′ y + ′2 y
dx ∂y ∂y ∂y ∂y ∂y ′ ∂y ∂y
2 2
 ∂f ∂ f ∂ f 
= y′ − ′ y ′ − ′2 y ′′ .
∂y ∂y ∂y ∂y
D’après l’équation d’Euler, le terme entre crochets est nul, ce qui démontre la
formule.

8
2.1.4 Exemple : La courbe brachistochrone
La courbe brachistochrone4 a été étudiée par Leibniz, Newton, Jacques et
Jean Bernoulli et Euler : En partant d’une position donnée sans vitesse initiale,
quelle est la trajectoire qui permet d’atteindre le plus rapidement possible une
position finale dans le champ de pesanteur ?

Pour préciser les notations, supposons que l’objet parte du point A de co-
ordonnées xA = 0 et yA = h > 0, sans vitesse initiale et glisse le long d’un
toboggan jusqu’au point B de coordonnées xB = d > 0 et yB = 0 (voir la figure
3). On doit donc exprimer la durée du parcours entre A et B, puis chercher la
courbe (la forme du toboggan) qui minimise cette durée.

h A

B
d x
Fig. 3 – Un objet part du point A sans vitesse initiale et glisse le long de la
courbe jusqu’au point B.

Soit s l’abscisse curviligne, un élément de courbe ds s’exprime en coordonnées

cartésiennes p p
ds = dx2 + dy 2 = 1 + y ′2 dx,
dy ds
où y ′ = dx . La vitesse de l’objet est donnée par v = dt . La conservation de
l’énergie mécanique implique
1
mv 2 + mgy = mgh.
2
Exprimons le temps de parcours T [y]. C’est une fonctionnelle qui dépend de la
forme y(x) du toboggan. On a
p
dx dx ds vdx v 2g(h − y)
= =p = p = p ,
dt ds dt 2
dx + dy 2 1+y ′2 1 + y ′2
soit s
1 + y ′2
dt = dx.
2g(h − y)
4 En grec, “Brakhisto” signifie ”le plus court” et “chronos”, le temps.

9
Le temps de parcours s’écrit donc :
Z d
T [y] = f (y, y ′ )dx,
0
s
′ 1 + y ′2
où f (y, y ) = ne dépend pas de x. On cherche donc la forme y(x)
2g(h − y)
qui minimise la fonctionnelle T [y]. Puisque f est indépendante de x, on peut
utiliser la formule de Beltrami
∂f
f − y′ = C,
∂y ′
où C est une constante. Donc
1 p
p = C 2g(h − y).
1 + y ′2
En élevant au carré on a donc :
R
1 + y ′2 = ,
h−y
où R = 1/2gC 2 . Posons y ′ = tan(θ/2) :
θ 1 R
1 + y ′2 = 1 + tan2 = θ
= ,
2 cos2 2
h−y
soit
θ R
y = h − R cos2 = h − (1 + cos θ).
2 2
Par ailleurs,
dy dy dθ R dθ θ θ dθ
y ′ = tan(θ/2) = = = sin θ = R sin cos ,
dx dθ dx 2 dx 2 2 dx
donc
θ
dx = R cos2 dθ,
2
soit en intégrant
R R
x= (θ + sin θ) + cte et y = h − (1 + cos θ),
2 2
que l’on peut écrire avec le changement de variable θ → θ − π et en imposant
x(0) = 0 :
R R
x = (θ − sin θ) et y = h − (1 − cos θ).
2 2
C’est l’équation d’une cycloı̈de. Reste à déterminer R pour que la trajectoire
passe par le point B. Sur la figure 4, on remarque que la pente à l’origine est
infinie pour obtenir une accélération initiale maximale5 .
5 On peut également s’intéresser à la courbe tautochrone (“tauto” le même) telle que tout

objet lâché sans vitesse initiale de n’importe quel point de la courbe arrive à son extrémité
en un temps indépendant du point de départ (mais pas à la même vitesse). On retrouve une
cycloı̈de.

10
 1

0.8

0.6

0.4

y
0.2

-0.2
0 0.5 1 1.5 2
x

Fig. 4 – La courbe brachistochrone entre le point A(0, 1) et le point B(2, 0) pour

R ≃ 1.035.

2.2 Approche variationnelle avec des contraintes

Dans un problème aux variations, il peut exister une (ou plusieurs) contrainte
sur la fonction y(x). Le problème isopérimétrique, en est une belle illustration :
Dans le plan, quelle est la courbe fermée qui a la plus grande aire pour un
périmètre donné ?
Dans le plan muni d’un répère cartésien, on cherche l’équation y(x) de cette
courbe telle que l’aire I
A[x, y] = ydx (2)

soit maximale avec la contrainte suivante : le périmètre

I p
l[x, y] = 1 + y ′2 dx (3)

est fixé.

2.2.1 Méthode des multiplicateurs de Lagrange

Pour simplifier, nous raisonnerons dans un premier temps à l’aide d’une
fonction de plusieurs variables. La généralisation à la fonctionelle ne posant pas
de problème.

Par exemple, un fabriquant souhaite produire un récipient cylindrique de

volume fixé V dont la surface S soit la plus faible possible (pour minimiser
la quantité de matériaux). Comment choisir la hauteur L et le rayon R du
cylindre ?
La surface est une fonction de L et R : S(L, R) = 2πR2 + 2πRL. On cherche
le minimum de S avec la contrainte V = πR2 L, où V est fixé. En remplaçant
L = V /πR2 dans S(L, R) on se ramène facilement à une fonction à une variable

11
dont il est facile de trouver le minimum obtenu pour R = (V /2π)1/3 , soit un
rapport L = 2R (qui se rapproche le plus d’une sphère).

Généralisons : Soit f (x1 , x2 , · · · , xn ) une fonction

de n variables. Pour que
f soit extrémale en un certain point x01 , x02 , ..., x0n , il faut que
X ∂f
df = dxi = 0. (4)
i
∂xi

Si toutes les variables x1 , x2 , ..., xn sont indépendantes, cette condition est vérifiée
pour
∂f
= 0, i = 1 · · · n. (5)
∂xi
Si les variables xi vérifient une contrainte supplémentaire sous la forme (con-
trainte de type holonôme)

g(x1 , x2 , · · · , xn ) = 0, (6)

les variables xi ne sont plus indépendantes et l’extremum de la fonction n’est

plus obtenu par l’équation (5).

En revanche, d’après l’équation (6), les variations dxi sont liées par :
X ∂g
dg = dxi = 0. (7)
i
∂xi

Puisque l’équation (7) relie les variables dxi entre elles, une façon de procéder
est d’exprimer dxn en fonction de dx1 , dx2 ,...,dxn−1 et de remplacer dxn dans
l’équation (4). On se ramène donc à un problème de n−1 variables indépendantes
avec n − 1 équations à résoudre (obtenues par l’annulation des coefficients des
dxi )...ce qui peut être très compliqué !

Une alternative efficace est d’introduire un paramètre µ qui sera calculé a

posteriori. L’ajout d’une variable pour lever une contrainte est assez courant en
physique. Voyons deux façons différentes d’introduire cette méthode :

i) Quel que soit µ, on peut écrire :

X ∂f ∂g

df + µdg = +µ dxi = 0.
i
∂xi ∂xi

Maintenant, choisissons µ tel que

∂f ∂g
+µ = 0.
∂xn ∂xn
On se ramène à nouveau à un problème de n − 1 variables indépendantes et
donc
∂f ∂g
+µ = 0, pour i = 1, 2, ..., n − 1.
∂xi ∂xi

12
Finalement, cette dernière équation est valable pour tout i = 1, 2, ..., n, comme
si toutes les variables xi étaient indépendantes. Elles ne le sont pas et s’expri-
ment toutes en fonction de µ. Reste à calculer le multiplicateur de Lagrange en
remplaçant les xi dans la contrainte (6).

ii) Adoptons un point de vue géométrique. Au cours d’une promenade en mon-

tagne, on cherche le point le plus élevé. Ici f (x, y) correspond à l’altitude du
point de coordonnées (x, y) et g(x, y) définit la trajectoire du promeneur qui est
imposée.

Fig. 5 – La fonction f (x, y) est l’altitude du point de coordonées (x, y) (sur une
carte) et le chemin parcouru (en rouge) est imposé.

Aux points où g coupe une ligne de niveau (f (x, y) = cte), le chemin est
en pente. Un extremum (stationnaire) correspond forcément à un point où le
chemin est tangent à une ligne de niveau. Autrement dit, au point extrémal les
vecteurs gradients sont colinéaires :
∇g = µ∇f,
où µ est le multiplicateur de Lagrange. On a donc
∂f ∂g
−µ = 0
∂x ∂x
∂f ∂g
−µ = 0.
∂y ∂y
On peut bien sûr ajouter des contraintes (en nombre inférieur au nombre de
variables).

Illustrons la méthode des multiplicateurs de Lagrange par deux exemples :

• Soit une fonction f (x, y, z) dont on cherche l’extremum avec la contrainte
g(x, y, z) = 0.
La fonction est extrémale pour
∂f ∂f ∂f
df = dx + dy + dz = 0. (8)
∂x ∂y ∂z

13
 Et la contrainte impose
∂g ∂g ∂g
dg = dx + dy + dz = 0. (9)
∂x ∂y ∂z

En exprimant dz en fonction de dx et dy à l’aide de l’équation (9) on a :

∂g ∂g
∂y
dz = − ∂x
∂g
dx − ∂g
dy
∂z ∂z

que l’on introduit dans l’équation (8) :

 ∂g   ∂g 
∂f ∂x ∂f ∂f ∂y ∂f
− ∂g
dx + − ∂g
dy = 0.
∂x ∂z
∂z ∂y ∂z
∂z

Et finalement :
∂g
∂f ∂x ∂f
= ∂g
∂x ∂z
∂z
∂g
∂f ∂y ∂f
= ∂g
.
∂y ∂z
∂z

Avec la méthode des multiplicateurs de Lagrange, on introduit le paramètre

µ et on recherche l’extrémum de f + µg soit :
∂f ∂g
+µ = 0
∂x ∂x
∂f ∂g
+µ = 0
∂y ∂y
∂f ∂g
+µ = 0
∂z ∂z
On choisit µ tel que :
∂f
µ = − ∂z
∂g
∂z
et on utilise cette expression de µ pour résoudre le système de deux
équations restant. On retrouve immédiatement les expressions obtenues
par la méthode directe6 .
6 Par exemple : f = x + y + z et g = x2 + y 2 + z 2 . On a donc
1 + 2µx = 0
1 + 2µy = 0
1 + 2µz = 0
1 1 2
qui implique x = y = z = − 2µ et la contrainte impose la valeur de µ : g = 3(− 2µ ) = 1 soit
√
x = y = z = ±1/ 3.

14
 • Reprenons l’exemple de la boı̂te cylindrique. Introduisons le multiplicateur
de Lagrange µ. On cherche les extrema de S + µV . Soit
∂S ∂V
+µ = 0,
∂L ∂L
∂S ∂V
+µ = 0.
∂R ∂R
Donc

2πR + µπR2 = 0,
4πR + 2πL + 2µπRL = 0.

Donc

R = −2/µ,
L = −4/µ.

On doit donc choisir L = 2R. On en déduit µ à l’aide de la contrainte :

V = πR2 L = −16π/µ3 donc µ = −( 16π
V )
1/3
.

Pour résumer : Trouver l’extremum d’une fonction à plusieurs variables liées

par une contrainte (6), revient à chercher l’extremum de la fonction f + µg
dont les variables sont supposées indépendantes. Le multiplicateur µ est cal-
culé a posteriori à l’aide de la contrainte. En général, on introduit autant de
multiplicateurs de Lagrange qu’il y a de contraintes.

2.2.2 Contraintes de type Holonôme

La méthode des multiplicateurs de Lagrange se généralise au cas d’une fonc-
tionnelle à plusieurs variables7 : Soit une fonction


f y1 (x), y2 (x), ...yn (x), y1′ (x), y2′ (x), ..., yn′ (x), x .

Si la contrainte sur les fonctions y1 , y2 ...yn peut s’écrire sous la forme

g y1 (x), y2 (x), ...yn (x) = 0, (10)

on dit que la contrainte est holonôme8 . L’extremum de la fonctionnelle

Z x2


S[y1 , y2 ..., yn ] = f y1 (x), y2 (x), ...yn (x), y1′ (x), y2′ (x), ..., yn′ (x), x dx, (11)
x1

soumise à la contrainte (10), est obtenu en recherchant l’extremum de la fonc-

tionnelle
Z x2  
′ ′ ′



f y1 , y2 , ...yn , y1 , y2 , ..., yn , x − µ(x)g y1 , y2 , ...yn dx.
x1
7 Consulter par exemple le chapitre 3.2) de “Mécanique, De la formulation Lagrangienne
au chaos hamiltonien” C. Gignoux et B. Silvestre-Brac (EDP Sciences,2002).
8 Un contre exemple : lorsque la contrainte s’écrit sous la forme d’une inégalité.

15
La fonction µ(x) est déterminée a posteriori par la contrainte (10).

Exemple : Les géodésiques sur la sphère

Quelle est le chemin le plus cours entre deux points à la surface d’une sphère ?

Soit une sphère de centre O et de rayon R et deux points A et B placés sur

sa surface. La position d’un point à la surface de la sphère obéit à la contrainte :
x2 (t) + y 2 (t) + z 2 (t) = R2 , (12)
où t est un paramètre (pas forcément le temps). L’arc élémentaire a une longueur
ds2 = dx2 + dy 2 + dz 2 = (x′2 + y ′2 + z ′2 )dt2 .
La longueur du chemin est donc
Z B p
l[x, y, z] = x′2 + y ′2 + z ′2 dt.
A

On recherche les fonctions x(t), y(t) et z(t) telle que l[x, y, z] soit minimale avec
la contrainte (12). On introduit donc la fonction µ(t) et on cherche à minimiser
RB
l’intégrale A h(x, y, z, x′ , y ′ , z ′ )dt où
p
h(x, y, z, x′ , y ′ , z ′ ) = x′2 + y ′2 + z ′2 − µ(t)(x2 + y 2 + z 2 − R2 )
L’équation d’Euler-Lagrange pour la fonction x donne
∂h d ∂h d x′
= −2µx = = p
∂x dt ∂x′ dt x′2 + y ′2 + z ′2
de même pour y et z. On peut donc écrire :
du
= −2µr, (13)
dt
où u est le vecteur unitaire tangent à la trajectoire :
1 dr
u= p .
x + y + z dt
′2 ′2 ′2

En prenant le produit vectoriel de l’équation (13) avec r on obtient

du
r× =0
dt
dr
et comme u et dt sont colinéaires
dr
× u = 0.
dt
On en déduit que
d
(r × u) = 0.
dt
Le vecteur r × u est donc constant entre les points A et B. Les vecteurs r et
u sont dans le plan équatorial passant par A et B. Le chemin le plus court est
donc l’arc de grand cercle (il y a deux chemins, mais un seul est le plus court).

16
2.2.3 Contraintes de forme intégrale
Dans les problèmes aux variations, la contrainte peut également s’exprimer
sous la forme d’une intégrale :
Z x2


C[y1 , y2 , ...yn ] = g y1 (x), y2 (x), ...yn (x), y1′ (x), y2′ (x), ..., yn′ (x), x dx. (14)
x1

Comme dans le cas précédent, l’extremum de la fonctionnelle (11) sera obtenu

en recherchant l’extremum de la fonctionnelle
Z x2  



f y1 , y2 , ...yn , y1′ , y2′ , ..., yn′ , x − µ g y1 , y2 , ...yn dx.
x1

Le multiplicateur de Lagrange µ est déterminé a posteriori par la contrainte

(14).

Reprenons le problème isopérimétrique : on cherche l’équation y(x) de la courbe

sous forme paramétrique x(t), y(t), où t estpun paramètre. p
Un élément de courbe
positionné en (x, y) a une longueur ds = dx2 + dy 2 = x′2 + y ′2 dt. La con-
trainte sur le périmètre est exprimée sous la forme d’une intégrale :
I p
l[x, y] = x′2 + y ′2 dt = cte. (15)

L’aire que l’on souhaite maximale est donnée par :

I I
A[x, y] = ydx = y(t)x′ (t)dt, (16)

ou encore, en intégrant par parties :

I
A[x, y] = − x(t)y ′ (t)dt. (17)

En utilisant la méthode des multiplicateurs de Lagrange, on recherche donc

l’extremum de I
S[x, y] = h1 (y, y ′ , x′ )dx,
p
où h1 (y, y ′ , x′ ) = yx′ − µ x′2 + y ′2 ne dépend pas explicitement de t et de x.
L’équation d’Euler-Lagrange pour x donné donc ( ∂h ∂x = 0, variable cyclique) :
1

d ∂h1
= 0,
dt ∂x′
donc
∂h1 x′
= y − µ = b, (18)
∂x′
p
x′2 + y ′2
où b est une constante. On pourrait exprimer l’équation d’Euler-Lagrange pour
y mais il est plus facile d’utiliser à présent la fonctionnelle (17) et donc h2 =

17
 p
−xy ′ − µ x′2 + y ′2 . L’équation d’Euler-Lagrange pour y donne donc ( ∂h
∂y = 0,
2

variable cyclique) :
d ∂h2
= 0,
dt ∂y ′
donc
∂h2 y′
= −x − µ p = −a, (19)
∂y ′ x′2 + y ′2
où a est une constante. En multipliant l’équation (18) par 2y ′ et l’équation (19)
par 2x′ on obtient en ajoutant ces deux équations :
d 2
x + y 2 − 2ax − 2by = 0.


2(xx′ + yy ′ − ax′ − by ′ ) =
dt
Soit
x2 + y 2 − 2ax − 2by = (x − a)2 + (y − b)2 − a2 − b2 = c2 ,
où c est une constante.
√ On trouve bien sur l’équation d’un cercle centré en (a, b)
et de rayon R = a2 + b2 + c2 . L’intérêt de la paramétrisation est que l’on peut
calculer facilement le multiplicateur de Lagrange. En effet, avec x(t) = a+R sin t
et y(t) = a + R cos t on trouve µ = R. On vérifie également que R = l/2π et
A = l2 /4π. Il est facile de se convaincre que cet extremum est bien un maximum.

2.2.4 Exemple : La chainette

On considère un fil pesant inextensible et homogène de masse linéique λ = 1
et de longueur l attaché à ses deux extrémités en x = 0 et en x = a. Déterminer
sa forme à l’équilibre9 .

Réponse :
Soit y(x) l’équation de cette forme. On a les conditions aux limites y(0) =
y(a) = 0 et la contrainte :
Z Z p Z ap
l = ds = dx2 + dy 2 = 1 + y ′2 dx = constante. (20)
0

Chaque élément de corde de longueur ds et de masse λds a pour coordonnées

(x, y(x)). Les coordonnées du centre de gravité sont

1 a p
Z
Xg (y ′ , x) = x 1 + y ′2 dx,
l 0
et a
1
Z p
Yg (y, y ′ ) = y 1 + y ′2 dx.
l 0
9 Cette
question a intéressé de nombreux savants comme Galilée, Leibniz, Jean Bernoulli
et Huygens.

18
La condition d’équilibre impose que le centre de gravité soit le plus bas possible :
yg minimale avec la contrainte (20). En utilisant la méthode des multiplicateurs
de Lagrange, on recherche donc l’extremum de
Z a
S[y] = h(y, y ′ )dx,
o
p p
où h(y, y ′ ) = y 1 + y ′2 + µ 1 + y ′2 qui est indépendante de x, on peut donc
utiliser la formule de Beltrami :
∂h p y ′2
h − y′ = (y + µ) 1 + y ′2 − (y + µ) p = C,
∂y ′
1 + y ′2

soit
y+µ
p = C.
1 + y ′2
Donc
y+µ 2
1 + y ′2 = ( )
C
y+µ
Posons z = C , on a

1p 2
1 + C 2 z ′2 = z 2 soit z′ = ± z − 1.
C
Une solution est de la forme z = cosh (αx + β) avec α = 1/C. Donc
x
y(x) = C cosh( + β) − µ.
C
a
Or y(0) = y(a) = 0 donc β = −a/2C et µ = C cosh( 2C ). Ainsi :

x − a/2 a

y(x) = C cosh( ) − cosh( ) .
C 2C
Reste à déterminer la constante C à l’aide de la contrainte (20) qui donne
l’équation implicite :
a
l = 2C sinh .
2C
On en déduit également

1 a p aC a a
Z
Xg = x 1 + y ′2 dx = sinh = ,
l 0 l 2C 2

où l’on a utiliser l’équation issue de la contrainte. On en déduit

a
Yg = C(1 − cosh ).
2C
La forme de la chainette est donc un cosinus hyperbolique (et non pas une
parabole comme le pensait Galilée).

19
3 Rudiments de mécanique analytique
3.1 De Newton à Lagrange : une reformulation de la mécanique
L’expérience montre que le mouvement d’un objet ponctuel est complètement
déterminé par sa position q(t) et sa vitesse q̇(t) à un instant donné t (on se place
pour l’instant dans un espace à une dimension). En effet, l’équation fondamen-
tale de la dynamique fournit l’expression de q̈(t) en fonction de q(t) et q̇(t) :


q̈(t) = f q(t), q̇(t) .

Connaissant les forces qui s’exercent sur le système et donc l’expression de la

fonction f , la trajectoire est construite localement, de proche en proche, en
passant de l’état mécanique à l’instant t à celui en t + dt : Ainsi, q̇(t) permet de
calculer q(t + dt) et q̈(t) permet de calculer q̇(t + dt). Le formalisme newtonien
est donc basé sur une approche différentielle.

Il existe une autre approche de la mécanique basée sur le principe de moin-

dre action10 : parmi tous les chemins possibles entre deux points donnés de
l’espace, la trajectoire physique (réelle) d’un objet est le chemin qui minimise
(ou en général “extrémise”) une certaine quantité définie sur toute la trajec-
toire, appelée action. On retrouve une idée très générale en physique qui a le
mérite d’être simple et élégante : les lois de la nature obéissent à un principe
d’économie. Contrairement à l’approche newtonienne, il s’agit donc d’un critère
global qui caractérise l’ensemble de la trajectoire.

L’enjeu n’est donc pas de trouver de nouvelles lois de la physique, mais de

reformuler l’approche newtonienne sous la forme du principe de moindre action.
Comment définir l’action de façon à retrouver les équations du mouvement ?
Le principe variationnel va nous permettre de répondre à cette question. Au-
delà de cette reformulation de la mécanique newtonienne, cette approche a des
développements fondamentaux en mécanique quantique, en électromagnétisme,
en relativité générale et plus généralement en théorie des champs.

3.2 Le formalisme Lagrangien

3.2.1 Le principe de moindre action
Un objet part d’un point A1 à l’instant t1 et arrive au point A2 à l’instant
t2 . Comme le montre la figure 6, il existe une infinité de chemins entre ces deux
points, pourtant l’objet en choisit un seul, la trajectoire physique.
Pour utiliser le principe variationnel, il faut introduire une fonctionnelle qui
dépende du chemin q(t), fictif ou réel, suivi par l’objet. Le plus simple est de
définir cette fonctionnelle comme une intégrale sur le temps le long du chemin
considéré. Il s’agit bien d’une quantité qui caractérise globalement le chemin. La
10 C’est Pierre-Louis de Maupertuis qui énonce ce principe “si digne de l’Être suprême” en

1744 et Joseph-Louis de Lagrange qui, en 1787, en donne la formulation complète que nous
connaissons aujourd’hui.

20
 q
A1 q1
q0

q
2
A2
q
3

t1 t2 t

Fig. 6 – Différents chemins qi (t) (i = 0, 1, 2, 3) en fonction du temps t. La

trajectoire réelle q0 représentant un mouvement de chute libre dans le champ de
pesanteur est représentée en traits pleins.

fonction à intégrer dépend naturellement du chemin considéré q(t), ainsi que de

sa dérivée par rapport au temps q̇(t). Plus généralement, cette fonction peut
dépendre explicitement du temps. On définit donc l’ action par
Z t2
S[q] = L(q(t), q̇(t), t)dt, (21)
t1

où la fonction L(q(t), q̇(t), t) est appelée le Lagrangien du système. Nous verrons
plus loin comment écrire cette fonction11 .
Remarquons que la fonction L dépend du temps explicitement, mais aussi
implicitement à travers q(t) et q̇(t). La dérivée de L par rapport à t est donc :
dL ∂L ∂L dq ∂L dq̇ ∂L ∂L ∂L
= + + = + q̇ + q̈.
dt ∂t ∂q dt ∂ q̇ dt ∂t ∂q ∂ q̇
Pour un chemin donné, tel que q(t1 ) et q(t2 ) sont fixés, S[q] prend une certaine
valeur. D’après le principe de moindre action, la trajectoire physique est le
chemin qui donne une valeur minimale12 à l’action. La trajectoire est donc
donnée par l’équation d’Euler-Lagrange :

∂L d ∂L
− = 0. (22)
∂q dt ∂ q̇

Plus généralement, le principe de moindre action implique que la dynamique

d’une quantité physique (position d’une particule, valeur d’un champ,...) se
11 Remarquons que le Lagrangien est défini à une fonction près, dépendant seulement du
dg(t)
temps, que l’on peut toujours écrire comme une dérivée totale, dt . La contribution à l’action
est donc une constante, g(t2 ) − g(t1 ), qui ne modifie pas l’extremum obtenu. On appelle cette
propriété l’invariance de jauge.
12 Il s’agit en général d’un minimum, d’où le nom de principe de moindre action, mais dans

certain cas l’extremum est un maximum.

21
déduit d’une unique grandeur, appelée action, qui doit être minimale pour
représenter la réalité.

3.2.2 Le Lagrangien
Comment écrire le Lagrangien d’un système ? Considérons pour commencer
une particule libre. D’après l’équation fondamentale de la dynamique13 , on doit
avoir
d
mq̈ = (mq̇) = 0. (23)
dt
Le Lagrangien d’une particule libre ne dépend a priori que de la vitesse q̇.
L’équation d’Euler-Lagrange impose donc
d ∂L(q̇)
= 0.
dt ∂ q̇
En comparant avec l’équation fondamentale de la dynamique (23), on peut
choisir
1
L = mq̇ 2 .
2
Dans le cas où le système est soumis à une force conservative14, il est raisonnable
d’écrire le Lagrangien sous la forme
1 2
L(q, q̇) = mq̇ + f (q),
2
où f est une fonction qui ne dépend que de la position. Retrouver l’équation
fondamentale de la dynamique à l’aide de l’équation d’Euler-Lagrange impose
f (q) = −U (q), où U est l’énergie potentielle d’interaction. Soit :
1 2
L(q, q̇, t) = mq̇ − U (q). (24)
2
Ainsi, dans le cas où U = U (q) l’énergie potentielle ne dépend que de la coor-
donnée q, l’équation d’Euler-Lagrange implique
∂U
mq̈ = − .
∂q
D’une façon générale, on écrira le Lagrangien comme

L = T − U,
13 Une méthode alternative est de raisonner à partir des symétries du système (Mécanique,

L. Landau et E. Lifchitz (Mir, Moscou, 1982)). Par symétrie, le Lagrangien d’une particule
libre ne peut dépendre explicitement de sa position q (invariance par translation d’espace),
du temps t (invariance par translation du temps) et de la direction de sa vitesse (isotropie
de l’espace). C’est donc une fonction de la norme de la vitesse uniquement (qui peut a priori
varier au cours du temps). A une dimension on a donc L = f (q̇ 2 ). Le principe de relativité de
Galilée impose ensuite la forme du Lagrangien.
14 L’utilisation exclusive de forces conservatives en mécanique analytique peut sembler être

un handicap, mais la dissipation de l’énergie n’intervient que de façon effective pour des
systèmes macroscopiques. Sur le plan fondamental, pour des systèmes microscopiques étudiés
dans le cadre de la physique quantique, cette restriction ne pose pas de problème.

22
où T est l’énergie cinétique et U l’énergie potentielle.

Considérons par exemple un pendule de longueur l décrit par l’angle θ (une

dimension). Son Lagrangien est (U (θ) = −mgl cos θ, où θ est l’angle avec la
verticale) :
m
L = l2 θ̇2 + mgl cos θ,
2
et l’équation d’Euler-Lagrange implique :
∂L d ∂L
= −mgl sin θ = = ml2 θ̈.
∂θ dt ∂ θ̇
r
g
On retrouve l’équation du mouvement du pendule θ̈ + sin θ = 0.
l

Dans le cas d’un système constitué de N particules dans un espace à trois

dimensions, on introduira les coordonnées généralisées qi , où i = 1, 2, ..., 3N et
les vitesses correspondantes q̇i . En coordonnées cartésiennes on a donc q1 = x1 ,
q2 = y1 , q3 = z1 , q4 = x2 , q5 = y2 , ... où (xj , yj , zj ) sont les coordonnées
spatiales de la j eme particule15 . On a alors 3N équations d’Euler-Lagrange :
∂L d ∂L
− =0 pour i = 1, 2, ...3N.
∂qi dt ∂ q̇i
Et le Lagrangien du système s’écrit :
3N
X 1
L= mq̇i2 − U (q1 , q2 , ..., q3N ).
i=1
2

15 Une contrainte sur les coordonnées réduit le nombre de degrés de liberté. Par exemple, une

particule qui se déplace à la surface d’une sphère n’a que deux degrés de liberté. Compte tenu
de la symétrie, on choisira les angles q1 = θ et q2 = φ plutôt que les coordonnées cartésiennes.
En général s’il y a m contraintes, il y a 3N − m degrés de liberté (c’est-à-dire de coordonnées
indépendantes qui fixent la position du système).

23
3.3 Le formalisme hamiltonien
3.3.1 Les équations de Hamilton
Si le Lagrangien ne dépend pas explicitement du temps ( ∂L
∂t = 0) :

3N
X ∂L
q̇i − L (25)
i
∂ q̇i

est une constante du mouvement d’après la formule de Beltrami16 . Autrement

dit, le Lagrangien n’est pas une constante du mouvement.

Dans le cas d’une force conservative (à une dimension), on voit que :

∂L 1 1
q̇ − L = mq̇ 2 − mq̇ 2 + U (q) = mq̇ 2 + U (q)
∂ q̇ 2 2
est une constante du mouvement. Autrement dit, l’énergie du système est bien
conservée. Ce résultat suggère l’utilisation d’une autre fonction que le Lagrang-
ien pour déduire les équations du mouvement. On définit le moment conjugué
(ou impulsion généralisée) pi associé à la coordonnée qi par

∂L
pi = . (26)
∂ q̇i

En général le moment conjugué est l’impulsion pi = mq̇i (cas des N particules

ponctuelles). Mais lorsque la coordonnée généralisée qi est un angle, son moment
conjugué est un moment cinétique. Ainsi, dans le cas du pendule de longueur l
décrit par l’angle θ. Le moment conjugué de θ est pθ = ml2 θ̇, c’est-à-dire son
moment cinétique.

Si le Lagrangien ne dépend pas explicitement de la variable qi on dit que

∂L
cette variable est cyclique ( = 0) et d’après les équations d’Euler-Lagrange
∂qi
∂L
son moment conjugué est une constante du mouvement (pi = = cte). C’est
∂ q̇i
le cas d’une particule libre dont le Lagrangien ne dépend pas de la position.
Ainsi, on relie la notion de symétrie à celle d’invariance.
16 Pour le montrer à nouveau dans le cadre de la mécanique analytique :
3N » –
dL X ∂L ∂L
= q̇i + q̈i
dt i
∂qi ∂ q̇i
3N » –
X d ∂L ∂L d
= ( )q̇i + (q̇i )
i
dt ∂ q̇i ∂ q̇i dt
3N
d X ∂L
= q̇i .
dt i ∂ q̇i

24
 A partir de l’expression (25), définissons une nouvelle fonction appelée Hamil-
tonien et notée H :
3N
X
H= q̇i pi − L = T + U. (27)
i=1

Le Hamiltonien se confond avec l’énergie du système dans le cas de forces con-

servatives et si ∂L
∂t = 0.

Reste à trouver les équations du mouvement à partir de cette nouvelle fonc-

tion. Pour cela différencions H et utilisons l’équation d’Euler-Lagrange :
3N
X  ∂L ∂L
 ∂L
dH = q̇i dpi + pi dq̇i − dqi + dq̇i − dt
i=1
∂qi ∂ q̇i ∂t
3N
X  ∂L  ∂L
= q̇i dpi − dqi − dt
i=1
∂qi ∂t
3N
X   ∂L
= q̇i dpi − ṗi dqi − dt.
i=1
∂t

On voit que les variables naturelles du Hamiltonien sont les coordonnées qi et

leurs moments conjugués pi . On en déduit les équations de Hamilton17 :

∂H ∂H
q̇i = et ṗi = − pour i = 1, 2, ...3N. (28)
∂pi ∂qi

Comme les 3N équations du second ordre d’Euler-Lagrange, les 6N équations du

premier ordre de Hamilton permettent d’exprimer les équations du mouvement.
Par ailleurs,
∂H ∂L
=− .
∂t ∂t
Ainsi, le Hamiltonien d’un système de N particules soumises à des forces con-
servatives s’écrit :
3N

X p2i
H {qi , pi } = + U (q1 , q2 , ..., q3N ).
i=1
2m
17 Une méthode alternative est d’utiliser directement le principe variationnel (δS = 0) :
Z Z
δS = δ (pq̇ − H)dt = (q̇δp + pδq̇ − δH)dt

dδq ∂H ∂H
Z
= (q̇δp + p − δq − δp)dt
dt ∂q ∂p
Z » –
∂H ∂H
= (q̇ − )δp − (ṗ + )δq dt + [pδq],
∂p ∂q
où l’on a intégré par parties. Le terme entre crochets étant nul, on retrouve les équations de
Hamilton.

25
Les équations de Hamilton impliquent pour i = 1, 2, ...3N :
pi ∂U
q˙i = et ṗi = − = Fi ,
m ∂qi
où Fi est la composante de la force associée à la coordonnée qi .

En général, on écrit le Lagrangien du système en fonction des qi et des q̇i ,

on calcule les moments conjugués pi , on en déduit le Hamiltonien en fonction
des qi et des pi , et on résout les équations de Hamilton.

Reprenons l’exemple du pendule de longueur l décrit par l’angle θ (une

dimension). Son Hamiltonien est :

p2
 
m 2 2
H = pθ θ̇ − L = ml2 θ̇2 − l θ̇ + mgl cos θ = θ 2 − mgl cos θ.
2 2ml

Les équations de Hamilton impliquent :

∂H pθ
θ̇ = =
∂pθ ml2
et
∂H
ṗθ = − = −mgl sin θ
∂θ
et on retrouve l’équation du mouvement.

Le portrait de phase

On appelle espace des phases, l’espace à 6N dimensions (pour N particules

qui évoluent dans un espace à trois dimensions) dont les coordonnées sont les
3N coordonnées généralisées qi et les 3N moments conjugués pi . Le portrait de
phase est alors la représentation de la dynamique du système dans l’espace des
phases pour une énergie donnée E.

Dans le cas du pendule simple, l’espace des phases a deux dimensions : q = θ

et p = pθ . Et le portrait de phase est donné par :
p
p = ± 2ml2 (E + mgl cos q).

On distingue deux cas selon la valeur de l’énergie :

• si −mgl < E < mgl : le mouvement n’est possible que pour des angles
compris dans un certain intervalle qui augmente avec l’énergie (libration
ou oscillations libres).
• si mgl < E : le mouvement est possible quelle que soit la valeur de l’angle
(rotation).

26
 Fig. 7 – Le portrait de phase du pendule simple.

3.3.2 Les crochets de Poisson

Considérons une grandeur physique f dépendant des 3N coordonnées qi et
des 3N moments conjugués pi . Sa dérivée totale par rapport au temps est donnée
par
3N
df ∂f X ∂f ∂f

({qi , pi }, t) = + q̇j + ṗj .
dt ∂t j=1 ∂qj ∂pj
D’après les équations de Hamilton,
3N  
df ∂f X ∂f ∂H ∂f ∂H ∂f 
= + − = + f, H , (29)
dt ∂t j=1 ∂qj ∂pj ∂pj ∂qj ∂t

où l’on a introduit les crochets de Poisson de f avec H. En général, les crochets
de Poisson de deux fonctions f et g s’écrivent :
3N  
 X ∂f ∂g ∂f ∂g
f, g = − .
j=1
∂qj ∂pj ∂pj ∂qj

Les équations de Hamilton s’écrivent alors :

 
q̇i = qi , H et ṗi = pi , H pour i = 1, 2, ...3N.
Les crochets de Poisson permettent d’associer des propriétés de symétrie et d’in-
variance. Ainsi, d’après l’équation (29), une quantité
 physique f qui ne dépend
pas explicitement du temps est conservée si f, H = 0. Autrement dit, l’u-
niformité du temps correspond à la conservation de l’énergie. Par ailleurs, les
équations de Hamilton montrent clairement que si H ne dépend pas d’une cer-
taine coordonnée q0 (variable cyclique), son moment conjugué p0 est une con-
stante du mouvement.

Pour ouvrir une porte sur la mécanique quantique, les crochets de Poisson
sont à l’origine de la quantification et du commutateur entre deux opérateurs.
Ainsi, les crochets de Poisson 
q, p = 1,

27
annonce la relation fondatrice de la mécanique quantique18 :

[q, p] = i~,

où ~ = h/2π, h est la constante de Planck.

3.4 Conclusion
Pour traiter un problème de mécanique, nous avons trois approches à notre
disposition :
– Newtonienne : 3N équations (scalaires) différentielles du 2eme ordre à
résoudre (et 6N conditions initiales). Mais la notion de force sur laquelle
repose la mécanique newtonienne ne se généralise pas à l’échelle quantique.
– Lagrangienne : 3N équations différentielles du 2eme ordre à résoudre (et
6N conditions initiales).
– Hamiltonienne : 6N équations différentielles du 1eme ordre à résoudre (et
6N conditions initiales)19 .
La différence principale entre les approches newtonienne et lagrangienne est que
la première est différentielle et donc locale (relation entre la cause et l’effet). Le
principe de moindre action implique lui une approche variationnelle et donc
globale20 .

Les deux dernières approches qui constituent la mécanique analytique ont

les avantages suivants :
– En mécanique analytique on manipule des coordonnées généralisées qui
sont des scalaires et non pas des vecteurs.
– La mécanique analytique se généralise dans le cadre de la mécanique quan-
tique et de la théorie des champs, ce qui a permis le développement de la
physique des particules et des interactions fondamentales.

18 On consultera avec profit Mécanique quantique, 1. Fondements et premières applications,

C. Aslangul (de Boeck, Bruxelles, 2007).

19 En général, il est plus facile d’écrire un Lagrangien qu’un Hamiltonien à l’aide des

symétries du système étudié. Quitte à utiliser ensuite le formalisme Hamiltonien. Mais, il

n’est pas toujours possible de déduire le moment conjugué du Lagrangien et donc d’écrire le
Hamiltonien d’un système (par exemple dans le cas de l’électrodynamique classique).
20 On retrouve cette distinction en mécanique quantique entre la formulation différentielle

de Schrödinger et celle en terme d’intégrales de chemin de Feynman.

28
3.5 Exercices
3.5.1 L’oscillateur harmonique
L’énergie potentielle d’un oscillateur harmonique (classique) à une dimension
s’écrit
1
U (q) = mω 2 q 2 .
2
1) En déduire l’équation du mouvement à l’aide du formalisme lagrangien.
2) En déduire l’équation du mouvement à l’aide du formalisme hamiltonien.
3) Montrer à l’aide des équations de Hamilton que le Lagrangien suivant
décrit le mouvement d’un oscillateur amorti :
1 αt 2
L= e (q̇ − ω 2 q 2 ).
2

Réponse :

1) Le Lagrangien s’écrit :
1 2 1
L= mq̇ − mω 2 q 2 .
2 2
Et d’après l’équation d’Euler-Lagrange :
∂L d ∂L
=
∂q dt ∂ q̇
donc
−mω 2 q = mq̈.

2) Le Hamiltonien s’écrit :

p2 1
H= + mω 2 q 2 .
2m 2
Les équations de Hamilton donnent
p
q̇ = et ṗ = −mω 2 q.
m
On en déduit l’équation du mouvement :

q̈ + ω 2 q = 0.

∂L
3) Le moment conjugué s’écrit p = ∂ q̇ = q̇eαt donc

1 2 −αt 1 2 2 αt
H= p e + ω q e .
2 2

29
Les équations de Hamilton donnent

q̇ = pe−αt et ṗ = −ω 2 qeαt .

On en déduit l’équation du mouvement :

q̈ + ω 2 q + αq̇ = 0.

3.5.2 Le cerceau à vitesse angulaire constante

Une bille glisse sans frottement le long d’un cerceau de rayon R animé d’un
mouvement de rotation autour de son axe à vitesse angulaire ω constante. La
position de la bille sur le cerceau est déterminée par un seul degré de liberté,
l’angle θ avec la verticale. La position du cerceau est elle décrite par l’angle φ
(voir la figure 8).

Fig. 8 – La bille sur le cerceau en rotation.

1) Écrire le Lagrangien de ce système et en déduire l’équation du mouvement.

2) Écrire le Hamiltonien du système et en déduire à nouveau l’équation du
mouvement.
3) Montrer que la bille est soumise à un potentiel effectif p V (θ) dont on
étudiera le comportement en fonction de θ. On posera ω0 = g/R.
4) Quel est le comportement de la bille pour des vitesses de rotation lente
(ω < ω0 ) ?
5) Même question dans le cas de la rotation rapide (ω > ω0 ).

Réponse :

1) En coordonnées sphériques, la position de la bille est donc déterminée par

R, θ(t) et φ(t). Sa vitesse est donc Rθ̇~eθ + R sin θφ̇~eφ . Son énergie potentielle

30
est U (θ) = mgR cos θ. On en déduit le Lagrangien du système :
1
L(θ, θ̇) = mR2 (θ̇2 + ω 2 sin2 θ) − mgR cos θ.
2
On en déduit à l’aide de l’équation d’Euler-Lagrange :
∂L d ∂L
= ,
∂θ dt ∂ θ̇
soit
d
mR2 ω 2 sin θ cos θ + mgR sin θ = (mR2 θ̇) = mR2 θ̈,
dt
donc
mRω 2 sin θ cos θ + mg sin θ = mRθ̈.
A l’aide de l’équation fondamentale de la dynamique projetée sur ~eθ pour
éliminer la réaction du cerceau sur la bille et en tenant compte de la force
centrifuge de norme mω 2 R sin θ on aurait bien :

mRθ̈ = mg sin θ + (mω 2 R sin θ) cos θ.

2) Le moment conjugué de θ s’écrit :

∂L
pθ = = mR2 θ̇,
∂ θ̇
donc
1 p2θ 1
H = θ̇pθ −L = mR2 (θ̇2 −ω 2 sin2 θ)+mgR cos θ = 2
− mR2 ω 2 sin2 θ+mgR cos θ.
2 2mR 2
D’après les équations de Hamilton :
∂H pθ
θ̇ = =
∂pθ mR2
et
∂H
p˙θ = − = mR2 ω 2 sin θ cos θ + mgR sin θ.
∂θ
On retrouve bien sûr l’équation du mouvement.

3) D’après le Hamiltonien :
1
V (θ) = mR2 (2ω02 cos θ − ω 2 sin2 θ).
2
La fonction est symétrique par rapport à π (V (π − θ) = V (π + θ)). Et

dV
= −mR2 sin θ(ω02 + ω 2 cos θ).
dθ

31
 2
pour w=0.7w0
pour w=w0
pour w=1.3w0
1.5

0.5

V(theta)
-0.5

-1

-1.5

-2

-2.5
0 0.5 1 1.5 2 2.5 3
theta

Fig. 9 – Le potentiel V (θ) entre 0 et π.

Entre 0 et π, la dérivée est négative si ω02 + ω 2 cos θ > 0, donc si

ω02
cos(θ) > − .
ω2

4) Dans le cas d’une rotation lente (ω < ω0 ), cette condition est remplie
quelle que soit la valeur de θ. La dérivée est donc négative pour tout θ et
s’annule en 0 (équilibre instable) et en π (équilibre stable), en bas du cerceau.
Au voisinage de ce q point d’équilibre, la bille a un mouvement harmonique avec
ω02
la pulsation ωl = ω ω2 − 1. En effet :

d2 V
= −mR2 (ω02 cos θ + ω 2 (cos2 θ − sin2 θ))
dθ2
et
1 d2 V 

1
2
= mR2 ωl2 .
2 dθ π  2

5) Dans le cas d’une rotation rapide (ω > ω0 ), cette condition est remplie si
ω2
θ < θs tel que cos(θs ) = − ω02 . L’angle θs est donc un point d”équilibre (stable)
et 0 et π deviennent des points d’équilibre instables. En ce point, l’énergie
potentielle vaut
1 ω4
V (θs ) = − mR2 ω 2 (1 + 04 ).
2 ω
Laqbille a un mouvement oscillant autour de l’angle θs avec la pulsation ωr =
ω04
ω 1− ω4 . Quand ω → ∞, θs → π/2.

32
3.5.3 L’atome hydrogénoı̈de
On considère l’électron d’un atome hydrogénoide21 subissant l’effet du po-
tentiel central :
Ze2
V (r) = − . (30)
4πǫ0 r
1) Donner l’expression du Lagrangien de ce système en coordonnées sphériques.
2) Établir les équations de Lagrange dans ce même système de coordonnées.
3) En déduire l’existence de quantités conservées. Interpréter.
4) Écrire le Hamiltonien de ce système et les équations de Hamilton.

Réponse :

1 2 1 Ze2 1 Z ẽ2
1) On a L = mṙ −V (r) = mṙ2 + ˆ mṙ2 +
= dont on note qu’il
2 2 4πǫ0 r 2 r
est invariant par rotation. On a introduit une charge effective pour simplifier les
notations
En coordonnées sphériques on a : r = r ur et ṙ = ṙ ur + r θ̇ uθ + r sin θ ϕ̇ uϕ
d’où :
1 Z ẽ2
L(r, θ, ϕ) = m(ṙ2 + r2 θ̇2 + r2 sin2 θ ϕ̇2 ) + . (31)
2 r

2) D’après l’équation de Lagrange :

 
d ∂L ∂L
= , (32)
dt ∂ q̇i ∂qi
on a, pour les coordonnées r, θ et ϕ :

d   Z ẽ2
2 2 2

 (m ṙ) = m r θ̇ + r sin θ ϕ̇ − 2
dt r






d  2 

mr θ̇ = mr2 sin θ cos θ ϕ̇2 (33)


 dt


 d

mr2 sin2 θ ϕ̇ = 0.



dt

3) On déduit immédiatement des équations précédentes que la quantité

mr2 sin2 θ ϕ̇ est conservée. La conservation de cette grandeur est liée au fait
que le Lagrangien est indépendant de la variable angulaire ϕ. Elle traduit donc
l’invariance par rotation du système autour de l’axe Oz. Vérifions que la quantité
conservée est la projection du moment angulaire sur l’axe Oz.
On a J = r ∧ mṙ et Jz = J.uz où uz est donné par : uz = cos θ ur − sin θ uθ
d’où : Jz = mr2 sin2 θϕ̇.
21 Un atome hydrogénoı̈de est un atome auquel on a arraché tous ses électrons sauf un. Il se

comporte donc comme un atome d’hydrogène avec un noyau de charge Z > 1.

33
 L’invariance complète par rotation du système implique que les autres com-
posantes Jx et Jy sont aussi conservées. C’est la dissymétrie du système de
coordonnées qui privilégie l’axe Oz et, par conséquent, la composante Jz du
moment cinétique. Mais comme le choix de l’axe Oz est arbitraire la conserva-
tion de Ji vaut pour tout i = x, y ou z.

4) Calculons d’abord les moments conjugués :

∂L
pr = = mṙ
∂ ṙ
∂L
pθ = = mr2 θ̇
∂ θ̇
∂L
pϕ = = mr2 sin2 θϕ̇.
∂ ϕ̇
En coordonnées sphériques, le moment conjugué est donc différent de l’impulsion
p = mṙ. Le Hamiltonien du système s’écrit donc

H = ṙpr + θ̇pθ + ϕ̇pϕ − L

1 Z ẽ2

= mṙ2 + mr2 θ̇2 + mr2 sin2 θϕ̇2 − m(ṙ2 + r2 θ̇2 + r2 sin2 θ ϕ̇2 ) +
2 r
1 2 2 2 2 2 2 Z ẽ2
= m(ṙ + r θ̇ + r sin θ ϕ̇ ) −
2 r
p2r p2θ p2ϕ Z ẽ2
= + + − .
2m 2mr2 2mr2 sin2 θ r
Les équations de Hamilton sont :

∂H pr ∂H p2 p2ϕ Z ẽ2
ṙ = = et ṗr = − = θ3 + −
∂pr m ∂r mr mr3 sin2 θ r2
2
∂H pθ ∂H pϕ cos θ
θ̇ = = et ṗθ = − =
∂pθ mr 2 ∂θ mr2 sin3 θ
∂H pϕ ∂H
ϕ̇ = = et ṗϕ = − = 0.
∂pϕ mr2 sin2 θ ∂ϕ

En utilisant les expressions des moments conjugués, on retrouve bien les équations
du mouvement obtenues à l’aide du formalisme lagrangien.

3.5.4 La particule plongée dans un champ électromagnétique

On considère le Lagrangien d’une particule de masse m et de charge q plongée
dans un champ électromagnétique (E, B) décrit par le potentiel scalaire φ(r, t)
et vectoriel A(r, t) :
1
L= mv2 + q v.A(r, t) − q φ(r, t) . (34)
2
1) Déterminer l’équation du mouvement de la particule.

34
 2) Déterminer le Hamiltonien de la particule.
3) En déduire les équations de Hamilton.

Réponse :

1) On a :
1
L= mv2 + q v.A(r, t) − q φ(r, t) (35)
2
où l’on note que A et φ font intervenir, via leurs arguments, la position de la
particule.
On en tire :
∂L


 = mẋi + qAi (r, t)
 ∂ ẋi

(36)

 ∂L ∂A j (r, t) ∂φ(r, t)
=q ẋj − q


∂xi ∂xi ∂xi
et :    
d ∂L ∂Ai (r, t) ∂Ai (r, t)
= mẍi + q + ẋj (37)
dt ∂ ẋi ∂t ∂xj
avec sommation sous-entendue sur l’indice répété (j). L’équation du mouvement
de la particule est alors donnée par :
 
∂φ(r, t) ∂Ai (r, t) ∂Ai (r, t) ∂Aj (r, t)
mẍi = −q −q + ẋj − ẋj . (38)
∂xi ∂t ∂xj ∂xi
On montre à la fin de cet exercice que l’on retrouve bien l’équation usuelle de
la particule chargée dans un champ électromagnétique :
dv
m = q (E + v ∧ B) . (39)
dt

2) On détermine d’abord le moment conjugué pi . On a :

∂L
pi = = mẋi + qAi (r, t) (40)
∂ ẋi
d’où :
pi qAi (r, t)
ẋi (xj , pj ) = − . (41)
m m
On remarque que le moment conjugué est différent de l’impulsion pi = mẋi . Le
Hamiltonien est donné par :
X
H = pi ẋi − L(xi , ẋi (xj , pj ), t)
i

2
p2
 X  pi 
q X 1 X pi qAi qAi
= − pi Ai − m − −q Ai − + qφ
m m i 2 i
m m i
m m
(42)

35
où l’on a, cette fois, explicité les sommes (sur l’indice i). On a, en simplifiant
l’expression précédente :
1
H= (p − qA)2 + q φ (43)
2m
où l’on remarque que l’impulsion de la particule p est remplacée par p − qA.

∂H ∂H
3) Les équations de Hamilton sont données par : ṗi = − et q̇i =
∂xi ∂pi
soit, de façon explicite :
 
∂H ∂ 1 2 q ∂Ak ∂φ
= (p − qA(r, t)) + q φ(r, t) = − (pk − qAk ) +q
∂xi ∂xi 2m m ∂xi ∂xi
(44)
et
∂H 1
= (pi − qAi ) . (45)
∂pi m

Montrons à présent que l’équation du mouvement (38) est bien l’équation usuelle
de la particule chargée dans un champ électromagnétique :
dv
m = q (E + v ∧ B) . (46)
dt
Le champ électrique E vérifie :
∂B ∂
∇∧E = − = − (∇ ∧ A) (47)
∂t ∂t
A étant le potentiel vecteur.
On a donc :
∂A(r, t)
E=− − ∇φ(r, t) (48)
∂t
et la composante Ei du champ électrique est donnée par :
∂Ai (r, t) ∂φ(r, t)
Ei = − − (49)
∂t ∂xi
qui correspond bien (à la charge q près) aux deux premiers termes du second
membre de l’expression (38). La composante i du produit vectoriel v ∧ B est
donnée par :
(v ∧ B)i = ǫijk ẋj Bk (50)
où, de nouveau, la sommation sur les indices répétés (j et k) est sous-entendue.
Le tenseur antisymétrique de Levi-Civita ǫijk est définit par : ǫ123 = 1 ainsi
que toute permutation circulaire de ces indices : ǫ123 = ǫ231 = ǫ312 = 1. L’an-
tisymétrie implique : ǫ213 = ǫ321 = ǫ132 = −1 et ǫiik = ǫijj = ǫiii =etc= 0.

36
 Finalement une propriété importante de ce tenseur est :
X
ǫijk ǫilm = ǫijk ǫilm = (δjl δkm − δjm δkl ) (51)
i

δij étant le symbole de Kronecker. Comme B est un rotationnel B = ∇ ∧ A sa

composante Bk est donnée par :

Bk = ǫklm ∂l Am . (52)

On a donc :
(v ∧ B)i = ǫijk ǫklm ẋj ∂l Am = ǫkij ǫklm ẋj ∂l Am

= (δil δjm − δim δjl ) ẋj ∂l Am = ẋj ∂i Aj − ẋj ∂j Ai

(53)
 
∂Ai ∂Aj
=− ẋj − ẋj
∂xj ∂xi

qui sont bien les troisième et quatrième termes figurant dans le second membre
l’expression (38).

37