李 铮 ， 李 冰 ， 高宇星 ， 苏来博 ， 徐武彬

（1.广西科技大学机械与汽车工程学院，广西 柳州 545006；2.广西土方机械协同创新中心，广西 柳州 545006）

0 引言

由于转子系统制造误差存在的客观性，文献[4]和文献[5]使用有限差分法求解存在圆度误差的转子动力学模型，求解油膜力、轴心轨迹以及系统临界转速。考虑尺寸和形位误差及轴瓦表面弹性变形量的影响，使用三维形状误差模型取代二维模型分析误差滑动轴承性能[6-10]。

油膜力是造成轴承转子系统非线性化的主要因素，同时对转子系统动静态特性和运行稳定性产生重要影响。在求解油膜力表达式时，Yaser Bastani等[11]将修正函数应用于无限长（IL）或无限短（IS）的力模型，并通过数值方法求解二维雷诺方程逼近有限长轴承的油膜力计算表达式，预测系统在静态和不平衡载荷下的稳态行为。数值方法求解油膜力是目前研究转子动力学特性的主要方法，但其计算时间长、工作量大的缺点明显不能适应快速求解油膜特性的要求。因此，Dimitris Sfyris 等[12]将雷诺方程拆分为四个线性常微分方程，并与边界条件一起构成四个二阶非恒定系数的线性常微分方程。使用变量分离的方法将雷诺方程的解变换为一组特解和通解的加法和乘法形式进行近似解析解求解。张永芳和黑棣等[13-14]在分离变量的基础上结合变分原理将有限长和无限长径向滑动轴承油膜压力用其轴向和周向分离函数的乘积表示，求得油膜力近似解析解来提高计算效率。……