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

yms_1

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

阶数	我的结果	ANSYS的结果
1	45.659	45.648
2	64.993	64.987
3	85.196	85.192
4	147.47	147.47
5	157.95	157.94
6	178.51	178.51
7	185.94	185.94
8	267.24	267.24
9	345.43	345.43
10	392.27	392.27
11	424.21	424.21
12	445.31	445.31
13	455.22	455.22
14	597.36	597.36
15	724.19	724.19
16	786.68	786.68
17	888.21	888.21
18	941.48	941.48
19	970.67	970.67
20	1027.3	1027.3

将有限元控制方程改为（（ω*ω）【M】+（1+0.5i）【K】）{x}=0   ，计算复特征值。结果得到如下内容：

阶数	实部	虚部
1	1.5589E+35	-1.5589e+035i
2	1.5886E+27	-2.9465e+026i
3	1.0918E+35	+1.0918e+035i
4	1.7619E+27	+1.8351e+027i
5	2.2644E+26	+1.3394e+027i
6	1.8012E+27	-1.471e+027i
7	2.7135E+27	+2.6093e+027i
8	9.9111E+26	-1.0358e+027i
9	8.4299E+26	+6.3189e+026i
10	1.2330E+27	-1.4786e+027i
11	8.8330E+26	1.263e+027i
12	1.3528E+27	+1.1022e+027i
13	6.2785E+26	-5.3406e+026i
14	2.9609E+26	-4.9818e+026i
15	2.5363E+26	-4.5582e+026i
16	4.3492E+26	+3.3144e+026i
17	1.0270E+27	+8.0443e+026i
18	9.7542E+26	+8.5914e+026i
19	6.9768E+26	-4.0533e+026i
20	1.9850E+26	8.312e+026i

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

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

ChaChing

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

yms_1

谢谢您的关注。哈哈

yms_1

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

yms_1

版主们都跑哪里去了啊，在线等，不信你不出来！:'(