多级转子装配优化ANSYS转子动力学仿真

关键词：ANSYS；转子动力学；装配优化；临界转速；遗传算法

一、文章简要介绍

航空发动机高压转子由多级转子装配而成，装配角度怎么排，直接决定整根转子在临界转速附近的振动大小。传统装配优化要么盯着几何偏心，要么盯着不平衡质量，要么只压工作转速下的振动，结果往往是临界转速过不去。这篇论文提出一个归一化目标函数，把一阶和二阶临界转速下六个节点的振动响应一起压下去，再用遗传算法找最优装配角度。四级转子的实验验证里，这个方法让一阶临界转速最大位移降了74.7%，二阶降了11.9%，两个临界转速同时降，综合效果超过四种传统方法。论文发表在Aerospace 2026年第13卷310篇，DOI：10.3390/aerospace13040310。

图1 四级转子有限元模型：前轴、压气机、高压涡轮、后轴四段装配，标注6个关键监测节点

图1是四级转子的有限元模型。转子按1比1简化为前轴、压气机、高压涡轮、后轴四段，材料是钢，弹性模量210吉帕，泊松比0.3。节点6和39是前后支撑边界，节点14、15、21、27分别对应压气机前后段、压气机出口和高压涡轮，这六个位置是临界转速穿越时振幅最高的地方。

二、仿真步骤

步骤一：建梁单元转子动力学模型。先在MATLAB里搭梁单元有限元模型，考虑陀螺效应，扫0到15000转每分钟的转速范围，画出Campbell图（图2）。前向涡动频率曲线和1倍频线的交点就是临界转速，算出来一阶7421转每分钟、二阶11810转每分钟，对应一阶和二阶弯曲模态。

图2 MATLAB梁单元模型Campbell图：考虑陀螺效应，前向涡动频率与1倍频交点给出临界转速7421与11810 rpm

图2是MATLAB梁单元模型的Campbell图。转速范围0到15000转每分钟，两个前向涡动分支和1倍频线的交点，就是两个临界转速。

步骤二：ANSYS Workbench建三维实体模型交叉验证。为了让临界转速算得踏实，论文在ANSYS Workbench 2022R2里重建了四级转子的三维实体模型，同样考虑陀螺效应做模态分析，转速范围一致。得到的Campbell图见图3，一阶临界转速7301转每分钟、二阶12045转每分钟，和MATLAB梁单元结果分别差1.6%和1.9%。两种建模方式，梁单元和三维实体单元，网格策略也不一样，这点差异在工程上完全可以接受，两个模型互相印证，可靠性就立住了。

图3 ANSYS Workbench三维实体模型Campbell图：临界转速7301与12045 rpm，与梁单元结果相差1.6%和1.9%

图3是ANSYS Workbench三维实体模型的Campbell图。实体模型考虑了完整的几何细节，临界转速落在7301和12045转每分钟，与梁单元模型高度一致。

步骤三：定目标函数。装配角是离散变量，压气机用12颗螺钉，可选的装配角是30度的整数倍；高压涡轮用24颗螺钉，可选角是15度的整数倍。以六个节点在一阶和二阶临界转速下的最大位移幅值为基础，构造归一化目标函数F，把两个临界转速的响应归一化后加起来，谁的综合振幅小，谁就是好装配。

步骤四：遗传算法求解。搜索空间里所有装配序列用遗传算法找最优解，100代终止，连续20代不改进也提前停。算法直接在全离散可行空间里搜索，不需要插值。另外论文还做了穷举遍历，把所有装配序列的目标函数值都算一遍，既验证了遗传算法的结果，也把最差装配序列找了出来。

三、关键结果与发现

第一个结果：归一化目标函数给出最优装配序列。以F为目标，遗传算法收敛到装配序列θz2等于150度、θz3等于180度、θz4等于150度，对应一阶临界转速最大位移0.0035毫米、二阶0.0091毫米。穷举搜索确认这是全局最优，最差装配序列是0度、30度、0度，一阶和二阶位移涨到0.0214和0.0261毫米，差了六倍以上。

图4 归一化目标函数优化结果：(a)遗传算法收敛曲线，(b)全部装配角度序列穷举搜索结果

图4是优化过程。左边是遗传算法收敛曲线，右边是穷举搜索的全景图，最优序列在曲面底部，一眼能看出来。

第二个结果：传统目标函数只压得住一头。仿真对比里，以几何偏心C为目标，二阶临界转速位移只降了0.03%，一阶反而升了4.4%。以综合不平衡U为目标，一阶降70.4%、二阶降59.3%。以双目标C与U为目标，降61.1%和55.3%。以工作转速17000转每分的振动为目标，降72.1%和52.0%。这些方法在单个指标上有表现，但都不是专门为跨临界转速设计的，归一化F在综合效果上最平衡。

图5 仿真对比：几何偏心C、综合不平衡U、双目标C与U、工作转速振动四种目标函数最优装配的临界转速位移

图5是四种传统目标函数和本文目标函数的仿真对比柱状图。一阶和二阶临界转速的位移幅值都列在一起，谁优谁劣很清楚。

第三个结果：实验验证，一阶降74.7%、二阶降11.9%。实验用真空高速动平衡机把转子从0升到15000转每分钟，按6个节点位置布置电涡流和速度传感器，八种目标函数各装一次、各测一轮升速。默认装配下实测临界转速是7000和12000转每分钟，和MATLAB、ANSYS两个仿真模型对得上，偏差在0.4%到6.0%之间。归一化F的实测结果是两个临界转速同时下降：一阶74.7%、二阶11.9%。八种目标函数里，只有F和工作转速f14000能做到两个临界转速同时降，而F的综合幅度最大。

图6 四级转子振动试验台：真空高速动平衡机驱动，电涡流与速度传感器按6节点布置

图6是实验台架。转子水平安装在真空高速动平衡机上，前后轴颈由摆式轴承座支撑，传感器位置和有限元模型的六个节点一一对应。

图7 实验对比：八种目标函数最优装配序列升速试验的一阶与二阶临界转速最大位移幅值

图7是八种目标函数的实验对比。灰色是默认装配的基准，F在两个临界转速上的柱子都明显矮下去，别的目标函数要么只矮一头，要么一头更低另一头更高。

第四个结果：柔性动挠度必须算。用不考虑柔性动挠度的模型算出来的F撇，一阶虽然降了72.5%，二阶反而升了34.2%。这说明装配优化如果只用静态误差传递模型，临界转速附近的动态响应根本管不住，柔性转子动力学模型是这套方法的地基。

四、我们提供的服务

我们可提供以下仿真服务：①ANSYS转子动力学分析（模态、Campbell图、临界转速、谐响应）；②多级转子装配误差传递建模与装配角度优化；③MATLAB转子动力学程序开发与交叉验证；④转子-轴承系统振动故障诊断与抑制方案。如果你的转子也有临界转速振动问题，欢迎私信聊聊。

如需相关服务，欢迎通过公众号“320科技工作室”联系我们。