A purely frequency based Floquet-Hill formulation for the efficient stability computation of periodic solutions of ordinary differential systems

Abstract

Since the founding theory established by G. Floquet more than a hundred years ago, computing the stability of periodic solutions has given rise to various numerical methods, mostly depending on the way the periodic solutions are themselves determined, either in the time domain or in the frequency domain. In this paper, we address the stability analysis of branches of periodic solutions that are computed by combining a pure Harmonic Balance Method (HBM) with an Asymptotic Numerical Method (ANM). HBM is a frequency domain method for determining periodic solutions under the form of Fourier series and ANM is continuation technique that relies on high order Taylor series expansion of the solutions branches with respect to a path parameter. It is well established now that this HBM-ANM combination is efficient and reliable, provided that the system of ODE is first of all recasted with quadratic nonlinearities, allowing an easy manipulation of both the Taylor and the Fourier series. In this context, Hill's method, a frequency domain version of Floquet theory, is revisited so as to become a by-product of the HBM applied to a quadratic system, allowing the stability analysis to be implemented in an elegant way and with good computing performances. The different types of stability changes of periodic solutions are all explored and illustrated through several academic examples, including systems that are autonomous or not, conservative or not, free or forced.