Principal parametric resonances of a slender cantilever beam subject to axial narrow-band random excitation of its base

The non-linear integro-differential equations of motion for a slender cantilever beam subject to axial narrow-band random excitation are investigated. The method of multiple scales is used to determine a uniform first-order expansion of the solution of equations. According to solvability conditions, the non-linear modulation equations for the principal parametric resonance are obtained. Firstly, The largest Lyapunov exponent which determines the almost sure stability of the trivial solution is quantificationally resolved, in which, the modified Bessel function of the first kind is introduced. Results show that the increase of the bandwidth facilitates the almost sure stability of the trivial response and stabilizes the system for a lower acceleration oscillating amplitude but intensifies the instability of the trivial response for a higher one. Secondly, the first and second order non-trivial steady state response of the system is obtained by perturbation method and the corresponding amplitude–frequency curves are calculated when the bandwidth is very small. Results show that the effective non-linearity of whether the amplitude expectation of the first order steady state response or the amplitude expectation of the second order steady state response is of the hardening type for the first mode, whereas for the second mode the effective non-linearity of whether the amplitude expectation of the first order steady state response or the amplitude expectation of the second order steady state response is of the softening type. Finally, the stochastic jump and bifurcation is investigated for the first and second modal parametric principal resonance. The basic jump phenomena indicate that, under the conditions of system parameters with a smaller bandwidth, the most probable motion is around the non-trivial branch of the amplitude response curve, whereas with a higher bandwidth, the most probable motion is around the trivial one of the amplitude response curve. However, the stochastic jump is sometimes more sensitive to the change of the bandwidth, in other words, a small change of bandwidth may induce a series of stochastic jump and bifurcation.     

一类非线性系统在随机激励下的分叉

本文研究一类阻尼为线性，弹性恢复力为非线性的振动系统在随机外部激励作用下的随机分叉。文中采用广义稳态势和方法，求解系统响应的稳态联合概率密度函数。在此基础上根据由不变测度定义的随机分叉，讨论了具有权式分叉的确定性非线性系统在随机扰动下分叉行为。     

Jump and bifurcation of Duffing oscillator under narrow-band excitation

The jump and bifurcation of Duffing oscillator with hardening spring subject to narrow-band random excitation are systematically and comprehensively examined. It is shown that, in a certain domain of the space of the oscillator and excitation parameters, there are two types of more probable motions in the stationary response of the Duffing oscillator and jumps may occur. The jump is a transition of the response from one more probable motion to another or vise versa. Outside the domain the stationary response is either nearly Gaussian or like a diffused limit cycle. As the parameters change across the boundary of the domain the qualitative behavior of the stationary response changes and it is a special kind of bifurcation. It is also shown that, for a set of specified parameters, the statistics are unique and they are independent of initial condition. It is pointed out that some previous results and interpretations on this problem are incorrect. The project supported by National Natural Science Foundation of China     

Necessary and sufficient conditions for existence of stationary solutions of some nonlinear stochastic systems

At the state of statistical stationarity, the response of a nonlinear system under multiplicative random excitations can be either trivial or non-trivial, depending on the spectral levels of the excitations and the values of certain system parameters. Assuming that the random excitations are Gaussian white noises, the two types of response may be investigated by way of their stationary densities, which are obtainable for first order dynamical systems and for higher order dynamical systems belonging to the class of generalized stationary potential. Alternatively, the Lyapunov exponents can be computed for perturbation from either the trivial or non-trivial solution, since a negative sign for the greatest Lyapunov exponent provides both the necessary and sufficient conditions for the stability of sample functions with probability one. It is shown in two specific examples, that the boundary at which the greatest Lyapunov exponent changes its sign coincides with the boundary for regularity (or being normalizable) for the probability density in both the trivial and non-trivial solutions. Thus, the stability conditions in the strong sense of probability one and the weak sense in distribution are identical in these cases.     

Simultaneous combination and 1:3:5 internal resonances in a parametrically excited beam-mass system

The nonlinear dynamics of a base-excited slender beam carrying a lumped mass subjected to simultaneous combination parametric resonance of sum and difference type along with 1:3:5 internal resonances is investigated. Method of normal form is applied to the governing nonlinear temporal differential equation of motion to obtain a set of first-order differential equations which are used to obtain the steady-state, periodic, quasi-periodic and chaotic responses for different control parameters viz., amplitude and frequency of external excitation and damping. Frequency response, phase portraits, time spectra and bifurcation diagram are plotted to visualize the system behaviour with variation in the control parameters. Here, two distinct zones of trivial instability, blue sky catastrophe phenomena, jump down phenomena, simultaneous occurrence of periodic and chaotic orbits, period doubling of the mixed-mode periodic orbits leading to chaos, attractor merging crisis, boundary crisis, type II and on-off intermittencies are observed. Bifurcation diagram is plotted to facilitate the designer to choose a safe operating zone.     

A new study of chaotic behavior and the existence of Feigenbaum’s constants in fractional-degree Yin–Yang Hénon maps

hing dynamics in a square tank are numerically investigated when the tank is subjected to horizontal, narrowband random ground excitation. The natural frequencies of the two predominant sloshing modes are identical and therefore 1:1 internal resonance may occur. Galerkin’s method is applied to derive the modal equations of motion for nonlinear sloshing including higher modes. The Monte Carlo simulation is used to calculate response statistics such as mean square values and probability density functions (PDFs). The two predominant modes exhibit complex phenomena including “autoparametric interaction” because they are nonlinearly coupled with each other. The mean square responses of these two modes and the liquid elevation are found to differ significantly from those of the corresponding linear model, depending on the characteristics of the random ground excitation such as bandwidth, center frequency and excitation direction. It is found that the direction of the excitation is a significant factor in predicting the mean square responses. The frequency response curves for the same system subjected to equivalent harmonic excitation are also calculated and compared with the mean square responses to further explain the phenomena. Changing the liquid level causes the peak of the mean square response to shift. Furthermore, the risk of the liquid overspill from the tank is discussed by showing the three-dimensional distribution charts of the mean square responses of liquid elevations.     

Effects of excitation probability distribution on system responses

For a system subjected to a random excitation, the probability distribution of the excitation may affect behaviors of the system responses. Such effects are investigated for a variety of dynamical systems, including a linear oscillator, an oscillator of cubic non-linearity in both damping and stiffness, and a non-linear oscillator of the van der Pol type. The random excitations are assumed to be stationary stochastic processes, sharing the same spectral density, but with different probability distributions. Each excitation process is generated by passing a Brownian motion process through a non-linear filter, which is governed by an Ito stochastic differential equation. Monte Carlo simulations are carried out to obtain the transient and stationary properties of the system response in each case. It is shown that, under different excitations, the transient behaviors of the system response can be markedly different. The differences tend to reduce, however, as time of exposure to the excitations increases and the system reaches the stationary state.     

Dynamical responses of airfoil models with harmonic excitation under uncertain disturbance

In this paper, we investigate nonlinear dynamical responses of two-degree-of-freedom airfoil (TDOFA) models driven by harmonic excitation under uncertain disturbance. Firstly, based on the deterministic airfoil models under the harmonic excitation, we introduce stochastic TDOFA models with the uncertain disturbance as Gaussian white noise. Subsequently, we consider the amplitude–frequency characteristic of deterministic airfoil models by the averaging method, and also the stochastic averaging method is applied to obtain the mean-square response of given stochastic TDOFA systems analytically. Then, we carry out numerical simulations to verify the effectiveness of the obtained analytic solution and the influence of harmonic force on the system response is studied. Finally, stochastic jump and bifurcation can be found through the random responses of system, and probability density function and time history diagrams can be obtained via Monte Carlo simulations directly to observe the stochastic jump and bifurcation. The results show that noise can induce the occurrence of stochastic jump and bifurcation, which will have a significant impact on the safety of aircraft.     

Nonlinearly damped systems under simultaneous broad-band and harmonic excitations

The probability distribution of the response of a nonlinearly damped system subjected to both broad-band and harmonic excitations is investigated. The broad-band excitation is additive, and the harmonic excitations can be either additive or multiplicative. The frequency of a harmonic excitation can be either near or far from a resonance frequency of the system. The stochastic averaging method is applied to obtain the Itô type stochastic differential equations for an averaged system described by a set of slowly varying variables, which are approximated as components of a Markov vector. Then, a procedure based on the concept of stationary potential is used to obtain the exact stationary probability density for a class of such averaged systems. For those systems not belonging to this class, approximate solutions are obtained using the method of weighted residuals. Application of the exact and approximate solution procedures are illustrated in two specific cases, and the results are compared with those obtained from Monte Carlo simulations.     

Responses of Duffing oscillator with fractional damping and random phase

This paper aims to investigate dynamic responses of stochastic Duffing oscillator with fractional-order damping term, where random excitation is modeled as a harmonic function with random phase. Combining with Lindstedt–Poincaré (L–P) method and the multiple-scale approach, we propose a new technique to theoretically derive the second-order approximate solution of the stochastic fractional Duffing oscillator. Later, the frequency–amplitude response equation in deterministic case and the first- and second-order steady-state moments for the steady state in stochastic case are presented analytically. We also carry out numerical simulations to verify the effectiveness of the proposed method with good agreement. Stochastic jump and bifurcation can be found in the figures of random responses, and then we apply Monte Carlo simulations directly to obtain the probability density functions and time response diagrams to find the stochastic jump and bifurcation. The results intuitively show that the intensity of the noise can lead to stochastic jump and bifurcation.     

Vertical dynamics of a pantograph carbon-strip suspension under stochastic contact-force excitation

In electric trains, a pantograph is mounted on the roof of the train to collect power through contact with an overhead catenary wire. The carbon-strip suspension of a pantograph, along which contact with the catenary occurs, is subjected to harmonic and stochastic contact-force excitations. In this paper, vertical dynamics of the carbon-strip suspension is studied with an aim of improving the reliability and safety of running trains. A single-degree-of-freedom model of the carbon-strip suspension with nonlinear stiffness is developed using parameter values of the DSA X pantograph. Through stochastic averaging, a Fokker–Planck–Kolmogorov equation governing the stationary response of the carbon-strip suspension is set up. Based on the transition probability density of the stationary response, it is found that random jumps and bifurcations in the carbon-strip motion can occur. The possibility of motion bifurcations and the frequency of random jumps warrant consideration in advanced design of carbon-strip suspensions.     

Conditional Fluid-Particle Accelerations in Turbulence

Simple closures for average fluid-particle accelerations, conditional on fixed local fluid velocity, are considered in isotropic, homogeneous and stationary turbulence using exact probability density transport equations and are compared with direct numerical simulations (DNS). Such accelerations are common ingredients in Lagrangian stochastic models for fluid-particle trajectories in turbulence. One-particle accelerations are essentially trivial, so the focus is on two-particle relative accelerations, which are important in the relative dispersion process. The closure is simply a quadratic form in the velocity variable and this special form also defines the Eulerian velocity probability density function (pdf), and comparisons with DNS (for grids up to 512³) of both the acceleration closure and velocity pdf's are encouraging. Received 2 June 1997 and accepted 29 December 1997     

Two degree-of-freedom oscillator coupled to a non-ideal source

In the paper the one-mass two degree-of-freedom system with non-ideal excitation is considered. The resonance motion of the system is investigated. The mathematical model of the system contains three coupled second order differential equations. In the paper an analytical solving procedure is developed. The steady-state motion and the criteria for stability of solutions are developed. Two special cases of motion depending on the frequency properties of the system are studied. When the frequency properties in both orthogonal direction are equal there is only one resonance. If the frequency in one direction is two times higher than in other two different resonances occur: one in x and the other in y direction. The conditions for jump phenomena and for Sommerfeld effect are presented. The analytically obtained solutions are compared with numerical ones. They show good agreement.     

Stochastic bifurcations in a nonlinear tri-stable energy harvester under colored noise

In this paper, the stochastic bifurcations and the performance analysis of a strongly nonlinear tri-stable energy harvesting system with colored noise are investigated. Using the stochastic averaging method, the averaged Fokker–Plank–Kolmogorov equation and the stationary probability density (SPD) of the amplitude are obtained, respectively. Meanwhile, the Monte Carlo simulations are performed to verify the effectiveness of the theoretical results. D-bifurcation is studied through the largest Lyapunov exponent calculations, which implies the system undergoes D-bifurcation twice with increasing the nonlinear stiffness coefficients. The effects of the nonlinear stiffness coefficients, noise intensity and correlation time on P-bifurcation are discussed by the qualitative changes of the SPD. Moreover, the relationship between D- and P-bifurcation is explored. If the strength of stochastic jump has obvious gap with respect to the two statuses before and after the occurrence of P-bifurcation, the D-bifurcation will happen, too. Finally, the performance and the capability of harvesting energy from ambient random excitation are analyzed.     

The Study of a Parametrically Excited Nonlinear Mechanical System with the Continuation Method

The dynamic behavior of a one-degree-of-freedom, parametrically excited nonlinear system is investigated. The Galerkin method is applied to the principal and fundamental parameteric resonance of the system. The continuation method is used to study the change of harmonic oscillation with respect to the variation of excitation frequency. The numerical stability analysis of the trivial solution is carried out and the stable and unstable regions of the trivial solution are given. They are found to agree with the results obtained by the analytical method of Galerkin. Periodic solutions are traced and the coexistence of multi-periodic solutions is observed With the change of excitation frequency the large amplitude periodic-2 oscillation is found to be in the same closed branch with the small amplitude periodic-2 solution. In addition, the bifurcation pattern of the trivial solution is found to change from subcritical Hopf bifurcation into supercritical Hopf bifurcation with the increase of excitation amplitude. Combined with the conventional numerical integration method, new complex dynamic behavior is detected.     

The nonstationary effects on a softening duffing oscillator

This paper presents a numerical study for the bifurcations of a softening Duffing oscillator subjected to stationary and nonstationary excitation. The nonstationary inputs used are linear functions of time. The bifurcations are the results of either a single control parameter or two control parameters that are constrained to vary in a selected direction on the plane of forcing amplitude and forcing frequency. The results indicate: 1. Delay (memory, penetration) of nonstationary bifurcations relative to stationary bifurcations may occur. 2. The nonstationary trajectories jump into the neighboring stationary trajectories with possible overshoots, while the stationary trajectories transit smoothly. 3. The nonstationary penetrations (delays) are compressed to zero with an increasing number of iterations. 4. The nonstationary responses converge through a period-doubling sequence to a nonstationary limit motion that has the characteristics of chaotic motion. The Duffing oscillator has been used as an example of the existence of broad effects of nonstationary (time dependent) and codimensional (control parameter variations in the bifurcation region) inputs which markedly modify the dynamical behavior of dynamical systems.     

Exact stationary probability density functions for non-linear systems under Poisson white noise excitation

The paper presents exact stationary probability density functions for systems under Poisson white noise excitation. Two different solution methods are outlined. In the first one, a class of non-linear systems is determined whose state vector is a memoryless transformation of the state vector of a linear system. The second method considers the generalized Fokker-Planck (Kolmogorov-forward) equation. Non-linear system functions are identified such that the stationary solution of the system admits a prescribed stationary probability density function. Both methods make use of the stochastic integro-differential equations approach. This approach seems to have some computational advantages for the determination of exact stationary probability density functions when compared to the stochastic differential equations approach.     

Response Statistics of Three-Degree-of-Freedom Nonlinear System to Narrow-Band Random Excitation

The principal resonance of a 3-DOF nonlinear system to narrow-band random external excitations is investigated. The method of multiple scales is used to derive the equations for modulation of amplitude and phase. The behavior, stability and bifurcation of steady-state responses are studied by means of qualitative analysis. The effects of damping, detuning, and excitation intensity on responses are analyzed. The theoretical analyses are verified by numerical results. Both theoretical analyses and numerical simulations show that when the intensity of the random excitation increases, the nontrivial steady state solution may change from a limit cycle to a diffused limit cycle. Under some conditions, co-existence of two kinds of stable steady-state solutions, saturation and jump phenomena may occur. The stationary probability density function of responses for the co-existence case is obtained approximately.     

Nonlinear dynamics of excited plate immersed in granular matter

The interaction between granular matter and the elastic body is a complex issue due to the complex properties of granular matter. An experiment involving a sinusoidally excited plate buried in glass bead particles contained in a box is conducted. The motion behavior of the plate is observed and recorded by the strain gauge. The amplitude–frequency and phase–frequency curves are recorded to study the natural property of the plate in granular matter. In this experiment, jump phenomena are found in both the amplitude–frequency and phase–frequency planes in circumstances with smaller particle sizes, lower buried depths, and larger amplitudes of the excitation force. Otherwise, the period-doubling bifurcation, especially 3T, is found with the increase in the excitation force. These bifurcations usually occur in specific buried depth and excitation frequency band and require smaller particle sizes. The experiments with random-shaped particles exhibit no-jump phenomenon, but period-doubling bifurcation and chaos. These phenomena are sensitive to parameters and closely related to the varying process of the excitation frequency and force. Reasonable mechanisms are summarized qualitatively through some of our recent researches in this paper.     

Deterministic and Stochastic Response of Nonlinear Coupled Bending-Torsion Modes in a Cantilever Beam

The purpose of this study is to understand the main differences between the deterministic and random response characteristics of an inextensible cantilever beam (with a tip mass) in the neighborhood of combination parametric resonance. The excitation is applied in the plane of largest rigidity such that the bending and torsion modes are cross-coupled through the excitation. In the absence of excitation, the two modes are also coupled due to inertia nonlinearities. For sinusoidal parametric excitation, the beam experiences instability in the neighborhood of the combination parametric resonance of the summed type, i.e., when the excitation frequency is in the neighborhood of the sum of the first bending and torsion natural frequencies. The dependence of the response amplitude on the excitation level reveals three distinct regions: nearly linear behavior, jump phenomena, and energy transfer. In the absence of nonlinear coupling, the stochastic stability boundaries are obtained in terms of sample Lyapunov exponent. The response statistics are estimated using Monte Carlo simulation, and measured experimentally. The excitation center frequency is selected to be close to the sum of the bending and torsion mode frequencies. The beam is found to experience a single response, two possible responses, or non-stationary responses, depending on excitation level. Experimentally, it is possible to obtain two different responses for the same excitation level by providing a small perturbation to the beam during the test.