摘要：本文建立了挤压油膜阻尼器-滑动轴承转子模型，采用Runge-Kutta法求解系统的动力学响应问题，进而研究了转速和挤压油膜阻尼器弹性刚度对系统动力响应的影响。研究表明：随着转速增大，系统经历了从周期运动、准周期运动到混沌的过程；挤压油膜阻尼器内支承刚度越大，系统响应越弱，刚度较小时，系统会进入准周期运动形态。研究结果可为挤压油膜阻尼器转子系统的结构参数设计提供参考。
关键词：挤压油膜阻尼器；转子；支承刚度；振动
中图分类号：O322 文献标志码：A 文章编号：
Research on Dynamic Response of SFD-Sliding Bearing Rotor System with Parameters
WANG Ji-yan1
(School of Science, Shandong Jiaotong University, Jinan 250357, China)
Abstract: This paper established a model of SFD-bearing rotor, solved dynamic response question of the rotor system by adopting Runge-Kutta method, and analyzed the influences of speed and stiffness of elastic support of SFD on dynamic response of the system. The study has shown: as the speed increases, the system is in the process changing from periodic motion, quasi-periodic motion, to chaos; the greater the support stiffness, the weaker the system response, and when the stiffness is relatively small, the system will go into the state of periodic motion. These results can provide the reference for designing structure parameters of SFD rotor systems.
Keywords: squeeze film damper; rotor; stiffness; vibration
0 引言
旋转机械在机械领域中应用非常广泛。特别地，如果旋转机械振动过大，会导致机器不能满足工作要求，甚至会造成重大事故。挤压油膜阻尼器具有体积小、重量轻等优点，并且有明显的减振效果，因此常用于轴承转子结构中。很多学者对带挤压油膜阻尼器的轴承转子系统进行了深入研究。吕晓光等[1]研究了挤压油膜阻尼器轴承转子系统的动力学响应，发现系统由周期运动到混沌的响应变化。王红瑾等[2]在考虑了盘的转角变化的情况下，建立了带挤压油膜阻尼器的转子系统模型，讨论了转子转速、两盘偏心相位差及质量比对系统响应的影响。赵项伟等[3]建立了静偏心挤压油膜阻尼器轴承转子动力学模型，运用数值积分获取了系统的非线性响
应，研究了静偏心对转子轴心轨迹的影响，得到静偏心时轴心轨迹近似为椭圆的结论。
但是，并非任意的挤压油膜阻尼器都能起到提高旋转机械系统稳定性和减振的作用，参数设计是否合理是一个非常重要的原因。影响转子动力响应的因素有很多，如质量因素、阻尼因素和刚度因素等，本文主要研究转速和挤压油膜阻尼器的弹性刚度对转子动力响应的影响。
本文对挤压油膜阻尼器滑动轴承轴承转子系统进行了数学建模，利用Matlab软件计算转子系统的动力响应，研究转子转速和挤压油膜阻尼器内支承刚度的变化对轴承转子系统动力响应的影响，为转子轴承系统的稳定性研究提供参考。
1建模
(a)转子结构简图 (b) 转子结构侧视图
1 挤压油膜阻尼器滑动轴承转子系统结构图
挤压油膜阻尼器滑动轴承转子系统的机构简图[4]见图1(a)，图1(b)为系统结构侧视图，考虑轴承质量，为SFD中心，建立坐标系，为轴承体中心，建立与轴承体的连体坐标系，为转轴中心，同样建立与转轴一起运动的连体坐标系，为转盘中心。建立如下动力学方程：
其中，方程组（1）中为轴承体质量，为SFD内弹性支承刚度，和分别为SFD作用在轴承上沿和上的油膜力，和为滑动轴承作用在转轴上沿和上的油膜力，为轴承体连体坐标系转过的极角，为轴颈连体坐标系转过的极角。
2动力响应
MATLAB 软件是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。该软件包含大量计算算法的集合拥有600多个工程中要用到的数学运算函数，可以方便的实现用户所需的各种计算功能，广泛应用于科学研究和工程技术的各个领域。
本文采用Runge-Kutta法，通过MATLAB编程，计算系统的动力响应，为保证获得稳态响应，需要选择合适的步长。在具体的算例计算过程中，采用的相关参数见表1。
表1 转子参数
2.1转速对系统响应的影响
系统的基本参数不变，阻尼器内刚度为K1=2.45×108 N/m，在转速改变的情况下，分析计算系统的动力响应，见图2-4。系统在低转速时期会发生稳定的周期响应，如图2表示角速度为300rad/s时转盘的轨迹图和庞加莱映射图。转速增大，系统运动逐渐从稳定的周期响应发展为准周期运动，图3显示，转速为600 rad/s时转盘的轨迹图和庞加莱映射图，说明响应为准周期运动。转速继续增加，系统响应进入混沌运动状态。图4(a)表示转速为1000 rad/s时，转盘的运动轨迹，图4(b)表示转盘中心的庞加莱映射图，很显然，系统响应进入混沌。
系统动力响应会受转速的影响，从稳定的周期运动进入准周期，最终到混沌状态。
2.2阻尼器刚度对系统响应的影响
图5-8表示转速为500 rad/s，阻尼器刚度在一定范围内变化时转盘的动力响应轨迹图。可以看出，在相同转速下，阻尼器内弹性支承刚度越大，系统响应越小，符合一般弹性支承变化规律。同时，弹性支承刚度发生改变时，转盘运动轨迹也相应发生变化。当K1=2.45×106 N/m，系统发生拟周期响应(见图5)。K1大于2.45×107 N/m时，系统进行周期运动(见图6-8)。
(a)转盘中心轨迹 (b)转盘中心poincare映射
图8 K1=2.45×109 N/m系统动力响应
3结论
本文研究了挤压油膜阻尼器轴承转子系统的动力响应，通过计算分析，可得到以下结论：
（1）转盘转速的变化会影响系统的动力响应，系统经历了从周期运动到准周期，最终到达混沌的复杂运动；
（2）挤压油膜阻尼器刚度变化会影响系统的动力响应。挤压油膜阻尼器刚度的增大，会抑制转子的振动，刚度较小时系统会出现拟周期响应。
可为挤压油膜阻尼器转子系统的研究提供参考。
参考文献：
[1] 吕晓光, 赵玉成, 卢纪.挤压油膜阻尼器-滑动轴承-柔性转子系统的动力响应分析[J] .应用力学学报, 2007, 24(3):460-463.
[2] 王红瑾，秦卫阳，杨树林，邓旺群.带挤压油膜阻尼器双盘转子的参数变化对系统响应的影响. [J] .航空动力学报, 2009, 24(11):2545-2550.
[3] 赵项伟，罗贵火，王飞. 静偏心对挤压油膜阻尼器减振特性影响的数值分析[J] .航空发动机，2018, 44(3): 42-48.
[4] 王继燕, 赵玉成, 徐惠. SFD-转子系统动力响应分析[J].煤矿机械, 2012, 33(4): 98-101.