预测转子超转破裂转速及破裂模式的方法及转子构型方法技术

本发明专利技术公开了一种预测转子超转破裂转速及破裂模式的方法及转子构型方法，预测转子超转破裂转速及破裂模式的方法包括以下步骤：S1、建立转子的二维轴对称有限元模型及轮盘鼓筒单元的三维有限元模型；S2、对二维轴对称有限元模型进行静力学参数设置和计算分析，获得不同轮盘鼓筒单元的三维有限元模型的稳态的初始应力场和初始位移场；S3、将获得的初始应力场和初始位移场导入对应的轮盘鼓筒单元的三维有限元模型中进行初始工况仿真，获得对应的轮盘鼓筒单元的初始应力状态和位移状态；S4、向步骤S3获得的三维有限元模型中施加随时间变化的转速谱，设置动力学参数，计算后进行数据后处理，获得对应轮盘的破裂转速及对应的破裂模式。破裂模式。破裂模式。

【技术实现步骤摘要】
预测转子超转破裂转速及破裂模式的方法及转子构型方法


[0001]本专利技术涉及一种预测转子超转破裂转速及破裂模式的方法及转子构型方法。

技术介绍

[0002]航空发动机及燃气轮机的涡轮、压气机、风扇和涡轮增压器转子必须具有足够的强度以便能承受相应的超转及破裂转速，适航条款CCAR
‑
33
‑
R2对航空发动机转子的超转及破裂安全性提出了明确且清晰的要求，是航空发动机适航取证必须满足的强度要求。航空史中已发生很多起因转子超转破裂引起机匣包容失效等进而导致重大人员伤亡的航空事故，因此在设计及试验验证阶段准确预测转子的超转破裂转速及破裂模式就显得尤为重要。其中，转子是指由轴承支撑的旋转体，包括轮盘和鼓筒；超转破裂转速指转子的转速超过最高转速直至轮盘破裂时对应的转速；破裂模式指轮盘破裂的方式，包括径向破裂和周向破裂两种。
[0003]现有超转破裂转速计算方法一般为对简化的二维有限元转子模型或者简化的单盘三维有限元模型进行静力学分析，然后基于平均应力法得到结果。实践中发现采用现有静力学计算方法，破裂转速的数值严重依赖于修正系数的选取，难以准确预测破裂转速，更无法获得破裂模式，这给转子部件、核心机及整机试验造成不利影响，给后续航空发动机或燃气轮机产品的服役带来一定的风险。

技术实现思路

[0004]本专利技术要解决的技术问题是现有技术中转子的超转破裂转速计算方法无法准确预测破裂转速且无法获得破裂模式的缺陷，提供一种预测转子超转破裂转速及破裂模式的方法及转子构型方法。/>[0005]本专利技术是通过下述技术方案来解决上述技术问题：本专利技术提供了一种预测转子超转破裂转速及破裂模式的方法，所述转子包括多个轮盘和鼓筒，所述鼓筒的两端分别连接于相邻两个所述轮盘，所述方法包括以下步骤：S1、建立所述转子的二维轴对称有限元模型及轮盘鼓筒单元的三维有限元模型，其中，每个所述轮盘鼓筒单元包括一个轮盘和与所述轮盘连接的鼓筒的其中一半的部分；S2、对所述转子的二维轴对称有限元模型进行静力学参数设置和计算分析，获得所述转子中不同所述轮盘鼓筒单元的三维有限元模型的稳态的初始应力场和初始位移场；S3、将步骤S2中获得的各个所述轮盘鼓筒单元的初始应力场和初始位移场导入步骤S1中建立的对应的所述轮盘鼓筒单元的三维有限元模型中，进行初始工况仿真，获得对应的所述轮盘鼓筒单元的初始应力状态和位移状态；S4、向步骤S3获得的所述轮盘鼓筒单元的三维有限元模型中施加随时间变化的转速谱，设置动力学参数，计算后进行数据后处理，获得对应轮盘的破裂转速及对应的破裂模式。
[0006]在本方案中，通过上述方法克服了现有技术对因模型简化带来的计算风险，提高
了转子超转破裂转速的预测精度，为航空发动机及燃气轮机转子的超转破裂设计提供更准确的依据和指导，为适航取证提供有力的支撑。
[0007]较佳地，在步骤S1中，所述转子的二维轴对称有限元模型通过将转子的轴对称部分采用轴对称单元进行网格离散，其中少量的非轴对称结构如孔及螺栓采用带厚度的平面应力单元进行离散，转子叶片通过等效质量单元离散以考虑转子叶片所施加给所述转子的离心力。
[0008]在本方案中，采用上述方法建模可以得到较为准确的模型，从而得到较为准确的轮盘鼓筒单元的三维有限元模型的初始应力场和初始位移场，进而获得更加准确的轮盘的破裂转速及破裂模式。
[0009]较佳地，在步骤S1中，所述轮盘鼓筒单元的三维有限元模型建立的单元类型为六面体单元，在所有相邻所述六面体单元之间建立无厚度的结合力单元，其中，所述结合力单元用以模拟轮盘破裂的萌生、扩展进而至所述轮盘的破裂。
[0010]在本方案中，采用上述方法建模可以进一步得到较为准确的模型，从而得到更加准确的轮盘鼓筒单元的三维有限元模型的初始应力场和初始位移场，进而获得更加准确的轮盘的破裂转速及破裂模式。
[0011]较佳地，在步骤S4中，通过数据后处理确定初始失效的结合力单元，即破裂裂纹萌生的位置，其对应的转速即为破裂转速；继续增加转速后，轮盘和鼓筒会出现破裂，根据破裂时轮盘的形态，即可判断出轮盘的破裂模式。
[0012]较佳地，在步骤S2中，所述静力学参数包括接触参数、材料参数、边界条件和载荷。
[0013]在本方案中，通过设置多种静力学参数获得轮盘鼓筒单元的三维有限元模型的稳态的初始应力场和初始位移场。
[0014]较佳地，所述接触参数包括接触刚度参数，所述接触刚度参数取值为0.5
‑
1。
[0015]在本方案中，由于转子系统中存在大量的非线性接触，接触刚度参数不宜过大或过小。
[0016]较佳地，在步骤S2中，在计算分析求解时要考虑几何大变形，并通过设置多个载荷以保证收敛性，然后采用基于图理论的分布式并行求解算法进行求解，获得所述初始应力场和初始位移场。
[0017]较佳地，采用基于图理论的分布式并行求解算法进行求解时，收敛设置采用力收敛，收敛值设为100。
[0018]较佳地，在步骤S3中，在所述轮盘鼓筒单元的三维有限元模型的内部单元导入所述初始应力场，在所述轮盘鼓筒单元的三维有限元模型的边界施加所述初始位移场，然后进行准静态求解，获得所述轮盘鼓筒单元的初始应力状态和位移状态。
[0019]较佳地，进行准静态求解时采用动力松弛法，其中，所述动力松弛法的收敛准则设置为0.001。
[0020]较佳地，在步骤S4中，应用非线性显式动力学方法进行并行计算以获得相应的待处理数据。
[0021]较佳地，在步骤S4中，转速谱采用线性增加且加速率为8
‑
10m/s的转速谱，所述转速谱的起始转速为步骤S2中进行稳态计算对应的转子转速。
[0022]本专利技术还提供一种转子构型方法，包括以下步骤：
S10、进行转子的初始结构设计；S20、采用上述的预测转子超转破裂转速及破裂模式的方法获得初始结构设计的转子中的轮盘破裂的转速及对应的破裂模式，并对所述转子中的轮盘超转破裂的安全裕度进行评价，判断所述转子是否满足裕度要求；如不满足裕度要求，则执行步骤S30，如满足裕度要求，则执行步骤S40；S30、重新进行转子的结构设计，然后继续执行步骤S20；S40、确定转子的构型，以待后续加工。
[0023]本专利技术的积极进步效果在于：本专利技术的预测转子超转破裂转速及破裂模式的方法克服了现有技术对因模型简化带来的计算风险，提高了转子超转破裂转速的预测精度，为航空发动机及燃气轮机转子的超转破裂设计提供更准确的依据和指导，为适航取证提供有力的支撑。
附图说明
[0024]图1为本专利技术实施例中预测转子超转破裂转速及破裂模式的方法的流程图。
[0025]图2为转子的二维轴对称有限元模型的局部示意图。
[0026]图3为单个轮盘鼓筒单元的三维有限元模型的局部扇区模型。
[0027]图4为本专利技术实施例中提供的结合力单元示意图。
[0028]图5为本专利技术较佳实施例中转速谱示意图。<本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种预测转子超转破裂转速及破裂模式的方法，所述转子包括多个轮盘和鼓筒，所述鼓筒的两端分别连接于相邻两个所述轮盘，其特征在于，所述方法包括以下步骤：S1、建立所述转子的二维轴对称有限元模型及轮盘鼓筒单元的三维有限元模型，其中，每个所述轮盘鼓筒单元包括一个轮盘和与所述轮盘连接的鼓筒的其中一半的部分；S2、对所述转子的二维轴对称有限元模型进行静力学参数设置和计算分析，获得所述转子中不同所述轮盘鼓筒单元的三维有限元模型的稳态的初始应力场和初始位移场；S3、将步骤S2中获得的各个所述轮盘鼓筒单元的初始应力场和初始位移场导入步骤S1中建立的对应的所述轮盘鼓筒单元的三维有限元模型中，进行初始工况仿真，获得对应的所述轮盘鼓筒单元的初始应力状态和位移状态；S4、向步骤S3获得的所述轮盘鼓筒单元的三维有限元模型中施加随时间变化的转速谱，设置动力学参数，计算后进行数据后处理，获得对应轮盘的破裂转速及对应的破裂模式。2.如权利要求1所述的预测转子超转破裂转速及破裂模式的方法，其特征在于，在步骤S1中，所述转子的二维轴对称有限元模型通过将转子的轴对称部分采用轴对称单元进行网格离散，其中少量的非轴对称结构如孔及螺栓采用带厚度的平面应力单元进行离散，转子叶片通过等效质量单元离散以考虑转子叶片所施加给所述转子的离心力。3.如权利要求2所述的预测转子超转破裂转速及破裂模式的方法，其特征在于，在步骤S1中，所述轮盘鼓筒单元的三维有限元模型建立的单元类型为六面体单元，在所有相邻所述六面体单元之间建立无厚度的结合力单元，其中，所述结合力单元用以模拟轮盘破裂的萌生、扩展进而至所述轮盘的破裂。4.如权利要求3所述的预测转子超转破裂转速及破裂模式的方法，其特征在于，在步骤S4中，通过数据后处理确定初始失效的结合力单元，即破裂裂纹萌生的位置，其对应的转速即为破裂转速；继续增加转速后，轮盘和鼓筒会出现破裂，根据破裂时轮盘的形态，即可判断出轮盘的破裂模式。5.如权利要求1所述的预测转子超转破裂转速及破裂模式的方法，其特征在于，在步骤S2中，所述静力学参数包括接触参数、材料参数、边界条件和载荷。6.如权利要求5所述的预测转子超转破裂转速及破裂...

【专利技术属性】
技术研发人员：王小庆，刘文，张辉，魏铭瑛，
申请(专利权)人：中国航发商用航空发动机有限责任公司，
类型：发明
国别省市：