CN107977536B - 一种轮盘超转试验转接盘优化设计方法 - Google Patents

一种轮盘超转试验转接盘优化设计方法 Download PDF

Info

Publication number
CN107977536B
CN107977536B CN201711414531.5A CN201711414531A CN107977536B CN 107977536 B CN107977536 B CN 107977536B CN 201711414531 A CN201711414531 A CN 201711414531A CN 107977536 B CN107977536 B CN 107977536B
Authority
CN
China
Prior art keywords
disc
wheel disc
test
adapter
test wheel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711414531.5A
Other languages
English (en)
Other versions
CN107977536A (zh
Inventor
胡绪腾
金小杰
宋迎东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201711414531.5A priority Critical patent/CN107977536B/zh
Publication of CN107977536A publication Critical patent/CN107977536A/zh
Application granted granted Critical
Publication of CN107977536B publication Critical patent/CN107977536B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design

Abstract

本发明公开了一种轮盘超转试验转接盘优化设计方法，在保证试验轮盘定心效果情况下，通过优化参数轮缘长度，使转接盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移等于试验轮盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移。从而减少转接盘对试验轮盘径向位移的影响，使超转试验的结果更加准确，并提高超转试验成功率。

Description

一种轮盘超转试验转接盘优化设计方法
技术领域
本发明涉及一种航空发动机轮盘超转试验转接盘优化设计方法，属于发动机轮盘强度结构设计领域。
背景技术
轮盘可能因为在加速过程中燃油调节器失灵，发动机供油量增大，导致轮盘转速超过最大允许工作转速，而转子或轴发生故障可能也会让轮盘发生超转。带有加力燃烧室的军用发动机上最容易发生轮盘超转现象。所有发生在轮盘和轮缘的破裂均属于非包容性事故，这就意味着破裂产生的轮盘碎片在击穿发动机包容机匣后极有可能切断密集的油路系统，损毁控制元件，甚至打穿油箱或座舱，导致灾难性的后果。因此在轮盘结构设计中，必须进行预估和验证。
目前，国内外针对轮盘破裂转速的预估与验证得研究开展了大量工作。在进行轮盘超转试验时，试验轮盘需要通过转接盘与超转试验驱动轴相连，超转试验对定心的要求非常高。现有超转试验中，试验轮盘主要依靠短螺栓连接的定心。但很多轮盘连接处在接近破裂转速时会有较大的径向位移，这会将短螺栓剪断，导致超转试验失败。而采用航空发动机中依靠圆柱面抵紧的定心方式，让短螺栓只起到传扭作用，但这种定心方式下转接盘将影响试验轮盘径向位移，改变超转试验结果。
发明内容
为了解决现有转接盘会影响超转试验破裂结果的问题，本发明的目的是提供一种轮盘超转试验转接盘优化设计方法，以提高试验成功率，减少转接盘对超转试验结果的影响。
为实现上述目的，本发明采用的技术方案是：
一种轮盘超转试验转接盘优化设计方法，包括以下步骤：
步骤1：将转接盘的旋转对称剖面需设计成沿径向中心线轴对称的形式，将转接盘的轮缘部位厚度设计成比转接端厚度大；
步骤2：建立试验轮盘与转接盘二维有限元模型，计算破裂转速，并将过程写成APDL语言；
步骤3：计算试验轮盘破裂转速，并提取试验轮盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移；
步骤4：计算转接盘破裂转速；
步骤5：判断转接盘破裂转速是否大于试验轮盘破裂转速；若否，则增加转接盘轮毂、辐板、转接端厚度，并返回步骤4，直至转接盘破裂转速大于试验轮盘破裂转速；若是，则提取转接盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移，并进行步骤6；
步骤6：判断转接盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移是否大于试验轮盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移；若否，则进行步骤7；若是，则进行步骤8；
步骤7：增加转接盘二维有限元模型中轮缘长度，并返回步骤4；
步骤8：设置目标函数，优化参数为轮缘长度，利用有限元程序进行优化，得到最终优化后轮缘长度。
所述步骤8中，优化轮缘长度，使转接盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移等于试验轮盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移；优化完成后得到的轮缘长度值即为优化后转接盘轮缘长度。
有益效果：本发明的方法在保证试验轮盘定心效果情况下，通过优化参数轮缘长度L，使试验轮盘与转接盘连接处径向位移X1与X2在试验轮盘破裂转速N1时相等。从而减少转接盘对试验轮盘径向位移的影响，使超转试验的结果更加准确，并提高超转试验成功率。具有以下优点：
(1)保证试验轮盘定心效果，提高超转试验成功率。
(2)使超转试验的结果更加准确。
附图说明
图1为转接盘优化设计思路流程图；
图2a和2b为原型转接盘的结构示意图；其中，图2a为主视图，图2b为局部侧视图；
图3为旋转对称剖面中心对称示意图；
图4为试验机-转接盘-试验轮盘连接示意图；
图5为转接盘优化前后对比图；
图6为高压涡轮盘连接处，未优化转接盘连接处与优化转接盘连接处的转速-径向位移曲线图。
图中，1-试验轮盘；2-转接盘，201-轮缘，202-转接端，203-辐板，204-转接盘轮毂，205-螺栓孔；3-试验机驱动轴；A-转接盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点，B-试验轮盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点；L-轮缘长度，D1-轮缘厚度，D2-转接端厚度，D3-辐板厚度，D4-转接盘轮毂厚度。
具体实施方式
下面结合附图对本发明做更进一步的解释。
如图1所示，本发明的一种轮盘超转试验转接盘优化设计方法，包括以下步骤：
步骤1：将转接盘的旋转对称剖面需设计成沿其径向中心线对称的形式，可以保证转接盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点的转速-径向位移曲线和试验轮盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点的转速-径向位移曲线的拟合度更高；将轮缘部位厚度设计成比转接端厚度大，保证试验轮盘在旋转过程中良好定心。
如图4所示，试验轮盘1通过转接盘2连接至试验机驱动轴3。根据超转试验机型号，试验轮盘尺寸和以往设计经验，确定原型转接盘初始设计参数，如图1所示，该转接盘不能保证良好定心，在接近破裂时会因转接盘与试验轮盘径向位移不协调而剪断螺栓。增加轮缘部分，将转接盘旋转对称剖面设计成关于径向中心线轴对称形式，如图2所示，轮缘201初始长度L设置为10-20mm。轮缘厚度D1设计成比转接端厚度D2大10-30mm。
步骤2：建立试验轮盘与转接盘二维有限元模型，计算破裂转速，并将过程写成APDL语言；
通过试验获取试验轮盘和转接盘材料本构模型，建立多线性本构模型。建立试验盘与转接盘二维模型。设置单元行为方式为轴对称形式。对轮盘与转接盘施加轴向位移约束。根据经验施加全局旋转角速度，施加的角速度需大于所计算盘的破裂转速。采用大变形准静态分析，激活自动时间步中的弧长法选项，收敛准则设置为0.5。载荷子步数设置3000步(过小的步数会导致曲线震荡)，每5次记录一次结果。设定轮盘外缘的径向位移达到0.015mm时自动停止计算。并将上述设定全部导出为APDL语言。
步骤3：计算试验轮盘破裂转速N1，并提取试验轮盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移X1；如图4，该节点为B处；
通过有限元计算，提取试验轮盘破裂转速N1。提取B处坐标的节点在N1时的径向位移X1。
步骤4：计算转接盘破裂转速N2；
通过有限元计算，提取试验轮盘破裂转速N2。根据经验增加导出程序中转接盘模型的轮缘长度L，对于增加的L值不易过大。
步骤5：判断N2是否大于N1；
若否，则增加转接盘轮毂204的厚度D4、辐板203的厚度D3、转接端202的厚度D2各1-20mm，并返回步骤4，直至转接盘破裂转速大于试验轮盘破裂转速。若是，则提取转接盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移，该节点如图3中A所示，并进行步骤6；
本步骤能够保证试验时转速上升过程中试验轮盘先发生破裂。
步骤6：判断转接盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移X2是否大于试验轮盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移X1；
若否，则进行步骤7。若是，则进行步骤8；这样在步骤8的优化过程中，轮缘长度只会减小，保证转接盘的破裂转速不会因为步骤8而小于试验轮盘破裂转速。
步骤7：增加转接盘二维有限元模型中轮缘长度，并返回步骤4；
每次增加的转接盘模型的轮缘201长度为1-10mm，并返回步骤4计算。
增加轮缘长度，可以增加转接盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移。
步骤8：优化轮缘长度L至X1＝X2；
在有限元软件Ansys中设置目标函数为ABS(X2-X1)，优化参数为轮缘201长度L。步骤8中的轮缘201长度L初始值为增加过后，保证N2＞N1、X2≥X1的值。利用Ansys自带优化程序进行优化。优化完成后得到的L值即为优化后转接盘轮缘长度；通过减少转接盘在试验中对试验轮盘产生径向约束，减少对试验轮盘破裂转速试验结果的影响。
通过该优化方法，优化前后轮盘如图5所示。试验轮盘连接处，未优化转接盘连接处与优化转接盘连接处(图3中节点A处)转速-径向位移曲线如图6所示。优化转接盘与试验轮盘安装后破裂转速与试验盘破裂转速相差0.07％，未优化转接盘与试验轮盘安装后破裂转速与试验盘破裂转速相差1.01％。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以作出若干改进，这些改进也应视为本发明的保护范围。

Claims (5)

1.一种轮盘超转试验转接盘优化设计方法，其特征在于：包括以下步骤：
步骤1：将转接盘的旋转对称剖面需设计成沿径向中心线轴对称的形式，将转接盘的轮缘部位厚度设计成比转接端厚度大；
步骤2：建立试验轮盘与转接盘二维有限元模型，计算破裂转速，并将过程写成APDL语言；
步骤3：计算试验轮盘破裂转速，并提取试验轮盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移；
步骤4：计算转接盘破裂转速；
步骤5：判断转接盘破裂转速是否大于试验轮盘破裂转速；若否，则增加转接盘轮毂、辐板、转接端厚度，并返回步骤4，直至转接盘破裂转速大于试验轮盘破裂转速；若是，则提取转接盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移，并进行步骤6；
步骤6：判断转接盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移是否大于试验轮盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移；若否，则进行步骤7；若是，则进行步骤8；
步骤7：增加转接盘二维有限元模型中轮缘长度，并返回步骤4；
步骤8：设置目标函数，优化参数为轮缘长度，利用有限元程序进行优化，得到最终优化后轮缘长度。
2.根据权利要求1所述的轮盘超转试验转接盘优化设计方法，其特征在于：所述步骤2中，计算破裂转速为：对轮盘与转接盘施加轴向位移约束，施加全局旋转角速度，施加的角速度需大于所计算盘的破裂转速，采用大变形准静态分析，激活自动时间步中的弧长法选项，收敛准则设置为0.5，载荷子步数设置3000步，每5次记录一次结果，设定轮盘外缘的径向位移达到0.015mm时自动停止计算，并将上述设定全部导出为APDL语言。
3.根据权利要求1所述的轮盘超转试验转接盘优化设计方法，其特征在于：所述步骤5中，转接盘轮毂、辐板、转接端增加的厚度为1-20mm。
4.根据权利要求1所述的轮盘超转试验转接盘优化设计方法，其特征在于：所述步骤7中，每次增加转接盘二维有限元模型中轮缘长度为1-10mm。
5.根据权利要求1所述的轮盘超转试验转接盘优化设计方法，其特征在于：所述步骤8中，优化轮缘长度，使转接盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移等于试验轮盘二维有限元模型上试验轮盘与转接盘连接处螺栓孔圆心坐标的节点在试验轮盘破裂转速时的径向位移；优化完成后得到的轮缘长度值即为优化后转接盘轮缘长度。
CN201711414531.5A 2017-12-25 2017-12-25 一种轮盘超转试验转接盘优化设计方法 Active CN107977536B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711414531.5A CN107977536B (zh) 2017-12-25 2017-12-25 一种轮盘超转试验转接盘优化设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711414531.5A CN107977536B (zh) 2017-12-25 2017-12-25 一种轮盘超转试验转接盘优化设计方法

Publications (2)

Publication Number Publication Date
CN107977536A CN107977536A (zh) 2018-05-01
CN107977536B true CN107977536B (zh) 2020-02-07

Family

ID=62007460

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711414531.5A Active CN107977536B (zh) 2017-12-25 2017-12-25 一种轮盘超转试验转接盘优化设计方法

Country Status (1)

Country Link
CN (1) CN107977536B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111238804B (zh) * 2018-11-28 2021-08-03 中国航发商用航空发动机有限责任公司 航空发动机转子的超转破裂试验的试验件构型方法
CN112395799B (zh) * 2020-11-30 2023-03-10 中国电建集团成都勘测设计研究院有限公司 用于拱坝体形设计合理性及运行期工作状态的评价方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101853317A (zh) * 2010-04-20 2010-10-06 北京航空航天大学 一种涡轮盘结构概率设计系统的构建方法
CN103604587A (zh) * 2013-08-27 2014-02-26 杭州海罗马科技有限公司 高速旋转试验机及测试转子高速旋转强度性能的方法
CN105488252A (zh) * 2015-11-22 2016-04-13 沈阳黎明航空发动机(集团)有限责任公司 一种无螺栓结构涡轮盘挡板的破裂转速仿真分析方法
CN105512408A (zh) * 2015-12-14 2016-04-20 中国航空工业集团公司金城南京机电液压工程研究中心 盘片轴一体化超高转速涡轮转子三等分开槽预制方法
CN106446367A (zh) * 2016-09-09 2017-02-22 南京航空航天大学 基于弧长法非线性有限元分析的轮盘破裂转速预测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020141871A1 (en) * 2001-01-30 2002-10-03 Lakshimikantha Medamaranahally Axial-flow fan having inner and outer blades

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101853317A (zh) * 2010-04-20 2010-10-06 北京航空航天大学 一种涡轮盘结构概率设计系统的构建方法
CN103604587A (zh) * 2013-08-27 2014-02-26 杭州海罗马科技有限公司 高速旋转试验机及测试转子高速旋转强度性能的方法
CN105488252A (zh) * 2015-11-22 2016-04-13 沈阳黎明航空发动机(集团)有限责任公司 一种无螺栓结构涡轮盘挡板的破裂转速仿真分析方法
CN105512408A (zh) * 2015-12-14 2016-04-20 中国航空工业集团公司金城南京机电液压工程研究中心 盘片轴一体化超高转速涡轮转子三等分开槽预制方法
CN106446367A (zh) * 2016-09-09 2017-02-22 南京航空航天大学 基于弧长法非线性有限元分析的轮盘破裂转速预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
粉末冶金涡轮盘破裂转速分析与验证;秦仕勇 等;《燃气涡轮试验与研究》;20171031;第30卷(第5期);第36-41页 *

Also Published As

Publication number Publication date
CN107977536A (zh) 2018-05-01

Similar Documents

Publication Publication Date Title
CN107977536B (zh) 一种轮盘超转试验转接盘优化设计方法
EP2940283B1 (en) Turbomachine fan casing assembly
EP2952695B1 (en) Fan case
EP3438420B1 (en) Aero-engine and fusing method thereof under fan blade out load
CN102422139A (zh) 用于测试车辆的一个或多个制动元件的制动能力的方法和系统
CN105488252A (zh) 一种无螺栓结构涡轮盘挡板的破裂转速仿真分析方法
CN116050194B (zh) 一种确定涡轮转子无螺栓挡板径向配合紧度的方法
CN114361536B (zh) 燃料电池系统故障处理方法及装置
CN104071151A (zh) 一种飞车保护方法、装置和系统
CN103971760A (zh) 控制棒包壳管、控制棒组件及控制棒缓冲结构
CN105000173A (zh) 飞机刹车接地保护系统及方法
CN105512408A (zh) 盘片轴一体化超高转速涡轮转子三等分开槽预制方法
CN114611370A (zh) 预测转子超转破裂转速及破裂模式的方法及转子构型方法
CN108071429B (zh) 可失效转子支承结构及航空发动机
CN105756723B (zh) 一种主动熔断装置、方法及航空发动机
CN105523176A (zh) 一种飞机刹车动静态综合控制系统及方法
Shmotin et al. Numerical analysis of aircraft engine fan blade-out
CN102101278A (zh) 一种用于取出喷油器孔中衬套的工装装置
CN103808363A (zh) 一种机场道面检测用埋地传感器组件及实施方法
US8042406B2 (en) Method of determining overspeed reliability
EP3219965B1 (en) Takeoff power management system and method for gas turbine engines
CN203051521U (zh) 一种旋转盘轴的可靠连接结构
CN109606017A (zh) 一种带有保护装置的重型汽车轮毂
CN102319910A (zh) 伺服刀架用同步型过载保护结构
CN207059640U (zh) 一种半轴断裂保护装置及车辆

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant