Abstract
A power series expansion of the characteristic parameter is obtained in terms of the natural frequency, which is defined as an implicit function in a transcendental equation. It plays a central role in calculating the Fourier coefficient of exact solutions. A general algorithm is proposed to obtain the solution in the power series x=Σ bnyn to the equation Σ cnxn=y. The method is particularly effective when m (the modulus parameter of elliptic functions) approaches zero, which corresponds with the quasi-linear small oscillations in those systems with a hard, a soft or snap-through (half-swing mode) spring. The accuracy of existing (simple) approximate solutions to the autonomous Duffing equation, i.e., the perturbation, the averaging and the (two-term) harmonic balance methods, is examined by comparison with the present exact solutions.