Une masse est attachée à une tige légère de[br] longueur fixe [math]\ell[/math] tournant autour d'un axe. Le système a un[br] unique degré de liberté, l'angle [math]y[/math] qu'il fait avec la[br] verticale. Les forces exercées sur le mobile sont: la tension dans[br] la tige (le long de la tige), la pesanteur ([math]g[/math] constante et vers le[br] bas) et le frottement (proportionnel et contraire à la vitesse [math]\dot y[/math]). Le[br] fait que la tige est rigide nous indique que la tension annule la[br] composante radiale de la force et que seule compte sa composante[br] angulaire, projection du poids perpendiculairement à la tige,[br] donnant finalement l'équation:[br] [math] \ddot y+\mu \dot y+\frac g\ell \sin(y)=0.[/math][br][br] Cette équation n'est pas linéaire, notez la fonction [math]\sin[/math], ce qui[br] fait qu'elle est difficile à résoudre. On peut cependant, en prenant[br] un développement limité pour de petits angles, linéariser[br] cette équation.
Nous visualisons dans l'espace des phases [math](y,\dot y)[/math] (angle, vitesse) le champ de vecteurs que doit suivre le pendule. Si on lance le pendule assez fort, il va faire plusieurs tours avant que le frottement ne contraigne les oscillations à s'amenuiser. Décrivons la courbe rouge: partant d'une vitesse quasi nulle (pendule quasi en équilibre instable), l'angle négatif -π augmente, gagnant de la vitesse vers le bas, jusqu'à un angle nul et une vitesse maximale, puis la vitesse diminue, l'angle augment encore en ralentissant, jusqu'à une amplitude de battement maximale, une vitesse qui s'annule puis devient négative, le pendule retombe, jusqu'à une vitesse négative maximale pour un angle nul, et il remonte de l'autre côté. Le frottement a fait son office et l'oscillation est de moindre amplitude, le pendule est piégé dans un bassin autour d'un point d'équilibre (un angle multiple de 2π, ici 0).[br][br]Vous pouvez modifier le point de départ, le coefficient de frottement µ et le pas de la discrétisation. Observez deux intégrations numériques, en bleu la méthode d'Euler, en rouge la méthode de Runge-Kutta.