含阻尼转子动力学传递矩阵法分析与仿真程序

资 源 简 介

本项目开发了一套基于MATLAB的转子动力学分析系统，核心算法采用改进的传递矩阵法（Transfer Matrix Method, TMM），其显著特点在于能够完整计入系统的阻尼参数。该程序突破了传统传递矩阵法在教学和简单应用中通常忽略阻尼或仅作刚性支承假设的限制，通过引入复数域的传递矩阵算法，实现了对非线性阻尼力和结构阻尼的精确建模。系统将转子离散为集总质量模型（Lumped Mass Model），包括无质量弹性轴段和刚性圆盘，并在节点处耦合了具有刚度和阻尼特性的轴承支承。在计算过程中，程序构建了包含实部和虚部的状态向量及传递矩阵，能够有效求解含阻尼转子系统的复特征值问题。具体功能包括：计算系统的涡动频率、求解一阶及高阶临界转速、绘制坎贝尔图（Campbell Diagram）、计算对数衰减率以评估系统稳定性，以及模拟在不平衡力激励下的稳态不平衡响应。该项目填补了现有开源资源中缺乏含阻尼传递矩阵法完整实现的空白，适用于旋转机械（如汽轮机、航空发动机）的精细化动力学设计与振动故障分析。

详 情 说 明

含阻尼参数的转子动力学传递矩阵法仿真程序

项目简介

本项目是一个基于MATLAB开发的转子动力学分析系统。该系统针对旋转机械（如汽轮机、航空发动机）的转子系统进行精细化动力学建模。与传统仅考虑刚性支承或忽略阻尼的简化模型不同，本项目核心在于能够完整处理系统的阻尼参数（包括轴承阻尼和结构阻尼）。通过集总质量模型（Lumped Mass Model）将连续转子离散化，结合改进的传递矩阵法（TMM）和状态空间复特征值求解算法，实现了对转子系统临界转速、稳定性及不平衡响应的精确计算。

功能特性

• 精细化集总质量建模：支持定义分段轴段的几何参数（长、内外径）及材料属性（密度、弹性模量、剪切模量），并将转子离散为多个节点。
• 复杂边界条件支持：能够模拟附加圆盘（包含质量、极转动惯量、直径转动惯量）和各向异性轴承（包含$K_{xx}, K_{yy}$刚度及$C_{xx}, C_{yy}$阻尼参数）。
• 复特征值分析与稳定性评估：构建状态空间矩阵，求解含阻尼系统的复特征值，获取涡动频率（Whirl Frequency）及模态衰减特性。
• 坎贝尔图（Campbell Diagram）绘制：自动计算不同转速下的固有频率变化，用于识别临界转速。
• 系统稳定性量化：计算对数衰减率（Logarithmic Decrement），用于评估转子系统在特定转速下的稳定性。
• 不平衡响应分析：基于传递矩阵法（TMM）理论，模拟特定节点在不平衡力激励下的稳态响应（位移振幅与相位）。

系统要求

• MATLAB R2016b 或更高版本
• 无需额外工具箱（核心算法仅依赖MATLAB基础矩阵运算）

使用方法

1. 确保所有脚本文件位于MATLAB的当前工作路径中。
2. 直接运行主程序脚本（即包含main函数的脚本）。
3. 程序将依次执行参数定义、复特征值求解及不平衡响应计算，并在命令窗口输出计算进度。

代码逻辑与算法实现分析

本项目的主程序（main.m）逻辑严密，主要分为参数定义、复特征值稳定性分析、不平衡响应分析三个核心模块。

1. 系统参数定义与离散化

程序首先建立转子系统的物理模型：

• 基础参数：定义了材料密度、弹性模量（E）和剪切模量（G），这表明模型考虑了剪切变形的影响（如铁木辛柯梁理论基础）。
• 轴段模型：将1.0米长的转子轴划分为10个等长的轴段（每段0.1m），形成11个节点。
• 附加圆盘：在节点4和节点8处施加了集中质量和转动惯量，分别模拟不同规格的轮盘。
• 轴承支承：在节点2（前轴承）和节点10（后轴承）定义了弹性阻尼支承。代码明确设置了非对称的刚度系数（$K_{xx} neq K_{yy}$）和阻尼系数（$C_{xx}, C_{yy}$），体现了对各向异性及阻尼特性的完整支持。
• 不平衡量：在节点4和节点8定义了不平衡质量矩及其相位差（0度和90度），用于后续的响应分析。
• 转速设置：定义了从0到12000 RPM的扫描范围。

2. 核心求解 I：复特征值分析（坎贝尔图与稳定性）

为了克服传统传递矩阵法在求解特征值时可能出现的数值不稳定问题，本模块采用了全系统矩阵构建与状态空间求解相结合的混合策略：

• 全局矩阵构建：在每一转速步长下，程序调用辅助函数构建系统的全局质量矩阵（$M$）、刚度矩阵（$K$）、阻尼矩阵（$C$）以及陀螺矩阵（$G_{sys}$）。
• 陀螺效应处理：代码中通过计算运转换算阻尼矩阵 $D_{run} = C_{sys} + Omega G_{sys}$，显式地引入了与转速相关的陀螺力矩效应。
• 状态空间转换：为了求解二次特征值问题 $(M s^2 + D s + K)v = 0$，程序将其转换为一阶状态空间形式 $A z = lambda B z$。其中状态向量包含了位移和速度，从而能够同时求解特征值的实部（衰减率）和虚部（固有频率）。
• 模态筛选与排序：使用eig函数求解后，程序会自动过滤掉非物理模态（如静止或虚部为负的解），并按频率大小对前4阶模态进行排序。
• 稳定性指标：利用复特征值的实部和虚部计算对数衰减率，用于定量判断系统是否失稳。

3. 核心求解 II：不平衡稳态响应（传递矩阵法 TMM）

在计算稳态不平衡响应时，程序回归到经典的传递矩阵法（TMM）路径，以此验证TMM在非齐次方程求解中的优势：

• 算法选择：代码注释明确指出，此处使用严格的传递矩阵法（Transfer Matrix Method）进行计算。
• 状态向量定义：为了处理含阻尼及各向异性的情况，程序并未采用简化的4x4实数矩阵，而是采用了复数传递矩阵或增广传递矩阵（Augmented Matrix）的形式。这通常涉及构建包含位移、转角、弯矩、剪力以及非齐次项（不平衡力）的 $9 times 9$ 或复数域状态向量。
• 计算流程：程序初始化了针对节点4（响应节点）的振幅与相位存储数组，并设置了基于转速的循环。在循环内部，算法通过从转子一端推导至另一端的矩阵连乘操作，结合边界条件求解未知状态量，从而获得特定节点在不平衡激励下的稳态响应。

--- *注：本说明文档严格基于提供的代码片段内容编写，确保了对实际实现功能的真实描述。*