传递矩阵法求临界转速

huangziyuan212

将转子集总为15个节点，计算结果不对啊，求大神们指点。。。。。另外对于解析的方法怎么求解这个等截面轴的临界转速？？？

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 l20];

M=[m1 m2m2 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 v1v1 v1 v2 0];

J=[J1 J2 J2 J2 J2 J2 J2 J2 J2 J2 J2J2 J2 J3 J4];

k=0;

Tit=['第一阶频率的振型和弯矩图';'第二阶频率的振型和弯矩图';'第三阶频率的振型和弯矩图'];

forw=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 01];

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

没有人做吗？

wkl1990

if F*(-1)^k < 0 %求解临界转速，我看公式上不是等于0吗，为什么你这样写？

gudanyuzhou

见识一下子

天天狂奔

幸好我心思缜密，火前留名！

                               
登录/注册后可看大图

用不了多久我就会升级加精！当上八五！ 出任吧主！ 迎娶度娘！走上贴吧的巅峰！ 嘿嘿，想想还有点小激动呢

lbtv

1713573225

v1 这是什么参数

yejet

等截面轴的临界转速解析解推导方法及计算公式参考下面的图片

yejet

简单看你的程序，逻辑关系非常混乱，而且没有必要的说明
估计大部分人都不会有时间去帮你调整程序
所以于计算结果还是需要你自己去跟踪分析
1. 查看输入参数是否有问题
2. 循环计算是否有问题
3. 判断条件if F*(-1)^k < 0 %求解临界转速，看不懂附加一份资料，你看看别人是怎么做的吧

yejet

1713573225 发表于 2015-11-19 09:23
v1 这是什么参数

从楼主的程序来看V1应该仅仅是编程过程中出现的一个过程变量