非亏损系统二次特征值问题的Rt阶法iz降.pdfVIP

振 动 与 冲 击 第 2l卷第 4期 JOURNAL OFVIBRATION AND SHOCK 非亏损系统二次特征值 问题 的 Ritz降阶法 邸 娜 陈力奋 王 皓 (复旦大学力学与工程科学系 ，上海 200433) 摘 要 本文提出了一种利用Ritz降阶对含有重特征值及刚体模态的非亏损阻尼系统进行模态分析的方法。本 方法利用保守系统的左 、右截断模态作为 Ritz基，综合 出非保守系统二次特征值问题 的待求左 、右特征向量。文章证 明了 在大多数情况下，本方法可用较少 的耗费，提取满足给定精度的待求二次特征值问题的模态解。文末 的算例表明了本方 法的有效性 。 关键词：非亏损系统，二次特征值 问题 ，Ritz降阶 中图分类号：0327 0 引 言 设 Mo，C0，Ko，∈RN 分别表示』、r维非亏损阻尼 系统的质量、阻尼和刚度矩阵 ，则系统 的动力学方程 』、r自由度非亏损阻尼系统的动力学方程在数学 为 +c0+Koq：．厂(t)，若外激励 ．厂(t)∈ 可表 上转化为求解一个二次特征值 问题。常用的解法是把 示为广义加速度 、广义速度 口和广义坐标 口的线性 它投影到状态空间，求解一组一阶的2N维方程。但对 函数．厂(t)：M1；+C1+ 1g，1，c1，1，∈ Ⅳ， 于高阶系统来说 ，这不仅意味着大量 的计算耗费，而 此时 ，动力学方程可改写为 ： 且特性矩阵的性态难 以保证 ，导致计算 的条件数差。 nq+ + ：0 (1) 尤其 ，工程上一般只对低阶模态感兴趣 ，因此就更是 其 中，M ：Mo—M1，C：C0一C1，K：Ko—K1，分另0 得不偿失。近几十年来 ，寻求一种迅速而有效地求解 是系统的广义质量、广义阻尼和广义刚度矩阵，一般 二次特征值 问题的方法，已成为结构动力学和计算方 不具有任何对称性 。这里假定 为可逆矩阵。其相 法领域共 同感兴趣 的课题。AdhikariS j利用伽辽金 应的二次特征值问题为： 误差最小方法和 Nemnann展开 ，提出了一种基于保守 ( + C+ ) R：0 (20) 系统的完全模态计算相应非保守系统特征向量 的方 ( Mr+ C +Kr)L：0 (2b) 法。然而 ，该方法实质上是一种坐标变换法 ，当 』、r很 这里 ∈C， ∈CN(t：L，R)分别为特征值和特 大时 ，这是不切实 际的做法 ，并没有减少计算量。此 征向量 (左、右)。通过求解下列特征方程可得系统的 外 ，Adhikari法只考虑了保守系统与非保守系统仅有 2N个特征值 (i=1，2，…，2N)： 简单