|

楼主 |
发表于 2012-4-18 10:47
|
显示全部楼层
两位坛友,还有问题需要请教。我现在做了一个柱的实例,节点数为1430个,自由度为4290。现在提取了他的整体质量矩阵和整体刚度矩阵进行计算,采用matlab中[EigVector,EigValue]=eigs(k,m,n,'SM') 这个函数求其前几阶的特征值,但是结果与有限元模态分析结果相差很大。有限元结果1阶只有0.5Hz,但是计算出来的确有将近10Hz,而且各阶分布较近,不知道这是什么原因,将我模型的各系统矩阵贴上,请坛友指点。另外就是在提取模型的各系统矩阵时要不要对模型施加约束。
有限元计算前10阶频率为[0.527,0.53,3.21,3.22,6.04,8.60,8.64,9.69,15.95,16.01];
通过矩阵计算结果为[9.44,18.14,18.33,18.53,18.73,18.96,19.07,19.26,19.31,19.52];
其中sitffness为整理的前100阶整体刚度矩阵,mass为前100阶质量矩阵。
|
|