STOCHASTIC AVERAGING OF STRONGLY NON-LINEAR OSCILLATORS UNDER BOUNDED NOISE EXCITATION

A stochastic averaging method for strongly non-linear oscillators under external and/or parametric excitation of bounded noise is proposed by using the so-called generalized harmonics functions. The method is then applied to study the primary resonance of Duffing oscillator with hardening spring under external excitation of bounded noise. The stochastic jump and its bifurcation of the system are observed and explained by using the stationary probability density of amplitude and phase. Subsequently, the method is applied to study the dynamical instability and parametric resonance of Duffing oscillator with hardening spring under parametric excitation of bounded noise. The primary unstable region is delineated by evaluating the Lyapunov exponent of linearized system, and the response and jump of non-linear system around the unstable region are examined by using the sample functions and stationary probability density of amplitude and phase.     

Stability Analysis of an Inverted Pendulum Subjected to Combined High Frequency Harmonics and Stochastic Excitations

Stability of vertical upright position of an inverted pendulum with its suspension point subjected to high frequency harmonics and stochastic excitations is investigated. Two classes of excitations, i.e., combined high frequency harmonic excitation and Gaussian white noise excitation, and high frequency bounded noise excitation, respectively, are considered. Firstly, the terms of high frequency harmonic excitations in the equation of motion of the system can be set equivalent to nonlinear stiffness terms by using the method of direct separation of motions. Then the stochastic averaging method of energy envelope is used to derive the averaged Ito stochastic differential equation for system energy. Finally, the stability with probability 1 of the system is studied by using the largest Lyapunov exponent obtained from the averaged Ito stochastic differential equation. The effects of system parameters on the stability of the system are discussed, and some examples are given to illustrate the efficiency of the proposed procedure.     

Asymptotic stability with probability one of MDOF nonlinear oscillators with fractional derivative damping

In this paper, the asymptotic stability with probability one of multi-degree-of-freedom (MDOF) nonlinear oscillators with fractional derivative damping parametrically excited by Gaussian white noises is investigated. A stochastic averaging method and the Khasminskii’s procedure are employed to evaluate the largest Lyapunov exponent, whose sign determines the stability of the system. As an example, two coupled nonlinear oscillators with fractional derivative damping is worked out to demonstrate the proposed procedure and to examine the effect of fractional order on the stochastic stability of system. In particular, the case of factional order more than 1 is studied for the first time.     

Stationary Response of Lotka-Volterra System with Real Noises

A stochastic version of Lotka-Volterra model subjected to real noises is proposed and investigated. The approximate stationary probability densities for both predator and prey are obtained analytically. The original system is firstly transformed to a pair of It o stochastic differential equations. The Ito formula is then carried out to obtain the It o stochastic differential equation for the period orbit function. The orbit function is considered as slowly varying process under reasonable assumptions. By applying the stochastic averaging method to the orbit function in one period, the averaged Ito stochastic differential equation of the motion orbit and the corresponding Fokker-Planck equation are derived. The probability density functions of the two species are thus formulated. Finally, a classical real noise model is given as an example to show the proposed approximate method. The accuracy of the proposed procedure is verified by Monte Carlo simulation.     

Higher-order stochastic averaging to study stability of a fractional viscoelastic column

Stochastic stability of a fractional viscoelastic column axially loaded by a wideband random force is investigated by using the method of higher-order stochastic averaging. By modelling the wideband random excitation as Gaussian white noise and real noise and assuming the viscoelastic material to follow the fractional Kelvin–Voigt constitutive relation, the motion of the column is governed by a fractional stochastic differential equation, which is justifiably and uniformly approximated by an averaged system of Itô stochastic differential equations. Analytical expressions are obtained for the moment Lyapunov exponent and the Lyapunov exponent of the fractional system with small damping and weak random fluctuation. The effects of various parameters on the stochastic stability of the system are discussed.     

Stochastic stability of a fractional viscoelastic column under bounded noise excitation

The stability of a viscoelastic column under the excitation of stochastic axial compressive load is investigated in this paper. The material of the column is modeled using a fractional Kelvin–Voigt constitutive relation, which leads to that the equation of motion is governed by a stochastic fractional equation with parametric excitation. The excitation is modeled as a bounded noise, which is a realistic model of stochastic fluctuation in engineering applications. The method of stochastic averaging is used to approximate the responses of the original dynamical system by a new set of averaged variables which are diffusive Markov vector. An eigenvalue problem is formulated from the averaged equations, from which the moment Lyapunov exponent is determined for the column system with small damping and weak excitation. The effects of various parameters on the stochastic stability and significant parametric resonance are discussed and confirmed by simulation results.     

A stochastic epidemic model on homogeneous networks

In this paper, a stochastic SIS epidemic model on homogeneous networks is considered. The largest Lyapunov exponent is calculated by Oseledec multiplicative ergodic theory, and the stability condition is determined by the largest Lyapunov exponent. The probability density function for the proportion of infected individuals is found explicitly, and the stochastic bifurcation is analysed by a probability density function. In particular, the new basic reproductive number R*, that governs whether an epidemic with few initial infections can become an endemic or not, is determined by noise intensity. In the homogeneous networks, despite of the basic productive number R₀>1, the epidemic will die out as long as noise intensity satisfies a certain condition.     

Stochastic responses of Duffing–Van der Pol vibro-impact system under additive colored noise excitation

A response analysis procedure is developed for a vibro-impact system excited by colored noise. The non-smooth transformation is used to convert the vibro-impact system into a new system without impact term. With the help of the modified quasi-conservative averaging, the total energy of the new system can be approximated as a Markov process, and the stationary probability density function （PDF） of the total energy is derived. The response PDFs of the original system are obtained using the analytical solution of the stationary PDF of the total energy. The validity of the theoretical results is tested through comparison with the corresponding simulation results. Moreover, stochastic bifurcations are also explored.     

乘性色噪声激励下三稳态van der Pol-Duffing振子随机P-分岔

研究了乘性色噪声作用下三稳态van der Pol-Duffing振子的随机P-分岔问题. 首先应用随机平均法得到系统振动幅值稳态概率密度函数的表达式, 进而应用奇异性理论, 得到刻画随机P-分岔发生的临界参数条件的转迁集以及系统存在的典型稳态概率密度曲线, 并通过Monte-Carlo数值模拟进行了验证. 以此为基础讨论了噪声强度、相关时间、系统线性阻尼系数对随机P-分岔和系统稳态响应行为的影响.     

Stochastic bifurcations of generalized Duffing–van der Pol system with fractional derivative under colored noise

The stochastic bifurcation of a generalized Duffing–van der Pol system with fractional derivative under color noise excitation is studied. Firstly, fractional derivative in a form of generalized integral with time-delay is approximated by a set of periodic functions. Based on this work, the stochastic averaging method is applied to obtain the FPK equation and the stationary probability density of the amplitude. After that, the critical parameter conditions of stochastic P-bifurcation are obtained based on the singularity theory. Different types of stationary probability densities of the amplitude are also obtained. The study finds that the change of noise intensity, fractional order, and correlation time will lead to the stochastic bifurcation.     

On stochastic complex beam-beam interaction models with Gaussian colored noise

This paper is to continue our study on complex beam-beam interaction models in particle accelerators with random excitations Y. Xu, W. Xu, G.M. Mahmoud, On a complex beam-beam interaction model with random forcing [Physica A 336 (2004) 347-360]. The random noise is taken as the form of exponentially correlated Gaussian colored noise, and the transition probability density function is obtained in terms of a perturbation expansion of the parameter. Then the method of stochastic averaging based on perturbation technique is used to derive a Fokker-Planck equation for the transition probability density function. The solvability condition and the general transforms using the method of characteristics are proposed to obtain the approximate expressions of probability density function to order ε.Also the exact stationary probability density and the first and second moments of the amplitude are obtained, and one can find when the correlation time equals to zero, the result is identical to that derived from the Stratonovich-Khasminskii theorem for the same model under a broad-band excitation in our previous work.     

STOCHASTIC AVERAGING OF DUHEM HYSTERETIC SYSTEMS

The response of Duhem hysteretic system to externally and/or parametrically non-white random excitations is investigated by using the stochastic averaging method. A class of integrable Duhem hysteresis models covering many existing hysteresis models is identified and the potential energy and dissipated energy of Duhem hysteretic component are determined. The Duhem hysteretic system under random excitations is replaced equivalently by a non-hysteretic non-linear random system. The averaged Ito's stochastic differential equation for the total energy is derived and the Fokker-Planck-Kolmogorov equation associated with the averaged Ito's equation is solved to yield stationary probability density of total energy, from which the statistics of system response can be evaluated. It is observed that the numerical results by using the stochastic averaging method is in good agreement with that from digital simulation.     

Delay induced transitions in an asymmetry bistable system and stochastic resonance

The delay Fokker-Planck equation is given for an asymmetry bistable system with correlated Gaussian white noises. The small delay approximation based on the probability density approach is used and the approximate stationary probability density function is obtained. The phenomenon of delay induced transitions is found. When a weak periodic signal is added, the phenomenon of stochastic resonance is investigated. Expression of the signal-to-noise ratio (SNR) is obtained by using the two-state theory. It is sh...     

IN-SERVICE IDENTIFICATION OF NON-LINEAR DAMPING FROM MEASURED RANDOM VIBRATION

A non-linearly damped single-degree-of-freedom (s.d.o.f.) system under broadband random excitation is considered. A procedure for in-service identification of the damping characteristic from measured stationary response is described. The procedure is based on the stochastic averaging method. The explicit analytical solution is obtained for the integral equation, which relates the desired damping characteristics to the apparent force in the shortened equation for the slowly varying response amplitude, and thus to the measured probability density of the amplitude. The approach is of a non-parametric nature, which makes it convenient for testing hypotheses of damping mechanisms from measured random vibration data. Extensive results of numerical tests for the procedure are presented.     

Maximal Lyapunov exponent and rotation numbers for two coupled oscillators driven by real noise

Asymptotic expansions for the exponential growth rate, known as the Lyapunov exponent, and rotation numbers for two coupled oscillators driven by real noise are constructed. Such systems arise naturally in the investigation of the stability of steady-state motions of nonlinear dynamical systems and in parametrically excited linear mechanical systems. Almost-sure stability or instability of dynamical systems depends on the sign of the maximal Lyapunov exponent. Stability conditions are obtained under various assumptions on the infinitesimal generator associated with real noise provided that the natural frequencies are noncommensurable. The results presented here for the case of the infinitesimal generator having a simple zero eigenvalue agree with recent results obtained by stochastic averaging, where approximate ItÔ equations in amplitudes and phases are obtained in the sense of weak convergence.Dedicated to Thomas K. Caughey on the occasion of his 65th birthday.     

Response of a Duffing—Rayleigh system with a fractional derivative under Gaussian white noise excitation

In this paper,we consider the response analysis of a Duffing-Rayleigh system with fractional derivative under Gaussian white noise excitation.A stochastic averaging procedure for this system is developed by using the generalized harmonic functions.First,the system state is approximated by a diffusive Markov process.Then,the stationary probability densities are derived from the averaged Ito stochastic differential equation of the system.The accuracy of the analytical results is validated by the results from the Monte Carlo simulation of the original system.Moreover,the effects of different system parameters and noise intensity on the response of the system are also discussed.     

Stochastic Analysis of Predator–Prey Models under Combined Gaussian and Poisson White Noise via Stochastic Averaging Method

In the present paper, the statistical responses of two-special prey–predator type ecosystem models excited by combined Gaussian and Poisson white noise are investigated by generalizing the stochastic averaging method. First, we unify the deterministic models for the two cases where preys are abundant and the predator population is large, respectively. Then, under some natural assumptions of small perturbations and system parameters, the stochastic models are introduced. The stochastic averaging method is generalized to compute the statistical responses described by stationary probability density functions (PDFs) and moments for population densities in the ecosystems using a perturbation technique. Based on these statistical responses, the effects of ecosystem parameters and the noise parameters on the stationary PDFs and moments are discussed. Additionally, we also calculate the Gaussian approximate solution to illustrate the effectiveness of the perturbation results. The results show that the larger the mean arrival rate, the smaller the difference between the perturbation solution and Gaussian approximation solution. In addition, direct Monte Carlo simulation is performed to validate the above results.     

Response of SDOF nonlinear oscillators with lightly fractional derivative damping under real noise excitations

This paper deals with the response of single-degree-of-freedom (SDOF) strongly nonlinear oscillator with lightly fractional derivative damping to external and (or) parametric real noise excitations. First, the state vector of the displacement and the velocity is approximated by one-dimensional time-homogeneous diffusive Markov process of amplitude through using the stochastic averaging method. Then, the stationary probability density of amplitude is obtained by solving the Fokker-Planck-Kolmogorov (FPK) equation associated with the averaged It? equation of amplitude, in which the Fourier series expansions are used to obtain the explicit expressions of the drift and diffusion coefficients. Finally, the response of a Duffing oscillator with lightly fractional derivative damping under external and parametric real noise excitations is evaluated by using the proposed procedure and compared with that from the Monte Carlo simulation of original system.     

Lyapunov exponent diagrams of a 4-dimensional Chua system

We report numerical results on the existence of periodic structures embedded in chaotic and hyperchaotic regions on the Lyapunov exponent diagrams of a 4-dimensional Chua system. The model was obtained from the 3-dimensional Chua system by the introduction of a feedback controller. Both the largest and the second largest Lyapunov exponents were considered in our colorful Lyapunov exponent diagrams, and allowed us to characterize periodic structures and regions of chaos and hyperchaos. The shrimp-shaped periodic structures appear to be malformed on some of Lyapunov exponent diagrams, and they present two different bifurcation scenarios to chaos when passing the boundaries of itself, namely via period-doubling and crisis. Hyperchaos-chaos transition can also be observed on the Lyapunov exponent diagrams for the second largest exponent.     

Chaos Control for a Single Pendulum Damping System with Proper Random Phase Under Forced Oscillation

This paper deals with the effect of phase of a Gaussian white noise for a single pendulum system with damping. Based on the Khasminskii’s formulation of spherical coordinate and the extension of Wedig’s algorithm for linear stochastic system, the largest Lyapunov exponent is computed. Due to the change of the sign for the largest Lyapunov exponent, the chaotic behavior of the system is suppressed. Also Poincaré surface of section, phase portrait and the time evolution are analyzed to confirm the stability of the system, which demonstrates the above control methods are effective.