Voir la source de Oscillateur harmonique

{{Voir homonymes|Oscillateur}}
Un '''oscillateur harmonique''' est un [[oscillateur (mécanique)|oscillateur]] idéal dont l'évolution au cours du temps est décrite par une [[fonction trigonométrique|fonction sinusoïdale]], dont la [[fréquence]] ne dépend que des caractéristiques du système et dont l'amplitude est constante.

Ce [[modèle mathématique]] décrit l'évolution de n'importe quel [[système physique]] au voisinage d'une [[Équilibre statique (mécanique)|position d'équilibre stable]], ce qui en fait un outil transversal utilisé dans de nombreux domaines : [[mécanique (industrie)|mécanique]], [[électricité]] et [[Électronique (technique)|électronique]], [[optique]]. Il néglige les forces dissipatives ([[frottement]] par exemple). Dans la pratique, pour que l'[[amplitude]] de l'oscillation reste constante, il faut entretenir les oscillations en fournissant de l'[[énergie]].

== Définition générale et propriétés ==
=== Approche expérimentale : exemples d'oscillateurs réels ===
[[Fichier:Simple harmonic oscillator.gif|vignette|gauche|upright=0.3|Oscillation d'une masse suspendue à un ressort.]]
De nombreux systèmes physiques réels, pour lesquels on rend négligeable l'influence des frottements, se comportent comme des systèmes oscillants, qui peuvent être modélisés comme des oscillateurs harmoniques.

{{Article détaillé|système masse-ressort}}
* '''[[Système masse-ressort]]''' : Il s'agit d'un système très simple à réaliser. Si la masse est écartée légèrement de sa position d'équilibre et relâchée sans vitesse initiale, l'expérimentateur constate que celle-ci se met à osciller autour de cette position d'équilibre. La mesure, à l'aide d'un capteur adapté, de l'évolution au cours du temps de l'élongation ''x''(''t'') du [[ressort]] (par rapport à la position où la masse est au repos) montre que les (petites) oscillations sont purement sinusoïdales, de pulsation <math>\omega_0</math>, du moins si l'influence des frottements est négligée<ref>Il suffit d'effectuer les mesures pendant une durée assez faible, de l'ordre de quelques périodes, pour ne constater qu'un amortissement négligeable.</ref>. La valeur de <math>\omega_0</math> ne dépend ''pas'' de l'amplitude initiale de déplacement de la masse, mais uniquement de la valeur de la masse et des propriétés du ressort, plus précisément de sa [[Raideur (mécanique)|raideur]] ''k''.
: On peut montrer que l'équation du mouvement de la masse s'écrit <math>\frac{\mathrm{d}^2x}{\mathrm{d}t^2} +\omega_0^2 \, x(t)=0</math>, avec <math>\omega_0=\sqrt{\frac{k}{m}}</math>.
:Par suite l'équation horaire est donnée par <math>x(t)=x_m\cos{(\omega_0 t +\phi_0)}</math>, où <math>x_m</math> est '''l'amplitude''' des oscillations et <math>\phi_0</math> la phase à l'origine, qui dépendent des conditions initiales. Le mouvement est qualifié '''d'harmonique''', de pulsation propre <math>\omega_0</math>, la caractéristique du système.
[[Fichier:HarmOsc1.png|cadre|centré|Représentation schématique de l'évolution dans le temps de l'élongation de la masse ''m''.]]

{{démonstration|contenu=
Soit (''R'') le référentiel, considéré comme galiléen, lié au support sur lequel est fixé le ressort, ''Ox'' l'axe parallèle à la direction de l'axe du ressort, orienté positivement vers le bas, d'origine ''O'' la position de la masse ''m'' considérée comme ponctuelle à l'équilibre, <math>\ell_0</math> l'élongation du ressort à vide, <math>\ell</math> celle du ressort à l'équilibre. Lorsque la masse ''m'' est au repos il vient (en projection sur ''Ox'') l'équation <math>mg-k(\ell-\ell_0)=0</math>.
L'élongation totale du ressort pour une position <math>x=x(t)</math> quelconque de la masse ''m'' est donnée par <math>\Delta\ell=(\ell-\ell_0)+x(t)</math>, par suite l'application à celle-ci du principe fondamental de la dynamique dans le référentiel (''R'') donne en projection sur ''Ox'' :
:<math>m\frac{\mathrm{d}^2x}{\mathrm{d}t^2}=mg-k\left[(\ell-\ell_0)+x(t)\right]</math>,
soit compte tenu de l'équation statique précédente :
:<math>m\frac{\mathrm{d}^2x}{\mathrm{d}t^2}=-kx</math>,
soit encore <math>\frac{\mathrm{d}^2x}{\mathrm{d}t^2}+\frac{k}{m}x=0</math>.}}

[[Fichier:Lc circuit.svg|100px|vignette|gauche|Les oscillations libres d'un circuit LC constituent un exemple d'oscillations harmoniques.]]
[[Fichier:Curvas Oscilador LC.gif|vignette|Courbe de la tension aux bornes du condensateur (en haut) et aux bornes de la bobine (en bas). Celles-ci sont en opposition de phase, montrant le transfert d'énergie entre le condensateur et la bobine.]]
* '''Oscillations libres d'un [[circuit LC]]''': un [[Condensateur (électricité)|condensateur]] de [[capacité électrique]] ''C'', initialement chargé avec une charge <math>q_m</math> commence à se décharger à un instant choisi comme origine dans une [[Bobine (électricité)|bobine]] d'[[inductance]] ''L'', dont la résistance est considérée comme négligeable (cf. figure ci-contre). Il est facile de visualiser à l’[[oscilloscope]] l'évolution de la tension <math>u_C(t)</math> aux bornes du condensateur, et donc celle de sa charge du fait de la relation <math>u_C(t)=\frac{q(t)}{C}</math><ref>En toute rigueur cette relation n'est valable en régime variable que dans l'approximation des régimes quasi stationnaires, qui est bien vérifiée ici vu la taille du circuit. Cf.Perez, ''Electromagnétisme'', {{2e|édition}}, Paris, Masson, 1994.</ref>. Les résultats de mesures montrent que <math>q(t)</math> a une évolution sinusoïdale pure au cours du temps, de fréquence <math>\omega_0</math>, dont la valeur ne dépend pas là encore de la charge initiale du condensateur, mais uniquement des valeurs de ''L'' et de ''C''.

: L'utilisation de la [[loi des mailles]] pour le circuit considéré permet de montrer que l'évolution de <math>q(t)</math> obéit à l'équation différentielle <math>\frac{\mathrm{d}^2q}{\mathrm{d}t^2} +\omega_0^2 \, q(t)=0</math>, avec <math>\omega_0=\frac{1}{\sqrt{LC}}</math> pulsation propre des oscillations libres du circuit LC. Il vient l'expression <math>q(t)=q_m\cos{(\omega_0 t +\phi_0)}</math> pour l'évolution de la charge du condensateur, les valeurs de l'amplitude <math>q_m</math> et de la phase à l'origine <math>\phi_0</math> dépendant des conditions initiales.

Les analogies entre l'oscillateur mécanique qu'est le système masse-ressort et l'oscillateur électrique qu'est le circuit LC permettent d'identifier entre elles certaines grandeurs mécanique et électrique (voir infra [[#Analogie électro-mécanique|analogie électromécanique]]). Elles montrent aussi que les grandeurs permettant de caractériser l'état du système obéissent à une [[équation différentielle]] de même forme, dont les solutions sont harmoniques sinusoïdales. Il est ainsi possible de définir de façon générale un oscillateur harmonique comme un système dont l'évolution libre (non amortie du moins) est caractérisée par une équation de cette forme.

=== Oscillateur harmonique libre unidimensionnel ===
[[Fichier:Simple harmonic motion animation.gif|vignette|Illustration animée du mouvement harmonique unidimensionnel.]]
Les considérations expérimentales précédentes permettent de définir de façon générale un oscillateur harmonique, pour un système physique quelconque, "unidimensionnel" (ou à un [[degré de liberté (physique)|degré de liberté]]), c'est-à-dire dont l'état peut être décrit à un instant ''t'' donné par une seule [[grandeur physique]], notée ''s''(''t'')<ref>Cette grandeur est générale, et ne représentera par nécessairement une coordonnée géométrique comme en mécanique, il pourra s'agir d'une grandeur électrique par exemple, cf. exemple du circuit ''LC'' donné plus haut.</ref>. La nature précise de cette grandeur physique dépend du système considéré (cf. exemples précédents).

Un tel système se comporte comme un oscillateur harmonique (idéal) si lorsque le système est "libre", c'est-à-dire pouvant être considéré comme isolé, la grandeur physique ''s''(''t'') obéit à l'équation différentielle d'évolution suivante:
:<math>\frac{\mathrm{d}^2s}{\mathrm{d}t^2}+\omega_0^2\, s(t)=0</math>, (1),
avec <math>\omega_0</math> constante réelle positive, ayant les dimensions d'une ''fréquence'', appelée ''pulsation propre''<ref>La pulsation <math>\omega </math>, qui s'exprime en rad.s{{exp|-1}}, est une fréquence ''angulaire'', liée à la fréquence ''f'' (en Hz) par la relation <math>\omega=2\pi f</math>. Si le terme de ''pulsation propre'' est le plus rigoureux pour désigner <math>\omega_0 </math>, celui de ''fréquence propre'' est couramment employé à sa place en pratique.</ref> du système, qui ne dépend que des caractéristiques propres de celui-ci.

La solution générale de (1) <math>s(t)=s_m \cos{(\omega_0 t+\phi_0)}</math>, avec <math>s_m = \mathrm{cte}</math> amplitude des oscillations et <math>\phi_0 = \mathrm{cte}</math> phase à l'origine, ces deux constantes étant déterminées par les conditions initiales sur ''s''(''t'') et sa dérivée première.

=== Portrait de phase ===
[[Fichier:Pendulum 45deg.gif|vignette|271x271px|Portrait de phase d'un pendule]]
L'état de l'oscillateur à un instant donné est totalement déterminé par la connaissance de la position et de la vitesse. Puisque l'oscillateur n'a qu'un degré de liberté, un seul paramètre représente la position, et un seul la vitesse. Le [[portrait de phase]] d'un oscillateur harmonique est donc l'ensemble des trajectoires décrites par le système à partir de toutes les conditions initiales, représenté dans le plan <math>\left( s, \frac{\mathrm{d}s}{\mathrm{d}t}\right)</math> .

Le portrait de phase d'un oscillateur harmonique libre est une [[Ellipse (mathématiques)|ellipse]].

=== Interprétation en termes énergétiques ===
==== Mise en évidence de la conservation de l'énergie ====
Par définition un système idéal est un système dans lequel l'énergie totale est conservée, c'est-à-dire qu'elle demeure constante à travers le temps. Pour mettre l'aspect idéal de ce système en évidence et la signification physique de l'équation (1), il est très intéressant de multiplier celle-ci membre à membre par <math>\frac{\mathrm{d}s}{\mathrm{d}t}</math>. Il vient alors :
:<math>\frac{\mathrm{d}^2s}{\mathrm{d}t^2}\frac{\mathrm{d}s}{\mathrm{d}t}+\omega_0^2\,\frac{\mathrm{d}s}{\mathrm{d}t}s(t)=0</math>, soit encore <math>\frac{\mathrm{d}}{\mathrm{d}t}\left[\frac{1}{2}\left(\frac{\mathrm{d}s}{\mathrm{d}t}\right)^2+\frac{1}{2}\, \omega_0^2 \, s^2(t)\right]=0</math>.
Par suite pour un oscillateur harmonique la quantité <math>\frac{1}{2}\left(\frac{\mathrm{d}s}{\mathrm{d}t}\right)^2+\frac{1}{2}\, \omega_0^2 \, s^2(t) = e</math> est '''conservée'''. Aux dimensions près, celle-ci correspond à l''''énergie du système''' (de fait il sera toujours possible en pratique de multiplier ''e'' par des combinaisons de paramètres du système pour obtenir les "bonnes" dimensions).

Cette "énergie" est constituée de deux termes :
* un terme "d'énergie potentielle" en <math>\frac{1}{2}\omega_0^2 s^2(t)</math>;
* un terme "d'[[énergie cinétique]]" en <math>\frac{1}{2}\left(\frac{\mathrm{d}s}{\mathrm{d}t}\right)^2</math>, (défini ici à une constante multiplicative près).
La [[forme quadratique]] de l'énergie "potentielle" est en fait caractéristique d'un oscillateur harmonique : en effet, ce n'est qu'avec cette forme que la prise en compte du caractère conservatif du système <math>\frac{\mathrm{d}e}{\mathrm{d}t}=0</math> redonne l'équation d'évolution précédente.

Bien entendu, la nature exacte et les noms de ces termes dépendent du système étudié :
* pour le système masse-ressort, <math>s(t)=x(t)</math>, <math>\frac{\mathrm{d}x}{\mathrm{d}t}=v(t)</math> est la vitesse instantanée de la masse à l'instant ''t'', et <math>\omega_0^2=\frac{k}{m}</math>. Par suite <math>e=\frac{1}{2}v^2(t)+\frac{k}{2m} x^2(t)</math> et il suffit de multiplier membre à membre par la masse ''m'' pour obtenir l'énergie mécanique du système <math>E=\frac{1}{2}mv^2(t)+\frac{1}{2} kx^2(t)</math>: dans ce cas il s'agit bien de la somme de l'énergie cinétique et de l'énergie potentielle élastique de la masse ''m''. Dans la limite où les frottements sont négligeables, cette énergie mécanique est conservée.
* pour le circuit ''LC'', la grandeur ''s''(''t'') est la charge ''q''(''t'') de l'oscillateur à l'instant ''t'', <math>\frac{\mathrm{d}q}{\mathrm{d}t}=i(t)</math> est l'intensité dans le circuit, et <math>\omega_0=\frac{1}{\sqrt{LC}}</math>. Par suite <math>e=\frac{1}{2}i^2(t)+\frac{1}{2LC} q^2(t)</math> et il suffit de multiplier membre à membre cette quantité par l'[[inductance]] ''L'' pour obtenir l'énergie du système <math>E=\frac{1}{2}Li^2(t)+\frac{1}{2} \frac{q^2(t)}{C}</math>. Le terme en <math>\frac{1}{2}Li^2(t)</math>, qui correspond à l'énergie contenue dans la [[Bobine (électricité)|bobine]], est l'analogue "électrique" de l'énergie cinétique du système masse-ressort, et celui en <math>\frac{1}{2} \frac{q^2(t)}{C}</math>, qui correspond à l'énergie stockée dans le [[Condensateur (électricité)|condensateur]], est l'analogue de l'énergie potentielle mécanique.
Comme <math>s(t)=s_m \cos{(\omega_0 t +\phi_0)}</math>, l'énergie du système peut alors s'écrire sous la forme <math>e=\frac {1} {2} \omega_0^2 s_m^2</math> . L'énergie du système pour un oscillateur harmonique est donc proportionnelle au carré de l'amplitude des oscillations.

==== Interprétations des oscillations en termes énergétiques ====
Le caractère conservatif de l'énergie d'un oscillateur harmonique idéal est bien sûr lié à l'absence de termes de dissipation. Il est en effet facile de voir que l'équation (1) est invariante par renversement du temps <math>t\rightarrow -t</math> : une telle symétrie implique que l'énergie du système physique soit conservée (cf. [[Théorème de Noether (physique)|théorème de Noether]]). Cette [[conservation de l'énergie]] totale du système implique qu'au cours de son évolution il y a transfert entre les formes "potentielle" et "cinétique" de l'énergie du système :

:en fait, le caractère « oscillant » du système peut s'interpréter comme la conséquence du « va-et-vient » entre les deux formes d'énergie, dont la somme demeure (idéalement) constante au cours du temps.

Physiquement, un oscillateur harmonique présente la propriété fondamentale d'''isochronisme des oscillations'', c'est-à-dire que celles-ci possèdent une période donnée <math display="inline">T_0</math> (période propre), constante, au bout de laquelle elles se répètent à l'identique. Sur le plan fondamental, cette propriété résulte de l''''invariance par translation dans le temps''', qui elle-même implique le caractère conservatif de l'énergie du système<ref>Il est important de souligner que la propriété d'isochronisme des oscillations n'est ''pas'' limitée au cas de l'oscillateur harmonique. Ainsi un pendule simple non amorti, dans le domaine non linéaire, exhibera des oscillations présentant cette propriété, bien que non sinusoïdales, et également en raison de l'invariance par translation dans le temps. En effet lorsque le pendule, partant d'une extrémité, arrive à une autre au bout d'une durée ''T''/2, il n'a aucune raison de mettre une durée différente pour revenir à son extrémité d'origine, puisque les lois physiques sont invariantes par translation dans le temps. Il nécessitera donc la même durée ''T''/2 pour "revenir", d'où la période de ''T''. Bien entendu en présence de frottements, la symétrie par translation dans le temps est brisée, et il n'y aura plus d'isochronisme.</ref>{{,}}<ref>Sur la relation entre loi de conservation et symétries continues des systèmes physiques, cf. l'article [[Théorème de Noether (physique)|théorème de Noether]].</ref>. Sur le plan mathématique, l'invariance par translation dans le temps se traduit par l'absence de termes d'ordre impair en ''t'', par exemple en <math display="inline">\frac{\mathrm{d}s}{\mathrm{d}t}</math>, dans l'équation différentielle d'évolution. En présence de frottements, la symétrie par translation dans le temps est brisée, et il n'y aura plus d'isochronisme des oscillations (cf. oscillateur harmonique amorti).

== L'oscillateur harmonique comme modèle - limites ==
=== Importance physique de l'oscillateur harmonique ===
Bien qu'il s'agisse d'une idéalisation, le modèle de l'oscillateur harmonique est d'une importance fondamentale en Physique. En effet, tout système à un seul [[degré de liberté (physique)|degré de liberté]] noté ''q''(''t''), pour lequel il est possible de définir une [[énergie potentielle]] (notée <math>V(q)</math> de façon générale) minimale en <math>q=q_m</math>, se comporte au voisinage de ce minimum comme un oscillateur harmonique.

En effet, au voisinage de <math>q=q_m</math> l'énergie potentielle s'écrit <math>V(q)= V(q_m)+V'(q_m)(q-q_m)+\frac{1}{2}\left(\frac{\mathrm{d}^2V}{\mathrm{d}q^2}\right)_{q=q_m}\left(q-q_m\right)^2+O(q^3)</math>, avec <math>V'(q_m)=0</math> puisque <math>V(q)</math> est minimale en <math>q_m</math>. De plus, <math>k\equiv \left(\frac{\mathrm{d}^2V}{\mathrm{d}q^2}\right)_{q=q_m} >0 </math> car l'équilibre doit être stable pour qu'il y ait des oscillations. En prenant l'origine de l'énergie potentielle en <math>q=q_m</math>, <math>V(q_m)=0</math>, et en posant <math>s(t)=q(t)-q_m</math>, l'énergie potentielle se met au voisinage de ''tout minimum'' sous la forme de celle d'un oscillateur harmonique idéal :
:<math>V(s)=\frac{1}{2}\, k\, s^2(t)</math>,
''k'' étant une "constante de raideur" (par analogie avec le système masse-ressort) de l'oscillateur harmonique<ref>Dans le cas du système masse-ressort, il s'agit bien de la raideur de celui-ci, mais dans le cas du circuit ''LC'' on a en fait <math>k\equiv \frac{1}{C}</math>.</ref>.
[[Image:Puit potentiel parabole.svg|vignette|250px|droite|Potentiel d'un pendule simple et son approximation parabolique]]
:'''Exemple : [[pendule simple]]'''
Dans le cas d'un pendule simple, le système est à un seul degré de liberté noté ''θ'' (angle entre la verticale et le fil), et l'énergie potentielle de la masse ''m'' (assimilée à un [[point matériel]]) est donnée par <math>V(\theta)=mg\ell\left(1-\cos{\theta}\right)</math>, ''g'' étant la pesanteur et <math>\ell</math> la longueur du pendule. Cette énergie potentielle est clairement minimale en <math>\theta_m=0</math> et au voisinage de cette position d'équilibre l'énergie potentielle du système s'écrit :
:<math>V(\theta)\approx \frac{1}{2}\, m\, g\, \ell \, \theta^2</math> puisque <math>\frac{\mathrm{d}^2V}{\mathrm{d}\theta^2}=m \, g \, \ell\cos{\theta}</math>, dans ce cas <math>k=m \, g \, \ell</math>. L'équation différentielle du mouvement est alors dans ce domaine celle d'un oscillateur harmonique : <math>\ddot{\theta}+\omega_0^2\, \theta=0</math>, avec <math>\omega_0=\sqrt{\frac{g}{\ell}}</math> pulsation propre du pendule simple. Elle correspond à une linéarisation, dans la limite des petites oscillations, de l'équation "complète" <math>\ddot{\theta}+\omega_0^2\sin{\theta}=0</math> donnant l'évolution du pendule<ref>L'approximation harmonique est en fait valable dans un domaine assez vaste, de l'ordre de 30° environ avec 5% d'erreur.</ref>.

Dans le cas du pendule simple, l'approximation harmonique, dite aussi des petites oscillations, revient à remplacer le potentiel réel par le potentiel parabolique, caractéristique de l'oscillateur harmonique (cf. figure ci-contre).

=== Limitations du modèle de l'oscillateur harmonique ===
Le modèle de l'oscillateur harmonique possède des limites de deux ordres.

En premier lieu, il correspond à la situation d'un système physique pris au voisinage d'un point d'équilibre stable comme indiqué précédemment. Si l'écart à l'équilibre est important, il apparaît en général des termes supplémentaires non quadratiques (dits ''anharmoniques'') dans le développement de l'énergie potentielle.

Par ailleurs, il s'agit d'un modèle conservatif, qui néglige tous les termes de dissipation d'énergie, ou alors pour lequel ceux-ci sont exactement compensés par un apport extérieur d'énergie. Dans le premier cas, l'oscillateur harmonique est dit ''amorti'', dans le second il est dit ''entretenu''.

On peut dresser un parallèle entre l'oscillateur harmonique et la [[corde de Melde]]. Cette dernière expérience permet de mettre en avant une infinité de pulsations possibles pour la corde tandis que le ressort lui ne peut osciller que selon une unique pulsation. Ceci est expliqué par l'inertie de la masse accrochée au ressort qui ne peut pas se subdiviser contrairement à la longueur de la corde à vide. Notons que le ressort seul, sans masse accrochée, pourrait à priori se comporter comme la [[Expérience de Melde|Corde de Melde]]<ref>{{Lien web|nom1=Tulloue|prénom1=Geneviève|titre=Corde de Melde (Animation Flash)|url=http://www.sciences.univ-nantes.fr/sites/genevieve_tulloue/Ondes/ondes_stationnaires/melde.php|site=www.sciences.univ-nantes.fr|consulté le=2016-09-27}}</ref>.

==== Limites du modèle pour les systèmes réels ====
Pour tout système physique réel, le comportement n'est celui d'un oscillateur harmonique que dans un domaine restreint, même en négligeant tout terme de dissipation.

Par exemple, pour le système masse-ressort le système ne se comporte comme un oscillateur harmonique que dans la limite des petites oscillations, c'est-à-dire si l'on reste dans le '''''domaine linéaire''''' du ressort, où la force de rappel est proportionnelle à son élongation<ref>Pour un ressort, la force de rappel est donné par <math>F=kx</math> ([[loi de Hooke]], où ''k'' est la constante de raideur du ressort. Il s'agit d'une [[loi de déformation]] approchée, valable pour une élongation (différence entre longueur du ressort à ''t'' et sa longueur au repos) "petite". En pratique, il est possible de tracer expérimentalement la courbe entre ''F'' et ''x'': cette courbe est un segment de droite pour ''x'' petit, la raideur ''k'' est la pente de celle-ci.</ref>. Pour des élongations plus élevées, il y a apparitions de termes non linéaires (dits ''anharmoniques'') dans la relation force de rappel / élongation. Ceux-ci se caractérisent en pratique par l'apparition de termes en ''x''{{3}} ou plus dans l'expression de l'énergie potentielle. Pour des élongations encore plus élevées, il y aura déformation permanente du ressort, qui quitte alors le domaine [[déformation élastique|''élastique'']] pour entrer dans celui de la [[déformation des matériaux|déformation]] [[déformation plastique|''plastique'']], voire se rompre. Dans ce cas-là la modélisation du système par une énergie potentielle ne sera pas possible.

Par ailleurs, il est fait l'hypothèse que le ressort revient toujours dans son état initial après déformation : en pratique, et même pour de petites élongations pour lesquelles le ressort reste dans le domaine élastique, ceci n'est pas strictement vérifié, du fait d'une certaine usure des matériaux du ressort. En fait, l'état réel du ressort dépend de son "histoire" (phénomène d'[[hystérésis]]) : cet effet est particulièrement visible sur un grand nombre d'oscillations et peut être modélisé par un terme d'amortissement dit ''structurel''<ref>Cf. Meirovitch, ''Elements of vibration analysis'', 2nd ed., {{ISBN|0-070-41342-8}}, chapitre 1.</ref>.

==== Prise en compte de la dissipation d'énergie : oscillateur harmonique amorti ====
De façon générale, l'énergie d'un oscillateur harmonique sera donnée par une expression de la forme :
:<math>E=\frac{1}{2}\kappa\left(\frac{\mathrm{d}s}{\mathrm{d}t}\right)^2+\frac{1}{2}\, \kappa \, \omega_0^2 \, s^2(t)</math>,
expression dans laquelle <math>\kappa</math> représente un paramètre multiplicatif positif, dépendant du système considéré, et tel que ''E'' ait effectivement la dimension d'une énergie : ainsi pour un système masse-ressort <math>\kappa=m</math>, pour un circuit ''LC'', <math>\kappa=L</math>, pour un pendule simple <math>\kappa=m \, \ell^2</math>, etc.

Si la dissipation d'énergie est prise en compte, la relation <math>\frac{\mathrm{d}E}{\mathrm{d}t}=0</math> traduisant le caractère conservatif du système n'est plus valide, et il faudra tenir compte de la puissance perdue par dissipation notée <math>\mathcal{P}_{\rm diss}(t)</math>, qui est toujours '''négative'''<ref>Sinon cela correspondrait à un terme d'entretien ou d'amplification des oscillations, par apport extérieur d'énergie.</ref>, et telle que :
:<math>\frac{\mathrm{d}E}{\mathrm{d}t}=\mathcal{P}_{\rm diss}(t)</math>.

En règle générale, cette puissance des forces de dissipation n'a pas d'expression simple, et peut dépendre d'un très grand nombre de paramètres, parfois extérieurs au système considéré. Par exemple, pour un oscillateur mécanique, ce terme correspondra aussi bien à des forces dissipatives ''intérieures'' au système, traduisant les "frottements" entre ses différentes sous-parties, qu'à des [[forces intérieures et forces extérieures|forces extérieures]], avec un transfert thermique entre le système et l'extérieur. En toute rigueur il faudrait alors utiliser une approche thermodynamique plutôt que purement mécanique<ref>Cf. {{Landau|tome 1}}.</ref>.

Toutefois dans le cas où la '''dissipation est faible''', et que les propriétés physiques du système sont considérées comme invariantes au cours du temps, il est possible de développer <math>\mathcal{P}_{\rm diss}(t)</math> en puissance de <math>\frac{\mathrm{d}s}{\mathrm{d}t}</math>.

A l'ordre le plus bas, l'expression de <math>\mathcal{P}_{\rm diss}(t)</math> est alors de la forme<ref>Physiquement, l'ordre le plus bas en <math>\frac{ds}{dt}</math> doit bien être quadratique, en effet le signe de <math>\frac{ds}{dt}</math> variant au cours du temps <math>\mathcal{P}_{\rm diss}(t)</math> pourrait devenir positive ce qui n'est pas possible en dehors d'un apport extérieur d'énergie.</ref> :
:<math>\mathcal{P}_{\rm diss}=-\kappa \,\alpha\left(\frac{\mathrm{d}s}{\mathrm{d}t}\right)^2</math>,
<math>\alpha</math> étant un paramètre positif, supposé constant, ayant nécessairement les dimensions de l'inverse d'une durée<ref>En effet, puisque ''E'' a les dimensions d'une énergie, il est de même de <math>\kappa\left(\frac{ds}{dt}\right)^2</math>, or une puissance ayant les dimensions d'une énergie sur une durée, <math>\alpha</math> devra bien avoir celle de l'inverse d'une durée pour que l'expression de la puissance dissipée soit homogène.</ref> (cf. ci-après pour une discussion sur son interprétation physique).

[[Fichier:HarmOsc2b.png|vignette|200 px|droite|Système masse-ressort amorti. La force d'amortissement visqueux est proportionnelle à la vitesse.]]
Sur le plan physique, une telle expression peut se justifier par le fait que l'introduction d'un terme de dissipation doit donner une équation différentielle qui ne soit plus invariante par renversement du temps, dont un terme en <math display="inline">\frac{\mathrm{d}s}{\mathrm{d}t}</math> à l'ordre le plus bas, ce que seul peut réaliser un terme de dissipation dépendant de <math display="inline">\left(\frac{\mathrm{d}s}{\mathrm{d}t}\right)^2</math> au moins.

* Dans le cas du système masse-ressort, cette expression correspond à celle de la puissance d'une force d'amortissement visqueux de la forme <math>\vec{f}=-\lambda\, \vec{v}</math>, avec <math>\lambda>0</math>.
:Un montage expérimental permettant d'obtenir cette forme peut être celui de la figure ci-contre, où la masse est relié à un dispositif dissipateur plongeant dans un liquide, dont la [[viscosité]] peut être plus ou moins forte (cf. figure ci-contre à droite).
[[Fichier:HarmOscLCR.png|vignette|gauche|Circuit LC avec introduction d'une résistance ''R'' pour modéliser les pertes.]]

* Dans le cas du circuit ''LC'', elle correspond à la prise en compte de la [[résistance (électricité)|résistance]] ''R'' de la bobine et du reste du circuit pour modéliser les sources de pertes (cf. figure ci-contre à gauche).
:La puissance dissipée par [[effet Joule]] dans celle-ci est donnée par <math>P_J=R \, i^2(t)=r\left(\frac{\mathrm{d}q}{\mathrm{d}t}\right)^2</math>, ce qui correspond bien à l'expression proposée précédente.

Par suite, en tenant compte de cette expression approchée pour l'oscillateur faiblement amorti, il vient compte tenu de <math>\frac{\mathrm{d}E}{\mathrm{d}t}=\mathcal{P}_{\rm diss}(t)</math> l'équation différentielle donnant l'évolution de la grandeur ''s''(''t'') :
:<math>\frac{\mathrm{d}^2s}{\mathrm{d}t^2}\frac{\mathrm{d}s}{\mathrm{d}t}+\omega_0^2\,\frac{\mathrm{d}s}{\mathrm{d}t} \, s(t)=-\alpha\left(\frac{\mathrm{d}s}{\mathrm{d}t}\right)^2</math>.

Contrairement au cas de l'oscillateur harmonique idéal, cette équation n'est pas symétrique par renversement du temps du fait du terme en <math>\frac{\mathrm{d}s}{\mathrm{d}t}</math>, ce qui traduit sur le plan mathématique la non-conservation de l'énergie.

L'équation différentielle donnant l'évolution de l'oscillateur harmonique (faiblement) amorti peut être réécrite sous une forme plus physique en posant <math>\alpha=\frac{2}{\tau}</math> ou <math>\tau</math> est une durée caractéristique des processus de dissipation appelé '''[[temps de relaxation]]'''. Il vient alors :
:<math>\frac{\mathrm{d}^2s}{\mathrm{d}t^2}+\frac{2}{\tau}\frac{\mathrm{d}s}{\mathrm{d}t}+\omega_0^2\,s(t)=0</math>.

La résolution mathématique de cette équation différentielle permet de distinguer trois régimes selon les valeurs relatives de <math>\tau</math> et de la pulsation propre <math>\omega_0</math> :
[[Fichier:HarmOsc4.png|vignette|gauche|Allure de ''s''(''t'') en régime sous-critique.]]
* '''régime dit pseudo-périodique''' : si <math>\omega_0>\frac{1}{\tau}</math>, il vient <math>s(t)=A e^{-\frac{t}{\tau}}\cos{(\omega t +\phi_0)}</math>, ce qui correspond à des '''oscillations sinusoïdales amorties''' de '''pseudo-fréquence'''<ref>En effet ''s''(''t'') n'est plus une fonction périodique du fait de la décroissance exponentielle.</ref> <math>\omega=\omega_0\sqrt{1-\frac{1}{\tau^2\omega_0^2}} <\omega_0</math>, ''A'' et <math>\phi_0</math> étant des constantes dépendant des conditions initiales. Il n'y a plus d'oscillations au sens strict mais s(t) prend la forme d'une sinusoïde dont l'enveloppe décroît exponentiellement au cours du temps (cf. figure ci-contre) : on parle alors d'oscillations libres ''amorties''.

:L'étude de la courbe représentative de <math>s(t)</math> permet de déterminer la pseudo-période <math>T=\frac{2\pi}{\omega}</math> des oscillations, ainsi que le '''[[décrément logarithmique]] <math>\delta\equiv\ln{\left(\frac{s(t)}{s(t+T)}\right)}</math>, qui est donné dans le cas du régime sous-critique par <math>\delta=\frac{T}{\tau}</math>. Ceci permet la mesure expérimentale du '''temps de relaxation''' <math>\tau</math>, caractéristique des phénomènes dissipatifs.

* '''régime critique''' : si <math>\omega_0=\frac{1}{\tau}</math>, il vient <math>s(t)=e^{-\frac{t}{\tau}}\left(At+B\right)</math>, ''A'' et ''B'' étant deux constantes dépendant des conditions initiales. Il n'y a pas d'oscillations amorties, le système revient exponentiellement vers son état de repos. Il s'agit d'un cas limite.
* '''régime apériodique''' : si <math>\omega_0<\frac{1}{\tau}</math>, il vient <math>s(t)=A e^{-\beta_+\frac{t}{\tau}}+B e^{-\beta_-\frac{t}{\tau}}</math>, avec <math>\beta_{\pm}=1\pm\sqrt{1-\omega_0^2\tau^2}\,>0</math>, ''A'' et ''B'' étant deux constantes dépendant des conditions initiales. Il n'y a pas d'oscillations amorties, le système revient exponentiellement vers son état de repos.

== Oscillateurs mécaniques ==
=== Oscillations de translation ===

On peut mettre en oscillation une masse soumise à l'action d'un ressort. On peut suivant les cas, réaliser des oscillations verticales ou des oscillations horizontales (en utilisant un dispositif permettant de minimiser les frottements sur le support).

Dans les deux cas, les oscillations sont harmoniques : la fonction du temps <math>x(t)</math> de la position de la masse de part et d'autre de la position d'équilibre (statique) est une fonction sinus.
La période est indépendante de l’amplitude (isochronisme des oscillations) : elle ne dépend que de l'inertie du système (masse ''m'') et de la caractéristique de la force de rappel (constante de raideur ''k'' du ressort) :
<math> T = 2\pi\sqrt\frac{m}{k} .</math>

La constante de raideur ''k'' est exprimée en N/m. Pour ''k''&nbsp;=&nbsp;1&nbsp;N/m, il faudrait un Newton pour allonger le ressort d'un mètre.
En réalité, la période dépend aussi de la masse du ressort, qui est négligée dans la formule ci-dessus. Voir l'article détaillé : [[système masse-ressort]], pour une théorie plus complète.

Si les oscillations sont amorties par une force de frottement fluide (type visqueux à faible vitesse, force en <math>-\alpha v</math>), l'équation différentielle du mouvement peut s'écrire :

<math> m\ddot{x} + \alpha\dot{x} + k x = 0 .</math>

=== Oscillations de rotation ===
[[Image:torsion.gif|vignette|gauche|Schéma d'un pendule de torsion]]
{{Article détaillé|pendule de torsion}}

Le dispositif est constitué d'une barre fixée à un support par l'intermédiaire d'un fil de torsion : ce fil d'acier exerce un couple de rappel proportionnel à l'angle de torsion qu'on lui impose : <math>-C\theta\,</math>. Sur la barre on peut positionner deux masselottes de façon symétrique de façon à modifier le moment d'inertie.

La période est indépendante de l’amplitude (isochronisme des oscillations). Elle est donnée par la relation ci-dessous où ''J'' désigne le [[moment d'inertie]] de la barre munie des masselottes.

<math> T = 2\pi\sqrt\frac{J}{C} .</math>

Équation différentielle :
<math> J\ddot{\theta} + \alpha\dot{\theta} + C \theta = 0 .</math>

== Oscillateurs électriques ==
=== Circuit LC ===
Un '''[[circuit LC]]''' en [[électrocinétique]] est un circuit théorique comportant une [[bobine (électricité)|bobine]] idéale : parfaitement inductive ([[inductance]] ''L'' et résistance nulle) et un [[Condensateur (électricité)|condensateur]] (capacité ''C''). Les deux [[dipôle électrique|dipôles]] sont en outre totalement linéaires, ce qui est le cas avec des bobines à noyau d'air mais qui ne l'est pas pour de bobines avec un noyau [[ferromagnétique]].

Un tel circuit se comporte alors comme un oscillateur dont la période propre est :
<math>T_0 = 2\pi \sqrt{LC}.</math>

De la loi des mailles : <math> u_L + u_C = 0 \,</math>, et des équations caractéristiques des deux dipôles : <math> u_L = L \, \frac{\mathrm{d}i}{\mathrm{d}t}</math> et <math>i = C \, \frac{\mathrm{d}u_C}{\mathrm{d}t}</math>, on déduit :
:<math> L\,C\,\frac {\mathrm{d}^2u_C}{\mathrm{d}t^2}+ u_C = 0 </math>, avec <math>  u_C = \frac{q}{C} \,.</math>

L'équation différentielle peut donc s'écrire :
:<math> L\, \ddot{q} + \frac{q}{C} = 0 .</math>

=== Circuit RLC ===
Dans un circuit LC réel, on ne peut s'affranchir de la résistance. Celle-ci dissipe de l'énergie par effet Joule. Dans ce cas, l'équation différentielle qui régit les oscillations (amorties) peut s'écrire :
<math> L\ddot{q} + R\dot{q} + q /C = 0 .</math>

Remarque : on peut entretenir les oscillations grâce à un montage dit à résistance négative.

== Analogie électro-mécanique ==
Les oscillations mécaniques avec amortissement fluide, et les oscillations électriques d'un circuit RLC conduisent à deux équations différentielles du second ordre formellement identiques.

<math> \ddot{z} + 2\lambda \dot{z} + \omega_0^2 z = 0 .</math>

{| class="wikitable"
! Oscillateur générique
! RLC
! Masse soumise à un ressort
|-
|   
|<math>V</math> = tension
|<math>F</math> = Force
|-
|<math> \ddot{s} + 2\lambda \dot{s} + \omega_0^2 s = 0 </math>
|<math> V_L + V_R + V_C = L\frac{\mathrm di}{\mathrm dt} + R i + \frac{q}{C} = 0 </math>
|<math> ma + F_{fr} + F_R = m\ddot{x} + \alpha \dot{x} + k x = 0 </math>
|-
| <math>s</math>
| <math>q </math> = charge électrique
| <math>x</math> = déplacement
|-{{ligne grise}}
| <math>\dot{s}</math>
| <math>\dot{q} = i </math> = intensité
| <math>\dot{x}</math> = vitesse
|-
| <math>\ddot{s}</math> 
|<math>\ddot{q} = \frac{\mathrm{d}i}{\mathrm{d}t}</math>
| <math>\ddot{x}</math> = accélération
|-
| <math> \beta</math>
| <math>L</math> = inductance propre
| ''m'' = masse du mobile
|-{{ligne grise}}
| <math> \rho </math>
| ''R'' = résistance
| <math> \alpha </math> = coef de frottement
|-
| <math> \gamma </math>
| <math>\frac{1}{C}</math> = inverse de la capacité
| ''k'' = constante de raideur
|-{{ligne grise}}
| <math> T = 2\pi\sqrt\frac{\beta}{\gamma} </math> = période propre
| <math> T = 2\pi\sqrt{LC} </math> = période propre
| <math> T = 2\pi\sqrt\frac{m}{k} </math> = période propre
|-
| <math> f </math>
|  <math> P = RI^2 </math> : effet Joule
|  <math> f = \alpha  \dot{x} </math> : force de frottement
|-{{ligne grise}}
|<math> Q  </math> = facteur de qualité
|  <math> Q = \frac{1}{R}\sqrt\frac{L}{C}</math> 
|   <math> Q = \frac{1}{\alpha}\sqrt{mk}</math>
|}

== Oscillateurs électriques entretenus ==

On obtient des oscillations électriques entretenues en associant un élément amplificateur, à un circuit LC.''''' L'énergie est prélevée sur le circuit LC, amplifiée puis réinjectée sur le même circuit''''', tout cela par des couplages adéquats. L'oscillation a lieu lorsque les pertes du circuit LC sont compensées par l'énergie réinjectée par l'amplificateur. L'étude théorique peut être réalisée en considérant deux quadripoles bouclés l'un sur l'autre: Le quadripôle amplificateur et le circuit LC considéré comme un quadripôle filtre.

Différents montages pratiques d'oscillateurs ont été proposés : [[oscillateur Colpitts]], [[oscillateur Clapp]], [[oscillateur Hartley]], [[oscillateur de Pierce]].

Ces oscillateurs entretenus peuvent inclure un [[quartz (électronique)|quartz]], qui permet d'obtenir une meilleure stabilité de fréquence.
Les oscillateurs entretenus sont très utilisés en électronique par exemple, pour les horloges des circuits numériques, ou pour les pilotes des appareils de télécommunications.

== Notes et références ==
{{Références}}

== Articles connexes ==
* [[Rétroaction]]
* [[Oscillateur]]
* [[Systèmes oscillants à un degré de liberté]]
* [[Oscillateur harmonique quantique]]
* [[Oscillations amorties]]

== Bibliographie ==
* [[Vladimir Damgov]], ''Nonlinear and parametric phenomena. Applications in radiometric and mechanical systems,'' World Scientific, Series on Nonlinear Sciences, 2004.
* Leonard Meirovitch, ''Elements of vibration analysis'', 2nd ed., Mc Graw Hill, 1986 {{ISBN|0-070-41342-8}}.
* {{Landau|tome 1}}

{{Portail|physique}}

[[Catégorie:Système oscillant]]