0
Journal of Computational and Nonlinear Dynamics | Volume 1 | Issue 3 | Research Paper
RESEARCH PAPERS

Subharmonic-Response Computation and Stability Analysis for a Nonlinear Oscillator Using a Split-Frequency Harmonic Balance Method

[+] Author and Article Information
J. F. Dunne

Department of Engineering and Design, School of Science and Technology,  The University of Sussex, Falmer, Brighton, BN1 9QT, UKj.f.dunne@sussex.ac.uk

J. Comput. Nonlinear Dynam 1(3), 221-229 (Feb 20, 2006) (9 pages) doi:10.1115/1.2198875 History: Received September 05, 2005; Revised February 20, 2006
Article
References
Figures
Tables
Citing Articles

A split-frequency harmonic balance method (SF-HBM) is developed to obtain subharmonic responses of a nonlinear single-degree-of-freedom oscillator driven by periodic excitation. This method is capable of generating highly accurate periodic solutions involving a large number of solution harmonics. Responses at the excitation period, or corresponding multiples (such as period 2 and period 3), can be readily obtained with this method, either in isolation or as combinations. To achieve this, the oscillator equation error is first expressed in terms of two Mickens functions, where the assumed Fourier series solution is split into two groups, nominally associated with low-frequency or high-frequency harmonics. The number of low-frequency harmonics remains small compared to the number of high-frequency harmonics. By exploiting a convergence property of the equation-error functions, accurate low-frequency harmonics can be obtained in a new iterative scheme using a conventional harmonic balance method, in a separate step from obtaining the high-frequency harmonics. The algebraic equations (needed in the HBM part of the method) are generated wholly numerically via a fast Fourier transform, using a discrete-time formulation to include inexpansible nonlinearities. A nonlinear forced-response stability analysis is adapted for use with solutions obtained with this SF-HBM. Period-3 subharmonic responses are obtained for an oscillator with power-law nonlinear stiffness. The paper shows that for this type of oscillator, two qualitatively different period-3 subharmonic response branches can be obtained across a broad frequency range. Stability analysis reveals, however, that for an increasingly stiff model, neither of these subharmonic branches are stable.
FIGURES IN THIS ARTICLE
<>
Copyright © 2006 by American Society of Mechanical Engineers
Your Session has timed out. Please sign back in to continue.

References

1
Malatkar, P., and Nayfeh, A. H., 2002, “Calculation of the Jump Frequencies in the Response of S.D.O.F. Non-Linear Systems,” J. Sound Vib., 254 (5), pp. 1005–1011.
2
Nayfeh, A. H., and Mook., D. T., 1979, "Nonlinear Oscillations", Wiley, New York.
3
Worden, K., and Manson, G., 2005, “A Volterra Series Approximation to the Coherence of the Duffing Oscillator,” J. Sound Vib., 286 , pp. 529–547.
4
Lau, S. L., Cheung, Y. K., and Wu, S. Y., 1983, “Incremental Harmonic Balance Method With Multiple Time Scales for Aperiodic Vibration of Nonlinear Systems,” ASME J. Appl. Mech., 50 , pp. 871–876.
5
Hamdan, M. N., and Burton, T. D., 1993, “On the Steady State Response and Stability of Non-Linear Oscillators Using Harmonic Balance,” J. Sound Vib., 166 (2), pp. 255–266.
6
Worden, K., 1996, “On Jump Frequencies in the Response of the Duffing Oscillator,” J. Sound Vib., 198 , pp. 522–525.
7
Sensoy, S., and Huseyin, K., 1998, “On the Application of IHB Technique to the Analysis of Non-Linear Oscillations and Bifurcations,” J. Sound Vib., 215 (1), pp. 35–46.
8
Xu, Y. X., Bao, W. B., Schiehlen, W., and Hu, H. Y., 2001, “A 1/3 Pure Subharmonic Solution and Transient Process for the Duffing’s Equation,” Appl. Math. Mech., 22 (5), pp. 586–592.
9
Narayanan, S., and Jayaraman, K., 1993, “Chaotic Oscillations of a Square Prism in Fluid Flow,” J. Sound Vib., 166 (1), pp. 87–101.
10
Senjanović, I., 1994, “Harmonic Analysis of Nonlinear Oscillations of Cubic Dynamical Systems,” J. Ship Res., 38 (3), pp. 225–238.
11
Peyton Jones, J. C., and Çankaya, I., 1996, “Generalised Harmonic Analysis of Nonlinear Ship Roll Dynamics,” J. Ship Res., 40 (4), pp. 316–325.
12
von Groll, G., and Ewins, D. J., 2001, “The Harmonic Balance Method With Arc-length Continuation in Rotor/Stator Contact Problems,” J. Sound Vib.  [CrossRef], 241 , pp. 223–233.
13
Moon, B. Y., Kang, B. S., and Kim, B. S., 2001, “Dynamic Analysis of Harmonically Excited Non-Linear Structure System Using Harmonic Balance Method,” KSME Int. J., 15 (11), pp. 1507–1516.
14
Petrov, E. P., and Ewins, D. J., 2003, “Analytical Formulation of Friction Interface for Analysis of Nonlinear Multi-Harmonic Vibrations of Bladed Disks,” ASME J. Turbomach.  [CrossRef], 125 , pp. 364–371.
15
Raghothama, A., and Narayanan, S., 2002, “Periodic Response and Chaos in Nonlinear Systems With Parametric Excitation and Time Delay,” Nonlinear Dyn., 27 , pp. 341–365.
16
Al-shyyab, A., and Kahraman, A., 2005, “Nonlinear Dynamic Analysis of a Multi-Mesh Gear Train Using Multi-Term Harmonic Balance Method: Sub-Harmonic Motions,” J. Sound Vib., 279 , pp. 417–451.
17
Maple, R. C., King, P. I., Orkwis, P. D., and Wolff, J. M., 2004, “Adaptive Harmonic Balance Method for Nonlinear Time-Periodic Flows,” J. Comput. Phys., 193 , pp. 620–641.
18
Mickens, R. E., 1987, “Iteration Procedure for Determining Approximate Solutions to Nonlinear Oscillator Equations,” J. Sound Vib., 116 , pp. 185–187.
19
Mickens, R. E., 2005, “A Generalised Iteration Procedure for Calculating Approximations to Periodic Solutions of “truly Nonlinear Oscillations,” J. Sound Vib., 287 , pp. 1045–1051.

Figures

Grahic Jump Location
+Phase-plane trajectories obtained with the split-frequency HBM for the oscillator model with power-law stiffness for Ω∕ωn=5 with m=5, m=9, and m=13
View Large  |  Save Figure  |  Download Slide (.ppt)
Phase-plane trajectories obtained with the split-frequency HBM for the oscillator model with power-law stiffness for Ω∕ωn=5 with m=5, m=9, and m=13
Figure 5

Phase-plane trajectories obtained with the split-frequency HBM for the oscillator model with power-law stiffness for Ω∕ωn=5 with m=5, m=9, and m=13

Grahic Jump Location
+Equation error rms as a function of iteration number (shown as •) with 4 iterations for each new low-frequency harmonic (total number=8), fourth iteration shown overlaid with ∘): (a) m=5, (b) m=9, (c) m=13; (d), (e), and (f), show corresponding absolute magnitudes of the total 64 response harmonics (low-plus high-frequency)
View Large  |  Save Figure  |  Download Slide (.ppt)
Equation error rms as a function of iteration number (shown as •) with 4 iterations for each new low-frequency harmonic (total number=8), fourth iteration shown overlaid with ∘): (a) m=5, (b) m=9, (c) m=13; (d), (e), and (f), show corresponding absolute magnitudes of the total 64 response harmonics (low-plus high-frequency)
Figure 6

Equation error rms as a function of iteration number (shown as •) with 4 iterations for each new low-frequency harmonic (total number=8), fourth iteration shown overlaid with ∘): (a) m=5, (b) m=9, (c) m=13; (d), (e), and (f), show corresponding absolute magnitudes of the total 64 response harmonics (low-plus high-frequency)

Grahic Jump Location
+Subharmonic frequency response amplitudes as a function of frequency ratio Ω∕ωn for an oscillator with power-law stiffness parameter m=7 showing different solution branches
View Large  |  Save Figure  |  Download Slide (.ppt)
Subharmonic frequency response amplitudes as a function of frequency ratio Ω∕ωn for an oscillator with power-law stiffness parameter m=7 showing different solution branches
Figure 8

Subharmonic frequency response amplitudes as a function of frequency ratio Ω∕ωn for an oscillator with power-law stiffness parameter m=7 showing different solution branches

Grahic Jump Location
+Subharmonic frequency response amplitudes as a function of frequency ratio Ω∕ωn for an oscillator with power-law stiffness parameter m=11 showing different solution branches
View Large  |  Save Figure  |  Download Slide (.ppt)
Subharmonic frequency response amplitudes as a function of frequency ratio Ω∕ωn for an oscillator with power-law stiffness parameter m=11 showing different solution branches
Figure 9

Subharmonic frequency response amplitudes as a function of frequency ratio Ω∕ωn for an oscillator with power-law stiffness parameter m=11 showing different solution branches

Grahic Jump Location
+(a) Phase-plane trajectory for the stable subharmonic response of the Duffing oscillator (m=3); (b) equation error rms as a function of iteration number (shown as •) with four iterations for each new low-frequency harmonic (total number=8), fourth iteration shown overlaid with ●); (c) absolute magnitudes of the total 64 response harmonics (low- plus high-frequency); and (d) complex eigenvalues associated with the (stable) subharmonic response
View Large  |  Save Figure  |  Download Slide (.ppt)
(a) Phase-plane trajectory for the stable subharmonic response of the Duffing oscillator (m=3); (b) equation error rms as a function of iteration number (shown as •) with four iterations for each new low-frequency harmonic (total number=8), fourth iteration shown overlaid with ●); (c) absolute magnitudes of the total 64 response harmonics (low- plus high-frequency); and (d) complex eigenvalues associated with the (stable) subharmonic response
Figure 1

(a) Phase-plane trajectory for the stable subharmonic response of the Duffing oscillator (m=3); (b) equation error rms as a function of iteration number (shown as •) with four iterations for each new low-frequency harmonic (total number=8), fourth iteration shown overlaid with ●); (c) absolute magnitudes of the total 64 response harmonics (low- plus high-frequency); and (d) complex eigenvalues associated with the (stable) subharmonic response

Grahic Jump Location
+(a) Phase-plane trajectory for the unstable subharmonic response of the Duffing oscillator; (b) equation error rms as a function of iteration number (shown as • symbol) with four iterations for each new low-frequency harmonic (total number=8), fourth iteration shown overlaid with ●); (c) absolute magnitudes of the total 64 response harmonics (low- plus high-frequency); and (d) complex eigenvalues associated with the (unstable) subharmonic
View Large  |  Save Figure  |  Download Slide (.ppt)
(a) Phase-plane trajectory for the unstable subharmonic response of the Duffing oscillator; (b) equation error rms as a function of iteration number (shown as • symbol) with four iterations for each new low-frequency harmonic (total number=8), fourth iteration shown overlaid with ●); (c) absolute magnitudes of the total 64 response harmonics (low- plus high-frequency); and (d) complex eigenvalues associated with the (unstable) subharmonic
Figure 2

(a) Phase-plane trajectory for the unstable subharmonic response of the Duffing oscillator; (b) equation error rms as a function of iteration number (shown as • symbol) with four iterations for each new low-frequency harmonic (total number=8), fourth iteration shown overlaid with ●); (c) absolute magnitudes of the total 64 response harmonics (low- plus high-frequency); and (d) complex eigenvalues associated with the (unstable) subharmonic

Grahic Jump Location
+Contour plot of the initial-final-state-error surface for the Duffing oscillator for Ω∕ωn=5 computed via time-domain integration showing initial conditions for periodic solutions obtained with a systematic grid search (highlighted with an asterisk)
View Large  |  Save Figure  |  Download Slide (.ppt)
Contour plot of the initial-final-state-error surface for the Duffing oscillator for Ω∕ωn=5 computed via time-domain integration showing initial conditions for periodic solutions obtained with a systematic grid search (highlighted with an asterisk)
Figure 3

Contour plot of the initial-final-state-error surface for the Duffing oscillator for Ω∕ωn=5 computed via time-domain integration showing initial conditions for periodic solutions obtained with a systematic grid search (highlighted with an asterisk)

Grahic Jump Location
+Phase-plane trajectories obtained by time domain integration for the Duffing oscillator with Ω∕ωn=5 associated with the initial conditions found in Fig. 3
View Large  |  Save Figure  |  Download Slide (.ppt)
Phase-plane trajectories obtained by time domain integration for the Duffing oscillator with Ω∕ωn=5 associated with the initial conditions found in Fig. 3
Figure 4

Phase-plane trajectories obtained by time domain integration for the Duffing oscillator with Ω∕ωn=5 associated with the initial conditions found in Fig. 3

Grahic Jump Location
+Complex eigenvalues associated with the subharmonic response of the oscillator with power-law stiffness for Ω∕ωn=5 (a) m=7 and (b) m=9
View Large  |  Save Figure  |  Download Slide (.ppt)
Complex eigenvalues associated with the subharmonic response of the oscillator with power-law stiffness for Ω∕ωn=5 (a) m=7 and (b) m=9
Figure 7

Complex eigenvalues associated with the subharmonic response of the oscillator with power-law stiffness for Ω∕ωn=5 (a) m=7 and (b) m=9

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Dynamic Behavior in a Singular Delayed Bioeconomic Model
International Conference on Instrumentation, Measurement, Circuits and Systems (ICIMCS 2011)
Parallel Computing Using Multilevel Hypergraph Partitioning for Three-Dimensional Scramjet
International Conference on Mechanical and Electrical Technology, 3rd, (ICMET-China 2011), Volumes 1–3
Aerodynamic Performance Analysis
Axial-Flow Compressors> Chapter 9
Introduction and Definitions
Handbook on Stiffness & Damping in Mechanical Design> Chapter 1
Modes of Loading and Stiffness of Structural Components
Handbook on Stiffness & Damping in Mechanical Design> Chapter 2
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In