奔腾6 发表于 2007-7-9 11:02

请教转子传递矩阵

帮忙看一个转子动力学传递矩阵程序，一直得不到标准的剩余量曲线，还出现振荡，如果哪位前辈以前编过此程序，麻烦给传一个，解燃眉之急啊，多谢谢谢谢谢！！！！！
clear;
l=;
m=;
Jp=;
Jd=;
I=;
K=;
E=2.000e+011;
v=;
S=;
s=1;
x=[];
for n=200:800
for i=1:5
u11=;
u12=;
u21=(l(i)/(E*I(i))).*;
u22=;
s=sign(det(u21*S+u22))*s;
S=*inv();
end
D=det(S)*s;
x=;
end
n=200:1:800;
grid on
plot(n,x)

xsxsxs 发表于 2007-7-12 10:46

riccat 传递矩阵法 计算得到的剩余l量就是振荡的i

奔腾6 发表于 2007-7-12 11:01

我用此程序验算书上的算例（闻邦椿的《转子动力学》），结果相差很大，我反复检查过很多遍程序实在发现不了新的问题，如果前辈对此方面有所研究，请不吝赐教啊，很困惑，完整程序如下：
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)

vibration 发表于 2007-7-14 06:36

加注释

wenzihui 发表于 2007-7-15 22:12

Jp=;
Jd=;
I=;
E=4393;
v=;
你能解释一下这些的含义吗？？
感觉不太对！！

奔腾6 发表于 2007-7-21 10:23

回复 #5 wenzihui 的帖子

Jp=;
Jd=;
以上两组数据全设为0，是因为不计圆盘的转动惯量及陀螺力矩

v=;
令v等于零是不计剪切变形的影响

I=;
E=4393;
因为书上给的算例是l/(EI),而各点的E*I的积是一样的，所以就写成上面的形式
页: [1]
查看完整版本: 请教转子传递矩阵