质量矩阵和刚度矩阵怎么求频率

可以使用MATLAB求解这个问题

用Matlab 数值工具箱中库函数命令eig 来实现。

[ x ，d] = eig(ik ，im) % ik 和im 分别为结构的刚度矩阵和质量矩阵

d = sqrt (d)

for i = 1 :cn 　% cn 为结构的层数，即质点数

[ dl (i) ，j ] = min (d)

Xgd( : ，i) = x ( : ，j)

d(j) = max (d) + 1

end %以此循环对所求频率和振型进行排序

w = dl %所求自振频率

x = xgd 　%所求结构主振型

举个例子，可能会更清楚。

某三层钢筋混凝土结构，结构的各层特性参数为:第一层到第

三层质量m 分别为2 762 kg ，2 760 kg ，2 300 kg ，第一层到第三层

刚度k 分别为2. 485 ×104 N/ m ， 1. 921 ×104 N/ m ， 1. 522 ×

104 N/ m。地震波采用200gal El Centro 波，采样周期为0. 02 s。

单自由度体系自振频率(即共振频率)计算公式： ω=(k/m)^0.5

ω——振动圆频率

k——体系刚度

m——质点质量。