A split-frequency harmonic balance method for nonlinear oscillators with multi-harmonic forcing

Abstract A new harmonic balance method (HBM) is presented for accurately computing the periodic responses of a nonlinear sdof oscillator with multi-harmonic forcing and non-expansible nonlinearities. The presence of multi-harmonic forcing requires a large number of solution harmonics with a substantial increase in computational demand for either the conventional or the incremental HBM. In this method, the oscillator equation-error is first defined in terms of two functions (originally proposed for obtaining free-vibration periods in: R.E. Mickens, Iteration procedure for determining approximate solutions to nonlinear oscillator equations, Journal of Sound and Vibration 116 (1987) 185–187; and more recently: R.E. Mickens, A Generalised iteration procedure for calculating approximations to periodic solutions of “truly nonlinear oscillations”, Journal of Sound and Vibration 287 (2005) 1045–1051). A Fourier series solution is assumed, in which the total number of harmonics is fixed by the chosen discrete-time interval—this series is split into two partial sums nominally associated with either low-frequency or high-frequency harmonics. By exploiting a convergence property of the equation-error functions, the total solution is obtained in a new iterative scheme in which the low-frequency components are computed via a conventional HBM using a small number of algebraic equations, whereas the high frequency components are obtained in a separate step by updating. By gradually increasing the number of harmonics in the low-frequency group, the equation-error can be progressively reduced. Efficient use is made of FFT-based algebraic equation generation which allows an important class of non-expansible nonlinearities to be handled. The proposed method is tested on a Duffing-type oscillator, and an oscillator with a non-expansible 7th power stiffness term, where in both cases up to 24 component multi-harmonic forcing is applied. As a comparison, a conventional HBM is also used on the Duffing model in which the algebraic equations are generated in symbolic form to totally avoid errors from entering the formulation through complicated expansion of the cubic stiffness term (as in: I. Senjanovic, Harmonic analysis of nonlinear oscillations of cubic dynamical systems, Journal of Ship Research 38 (3) (1994) 225–238; and in: A. Raghothama, S. Narayanan, Periodic response and chaos in nonlinear systems with parametric excitation and time delay, Nonlinear dynamics 27 (2002) 341–365). The paper shows that in obtaining period-1 solutions, the computational accuracy and efficiency of the proposed method is very good.

[1]  Y. K. Cheung,et al.  Amplitude Incremental Variational Principle for Nonlinear Vibration of Elastic Systems , 1981 .

[2]  Soon-Yi Wu,et al.  Incremental harmonic balance method with multiple time scales for aperiodic vibration of nonlinear systems , 1983 .

[3]  J. C. Peyton Jones,et al.  Polyharmonic Balance Analysis of Nonlinear Ship Roll Response , 2004 .

[4]  T. D. Burton,et al.  On the Steady State Response and Stability of Non-Linear Oscillators Using Harmonic Balance , 1993 .

[5]  Wanda Szemplińska-Stupnicka.  Secondary resonances and approximate models of routes to chaotic motion in non-linear oscillators , 1987 .

[6]  Keith Worden.  On jump frequencies in the response of the Duffing oscillator , 1996 .

[7]  X. Wei,et al.  RESPONSE OF A DUFFING OSCILLATOR TO COMBINED DETERMINISTIC HARMONIC AND RANDOM EXCITATION , 2001 .

[8]  D. Ewins,et al.  The Harmonic Balance Method with arc-length continuation in rotor/stator contact problems , 2001 .

[9]  A. H. Nayfeh,et al.  Calculation of the jump frequencies in the response of s.d.o.f. non-linear systems , 2002 .

[10]  Ahmet Kahraman,et al.  Non-linear dynamic analysis of a multi-mesh gear train using multi-term harmonic balance method: sub-harmonic motions , 2005 .

[11]  I. Çankaya,et al.  Generalized Harmonic Analysis of Nonlinear Ship Roll Dynamics , 1996 .

[12]  Y. K. Cheung,et al.  Incremental time—space finite strip method for non‐linear structural vibrations , 1982 .

[13]  P. Orkwis,et al.  Adaptive harmonic balance method for nonlinear time-periodic flows , 2004 .

[14]  Leon O. Chua,et al.  Algorithms for computing almost periodic steady-state response of nonlinear systems to multiple input frequencies , 1981 .

[15]  Ronald E. Mickens,et al.  A generalized iteration procedure for calculating approximations to periodic solutions of “truly nonlinear oscillators” , 2005 .

[16]  Kay Hameyer,et al.  Frequency domain finite element approximations for saturable electrical machines under harmonic driving conditions , 1999 .

[17]  John E. T. Penny,et al.  The Accuracy of Jump Frequencies in Series Solutions of the Response of a Duffing Oscillator , 1994 .

[18]  S. Narayanan,et al.  Chaotic Oscillations of a Square Prism in Fluid Flow , 1993 .

[19]  S. Narayanan,et al.  Periodic Response and Chaos in Nonlinear Systems with Parametric Excitation and Time Delay , 2002 .

[20]  Ivo Senjanović.  HARMONIC ANALYSIS OF NONLINEAR OSCILLATIONS OF CUBIC DYNAMICAL SYSTEMS , 1994 .

[21]  Byeong Soo Kim,et al.  Dynamic analysis of harmonically excited non-linear structure system using harmonic balance method , 2001 .

[22]  B. Wu,et al.  MODIFIED MICKENS PROCEDURE FOR CERTAIN NON-LINEAR OSCILLATORS , 2002 .

[23]  Ronald E. Mickens,et al.  Iteration procedure for determining approximate solutions to non-linear oscillator equations , 1987 .

[24]  K. Huseyin,et al.  On the application of IHB technique to the analysis of non-linear oscillations and bifurcations , 1998 .