2D梁的模态分析

以下面的问题为例，讨论FEM进行梁的模态分析：

单跨梁受约束如图所示，梁的密度ρ=3e3kg/m3，试求梁的自振频率及振型。设梁的弹性模量为E=200GPa,横截面积A=1.0e-2㎡，长度L=2，界面惯性矩为I=1E-4m⁴。

然后主程序fem.m

主程序调用下面两个代码：

生成单元刚度矩阵的程序Beam2D2Node_Stima.m：

生成单元质量矩阵的程序:

对于上面的程序，给出了前五阶固有频率：

a=sqrt(diag(V));

a(1:5)

1.0e+04 *

    0.0718

    0.4498

    1.2594

    2.4679

    4.0796

并且给出了其对应的振型：

挠度的振型：

subplot(5,1,1), plot(coordinates,[V(1:2:end,1);0])
subplot(5,1,2), plot(coordinates,[V(1:2:end,2);0])
subplot(5,1,3), plot(coordinates,[V(1:2:end,3);0])
subplot(5,1,4), plot(coordinates,[V(1:2:end,4);0])
subplot(5,1,5), plot(coordinates,[V(1:2:end,5);0])

相关资源：连续梁的弯矩计算软件V1.0绿色版_连续梁-其它代码类资源-CSDN文库