简支梁非线性响应的三维路径积分法分析

采用基于Gauss-Legendre积分公式的三维路径积分法,分析了在过滤高斯白噪声激励下的简支梁非线性随机振动响应的概率密度函数;联立一阶滤波方程与简支梁一阶模态的振动模型,得到在过滤高斯白噪声激励下的简支梁随机振动模型,基于Gauss-Legendre积分公式的积分法和短时高斯近似法求解响应的概率密度函数值。结果表明,三维路径积分法计算值与蒙特卡洛模拟值符合良好,即使在尾部区域也符合良好。三维路径积分法比等效线性化法的计算精度更高。     

三维路径积分法的积分区间以及子区间数目选取方法研究

对基于Gauss-Legendre公式的路径积分法在二维、三维模型上的运用进行了研究。以高斯白噪声激励下的附加集中质量块悬臂梁、单稳态、双稳态、三稳态能量收集器为例,采用路径积分法计算平稳或非平稳响应的概率密度函数解,重点研究了积分区间以及子区间数目的选取对于路径积分法计算精度和效率的影响,采用蒙特卡洛数值模拟法的计算结果作为标准值,与路径积分法计算的结果进行了对比。提出了先对积分区间以及子区间数目进行粗选,再根据粗选结果进行负反馈调节积分区间的方法,提高了路径积分法的计算精度。     

基于路径积分法的悬索非线性随机振动响应分析

随机荷载激励下悬索过大的动力响应将影响其正常使用与安全,对其响应概率密度函数的求解与分析是评估悬索随机动力响应的重要途径之一。针对悬索在高斯白噪声激励下的随机振动模态响应,利用基于Gauss-Legendre积分和短时高斯转移概率密度假定的路径积分法,研究了模态振动响应的概率密度函数的平稳数值解与非平稳数值解,并进一步开展了参数研究,揭示了不同参数影响下概率密度函数的分布规律。将路径积分法所得的平稳解和非平稳解,分别与FPK方程的精确平稳解、等效线性化法所得平稳解及蒙特卡罗模拟非平稳解进行对比,结果表明,路径积分法所得的概率密度函数解分别与精确平稳解及蒙特卡罗模拟非平稳解符合良好。对于平稳响应,由于位移二次非线性项的存在,位移概率密度函数分布呈非对称分布形式,但速度概率密度函数并不受其影响,仍服从对称分布;非平稳响应概率密度函数初始时刻峰值较大,且在初始阶段峰值是随着时间不断变化的,波动较明显,随着时间推移逐渐平稳。研究结果对于悬索非平稳响应研究具有重要的工程意义。     

空间相关过滤白噪声激励下结构的随机振动离散分析方法

以随机振动的离散分析方法为基础,讨论了在空间相关过滤白噪声激励下结构的随机振动分析,给出了计算结构均值和均方响应的具体公式,并对一空间索网结构进行了风振响应计算。     

加性非平稳激励下结构速度响应的FPK方程降维

结构在随机激励下的非线性响应分析是具有高度挑战性的困难问题. 对于白噪声或过滤白噪声激励,求解FPK方程将获得结构响应 的精确解. 遗憾的是,对于非线性多自由度系统,FPK方程难以直接求解. 事实上,其数值解法严重受限于方程维度,而解析求解 则仅适用于少数特定的系统,且多是稳态解. 因此,将FPK方程进行降维,是求解高维随机动力响应分析问题的重要途径. 本文针 对幅值调制的加性白噪声激励下多自由度非线性结构的非平稳随机响应分析问题,将联合概率密度函数满足的高维FPK方程进行降 维. 针对结构速度响应概率密度函数求解,通过引入等价漂移系数,原FPK方程可转化为一维FPK型方程. 建议了构造等价漂移系数 的条件均值函数方法. 进而,采用路径积分方法求解降维FPK型方程,得到速度概率密度函数的数值解答. 结合单自由度Rayleigh 振子、十层线性剪切型框架和非线性剪切型框架结构在幅值调制的加性白噪声激励下的非平稳速度响应求解,讨论了本文方法的精 度和效率,验证了其有效性.      

加性非平稳激励下结构速度响应的FPK方程降维

结构在随机激励下的非线性响应分析是具有高度挑战性的困难问题.对于白噪声或过滤白噪声激励,求解FPK方程将获得结构响应的精确解.遗憾的是,对于非线性多自由度系统,FPK方程难以直接求解.事实上,其数值解法严重受限于方程维度,而解析求解则仅适用于少数特定的系统,且多是稳态解.因此,将FPK方程进行降维,是求解高维随机动力响应分析问题的重要途径.本文针对幅值调制的加性白噪声激励下多自由度非线性结构的非平稳随机响应分析问题,将联合概率密度函数满足的高维FPK方程进行降维.针对结构速度响应概率密度函数求解,通过引入等价漂移系数,原FPK方程可转化为一维FPK型方程.建议了构造等价漂移系数的条件均值函数方法.进而,采用路径积分方法求解降维FPK型方程,得到速度概率密度函数的数值解答.结合单自由度Rayleigh振子、十层线性剪切型框架和非线性剪切型框架结构在幅值调制的加性白噪声激励下的非平稳速度响应求解,讨论了本文方法的精度和效率,验证了其有效性.     

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

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

剪切模量是幂函数的成层土的地震随机反应

基于改进的一维剪切梁模型，对剪切模量是其深度的某一幂函数的成层非均质土层，得到其稳态动力响应的封闭型解析表达式。首次证明了这种土层振型函数的正交性，然后利用随机振动理论，并基于基岩输入地震加速度的功率谱密度函数：白噪声谱和过滤白噪声谱。研究了该土层对地震的随机动力响应问题。计算结果表明，1)在基岩输入地震加速度的功率谱为白噪声谱的情况下，土层的最大期望反应均有别于过滤白噪声谱时的相应值；2)平稳输入与输出过高地估计了土层的随机响应。     

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

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

基于state-space-split法的滞回非线性系统随机动力响应分析

采用SSS(state-space-split)法,建立了引入Bouc-Wen滞回模型的杜芬非线性系统在高斯白噪声激励下的概率密度函数(PDF)的近似求解方法,分析了其随机动力响应变化规律.首先,将Bouc-Wen滞回模型引入杜芬非线性系统,分别考虑非线性系统中的几何非线性和材料非线性对动力响应的影响.随后,对该模型进...     

Path integral solution of vibratory energy harvesting systems

A transition Fokker-Planck-Kolmogorov(FPK) equation describes the procedure of the probability density evolution whereby the dynamic response and reliability evaluation of mechanical systems could be carried out. The transition FPK equation of vibratory energy harvesting systems is a four-dimensional nonlinear partial differential equation. Therefore, it is often very challenging to obtain an exact probability density. This paper aims to investigate the stochastic response of vibration energy harvesters(VEHs)under the Gaussian white noise excitation. The numerical path integration method is applied to different types of nonlinear VEHs. The probability density function(PDF)from the transition FPK equation of energy harvesting systems is calculated using the path integration method. The path integration process is introduced by using the GaussLegendre integration scheme, and the short-time transition PDF is formulated with the short-time Gaussian approximation. The stationary probability densities of the transition FPK equation for vibratory energy harvesters are determined. The procedure is applied to three different types of nonlinear VEHs under Gaussian white excitations. The approximately numerical outcomes are qualitatively and quantitatively supported by the Monte Carlo simulation(MCS).     

Non-stationary response of a van der Pol-Duffing oscillator under Gaussian white noise

This paper investigates the probability density evolution process of a van der Pol-Duffing oscillator under Gaussian white noise. A path integration method is employed with the Gauss–Legendre integration scheme. In the path integration method, the short-time Gaussian approximation scheme is used for computing the one-step transition probability density. Two cases are considered with slight nonlinearity or strong nonlinearity in displacement. The stationary and non-stationary responses of the oscillator are studied. Compared with the simulation result, the path integration method can present a satisfactory probability density function (PDF) solution for each case. Different probability density evolution processes are observed correspondingly. In the case of slight nonlinearity, the PDF undergoes a clockwise motion around the origin. The peak region gradually expands and the PDF eventually forms a circle. By contrast, the strong nonlinearity drives the oscillator to oscillate around the limit cycle. In such a case, the PDF rapidly forms a circle. The circle keeps its shape and develops until the oscillator becomes stationary. More complicated phenomena can be studied by the adopted path integration method.     

Numerical path integration of a non-homogeneous Markov process

The numerical path integration method, based on Gauss-Legendre integration scheme, is applied to a Duffing oscillator subject to both sinusoidal and white noise excitations. The response of the system is a Markov process with one of the drift coefficients being periodic. It is a non-homogeneous Markov process that does not have a stationary probability distribution. When applying the numerical procedure, the values of transition probability density at the Gaussian-Legendre quadrature points need only be calculated for time steps of the first period of the sinusoidal excitation, and they can be saved for use in all subsequent periods. The numerical procedure is capable of capturing the evolution of the probability density from an initial distribution to one that is changing and rotating periodically in the phase space.     

泊松白噪声激励下强非线性系统的半解析瞬态解

自然界与工程中都普遍存在着随机扰动,且大多数呈现出固有的非高斯性质,若采用高斯激励建模可能会导致巨大的误差.泊松白噪声作为一种典型且重要的非高斯激励模型,已引起了广泛的关注.目前,泊松白噪声激励下系统的动态特性分析主要集中于稳态响应的研究,而针对瞬态响应的求解难度仍较大,需进一步发展.本文引入径向基神经网络,提出了一种泊松白噪声激励下单自由度强非线性系统瞬态响应预测的高效半解析方法.首先将广义Fokker-Plank-Kolmogorov (FPK)方程的瞬态解表示为一组含时变待定权值系数的高斯径向基神经网络;然后采用有限差分法离散时间导数项,并结合随机取样技术构造含时间递推式的损失函数;最后通过拉格朗日乘子法使得损失函数最小化获得时变最优权值系数.作为算例,探究了两个经典强非线性系统,并采用蒙特卡罗模拟方法对解析结果加以验证.结果表明:本文方法所获得的瞬时概率密度函数与蒙特卡罗模拟数据吻合地较好,并且算法具备较高的计算效率.在系统响应的整个演化过程中,本文所提方法能够非常有效地捕捉到系统响应在各个时刻下的复杂非线性特征.此外,本文方法所获得的高精度半解析瞬态解,不仅可作为基准解检验其...     

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

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

Probabilistic solution of non-linear random ship roll motion by path integration

This paper investigates the probability density function (PDF) of non-linear random ship roll motion using a previously developed path integration method. The mathematical model of ship rolling motion consists of a linear-plus-cubic damping and a non-linear restoring moment in the form of odd-order polynomials up to fifth-order terms. In the path integration method, the interpolation scheme is based on the Gauss–Legendre quadrature integration rule and the short-time transition probability density function is formulated by short-time Gaussian approximation. The present work extends the path integration method to the case of non-linear random ship roll motion. Different values of non-linearity coefficient and excitation intensity are used to examine the effectiveness of the path integration method. Numerical analysis shows that the results of the path integration method agree well with the simulation results, even in the tail region. The path integration method is effective and it is simply implemented in the examined cases. Due to the presence of non-linear damping terms and non-linear restoring moment terms, the PDFs of roll angle and angular velocity exhibit highly non-Gaussian behaviors.     

Vibration reduction evaluation of a linear system with a nonlinear energy sink under a harmonic and random excitation

The nonlinear behaviors and vibration reduction of a linear system with a nonlinear energy sink(NES)are investigated.The linear system is excited by a harmonic and random base excitation,consisting of a mass block,a linear spring,and a linear viscous damper.The NES is composed of a mass block,a linear viscous damper,and a spring with ideal cubic nonlinear stiffness.Based on the generalized harmonic function method,the steady-state Fokker-Planck-Kolmogorov equation is presented to reveal the response of the system.The path integral method based on the Gauss-Legendre polynomial is used to achieve the numerical solutions.The performance of vibration reduction is evaluated by the displacement and velocity transition probability densities,the transmissibility transition probability density,and the percentage of the energy absorption transition probability density of the linear oscillator.The sensitivity of the parameters is analyzed for varying the nonlinear stiffness coefficient and the damper ratio.The investigation illustrates that a linear system with NES can also realize great vibration reduction under harmonic and random base excitations and random bifurcation may appear under different parameters,which will affect the stability of the system.     

分数导数型本构关系描述粘弹性梁的振动分析

本文研究粘弹性梁在周期激励作用下的受迫振动问题。梁的材料满足Kelvin-Volgt分数导数型本构关系。基于动力学方程、本构关系和应变－位移关系建立了小变形粘弹性梁的振动方程。采用分离变量法分析粘弹性梁的自由振动,导出模态坐标满足的常微分－积分方程和模态函数满足的常微分方程,对于两端简支的截面梁给出了固有频率和模态函数。对于简谐激励作用下粘弹性梁的受迫振动,利用模态叠加得到了稳态响应。最后给出数值算例说明本文方法的应用。     

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.