okok12 发表于 2007-8-12 11:51

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

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

pengweicai 发表于 2007-8-12 12:43

ABLE方法
QR算法

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

欧阳中华 发表于 2007-8-12 13:44

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

appleseed05 发表于 2007-8-12 18:38

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

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

欧阳中华 发表于 2007-8-12 20:30

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

okok12 发表于 2007-8-12 20:56

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

欧阳中华 发表于 2007-8-13 08:49

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

zhaowu 发表于 2007-8-13 19:41

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

pengweicai 发表于 2007-8-13 21:41

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

xiaoyao111 发表于 2007-8-27 21:34

QR不行么,
可以吧

billions1943 发表于 2015-6-30 21:26

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

z1z11z111 发表于 2016-2-3 19:44

用数值求特征值没有问题,但无法用模特分析。因为模态分析的整个原来可是基于模态正交的基础上的,现在非对称,没有正交性。
页: [1]
查看完整版本: 刚度、阻尼矩阵非对称时的特征值问题