Équations différentielles · Amortissement

Amortissement

L’amortissement retire de l’énergie à une oscillation. Réglez ζ pour comparer les retours oscillants et ceux qui ne franchissent pas l’équilibre. Le temps est mesuré en T₀, la période sans amortissement.

01

L’amortissement dans cette visualisation

L’amortissement est la perte d’énergie qui fait disparaître une oscillation. Pour une masse sur un ressort muni d’un amortisseur, mx′′+cx′+kx=0m x'' + c x' + k x = 0, ce qui s’écrit x′′+2ζω0x′+ω02x=0x'' + 2\zeta\omega_0 x' + \omega_0^2 x = 0 avec la pulsation propre ω0=k/m\omega_0 = \sqrt{k/m} et le taux d’amortissement ζ=c/(2mk)\zeta = c/(2\sqrt{mk}). En dessous de 1, le mouvement est sous-amorti (il oscille) ; à 1, l’amortissement est critique ; au-dessus de 1, le mouvement est suramorti (il revient lentement). Pour une masse suspendue verticalement, la gravité ne fait qu’abaisser la position d’équilibre de mg/k ; x est mesuré à partir de cette position. La page lâche la masse sans vitesse initiale en x = 1 et cherche si davantage d’amortissement l’arrête plus tôt.

Les racines déterminent le mouvement ; une racine double donne aussi un terme t·e^(λt)
λ2+2ζω0λ+ω02=0,λ=ω0(−ζ±ζ2−1),σ=−max⁡Re⁡λ\lambda^2 + 2\zeta\omega_0\lambda + \omega_0^2 = 0, \qquad \lambda = \omega_0\left(-\zeta \pm \sqrt{\zeta^2 - 1}\right), \qquad \sigma = -\max \operatorname{Re}\lambda

À droite se trouve le plan des racines. Le produit des deux racines vaut ω02\omega_0^2 : tant que le mouvement oscille, elles sont donc sur le cercle de rayon ω₀. La partie réelle −ζω0-\zeta\omega_0 fixe la décroissance, la partie imaginaire ωd=ω01−ζ2\omega_d = \omega_0\sqrt{1-\zeta^2} la pulsation : l’amortissement ralentit donc aussi chaque oscillation. À ζ = 1, les racines se rejoignent en −ω₀. Au-delà, l’une part vers la gauche et l’autre revient lentement vers 0, à environ −ω0/(2ζ)=−k/c-\omega_0/(2\zeta) = -k/c pour ζ grand. La composante la plus lente du mouvement décroît comme e−σte^{-\sigma t}, et ce taux de décroissance σ, distance de la racine la plus à droite à l’axe imaginaire, est maximal exactement à ζ = 1. Sur x(t), la courbe bleue en pointillés le montre : l’enveloppe ±e−ζω0t/1−ζ2\pm e^{-\zeta\omega_0 t}/\sqrt{1-\zeta^2} en régime sous-amorti (la courbe la touche entre ses pics), le terme de la racine lente en régime suramorti.

« S’arrêter » demande une règle. La page utilise T1%T_{1\%}, le premier instant après lequel |x| reste définitivement à moins de 1 % du déplacement initial. La bande sous x(t) agrandit cette zone. Les extrema se trouvent en tn=nπ/ωdt_n = n\pi/\omega_d avec ∣x(tn)∣=e−ζω0tn|x(t_n)| = e^{-\zeta\omega_0 t_n}, ce qui détermine exactement quels pics sortent de la bande ; le dernier retour est ensuite trouvé par dichotomie sur la solution explicite, donc la valeur affichée est numérique. Un ζ plus petit donne un temps d’établissement légèrement plus tardif dès qu’un pic supplémentaire sort de la bande : la courbe de T1%T_{1\%} en fonction de ζ présente donc des sauts. Son minimum n’est pas à l’amortissement critique, mais à ζ∗=ln⁡100/π2+ln⁡2100≈0.826\zeta^* = \ln 100/\sqrt{\pi^2 + \ln^2 100} \approx 0.826, où l’unique dépassement touche exactement −1 % : 0.667 T₀ contre 1.057 T₀. L’amortissement critique n’est le plus rapide que parmi les mouvements qui ne franchissent jamais x = 0, et seulement sans vitesse initiale ; une autre bande déplace ζ* (une bande de 2 % donne environ 0.780).

02

À observer

  • Deux sens différents de « plus rapide » : le taux de décroissance σ est maximal exactement à ζ = 1 ; le temps d’établissement est minimal à ζ*, un peu en dessous. La valeur littérale 0.826 n’est pas ζ* : elle dépasse de 1.0015 %, et son temps d’établissement saute à 0.896 T₀. La page garde ζ* exact.
  • Doubler un taux d’amortissement élevé double approximativement le temps du retour lent : la constante de temps de la racine lente est proche de c/k.
  • x(t) est la solution exacte dans les trois régimes ; pour un taux saisi, le régime est déterminé sur le nombre tel qu’il est écrit : 0.99999999999999999 est donc sous-amorti. Le fluide dans le cylindre n’est qu’un schéma.
03

Commandes

  • Déplacez le curseur ζ (il s’aligne sur 0.1, 1 et 10), saisissez une fraction ou un nombre décimal, ou appuyez sur ζ = 1 ou ζ*. Une racine peut être déplacée le long de sa trajectoire ; elle s’aligne en −ω₀ et ne change de branche qu’à cet endroit. Faites glisser le curseur temporel rouge pour déplacer la masse.
  • Le Mode présentation ouvre une vue en plein écran : ← → modifient ζ d’un pas du curseur, H masque le taux de décroissance et les temps d’établissement.
04

Voir aussi

Sur la carte trace-déterminant de la page du portrait de phase, l’oscillateur amorti est la droite horizontale Δ = ω₀².

Pour aller plus loin: Wikipédia: Damping (Anglais); OpenStax, University Physics Volume 1, 15.5 Damped Oscillations (Anglais).