求转子扭转振动临界转速计算的matlab程序
大神们，帮一下小弟吧 本帖最后由 weilong 于 2012-6-6 09:55 编辑有没有好心人帮忙？在线等 本帖最后由 weilong 于 2012-6-6 15:30 编辑
clc
clear
%等截面轴参数
l1=0.062;l2=0.07;l3=0.263;l4=0.016;l5=0.164;l6=0.109;
l7=0.169;l8=0.095;l9=0.016;l10=0.06;l11=0.185;l12=0.085;l13=0.105;l14=0.015;l15=0.04;
d=0.05;A=pi*d*d/4;
%轮盘参数
D=0.5;
h=0.025;
%盘轴材料参数（忽略轴的质量）
a=1;
u=0.3;
rou=7800;
E=2.0e11;
G=E/(2*(1+u));
I=pi*(d^4)/64;
K1=pi*Gd^4/(32*l1);K2=pi*Gd^4/(32*l2);K3=pi*Gd^4/(32*l3);K4=pi*Gd^4/(32*l4);K5=pi*Gd^4/(32*l5);
K6=pi*Gd^4/(32*l6);K7=pi*Gd^4/(32*l7);K8=pi*Gd^4/(32*l8);K9=pi*Gd^4/(32*l9);K10=pi*Gd^4/(32*l10);
K11=pi*Gd^4/(32*l11);K12=pi*Gd^4/(32*l12);K13=pi*Gd^4/(32*l13);K14=pi*Gd^4/(32*l14);K15=pi*Gd^4/(32*l15);
v1=6*E*I/(a*G*A*l1*l1);v2=6*E*I/(a*G*A*l2*l2);v3=6*E*I/(a*G*A*l3*l3);v4=6*E*I/(a*G*A*l4*l4);v5=6*E*I/(a*G*A*l6*l6);v6=6*E*I/(a*G*A*l6*l6);v7=6*E*I/(a*G*A*l7*l7);v8=6*E*I/(a*G*A*l8*l8);v9=6*E*I/(a*G*A*l9*l9);v10=6*E*I/(a*G*A*l10*l10);v11=6*E*I/(a*G*A*l11*l11);v12=6*E*I/(a*G*A*l12*l12);v13=6*E*I/(a*G*A*l13*l13);v14=6*E*I/(a*G*A*l14*l14);
v15=6*E*I/(a*G*A*l15*l15);
mi=rou*pi*D^2/4*h;%轮盘的集质量
Jp=mi*D^2/8; Jd=Jp/2;
Ji=Jp-Jd;
%参数的数组形式
L=;
M=;
K=;
v=;
J=;
{:{39}:} 本帖最后由 weilong 于 2012-6-6 18:49 编辑
求助啊,,,, 本帖最后由 weilong 于 2012-6-7 22:55 编辑
这就是要计算的矩阵，哪位大神帮帮小弟？ 回复 6 # weilong 的帖子
你的问题是什么啊？ 就是如何实现矩阵运算？？？
大哥帮一下小弟 回复 8 # weilong 的帖子
提问时，尽量让别人知道你要问什么。
矩阵运算？加减乘除都是运算，要做什么？
我猜你要要用传递矩阵法算扭转频率，是么？ 是啊，能帮我弄一下吗？大哥 回复 7 # yyxt007 的帖子
是啊，能帮我弄一下吗？大哥 回复 10 # weilong 的帖子
根据边界条件，可求得系统的频率方程式。方程的阶次较高，一般采用频率搜索法求解。该方法的具体描述，请参考闻邦椿院士等编写的《高等转子动力学》。 回复 12 # yyxt007 的帖子
会者不难啊，您如果方便的话，帮我做一下，，，好人好事啊。 大哥，您教育的很对啊，但是，我也有纠结之处啊，还有很多事要处理，您帮帮忙，小弟感激不尽啊 回复 12 # yyxt007 的帖子
大哥，您教育的很对啊，但是，我也有纠结之处啊，还有很多事要处理，您帮帮忙，小弟感激不尽啊
页:
[1]
2