刚度、阻尼矩阵非对称时的特征值问题
各位,我最近在做一个系统运动方程的刚度、阻尼矩阵非对称时的 模态参数识别问题,遇到如下问题:在结构本身的刚度、阻尼矩阵都是非对称矩阵的时候,怎么求解其特征值,谢谢@ ABLE方法
QR算法
Ritz向量 等
建议到 期刊或者学位论文网搜搜看。 .
系统出现非对称矩阵时,往往简单地反应为非对称矩阵的特征值如何求得,不过建议求非对称矩阵系统特征值同时,考虑一下,非对称出现的物理背景是什么?计算出的非对称矩阵系统的特征值表示的是什么? 这个在计算方法上我作过一些比较,对于非对称的特征值问题或者非正定矩阵的特征值问题,我采用自己编写的计算方法(可能是自己只是采用常规的公式或者用的不太明白)比如子空间迭代,QR方法等等,计算结果都没有matlab自带的QZ方法计算结果好。看到过一些文献说对于这种问题,采用Lanczos或者ritz这些方法结果会比较好一些,这个没自己写过程序,但是用过被人写的程序算过,结果比QR和子空间迭代要好。不过我个人觉得采用matlab的QZ方法计算结果还是不错的:lol :lol
[ 本帖最后由 appleseed05 于 2007-8-12 18:40 编辑 ] .
matlab方法仅仅只能对计算规模比较小的问题,而且计算速度较慢、占用内存也相对大得多.. .. to 欧阳,我做风工程,最明显的就是气动导数改变了系统的刚度矩阵以及阻尼矩阵,而且气动导数对这些矩阵的改变是非对称的。 .
涉及流体时,往往存在升力,由此导致“负阻尼”效应,物理上还要考虑是不是会出现颤振现象.. .. 用matlab中eig命令可能会容易些。我对比过,精度也可以。
另外,还有子空间迭代、Lanczos法。 楼主可以根据情况而定,如果问题不是太大,精度要求也不是特高。 为了简单,高效率起见,可以采用MATLAB计算。 要自己编一个计算特征值的程序要费不少功夫。 QR不行么,
可以吧 我用的mathcad,命令不太一样··· 用数值求特征值没有问题,但无法用模特分析。因为模态分析的整个原来可是基于模态正交的基础上的,现在非对称,没有正交性。
页:
[1]