2010-2011 Mécanique Analytique , Corrigé 5 Assistants : jaap.kroes@ep .ch & ben

2010-2011 Mécanique Analytique , Corrigé 5 Assistants : jaap.kroes@ep .ch & benjamin.audren@ep .ch Exercice 1 : Problème de la brachistochrone Dans cet exercice, nous cherchons la courbe dans le plan (x, y) minimisant le temps de parcours entre l'origine O et le point P(x0, y0) dans un champ gravi que constant. Il nous faut tout d'abord utiliser la conservation de l'énergie pour exprimer la vitesse de la particule en fonction de y. En eet, en prenant y = 0 pour niveau de référence du potentiel gravi que on a 1 2mv2 −mgy = 0 → v(y) = p 2gy La stratégie est maintenant de minimiser le temps de parcours entre O et P. Le temps nécessaire pour parcourir une longeur ds = p dx2 + dy2 est dt = ds/v(y). La fonction à minimiser dépendra donc explicitement de y. Le point clé ici est de bien comprendre l'analogie avec la minimisation de l'action comme vue en cours. Lorsque l'on dé ni l'action comme intégrale sur le temps du Lagrangien, c'est en cherchant à minimiser l'action que l'on trouve les équations d'Euler-Lagrange. Ici, on cherche à minimiser le temps de parcours, que l'on écrit comme une intégrale sur une distance (dx ou dy selon le choix) d'une quantité. Chercher à minimiser ce temps de parcours revient à trouver des equations d'Euler-Lagrange pour la quantité intégrée. Nous devons maintenant choisir quelle paramétrisation utiliser : soit (x, y(x)) soit (x(y), y). Nous allons voir que, dans le cadre ce cet exercice, les deux choix mènent à une quantité conservée qui nous permettra de simpli er la résolution du système. En eet nous allons voir que dans le premier cas la fonction hamil- tonienne est une quantité conservée, alors que dans le second une variable cyclique apparaîtra. Commençons par le choix (x(y), y). L'élément de ligne s'écrit alors comme ds = √ 1 + x′2 dy et le temps à minimiser comme T = 1 √2g Z y0 0 dy s 1 + x′2 y ≡ 1 √2g Z dy F On voit qu'il s'agit de trouver la courbe réalisant le meilleur compromis entre une courte distance et une grande vitesse. Comme discuté ce choix fait apparaître une variable cyclique qui nous permet de résoudre le système plus facilement : d dy  ∂F ∂x′  = ∂F ∂x = 0 → ∂F ∂x′ = ∂ ∂x′ s 1 + x′2 y = const = a En eectuant la dérivée par rapport à x′, l'équation d'Euler-Lagrange prend la forme x′ √ 1 + x′2 = a√y → x′ = r y b −y avec b = a−2. La forme de la courbe est donc donnée par l'intégrale suivante : x(y) = Z dx + c = Z dy dx dy + c = Z dy r y b −y + c Voyons maintenant que l'on retrouve le même résultat en utilisant la paramétrisation (x, y(x)). L'élé- ment de ligne serait alors ds = p 1 + y′2 dx et l'intégrale à minimiser T = 1 √2g Z x0 0 dx s 1 + y′2 y ≡ 1 √2g Z dx G 1 L'équation d'Euler-Lagrange serait alors de la forme d dx ∂G ∂y′  = ∂G ∂y qui n'admet pas de variable cyclique. Le point important est de noter que la fonction G ne dépend pas explicitement de x et donc que la fonction hamiltonienne est conservée : h = ∂G ∂y′ y′ −G = − 1 p y(1 + y′2) = const En séparant les variables on obtient y′ = dy dx = s b −y y → x(y) = Z dx = Z dy r y b −y + c où b = h−2. On a donc bien retrouvé la même équation qu'auparavant. Des deux paramétrisations découle une quantité conservée qui conduit à une simpli cation du problème. La manifestation de cette quantité conservée est propre au choix de la paramétrisation : pour l'une une variable cyclique, pour l'autre la fonction hamiltonienne. Nous pouvons maintenant résoudre la dernière intégrale a n de trouver la forme explicite du rail. En faisant le changement de variables y = b sin2 θ on obtient x(θ) = 2b Z dθ sin2 θ + c = b 2(2θ −sin 2θ) + c Etant donné que la courbe passe par le point O, la constante c est nulle. En posant ϕ = 2θ et R = b/2, on a la forme paramétrique suivante de la courbe recherchée : ( x(ϕ) = R(ϕ −sin ϕ) y(ϕ) = R(1 −cos ϕ) qui est l'équation paramétrique d'une cycloïde de rayon R. Noter qu'en calculant dt = ds/v en utilisant x′ = dx dy = dx dϕ  dy dϕ −1 on a dt = s 1 + x′2 2gy dy = 1 √2gR s 2 −2 cos ϕ sin2 ϕ 1 √1 −cos ϕR sin ϕdϕ = s R g dϕ et donc T = p R/g ϕ0 où ϕ0 est l'angle nal, c'est à dire l'angle pour lequel x(ϕ0) = x0 et y(ϕ0) = y0. On peut maintenant remarquer une chose a priori surprenante : il existe des situations pour lesquelles la forme idéale du rail passe au dessous de P avant de remonter. En eet, ce genre de trajectoire est caractérisé par 0 ! = dy dx = sin ϕ 1 −cos ϕ Il faut donc que ϕ0 > π pour que cela puisse se produire. Cette condition peut être traduite sur la pente entre O et P. En eet on a ϕ0 > π ⇐ ⇒ y0 x0 = 1 −cos ϕ0 ϕ0 −sin ϕ0 < 2 π 2 -2 -1.5 -1 -0.5 0 0 1 2 3 4 5 6 t-sin(t), -1+cos(t) Fig. 1  Exemple de rail admettant un minimum : ϕ0 = 3π/2 Exercice 2 : Problème de la surface minimale de révolution Tout comme dans l'exercice précédent, nous avons la liberté de paramétriser la courbe en fonction de x ou de y. Nous allons tout d'abord prendre l'option (x, y(x)) a n d'obtenir une variable cyclique puis montrer que l'autre paramétrisation conduit à une autre quantité conservée : la fonction hamiltonienne. L'aire d'une bande de la surface vaut 2πxds = 2πx p 1 + y′(x)2 dx donc la surface totale mesure A = 2π Z x2 x1 dx x q 1 + y′2 Pour que cette dernière soit extrémale, il faut résoudre l'équation d'Euler-Lagrange d dx ∂F ∂y′  = ∂F ∂y → ∂F ∂y′ = ∂ ∂y′ x p 1 + y′2 = const = a On obtient donc en eectuant la dérivée par rapport à y′ : xy′(x) = a p 1 + y′(x)2 → y′(x) = r a2 x2 −a2 En eectuant le changement de variables x = a cosh θ on trouve y(θ) = aθ + b → x(y) = a cosh y −b a  On détermine a et b à partir des conditions x(0) = 1 et x(2) = cosh(2), ce qui donne a = 1 et b = 0. Voyons maintenant que de l'autre paramétrisation découle la conservation de la fonction hamiltonienne. L'intégrale à minimiser s'écrit cette fois-ci comme A = 2π Z y2 y1 dy x p 1 + x′2 ≡2π Z y2 y1 dy A où la fonction A ne dépend pas explicitement de y. La fonction hamiltonienne est donc conservée : h = ∂A ∂x′ x′ −A = − x √ 1 + x′2 = const 3 En séparant les variables on trouve x′ = dx dy = r x2 −a2 a2 → y(x) = Z dy = Z dx r a2 x2 −a2 où a2 = h2, qui est exactement la même expression que l'on a trouvée avec la première paramétrisation. Exercice 3 : Rotateurs couplés Le Lagrangien du système est donné par L = 1 2m1a2 ˙ θ2 1 + 1 2m2b2 ˙ θ2 2 −1 2kl2 où l est la distance entre les deux points matériels. Les positions des deux masses sont : ⃗ x1 =   a cos θ1 a sin θ1 0   ⃗ x2 =   b cos θ2 b sin θ2 z   La distance l est alors : l2 = ||⃗ x1 −⃗ x2||2 = a2 + b2 + z2 −2ab cos(θ1 −θ2). Le Lagrangien ne dépendant pas explicitement du temps, la fonction hamiltonienne est conservée : h = 2 X i=1 ∂L ∂˙ θi ˙ θi −L = 1 2m1a2 ˙ θ2 1 + 1 2m2b2 ˙ θ2 2 + 1 2kl2 On remarque qu'il n'y a pas de variable cyclique, toutefois physiquement, on s'attend à la conservation du moment cinétique total. On peut le trouver en appliquant le théorème de Noether. Le Lagrangien est invariant par rotation, c'est-à-dire sous la transformation θi →θi(s) = θi + s La quantité suivante est donc conservée : 2 X i=1 ∂L ∂˙ θi ∂θi(s) ∂s = m1a2 ˙ θ1 + m2b2 ˙ θ2 Le moment cinétique est donc bien conservé ! On peut aussi le voir en eectuant le changement de variable suivant uploads/s3/ corrige-05.pdf

Tags
Administrationfonction conservée minimiser comme quantité
Documents similaires
  • 61
  • 0
  • 0
Afficher les détails des licences
Licence et utilisation
Gratuit pour un usage personnel Attribution requise
Partager