你阻尼矩阵不对称估计是不是MATLAB里的求特征函数解不了... 回复 15 # dingxinran 的帖子
vbr sf(m,d,k)不对，少了个下划线，vbr_sf(m,d,k)就行了。 kukumalu180614 发表于 2012-2-26 09:20 static/image/common/back.gif
回复 15 # dingxinran 的帖子
vbr sf(m,d,k)不对，少了个下划线，vbr_sf(m,d,k)就行了。
ok 谢谢了 哈哈 固有频率和振型时是不考虑阻尼的
用=eig(inv(M)*K);
计算出特征向量阵vec和特征值val，再排序即可得到振型矩阵和频率
agu19851116 发表于 2012-3-1 22:00 static/image/common/back.gif
固有频率和振型时是不考虑阻尼的
用=eig(inv(M)*K);
计算出特征向量阵vec和特征值val，再排序 ...
得到特征向量后怎么画振型图啊 ME! 发表于 2013-1-8 19:29
得到特征向量后怎么画振型图啊
特征值对应特征向量
V为特征向量的话
n为其中一列的行数,i为某阶固有频率对应的阶数
plot(1:n,V(i，1:n))%第一阶的振型图
本帖最后由 ME! 于 2013-1-20 14:18 编辑
伤痕累累 发表于 2013-1-10 18:15 http://www.chinavib.com/static/image/common/back.gif
特征值对应特征向量
V为特征向量的话
n为其中一列的行数,i为某阶固有频率对应的阶数
不对吧，特征向量里面既有x,y,还有转角的位移参数，你这样画图，横坐标是1：自由度数，纵坐标为某阶特征向量中所有的值，而且我也试了，画的图不对
x坐标应为模型的坐标，即单元节点处的坐标，我想问的是y坐标是振型向量中y位移的单元节点处的值吗 forj=1:3:48;plot(x(j),V(j,48));
本帖最后由 伤痕累累 于 2013-1-20 18:48 编辑
ME! 发表于 2013-1-20 14:03 http://forum.chinavib.com/static/image/common/back.gif
不对吧，特征向量里面既有x,y,还有转角的位移参数，你这样画图，横坐标是1：自由度数，纵坐标为某阶特征 ...
对应的。。。
根据你选取的对应来。东西是死的，人是活的。
我所说的只是普遍的，通用的，就如楼主所说的那种质量刚度模型，是没有考虑转角的。见郑兆昌《机械振动》上册。
你要是得XY的，就该只提取XY 的特征向量，如果不知道变通，那肯定是不对的了。
再给你看看我的例子，
我的里面也有转角部分，出质量的振型图时已经将转角那部分去掉啦。
=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
振型为一阶弯曲，二阶锥型。和理论上的结论完全一致！
伤痕累累 发表于 2013-1-20 18:33 static/image/common/back.gif
对应的。。。
根据你选取的对应来。东西是死的，人是活的。
我所说的只是普遍的，通用的，就如楼主所 ...
你这不是第三阶和第四阶振型吗？怎么是1,2阶，还有你哼坐标的表示方式和你的模型是对不上的，只能说你的是画出大概的样子，不过你的意思我理解了！ ME! 发表于 2013-1-21 15:18 static/image/common/back.gif
你这不是第三阶和第四阶振型吗？怎么是1,2阶，还有你哼坐标的表示方式和你的模型是对不上的，只能说你的是 ...
你还是好好多看点书吧。横坐标我只是没有处理，偶数表示的是转角，我给那部分去掉了。这只是从程序里随意做的一个例子，我做这个只是用来验证模型是否正确，不是去发表论文的，所以没有做那么标准。还有3,4并不表示3,4,阶。只是对应的特征值的特征向量。如果你再无法理解，我实在没有办法了，表达能力欠缺。这都是最最基础的东西了。 伤痕累累 发表于 2013-1-21 17:08 static/image/common/back.gif
你还是好好多看点书吧。横坐标我只是没有处理，偶数表示的是转角，我给那部分去掉了。这只是从程序里随意 ...
意思是你的第一阶固有频率出现在第三个特征值上？不过我自己的振型图，我已经画出来了！ 伤痕累累 发表于 2013-1-21 17:08 static/image/common/back.gif
你还是好好多看点书吧。横坐标我只是没有处理，偶数表示的是转角，我给那部分去掉了。这只是从程序里随意 ...
我又出现问题了，http://forum.chinavib.com/forum.php?mod=viewthread&tid=124213&page=1&extra=
页:
1
[2]