基于Prohl与Ricatti传递矩阵法的转子固有频率计算平台

资 源 简 介

本项目是专门针对旋转机械转子动力学设计的MATLAB计算程序，主要功能是应用传递矩阵法（TMM）精确求解复杂转子系统的固有频率、临界转速以及相应的振型。 程序完整实现了Prohl传递矩阵理论，通过构建状态向量和单步传递矩阵，将转子物理模型离散化为包含圆柱轴段、圆锥轴段、轮盘及柔性支承的耦合系统。 为了克服传统传递矩阵法在计算高阶频率时由于截断误差累积导致的数值不稳定问题，本项目重点引入了Ricatti传递矩阵算法。 Ricatti法通过将状态向量之间的线性关系转化为非线性递推关系，有效避免了行列式搜索过程中的病态矩阵问题，显著提升了对于多支承、变截面长转子系统的计算精度和数值稳定性。 系统支持多自由度支承模拟，能够分析包括陀螺力矩、剪切变形以及转动惯量在内的多种动力学因素对系统特性的影响。 该平台适用于航空、电力、机械制造等领域中对各类高速泵、压缩机、燃气轮机转子系统的振动特性预估与故障分析。

详 情 说 明

基于Prohl与Ricatti传递矩阵法的转子系统固有频率计算平台

项目介绍

本项目是一个基于MATLAB开发的转子动力学计算工具，旨在通过传递矩阵法（Transfer Matrix Method, TMM）对复杂转子系统进行建模与振动特性分析。程序集成了传统的Prohl传递矩阵法与数值稳定性更高的Ricatti传递矩阵法，通过对计算结果的对比，直观展示不同算法在处理多支承、变截面转子系统时的数值表现。该平台能够精确计算系统的固有频率、识别临界转速，并实现转子结构可视化与前三阶振型的仿真绘制。

功能特性

• 双算法并行计算：同步实现Prohl法与Ricatti法，用于对比分析长轴系计算过程中的数值稳定性问题。
• 复杂结构建模：支持由多个圆柱/圆锥轴段、集中质量轮盘及弹性支承组成的耦合系统建模。
• 动力学因素完备：在建模中考虑了Timoshenko梁模型的旋转惯量、剪切变形，以及轮盘的陀螺效应。
• 自动搜索与提取：通过频率扫掠与过零点检测算法，自动提取指定范围内的系各阶固有频率。
• 多维度可视化报告：

- 算法数值稳定性对比曲线（归一化行列式值）。 - 转子物理结构离散化示意图。 - 前三阶标准化振型图。

系统要求

• 软件环境：MATLAB R2016b 或更高版本。
• 硬件环境：兼容主流计算设备，无需特殊显卡支持。

实现逻辑与功能模块说明

程序的设计遵循典型的转子动力学分析流程，逻辑分为参数定义、矩阵构建、频率搜索、结果提取与可视化五个核心步骤：

1. 转子系统参数定义

程序预设了标准的材料属性（如钢材的弹性模量、密度、剪切模量）和截面剪切系数。用户通过数组定义轴段的几何尺寸（长度、内外径）、轮盘的物理属性（安装节点、质量、极转动惯量、直径转动惯量）以及支承点的径向刚度。

2. 单元基础矩阵构建

程序根据状态向量 [位移 $w$, 转角 $phi$, 弯矩 $M$, 剪力 $Q$] 构建以下单元矩阵：

• 轴段矩阵：基于Timoshenko梁理论，包含了考虑剪切变形与旋转惯量的传递项，描述轴段两端状态向量的线性转换。
• 轮盘矩阵：引入质量项与陀螺力矩项，反映旋转状态下轮盘对系统弯曲振动的影响。
• 支承矩阵：将支承离散为作用于特定节点的点刚度。

3. Prohl 传递矩阵算法实现

该算法通过将所有单元矩阵按顺序累乘，得到整个转子系统的总传递矩阵。根据边界条件（如两端自由，即 $M=0, Q=0$），提取总矩阵的子块并计算其行列式。当行列式值为零时，对应的频率即为固有频率。

4. Ricatti 传递矩阵算法实现

为了解决高阶频率计算时可能出现的数值溢出或病态矩阵问题，程序引入了Ricatti变换。其核心逻辑是将状态向量分为位移项和受力项，通过非线性递推公式更新关联矩阵 $S$： $S_{i+1} = (u_{21} + u_{22} S_i)(u_{11} + u_{12} S_i)^{-1}$ 该方法避免了直接累乘导致的量级失控，显著提升了计算长轴系时的精度。

5. 频率搜索与振型恢复

• 根搜索：在设定的频率范围内进行步进扫掠，通过检测行列式值的符号变化和插值补偿，精确锁定固有频率点。
• 振型恢复：在确定固有频率后，程序回到初始节点，利用已知的传递关系逐节点恢复位移分布，并进行归一化处理。

6. 结果展示

程序自动生成绘图窗口，对比Prohl法与Ricatti法行列式的收敛情况，并给出转子的轴向剖面结构图。最后，通过彩色曲线展示转子在各阶固有频率下的相对振幅分布。

使用方法

1. 打开MATLAB，将当前工作路径设置为程序所在目录。
2. 运行计算脚本，程序将开始执行频率扫掠。
3. 在命令行窗口查看提取到的各阶固有频率（rad/s）。
4. 观察自动生成的图形窗口，分析转子的结构布局、算法稳定性对比以及振型仿真结果。
5. 根据实际需求，修改脚本开头的几何参数和材料属性，即可对不同的转子模型进行分析。