Construction of Invariant Torus Using Toeplitz Jacobian Matrices/Fast Fourier Transform Approach

The invariant torus is a very important case in the study of nonlinear autonomous systems governed by ordinary differential equations (ODEs). In this paper a new numerical method is provided to approximate the multi-periodic surface formed by an invariant torus by embedding the governing ODEs onto a set of partial differential equations (PDEs). A new characteristic approach to determine the stability of resultant periodic surface is also developed. A system with two strongly coupled van der Pol oscillators is taken as an illustrative example. The result shows that the Toeplitz Jacobian Matrix/Fast Fourier Transform (TJM/FFT) approach introduced previously is accurate and efficient in this application. The application of the method to normal multi-modes of nonlinear Euler beam is given in [1].     

Parametric Excitation for Two Internally Resonant van der Pol Oscillators

The response of a system of two nonlinearly coupled van der Poloscillators to a principal parametric excitation in the presence ofone-to-one internal resonance is investigated. The asymptoticperturbation method is applied to derive the slow flow equationsgoverning the modulation of the amplitudes and the phases of the twooscillators. These equations are used to determine steady-stateresponses, corresponding to a periodic motion for the starting system(synchronisation), and parametric excitation-response andfrequency-response curves. Energy considerations are used to studyexistence and characteristics of limit cycles of the slow flowequations. A limit cycle corresponds to a two-period amplitude- andphase-modulated motion for the van der Pol oscillators. Two-periodmodulated motion is also possible for very low values of the parametricexcitation and an approximate analytic solution is constructed for thiscase. If the parametric excitation increases, the oscillation period ofthe modulations becomes infinite and an infinite-period bifurcationsoccur. Analytical results are checked with numerical simulations.     

Constructing Galerkin's approximations of invariant tori using MACSYMA

Invariant tori of solutions for nonlinearly coupled oscillators are generalizations of limit cycles in the phase plane. They are surfaces of aperiodic solutions of the coupled oscillators with the property that once a solution is on the surface it remains on the surface. Invariant tori satisfy a defining system of nonlinear partial differential equations. This case study shows that with the help of a symbolic manipulation package, such as MACSYMA, approximations to the invariant tori can be developed by using Galerkin's variational method. The resulting series must be manipulated efficiently, however, by using the Poisson series representation for multiply periodic functions, which makes maximum use of the list processing techniques of MACSYMA. Three cases are studied for the single van der Pol oscillator with forcing parameter =0.5, 1.0, 1.5, and three cases are studied for a pair of nonlinearly coupled van der Pol oscillators with forcing parameters =0.005, 0.5, 1.0. The approximate tori exhibit good agreement with direct numerical integrations of the systems.Contribution of the National Institute of Standards and Technology, a Federal agency. Not subject to copyright.     

The Response of a Parametrically Excited van der Pol Oscillator to a Time Delay State Feedback

We investigate the parametric resonance of a van der Pol oscillator under state feedback control with a time delay. Using the asymptotic perturbation method, we obtain two slow-flow equations on the amplitude and phase ofthe oscillator. Their fixed points correspond to a periodic motion forthe starting system and we show parametric excitation-response andfrequency-response curves. We analyze the effect of time delay andfeedback gains from the viewpoint of vibration control and use energyconsiderations to study the existence and characteristics of limit cycles of the slow-flow equations. A limit cycle corresponds to a two-periodmodulated motion for the van der Pol oscillator. Analytical results areverified with numerical simulations. In order to exclude the possibilityof quasi-periodic motion and to reduce the amplitude peak of theparametric resonance, we find the appropriate choices for the feedbackgains and the time delay.     

Energy harvesting from quasi-periodic vibrations using electromagnetic coupling with delay

In this paper, we study quasi-periodic vibrational energy harvesting in a delayed self-excited oscillator with a delayed electromagnetic coupling. The energy harvester system consists in a delayed van der Pol oscillator with delay amplitude modulation coupled to a delayed electromagnetic coupling mechanism. It is assumed that time delay is inherently present in the mechanical subsystem of the harvester, while it is introduced in the electrical circuit to control and optimize the output power of the system. A double-step perturbation method is performed near a delay parametric resonance to approximate the quasi-periodic solutions of the harvester which are used to extract the quasi-periodic vibration-based power. The influence of the time delay introduced in the electromagnetic subsystem on the performance of the quasi-periodic vibration-based energy harvesting is examined. In particular, it is shown that for appropriate values of amplitudes and frequency of time delay the maximum output power of the harvester is not necessarily accompanied by the maximum amplitude of system response.     

Vibration control for the primary resonance of the van der Pol oscillator by a time delay state feedback

We investigate the primary resonance of an externally excited van der Pol oscillator under state feedback control with a time delay. By means of the asymptotic perturbation method, two slow-flow equations on the amplitude and phase of the oscillator are obtained and external excitation-response and frequency-response curves are shown. We discuss how vibration control and high amplitude response suppression can be performed with appropriate time delay and feedback gains. Moreover, energy considerations are used in order to investigate existence and characteristics of limit cycles of the slow-flow equations. A limit cycle corresponds to a two-period modulated motion for the van der Pol oscillator. We demonstrate that appropriate choices for the feedback gains and the time delay can exclude the possibility of modulated motion and reduce the amplitude peak of the primary resonance. Analytical results are verified with numerical simulations.     

多尺度法二义性的一种解释

多尺度法是为解决含小参数系统发展起来的应用最广泛的摄动法之一. 在求解高阶近 似方程时,多尺度法一般只求特解. 用多尺度法求解van der Pol 方程的三阶解时 将出现矛盾. 以van der Pol方程为例,证明了忽略一阶修正量中的一阶谐波 项使得混合偏导数不能交换顺序,从而导致了多尺度法的二义性和另一个数学矛盾. 在求解一阶修正量时采用含有一阶谐波项的全解,消除了二义性和该矛盾. 该 方法所求得的近似解与数值解进行了比较,结果非常吻合,验证了其合理性.     

Frequency lock-in during nonlinear vibration of an airfoil coupled with van der Pol Oscillator

The frequency lock-in during the nonlinear vibration of a turbomachinery blade is modeled using a spring-mounted airfoil coupled with a van der Pol Oscillator (VDP) oscillator. The proposed reduced-order model uses the nonlinear VDP oscillator to represent the oscillatory nature of wake dynamics caused by the vortex shedding. The damping term in the VDP oscillator is assumed to be nonlinear. The coupled equations governing the pitch and plunge motion of an airfoil are used to approximate the vibration of a turbomachinery blade. Springs having cubic-order nonlinearity for their stiffnesses are used to mount the airfoil. The unsteady lift acting on the blade is modeled using a self-excited nonlinear wake oscillator. The model for wake dynamics takes into account the influence of blade inertia. The nonlinear coupled three degrees of freedom (dof) aeroelastic system is studied for instability resulting in the frequency lock-in phenomenon. The equations are transformed into non-dimensional form, and then the frequencies of the coupled system are plotted to demonstrate the frequency lock-in. Further, the method of multiple scales is used to derive modulation equations which represent the amplitude and phase of the oscillation. The results obtained using the method of multiple scales are compared with direct numerical solutions to verify the present modeling method. The steady-state amplitudes of the response are plotted against the detuning parameter, which represents the frequency response curve. Further, the sensitivity of non-dimensional parameters such as coupling coefficients, mass ratio, reduced velocity, static unbalance, structural damping coefficient and the ratio of uncoupled pitch and plunge natural frequencies on the frequency response is investigated. The study revealed that parameters such as mass ratio, reduced velocity, structural damping coefficient, and coupling coefficients have a stronger influence in suppressing the amplitude of vibration. Meanwhile, parameters such as the frequency ratio, static unbalance, reduced velocity, and mass ratio significantly affect the range of frequency in which the lock-in phenomenon happens. Further, linear perturbation analysis is done to understand the qualitative effect of the system parameters such as coupling coefficients, mass ratio, frequency ratio, and static unbalance on the range of lock-in.     

Bifurcation and dynamic response analysis of rotating blade excited by upstream vortices

A reduced model is proposed and analyzed for the simulation of vortex-induced vibrations (VIVs) for turbine blades. A rotating blade is modelled as a uniform cantilever beam, while a van der Pol oscillator is used to represent the time-varying characteristics of the vortex shedding, which interacts with the equations of motion for the blade to simulate the fluid-structure interaction. The action for the structural motion on the fluid is considered as a linear inertia coupling. The nonlinear characteristics for the dynamic responses are investigated with the multiple scale method, and the modulation equations are derived. The transition set consisting of the bifurcation set and the hysteresis set is constructed by the singularity theory and the effects of the system parameters, such as the van der Pol damping. The coupling parameter on the equilibrium solutions is analyzed. The frequency-response curves are obtained, and the stabilities are determined by the Routh-Hurwitz criterion. The phenomena including the saddle-node and Hopf bifurcations are found to occur under certain parameter values. A direct numerical method is used to analyze the dynamic characteristics for the original system and verify the validity of the multiple scale method. The results indicate that the new coupled model is useful in explaining the rich dynamic response characteristics such as possible bifurcation phenomena in the VIVs.     

Regular oscillations,chaos, and multistability in a system of two coupled van der Pol oscillators: numerical and experimental studies

In this paper, the dynamics of a system of two coupled van der Pol oscillators is investigated. The coupling between the two oscillators consists of adding to each one’s amplitude a perturbation proportional to the other one. The coupling between two laser oscillators and the coupling between two vacuum tube oscillators are examples of physical/experimental systems related to the model considered in this paper. The stability of fixed points and the symmetries of the model equations are discussed. The bifurcations structures of the system are analyzed with particular attention on the effects of frequency detuning between the two oscillators. It is found that the system exhibits a variety of bifurcations including symmetry breaking, period doubling, and crises when monitoring the frequency detuning parameter in tiny steps. The multistability property of the system for special sets of its parameters is also analyzed. An experimental study of the coupled system is carried out in this work. An appropriate electronic simulator is proposed for the investigations of the dynamic behavior of the system. Correspondences are established between the coefficients of the system model and the components of the electronic circuit. A comparison of experimental and numerical results yields a very good agreement.     

Multiple scales analysis for double Hopf bifurcation with?1:3?resonance

The dynamical behavior of a general n-dimensional delay differential equation (DDE) around a 1:3 resonant double Hopf bifurcation point is analyzed. The method of multiple scales is used to obtain complex bifurcation equations. By expressing complex amplitudes in a mixed polar-Cartesian representation, the complex bifurcation equations are again obtained in real form. As an illustration, a system of two coupled van der Pol oscillators is considered and a set of parameter values for which a 1:3 resonant double Hopf bifurcation occurs is established. The dynamical behavior around the resonant double Hopf bifurcation point is analyzed in terms of three control parameters. The validity of analytical results is shown by their consistency with numerical simulations.     

A refined asymptotic perturbation method for nonlinear dynamical systems

In this paper, a refined asymptotic perturbation method for general nonlinear dynamical systems is proposed for the first time. This method can be considered as an alternative means for the traditional multiple scales method. Moreover, it is easier to be understood and used to carry out higher-order perturbation analysis. In addition, three examples including the Duffing equation, a system with quadratic and cubic nonlinearities to a subharmonic excitation, as well as the coupled van der Pol oscillator with parametrical excitations are investigated to illustrate the validity and usefulness of the proposed technique. The analytical and numerical results show good agreement.     

The transition from phase locking to drift in a system of two weakly coupled van der pol oscillators

We investigate the slow flow resulting from the application of the two variable expansion perturbation method to a system of two linearly coupled van der Pol oscillators. The slow flow consists of three non-linear coupled odes on the amplitudes and phase difference of the oscillators. We obtain regions in parameter space which correspond to phase locking, phase entrainment and phase drift of the coupled oscillators. In the slow flow, these states correspond respectively to a stable equilibrium, a stable limit cycle and a stable libration orbit. Phase entrainment, in which the phase difference between the oscillators varies periodically, is seen as an intermediate state between phase locking and phase drift. In the slow flow, the transitions between these states are shown to be associated with Hopf and saddle-connection bifurcations.     

Response statistics of van der Pol oscillators excited by white noise

The joint probability density function of the state space vector of a white noise exoited van der Pol oscillator satisfics a Fokker-Planck-Kolmogorov (FPK) equation. The paper describes a numerical procedure for solving the transient FPK equation based on the path integral solution (PIS) technique. It is shown that by combining the PIS with a cubic B-spline interpolation method, numerical solution algorithms can be implemented giving solutions of the FPK equation that can be made accurate down to very low probability levels. The method is illustrated by application to two specific examples of a van der Pol oscillator.     

Chaotic flexural oscillations of a spinning nanoresonator

The paper investigates the chaotic flexural oscillations of the spinning nanoresonator. The influence of cubic nonlinearity arising from the van der Waals interactions between two neighboring layers of carbon nanotubes on the structural oscillations of the system is considered. The integral–differential equations describing the flexural displacements of the nanoresonator are discretized into two coupled Duffing-type equations using the Galerkin–Ritz procedures. The linear stiffness can be either positive or negative, depending on the amplitudes of the linear trap rigidity arising from both the van der Waals interactions and the axial tensile loads. The chaotic flexural oscillations of the appropriately excited spinning nanoresonator are predicted theoretically. Using the Nayfeh–Mook multiscale perturbation algorithms, the coupled Duffing-type equations with linear positive stiffness may be transformed into autonomous equations of slowly modulated amplitudes whose equilibrium points and chaotic dynamics are investigated numerically. The potential chaotic oscillations of the elastic nanoresonator can be determined by the Melnikov–Holmes–Marsden (MHM) integral associated with the homoclinic/heteroclinic solutions of the disturbed Hamiltonian systems with linear negative stiffness. The findings are validated through the Poincare sections and Lyapunov exponents.     

Study of mixed-mode oscillations in a parametrically excited van der Pol system

This paper investigates the effects of slowly varying parametric excitation on the dynamics of van der Pol system. Periodic bifurcation delay behaviors are exhibited when the parametric excitation slowly passes through Hopf bifurcation value of the controlled van der Pol system. The first bifurcation delay behavior relies on initial conditions, while the bifurcation delay behaviors that follow the first one are immune to initial conditions. These bifurcation delay behaviors result in a hysteresis loop between the spiking attractor and the rest state, which is responsible for the generation of mixed-mode oscillations. Then an approximate calculation for the number of spikes in each cluster of repetitive spiking of mixed-mode oscillations is explored based on bifurcation delay behaviors. Theoretical results agree well with numerical simulations.     

Using delay to quench undesirable vibrations

A van der Pol type system with delayed feedback is explored by employing the two variable expansion perturbation method. The perturbation scheme is based on choosing a critical value for the delay corresponding to a Hopf bifurcation in the unperturbed ε=0 system. The resulting amplitude–delay relation predicts two Hopf bifurcation curves, such that in the region between these two curves oscillations will be quenched. The perturbation results are verified by comparison with numerical integration.     

Stability and bifurcation analysis in the delay-coupled nonlinear oscillators

This paper investigates the dynamical behavior of two oscillators with nonlinearity terms, which are coupled with finite delay parameters. Each oscillator is a general class of second-order nonlinear delay-differential equations. The system of delay differential equations is analyzed by reducing the delay equations to a system of ordinary differential equations on a finite-dimensional center manifold, the corresponding to an infinite-dimensional phase space. In addition, the characteristic equation for the linear stability of the trivial equilibrium is completely analyzed and the stability region is illustrated in the parameters space. Our analysis reveals necessary coefficients of the reduced vector field on the center manifold for studying the bifurcations of the trivial equilibrium such as transcritical, pitchfork, and Hopf bifurcation. Finally, we consider the delay-coupled van der Pol equations.     

Dynamical analysis of Mathieu equation with two kinds of van der Pol fractional-order terms

In this paper the dynamics of Mathieu equation with two kinds of van der Pol (VDP) fractional-order terms is investigated. The approximately analytical solution is obtained by the averaging method. The steady-state solution, existence conditions and stability condition for the steady-state solution are presented, and it is found that the two kinds of VDP fractional coefficients and fractional orders remarkably affect the steady-state solution, which is characterized by the additional damping coefficient (ADC) and additional stiffness coefficient (ASC). The comparisons between the analytical and numerical solutions verify the correctness and satisfactory precision of the approximately analytical solution. The presented typical amplitude–frequency curves illustrate the important effects of two kinds of VDP fractional-order terms on system dynamics. The application of two VDP fractional-order terms in vibration control is discussed. At last, the detailed results are summarized and the conclusions are made.     

Multi-frequency oscillations in self-excited systems

A self-excited three-mass chain system is considered here. For a self-excitation of van der Pol type, the possibility of multi-frequency oscillations is investigated. Both analytical approximate solutions and numerical simulation are used. The averaging method is used to establish existence and stability of the normal modes, the two-frequency modes as well as the three-frequency oscillations solutions. We found at first that the single mode seems to prevail. However a three-frequency solution can be stabilised by adapting the system slightly. A generic bifurcation diagram is given where all the possible phase portraits are sketched. The flow turns out to be quite predictable. There is no “room” for chaos or strange attractors. This behaviour is not typical for systems of coupled oscillators but turns out to be partly related to the involved symmetries as well as the particular choice of the system parameters.