计算的固有频率为负值？

king · 发表于 2006-9-19 18:50

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

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

x

我在用Matlab计算一个振动系统的固有频率时，得出的结果有一个是负值（特征值矩阵对角元素有一个为负值），程序检查过多遍了，不知道是否能忽略负号的影响，请高手指点。

[ 本帖最后由 lxq 于 2006-9-19 21:52 编辑 ]

xinyuxf · 发表于 2006-9-19 18:54

计算的固有频率里不应该出现负值吧，这应该还是程序的问题，是不是装配矩阵时有错误啊？建议再查查程序

AaronSpark · 发表于 2006-9-20 00:28

把你的程序贴一下，大家帮你看看

king · 发表于 2006-9-22 08:43

都搞了好几天了，最后发现是书上的公式错了，现在终于解决了。
哎，现在的书籍真是不能相信，作者都太不负责任了。

king · 发表于 2006-9-22 08:47

还是非常感谢大家的提示

ME! · 发表于 2012-11-23 21:47

我计算出来的固有频率有复数w=eig(inv(M)*K)
wn=sqrt(w)
f=wn/2/pi
w有负数，所以wn有复数，如果是因为刚度矩阵不对称导致，那么要忽略那些小于0和为复数的值吗，是不是取固有频率可以直接从非负的实数开始

xuning0418 · 发表于 2012-11-26 19:00

请问一下，matlab如何画板的振动频率。。。

happy · 发表于 2012-11-27 09:52

ME! 发表于 2012-11-23 21:47

登录/注册后可看大图

我计算出来的固有频率有复数w=eig(inv(M)*K)
wn=sqrt(w)
f=wn/2/pi

建议专门学一下复模态理论
复模态的固有频率为复特征的模

happy · 发表于 2012-11-27 09:53

xuning0418 发表于 2012-11-26 19:00

登录/注册后可看大图

请问一下，matlab如何画板的振动频率。。。

振动频率往往都是独立的数值形式出现的，你想问的是模态怎么画吧？

ME! · 发表于 2012-11-27 10:30

请问8楼“happy“楼主，将计算出来的固有频率取模就可以了吗？，我总觉得w不可能为负数，还有什么其它原因造成这个结果吗

happy · 发表于 2012-11-27 10:50

ME! 发表于 2012-11-27 10:30

登录/注册后可看大图

请问8楼“happy“楼主，将计算出来的固有频率取模就可以了吗？，我总觉得w不可能为负数，还有什么其它原因造 ...

估计是你的质量矩阵和刚度矩阵有问题造成
复模态一般出现在非正定阻尼的情况下

ME! · 发表于 2012-11-27 18:54

本帖最后由 ME! 于 2012-11-27 19:18 编辑

E=2.11e11;
rho=7800;
A=pi*0.1^2;
I=pi*0.2^4/64;
L=0.01;
W1=(A*L*rho)/3 ;
W2=(A*L*rho)/6;
W3=(156*A*L*rho)/420;
W4=(22*A*L^2*rho)/420;
W5=(54*A*L*rho)/420 ;
W6=(13*A*L^2*rho)/420;
W7=(3*A*L^3*rho)/420;
W8=(4*A*L^3*rho)/420;
M=[W1  0 0 W2 0 0;
0 W3 W4  0 W5 -W6;
0 W4 W8 0 W6 -W7
W2  0 0 W1 0 0;
0 W5 W6 0 W3 -W4
0 -W6 -W7 0 -W4 W8];
R1=A*E/L;
R2=12*I*E/(L^3);
R3=(6*E*I)/(L*L);
R4=4*I*E/L;
R5=2*I*E/L;
K=[R1,0,0,-R1,0,0;
0,R2,R3,0,-R2,R3;
0,R3,R4,0,-R3,R5;
-R1,0,0,R1,0,0;
0,-R2,-R3,0,R2,-R3;
0,R3,R5,0,-R3,R4];
D=K/M;
[ZX W9]=eig(D);
W=eig(D);
wn=sqrt(W)
求出来的W总有负值

ChaChing · 发表于 2012-11-27 21:25

ME! 发表于 2012-11-27 18:54

登录/注册后可看大图

E=2.11e11;
rho=7800;
A=pi*0.1^2;

M是对称矩阵, 也是正定矩阵
K是对称矩阵, 但并非是正定矩阵
Ref :http://forum.chinavib.com/thread-108027-1-1.html

xuning0418 · 发表于 2012-11-27 22:06

happy 发表于 2012-11-27 09:53

登录/注册后可看大图

振动频率往往都是独立的数值形式出现的，你想问的是模态怎么画吧？

嗯，是模态，，问一个小问题哈，板的振动解一般常用哪种形式，有带双曲函数形式的吗？谢谢，谢谢！！！

happy · 发表于 2012-11-28 14:32

xuning0418 发表于 2012-11-27 22:06

登录/注册后可看大图

嗯，是模态，，问一个小问题哈，板的振动解一般常用哪种形式，有带双曲函数形式的吗？谢谢，谢谢！！！

具体解的形式和你板的结构以及约束方式有关
一些简单的板模态表达式可以在一些板壳理论的书中找到

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

[综合讨论] 计算的固有频率为负值？

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

麻烦帮我看一下