thumb|Une section de Poincaré de l'équation de Duffing forcée suggérant un comportement chaotique .
Léquation de Duffing (ou oscillateur de Duffing), du nom de (1861–1944), est une équation différentielle non linéaire du second ordre utilisée pour modéliser certains oscillateurs amortis et forcés. L'équation s'écrit
qui décrit le déplacement en fonction du temps .
L'équation décrit le mouvement d'un oscillateur amorti avec un potentiel plus complexe qu'un mouvement harmonique simple (cas correspondant à ) ; un modèle physique serait un pendule pesant où la raideur du ressort ne suivrait pas la loi de Hooke.
L'équation de Duffing est un exemple de système dynamique simple pouvant présenter un comportement chaotique, comme l'oscillateur de Van der Pol. Plus encore, le système de Duffing présente en réponse fréquentielle le phénomène de résonance de saut qui se rapproche d'un comportement d'hystérésis en fréquence.
Paramètres Les paramètres donnés dans l'équation caractérisent les différents effets : contrôle le taux d'amortissement, contrôle la raideur linéaire, contrôle le taux de non-linéarité dans la force restauratrice ; le cas correspondant à un oscillateur harmonique amorti et forcé simple, est l'amplitude de la force conductrice périodique, est la fréquence angulaire de cette force.
L'équation de Duffing peut être vue comme décrivant les oscillations d'une masse attachée à un ressort non linéaire et un amortisseur linéaire. La force restauratrice fournie par le ressort s'écrit .
Pour et , le ressort est dit raide ; si , on parle de ressort souple (toujours pour ). Les mêmes qualificatifs sont alors appliqués à l'équation de Duffing, selon les valeurs de (et )..
Le nombre de paramètres dans l'équation de Duffing peut être réduit à 2 par changement d'échelle. Un exemple consiste à changer et en. et en supposant (d'autres changements de variables sont possibles, selon les valeurs des paramètres ou l'importance portée à certains termes). Dans ce cas.:
avec
L'équation ne dépend plus que de trois coefficients et des deux conditions initiales.
Méthodes de résolution Il n'est pas possible de donner une solution exacte à l'équation de Duffing sous forme symbolique. Cependant, plusieurs méthodes d'approximation fonctionnent bien : un développement en série de Fourier permet d'avoir une description du mouvement pour un ordre de précision donné ; le terme , parfois appelé terme de Duffing, peut être approché aussi près que souhaité et le système est vu comme un oscillateur harmonique simple perturbé ; la méthode de Frobenius ; des méthodes de résolution numérique classique (méthode d'Euler, Runge-Kutta...) ; la méthode par homotopie est utilisée avec de bons résultats, même dans le cas de forte non-linéarité...
Dans le cas spécial de l'équation sans amortissement () et sans forçage (), on peut écrire la solution exacte avec les fonctions elliptiques de Jacobi..
Bornes de la solution pour l'oscillateur non forcé Oscillateur non amorti En multipliant l'équation de Duffing non amorti et non forcé (), avec , on obtient. :
où est une constante déterminée par les conditions initiales et
La substitution dans montre que le système est hamiltonien :
Quand et sont positifs, la solution est bornée:
avec le hamiltonien positif.
Oscillateur amorti De façon similaire, pour l'oscillateur amorti.
avec pour l'amortissement. Sans forçage, l'oscillateur de Duffing amorti va tendre vers un de ses points d'équilibre stable. Ces points d'équilibre, stable ou non, correspondent aux cas où . Si , le seul point d'équilibre stable est en ; si et , l'équilibre stable peut être atteint en .
Réponse fréquentielle right|thumb|upright=1.5|Réponse fréquentielle en fonction de pour l'équation de Duffing, avec et un amortissement à . Les parties en pointillés de la réponse fréquentielle sont instables. L'oscillateur de Duffing forcé avec non-linéarité forcée est décrite par l'équation différentielle ordinaire :
La réponse fréquentielle de cet oscillateur décrit l'amplitude d'une réponse en régime permanent de l'équation () à une fréquence donnée pour une excitation . Pour un oscillateur linéaire avec , la réponse fréquentielle est également linéaire. Cependant, pour un terme de Duffing non nul, la réponse fréquentielle devient non linéaire. Selon le type de non-linéarité, l'oscillateur de Duffing peut montrer une réponse fréquentielle raide, souple ou mixte. Par la méthode d'analyse par homotopie ou d'équilibrage harmonique, on peut dériver une équation de la réponse fréquentielle de la forme.:
Pour les paramètres de l'équation de Duffing, cette équation algébrique donne l'amplitude d'oscillation en régime permanent à une excitation fréquentielle donnée.
Sauts thumb|right|Sauts dans la réponse fréquentielle pour , , . Les variations de fréquence créent un hystérésis entre les points A, B, C et D. Pour certains ordres de valeurs des paramètres de l'équation de Duffing, la réponse fréquentielle peut ne plus une fonction à valeurs réelles de la fréquence de forçage . Pour un oscillateur à ressort raide ( et pris assez grand), la réponse fréquentielle penche vers les hautes fréquences, et vers les basses fréquences pour celui à ressort souple ( et pris assez grand). Le côté surplombant est instable et la réponse correspondante ne peut être maintenue sur une longue durée. Par conséquent, le phénomène de saut apparait : quand la fréquence angulaire est augmentée lentement (avec les autres paramètres fixés), l'amplitude de réponse chute violemment, si la fréquence angulaire est lentement réduite, l'amplitude fait un bond vers la branche haute de la réponse fréquentielle. Les deux sauts ne correspondent pas, créant ainsi un hystérésis dans le système dépendant de l'évolution de la fréquence.
Exemples
Certains exemples typiques de séries temporelles et portraits de phase de l'équation de Duffing, montrant l'apparition de sous-harmoniques par des bifurcations à doublement de période – ou des comportements chaotiques. L'amplitude forcée augmente de à . Les autres paramètres ont les valeurs et . Les conditions initiales sont et Les points rouges dans les portraits de phase sont aux temps pour un multiple entier de la période .
Références
Voir aussi Articles connexes Oscillateur de Van der Pol Théorie du chaos
Bibliographie
Liens externes Duffing oscillator sur Scholarpedia
Catégorie:Équation différentielle Catégorie:Théorie du chaos