基于Hamilton原理的转子-汽封-轴承系统非线性动力学问题研究.docxVIP

基于Hamilton原理的转子-汽封-轴承系统非线性动力学问题研究

一、研究背景与意义

在旋转机械领域，转子-汽封-轴承系统作为核心组成部分，其运行稳定性直接决定设备的安全与效率。随着机组向高转速、大功率方向发展，系统中非线性因素（如汽封间隙的流体激振力、轴承油膜的非线性支撑力、转子材料的非线性变形等）的影响愈发显著，易引发系统共振、碰摩等故障。

Hamilton原理作为分析动力学问题的重要工具，通过能量变分形式建立系统运动方程，可有效整合多部件的非线性能量项，为复杂系统动力学建模提供统一框架。开展基于该原理的系统非线性动力学研究，能精准揭示系统动态行为演化机制，为旋转机械的优化设计与故障预警提供理论支撑。

二、Hamilton原理在系统建模中的适配性

Hamilton原理的核心是：完整、保守系统的真实运动轨迹，满足在给定时间区间内，系统的作用量（动能与势能差值的积分）取极值，其数学表达式为：

\delta\int_{t_1}^{t_2}(T-V)dt=0

其中，T为系统总动能，V为总势能。

对于转子-汽封-轴承系统，该原理具有天然适配性：

动能整合：可统一描述转子的旋转动能、平动动能，以及汽封流体、轴承油膜的运动动能，避免多部件动能分析的碎片化。

势能覆盖：能涵盖转子的弹性势能、轴承油膜的挤压势能、汽封流体的压力势能等非线性势能项，完整反映系统能量存储与转换关系。

非线性项嵌入：通过变分过程，可直接将汽封间隙的非线性流体力、轴承油膜的非线性刚度阻尼等因素，以能量导数形式融入运动方程，简化非线性项的耦合过程。

三、系统非线性因素分析

（一）汽封系统的非线性

汽封间隙随转子振动呈周期性变化，导致流体激振力呈现非线性特征。基于流体力学理论，汽封间隙内的压力分布可通过雷诺方程描述，而压力梯度产生的流体激振力随间隙变化呈现非线性关系，需在能量建模中引入非线性势能项。

（二）轴承系统的非线性

滑动轴承的油膜支撑力具有强非线性，其刚度与阻尼系数随转子偏心率、转速变化而动态调整。采用短轴承理论或有限元法求解油膜压力时，油膜力与转子位移、速度的关系需通过非线性函数拟合，在Hamilton原理的势能计算中需重点考虑油膜力的非线性能量耗散与存储特性。

（三）转子系统的非线性

转子的材料非线性（如弹性模量随应力变化）、几何非线性（如大变形引发的附加弯矩），以及转子与汽封、轴承的碰摩非线性（如碰撞时的接触力与能量损耗），均会导致系统动能与势能的非线性变化，需通过分段函数或非线性本构关系在能量模型中体现。

四、基于Hamilton原理的动力学建模流程

系统能量分解

动能T：包含转子的旋转动能（\frac{1}{2}J\omega^2，J为转动惯量，\omega为角速度）、平动动能（\frac{1}{2}m\dot{x}^2+\frac{1}{2}m\dot{y}^2，m为转子质量，x,y为横向位移），以及汽封流体、轴承油膜的运动动能。

势能V：涵盖转子弹性势能（\frac{1}{2}k_xx^2+\frac{1}{2}k_yy^2，k_x,k_y为转子刚度）、轴承油膜势能、汽封流体压力势能，以及非线性因素引发的附加势能（如碰摩势能、材料非线性势能）。

变分运算与运动方程推导

根据Hamilton原理，对作用量\int_{t_1}^{t_2}(T-V)dt求变分并令其为零，结合拉格朗日方程\frac{d}{dt}(\frac{\partialL}{\partial\dot{q}})-\frac{\partialL}{\partialq}=Q（L=T-V为拉格朗日量，q为广义坐标，Q为非保守力），推导系统非线性运动微分方程，形式如下：

m\ddot{x}+c_x(\dot{x},\dot{y})\dot{x}+k_x(x,y)x+F_{x,\text{seal}}(x,y)+F_{x,\text{bearing}}(x,y,\dot{x},\dot{y})=0

m\ddot{y}+c_y(\dot{x},\dot{y})\dot{y}+k_y(x,y)y+F_{y,\text{seal}}(x,y)+F_{y,\text{bearing}}(x,y,\dot{x},\dot{y})=0

其中，c_x,c_y为非线性阻尼系数，F_{\text{seal}},F_{\text{bearing}}分别为汽封流体激振力与轴承油膜力。

边界条件与参数确定

结合转子的支撑方式（如固定端、滚动轴承支撑）、汽封间隙设计值、轴承结构参数（如油膜厚度、供油压力），确定方程的初始条件（如初始位移、角速度）与边界