帮忙看下程序,关于转子固有频率的
clcclear
%等截面轴参数
l=0.45;%轴的总长
d=0.01;
A=pi*d*d/4;
%轮盘参数
D=0.04;
h=0.02;
%盘轴材料参数(忽略轴的质量)
a=1;%截面剪切系数
u=0.3;%泊松比
rou=7800;%密度
E=2.1e11;%弹性模量
G=E/(2*(1+u));%切变模量
I=pi*(d^4)/64;%圆截面的惯性矩
K=2.1e7;%弹性支撑刚度
v=6*E*I/(a*G*A*l*l);%
mi=rou*pi*D^2/4*h;%轮盘的集质量
Jp=mi*D^2/8; Jd=Jp/2;%转动惯量
Ji=Jp-Jd;
%参数的数组形式
L=;%轴的分段
M=;
K=;
v=;
J=;
k=0;
for w=0:0.01:100;
for i=1:9;
T(:,:,i)=[1 1/K(i);
-J(i)*w^2 1;];
end
H=T(:,:,1);
for i2=2:9;
H=T(:,:,i2)*H;
end
F=H(1,1)*H(2,2)-H(1,2)*H(2,1);
if F*(-1)^k < 0 %求解临界转速
k=k+1;
wi(k)=w;
w=wi(k);
ni(k)=wi(k)*30/pi;
end
end
ni;
wi;
为什么出不了结果? 出错提示!?
take a look
5.建议提问的网友分清 编程问题 和 专业问题 http://www.chinavib.com/thread-36746-1-1.html
6.提问的智慧!!!!(发帖前请认真阅读) http://www.chinavib.com/thread-21991-1-1.html
from http://home.vibunion.com/blog-63979-18250.html 你问的问题太泛了,具体点别人才好给你建议
页:
[1]