含分数阶导数项随机Duffing振子稳态响应分析.docVIP

含分数阶导数项随机Duffing振子稳态响应分析 　　摘要： 对一个含有分数阶导数项阻尼的、Gaussian白噪声激励下的Duffing振子进行了稳态响应分析。首先，基于能量平衡理论，运用等效线性化方法，计算等效系统的线性阻尼及自然频率，建立统计意义下的等效线性化系统。然后，利用平均法建立随机It方程，得到随机响应的Markovian近似；给出描述振子振幅概率密度函数演化的FokkerPlanck方程，并得到它的稳态解。进一步，对于含有响应振幅的等效线性系统，借助由Laplace变换得到的转换函数，得到原系统的条件功率谱密度，结合振幅的稳态概率密度作为权重函数，给出原系统功率谱密度的估计，以及响应的统计量的估计。数值模拟的结果说明所提出的功率谱密度的近似解析表达式是可靠的，它甚至适用于Duffing振子具有强非线性回复力的情形，因为它可以较好地表现出功率谱密度共振频谱加宽及多峰现象的出现。 　　关键词： 分数阶导数； 等效线性化法； 随机平均法； 条件功率谱密度； 响应的功率谱密度估计 　　中图分类号： O322； O324 文献标志码： A 文章编号：2015 　　DOI：10.16385/j.cnki.issn2015.03.006 　　引 言 　　随机响应分析是随机动力学研究的重要方面，已有了一些经典的方法和结果[13]。对于响应的功率谱密度，文献[4]利用经典的等效线性化方法表现出缺陷，它给出的估计得到了正确的共振频率，却高估了共振处的峰值而且低估了频谱的宽度。一种改进的等效线性化方法在文献[5]中被提出，它与传统方法的主要区别在于将等效线性系统的待定参数设定为响应振幅的函数，从而建立了一个求解等效的线性刚度的计算方法。在文献[6]中，首次出现了术语“条件功率谱密度”，并通过与随机平均法的结合，以概率密度函数为权重函数，将响应的功率谱密度作为条件功率谱密度的加权和，给出了随机响应功率谱密度的有效估计。在文献[4]中被正式提出后，条件功率谱密度的概念开始被广泛地应用于响应功率谱密度的估计，给出了理想的近似解析结果，但它的数学严密性及合理性的说明却一直是一个空白，直到在文献[7]中，Spanos等才给出了严格的证明。 　　近年来，分数阶微积分在结构动力学及工程力学领域，越来越多的研究表明分数阶导数可以用较少的参数来模拟某些黏弹性材料的本构关系[810]，在确定性的响应分析中，基于特征向量展开，Laplace变换[89]和Fourier变换[11]被用来得到精确的解析解。文献[12]中，平均法被用来分析响应的振幅。在随机响应的情形，已有的方法主要适用于求解弱非线性系统。文献[13]中，基于广义谐和函数的随机平均法被用于对高斯白噪声激励下的含有分数阶导数的强非线性系统进行随机响应分析。此外，数值方法方面，已有的结果如有限差分方法、有限元方法及无网格方法等[1416]，都一定程度上受到分数阶导数全局依赖性这一本质特征对计算速度的影响。文献[17]给出的算法在一定程度上减低了分数阶导数项对历史数据的依赖和记忆性，提高了数值计算的速度。 　　本文旨在建立一个对含有分数阶导数项的随机Duffing振子进行随机响应分析的方法。利用改进的等效线性化方法，首先得到一个含分数阶导数项的线性系统，其中的线性阻尼及自然频率都是振幅的函数。将分数阶导数项作为一个弱阻尼项，对得到的系统运用平均法，建立随机It方程，得到响应的FokkerPlanck方程并求得其稳态解。然后，借助Laplace变换，得到等效线性系统的转换函数及条件功率谱密度，并综合之前的结果，对条件功率谱密度利用振幅的概率密度函数进行加权求和，得到原系统响应功率谱密度的估计。最后，通过数值模拟的结果来验证响应的稳态概率密度及功率谱密度轨迹的合理性。 　　6 结 论 　　对于含分数阶导数项的强非线性随机Duffing振子，本文给出了一个响应功率谱密度估计的程序。在假设分数阶导数项表现为结构阻尼的前提下，本文给出的功率谱密度估计程序是可靠的，即使是在较强的非线性系统情形。首先，利用改进的等效线性化法，得到一个与原系统等效的含分数阶导数项的线性随机Duffing振子，通过计算等价系统的线性阻尼和自然频率，原系统的非线性回复力被一个线性的回复力代替，而且它是振幅响应的函数。然后，随机平均法被作用于得到的线性振子，通过系统振幅响应的Markovian近似，建立了振幅响应的随机It方程；写出描述振幅概率密度函数演化的FokkerPlanck方程，并给出其稳态解。最后，利用等效线性系统的转换函数，得到随机响应的条件功率谱密度，并结合得到的振幅稳态概率密度函数，得到响应功率谱密度的估计。数值结果表明，文中所给出的近似解析结