Exercices de cours du chapitre IV : systèmes à N DDL 19 Exercice IV-6: Double p

Exercices de cours du chapitre IV : systèmes à N DDL 19 Exercice IV-6: Double pendule Objectif : Comparaison des réponses linéarisée et non linéarisée (sous Matlab). Le Double pendule est constitué de deux masses 1 2 ( , ) m m soumises à leur poids propre. Elles sont reliées par deux fils inextensibles de longueur respective 1 2 ( , ) A A . Effectuez la mise en équations, sans linéarisation. Pour simplifier la programmation on prendra : 1 2 m m m = = , et 1 2 = = A A A Linéariser les équations et déterminez la réponse à des conditions initiales données. Sous Matlab comparer les solutions linéarisée et non linéarisée. θ1 θ2 y G o G g G xo 1 A 2 A 1 m 2 m Corrigé de l’exercice IV-6: Double pendule Mise en équations : • ( ) 2 2 2 1 1 1 2 1 1 1 2 2 2 2Ec m m y y θ θ θ = + + G G    A A A 2 2 2 2 1 2 1 1 2 2 2 2 1 2 1 2 2 1 2 ( ) 2 cos( ) Ec m m m m θ θ θ θ θ θ = + + + −     A A A A • ( ) 1 1 1 2 1 1 2 2 cos cos cos Ep m g m g θ θ θ = − − + A A A • 0 T δ = Équations de Lagrange, non linéarisées : θ1 θ2 y G o G g G xo 1 A 2 A 1 m 2 m 2 2 1 2 1 1 2 1 2 2 2 1 2 1 2 2 2 1 1 2 1 1 ( ) cos( ) sin( ) ( ) sin 0 m m m m m m g θ θ θ θ θ θ θ θ + + − − − + + =    A A A A A A 2 2 2 2 2 2 1 2 1 2 1 2 1 2 1 2 1 2 2 2 cos( ) sin( ) sin 0 m m m m g θ θ θ θ θ θ θ θ + − + − + =    A A A A A A Avec 1 2 m m m = = et 1 2 = = A A A 2 1 2 2 1 2 2 1 1 2 cos( ) sin( ) 2 sin / 0 g θ θ θ θ θ θ θ θ + − − − + =    A 2 2 1 2 1 1 2 1 2 cos( ) sin( ) sin / 0 g θ θ θ θ θ θ θ θ + − + − + =    A Ce que nous pouvons écrire sous la forme : 2 2 1 1 2 2 1 1 2 2 1 2 1 2 1 2 2 cos( ) sin( ) 2 sin / cos( ) 1 sin( ) sin / g g θ θ θ θ θ θ θ θ θ θ θ θ θ θ − ⎧ ⎫ ⎧ ⎫ − − ⎡ ⎤ = ⎨ ⎬ ⎨ ⎬ ⎢ ⎥ − − − − ⎣ ⎦⎩ ⎭ ⎩ ⎭   A   A Pour résoudre ces équations différentielles couplées non linéaire d’ordre deux, nous allons nous ramener à un système différentiel d’ordre un en effectuant un changement de variables. Posons { } 1 1 2 2 1 2 3 4 T X x x x x θ θ θ θ =< > = < >   En notant : 31 3 1 2 1 cos( ) cos( ) C x x θ θ = − = − , 31 3 1 sin( ) S x x = − etc Nous obtenons : { } 2 2 31 4 31 1 4 2 31 1 31 3 1 0 0 0 0 2 0 2 / 0 0 1 0 0 0 1 / x C x S gS X x C x S gS ⎡ ⎤ ⎧ ⎫ ⎢ ⎥ ⎪ ⎪ − ⎪ ⎪ ⎢ ⎥ = ⎨ ⎬ ⎢ ⎥ ⎪ ⎪ ⎢ ⎥ ⎪ ⎪ − − ⎢ ⎥ ⎣ ⎦ ⎩ ⎭ A  A Système différentiel à quatre équations que nous pouvons résoudre avec Matlab en utilisant le solveur « ODE45 » Linéarisation : 2 2 2 2 2 1 2 1 2 2 2 2 Ec m m θ θ θ θ ≅ + +     A A A Î[ ] 2 2 1 1 1 M m ⎡ ⎤ = ⎢ ⎥ ⎣ ⎦ A Exercices de cours du chapitre IV : systèmes à N DDL 20 2 2 1 2 2 2 Ep mg mg θ θ ≅ + A A Î[ ] 2 0 0 1 K mg ⎡ ⎤ = ⎢ ⎥ ⎣ ⎦ A Fréquences et modes propres : [ ] [ ] ( ) det 0 K M λ − = Î: 2 1 2 1 2 g ω = + A , 2 2 2 2 1 g ω = − A Modes propres : … Î{ } 1 2 2 T Z =< > et { } 2 2 2 T Z =< −> non normés [ ] 2 2 2 2 Z ⎡ ⎤ = ⎢ ⎥ − ⎣ ⎦ Î [ ] 1 2 1 1 4 2 1 Z − ⎡ ⎤ = ⎢ ⎥ − ⎢ ⎥ ⎣ ⎦ Oscillations libres : 2 i 0 i i q q ω + =  Î cos sin io i io i i i q q q t t ω ω ω = +  Lâcher sans vitesse initiale à t=0 Î 0 0 i q =  Dans une position donnée Î { } [ ] 1 1 1 2 2 1 2 2 1 4 2 o o o o o o o q Z θ θ θ θ θ θ − ⎧ ⎫ + ⎧ ⎫ ⎪ ⎪ = = ⎨ ⎬ ⎨ ⎬ − ⎩ ⎭ ⎪ ⎪ ⎩ ⎭ Puis [ ]{ } 1 1 1 2 2 2 1 1 2 2 2( cos cos ) 2( cos cos ) o o o o q t q t Z q q t q t θ ω ω θ ω ω ⎧ ⎧ ⎫ + ⎪ = = ⎨ ⎬ ⎨ − ⎪ ⎩ ⎭ ⎩ Solution que nous pouvons directement programmer dans Matlab La comparaison des deux solutions sera réalisée en simulant le mouvement pour différentes positions initiales. A titre d’exercice Vous pouvez aussi compléter le script Matlab pour prendre en compte des vitesses initiales non nulles dans la solution linéarisée. (Simple) Vous pouvez reprendre l’exercice pour deux longueurs de pendule différentes et modifier en conséquence le script Matlab. (Un peu plus long) uploads/s3/ pendule-double-lagrange-corrige.pdf

Tags
Administrationéquations pendule matlab double linéarisée
Documents similaires
  • 27
  • 0
  • 0
Afficher les détails des licences
Licence et utilisation
Gratuit pour un usage personnel Attribution requise
Partager