求一个求四自由度振动系统固有频率的MATLAB程序

hlhl · 发表于 2008-10-5 14:16

请问楼主是怎么解出振动频率的？是解关于s 的方程det(s^2*[M]+s*[C]+[K])=0求得的吗？

jww100 · 发表于 2008-10-5 23:58

对，具体的你可以看看8楼chaching老师给我回复的帖子

咕噜噜 · 发表于 2008-10-6 08:38

原帖由 jww100 于 2008-10-3 22:00 发表

 登录/注册后可看大图

谢谢院长指点，下图是我加入阻尼算出来的特征值
您的意思是这时的振动的频率取虚部就是
w1=0 w2=0 w3=0.5074 w4=-0.5074 w5=0.9977 w6=-0.9977 w7=0 w8=0
您看这样对吗？

不太对，主要是第一、第二阶固有频率有问题，咋出现负实根了呢，你改变一下参数试试看
也不是说一定不能有负实根，但是通常负实根不会出现

咕噜噜 · 发表于 2008-10-6 08:40

你利用方程det(s^2*[M]+s*[C]+[K])=0求得特征值？
这样做不太正规的，可能会出现问题
你最好还是计算特征值，给你一个方法，我不知道你学过想空间没有，把这个方程转换到相空间，然后用eig求特征值

ChaChing · 发表于 2008-10-6 08:43

不好意思, 我的语文水平很差
有些词英文不清楚, 有些词中文不清楚, 总是用习惯的词
rigid body mode就是频率为0的模态, 即结构没有变形的模态
(strain energy应变能=0, 4个弹簧处皆没有变形)
还有可否直呼我名字即可, 我仅是小工程师, 怕害了老师名号

jww100 · 发表于 2008-10-6 09:21

谢谢院长，我已经用eig求出了固有频率，但只是无阻尼的情况。
有阻尼的我还是没辙。

jww100 · 发表于 2008-10-6 09:29

呵呵，老师这么谦虚，那我就恭敬不如从命了，以后我就呼您名字。

ChaChing · 发表于 2008-10-6 09:56

试过10楼的提议吗? state space应该与小咕院长提到的空间相同
运动方程Mx"+Cx'+Kx=f, 设q1=x, q2=x' => [q1; q2]'=A*[q1;q2]+B*f
其中A即为特徵矩阵 A=[zeros(n,n) ones(n,n); -inv(M)*K, -inv(M)*C], 求eig(A)即可得
Ev=eig(A); fn = imag(Ev)/pi/2;

[ 本帖最后由 ChaChing 于 2008-10-6 10:03 编辑 ]

jww100 · 发表于 2008-10-6 10:12

谢谢ChaChing,我赶紧试试去
当时没试，呵呵，因为当时没看明白

咕噜噜 · 发表于 2008-10-6 11:36

原帖由 jww100 于 2008-10-6 09:21 发表

 登录/注册后可看大图

谢谢院长，我已经用eig求出了固有频率，但只是无阻尼的情况。
有阻尼的我还是没辙。

无阻尼的你竟然算出来复数？:@L算错了

jww100 · 发表于 2008-10-6 14:36

被BS了:@L
院长，复数的是我求带阻尼的

咕噜噜 · 发表于 2008-10-6 16:30

呵呵，还以为你说的是无阻尼的呢
不好意思，我理解错了，没有看不起你哟
慢慢来，做出来和大家分享一下

jww100 · 发表于 2008-10-6 22:36

谢谢院长鼓励，有这么多高人指导，加上咱的努力，咱一定能激流勇进，冲破万难啊，哈哈

[ 本帖最后由 jww100 于 2008-10-6 22:49 编辑 ]

甜炒栗子 · 发表于 2011-8-27 15:52

本帖最后由甜炒栗子于 2011-8-27 15:52 编辑

这类问题主要是eig函数的使用。
不过有阻尼的确实没遇到过。学习学习

李佳期 · 发表于 2013-11-25 08:25

楼主，这个问题你最后怎么解决的啊？我也遇到了，求指导~

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

[综合讨论] 求一个求四自由度振动系统固有频率的MATLAB程序

回复 16楼 hlhl 的帖子

回复 15楼 jww100 的帖子

回复 18楼咕噜噜的帖子

回复 20楼 ChaChing 的帖子

回复 22楼 jww100 的帖子

回复 23楼 ChaChing 的帖子

回复 25楼咕噜噜的帖子

回复 27楼咕噜噜的帖子

点评