刚度、阻尼矩阵非对称时的特征值问题

okok12

各位，我最近在做一个系统运动方程的刚度、阻尼矩阵非对称时的 模态参数识别问题，遇到如下问题：
   在结构本身的刚度、阻尼矩阵都是非对称矩阵的时候，怎么求解其特征值，谢谢@

pengweicai

ABLE方法
QR算法

Ritz向量 等
建议到 期刊或者学位论文网  搜搜看。

欧阳中华

.
   系统出现非对称矩阵时，往往简单地反应为非对称矩阵的特征值如何求得，不过建议求非对称矩阵系统特征值同时，考虑一下，非对称出现的物理背景是什么?计算出的非对称矩阵系统的特征值表示的是什么？

appleseed05

这个在计算方法上我作过一些比较，对于非对称的特征值问题或者非正定矩阵的特征值问题，我采用自己编写的计算方法（可能是自己只是采用常规的公式或者用的不太明白）比如子空间迭代，QR方法等等，计算结果都没有matlab自带的QZ方法计算结果好。看到过一些文献说对于这种问题，采用Lanczos或者ritz这些方法结果会比较好一些，这个没自己写过程序，但是用过被人写的程序算过，结果比QR和子空间迭代要好。不过我个人觉得采用matlab的QZ方法计算结果还是不错的:lol :lol

[ 本帖最后由 appleseed05 于 2007-8-12 18:40 编辑 ]

欧阳中华

.
   matlab方法仅仅只能对计算规模比较小的问题，而且计算速度较慢、占用内存也相对大得多.. ..

okok12

to 欧阳，我做风工程，最明显的就是气动导数改变了系统的刚度矩阵以及阻尼矩阵，而且气动导数对这些矩阵的改变是非对称的。

欧阳中华

.
   涉及流体时，往往存在升力，由此导致“负阻尼”效应，物理上还要考虑是不是会出现颤振现象.. ..

zhaowu

用matlab中eig命令可能会容易些。我对比过，精度也可以。
另外，还有子空间迭代、Lanczos法。

pengweicai

楼主可以根据情况而定，如果问题不是太大，精度要求也不是特高。 为了简单，高效率起见，可以采用MATLAB计算。 要自己编一个计算特征值的程序要费不少功夫。

xiaoyao111

QR不行么，
可以吧

billions1943

我用的mathcad  ，命令不太一样···

z1z11z111

用数值求特征值没有问题，但无法用模特分析。因为模态分析的整个原来可是基于模态正交的基础上的，现在非对称，没有正交性。