利用传递矩阵和Riccati传递矩阵法分析转子临界转速.docVIP

利用传递矩阵法和Riccati传递矩阵法分析转子临界转速 所需求解转子参数 将转子简化为如下所示： 三个盘的参数为： QUOTE m1=m2= 另，阶梯轴的三段轴的截面惯性矩分别为： QUOTE J1=1.7 cm4J2 三段轴的单位长度轴段的质量分别为： QUOTE m1=2.47 kgmm 试算转轴的传递矩阵 取试算转速 QUOTE p=ω=1000 rads ; 则，各轴段的传递矩阵分别为： 第1段 QUOTE l=0.06 mJ=1.7×10-8 m4m=2.47 第2段 QUOTE l=0.15 mJ=3.2×10-8 m4 第3段 第4段 第5段 第6段 此6段传递矩阵均采用MATLAB编程求解，MATLAB的源文件为HYPERLINK \l OLE_LINK1H.m 采用传递矩阵法进行各段轴的状态参数的传递 初始参数列阵为: 令 QUOTE X01=1 ，则初始矩阵可化为： QUOTE 1θ0125000θ017e6 以初始矩阵乘第一轴段的传递矩阵，则可得第一段轴的终端状态参数： 由于考虑支座的支撑刚度系数变化从，先取，那么，此处，则可得支座A后第2段的起始端参数阵为： 用第2段的传递矩阵乘此矩阵，可得第2段终端参数： 用中间圆盘的传递矩阵乘第2段终端参数阵，即可得第3段起始端参数： 用第3段传递矩阵乘其始端参数矩阵： 用上式乘以支座刚