huangziyuan212 发表于 2014-6-18 16:19

传递矩阵法求临界转速


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

clcclearl1=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=;M=;K=;v=;J=; k=0;Tit=['第一阶频率的振型和弯矩图';'第二阶频率的振型和弯矩图';'第三阶频率的振型和弯矩图'];forw=0:0.01:4000;for i=1:15;T(:,:,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];endH=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;endend

huangziyuan212 发表于 2014-6-20 10:21

没有人做吗？{:{19}:}

wkl1990 发表于 2014-8-28 19:44

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

gudanyuzhou 发表于 2014-9-1 18:36

见识一下子

天天狂奔 发表于 2014-9-3 09:13

幸好我心思缜密，火前留名！www.epzw.la/files/article/html/57/57599/用不了多久我就会升级加精！当上八五！ 出任吧主！ 迎娶度娘！走上贴吧的巅峰！ 嘿嘿，想想还有点小激动呢

lbtv 发表于 2014-9-3 15:09

{:{39}:}{:{39}:}

1713573225 发表于 2015-11-19 09:23

v1 这是什么参数

yejet 发表于 2015-11-20 12:44

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




yejet 发表于 2015-11-20 12:50

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

yejet 发表于 2015-11-20 12:55

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

从楼主的程序来看V1应该仅仅是编程过程中出现的一个过程变量
页: [1]
查看完整版本: 传递矩阵法求临界转速