Global residue harmonic balance method to periodic solutions of a class of strongly nonlinear oscillators

A new approach, namely the global residue harmonic balance method, was advanced to determine the accurate analytical approximate periodic solution of a class of strongly nonlinear oscillators. A class of nonlinear jerk equation containing velocity-cubed and velocity times displacements-squared was taken as a typical example. Unlike other harmonic balance methods, all the former residual errors are introduced in the present approximation to improve the accuracy. Comparison of the result obtained using this approach with the exact one and simplicity and efficiency of the proposed procedure. The method can be easily extended to other strongly nonlinear oscillators.     

Efficient determination of higher-order periodic solutions using n-mode harmonic balance

This paper presents a systematic procedure to explicitly determinethe algebraic equations arising from the method of harmonicbalance with an arbitrary number of modes in the assumed solutions.The technique can be used for a wide variety of nonlinear oscillators(including systems of ordinary differential equations). Themethod is illustrated in the case of second-order differentialequations with nonlinear restoring force. Although numericalmethods have been employed to solve the resulting systems ofalgebraic equations, the general approach is analytic. As such,this study confirms independently (i.e. nonsimulation) the period-doublingcascade of an escape equation including the bifurcation universalscaling laws.     

Relationship between the homotopy analysis method and harmonic balance method

This paper presents a study of the relationship between the homotopy analysis method (HAM) and harmonic balance (HB) method. The HAM is employed to obtain periodic solutions of conservative oscillators and limit cycles of self-excited systems, respectively. Different from the usual procedures in the existing literature, the HAM is modified by retaining a given number of harmonics in higher-order approximations. It is proved that as long as the solution given by the modified HAM is convergent, it converges to one HB solution. The Duffing equation, the van der Pol equation and the flutter equation of a two-dimensional airfoil are taken as illustrations to validate the attained results.     

Nonlinear transversely vibrating beams by the improved energy balance method and the global residue harmonic balance method

Mathematical modeling of many engineering systems such as beam structures often leads to nonlinear ordinary or partial differential equations. Nonlinear vibration analysis of the beam structures is very important in mechanical and industrial applications. This paper presents the high order frequency-amplitude relationship for nonlinear transversely vibrating beams with odd and even nonlinearities using the improved energy balance method and the global residue harmonic balance method. The accuracy of the energy balance method is improved based on combining features of collocation method and Galerkin–Petrov method, and an improved harmonic balance method is proposed which is called the global residue harmonic balance method. Unlike other harmonic balance methods, all the former global residual errors are introduced in the present approximation to improve the accuracy. Finally, preciseness of the present analytic procedures is evaluated in contrast with numerical calculations methods, giving excellent results.     

Elliptic harmonic balance method for two degree-of-freedom self-excited oscillators

Elliptic harmonic balance (EHB) method as an analytical method is widely used for strongly non-linear oscillators with a single degree-of-freedom (DOF). The oscillation equations are expressed by elliptic functions, and then the expressions are expanded as harmonics of elliptic type while only the first harmonic is retained. To the best of our knowledge, however, it seems that the EHB method has not been found applications in two or multiple DOFs systems. One possible reason is that the EHB method may cause a difficult problem that the number of equations obtained by harmonic balancing is not equal to that of unknowns. To this end, in the present paper, the EHB method is therefore extended to study a class of strongly self-excited oscillators with two DOFs. Prior to harmonic balancing, an additional equation is established to tackle the aforementioned problem. Illustrative examples show that solutions of limit cycles obtained by the proposed method are in good agreement with the numerical solutions.     

Application of the method of Lyapunov periodic functions

The paper is concerned with asymptotic behavior of continuous and discrete phase control systems involving periodic differentiable nonlinear vector functions and featuring nonunique equilibrium state.     

Localization of hidden attractors of the generalized Chua system based on the method of harmonic balance

Chua’s circuits, which were introduced by Leon Ong Chua in 1983, are simplest electric circuits operating in the mode of chaotic oscillations. Systems of differential equations describing the behavior of Chua’s circuits are three-dimensional autonomous dynamical systems with scalar nonlinearity. In the standard Chua system, chaotic oscillations are excited in the classical manner, namely, starting from the neighborhood of the unstable zero equilibrium, after the transient process, the system trajectory tends to a Chua attractor.     

Application of the discontinuous Galerkin finite element method in small deformation regimes

In this paper, a finite element formulation is defined in the framework of the discontinuous Galerkin method. Discontinuous Galerkin (dG) methods are classically used in fluid mechanics, however recently their application in solid mechanics has become more vivid among scientists. Of special interest is their application in elliptic problems with constraints such as incompressibility which leads to volumetric locking phenomenon and also in some structural models of shells, plates and beams with compatibility constraints, which brings about shear locking [1]. While classical standard Galerkin methods must be continuous, dG methods can be applied for discontinuities across element boundaries, where a jump of a value (displacement) can be observed. In the present work, a dG method is applied to a linear elastic bar, where a weak discontinuity is allowed in the bar. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Incremental harmonic balance method for nonlinear flutter of an airfoil with uncertain-but-bounded parameters

The limit cycle oscillation of a two-dimensional airfoil with parameter variability in an incompressible flow is investigated using the incremental harmonic balance (IHB) method. The variable parameters, such as the wind speed, the cubic plunge and pitch stiffness coefficients, are modeled as either bounded uncertain or stochastic parameters. In the solution process of the IHB method, the bounded parameters are considered as an active increment. Taking all values over the considered bounded regions of the active parameters provides us with a series of IHB solutions of limit cycle oscillations of the airfoil. With the aid of the attained solutions, the bounds and some statistical properties of the limit cycle oscillations are determined and compared with Monte Carlo simulation (MCS) results. Numerical examples show that the proposed approach is valid and effective for analyzing strongly nonlinear vibration problems with bounded uncertainties.     

An easy trick to a periodic solution of relativistic harmonic oscillator

In this paper, the relativistic harmonic oscillator equation which is a nonlinear ordinary differential equation is investigated by Homotopy perturbation method. Selection of a linear operator, which is a part of the main operator, is one of the main steps in HPM. If the aim is to obtain a periodic solution, this choice does not work here. To overcome this lack, a linear operator is imposed, and Fourier series of sines will be used in solving the linear equations arise in the HPM. Comparison of the results, with those of resulted by Differential Transformation and Harmonic Balance Method, shows an excellent agreement.     

Stable and unstable flow regimes for active fluids in the periodic setting

Depending on the involved physiobiological parameters, stable or unstable behavior in active fluids is observed. In this paper a rigorous analytical justification of (in-)stability within the corresponding regimes is given. In particular, occurring instability for the manifold of ordered polar states caused by self-propulsion is proved. This represents the prerequisite for active turbulence patterns as observed in a number of applications. The approach is carried out in the periodic setting and is based on the generalized principle of linearized (in)-stability related to normally stable and normally hyperbolic equilibria.     

Application of the invariant embedding method for computing nonlinear periodic oscillations

The invariant embedding method is applied to compute periodic solutions of nonlinear systems that depend continuously on the perturbation frequency. Analysis of such systems in the region of multivalued solutions is demonstrated. An example of an essentially nonlinear system is examined.Khar'kov Polytechnical Institute. Translated from Vychislitel'naya i Prikladnaya Matematika, No. 68, pp. 136–140, 1989.     

A periodic functional approach to the calculus of variations and the problem of time-dependent damped harmonic oscillators

I discuss the problem of time-dependent harmonic oscillators on the basis of a periodic functional approach to the calculus of variations. Both the Lagrangian and Hamiltonian formulations are explored and discussed in some detail. Some interesting consequences are revealed.     

Exact front, soliton and hole solutions for a modified complex Ginzburg-Landau equation from the harmonic balance method

A novel approach of using harmonic balance (HB) method is presented to find front, soliton and hole solutions of a modified complex Ginzburg-Landau equation. Three families of exact solutions are obtained, one of which contains two parameters while the others one parameter. The HB method is an efficient technique in finding limit cycles of dynamical systems. In this paper, the method is extended to obtain homoclinic/heteroclinic orbits and then coherent structures. It provides a systematic approach as various methods may be needed to obtain these families of solutions. As limit cycles with arbitrary value of bifurcation parameter can be found through parametric continuation, this approach can be extended further to find analytic solution of complex quintic Ginzburg-Landau equation in terms of Fourier series.     

Nonlinear analysis of a parametrically excited beam with intermediate support by using Multi-dimensional incremental harmonic balance method

In this paper, a nonlinear Euler-Bernoulli beam under a concentrated harmonic excitation with intermediate nonlinear support is investigated. Continuous expression for the kinetic energy, potential energy and dissipation function are constructed. An energy method based on the Lagrange equation combined with the Galerkin truncation is used for discretizing the governing equation. The Multi-dimensional incremental harmonic balance method (MIHBM) is derived, and the comparisons between the numerical results and the approximate analytical solutions based on the MIHBM verify the excellent accuracy of the MIHBM. The steady state dynamic of the beam is investigated by MIHBM. In order to investigate the energy transmission and understand the vibration response of the Euler-Bernoulli beam, the effects of the key parameters on the dynamic behaviors are studied and discussed, individually. The results show that the amplitude-frequency curves exhibits softening nonlinear behavior in the super-harmonic resonance region, and near resonant region the hardening nonlinear behavior is observed depending on the different parameters. Nonlinear dynamic analysis, such as bifurcation, 3-D frequency spectrum, waveform, frequency spectrum, phase diagram and Poincaré map, are also presented in order to study the influences of the key parameters on the vibration behaviors for the beam in a more accurate manner. In addition, the path to chaotic motion is observed to be through a sequence of the periodic motion and quasi-periodic motion.     

Parameters of the periodic response of a system with friction-affected constraints and harmonic excitation

Summary The periodic response of a single-freedom multi-body system with friction-affected constraints acted upon by a harmonic excitation is determined by numerical simulation. The contribution of the constraints to the generalized friction force, the influence of gravity and zero-gravity environment, and the consequences, when several constraints are considered without friction, are investigated. The effects of a small artificial damping on the response, and the possibilities of an equivalent damping for the friction-affected constraints are also examined.

---

Zusammenfassung Die periodische Bewegung eines Mehrkörpersystems mit einem Freiheitsgrad und reibungsbehafteten Bindungen bei harmonischer Erregung wird durch numerische Simulation bestimmt. Untersucht werden der Beitrag der Bindungen zur verallgemeinerten Reibungskraft, der Einfluß des Schwerefeldes und der Schwerelosigkeit und die Folgen, wenn einige Bindungen als reibungsfrei angenommen werden. Die Auswirkungen einer kleinen künstlichen Dämpfung auf die Bewegung, und die Möglichkeiten einer äquivalenten Dämpfung an Stelle der reibungsbehafteten Bindungen wird auch untersucht.

     

Applying carvalho's method to find periodic solutions of difference equations

In this paper we apply the method initially developed in [1] for differential-difference equations, to the case of difference equations, in order to find 2 and 3-periodic solutions of some equations that often appear in the literatures as are for instance the case of Applications 2,5 which are examples of population growth models, and Application 4, which is a standard example of nonlinear higher order scalar difference equation depending on two parameters (see, Kocik and Ladas [3]).