请问楼主是怎么解出振动频率的？是解关于s 的方程det(s^2*+s*+)=0求得的吗？
回复 16楼 hlhl 的帖子
对，具体的你可以看看8楼chaching老师给我回复的帖子 原帖由 jww100 于 2008-10-3 22:00 发表 http://www.chinavib.com/forum/images/common/back.gif谢谢院长指点，下图是我加入阻尼算出来的特征值
您的意思是这时的振动的频率取虚部就是
w1=0 w2=0 w3=0.5074 w4=-0.5074 w5=0.9977 w6=-0.9977 w7=0 w8=0
您看这样对吗？
不太对，主要是第一、第二阶固有频率有问题，咋出现负实根了呢，你改变一下参数试试看
也不是说一定不能有负实根，但是通常负实根不会出现 你利用方程det(s^2*+s*+)=0求得特征值？
这样做不太正规的，可能会出现问题
你最好还是计算特征值，给你一个方法，我不知道你学过想空间没有，把这个方程转换到相空间，然后用eig求特征值
回复 15楼 jww100 的帖子
不好意思, 我的语文水平很差有些词英文不清楚, 有些词中文不清楚, 总是用习惯的词
rigid body mode就是频率为0的模态, 即结构没有变形的模态
(strain energy应变能=0, 4个弹簧处皆没有变形)
还有可否直呼我名字即可, 我仅是小工程师, 怕害了老师名号
回复 18楼 咕噜噜 的帖子
谢谢院长，我已经用eig求出了固有频率，但只是无阻尼的情况。有阻尼的我还是没辙。
回复 20楼 ChaChing 的帖子
呵呵，老师这么谦虚，那我就恭敬不如从命了，以后我就呼您名字。回复 22楼 jww100 的帖子
试过10楼的提议吗? state space应该与小咕院长提到的空间相同运动方程Mx"+Cx'+Kx=f, 设q1=x, q2=x' => '=A*+B*f
其中A即为特徵矩阵 A=, 求eig(A)即可得
Ev=eig(A); fn = imag(Ev)/pi/2;
[ 本帖最后由 ChaChing 于 2008-10-6 10:03 编辑 ]
回复 23楼 ChaChing 的帖子
谢谢ChaChing,我赶紧试试去当时没试，呵呵，因为当时没看明白 原帖由 jww100 于 2008-10-6 09:21 发表 http://www.chinavib.com/forum/images/common/back.gif
谢谢院长，我已经用eig求出了固有频率，但只是无阻尼的情况。
有阻尼的我还是没辙。
无阻尼的你竟然算出来复数？:@L算错了
回复 25楼 咕噜噜 的帖子
被BS了:@L院长，复数的是我求带阻尼的 呵呵，还以为你说的是无阻尼的呢
不好意思，我理解错了，没有看不起你哟
慢慢来，做出来和大家分享一下
回复 27楼 咕噜噜 的帖子
谢谢院长鼓励，有这么多高人指导，加上咱的努力，咱一定能激流勇进，冲破万难啊，哈哈[ 本帖最后由 jww100 于 2008-10-6 22:49 编辑 ] 本帖最后由 甜炒栗子 于 2011-8-27 15:52 编辑
这类问题主要是eig函数的使用。
不过有阻尼的确实没遇到过。学习学习 楼主，这个问题你最后怎么解决的啊？我也遇到了，求指导~