The Duffing equation containing a cubic nonlinearity is probably the most popular example of a nonlinear oscillator. For its harmonically excited, slightly damped, and softening version, stationary large amplitude solutions at subcritical excitation frequencies are obtained when standard semi-analytical methods like Harmonic Balance or Perturbation Analysis are applied. These solutions have the shape of a nose in the amplitude-frequency diagram. In prior work, it has been observed that these solutions may contain large errors and that high ansatz orders may be necessary when applying the Harmonic Balance or other semi-analytical methods to make them converge. Some of these solutions are observed to be asymptotically stable, while in most cases, they are unstable. The current paper aims to give a descriptive explanation for this behavior of the nose solutions, which is mainly related to the exact solution of the free undamped vibrations. Based on this, approximations of the nose solutions are calculated with a procedure combining properties of Perturbation Analysis and Harmonic Balance. Therein, the exact solution of the free undamped vibrations is taken as the zeroth approximation, while higher-order solution parts are calculated by balancing the harmonics, and the phase shift of the zeroth approximation is calculated by a residuum minimization. This method just requires the solution of a system of linear algebraic equations, while systems of nonlinear algebraic equations have to be solved in the case of directly applying Harmonic Balance.
Read full abstract