奔腾6 发表于 2007-7-13 09:57

转子动力学中编程问题讨论

我用此程序验算书上的算例（闻邦椿的《转子动力学》），结果相差很大，我反复检查过很多遍程序实在发现不了新的问题，如果前辈对此方面有所研究，请不吝赐教啊，很困惑，完整程序如下：
clear;
l=;
m=;
Jp=;
Jd=;
I=;
E=4393;
v=;
k=;
kb=;
mb=;
S=;
s=1;
x=[];
for n=1864:1864
for i=1:13
K=k(i)*(kb(i)-mb(i)*n^2)/(k(i)+kb(i)-mb(i)*n^2)
%l(i)/(E*I(i))
u11=;
u12=;
u21=(l(i)/(E*I(i))).*;
u22=;
%s=sign(det(u21*S+u22))*s;
S=*inv()
end
D=det(S);
x=;
end
n=1864:1:1864;
grid on
plot(n,x)

[ 本帖最后由 eight 于 2007-7-13 13:34 编辑 ]

eight 发表于 2007-7-13 13:33

原帖由 奔腾6 于 2007-7-13 09:57 发表 http://www.chinavib.com/forum/images/common/back.gif
我用此程序验算书上的算例（闻邦椿的《转子动力学》），结果相差很大，我反复检查过很多遍程序实在发现不了新的问题，如果前辈对此方面有所研究，请不吝赐教啊，很困惑，完整程序如下：
clear;
l=

请把你的问题叙述清楚，别人的时间也是很宝贵的。建议先看看各个置顶贴，特别是：聚宝盆

奔腾6 发表于 2007-7-13 15:26

回复 #2 eight 的帖子

上面的程序运行没有问题，但是得出的曲线不是标准的剩余量曲线，不知是哪里的问题，请对转子动力学有所研究的前辈给与指点。
页: [1]
查看完整版本: 转子动力学中编程问题讨论