A high dimensional harmonic balance approach for an aeroelastic airfoil with cubic restoring forces

This is a study of a two-dimensional airfoil including a cubic spring stiffness placed in an incompressible flow. A new formulation of the harmonic balance method is employed for the aeroelastic airfoil to investigate the amplitude and frequency of the limit cycle oscillations. The results are compared with the results from the classical harmonic balance approach and from the conventional time marching integration method.     

Localized Basis Function Method for Computing Limit Cycle Oscillations

An alternate approach to the standard harmonic balance method (based on Fourier transforms) is proposed. The proposed method begins with an idea similar to the harmonic balance method, i.e. to transform the initial set of differential equations of the dynamics to a set of discrete algebraic equations. However, as distinct from previous harmonic balance techniques, the proposed method uses a set of basis functions which are localized in time and are not necessarily sinusoidal. Also as distinct from previous harmonic balance methods, the algebraic equations obtained after the transformation of the differential equations of the dynamics are solved in the time domain rather than the frequency domain. Numerical examples are provided to demonstrate the performance of the method for autonomous and forced dynamics of a Van der Pol oscillator.     

The Simplest Resonance Capture Problem,Using Harmonic Balance Based Averaging

We study the dynamics of capture into, or escape from, resonance in a strongly nonlinear oscillator with weak damping and forcing, using harmonic balance based averaging (HBBA). This system provides the simplest example of resonance capture that we know of. The HBBA technique, here adapted to tackle nonlinear resonances, provides a harmonic balance assisted approximation to the underlying, asymptotically correct, averaged dynamics. Allowing the harmonic balance approximation makes a variety of systems analytically tractable which might otherwise be intractable. The evolution equations for amplitude and phase of oscillations are derived first. Restricting attention near the primary resonance, the slow flow equations are approximately averaged. The resulting flow transparently shows the stable and unstable primary resonant solutions, as well as the trajectories that get captured into resonance and the ones that escape. Good agreement with numerics is obtained, showing the utility of HBBA near resonance manifolds.     

An incremental method for limit cycle oscillations of an airfoil with an external store

This paper proposes an incremental method, which is based on the harmonic balance method, to analyze the nonlinear aeroelastic problem of an airfoil with an external store. The governing equations of limit cycle oscillations (LCOs) of the airfoil are deduced by the harmonic balancing procedure. Different from usual procedures, the harmonic balance equations are not solved directly but instead transformed into an equivalent minimization problem. The minimization problem is solved using the Levenberg–Marquardt method. Numerical examples show that the LCOs obtained by the presented method are in excellent agreement with numerical solutions. The bifurcation of the LCOs is further analyzed using the Floquet theory. It is found that the LCOs exhibit saddle-node, symmetry breaking and period-doubling bifurcations with the wind speed as control parameter. Compared with the harmonic balance method, the presented method has a wider convergence region and hence makes it easier to choose a proper initial guess for iterations.     

Recent advances and comparisons between harmonic balance and Volterra-based nonlinear frequency response analysis methods

Harmonic balance and Volterra-based analysis methods are well known, but the capabilities of these methods have been limited by significant issues of complexity which either constrain their application to relatively simple cases, or limit the accuracy of analysis in more complex cases. This study briefly summarizes recent results which effectively extend the capabilities of both harmonic balance and Volterra-based analysis by making complex analyses much more feasible. The new capabilities and performance of the two approaches are then evaluated and compared using benchmark case studies of a Duffing oscillator and a nonlinear automotive damper. The results offer new insights and lead to different conclusions on the relative merits of harmonic balance versus Volterra-based analysis relative to prior studies and similar benchmark analyses.     

Iterative harmonic balance for period-one rotating solution of parametric pendulum

In this study, an iterative method based on harmonic balance for the period-one rotation of parametrically excited pendulum is proposed. Based on the definition of the period-one rotating orbit, the exact form of the solution can be obtained using the Fourier series. An iterative harmonic balance process is proposed to estimate the coefficients in the exact solution form. The general formula for each iteration step is presented. The method is evaluated using two criteria, which are the system energy error and the global residual error. The performance of the proposed method is compared with the results from multiscale method and perturbation method. The numerical results obtained with the Dormand?CPrince method (ODE45 in MATLAB?) are used as the baseline of the evaluation.     

基于增量谐波平衡的参激系统非线性识别法

将增量谐波平衡法应用到非线性系统的建模和参数识别中,针对Mathieu-Duffing方程,推导了利用增量谐波平衡原理识别参数激励非线性系统参数的方法. 该方法改进了增量谐波平衡方法的推导过程,通过数值模拟对比研究了谐波平衡非线性识别(harmonic balance nonlinearity identification, HBNID)和增量谐波平衡非线性识别(incremental harmonic balance nonlinearity identification, IHBNID)的效果,验证了增量谐波平衡非线性识别的有效性. 结果表明,增量谐波平衡非线性识别的计算效率较高,计算精度和抗噪能力都优于谐波平衡非线性识别.      

The application of reduced space harmonic balance method for the nonlinear vibration problem in rotor dynamics

The reduced space harmonic balance method is demonstrated to find the maximum vibration responses of rotor systems. Within the reduced space SQP method, transition from the high dimensional optimization space to the desired reduced space is accomplished by resorting to the null space decomposition technique, resulting in the elimination of the harmonic balance constraints. Numerical examples of rotor systems are presented to show the applicability of the proposed methodology.     

An analytical approximate technique for a class of strongly non-linear oscillators

An analytical approximate technique for large amplitude oscillations of a class of conservative single degree-of-freedom systems with odd non-linearity is proposed. The method incorporates salient features of both Newton's method and the harmonic balance method. Unlike the classical harmonic balance method, accurate analytical approximate solutions are possible because linearization of the governing differential equation by Newton's method is conducted prior to harmonic balancing. The approach yields simple linear algebraic equations instead of non-linear algebraic equations without analytical solution. With carefully constructed iterations, only a few iterations can provide very accurate analytical approximate solutions for the whole range of oscillation amplitude beyond the domain of possible solution by the conventional perturbation methods or harmonic balance method. Three examples including cubic-quintic Duffing oscillators are presented to illustrate the usefulness and effectiveness of the proposed technique.     

Limit cycle oscillations of rectangular cantilever wings containing cubic nonlinearity in an incompressible flow

Limit cycle oscillations (LCO) as well as nonlinear aeroelastic analysis of rectangular cantilever wings with a cubic nonlinearity are investigated. Aeroelastic equations of a rectangular cantilever wing with two degrees of freedom in an incompressible potential flow are presented in the time domain. The harmonic balance method is modified to calculate the LCO frequency and amplitude for rectangular wings. In order to verify the derived formulation, flutter boundaries are obtained via a linear analysis of the derived system of equations for five different cases and compared with experimental data. Satisfactory results are gained through this comparison. The problem of finding the LCO frequency and amplitude is solved via applying the two methods discussed for two different cases with hardening cubic nonlinearities. The results from first-, third- and fifth-order harmonic balance methods are compared with the results of an exact numerical solution. A close agreement is obtained between these harmonic balance methods and the exact numerical solution of the governing aeroelastic equations. Finally, the nonlinear aeroelastic analysis of a rectangular cantilever wing with a softening nonlinearity is studied.     

On Constructing the Unique Solution for the Necking in a Hyper-Elastic Rod

In this paper, a novel approach which considers gradient effects and uses non-deforming boundary conditions is adopted to construct the unique solution for necking in a hyper-elastic rod. We study the problem of the large axially symmetric deformations of a rod composed of an incompressible Ogden’s hyper-elastic material subject to a tensile stress (or a given displacement) when its two ends are fixed to rigid bodies. The attention is on the class of energy functions for which the stress–strain curve in the case of the uniaxial tension has a peak and valley combination. A phase-plane analysis is introduced to study the qualitative behaviour of the solutions. Then, by using the non-deforming conditions at two ends, the solutions corresponding to trajectories in different phase planes are obtained. It turns out that the non-deforming conditions play an important role in selecting the solutions. Further, by converting the problem into a displacement-controlled problem, the unique solution is obtained. The engineering strain and engineering stress curve plotted from our solution exhibits two interesting phenomena: (i) After the stress reaches the peak value there is a sudden stress drop; (ii) Afterwards it is followed by a stress plateau. Some mathematical explanations on these two phenomena are then given.     

Nonlinear vibration of a two-mass system with nonlinear stiffnesses

An analytical approach is developed for nonlinear free vibration of a conservative, two-degree-of-freedom mass–spring system having linear and nonlinear stiffnesses. The main contribution of the proposed approach is twofold. First, it introduces the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation and, more significantly, the treatment a nonlinear differential system by linearization coupled with Newton’s method and harmonic balance method. New and accurate higher-order analytical approximate solutions for the nonlinear system are established. After solving the nonlinear differential equation, the displacement of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, this higher-order Newton–harmonic balance (NHB) method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation unlike the classical harmonic balance method which results in complicated algebraic equations requiring further numerical analysis. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and the classical harmonic balance methods. Two examples of nonlinear two-degree-of-freedom mass–spring system are analyzed and verified with published result, exact solutions and numerical integration data.     

Forced response of quadratic nonlinear oscillator: comparison of various approaches

The primary resonances of a quadratic nonlinear system under weak and strong external excitations are investigated with the emphasis on the comparison of different analytical approximate approaches. The forced vibration of snap-through mechanism is treated as a quadratic nonlinear oscillator. The Lindstedt-Poincaré method, the multiple-scale method, the averaging method, and the harmonic balance method are used to determine the amplitude-frequency response relationships of the steady-state responses. It is demonstrated that the zeroth-order harmonic components should be accounted in the application of the harmonic balance method. The analytical approximations are compared with the numerical integrations in terms of the frequency response curves and the phase portraits. Supported by the numerical results, the harmonic balance method predicts that the quadratic nonlinearity bends the frequency response curves to the left. If the excitation amplitude is a second-order small quantity of the bookkeeping parameter, the steady-state responses predicted by the second-order approximation of the LindstedtPoincaré method and the multiple-scale method agree qualitatively with the numerical results. It is demonstrated that the quadratic nonlinear system implies softening type nonlinearity for any quadratic nonlinear coefficients.     

Epitaxy of Binary Compounds and Alloys

The present work aims at constructing a theoretical framework within which to address the issues of morphological instabilities (one-dimensional step bunching and two-dimensional step meandering), alloying, and phase segregation in binary systems in the context of (physical or chemical) vapor deposition. The length scale of interest, although nanoscopic, is sufficiently large that the steps on a vicinal surface can be viewed as smooth curves and, correspondingly, the theory is a continuum one. In a departure from theories inaugurated by Burton, Cabrera, and Frank [The growth of crystals and the equilibrium structure of their surfaces. Phil. Trans. Roy. Soc. A 243 (1951) 299–358] the steps are endowed with a thermodynamic structure whose main ingredients are a step free-energy density and edge species chemical potentials. Moreover, crystal anisotropy, with its altering of the dynamics of steps and the associated morphological instabilities, is accounted for – in a manner consistent with the second law – both in the thermodynamic and kinetic properties of terraces and, more importantly, of steps. Additionally, in contrast with most of the literature on the subject (cf. [J. Krug, Introduction to step dynamics and step instabilities. In: A. Voigt (ed.) Multiscale Modeling in Epitaxial Growth. Birkhäusser, Berlin (2005)]), adsorption–desorption along the steps, bulk atomic diffusion, and chemical reactions (both on the terraces and along the step edges) are incorporated and coupled to the other mechanisms, e.g., terrace adatom diffusion and step attachment–detachment kinetics, whose interplay governs the evolution of steps on vicinal surfaces. Importantly, aided by the concept of configurational forces for which a separate balance law is postulated Configurational Forces as Basic Concepts of Continuum Physics. Springer, Berlin Heidelberg New York (2000)]), the proposed theory allows the steps to evolve away from local equilibrium thus contributing to a general treatment of the dynamics of steps. Finally, a specialization to the epitaxy of binary compounds and alloys is afforded, yielding a generalization of the classical Gibbs–Thomson relation in the former and novel evolution equations for an individual step in the latter.     

Constrained parameter-splitting perturbation method for the improved solutions to the nonlinear vibrations of Euler–Bernoulli cantilevers

In this paper, a new method named constrained parameter-splitting perturbation method for improving the solutions obtained from the parameter-splitting perturbation method is proposed for solving the problems in some extremal cases, such as the strongly nonlinear vibration of an Euler–Bernoulli cantilever. The proposed method takes the advantages of both the perturbation method and the harmonic balance method. The idea is that the solution obtained by the parameter-splitting perturbation method is substituted into the equation of motion and then the accumulative error of the equation is minimized for determining the unknown splitting parameters under the constraints constructed under the frame of harmonic balance method. The forced vibration of an oscillator with cubic geometric nonlinearity and inertia nonlinearity and the forced vibration of a planar microcantilever beam with a lumped tip mass are studied as examples to reveal the efficacy of the proposed method. The inspection of the steady-state response including its stability is conducted by means of comparing the frequency-response curves obtained by the proposed method with those obtained by the numerical continuation method and harmonic balance method, respectively, to show the efficacy and the advantages of the proposed method. Meanwhile, the nonlinear ordering effect on the solutions of the proposed method is also studied by comparing the results obtained by using different nonlinear orderings in the systems. In the last, we found through convergence examinations that it is necessary to have corrections to the erroneous solution which are obtained by harmonic balance method and Floquet theory in stability analysis.

     

On the harmonic balance with linearization for asymmetric single degree of freedom non-linear oscillators

This paper deals with analytical approximation of non-linear oscillations of conservative asymmetric single degree of freedom systems, using the method of harmonic balance with linearization. This technique which consists of linearizing the governing equations prior to harmonic balance permits us to avoid solving complicated non-linear algebraic equations. But it could be applied only to symmetric oscillations for which it proves to be very simple and effective. This restriction is due to the fact that the method requires an appropriate initial approximate solution as input. Such a solution could not be readily identified for nonsymmetric oscillations, contrary the symmetric case where the fundamental harmonic works well. For these nonsymmetric oscillations, we propose in this paper to consider an initial approximation which consists of a small bias plus the fundamental harmonic. By expanding the corresponding harmonic balance equations respectively to first and second order in the bias, we are able to easily determine the bias and thus the required initial approximate solution that yields consistent solution at higher order. We use three examples to illustrate the proposed approach and reveal its simplicity and its very good convergence.     

Analytical approximation to large-amplitude oscillation of a non-linear conservative system

This paper deals with non-linear oscillation of a conservative system having inertia and static non-linearities. By combining the linearization of the governing equation with the method of harmonic balance, we establish analytical approximate solutions for the non-linear oscillations of the system. Unlike the classical harmonic balance method, linearization is performed prior to proceeding with harmonic balancing, thus resulting in a set of linear algebraic equations instead of one of non-linear algebraic equations. Hence, we are able to establish analytical approximate formulas for the exact frequency and periodic solution. These analytical approximate formulas show excellent agreement with the exact solutions, and are valid for small as well as large amplitudes of oscillation.     

一种单元谐波平衡法

基于有限元离散,对于工程中的非线性响应问题,提出一种单元谐波平衡法．与常规的谐波平衡法不同,本文将谐波平衡方程建立在有限元素上,从而兼顾了有限元素法和常规谐波平衡法两大优势．有限元技术的应用能使得求解问题的范围扩大到复杂工程结构,而谐波平衡概念的使用将使得含有复杂变形和复杂本构关系的动力学响应问题得到有效解决．所提方法能适用于工程结构中具有复杂非线性关系的动力学响应问题．由于谐波平衡法的实施依赖于谐波系数方程及其切线刚度矩阵的解析推导,尽管已经局限到有限元素上,但对于较为复杂一些的本构关系,推导仍非易事．为解决这些问题,放弃通常对于变形梯度和应变张量所作的向量假设,而是从连续介质力学中基本的几何关系入手,提出一种矩阵分解形式．通过利用张量的内蕴导数定义以及关于迹函数的有关性质,给出应力增量的一种新的表现形式．当它与变形梯度的矩阵分解相结合时,使得切线刚度矩阵的导出变得十分简单,而且所得计算形式也比通常紧凑和方便许多．     

Incremental harmonic balance method for periodic forced oscillation of a dielectric elastomer balloon

Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices call for an efficient and reliable method to analyze the dynamic response of DE. This remains to be a challenge since the resultant vibration equation of DE, for example, the vibration of a DE balloon considered here is highly nonlinear with higher-order power terms and time-dependent coefficients. Previous efforts toward this goal use largely the numerical integration method with the simple harmonic balance method as a supplement. The numerical integration and the simple harmonic balance method are inefficient for large parametric analysis or with difficulty in improving the solution accuracy. To overcome the weakness of these two methods,we describe formulations of the incremental harmonic balance(IHB) method for periodic forced solutions of such a unique system. Combined with an arc-length continuation technique, the proposed strategy can capture the whole solution branches, both stable and unstable, automatically with any desired accuracy.     

A 1/3 Pure Subharmonic Solution and Transient Process for the Duffing's Equation

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｍｏｄｅｌｉｎｇｆｏｒｄｙｎａｍｉｃｓｙｓｔｅｍ ,ｗｅｏｆｔｅｎｇｅｔｎｏｎｌｉｎｅａｒｉｔｙｍａｔｈｅｍａｔｉｃａｌｅｑｕａｔｉｏｎ .Ｉｔｉｓｄｉｆｆｉｃｕｌｔｔｏｇｉｖｅｏｕｔｐｕｒｅａｎａｌｙｔｉｃｒｅｓｏｌｕｔｉｏｎｔｏｔｈｅｅｑｕａｔｉｏｎ ,ｍｏｓｔａｎｓｗｅｒｓｗｏｒｋｅｄｏｕｔａｒｅｄｅｐｅｎｄｅｎｔｕｐｏｎｎｕｍｅｒｉｃａｌｓｏｌｕｔｉｏｎ .Ｈｅｎｃｅ,ｉｔｉｓｅｘｔｒｅｍｅｌｙｉｍｐｏｒｔａｎｔｔｏｆｉｎｄｏｕｔａｓｕｉｔａｂｌｅａｐｐｒｏｘｉｍａｔ…