[求助]请教大虾用MATLAB计算发动机固有频率和主振型的问题

bulobuloo · 发表于 2006-6-2 12:37

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

最近毕业设计要答辩了但是matlab计算这一块的问题还没有解决急上眉梢了都 有没有哪位前辈搞过类似的课题帮我写个程序 程序主要是通过已知的质量矩阵和刚度矩阵来求固有频率和主振型的 跪求了急急急！！

[此贴子已经被cdwxg于2006-6-3 0:41:22编辑过]

lxq · 发表于 2006-6-2 12:54

具体没有什么公式吗? 没学过这方面的不懂咯

wangzi629 · 发表于 2006-6-2 13:09

[ZX,W2]=eig(inv(M)*K); %M质量矩阵（对角阵），K刚度矩阵（对称阵） W=sqrt(W2) %固有频率 ZX %振型

yangzj · 发表于 2006-6-2 13:12

看看振动的书,大概是求inv(M)*K的特征值和特征向量就行了

wangzi629 · 发表于 2006-6-2 13:13

W阵主对角线元素即为固有频率。对应的ZX阵列向量即为该频率下的振型。当然你也可以按频率大小排序，但要注意对应关系。

bulobuloo · 发表于 2006-6-2 13:18

楼上的能把控制画出固有频率和主振形图形的那部分程序给具体写下么
你说的这个程序是计算这个的

恩我看过曾令贤的一篇论文也是用的这个程序导师说这样做的精度不够但是没办法了时间紧太紧！！
呵呵谢谢大家回帖这个论坛热心的前辈真不少
还请有能力的各位把后面排序和PLOT曲线的那部分帮下忙嘿嘿惭愧了

[ 本帖最后由 ChaChing 于 2010-3-31 00:55 编辑 ]

bulobuloo · 发表于 2006-6-2 13:51

M=[2000 0 0 0 0 0 0 2000 0 0 0 0 0 0 2000 0 0 0 0 0 0 7056000 0 0 0 0 0 0 2587000 0 0 0 0 0 0 8428000]; K=[235200 0 0 0 0 0 0 165800 0 0 0 -1028000 0 0 58200 0 3609000 0 0 0 0 370700000 0 0 0 0 3609000 0 479900000 0 0 -10280000 0 0 0 1405000000];
p=[]; [ZX,W2]=eig(inv(M)*K); W=sqrt(W2); n=size(inv(M)*K); for i=1:n p(i)=W(i,i); ZX(:,i)=ZX(:,i)/ZX(1,i); end for j=1:n-1 for i=j+1:n if p(j)>p(i) t=p(i); p(i)=p(j); p(j)=t; q=ZX(:,i); ZX(:,i)=ZX(:,j); ZX(:,j)=q; end end end for i=1:n subplot(6,1,i) plot(ZX(:,i)) ylabel('主振型') xlabel('x') grid on end
这个是我按照提示编写的程序我在软件上运行时画出的图特征不明显不知道是哪里出了错误请前辈门帮看看哪有错误！1

wangzi629 · 发表于 2006-6-6 17:59

你的K阵写错了吧，它不是对称阵，2行6列和6行2列的数差一个数量级。

当然，还可以用现成的函数sort进行排序，这里用起泡法进行的排序。不知你的x轴单位是什么？这个程序你运行下吧，看看符合你的要求不？

[ 本帖最后由 ChaChing 于 2010-3-31 00:58 编辑 ]

wangzi629 · 发表于 2006-6-6 18:43

clc clear M=[2000 0 0 0 0 0 0 2000 0 0 0 0 0 0 2000 0 0 0 0 0 0 7056000 0 0 0 0 0 0 2587000 0 0 0 0 0 0 8428000]; K=[235200 0 0 0 0 0 0 165800 0 0 0 -1028000 0 0 58200 0 3609000 0 0 0 0 370700000 0 0 0 0 3609000 0 479900000 0 0 -1028000 0 0 0 1405000000];
[ZX,W2]=eig(inv(M)*K); W=sqrt(W2); n0=size(inv(M)*K); n=n0(2);
for i=1:n p(i)=W(i,i); ZX_max=max(ZX(:,i)); ZX(:,i)=ZX(:,i)/ZX_max; end
for j=1:n-1 for i=1:n-j if p(i)>p(i+1) tmp=p(i); p(i)=p(i+1); p(i+1)=tmp; ZX_tmp=ZX(:,i); ZX(:,i)=ZX(:,i+1); ZX(:,i+1)=ZX_tmp; end end end
for i=1:n subplot(6,1,i) plot(ZX(:,i)) ylabel_1=['第' num2str(i) '阶主振型']; ylabel(ylabel_1) xlabel('x') grid on end

[此贴子已经被作者于2006-6-6 18:48:04编辑过]

bulobuloo · 发表于 2006-6-6 22:44

谢谢您花时间帮我解决这个问题在你的程序中把p(i)=W(i,i);句改为p(i)=10*W(i,i)/(2*pi);所得到的固有频率与书上答案一致不知道是为什么不管怎么样还是十分感谢你的帮助有机会还希望得到您的指导和帮助！~我QQ52976010 希望我门能成为朋友

wangzi629 · 发表于 2006-6-7 09:16

p(i)=10*W(i,i)/(2*pi)这是什么意思？W(i,i)就是固有频率啊。还有你的K阵，M阵数值看一下，和书上一致不？

michael111 · 发表于 2006-6-7 22:12

刚度矩阵是不是有问题？？求自振频率要处理边界条件，在质量矩阵和刚度矩阵中划去相应的行列，再用eig命令求解不会有错的！！！！！！

realyyy · 发表于 2006-6-17 17:19

<DIV class=quote>以下是引用bulobuloo在2006-6-6 22:44:52的发言： 
谢谢您花时间帮我解决这个问题在你的程序中把p(i)=W(i,i);句改为p(i)=10*W(i,i)/(2*pi);所得到的固有频率与书上答案一致不知道是为什么不管怎么样还是十分感谢你的帮助有机会还希望得到您的指导和帮助！~我QQ52976010 希望我门能成为朋友 </DIV>
 不会吧？！W是在这句命令之前算出来的啊！p(i)对W没有影响的。 我试过了：把p(i)=W(i,i);改为p(i)=10*W(i,i)/(2*pi);对结果没有任何影响！（不影响W也不影响排序！）

xuebx · 发表于 2006-6-18 08:32

感谢wangzi629坛友的精彩回答

多情清秋 · 发表于 2006-6-19 08:16

wangzi629加威望3点，bulobuloo加威望2点，realyyy、michael111、yangzj加威望1点
多情清秋 2006年6月19日

账号		自动登录	找回密码
密码			我要加入

[综合讨论] [求助]请教大虾用MATLAB计算发动机固有频率和主振型的问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（bulobuloo）谢谢您花时间帮我解决这个问题在...

回复：（bulobuloo）[求助]请教大虾用MATLAB计算发动机...

[综合讨论] [求助]请教大虾用MATLAB计算发动机固有频率和主振型的问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（bulobuloo）谢谢您花时间帮我解决这个问题 在...

回复：（bulobuloo）[求助]请教大虾用MATLAB计算发动机...

回复：（bulobuloo）谢谢您花时间帮我解决这个问题在...