刚性轴颈轴承特性数值分析工具

资 源 简 介

该程序是一个专门用于计算滑动轴颈轴承特征参数的MATLAB工具，其核心设定是在计算过程中忽略轴承衬套和轴颈的弹性变形，将其视为理想的刚性体。程序的核心流程通过数值求解流体动力润滑的雷诺方程实现，主要应用于计算在稳态运行条件下的轴承性能指标。实现方法上，系统采用有限差分法对轴承展开面的二维油膜压力场进行离散建模，通过迭代算法直至压力场收敛。其核心功能包括自动计算偏心率、偏心角、油膜厚度分布、油膜压力分布等基本物理量，并在此基础上进一步导出轴承的承载力、摩擦力、侧漏流量以及在微小扰动下的八个动力特性系数，包括

详 情 说 明

刚性轴颈轴承特性数值分析工具

项目介绍

本工具是一个基于 MATLAB 开发的数值仿真程序，旨在对滑动轴颈轴承进行稳态性能分析与动力特性评估。程序核心基于流体动力润滑理论，将轴承衬套和轴颈视为刚性体，不计其数其受压后的弹性变形。通过数值方法求解二维雷诺方程（Reynolds Equation），本工具能够精确模拟轴承在特定工况下的油膜压力分布，并导出工程设计所需的关键性能参数及转子动力学分析必备的线性刚度和阻尼系数。

功能特性

1. 平衡位置自动搜寻：程序能够根据设定的外部载荷，通过内置的迭代算法自动寻找对应的平衡偏心率和偏心角，确保计算结果符合实际受力工况。
2. 流场压力数值模拟：采用有限差分法（FDM）对轴承展开面的雷诺方程进行离散化求解，支持自定义周向与轴向的网格密度。
3. 全面性能参数导出：除压力分布外，还可计算最小油膜厚度、侧漏流量、摩擦功耗及轴承实际承载力。
4. 动力特性分析：利用微小扰动法计算轴承在平衡位置处的八个动力特性系数，包括四个线性刚度系数（Kxx, Kxy, Kyx, Kyy）和四个线性阻尼系数（Cxx, Cxy, Cyx, Cyy）。
5. 多维度数据可视化：提供压力分布三维图、等值线图、中剖面压力曲线以及油膜厚度极坐标分布图，方便直观分析物理场。

实现逻辑与系统流程

1. 参数定义与网格划分

程序首先配置轴承的几何特征（直径、长度、间隙）、工况参数（转速、润滑油粘度）以及目标载荷。随后，将轴承内表面在周向（theta）和轴向（z）进行网格离散，构建用于计算的坐标矩阵。

2. 平衡位置寻优算法

由于偏心率和偏心角不是直接给定的，程序启动一个双变量迭代过程：

• 给定初始偏心率和偏心角。
• 调用雷诺方程求解模块计算当前的油膜力。
• 计算计算力与目标载荷的残差，以及力的合力方向与偏心方向的偏差。
• 采用步长控制（alpha, beta）更新偏心率和偏心角，直至载荷平衡条件满足收敛精度。

3. 雷诺方程数值求解 (算法核心)

核心求解器实现了以下逻辑：

• 离散格式：将经典的 Reynolds 方程通过有限差分转换为代数方程组。
• 挤压项处理：考虑轴颈中心的速度项（x_dot, y_dot），以支持动力系数计算中的挤压效应模拟。
• 超松弛迭代 (SOR)：采用超松弛法提高压力场求解的收敛速度，并结合雷诺边界条件（Swift-Stieber 条件），通过在迭代过程中强制负压区域归零来模拟油膜破裂（空穴效应）。
• 边界条件：轴向两端压力为零，周向满足连续性周期边界条件。

4. 动力特性系数计算

程序在平衡位置附近施加微小的位移扰动（delta*C）和速度扰动（delta*omega*C）：

• 通过计算位移扰动后的压力场变化，利用差分近似求导，得到四个刚度系数。
• 通过计算速度扰动带来的挤压油膜力变化，求得四个阻尼系数。
• 最后将各系数进行无量纲化处理，便于不同规格轴承的性能对比。

5. 后处理与集成计算

在获得压力场后，通过数值积分法（矩形或梯形积分）计算轴向各截面的受力方向分量。同时，根据膜厚梯度和压力梯度进一步导出摩擦力和侧漏量。

关键算法说明

• 有限差分法 (Finite Difference Method)：程序将偏微分方程转化为包含中心差分格式的代数方程，利用网格节点间的相互关系描述压力场。
• 微小扰动法：这是一种线性化技术，假定在平衡点附近的微小扰动下，轴承力与位移/速度呈线性关系，从而获得刚度和阻尼。
• 迭代收敛判定：采用相对残差作为压力场和平衡位置迭代的终止准则，确保了数值结果的稳定性。

系统要求

• 环境要求：MATLAB R2016a 或更高版本。
• 工具箱：本程序完全基于标准 MATLAB 内置函数编写，无需安装额外的工具箱（如 Optimization Toolbox）。
• 硬件建议：由于包含多层嵌套迭代，建议内存不低于 8GB 以保证大规模网格划分时的计算效率。

使用方法

1. 打开 MATLAB 环境，将程序源代码所在的文件夹设置为当前工作目录。
2. 在脚本顶部的“输入参数设置”区域，根据实际轴承设计需求修改直径（D）、长度（L）、半径间隙（C）、转速（N）以及外部载荷（W_target）等参数。
3. 直接点击运行（Run）脚本。
4. 程序将在“命令行窗口”实时打印寻优后的偏心率、偏心角以及无量纲动力特性系数。
5. 程序运行结束会自动弹出四个子图组成的图形窗口，展示轴承详尽的物理特性分布。