sunjinyuan 发表于 2009-7-1 15:05

求助：关于动力总成的特征值求解问题

我编了些语句，计算动力总成的6阶特征值，但是算下来就是不对啊，现在基本确定K和M矩阵没问题，我忽略了阻尼，所以写的语句如下：
format long
M=Mmotor(207.5,5.0892,15.8066,14.2851,-0.1505,-0.0590,-3.0377);
K1=Kmount(135000,95000,190000,0,pi/2,pi/2,pi/2,(18/180)*pi,(72/180)*pi,pi/2,(108/180)*pi,(18/180)*pi,-0.13138,0.24431,-0.09928);
K2=Kmount(135000,95000,190000,0,pi/2,pi/2,pi/2,(18/180)*pi,(108/180)*pi,pi/2,(72/180)*pi,(18/180)*pi,-0.13868,-0.21069,-0.10028);
K3=Kmount(150000,65000,140000,0,pi/2,pi/2,pi/2,0,pi/2,pi/2,pi/2,0,0.65702,0.00631,-0.30858);
K=K1+K2+K3;
ww=(inv(M))*K;
w=sqrtm(ww);
F=w/(2*pi);
eig(F)
不知道最后几行哪里不对，请做过的高手帮下忙，不胜感谢！

ChaChing 发表于 2009-7-1 18:47

ww=(inv(M))*K; w=sqrtm(ww); F=w/(2*pi); eig(F)
改为下式试试!
ww=(inv(M))*K;w=eig(ww); f=w/(2*pi)

sunjinyuan 发表于 2009-7-2 08:32

谢谢楼上的回复。
我试了下，还是不对，没有对ww开方啊！
另外问下，忽略阻尼和没有忽略阻尼的结果相差大吗？

是直接求解-w^2*=0，得出w的吗？还是需要迭代计算？

[ 本帖最后由 sunjinyuan 于 2009-7-2 08:34 编辑 ]

ChaChing 发表于 2009-7-2 14:13

sorry! forget sqrt
ww=(inv(M))*K;w=eig(ww); f=sqrt(w)/(2*pi)

sunjinyuan 发表于 2009-7-3 15:33

特别谢谢楼上的回复，我试了，还是和以前的结果一样，看来真邪门了，还是要谢谢楼上的热心回答！

xxm1215 发表于 2009-12-29 13:25

回复 楼主 sunjinyuan 的帖子

楼主所用的是MATLAB编写的程序吗？如果是，我觉得你的质量矩阵不对啊，对于动力总成质量矩阵而言，其矩阵应当是6*6的方阵啊．刚度矩阵同样也是6*6的方阵．另外，ww=(inv(M))*K最好改为ww=K/M,因为mablab在计算(inv(M))*K的精度远不如K/M，二者有差别，你可以试试！
页: [1]
查看完整版本: 求助：关于动力总成的特征值求解问题