谁说我不会 发表于 2016-3-31 17:17

帮忙看下程序,关于转子固有频率的

clc
clear
%等截面轴参数
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;



为什么出不了结果?

ChaChing 发表于 2016-3-31 18:44

出错提示!?

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

hustxyong 发表于 2016-4-1 09:06

你问的问题太泛了,具体点别人才好给你建议
页: [1]
查看完整版本: 帮忙看下程序,关于转子固有频率的