光轴临界转速ansys和理论解差别很大

huangziyuan212

用什么方法能精确计算转子临界转速，Riccati方法、解析解、ansys、有限元方法，怎么才能吻合较好？求教

上课睡觉

建议把你的具体问题贴出来，包括模型及计算代码等
只要模型是合理的，不管采用什么算法，前几阶的临界转速误差都不会太大

huangziyuan212

上课睡觉 发表于 2014-6-6 06:42
建议把你的具体问题贴出来，包括模型及计算代码等
只要模型是合理的，不管采用什么算法，前几阶的临界转速 ...

clc
clear
l1=0.007;
l2=0.008;
d=0.0079;

m1=0.00137;
m2=0.00274;
m3=0.00293;
m4=0.00156;
K1=2.5e7;
A=pi*d*d/4;
a=0.886;
u=0.3;
rou=8000;
E=1.93e11;
G=E/(2*(1+u));
I=pi*(d^4)/64;
v1=6*E*I/(a*G*A*l1*l1);
v2=6*E*I/(a*G*A*l2*l2);

Jp1=10.6487e-3;  
Jp2=21.2974e-3;
Jp3=22.81856e-3;
Jp4=12.16995e-3;

Jd1=-5.823e-3;
Jd2=-11.646e-3;
Jd3=-16.378e-3;
Jd4=-10.555e-3;

J1=Jp1-Jd1;
J2=Jp2-Jd2;
J3=Jp3-Jd3;
J4=Jp4-Jd4;

%数组参数
L=[l1 l1 l1 l1 l1 l1 l1 l1 l1 l1 l1 l1 l1 l2 0];
M=[m1 m2 m2 m2 m2 m2 m2 m2 m2 m2 m2 m2 m2 m3 m4];
K=[K1 0 0 0 0 0 0 0 0 0 0 0 0 0 K1];
v=[v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 v2 0];
J=[J1 J2 J2 J2 J2 J2 J2 J2 J2 J2 J2 J2 J2 J3 J4];
k=0;
Tit=['第一阶频率的振型和弯矩图';'第二阶频率的振型和弯矩图';'第三阶频率的振型和弯矩图'];
for w=0:0.01:4000;
for i=1:15;
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:15;
       H=T(:,:,i2)*H;
end
    F=H(3,1)*H(4,2)-H(3,2)*H(4,1);
    if F*(-1)^k < 0 %求解临界转速
       k=k+1;
       wi(k)=w;
       w=wi(k)
        ni(k)=wi(k)*30/pi;
end
end

huangziyuan212

上课睡觉 发表于 2014-6-6 06:42
建议把你的具体问题贴出来，包括模型及计算代码等
只要模型是合理的，不管采用什么算法，前几阶的临界转速 ...

file:///C:/Users/hp/AppData/Local/Temp/msohtml1/01/clip_image002.gif
这个轴的参数，为什么程序没有解

huangziyuan212

huangziyuan212 发表于 2014-6-12 10:54
这个轴的参数，为什么程序没有解

一根长106mm，直径为7.9mm的光轴，我集总质量后为15个圆盘，前13轴段7mm长，第14轴段为8mm

dollfish000

振动理论书上计算的第一第二阶固有频率，看起来非常繁琐，还有列表求解的，不过是个学习的过程需要了解。误差大和约束方法可能有关系，也和ANSYS中采用的一维梁或三维实体单元有关系。