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

《利用传递矩阵法和Riccati传递矩阵法分析转子临界转速》由会员分享，可在线阅读，更多相关《利用传递矩阵法和Riccati传递矩阵法分析转子临界转速(13页珍藏版)》请在人人文库网上搜索。

1、利用传递矩阵法和Riccati传递矩阵法分析转子临界转速一、 所需求解转子参数将转子简化为如下所示：三个盘的参数为：另，阶梯轴的三段轴的截面惯性矩分别为： 三段轴的单位长度轴段的质量分别为：二、 试算转轴的传递矩阵取试算转速 ; 则，各轴段的传递矩阵分别为：第1段第2段第3段第4段第5段第6段此6段传递矩阵均采用MATLAB编程求解，MATLAB的源文件为H.m三、 采用传递矩阵法进行各段轴的状态参数的传递初始参数列阵为:令，则初始矩阵可化为：以初始矩阵乘第一轴段的传递矩阵，则可得第一段轴的终端状态参数：由于考虑支座的支撑刚度系数变化从，先取，那么，此处，则可得支座A后第2段的起始端参数阵为：。

2、用第2段的传递矩阵乘此矩阵，可得第2段终端参数：用中间圆盘的传递矩阵乘第2段终端参数阵，即可得第3段起始端参数：用第3段传递矩阵乘其始端参数矩阵：用上式乘以支座刚度矩阵，得其终端参数：则，根据可得： ,则可得支座B后第4段的起始端参数阵为：同上，用此段轴的传递函数乘其起始端的状态参数，可得：则，根据可得： 则，可得第5段的起始参数矩阵：其中，为铰链处的转角。用第5段的传递矩阵乘此参数矩阵，即得第5段的终端参数：用上式乘以支座刚度矩阵可得第6段的初始状态参数阵：则，用第6段的传递矩阵乘此状态参数即可得其终端的参数阵：根据最右边盘得传递矩阵，可得转子终端的状态参数：则根据终端的自由状态，则应该;通。

3、过令解出，并将其带入到的表达式中，可得：；此处使用的MATLAB源程序为calculate.m在MATLAB中使用线性插值法寻找最佳p值使得逼近于0。其程序为rotor.m经计算，考虑支撑刚度变化之间时，取时，一阶临界转速值为取时，一阶临界转速值为取时，一阶临界转速值为取时，一阶临界转速值为取时，一阶临界转速值为取时，一阶临界转速值为因此随着刚度的增加，一阶临界转速的值越来越大，而当不考虑支座的刚度变化，