转子动力学作业（求解转子系统前三个临界转速和主振型的传递矩阵法）.doc

节选段落一：
%求解转子系统前三个临界转速和主振型的传递矩阵法
clc
clear
%等截面轴参数
l1=0.12;
d=0.04;
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=2.0e7;
v1=6*E*I/(a*G*A*l1*l1);
mi=rou*pi*D^2/4*h;%轮盘的集质量
Jp=mi*D^2/8; Jd=Jp/2;
Ji=Jp-Jd;
%参数的数组形式
L=[l1 l1 l1 l1 l1 l1 l1 l1


节选段落二：
+003 * 0.0150 1.7662
00.20.40.60.811.21.41.6
-1
0
1
第一阶频率的振型和弯矩图
轴长
不平衡值
00.20.40.60.811.21.41.6
-1
0
1
第二阶频率的振型和弯矩图
轴长
不平衡值
00.20.40.60.811.21.41.6
-1
0
1
第三阶频率的振型和弯矩图
轴长
不平衡值
振型
弯矩
00.511.5
-1
0
1
第一阶频率的振型和弯矩图
轴长
不平衡值
00.511.5
-1
0
1
第二阶频率的振型和弯矩图
轴长
不平衡值
00.511.5
-1
0
1
第三阶频率的振型和弯矩图
轴长
不平衡值
振型


节选段落三：
弯矩
00.511.5
-1
0
1
第一阶频率的振型和弯矩图
轴长
不平衡值
00.511.5
-1
0
1
第二阶频率的振型和弯矩图
轴长
不平衡值
00.511.5
-1
0
1
第三阶频率的振型和弯矩图
轴长
不平衡值
振型
弯矩
00.20.40.60.811.21.41.6
-1
0
1
2
3
x 10
-6
wi(1)时的振动响应图
轴长
不平衡值
00.20.40.60.811.21.41.6
-5
0
5
10
x 10
-6
wi(2)时的振动响应图
轴长
不平衡值
00.20.40.60.811.21.41.6
-1
0
1
2
x 10
-6
wi(1)时的振动响应图