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.
This paper is concerned with a Taylor series-based continuation algorithm, the so-called Asymptotic Numerical Method (ANM). It describes a generic continuation procedure to apply the ANM principle at best, in other words, that presents a high level of genericity without paying the price of this genericity by low computing performances. The way to quadratically recast a system of equation is now part of the method itself, and the way to handle elementary transcendental function is detailed with great attention. A sparse tensorial formalism is introduced for the internal representation of the system, which, when combined with a block condensation technique, provides a good computational efficiency of the ANM. Three examples are developed to show the performance and the versatility of the implementation of the continuation tool. Its robustness and its accuracy are explored. Finally, the potentiality of this method for complex nonlinear finite element analysis is enlightened by treating 2D elasticity problems with geometrical nonlinearities. KEYWORDS asymptotic numerical method, continuation, finite element method, nonlinear systems, quadratic recast, Taylor series *The term "numerical" stands in the name of the method because of this finite element discretization. Int J Numer Methods Eng. 2019;119:261-280.wileyonlinelibrary.com/journal/nme
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.