求助：有阻尼时MATLAB中求固有频率以及振型

欧阳中华 · 发表于 2012-2-25 20:42

.
你阻尼矩阵不对称估计是不是MATLAB里的求特征函数解不了...

kukumalu180614 · 发表于 2012-2-26 09:20

回复 15 # dingxinran 的帖子

vbr sf(m,d,k)不对，少了个下划线，vbr_sf(m,d,k)就行了。

dingxinran · 发表于 2012-2-29 15:12

kukumalu180614 发表于 2012-2-26 09:20

登录/注册后可看大图

回复 15 # dingxinran 的帖子

vbr sf(m,d,k)不对，少了个下划线，vbr_sf(m,d,k)就行了。

ok 谢谢了哈哈

agu19851116 · 发表于 2012-3-1 22:00

固有频率和振型时是不考虑阻尼的
用[vec,val]=eig(inv(M)*K);
计算出特征向量阵vec和特征值val，再排序即可得到振型矩阵和频率

ME! · 发表于 2013-1-8 19:29

agu19851116 发表于 2012-3-1 22:00

登录/注册后可看大图

固有频率和振型时是不考虑阻尼的
用[vec,val]=eig(inv(M)*K);
计算出特征向量阵vec和特征值val，再排序 ...

得到特征向量后怎么画振型图啊

伤痕累累 · 发表于 2013-1-10 18:15

ME! 发表于 2013-1-8 19:29
得到特征向量后怎么画振型图啊

特征值对应特征向量
V为特征向量的话
n为其中一列的行数,i为某阶固有频率对应的阶数
plot(1:n,V(i，1:n))%第一阶的振型图

ME! · 发表于 2013-1-20 14:03

本帖最后由 ME! 于 2013-1-20 14:18 编辑

伤痕累累发表于 2013-1-10 18:15

登录/注册后可看大图

特征值对应特征向量
V为特征向量的话
n为其中一列的行数,i为某阶固有频率对应的阶数

不对吧，特征向量里面既有x,y,还有转角的位移参数，你这样画图，横坐标是1：自由度数，纵坐标为某阶特征向量中所有的值，而且我也试了，画的图不对
x坐标应为模型的坐标，即单元节点处的坐标，我想问的是y坐标是振型向量中y位移的单元节点处的值吗 forj=1:3:48; plot(x(j),V(j,48));

伤痕累累 · 发表于 2013-1-20 18:33

本帖最后由伤痕累累于 2013-1-20 18:48 编辑

ME! 发表于 2013-1-20 14:03

登录/注册后可看大图

不对吧，特征向量里面既有x,y,还有转角的位移参数，你这样画图，横坐标是1：自由度数，纵坐标为某阶特征 ...

对应的。。。
根据你选取的对应来。东西是死的，人是活的。
我所说的只是普遍的，通用的，就如楼主所说的那种质量刚度模型，是没有考虑转角的。见郑兆昌《机械振动》上册。
你要是得XY的，就该只提取XY 的特征向量，如果不知道变通，那肯定是不对的了。
再给你看看我的例子，
我的里面也有转角部分，出质量的振型图时已经将转角那部分去掉啦。

[V,D]=eig(K,M);
n=length(K);
W=zeros(n,1);
for i=1:n
W(i)=sqrt(D(i,i));%固有频率
end
plot(1:2:n,V(1:2:n,4),1:2:n,V(1:2:n,3),'k')%前两阶振型
grid on
未命名1.jpg

振型为一阶弯曲，二阶锥型。和理论上的结论完全一致！

ME! · 发表于 2013-1-21 15:18

伤痕累累发表于 2013-1-20 18:33

登录/注册后可看大图

对应的。。。
根据你选取的对应来。东西是死的，人是活的。
我所说的只是普遍的，通用的，就如楼主所 ...

你这不是第三阶和第四阶振型吗？怎么是1,2阶，还有你哼坐标的表示方式和你的模型是对不上的，只能说你的是画出大概的样子，不过你的意思我理解了！

伤痕累累 · 发表于 2013-1-21 17:08

ME! 发表于 2013-1-21 15:18

登录/注册后可看大图

你这不是第三阶和第四阶振型吗？怎么是1,2阶，还有你哼坐标的表示方式和你的模型是对不上的，只能说你的是 ...

你还是好好多看点书吧。横坐标我只是没有处理，偶数表示的是转角，我给那部分去掉了。这只是从程序里随意做的一个例子，我做这个只是用来验证模型是否正确，不是去发表论文的，所以没有做那么标准。还有3,4并不表示3,4,阶。只是对应的特征值的特征向量。如果你再无法理解，我实在没有办法了，表达能力欠缺。这都是最最基础的东西了。

ME! · 发表于 2013-1-22 14:20

伤痕累累发表于 2013-1-21 17:08

登录/注册后可看大图

你还是好好多看点书吧。横坐标我只是没有处理，偶数表示的是转角，我给那部分去掉了。这只是从程序里随意 ...

意思是你的第一阶固有频率出现在第三个特征值上？不过我自己的振型图，我已经画出来了！

ME! · 发表于 2013-3-6 16:32

伤痕累累发表于 2013-1-21 17:08

登录/注册后可看大图

你还是好好多看点书吧。横坐标我只是没有处理，偶数表示的是转角，我给那部分去掉了。这只是从程序里随意 ...

我又出现问题了，http://forum.chinavib.com/forum. ... p;page=1&extra=

402144999 · 发表于 2014-7-26 14:01

提示: 作者被禁止或删除内容自动屏蔽

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

402144999 402144999 当前离线积分 22 头像被屏蔽	发表于 2014-7-26 14:01 \| 显示全部楼层提示: 作者被禁止或删除内容自动屏蔽
402144999 402144999 当前离线积分 22 头像被屏蔽
	回复支持反对使用道具举报显身卡

[转子动力学] 求助：有阻尼时MATLAB中求固有频率以及振型

点评

评分

点评