变分损伤模型LS-DYNA子程序开发

关键词：ANSYS；转子动力学；装配优化；临界转速；遗传算法一、文章简要介绍航空发动机高压转子由多级转子装配而成，装配角度怎么排，直接决定整根转子在临界转速附近的振动大小。传统装配优化要么盯着几何偏心，要么盯着不平衡质量，要么只压工作转速下的振动，结果往往是临界转速过不去。这篇论文提出一个归一化目标函数，把一阶和二阶临界转速下六个节点的振动响应一起压下去，再用遗传算法找最优装配角度。四级转子的实验验证里，这个方法让一阶临界转速最大位移降了74.7%，二阶降了11.9%，两个临界转速同时降，综合效果超过四种传统方法。论文发表在Aerospace 2026年第13卷310篇，DOI：10.3390/aerospace13040310。图1是四级转子的有限元模型。转子按1比1简化为前轴、压气机、高压涡轮、后轴四段，材料是钢，弹性模量210吉帕，泊松比0.3。节点6和39是前后支撑边界，节点14、15、21、27分别对应压气机前后段、压气机出口和高压涡轮，这六个位置是临界转速穿越时振幅最高的地方。二、仿真步骤步骤一：建梁单元转子动力学模型。先在MATLAB里搭梁单元有限元模型，考虑陀螺效应，扫0到15000转每分钟的转速范围，画出Campbell图（图2）。前向涡动频率曲线和1倍频线的交点就是临界转速，算出来一阶7421转每分钟、二阶11810转每分钟，对应一阶和二阶弯曲模态。图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转每分钟，与梁单元模型高度一致。步骤三：定目标函数。装配角是离散变量，压气机用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是优化过程。左边是遗传算法收敛曲线，右边是穷举搜索的全景图，最优序列在曲面底部，一眼能看出来。第二个结果：传统目标函数只压得住一头。仿真对比里，以几何偏心C为目标，二阶临界转速位移只降了0.03%，一阶反而升了4.4%。以综合不平衡U为目标，一阶降70.4%、二阶降59.3%。以双目标C与U为目标，降61.1%和55.3%。以工作转速17000转每分的振动为目标，降72.1%和52.0%。这些方法在单个指标上有表现，但都不是专门为跨临界转速设计的，归一化F在综合效果上最平衡。图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是实验台架。转子水平安装在真空高速动平衡机上，前后轴颈由摆式轴承座支撑，传感器位置和有限元模型的六个节点一一对应。图7是八种目标函数的实验对比。灰色是默认装配的基准，F在两个临界转速上的柱子都明显矮下去，别的目标函数要么只矮一头，要么一头更低另一头更高。第四个结果：柔性动挠度必须算。用不考虑柔性动挠度的模型算出来的F撇，一阶虽然降了72.5%，二阶反而升了34.2%。这说明装配优化如果只用静态误差传递模型，临界转速附近的动态响应根本管不住，柔性转子动力学模型是这套方法的地基。来源：320科技工作室