|
|
马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。
您需要 登录 才可以下载或查看，没有账号？我要加入
x
求救，自己用传递矩阵法算的轴的固有频率，和ANSYS算出的固有频率不一致！找不出哪里错了.这是我编的MATLAB程序：clc
clear
L=[0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035]; %轴段长度
d=0.05; %直径
A=pi*d*d/4; %截面积
%盘轴材料参数
a=0;
k=0;
u=0.3; %泊松比
rou=7850; %密度
E=2.06e11; %弹性模量
K1=2.0e16; %支撑处刚度
G=E/(2*(1+u));
I=pi*(d^4)/64;
for i=1:20;
M(i)=rou*pi*d^2/4*L(i);%轮盘的集质量
I1(i)=0.78*L(i)*(d^4)/1.0e6;
I3(i)=M(i)*(d^2)/8;
v(i)=0;%6*E*I/(a*G*A*L(i)*L(i));%剪切影响系数
Jp=M(i)*d*d/8; %极转动惯量
Jd=Jp/2; %直径转动惯量
J(i)=0;%Jp-Jd;
end
K=[K1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 K1]; %支撑刚度
for w=0:1:20000;
for i=1:20;
T(:,:,i)=[1+(L(i)^3)*(1-v(i))*(M(i)*w^2-K(i))/(6*E*I) L(i)+L(i)^2*J(i)*w^2/(2*E*I) L(i)^2/(2*E*I) L(i)^3*(1-v(i))/(6*E*I);
(L(i)^2)*(M(i)*w^2-K(i))/(2*E*I) 1+L(i)*J(i)*w^2/(E*I) L(i)/(E*I) L(i)^2/(2*E*I);
L(i)*(M(i)*w^2-K(i)) J(i)*w^2 1 L(i);
M(i)*w^2-K(i) 0 0 1];%传递矩阵
end
H=T(:,:,1);
for i2=2:20;
H=T(:,:,i2)*H;
end
F=H(3,1)*H(4,2)-H(4,1)*H(3,2);
if F*(-1)^k < 0 %求解临界转速
k=k+1;
wi(k)=w %固有频率
ni(k)=wi(k)*30/pi;%临界转速
end
end
补充内容 (2013-1-5 09:40):
转子的固有频率，我用ANSYS算的固有频率和论坛算例算一算出的临界转速再乘以pi/30完全不一样啊，是不是固有频率不能用临界转速来换算呢？
这是算列一
算例一：如图所示两端简支的光轴系统，轴长为2m，轴直径为0.1m，转子弹性模量为2.095e11Pa，转子材料密度为7.85e3kg/m3。
用论坛程序算出该转子可以求得临界转速的理论解，理论解为：
一阶临界转速：3043 转/分
那么换算成固有频率为w=3043*pi/30=318.5
但我用ANSYS计算这根轴的固有频率为634，为什么呢？？？？ |
|
[复制链接]
发表于 2012-12-11 17:32
提升卡
变色卡
显身卡
楼主
