非对称非线性系统非平稳随机响应中心差分法.pdfVIP

第22卷第 2期 振 动 工 程 学 报 Vo1．22No．2 2009年 4月 JournalofVibrationEngineering Apr．2009 非对称非线性系统非平稳随机响应的中心差分法 张 明，聂 宏 (南京航空航天大学飞行器先进设计技术国防重点学科实验室，江苏 南京 210016) 摘要：利用随机中心差分法离散化格式，结合时变等价线性化步骤，建立了多 自由度非线性系统受非 白噪声随机激 励作用下非平稳响应协方差矩阵的递推算法。在考虑到非对称项引起的非零均值响应条件下，对等价线性化方法 进行修正，使之适用于非对称系统 。算例表明：该方法计算结果和MonteCarlo模拟结果符合较好 ，且对所给算例 ， 此方法 比原中心差分法精确。该方法计算步骤简单，便于在计算机上实现，计算效率高。 关键词：随机振动；非线性；非对称；中心差分法；等价线性化 中图分类号 ：0322；0324 文献标识码 ：A 文章编号：1004—4523(2009)02—0128—06 引 言 1 随机动力学模型的建立 在非线性系统 随机响应 的预测方法中，等价线 考虑以下 自由度非线性系统在随机激励下的 性化方法仍然是最简单可行 ，也是工程 中应用最为 振动 问题 广泛的一种方法[1]。采用等价线性化方法 ，可与时域 MX+ CX+ KX+ g( ，)= ，(f) (1) 直接积分法 、模态分析法 、虚拟激励法、双频谱法等 式 中 M，C，K分别为系统的质量、阻尼和刚度矩 结合[2]，求解非线性系统受非平稳随机激励的问 阵；f()为零均值平稳或非平稳高斯随机过程 ， 题。其中，时域直接积分法 中的中心差分法由于考虑 g(X，)为非线性的刚度和阻尼 ，是对称的，即为奇 了等效参数 的时变性，且计算步骤简单，便于在计算 函数，满足g(X，)一一g(，)，是关于 和 函数 机上实现 ，得到了比较好的发展[8叫 。但是 ，对于刚 构成 的列阵 度和阻尼同时含非线性的多 自由度系统，尤其是激 g(X，)一 [g (，)，g (，)，…，g (，x)3 励为非白噪声激励 的情形 ，已有 的文献并没有建立 (2) 响应协方差矩阵计算 的一般完整递推算法 。 由等价线性化理论 ，当假定响应为零均值 的高 上述采用等价线性化方法的文献 中非线性项大 斯过程时，相应 的等价线性化系统为 多为对称的，即为奇函数 [3]。对于非对称系统 ，响 MX+ C + K， 一 ，(f) (3) 应的均值非零 ，等价线性化方法的基本假设不成立 ， 式中 系数C 和 的元素为 直接应用会带来误差。文献[2]虽然计及了非对称项 引起 的响应均值变化 ，但是并未考虑其对等效参数 一 E[ ]+c， ㈤ 及二阶矩的影响。 K===E[ ]+K (5) Spanos和朱位秋提 出了处理非对称非线性和 E(·)表示取期望。显然，C嘶和 嘶是响应 和 非零均值响应的等价线性化方法[1卜¨]，将非零均值 矩的函数，对于非平稳系统 ， “和 K 均是时变 响应从随机微分方程中分离 ，使其响应均值为零 ，然 的。考虑等效系数的时变性，将时间t按间隔 离散 后采用等价线性化方法求解 。本文将