cqucivil 发表于 2014-3-28 11:05

matlab已知刚度、质量、阻尼矩阵怎样求频率?

已知一个空间框架的刚度矩阵、质量矩阵和阻尼矩阵 ,能否由matlab求出其频率? 求大神指点,有没有相关程序

cqucivil 发表于 2014-3-28 12:47

{:{39}:}

歪猴 发表于 2014-3-28 12:55

多自由度的求解公式清楚吧,,,把他们编程一下,,就可以了   

cqucivil 发表于 2014-3-28 14:13

歪猴 发表于 2014-3-28 12:55
多自由度的求解公式清楚吧,,,把他们编程一下,,就可以了

谢谢那个清楚但是根据公式要求响应啊能否具体指点一二

yyxt007 发表于 2014-3-29 17:42

eig函数。

歪猴 发表于 2014-4-19 10:31

cqucivil 发表于 2014-3-28 14:13
谢谢那个清楚但是根据公式要求响应啊能否具体指点一二

同济大学 ji小熊版的 《 汽车振动学》 有相关程序   你你可以查一下

suffer 发表于 2014-4-19 21:21

当有阻尼时(M——质量矩阵;C—— 阻尼矩阵;K——刚度矩阵)
M*X''+CX'+KX=0

有一个近似算法是现在常用的,就是先求出无阻尼系统的固有频率及相应振型矩阵,用命令eig(K,M,'qz')或者eig(K,M),将系统正则化,求出正则阻尼矩阵,其可能不是对角矩阵,但取其对角元素即可,那么阻尼系统固有频率为wr=w*sqrt(1-c),其中c为正则阻尼矩阵对角线上元素,而系统的振型不变,因此原来地振型矩阵中各列向量仍对应阻尼系统各固有频率

function =vbr_sf(m,d,k)
%vbr_sf   vbr_sf(m,d,k)
%      =vbr4(m,d,k)
%      function vbr_sf finds the mode shapes and natural frequencies of
%      a linear second order matrix equation.
%      =vbr_sf(m,k) finds the mode shapes and natural frequencies
%      for the undamped case.
if nargin==2
   k=d;
   =eig(m\k);
   w=sqrt(w);
end
if nargin==3
   if norm(d/m*k-k/m*d) < 1e-8*norm(k/m*d)
%disp('Damping is proportional, eigenvectors are real.')
=eig(m\k);
w=sqrt(w);
zeta=(v'*m*v)\(v'*d*v)/2/w;
else
%disp('Damping is non-proportional, eigenvectors are complex.')
a=;
=eig(a);
w=abs(w1);
zeta=-real(w1)/w;
   end
end
w=diag(w);zeta=diag(zeta);

来自:http://blog.sina.com.cn/s/blog_9f81fb09010157c1.html
页: [1]
查看完整版本: matlab已知刚度、质量、阻尼矩阵怎样求频率?