Technical Brief

A New Spatial and Temporal Harmonic Balance Method for Obtaining Periodic Steady-State Responses of a One-Dimensional Second-Order Continuous System

[+] Author and Article Information
X. F. Wang

Department of Mechanical Engineering,
University of Maryland,
Baltimore County,
1000 Hilltop Circle,
Baltimore, MD 21250
e-mail: xuefeng1@umbc.edu

W. D. Zhu

Fellow ASME
Division of Dynamics and Control,
School of Astronautics,
Harbin Institute of Technology,
P.O. Box 137,
Harbin 150001, China;
Department of Mechanical Engineering,
University of Maryland Baltimore County,
1000 Hilltop Circle,
Baltimore, MD 21250
e-mail: wzhu@umbc.edu

1Corresponding author.

Manuscript received October 21, 2015; final manuscript received June 23, 2016; published online October 18, 2016. Assoc. Editor: Alexander F. Vakakis.

J. Appl. Mech 84(1), 014501 (Oct 18, 2016) (6 pages) Paper No: JAM-15-1567; doi: 10.1115/1.4034011 History: Received October 21, 2015; Revised June 23, 2016

A new spatial and temporal harmonic balance (STHB) method is developed for obtaining periodic steady-state responses of a one-dimensional second-order continuous system. The spatial harmonic balance procedure with a series of sine and cosine basis functions can be efficiently conducted by the fast discrete sine and cosine transforms, respectively. The temporal harmonic balance procedure with basis functions of Fourier series can be efficiently conducted by the fast Fourier transform (FFT). In the STHB method, an associated set of ordinary differential equations (ODEs) of a governing partial differential equation (PDE), which is obtained by Galerkin method, does not need to be explicitly derived, and complicated calculation of a nonlinear term in the PDE can be avoided. The residual and the exact Jacobian matrix of an associated set of algebraic equations that are temporal harmonic balanced equations of the ODEs, which are used in Newton–Raphson method to iteratively search a final solution of the PDE, can be directly obtained by STHB procedures for the PDE even if the nonlinear term is included. The relationship of Jacobian matrix and Toeplitz form of the system matrix of the ODEs provides an efficient and convenient way to stability analysis for the STHB method; bifurcations can also be indicated. A complex boundary condition of a string with a spring at the boundary can be handled by the STHB method in combination with the spectral Tau method.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.


Meirovitch, L. , 1997, Principles and Techniques of Vibrations, Vol. 1, Prentice Hall, Upper Saddle River, NJ.
Szemplińska-Stupnicka, W. , 1978, “ The Generalized Harmonic Balance Method for Determining the Combination Resonance in the Parametric Dynamic Systems,” J. Sound Vib., 58(3), pp. 347–361. [CrossRef]
Tamura, H. , Tsuda, Y. , and Sueoka, A. , 1981, “ Higher Approximation of Steady Oscillations in Nonlinear Systems With Single Degree of Freedom: Suggested Multi-Harmonic Balance Method,” Bull. JSME, 24(195), pp. 1616–1625. [CrossRef]
Lau, S. , Cheung, Y. , and Wu, S. , 1982, “ A Variable Parameter Incrementation Method for Dynamic Instability of Linear and Nonlinear Elastic Systems,” ASME J. Appl. Mech., 49(4), pp. 849–853. [CrossRef]
Ferri, A. , 1986, “ On the Equivalence of the Incremental Harmonic Balance Method and the Harmonic Balance–Newton Raphson Method,” ASME J. Appl. Mech., 53(2), pp. 455–457. [CrossRef]
Ling, F. , and Wu, X. , 1987, “ Fast Galerkin Method and Its Application to Determine Periodic Solutions of Non-Linear Oscillators,” Int. J. Non-Linear Mech., 22(2), pp. 89–98. [CrossRef]
Guillen, J. , and Pierre, C. , 1999, “ An Efficient, Hybrid, Frequency-Time Domain Method for the Dynamics of Large-Scale Dry-Friction Damped Structural Systems,” IUTAM Symposium on Unilateral Multibody Contacts, Springer, Munich, Germany, pp. 169–178.
Wang, X. , and Zhu, W. , 2015, “ A Modified Incremental Harmonic Balance Method Based on the Fast Fourier Transform and Broyden's Method,” Nonlinear Dyn., 81(1–2), pp. 981–989. [CrossRef]
Hsu, C. , and Cheng, W.-H. , 1973, “ Applications of the Theory of Impulsive Parametric Excitation and New Treatments of General Parametric Excitation Problems,” ASME J. Appl. Mech., 40(1), pp. 78–86. [CrossRef]
Wereley, N. M. , 1990, “ Analysis and Control of Linear Periodically Time Varying Systems,” Ph.D. thesis, Massachusetts Institute of Technology, Cambridge, MA.
Chivukula, R. K. , and Reznik, Y. A. , 2011, “ Fast Computing of Discrete Cosine and Sine Transforms of Types VI and VII,” Proc. SPIE 8135, p. 812505.
Johnson, D. , 1996, “ Chebyshev Polynomials in the Spectral Tau Method and Applications to Eigenvalue Problems,” Report No. 198451.
Mote, C. , 1966, “ On the Nonlinear Oscillation of an Axially Moving String,” ASME J. Appl. Mech., 33(2), pp. 463–464. [CrossRef]


Grahic Jump Location
Fig. 2

Frequency-response curves for the case of strong nonlinearity with kd = 30; N = 5 and 10. Unstable solutions for N = 5 are indicated.

Grahic Jump Location
Fig. 1

Frequency-response curves for the case of weak nonlinearity with kd = 10; N = 5 and 10

Grahic Jump Location
Fig. 3

Frequency-response curves for cases of the string with fixed–fixed boundary conditions and fixed-spring boundary conditions; k0 = 10,000 for the large spring stiffness case and k0 = 10 for the small spring stiffness case



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
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