In this Note, we present a harmonic-based numerical method to determine the local stability of periodic solutions of dynamical systems. Based on the Floquet theory and the Fourier series expansion (Hill method), we propose a simple strategy to sort the relevant physical eigenvalues among the expanded numerical spectrum of the linear periodic system governing the perturbed solution. By mixing the harmonic-balance method and asymptotic numerical method continuation technique with the developed Hill method, we obtain a purely-frequency based continuation tool able to compute the stability of the continued periodic solutions in a reduced computation time. To validate the general methodology, we investigate the dynamical behavior of the forced Duffing oscillator with the developed continuation technique.

Dans cette Note, nous présentons une méthode numérique fréquentielle pour déterminer la stabilité des solutions périodiques d'un système dynamique. La méthode, basée sur la théorie de Floquet et le développement en série de Fourier (méthode de Hill), consiste à extraire les valeurs propres physiques de l'ensemble des valeurs propres numériques du système perturbé étendu dans le domaine fréquentiel. En combinant alors la méthode de l'équilibrage harmonique et la méthode asymptotique numérique avec la précédente méthode de Hill, on obtient un outils de continuation purement fréquentiel où le calcul de la stabilité des solutions suivies est quasiment immédiat. Afin de valider la méthode, nous appliquons la méthode de continuation à un oscillateur de Duffing forcé.