Stability of two-mode internal resonance in a nonlinear oscillator

We analyze the stability of synchronized periodic motion for two coupled oscillators, representing two interacting oscillation modes in a nonlinear vibrating beam. The main oscillation mode is governed by the forced Duffing equation, while the other mode is linear. By means of the multiple-scale approach, the system is studied in two situations: an open-loop configuration, where the excitation is an external force, and a closed-loop configuration, where the system is fed back with an excitation obtained from the oscillation itself. The latter is relevant to the functioning of time-keeping micromechanical devices. While the accessible amplitudes and frequencies of stationary oscillations are identical in the two situations, their stability properties are substantially different. Emphasis is put on resonant oscillations, where energy transfer between the two coupled modes is maximized and, consequently, the strong interdependence between frequency and amplitude caused by nonlinearity is largely suppressed.     

On the response of a harmonically excited two degree-of-freedom system consisting of a linear and a nonlinear quasi-zero stiffness oscillator

There are many systems which consist of a nonlinear oscillator attached to a linear system, examples of which are nonlinear vibration absorbers, or nonlinear systems under test using shakers excited harmonically with a constant force. This paper presents a study of the dynamic behaviour of a specific two degree-of-freedom system representing such a system, in which the nonlinear system does not affect the vibration of the forced linear system. The nonlinearity of the attachment is derived from a geometric configuration consisting of a mass suspended on two springs which are adjusted to achieve a quasi-zero stiffness characteristic with pure cubic nonlinearity. The response of the system at the frequency of excitation is found analytically by applying the method of averaging. The effects of the system parameters on the frequency-amplitude response of the relative motion are examined. It is found that closed detached resonance curves lying outside or inside the continuous path of the main resonance curve can appear as a part of the overall amplitude-frequency response. Two typical situations for the creation of the detached resonance curve inside the main resonance curve, which are dependent on the damping in the nonlinear oscillator, are discussed.     

Suppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber

In a single degree-of-freedom weakly nonlinear oscillator subjected to periodic external excitation, a small-amplitude excitation may produce a relatively large-amplitude response under primary resonance conditions. Jump and hysteresis phenomena that result from saddle-node bifurcations may occur in the steady-state response of the forced nonlinear oscillator. A simple mass-spring-damper vibration absorber is thus employed to suppress the nonlinear vibrations of the forced nonlinear oscillator for the primary resonance conditions. The values of the spring stiffness and mass of the vibration absorber are significantly lower than their counterpart of the forced nonlinear oscillator. Vibrational energy of the forced nonlinear oscillator is transferred to the attached light mass through linked spring and damper. As a result, the nonlinear vibrations of the forced oscillator are greatly reduced and the vibrations of the absorber are significant. The method of multiple scales is used to obtain the averaged equations that determine the amplitude and phases of the first-order approximate solutions to primary resonance vibrations of the forced nonlinear oscillator. Illustrative examples are given to show the effectiveness of the dynamic vibration absorber for suppressing primary resonance vibrations. The effects of the linked spring and damper and the attached mass on the reduction of nonlinear vibrations are studied with the help of frequency response curves, the attenuation ratio of response amplitude and the desensitisation ratio of the critical amplitude of excitation.     

An adaptive harmonic balance method for predicting the nonlinear dynamic responses of mechanical systems—Application to bolted structures

Aeronautical structures are commonly assembled with bolted joints in which friction phenomena, in combination with slapping in the joint, provide damping on the dynamic behavior. Some models, mostly nonlinear, have consequently been developed and the harmonic balance method (HBM) is adapted to compute nonlinear response functions in the frequency domain. The basic idea is to develop the response as Fourier series and to solve equations linking Fourier coefficients. One specific HBM feature is that response accuracy improves as the number of harmonics increases, at the expense of larger computational time. Thus this paper presents an original adaptive HBM which adjusts the number of retained harmonics for a given precision and for each frequency value. The new proposed algorithm is based on the observation of the relative variation of an approximate strain energy for two consecutive numbers of harmonics. The developed criterion takes the advantage of being calculated from Fourier coefficients avoiding time integration and is also expressed in a condensation case. However, the convergence of the strain energy has to be smooth on tested harmonics and this constitutes a limitation of the method. Condensation and continuation methods are used to accelerate calculation. An application case is selected to illustrate the efficiency of the method and is composed of an asymmetrical two cantilever beam system linked by a bolted joint represented by a nonlinear LuGre model. The practice of adaptive HBM shows that, for a given value of the criterion, the number of harmonics increases on resonances indicating that nonlinear effects are predominant. For each frequency value, convergence of approximate strain energy is observed. Emergence of third and fifth harmonics is noticed near resonances both on vibratory responses and on approximate strain energy. Parametric studies are carried out by varying the excitation force amplitude and the threshold value of the adaptive algorithm. Maximal amplitudes of vibration and frequency response functions are plotted for three different points of the structure. Nonlinear effects become more predominant for higher force amplitudes and consequently the number of retained harmonics is increased.     

Intraband discrete breathers in disordered nonlinear systems. I. Delocalization

The existence of Discrete Breathers or DBs (also called Intrinsic Localized Modes or ILMs) and multibreathers, is investigated in a simple one-dimensional chain of random anharmonic oscillators with quartic potentials coupled by springs. When the breather frequency is outside and above the linearized (phonon) spectrum, the existence theorems and numerical methods previously used in periodic nonlinear models for finding time-periodic and spatially localized solutions, hold identically in random nonlinear systems. These solutions are extraband discrete breathers (EDBs). When the frequencies penetrate inside the linearized spectrum, the existence theorems do not hold. Our numerical investigations demonstrate that the strict continuation of (localized) EDBs as intraband discrete breathers (IDBs) is impossible because of cascades of bifurcations generating many discontinuities. A detailed analysis of these bifurcations for small systems with increasing sizes, shows that only a relatively small subset of the spatially extended multibreathers can be strictly continued while their frequency varies inside the phonon spectrum. We propose an ansatz for finding the coding sequences of these solutions and continuing safely these multibreathers in finite systems of any size. This continuation ends at a lower limit frequency where the solution annihilates through a bifurcation with another multibreather. A smaller subset of these multibreather solutions can be continued to amplitude zero and become linear localized modes at this limit. Conversely, any linear localized mode can be continued when increasing its frequency as an extended multibreather. Extrapolation of these results to infinite systems yields the main conclusion of this first part which is that nonlinearity in disordered systems (with localized eigenmodes only) restores their capability of energy transportation by generating infinitely many spatially extended time-periodic solutions. This approach yields mainly spatially extended solutions, except sometimes at their bifurcation points. In the second part of this work, which is presented in our next article, we develop an accurate method for calculating in situ localized IDBs.     

Multi-mode nonlinear resonance ultrasound spectroscopy for defect imaging: an analytical approach for the one-dimensional case

A nonlinear version of the resonance ultrasound spectroscopy (RUS) theory is presented as an extension of the RUS formalism to the treatment of microdamage characterized by nonlinear constitutive equations. General analytical equations are derived for the one-dimensional case, describing the excitation amplitude dependent shift in the resonance frequency and the generation of harmonics resulting from the interaction between bar modes due to the presence of either localized or volumetrically distributed nonlinearity. Solutions are obtained for classical cubic nonlinearity, as well as for the more interesting case of hysteresis nonlinearity. The analytical results are in excellent quantitative agreement with numerical calculations from a multiscale model. Finally, the analytical formulas are exploited to infer critical information about damage position, degree of nonlinearity, and width of the damage zone either from the shifts in resonance frequency occurring at different excitation modes, or from the shift and the harmonics predicted at a single mode. Unlike other techniques, the multi-mode-nonlinear RUS method does not require a spatial scan to locate the defect, as it lets different excitation modes, with different vibration patterns, probe the structure. Two general methods are suggested for inverting experimental data.     

The effect of nonlinearity on unstable zones of Mathieu equation

Mathieu equation is a well-known ordinary differential equation in which the excitation term appears as the non-constant coefficient. The mathematical modelling of many dynamic systems leads to Mathieu equation. The determination of the locus of unstable zone is important for the control of dynamic systems. In this paper, the stable and unstable regions of Mathieu equation are determined for three cases of linear and nonlinear equations using the homotopy perturbation method. The effect of nonlinearity is examined in the unstable zone. The results show that the transition curves of linear Mathieu equation depend on the frequency of the excitation term. However, for nonlinear equations, the curves depend also on initial conditions. In addition, increasing the amplitude of response leads to an increase in the unstable zone.     

The fractional discrete nonlinear Schrödinger equation

We examine a fractional version of the discrete nonlinear Schrödinger (dnls) equation, where the usual discrete laplacian is replaced by a fractional discrete laplacian. This leads to the replacement of the usual nearest-neighbor interaction to a long-range intersite coupling that decreases asymptotically as a power-law. For the linear case, we compute both, the spectrum of plane waves and the mean square displacement of an initially localized excitation in closed form, in terms of regularized hypergeometric functions, as a function of the fractional exponent. In the nonlinear case, we compute numerically the low-lying nonlinear modes of the system and their stability, as a function of the fractional exponent of the discrete laplacian. The selftrapping transition threshold of an initially localized excitation shifts to lower values as the exponent is decreased and, for a fixed exponent and zero nonlinearity, the trapped fraction remains greater than zero.     

Nonlinear model order reduction of jointed structures for dynamic analysis

Assembled structures generally show weak nonlinearity, thus it is rather commonplace to assume that their modes are both linear and uncoupled. At small to modest amplitude, the linearity assumption remains correct in terms of stiffness but, on the contrary, the dissipation in joints is strongly amplitude-dependent. Besides, the modes of any large structure may be LOCALLY collinear in the localized region of a joint. As a result the projection of the structure on normal modes is not appropriate since the corresponding generalized coordinates may be strongly coupled. Instead of using this global basis, the present paper deals with the use of a local basis to reduce the size of the problem without losing the nonlinear physics. Under an appropriate set of assumptions, the method keeps the dynamic properties of joints, even for large amplitude, which include coupling effects, nonlinear damping and softening effects. The formulation enables us to take into account FE models of any realistic geometry. It also gives a straightforward process for experimental identification. The formulation is detailed and investigated on a jointed structure.     

Multiple internal resonances and non-planar dynamics of shallow suspended cables to the harmonic excitations

In the present study, the nonlinear response of a shallow suspended cable with multiple internal resonances to the primary resonance excitation is investigated. The method of multiple scales is applied directly to the nonlinear equations of motion and associated boundary conditions to obtain the modulation equations and approximate solutions of the cable. Frequency–response curves and force–response curves are used to study the equilibrium solution and its stability. The effects of the excitation amplitude on the frequency–response curves of the cable are also analyzed. Moreover, the chaotic dynamics of the shallow suspended cable is investigated by means of numerical simulations.     

Localization of nonlinear waves in elastic bodies with inclusions

By the example of the problem of the motion of a semi-infinite string lying on an elastic base, a method for describing wave localization near inclusions is proposed for the case of a cubic nonlinearity of the base. The method applies the perturbation technique to the amplitude of a localized mode. The nature of the divergences is revealed, and the secular terms are found to belong to one of two types: inphase or antiphase with the localized wave. It is shown that a combination of the renormalization method and multiscale method provides a convergence of the solutions, which are sought for in the form of power series in the amplitude of the localized mode. It is found that the localization process is determined by the type of the discrete spectrum, type of the nonlinearity, and type of dispersion. The nonlinearity of the elastic base produces two characteristic effects. First, the frequency of the localized wave becomes dependent on the wave amplitude. Second, the system can generate traveling waves at multiple frequencies, which withdraw energy from the localized wave and cause it to decay. The decay behavior is determined by the minimum frequency of these traveling waves (because it must be higher than the cutoff frequency). The lifetime of the localized wave as a function of the mass of a dynamic inclusion exhibits a number of maxima. In particular, the first maximum corresponds to the minimum amplitude of the traveling wave at the triple frequency.     

Theory and experiments for large-amplitude vibrations of empty and fluid-filled circular cylindrical shells with imperfections

The large-amplitude response of perfect and imperfect, simply supported circular cylindrical shells to harmonic excitation in the spectral neighbourhood of some of the lowest natural frequencies is investigated. Donnell's non-linear shallow-shell theory is used and the solution is obtained by the Galerkin method. Several expansions involving 16 or more natural modes of the shell are used. The boundary conditions on the radial displacement and the continuity of circumferential displacement are exactly satisfied. The effect of internal quiescent, incompressible and inviscid fluid is investigated. The non-linear equations of motion are studied by using a code based on the arclength continuation method. A series of accurate experiments on forced vibrations of an empty and water-filled stainless-steel shell have been performed. Several modes have been intensively investigated for different vibration amplitudes. A closed loop control of the force excitation has been used. The actual geometry of the test shell has been measured and the geometric imperfections have been introduced in the theoretical model. Several interesting non-linear phenomena have been experimentally observed and numerically reproduced, such as softening-type non-linearity, different types of travelling wave response in the proximity of resonances, interaction among modes with different numbers of circumferential waves and amplitude-modulated response. For all the modes investigated, the theoretical and experimental results are in strong agreement.     

Contribution to harmonic balance calculations of self-sustained periodic oscillations with focus on single-reed instruments

The harmonic balance method (HBM) was originally developed for finding periodic solutions of electronical and mechanical systems under a periodic force, but has been adapted to self-sustained musical instruments. Unlike time-domain methods, this frequency-domain method does not capture transients and so is not adapted for sound synthesis. However, its independence of time makes it very useful for studying any periodic solution, whether stable or unstable, without care of particular initial conditions in time. A computer program for solving general problems involving nonlinearly coupled exciter and resonator, HARMBAL, has been developed based on the HBM. The method as well as convergence improvements and continuation facilities are thoroughly presented and discussed in the present paper. Applications of the method are demonstrated, especially on problems with severe difficulties of convergence: the Helmholtz motion (square signals) of single-reed instruments when no losses are taken into account, the reed being modeled as a simple spring.     

Nonlinear dynamics of harmonically excited circular cylindrical shells containing fluid flow

In the present study, the geometrically nonlinear vibrations of circular cylindrical shells, subjected to internal fluid flow and to a radial harmonic excitation in the spectral neighbourhood of one of the lowest frequency modes, are investigated for different flow velocities. The shell is modelled by Donnell's nonlinear shell theory, retaining in-plane inertia and geometric imperfections; the fluid is modelled as a potential flow with the addition of unsteady viscous terms obtained by using the time-averaged Navier-Stokes equations. A harmonic concentrated force is applied at mid-length of the shell, acting in the radial direction. The shell is considered to be immersed in an external confined quiescent liquid and to contain a fluid flow, in order to reproduce conditions in previous water-tunnel experiments. For the same reason, complex boundary conditions are applied at the shell ends simulating conditions intermediate between clamped and simply supported ends. Numerical results obtained by using pseudo-arclength continuation methods and bifurcation analysis show the nonlinear response at different flow velocities for (i) a fixed excitation amplitude and variable excitation frequency, and (ii) fixed excitation frequency by varying the excitation amplitude. Bifurcation diagrams of Poincaré maps obtained from direct time integration are presented, as well as the maximum Lyapunov exponent, in order to classify the system dynamics. In particular, periodic, quasi-periodic, sub-harmonic and chaotic responses have been detected. The full spectrum of the Lyapunov exponents and the Lyapunov dimension have been calculated for the chaotic response; they reveal the occurrence of large-dimension hyperchaos.     

Finite-amplitude standing acoustic waves in a cubically nonlinear medium

The acoustic field in a resonator filled with a cubically nonlinear medium is investigated. The field is represented as a linear superposition of two strongly distorted counterpropagating waves. Unlike the case of a quadratically nonlinear medium, the counterpropagating waves in a cubically nonlinear medium are coupled through their mean (over a period) intensities. Free and forced standing waves are considered. Profiles of discontinuous oscillations containing compression and expansion shock fronts are constructed. Resonance curves, which represent the dependences of the mean field intensity on the difference between the boundary oscillation frequency and the frequency of one of the resonator modes, are calculated. The structure of the profiles of strongly distorted “forced” waves is analyzed. It is shown that discontinuities are formed only when the difference between the mean intensity and the detuning takes certain negative values. The discontinuities correspond to the jumps between different solutions to a nonlinear integro-differential equation, which, in the case of small dissipation, degenerates into a third-degree algebraic equation with an undetermined coefficient. The dependence of the intensity of discontinuous standing waves on the frequency of oscillations of the resonator boundary is determined. A nonlinear saturation is revealed: at a very large amplitude of the resonator wall oscillations, the field intensity in the resonator ceases depending on the amplitude and cannot exceed a certain limiting value, which is determined by the nonlinear attenuation at the shock fronts. This intensity maximum is reached when the frequency smoothly increases above the linear resonance. A hysteresis arises, and a bistability takes place, as in the case of a concentrated system at a nonlinear resonance.     

Stationary dark localized modes: discrete nonlinear Schrödinger equations

Various kinds of stationary dark localized modes in discrete nonlinear Schr?dinger equations are considered. A criterion for the existence of such excitations is introduced and an estimation of a localization region is provided. The results are illustrated in examples of the deformable discrete nonlinear Schr?dinger equation, of the model of Frenkel excitons in a chain of two-level atoms, and of the model of a one-dimensional Heisenberg ferromagnetic in the stationary phase approximation. The three models display essentially different properties. It is shown that at an arbitrary amplitude of the background it is impossible to reach strong localization of dark modes. In the meantime, in the model of Frenkel excitons, exact dark compacton solutions are found.     

Equilibrium distribution of the wave energy in a carbyne chain

The steady-state energy distribution of thermal vibrations at a given ambient temperature has been investigated based on a simple mathematical model that takes into account central and noncentral interactions between carbon atoms in a one-dimensional carbyne chain. The investigation has been performed using standard asymptotic methods of nonlinear dynamics in terms of the classical mechanics. In the first-order nonlinear approximation, there have been revealed resonant wave triads that are formed at a typical nonlinearity of the system under phase matching conditions. Each resonant triad consists of one longitudinal and two transverse vibration modes. In the general case, the chain is characterized by a superposition of similar resonant triplets of different spectral scales. It has been found that the energy equipartition of nonlinear stationary waves in the carbyne chain at a given temperature completely obeys the standard Rayleigh–Jeans law due to the proportional amplitude dispersion. The possibility of spontaneous formation of three-frequency envelope solitons in carbyne has been demonstrated. Heat in the form of such solitons can propagate in a chain of carbon atoms without diffusion, like localized waves.     

Approximate step response of a nonlinear hydraulic mount using a simplified linear model

A simplified dynamic stiffness type linear model is used to analytically find the step responses of a nonlinear hydraulic mount in terms of the transmitted force and top chamber pressure. The closed form solution could be efficiently implemented with effective mount parameters, and peak value and the decay curve predictions could provide some insight into the nonlinear behavior. The analytical solutions to an ideal step input correlate well with both numerical simulations (of the same linear model) and measurements when a step-like displacement excitation is applied to fixed and free decoupler mounts.     

Localized and periodic wave patterns for a nonic nonlinear Schrödinger equation

Propagating modes in a class of ‘nonic’ derivative nonlinear Schrödinger equations incorporating ninth order nonlinearity are investigated by application of two key invariants of motion. A nonlinear equation for the squared wave amplitude is derived thereby which allows the exact representation of periodic patterns as well as localized bright and dark pulses in terms of elliptic and their classical hyperbolic limits. These modes represent a balance among cubic, quintic and nonic nonlinearities.     

Instabilities in the main parametric resonance area of a mechanical system with a pendulum

Vibrations of an autoparametric system, composed of a nonlinear mechanical oscillator with an attached damped pendulum, around the principal resonance region, are investigated in this paper. Approximate analytical solutions of the model are determined on the basis of the Harmonic Balance Method (HBM). Correctness of the analytical results is verified by numerical simulations and experimental tests performed on an especially prepared experimental rig. The influence of all essential parameters such as damping, excitation amplitude and frequency, nonlinear stiffness of the spring, on the localisation of the instability region and the system dynamics is presented in the work. Regions of regular system oscillations, chaotic motions, and full rotation of the pendulum are confirmed experimentally.