yms_1 发表于 2009-10-12 15:14

【求助】复特征值求解时的一个诡异现象。

自编有限元，计算结构的质量矩阵【M】和刚度矩阵【K】，然后求解（（ω*ω）【M】+【K】）{x}=0的广义特征值。固有频率计算结果与ANSYS的结果对比如下：
阶数我的结果ANSYS的结果145.65945.648264.99364.987385.19685.1924147.47147.475157.95157.946178.51178.517185.94185.948267.24267.249345.43345.4310392.27392.2711424.21424.2112445.31445.3113455.22455.2214597.36597.3615724.19724.1916786.68786.6817888.21888.2118941.48941.4819970.67970.67201027.31027.3

将有限元控制方程改为（（ω*ω）【M】+（1+0.5i）【K】）{x}=0   ，计算复特征值。结果得到如下内容：
阶数实部虚部1 1.5589E+35-1.5589e+035i2 1.5886E+27-2.9465e+026i3 1.0918E+35+1.0918e+035i4 1.7619E+27+1.8351e+027i5 2.2644E+26+1.3394e+027i6 1.8012E+27-1.471e+027i7 2.7135E+27+2.6093e+027i8 9.9111E+26-1.0358e+027i9 8.4299E+26+6.3189e+026i10 1.2330E+27-1.4786e+027i11 8.8330E+261.263e+027i12 1.3528E+27+1.1022e+027i13 6.2785E+26-5.3406e+026i14 2.9609E+26-4.9818e+026i15 2.5363E+26-4.5582e+026i16 4.3492E+26+3.3144e+026i17 1.0270E+27+8.0443e+026i18 9.7542E+26+8.5914e+026i19 6.9768E+26-4.0533e+026i20 1.9850E+268.312e+026i


结果根本不对，这是为什么呢？
我分别采用了fortran的IMSL库中的DGVLCG函数和Matlab中的eig函数，结果都差不多这样。请大家指点！！

[ 本帖最后由 yms_1 于 2009-10-12 15:16 编辑 ]

ChaChing 发表于 2009-10-13 07:51

Matlab中的eig函数? LZ怎么用!?
方便给齐资料与程式?

yms_1 发表于 2009-10-13 09:54

回复 沙发 ChaChing 的帖子

谢谢您的关注。哈哈

yms_1 发表于 2009-10-13 09:56

由于矩阵数据较大，必须压缩。要看矩阵数据的话，请将work.dat改为work.rar

yms_1 发表于 2009-10-14 09:07

版主们都跑哪里去了啊，在线等，不信你不出来！:'(
页: [1]
查看完整版本: 【求助】复特征值求解时的一个诡异现象。