Sol 1
Sol 1
Sol 1
TD de Mécanique Analytique
Corrigé de la série N˚ 1
Parcours Physique Moderne : Equations de Lagrange
Corrigé de l’exercice 1
RO l mg l RB
θ
O x
F B
1. Voir cours. Quand le système est statique ou il se déplace d’un mouvement uniforme, le travail de
toutes les forces est nul, pas seulement des forces intérieures car la résultante des forces extérieures
est nulle.
2. (a) La masse m se déplace dans un plan. Comme, OM = l, alors elle a un seul dégré de libérté.
Le mouvement de m peut être bien repéré par la variable θ, qui sera utilisée comme coordonnée
généralisée.
(b) On dénombre quatre forces : les réactions normales, puisqu’il n’y a pas de frottement, aux points
O et B R ~ o et R
~ B , le poids m~g et la force F
~ appliquée au point B.
Le principe de d’Alembert stipule que le travail des forces intérieures lors d’un déplacement
virtuel est nul. Considérons le déplacement virtuel δθ et calculons la force généralisée selon cette
coordonnée :
X
Qθ = ~i ∂~ri
F
i
∂θ
où ~ri est le vecteur qui repère le point d’application de la force. Ce qui donne, en utilisant la base
cartésienne
−→
OA = lcosθ~i + lsinθ~j
−−→
OB = 2lcosθ~i
−→ −−→
∂ OA ~ ∂ OB
Qθ = −mg~j +F
∂θ ∂θ
= mglcosθ − 2lF sinθ
m~g = −mg~j et F~ = −F~i et sachant que R0 et RB ne travaillent pas le déplacement leur est
perpendiculaire. Or le principe de d’Alembert donne
mg
Qθ δθ = 0 =⇒ F = cotθ
2
donc pour cette valeur, le système sera statique.
1
(c) Pour déterminer la valeur de la réaction en B R ~ B , il suffit de prendre comme déplacement virtuel
du point B un cerle de rayon 2lcosθ avec θ constant. Ainsi, F ~ ne travaille pas, de même pour
~ ~
R0 . Les deux forces qui travaillent sont le poids et RB . Dans cette configuration, les vecteurs
(~i, ~j) deviennent mobiles. Soit R0 le référentiel par rapport auquel on effectue cette rotation
et soit (~i0 , ~j0 , ~k0 ) la base orthonormée liée à R0 . Repérons la rotation par l’angle ψ telle que
d d
ψ = (~i,~i0 ) = (~j, ~j0 ), ce qui donne d~i/dψ = ~j et d~j/dψ = −~i et ψ̇ = constante.
Calculons la coordonnée généralisée associée à ψ
−→ −−→
∂ OA ∂ OB
Qψ = · m~g + · (RB~j)
∂ψ ∂ψ
= l(cosθ~j − sinθ~i) · (−mg~j0 ) + 2lcosθ~j · (RB~j)
= −mgl(cosθcosψ − sinθsinψ) + 2lRB cosθ
= −mglcos(θ − ψ) + 2lRB cosθ.
Or comme ψ̈ = 0, cela implique que le module de la vitesse par rapport à R0 de tous les points
du treillis est constant et comme le mouvement est circulaire alors l’accélération par rapport à
V2 VB2
R0 est centrale ~γA = lA ~nA et ~γB = 2lcosθ ~nB ce qui implique que l’accélération généralisée selon
ψ est donnée par
−→ −−→
∂ OA ∂ OB
Aψ = · ~γA + · ~γB = 0
∂ψ ∂ψ
−→ −→ −−→ −−→ −→ −−→
puisque ~nA = OA/kOAk et ~nB = OB/kOBk, d’une part, et OA ⊥ ~nA et OB ⊥ ~nB , d’autre part.
Le principe de d’Alembert Qψ δψ = Aψ δψ = 0 permet d’écrire
mgcos(θ − ψ)
Qψ = 0 =⇒ RB =
2cosθ
et cette relation est valable quelque soit la valeur de ψ et donc en particulier pour ψ = 0 qui nous
ramène à la situation du treillis statique
mg
RB = .
2
1 Corrigé de l’exercice 2
ur
θ M
uϕ
r Y
O uθ
ρ
ϕ
En raison de la symétrie sphérique du problème, nous allons utiliser les coordonnées sphériques (r, θ, ϕ).
Soit (~er , ~eθ , ~eϕ ) la base sphérique. Rappelons que les vecteurs de la base sont en rotation dans le repère
2
~ = θ̇~eϕ + ϕ̇~k, ~k étant le vecteur selon Oz de la base
R(O, xyz), qui est galiléen, avec le vecteur de rotation Ω
cartésienne. Rappelons aussi que
d~ei ~ ∧ ~ei
= Ω
dt
~ei étant l’un des vecteurs de la base sphérique.
1. La position de m est repérée dans le système de coordonnées sphériques par trois coordonnées,
−−→
OM = r. Comme, m doit se déplacer à l’intérieur de la sphère, alors le nombre de degré de liberté
de m est 2. Ainsi seules θ et ϕ décrivent complètement le mouvement de m.
2. La vitesse de m est
d
~v = (a~er )
dt
= aΩ~ ∧ ~er
= a(θ̇~eθ + ϕ̇sinθ~eϕ )
et
−−→ −−→
∂ OM ∂ OM
Qϕ = (−mg~k) · + (−N~er ) ·
∂ϕ ∂ϕ
∂~er
= −mg(a~k − N~er ) ·
∂ϕ
= −mg(a~k − N~er ) · (~k ∧ ~er )
= −mgsinθ(r~k − N~er ) · ~eϕ
= 0
3
4. Si θ et ϕ̇ = Ω sont constants, les équations du mouvement deviennent
g
ϕ̇2 =
acosθ
qui n’est valide que si cosθ > 0 =⇒ −π/2 < θ < π/2, ce qui est le cas. Ainsi les équations du
mouvement sont
r
g
ϕ̇ =
acosθ
θ = constante.
Exercice 3
Soit une perle de masse m qui coulisse sans frottement sur un cerceau de rayon R, figure 3.
φ M
θ R
O y
1. En raison de la symétrie sphérique du problème, nous utilisons les coordonnées sphériques pour décrire
la position de la perle.
La perle se déplace sur le cerceau alors p
r = R et comme son mouvement dans le plas Oxy est circulaire
de rayon Rsinθ alors ρ = Rsinθ ; ρ = x2 + y 2 . Ainsi, les deux contraintes sur le mouvement de la
perle réduisent le nombre de degrés de liberté à 1. On utilise donc θ comme coordonnée généralisée
du mouvement de la perle.
2. Calculons l’énergie cinétique,
−−→ ~ d~er
OM = R~er =⇒ V = R
dt
~
= RΩ ∧ ~er
= R(ϕ̇~k + θ̇~eϕ ) ∧ ~er
= R ϕ̇sinθ~eϕ + θ̇~eθ
Calculons l’énergie potentielle. Pour ce faire calculons le travail des forces appliquées à la perle
qui sont son poids m~g = −mg~k et R ~ = −RN ~er . La réaction normale ne travaille pas puisque
~
δW (RN ) = −RN ~er · d~er = 0 puisque ~er ⊥ d~er . Le déplacement élémentaire est
→
− −−→
dl = dM = Rd~er = RdtΩ ∧ ~er = R dϕ~k + dθ~eϕ ∧ ~er = Rdϕsinθ~eϕ + Rdθ~eθ
4
ce qui donne pour le travail élémentaire du poids
Rappelons que ~k · ~eϕ = 0 et on prend la constante du potentiel égale à 0 puisqu’elle n’a pas d’effet
sur le mouvement de la perle.
Le lagrangien du système peut s’exprimer ainsi comme
1
L = T −V = mR2 Ω2 sin2 θ + θ̇ 2 + mgRcosθ.
2
3. Calculons le moment conjugué p de la perle :
∂L
p = pθ = = mR2 θ̇.
∂ θ̇
On note que le moment conjugué n’est d’autre que le moment cinétique de la perle, ce qui est prévu
étant donné que la coordonnée généralisée est un angle. On peut écrire p = I θ̇ avec I = mR2 le
moment d’inertie de la perle par rapport à Oz.
L’expression du Hamiltonien est
H = pθ̇ − L
1
= I θ̇ 2 − I Ω2 sin2 θ + θ̇ 2 − mgRcosθ
2
1 2 1 2 2
= I θ̇ − IΩ sin θ − mgRcosθ
2 2
p2 1 2 2
= − IΩ sin θ − mgRcosθ
2I 2
2
où le terme p2I peut être identifié à une énergie cinétique de rotation de la perle autour de ~eϕ , ce qui
nous permet d’interpréter le terme restant comme une énergie potentielle effective Ũ (θ). Aussi
p2 1
H = + Ũ (θ) avec Ũ (θ) = − IΩ2 sin2 θ − mgRcosθ
2I 2
4. Cherchons les extremums de Ũ (θ) :
dŨ (θ)
= −IΩ2 sinθcosθ + mgRsinθ
dθ
= 0 =⇒ −sinθ IΩ2 cosθ − mgR = 0
(
sinθ0 = 0 =⇒ θ0 = 0 puisque θ ∈ [−π/2, π/2]
=⇒
cosθ1 = mgR g g
IΩ2 = RΩ2 =⇒ θ1 = ±arccos RΩ2 .
d2 Ũ (θ)
= −IΩ2 cos(2θ) + mgRcosθ
dθ 2
5
on en déduit que
d2 Ũ (θ) 2 2 2 Ωc
2
= mgR − IΩ = mR Ω ( − 1) < 0 =⇒ θ0 = 0 est un maximum
dθ 2 θ0 =0 Ω2
d2 Ũ (θ)
= −2IΩ2 cos2 θ1 + IΩ2 + mgRcosθ1
dθ 2 θ1
= cosθ1 (IΩ2 cosθ1 − 1) + IΩ2 (1 − cos2 θ)
= IΩ2 (1 − cos2 θ) > 0 =⇒ θ1 est un minimum
Pour retrouver la trajectoire de la perle si les conditions initiales sont θ(t = 0) = 0 et θ̇(t = 0) = 0, il suffit
de résoudre l’équation du mouvement autour de θ = 0.
L’équation du mouvement est obtenue en utilisant soit les équations de Lagrange soit les équations de
Hamilton. En faisant usage de ces dernières, on
∂H p
θ̇ = =⇒ θ̇ =
∂p I
∂H 1
ṗ = − =⇒ ṗ = − IΩ2 sin2θ + mgRsinθ
∂θ 2
or en utilisant la première équation on a ṗ = I θ̈, ce qui donne finalement
1
θ̈ + Ω2 sin2θ − Ω2c sinθ = 0
2
où nous avons utilisé Ωc = g/R et I = mR2 . Comme θ → 0 alors sin2θ = 2θ + O(θ 2 ) et sinθ = θ + O(θ 2 ),
2
ce qui donne
θ̈ + (Ω2 − Ω2c )θ = 0
qui est une équation différentielle de second ordre sans second membre
p à coefficients constants dont l’équation
caractérestique 2 2 2 2 2
√ est
2
r +(Ω −Ω
2
√ c2 ) =2 0 ayant comme solutions r = ±i Ω − Ωc car Ω > ΩC . Aussi, la solution
est θ = Ae−i Ω p −Ωc t
+ Bei Ω −Ωc t A et B étant déterminées par les conditions initiales θ(0) = 0 = A + B
et θ̇(0) = 0 = −i Ω − Ω2c (A − B) ce qui donne A = B = 0 et donc la solution se réduit à θ = 0. Ce résultat
2
est bien prévisible puisque la position θ = 0 est une position d’équilibre donc si θ̇ = 0, l’anneau campe sur
sa position d’équilibre.
Corrigé de l’exercice 4
Une sphère S1 pleine de rayon a roule sans glisser sur une deuxième sphère S2 de rayon b. A l’instant t = 0,
S1 est située au sommet de S2 , voir figure ci-dessous. On note par R(O, xyz) muni de la base orthonormée
directe cartésienne (~i, ~j, ~k) et soit C le centre de S1 . I est le point de contact de S1 avec S2 .
z (S1)
RN
(S2) RT ψ
C
a
θ I
b mg
j eρ
O x
i
eϕ
6
1. Cas où S1 reste en contanct avec S2 :
(a) Nous avons deux méthodes pour établir la condition de roulement sans glissement.
(i) Considérons un point M confondu avec le point de contact entre les deux sphères à l’instant
initial. Nous avons roulement sans glissement si la distance parcourue par M à l’instant t est
égale à celle parcourue par I ; ainsi
aψ = bθ.
(ii) La deuxième approche consiste à calculer la vitesse du point de contact I et d’exprimer qu’il
y a roulement sans glissement et donc V ~ (I ∈ S1 ) = V~ (I ∈ S2 ), ce qui donne sachant que S2
est immobile
~ (I ∈ S1 ) = 0
V
−→
Pour la vitesse de I, on décompose OI en passant par C pour faire apparaı̂tre le mouvement
−→ −−→ −→
de I en tant que point de S1 , OI = OC + CI sachant que
−−→
dOC ~c = θ̇~j ∧ −
−→
= V OC = (a + b)θ̇~k ∧ ~eρ
dt R
−→
dCI
= −a(θ̇ + ψ̇)~k ∧ ~eρ
dt R
−−→ −−→
avec ~eρ = OC/kOCk. Comme
−−→ −→
~ (I ∈ S1 ) = dOC dCI
V +
dt R dt R
= (a + b)θ̇ − a(θ̇ + ψ̇) ~k ∧ ~eρ
= 0 =⇒ bθ̇ = aψ̇ =⇒ bθ = aψ
puisque θ(0) = ψ(0) = 0. On retrouve ainsi le résultat établi dans la première approche.
(b) Calculons l’énergie cinétique de S1 en utilisant le théorème de Koenig, sachant que le vecteur
~ 1 /R) = (θ̇ + ψ̇)~j et que le moment d’inertie de la sphère par rapport à son centre est
rotation Ω(S
I = 2/5 ma2 ,
IΩ2 1 ma2 1
T = + mVc2 = (θ̇ + ψ̇)2 + m(a + b)θ̇ 2
2 2 5 2
Pour établir l’expression de l’énergie potentielle de S1 , remarquons que les forces qui travaillent
lors d’un déplacement élémentaire sont le poids m~g et R ~ T alors que R
~ N ne travaille pas. On
calculera ainsi le lagrangien en incluant l’énergie potentielle du poids alors que R ~ T interviendra
dans les équations de Lagrange par ses composantes généralisées.
→
−
Calculons alors le travail de m~g = −mg~k sachant que dl = (a+b)d~eρ avec d~eρ = dθ −sinθ~k + cosθ~i
→
−
δW (m~g ) = −mg~k · dl
= −mg(a + b)~k · d~eρ
= −mg(a + b)dθ~k · −sinθ~k + cosθ~i
= +mg(a + b)dθsinθ =⇒ dV = mg(a + b)dθsinθ =⇒ V = −mg(a + b)cosθ + K(= 0)
7
que θ et ψ sont inter-dépendante.
Le lagrangien est ainsi donné par
(d) Les équations du mouvement sont données par les équations de Lagrange en tenant compte de
R~ T par ses coordonnées généralisées Qθ et Qψ
d ∂L ∂L
− = Qθ
dt ∂ θ̇ ∂θ
d ∂L ∂L
− = Qψ
dt ∂ ψ̇ ∂ψ
avec
−→
∂ OI ~
Qθ = · RT = ~eϕ · (−RT ~eϕ ) = −RT
∂θ
−→
∂ OI ~
Qψ = · RT = 0
∂ψ
et
∂L
∂θ = mg(a + b)sinθ
2
∂L
∂ θ̇
= 2ma 2
5 (θ̇ + ψ̇) + m(a + b) θ̇
∂L
∂ψ = 0
∂L
= 2ma
2
∂ ψ̇ 5 (θ̇ + ψ̇)
2ma2
(θ̈ + ψ̈) + m(a + b)2 θ̈ − mg(a + b)sinθ = −RT
5
2ma2
(θ̈ + ψ̈) = 0
5
ce qui donne
2ma2
(θ̈ + ψ̈) + m(a + b)2 θ̈ − mg(a + b)sinθ = bλ
5
2ma2
(θ̈ + ψ̈) = −aλ
5
bθ̇ − aψ̇ = 0
8
en calculant λ à partir de la deuxième équation, λ = −2ma/5 (θ̇ + ψ̇), et en substituant son
expression dans la première équation, on obtient
2ma2 2mab
(θ̈ + ψ̈) + m(a + b)2 θ̈ − mg(a + b)sinθ = − (θ̈ + ψ̈)
5 5
2ma
=⇒ (a + b)(θ̈ + ψ̈) + m(a + b)2 θ̈ − mg(a + b)sinθ = 0
5
On en déduit que
b
ψ̈ = θ̈
a
5bg
= sin(π − ψ)
7a(a + b)
5bg
= sin(ψ)
7a(a + b)
5g
θ̈ = sinθ
7(a + b)
5bg
ψ̈ = sin(ψ)
7a(a + b)
2mg
RT = b sinθ
7
2. Cas où S1 peut perdre le contact avec S2 :
(a) Comme la sphère S1 peut quitter S2 la distance entre O et C peut varier et diffère de a + b
au moment où S1 quitte S2 alors qu’en amont de ce moment les équations du mouvement sont
identiques à celles établies dans la question précédente avec la contrainte
−−→
kOCk = r = a + b
ma2 1
T = (θ̇ + ψ̇)2 + m(ṙ 2 + r 2 θ̇ 2 ).
5 2
alors que l’énergie potentielle reste la même. Le lagrangien est alors
9
Seule l’énergie potentielle du poids est inclue dans l’expression du Lagrangien. Comme nous
avons cette fois-ci trois coordonnées et deux contraintes les équations de Lagrange sont obtenues
en utilisant deux multiplicateurs de Lagrange λ1 et λ2 comme suit
−→
d
dt
∂L
∂ ṙ − ∂L
∂r = a1r λ1 + a2r λ2 = ∂ OI
∂r · (−RT ~
eϕ + RN ~eρ ) = RN
−→
d ∂L
− ∂L
= a1θ λ1 + a2θ λ2 = ∂ OI
∂θ · (−RT ~
eϕ + RN ~eρ ) = −rRT
dt ∂ θ̇ ∂θ
−→
OI · (−R ~e + R ~e ) = 0
d ∂L
− ∂L
= a1ψ λ1 + a2ψ λ2 = ∂∂ψ T ϕ N ρ
dt ∂ ψ̇ ∂ψ
avec δr = 0 =⇒ a1r = 1 et a1θ = a1ψ = 0 et bδθ − aδψ = 0 =⇒ a2r = 0, a2θ = b et a2ψ = −a. Or
∂L
= mr θ̇ 2 − mgcosθ
∂r
∂L
= mṙ
∂ ṙ
ce qui donne pour l’ensemble des équations, en utilisant les resultats de la question précédente
mr̈ − mr θ̇ 2 + mgcosθ = λ1
2ma2
(θ̈ + ψ̈) + mr 2 θ̈ − mgrsinθ = bλ2
5
2ma2
(θ̈ + ψ̈) = −aλ2
5
r = a+b
bθ − aψ = 0
10g g 17 10
(1 − cosθ) − cosθ = 0 =⇒ cosθ =
7(a + b) a+b 7 7
10
cosθ =
17
10
θ = arccos( ) = 54˚.
17
Corrigé de l’exercice 5
Dans un espace à deux dimensions (x, z), on considère un milieu matériel d’indice de réfraction n=n(z).
La distance parcourue ds est liée à l’indice de réfraction par ds = cdt/n, où c est la vitesse de la lumière
dans le vide, voir figure ci-dessous.
10
Z
θ
z+dz ds= dx2+dz2
z
O
x x+dx X
1. Comme indiqué ci-dessus, si l’on prend un élément de distance dans le plan (x, z) dans un milieu
d’indice de réfraction n = n(z), nous avons
p c n(z) p 2
ds = dx2 + dz 2 = vdt = dt =⇒ dt = dx + dz 2 .
n(z) c
Si l’on définit la quantité
s
Z Z Z Z
B B n(z) p 2 B n(z) dx2 B n(z) p
S = dt = dx + dz 2 = 1 + 2 dz = 1 + x′2 dz
A A c A c dz A c
où x′ = dx/dz et minimiser le chemin optique consiste à minimiser le temps que mettra la lumière
pour parcourir un chemin d’un point A à un point B et donc cela reviendra à minimiser la quantité
S, qui jouera le rôle d’une action. Nous pouvons déduire alors que la fonctionnelle
n(z) √
L = 1 + x′ 2 = L(x′ , z)
c
joue le rôle d’un lagrangien. Notons que z joue le rôle du temps dans le formulation de Lagrange. On
note que L ne dépend pas de x et donc cette dernière est une variable cyclique ce qui implique que
son moment conjugué est une intégrale première
11
3. Trouver la trajectoire du rayon de lumière consiste à trouver l’équation z = z(x). On a le choix de
partir soit du résultat de la première question ou de celui de la deuxième question.
Prenons le second résultat,
n(z)
H= √ =K
c 1 + z ′2
et la constante K est déterminée à partir de z ′ (x = 0) = 0 et z(x = 0) = 0 ce qui donne, sachant que
n(z) = n0 + λz =⇒ n(0) = n0
n0
K = .
c
Ainsi, nous pouvons écrire
n0 + λz n0
√ =
c 1 + z ′2 c
s
dz λz 2
=⇒ z ′ = = ± (1 + ) − 1.
dx n0
n0 p du λ
du = ± u2 − 1dx =⇒ √ = ± dx.
λ 2
u −1 n0
12