Random excitation of coupled oscillators with single and simultaneous internal resonances

The primary objective of this paper is to examine the random response characteristics of coupled nonlinear oscillators in the presence of single and simultaneous internal resonances. A model of two coupled beams with nonlinear inertia interaction is considered. The primary beam is directly excited by a random support motion, while the coupled beam is indirectly excited through autoparametric coupling and parametric excitation. For a single one-to-two internal resonance, we used Gaussian and non-Gaussian closures, Monte Carlo simulation, and experimental testing to predict and measure response statistics and stochastic bifurcation in the mean square. The mean square stability boundaries of the coupled beam equilibrium position are obtained by a Gaussian closure scheme. The stochastic bifurcation of the coupled beam is predicted theoretically and experimentally. The stochastic bifurcation predicted by non-Gaussian closure is found to take place at a lower excitation level than the one predicted by Gaussian closure and Monte Carlo simulation. It is also found that above a certain excitation level, the solution obtained by non-Gaussian closure reveals numerical instability at much lower excitation levels than those obtained by Gaussian and Monte Carlo approaches. The experimental observations reveal that the coupled beam does not reach a stationary state, as reflected by the time evolution of the mean square response. For the case of simultaneous internal resonances, both Gaussian and non-Gaussian closures fail to predict useful results, and attention is focused on Monte Carlo simulation and experimental testing. The effects of nonlinear coupling parameters, internal detuning ratios, and excitation spectral density level are considered in both investigations. It is found that both studies reveal common nonlinear features such as bifurcations in the mean square responses of the coupled beam and modal interaction in the neighborhood of internal resonances. Furthermore, there is an upper limit for the excitation level above which the system experiences unbounded response in the neighborhood of simultaneous internal resonances.     

Planar and non-planar non-linear dynamics of suspended cables under random in-plane loading-I. Single internal resonance

The non-linear interaction of the in-plane and out-of-plane motions of a suspended cable in the neighbourhood of 2:1 internal resonance under random loading is studied. The random loading acts externally on the in-plane mode, while the out-of-plane mode is non-linearly coupled with the in-plane mode. Any non-trivial motion of the out-of-plane mode is mainly due to this non-linear coupling, which becomes significant in the neighbourhood of internal resonance. The response statistics are estimated by employing the Fokker-Planck equation together with Gaussian and non-Gaussian closures. Monte-Carlo simulation is also used for numerical verification. Away from the internal resonance condition, the response is governed by the inplane motion, and the non-Gaussian closure solution is found to be in good agreement with numerical simulation results. The stochastic bifurcation of the out-of-plane mode is predicted by Gaussian and non-Gaussian closures, and by Monte-Carlo simulation. The non-Gaussian closure can only predict bounded solutions within a limited region. The on-off intermittency of the second mode is observed in the Monte-Carlo simulation over a small range of excitation level. The influence on response statistics of excitation level and cable parameters, such as internal detuning, damping ratios, and sag-to-span ratio, is reported.     

Random excitation of nonlinear coupled oscillators

This paper presents the experimental results of random excitation of a nonlinear two-degree-of-freedom system in the neighborhood of internal resonance. The random signals of the excitation and response coordinates are processed to estimate the mean squares, autocorrelation functions, power spectral densities, and probability density functions. The results are qualitatively compared with those predicted by the Fokker-Planck equation together with a non-Gaussian closure scheme. The effects of system damping ratios, nonlinear coupling parameter, internal detuning ratio, and excitation spectral density level are considered in both studies except the effect of damping ratios is not considered in the experimental investigation. Both studies reveal similar dynamic features such as autoparametric absorber effect and stochastic instability of the coupled system. The experimental results show that the autocorrelation function of the coupled system has the feature of ultra narrow band process and degenerates to a periodic one as the internal detuning departs from the exact internal resonance condition. The measured probability density functions of the response of the main system suggests that the Gaussian representation is sufticient as long as the excitation level is relatively low in the neighborhood of the system internal resonance condition.     

Parametric random excitation of nonlinear coupled oscillators

The stochastic bifurcation and response statistics of nonlinear modal interaction under parametric random excitation are studied analytically, numerically and experimentally. Two basic definitions of stochastic bifurcation are first introduced. These are bifurcation in distribution and bifurcation in moments. bifurcation in moments is examined for the case of a coupled oscillator subjected to parametric filtered white noise. The center frequency of the excitation is selected to be close to either twice the first mode or second mode natural frequencies or the sum of the two. The stochastic bifurcation in moments is predicted using the Fokker-Planck equation together with gaussian and non-Gaussian closures and numerically using the Monte Carlo simulation. When one mode is parametrically excited it transfers energy to the other mode due to nonlinear modal interaction. The Gaussian closure solution gives close results to those predicted numerically only in regions well remote from bifurcation points. However, bifurcation points predicted by the non-Gaussian closure are in good agreement with those estimated by numerical simulation. Depending on the excitation level, the probability density of the excited mode is strongly non-Gaussian and exhibits multi-maxima as predicted by Monte Carlo simulation. Experimental tests are carried out at relatively low excitation levels. In the neighborhood of stochastic bifurcation in mean square the measured results exhibit different regimes of response characteristics including zero motion and occasional small random motion regimes. These two regimes are characterized by the phenomenon of on-off intermittency. Both regimes overlap and thus it is difficult to locate experimentally the bifurcation point.     

非高斯随机激励下滞迟系统响应与首次穿越失效

针对由有界噪声、泊松白噪声和高斯白噪声共同构成的非高斯随机激励,通过Monte Carlo数值模拟方法研究了此激励作用下双线性滞迟系统和Bouc-Wen滞迟系统这两类经典滞迟系统的稳态响应与首次穿越失效时间。一方面,分析了有界噪声和泊松白噪声这两种分别具有连续样本函数和非连续样本函数的非高斯随机激励,在不同激励参数条件下对双线性滞迟系统和Bouc-Wen滞迟系统的稳态响应概率密度、首次穿越失效时间概率密度及其均值的不同影响;另一方面,揭示了在这类非高斯随机激励荷载作用下,双线性滞迟系统的首次穿越失效时间概率密度将出现与Bouc-Wen滞迟系统的单峰首次穿越失效时间概率密度截然不同的双峰形式。     

Nonlinear Random Response of Ocean Structures Using First- and Second-Order Stochastic Averaging

This paper deals with the dynamic response of nonlinear elastic structure subjected to random hydrodynamic forces and parametric excitation using a first- and second-order stochastic averaging method. The governing equation of motion is derived by using Hamilton's principle, taking into account inertia and curvature nonlinearities and work done due to hydrodynamic forces. Within the framework of first-order stochastic averaging, the system response statistics and stability boundaries are obtained. Unfortunately, the effects of nonlinear inertia and curvature are not reflected in the final results, since the contribution of these nonlinearities is lost during the averaging process. In the absence of hydrodynamic forces, the method fails to give bounded response statistics, and the analysis yields stability conditions. It is the second-order stochastic averaging which can capture the influence of stiffness and inertia nonlinearities that were lost in the first-order averaging process. The results of the second-order averaging are compared with those predicted by Gaussian and non-Gaussian closures and by Monte Carlo simulation. In the absence of parametric excitation, the non-Gaussian closure solutions are in good agreement with Monte Carlo simulation. On the other hand, in the absence of hydrodynamic forces, second-order averaging gives more reliable results in the neighborhood of stochastic bifurcation. However, under pure parametric random excitation, the stochastic averaging and Monte Carlo simulation predict the on-off intermittency phenomenon near bifurcation point, in addition to stochastic bifurcation in probability.     

Nonlinear interactions in the planar dynamics of cable-stayed beam

An analytical model is proposed to study the nonlinear interactions between beam and cable dynamics in stayed-systems. The integro-differential problem, describing the in-plane motion of a simple cable-stayed beam, presents quadratic and cubic nonlinearities both in the cable equation and at the boundary conditions. Mainly studied are the effects of quadratic interactions, appearing at relatively low oscillation amplitude. To this end an analysis of the sensitivity of modal properties to parameter variations, in intervals of technical interest, has evidenced the occurrence of one-to-two and two-to-one internal resonances between global and local modes. The interactions between the resonant modes evidences two different sources of oscillation in cables, illustrated by simple 2dof discrete models.In the one-to-two global–local resonance, a novel mechanism is analyzed, by which cable undergoes large periodic and chaotic oscillations due to an energy transfer from the low-global to high-local frequencies.In two-to-one global–local resonance, the well-known parametric-induced cable oscillation in stayed-systems is correctly reinterpreted through the autoparametric resonance between a global and a local mode. Increasing the load the saturation of the global oscillations evidences the energy transfer from high-global to low-local frequencies, producing large cable oscillations. In both cases, the effects of detuning from internal and external resonance are presented.     

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.     

Stochastic non-linear flutter of a panel subjected to random in-plane forces

—An analysis of non-linear flutter of a simply-supported panel exposed to supersonic gas flow and random in-plane forces is presented for two- and three-mode interactions. A first order quasi-steady state aerodynamic piston theory is used to model the aerodynamic loading. The Fokker-Planck equation is used to derive a general moment equation for two- and three-mode interactions. For stability analysis the moment equation is consistent and the mean square stability boundaries of the equilibrium are obtained in terms of the system parameters. The stability boundaries reveal common features to those predicted by the deterministic theory of panel nutter. For the non-linear response the moment equation is found inconsistent and a cumulant-neglect closure is used by setting cumulants of fifth and sixth orders to zero. This first order non-Gaussian closure is carried out to solve for the response statistics in terms of the air-to-plate mass ratio, aerodynamic pressure, modal damping, and in-plane random force spectral density. It is found that the non-Gaussian solution yields higher levels for the response statistics than those obtained by the Gaussian solution. The inclusion of more modes results in a reduction of the response levels and expands the stability region.     

悬索在外激励作用下的1∶3内共振分析Ⅰ∶离散法

研究了悬索在受到外激励作用和考虑1∶3内共振情况下的两模态非线性响应。对于一定范围内的悬索弹性-几何参数而言,悬索第三阶面内对称模态的固有频率接近于第一阶面内对称模态的固有频率的3倍,从而导致1∶3内共振的存在。首先利用Galerkin方法把悬索的面内运动方程进行离散,然后利用多尺度法对离散的运动方程进行摄动,可得到两组不同主共振情况下的平均方程。     

斜拉桥拉索在轴向窄带随机激励下的振动响应

导出了拉索在考虑垂度以及索张力沿索长变化时的参激随机微分方程,进一步给出了预测拉索在窄带随机激励下响应的近似理论解------用统计矩截断法求解矩方程,获得高斯闭合解和一阶非高斯闭合解. 以南京长江二桥约330米长的A20拉索为研究对象,对以上高斯闭合解和一阶非高斯闭合解进一步进行数值求解以获得拉索的响应,并采用Monte-Carlo数值方法对求解进行验证. 分析了拉索振动的一般特征,特别分析了激励中心频率和拉索频率比为1和2时的响应随激励带宽的变化特征,得到了一些新的结论.      

On statistical quasi-linearization

In carrying out the statistical linearization procedure to a non-linear system subjected to an external random excitation, a Gaussian probability distribution is assumed for the system response. If the random excitation is non-Gaussian, however, the procedure may lead to a large error since the response of bother the original non-linear system and the replacement linear system are not Gaussian distributed. It is found that in some cases such a system can be transformed to one under parametric excitations of Gaussian white noises. Then the quasi-linearization procedure, proposed originally for non-linear systems under both external and parametric excitations of Gaussian white noises, can be applied to these cases. In the procedure, exact statistical moments of the replacing quasi-linear system are used to calculate the linearization parameters. Since the assumption of a Gaussian probability distribution is avoided, the accuracy of the approximation method is improved. The approach is applied to non-linear systems under two types of non-Gaussian excitations: randomized sinusoidal process and polynomials of a filtered process. Numerical examples are investigated, and the calculated results show that the proposed method has higher accuracy than the conventional linearization, as compared with the results obtained from Monte Carlo simulations.     

Application of an Itô-based approximation method to random vibration of a pinching hysteretic system

In this paper, an extension of the Cumulant-Neglect closure scheme is utilized for the random vibration analysis of a single degree of freedom system with a general pinching hysteresis restoring force. The hysteresis element used in the system model can simulate commonly observed forms of stiffness, strength and pinching degradations. The second order statistics of the system response to a stationary Gaussian white noise input are derived using an Itô-based approximation technique. The validity of these response statistics are then verified by comparing them to Monte Carlo simulation results. The numerical studies performed for different combinations of degradation parameters and excitation levels show that the response estimates obtained by this solution method are in good agreement with Monte Carlo simulation. These studies also indicate the applicability of this technique for response analysis of complicated forms of non-linearities.     

Dynamical modeling and non-planar coupled behavior of inclined CFRP cables under simultaneous internal and external resonances

A dynamic model for an inclined carbon ?ber reinforced polymer(CFRP)cable is established, and the linear and nonlinear dynamic behaviors are investigated in detail. The partial differential equations for both the in-plane and out-of-plane dynamics of the inclined CFRP cable are obtained by Hamilton's principle. The linear eigenvalues are explored theoretically. Then, the ordinary differential equations for analyzing the dynamic behaviors are obtained by the Galerkin integral and dimensionless treatments.The steady-state solutions of the nonlinear equations are obtained by the multiple scale method(MSM) and the Newton-Raphson method. The frequency-and force-response curves are used to investigate the dynamic behaviors of the inclined CFRP cable under simultaneous internal(between the lowest in-plane and out-of-plane modes) and external resonances, i.e., the primary resonances induced by the excitations of the in-plane mode,the out-of-plane mode, and both the in-plane mode and the out-of-plane mode, respectively. The effects of the key parameters, e.g., Young's modulus, the excitation amplitude,and the frequency on the dynamic behaviors, are discussed in detail. Some interesting phenomena and results are observed and concluded.     

Variations in steepness of the probability density function of beam random vibration

Dynamic behaviour of a beam, subjected to stationary random excitation, has been investigated for the situation in which the response is different from the model of a Gaussian random process. The study was restricted to the case of symmetric non-Gaussian probability density functions of beam vibrations. There are two possible causes of deviations of the system response from the Gaussian model: the first, nonlinear behaviour, concerns the system itself and the second is external when the excitation is not Gaussian. Both cases have been considered in the paper. To clarity the conclusions for each case and to avoid interference of these different types of system behaviour, two beam structures, clamped-clamped and cantilevered, have been studied. A numerical procedure for prediction of the nonlinear random response of a clamped-clamped beam under the Gaussian excitations was based on a linear modal expansion. Monte Carlo simulation was undertaken using Runge–Kutta integration of the generalised coordinate equations. Probability density functions of the beam response were analysed and approximated making use of different theoretical models. An experimental study has been carried out for a linear system of a cantilevered beam with a point mass at the free end. A pseudo-random driving signal was generated digitally in the form of a Fourier expansion and fed to a shaker input. To generate a non-Gaussian excitation a special procedure of harmonic phase adjustment was implemented instead of the random choice. In so doing, the non-Gaussian kurtosis parameter of the beam response was controlled.     

Internal-external resonance of a curved pipe conveying fluid resting on a nonlinear elastic foundation

In this study, the forced vibration of a curved pipe conveying fluid resting on a nonlinear elastic foundation is considered. The governing equations for the pipe system are formed with the consideration of viscoelastic material, nonlinearity of foundation, external excitation, and extensibility of centre line. Equations governing the in-plane vibration are solved first by the Galerkin method to obtain the static in-plane equilibrium configuration. The out-of-plane vibration is simplified into a constant coefficient gyroscopic system. Subsequently, the method of multiple scales (MMS) is developed to investigate external first and second primary resonances of the out-of-plane vibration in the presence of three-to-one internal resonance between the first two modes. Modulation equations are formed to obtain the steady state solutions. A parametric study is carried out for the first primary resonance. The effects of damping, nonlinear stiffness of the foundation, internal resonance detuning parameter, and the magnitude of the external excitation are investigated through frequency response curves and force response curves. The characteristics of the single mode response and the relationship between single and two mode steady state solutions are revealed for the second primary resonance. The stability analysis is carried out for these plots. Finally, the approximately analytical results are confirmed by the numerical integrations.     

Nonlinear random coupled motions of structural elements with quadratic nonlinearities

The second-order closure method is used to analyze the nonlinear response of two-degree-of-freedom systems with quadratic nonlinearities. The excitation is assumed to be the sum of a deterministic harmonic component and a random component. The case of primary resonance of the second mode in the presence of a two-to-one internal (autoparametric) resonance is investigated. The method of multiple scales is used to obtain four first-order ordinary-differential equations that describe the modulation of the amplitudes and phases of the two modes. Applying the second-order closure method to the modulation equations, we determine the stationary mean and mean-square responses. For the case of a narrow-band random excitation, the results show that the presence of the nonlinearity causes multi-valued mean-square responses. The multi-valuedness is responsible for a jump phenomenon. Contrary to the results of the linear analysis, the nonlinear analysis reveals that the directly excited second mode takes a small amount of the input energy (saturates) and spills over the rest of the input energy into the first mode, which is indirectly excited through the autoparametric resonance.     

随机激励下滞迟系统的稳态响应闭合解

滞迟系统属于一类典型的强非线性系统,滞迟力不仅取决于系统的瞬时变形,还与变形历程有关.虽然滞迟系统的随机振动问题已被广泛研究,但至今尚未得到滞迟系统随机响应概率密度函数的精确闭合解.本文运用迭代加权残值法获得了高斯白噪声激励下Bouc-Wen滞迟系统稳态响应概率密度函数的近似闭合解.首先,运用等效线性化法求出系统的稳态高斯概率密度函数;然后以此构造权函数,应用加权残值法求得了系统指数多项式形式的非高斯概率密度函数;最后引入迭代的过程,逐步优化权函数,提高计算所得结果的精度.以随机地震激励下钢纤维陶粒混凝土结构的稳态响应作为算例,其中Bouc-Wen模型的参数是基于拟静力学试验数据,并应用最小二乘法辨识获得.与Monte Carlo模拟结果相比,等效线性化法得到的结果精度较差;由加权残值法得到的结果能够表现出非线性特征,但其精度依然无法令人满意;采用迭代加权残值法得到的近似闭合解与Monte Carlo模拟的结果吻合非常好;对于较强随机激励情形,采用渐进迭代加权残值法具有较高的求解效率,所获得的理论解析解具有较高的精度.结果表明,所获得的近似闭合解不仅对于土木工程领域具有重要的实际应用价值,而且还可作为检验其他非线性系统随机响应预测方法的精度的标准.     

单自由度超弹性SMA减振结构随机振动控制理论研究

利用形状记忆合金（Shape Memory Alloy,简称SMA）丝的超弹性,提出了一种具有复位功能的阻尼器。在SMA丝的Graesser本构模型基础上,建立了阻尼器恢复力的滑移双线性模型;假定滞回面积相等,提出了恢复力的滑移刚塑性模型以近似简化滑移双线性模型。采用等价线性化法建立了单自由度超弹性SMA减振结构在高斯白噪声激励下的平稳随机振动分析公式。通过一算例,考虑不同激励谱密度和结构阻尼比：比较了等价线性法和蒙特卡罗（Monte Carlo）模拟法计算的结构振动响应（位移标准差和速度标准差）,证明了SMA减振结构随机振动控制理论的有效性;比较了等价线性减振结构和无控结构的动力特性（刚度和阻尼比）和振动响应,说明了SMA阻尼器能提高结构的刚度和阻尼比,因而可有效抑制结构的振动。     

Super-Harmonic and Internal Resonances of a Suspended Cable with Nearly Commensurable Natural Frequencies

In this paper, super-harmonic and internal resonance characteristics ofa viscously damped cable with nearly commensurable natural frequenciesare investigated by use of a novel method. The proposed frequency-domainsolution method is based on the combined use of a three-dimensionalnonlinear finite element approach and the incremental harmonic balancetechnique. It is an accurate algorithm in the sense that it accommodatesmulti-harmonic components and no mode-based model reduction is utilizedin the solution process. The alternating frequency/amplitude-controlledalgorithm enables complete solution to the frequency-response curvesincluding unstable branches, sub- and super-harmonic resonance andinternal resonance. A suspended cable paradigm under internal resonancecondition is studied using the proposed method. Nonlinear response andmodal interaction characteristics of the cable at different frequencyregions are identified from analysis of response profiles and harmoniccomponent features. The super-harmonic and internal resonance responsesare respectively characterized based on the harmonic distribution. Underan in-plane harmonic excitation, the two-to-one internal resonancebetween the in-plane and out-of-plane modes and the super-harmonicresonance around the second symmetric in-plane mode are revealed. Strongnonlinear interaction among different modes in the parameter spaceranging from primary resonance to super-harmonic resonance is observed.