非经典系统的模态降阶与动柔度阵的新算法.pdfVIP

行在海面上船只的运动脚，或者飞行器的振动f6】。不对称的阻尼和冈Ⅱ度阵经常出现在陀螺 问题及含有随动力的系统中，以及有的系统质量阵也是不对称的(soom和硒mp】)。研究 人员已经着手解决非经典线性系统的二次特征值问题。FawzS，和Bishop州提出了非经典系统 中的左、右特征向量所满足的一些关系式．进而，给出了一个规范化左、右特征向量的方法。 Inman p1研究了一类特殊的非经典系统，这类系统可由线性变换转化为等价的系数矩阵为对 称的经典系统，保持特征值不变。后来Ahraadian和hm∞…给出了这类可对称化的非经典 系统存在正交规范模态的充要条件。之后，Caughey和Mapll给出了更一般的非经典非保守 系统可以解耦的充要条件。在他们接下来一篇文章中，Ma Caughey[12】应用了等价变换方 法求解这类可以解耦的非经典系统。然而，他们的工作局限于保守系统或者特殊的可以通过 等价变换解耦的非经典非保守系统。 对于一般的非经典系统，则通常把问题投影到状态空间求解113-l“，从而把=次特征值 问题降阶为一次特征值问题。但却把系数矩阵的尺寸变成了原来的两倍。对于高阶系统来说， 这不仅意味着大量的计算耗费，而且特性矩阵的性态难以保证，比如矩阵原有的稀疏性被破 坏了。 在工程计算中，也发展出一些有效的近似算法。应用摄动方法p”1】可解决只有轻微阻 尼的二次特征值问题。 因为是小阻尼系统，所以这种方法可近似采用相应的保守系统作为 未受扰动的系统。但是摄动法局限于小阻尼系统，而且为了求解特征值或特征向量的任意 lj}传2)阶近似项，需要知道前面所有征一1)阶的近似项，并且它们都是相关的。另～种 算法都局限于经典系统。矩阵迭代法适用于计算系统的最低几阶固有频率和模态；子空间迭 代法是矩阵迭代法与Ritz法的结合，是直接针对二次特征值问题的算法．而不需要把问题 转化到状态空间；I．anczos算法是基于K．,ylov序列的一种递归算法，最早是针对保守系统提 8l，在商业软件中得到应用B9一“。但 出的，80年代后期被推广到解决二次特征值问题pl-3 是无论摄动法还是迭代法都有一个缺点，就是对于一个普通的结构系统需要付出大量的计算 耗费。 2 系统的完全模态集计算相应非保守系统特征向量的方法。这种方法可以解决含有一般的阻尼 矩阵的系统，而且不需要化到状态空间求解。然而．该方法实质上是一种坐标变换法，当自 由度很大时，这是不切实际的做法，并没有减少计算量。此外，A【mikari法只考虑了保守系 统与非保守系统仅有成对出现的共轭复根的简单特征值的情形，而忽略了工程实际中经常碰 L．F．ch％[42I等将Ad．bikari法扩展 到的重特征值情形，也没讨论零特征值情形。H，W．Song 到保守系统有重特征值的情况，并对特征值为实数的情况作了讨论，但是依然基于保守系统 的完全模态集。没有减少计算量，而且对于非保守系统含有重特征值的情形也未加讨论。 本文的第二章将提出一种利用Ritz降阶对含有重特征值及零特征值的非经典系统进行 模态分析的方法。本方法利用保守系统的左、右截断模态作为Ritz基，综合出非保守系统 二次特征值问题的待求左、右特征向量。这样既不需要把问题转化到状态空间求解，又通过 Ritz降阶进一步减少自由度数，因此在大多数情况下，本方法可用较少的耗费，提取满足给 定精度的待求二次特征值问题的模态解。另外，本方法的通用性覆盖了含有重特征值及零特 征值的系统。 另一方面，工程振动问题中，动刚度阵联系着系统的响应和激励。在求响应时。需要计 算动刚度阵的逆阵一一动柔度阵。非经典系统的动刚度阵一般为二阶九阵： D2轨)=矿肘+kC+五。其中』I，，C，置∈R””分别为非经典系统的质量、阻尼和刚度阵， 可以既非对称又非反称。本文中研究的非经典系统动刚度阵就是数学上的九矩阵，阶数不限， 涵概了一般的二阶(一阶)动刚度阵。由于九矩阵D伉)的各元素是特征值九的多项式，直 接求解其逆阵——动柔度阵并非易事，促使一些学者做了这方面的研究。LcIllIgp3】早期仅研 究了保守系统动柔度阵的算法。胡海昌阻】针对非保守系统提出了预解式法，利用部分分式 展开得到了动柔度阵的表达式，但要求质量阵对称正定，同时假定系统不存在重特征值。 Pa