matlab求固有频率,eig()函数求固有频率和有限元软件求固有频率结果不相等

本帖最后由 爱模型啊 于 2017-12-26 16:54 编辑

我现在在用Abaqus有限元软件做模态分析，已经求出了各阶的模态频率，也知道怎么输出系统的质量矩阵、系统的刚度矩阵。

在有限元软件中，对一个2D平面正方形做模态分析，求各阶的模态频率，边界条件是对正方形下边全部约束自由度，划分了4个网格，如下所示：

201712261.PNG (8.36 KB, 下载次数: 31)

2017-12-26 16:17 上传

然后我想在matlab里面，通过eig()函数对系统的质量矩阵、系统的刚度矩阵计算来求出模态频率，来和Abaqus有限元软件的结果进行对比，这是我的代码：

clc

clear

K13=[1.00000000000000e+36        36785714285.7143        -36835233516.4835        -2829670329.67033        0        0        36736195054.9451        2829670329.67033        -76351579670.3297        -36785714285.7143        0        0        0        0        0        0        0        036785714285.7143        1.00000000000000e+36        2829670329.67033        36736195054.9451        0        0        -2829670329.67033        -36835233516.4835        -36785714285.7143        -76351579670.3297        0        0        0        0        0        0        0        0

-36835233516.4835        2829670329.67033        1.00000000000000e+36        0        -36835233516.4835        -2829670329.67033        -76351579670.3297        36785714285.7143        73472390109.8901        0        -76351579670.3297        -36785714285.7143        0        0        0        0        0        0

-2829670329.67033        36736195054.9451        0        1.00000000000000e+36        2829670329.67033        36736195054.9451        36785714285.7143        -76351579670.3297        0        -73670467032.9670        -36785714285.7143        -76351579670.3297        0        0        0        0        0        0

0        0        -36835233516.4835        2829670329.67033        1.00000000000000e+36        -36785714285.7143        0        0        -76351579670.3297        36785714285.7143        36736195054.9451        -2829670329.67033        0        0        0        0        0        0

0        0        -2829670329.67033        36736195054.