CN110134990A - 航空发动机叶片前缘危险部位的检测方法 - Google Patents

航空发动机叶片前缘危险部位的检测方法 Download PDF

Info

Publication number
CN110134990A
CN110134990A CN201910264098.4A CN201910264098A CN110134990A CN 110134990 A CN110134990 A CN 110134990A CN 201910264098 A CN201910264098 A CN 201910264098A CN 110134990 A CN110134990 A CN 110134990A
Authority
CN
China
Prior art keywords
blade
stress
leading edge
aviation engine
finite element
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.)
Granted
Application number
CN201910264098.4A
Other languages
English (en)
Other versions
CN110134990B (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 CN201910264098.4A priority Critical patent/CN110134990B/zh
Publication of CN110134990A publication Critical patent/CN110134990A/zh
Application granted granted Critical
Publication of CN110134990B publication Critical patent/CN110134990B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • 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
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种航空发动机叶片前缘危险部位的检测方法。通过航空发动机工作状态试车中进行动应力测试、对叶片进行有限元模态计算、提取与动应力测试结果对应模态下的振动应力仿真结果等步骤，能够检测出航空发动机叶片前缘各个部位工作载荷下的应力水平，可以准确得到航空发动机叶片前缘危险部位及服役时可能的损伤尺寸，可以为航空发动机叶片抗外物损伤设计以及外场外物损伤修理原则制定提供数据支撑。

Description

航空发动机叶片前缘危险部位的检测方法
技术领域
本发明涉及航空发动机叶片检测领域。
背景技术
航空发动机在近地面服役过程中，会不可避免地吸入石子、金属碎片等硬物。这些外来硬物会以100m/s～350m/s的相对速度撞击航空发动机风扇或压气机叶片前缘，造成严重的冲击损伤。这种外来硬物的冲击损伤被称作“外物损伤”(简称FOD，Foreign ObjectDamage)。
受损伤的叶片疲劳强度会大大降低，在发动机服役过程中，会发生不可预估的疲劳断裂，造成严重的飞行事故。在现阶段的航空发动机叶片设计中，对外物损伤的损伤容限设计已经成为了不可缺少的一部分。航空发动机叶片的工作载荷是叶片高速旋转带来的离心力载荷产生的稳态应力以及叶片在旋转状态下的振动产生的振动应力的耦合作用。可定义叶片前缘的危险部位要求为前缘的高振动应力和高稳态应力的部位。由于叶片前缘在工作载荷下的应力分布不均匀，会使得前缘各部位的应力水平有很大的差异，当外物损伤发生在前缘的危险部位时，叶片极易发生疲劳破坏，需要针对该部位进行模拟外物损伤及高周疲劳试验来考核叶片的抗外物损伤性能。在制定外物损伤叶片的修理原则时，也需要对叶片前缘不同应力水平的部位进行划分，以制定出不同部位需要修理的损伤尺寸，对于前缘危险部位出现的损伤，应严格控制损伤的维修尺寸，较严重的损伤则应予以报废处理。可见无论是在叶片的设计阶段还是外场服役的维修阶段，都需要分析叶片前缘的应力水平，确定叶片前缘的危险部位。然而在现阶段仍缺少一种严格的航空发动机叶片前缘危险部位的检测方法。
发明内容
发明目的：本发明的目的在于提供一种可以准确检测到航空发动机叶片前缘危险部位及服役时可能的损伤尺寸的检测方法。
技术方案：本发明可采用以下技术方案：
一种航空发动机叶片前缘危险部位的检测方法，包括以下步骤：
(1)、在航空发动机叶片的测点通过测点应变片监测的应变数据，获取叶片工作状态下的实际固有频率fd，以及共振状态下的测点的真实振动应力σt；
(2)、使用ANSYS有限元分析软件，按照叶片工作状态下的边界条件及转速，建立叶片的有限元模型，并得到叶片的各阶动频，若计算的某阶动频结果与步骤(1)动应力测试实际测得的固有频率fd误差<5％，则可以确定该阶模态对应动应力测试得到的实际固有频率及动应力结果；并在ANSYS有限元分析软件后处理中提取叶片在该阶模态下的振动应力计算结果；
(3)、叶片前缘动应力计算：在(2)的模态应力结果中提取(1)动应力测试测点的振动应力σt0以及前缘部位各节点的振动应力σa0，通过公式
计算得到前缘部位各节点的真实振动应力σa；
(4)、按照与公式(2)中模态计算相同的边界条件和转速，使用ANSYS有限元分析软件对叶片的有限元模型进行静应力计算，在ANSYS有限元分析软件后处理中得到叶片在工作状态下的静应力结果，并提取叶片前缘部位各节点的静应力σm；
(5)、将步骤(4)中的静应力σm与步骤(3)中的振动应力σa转换为应力比为R＝-1的等效应力；
(6)、对比叶片前缘各个部位节点的等效应力，等效应力最大的节点部位即为该叶片前缘的危险部位。
进一步的，步骤(2)中，使用ANSYS有限元分析软件对航空发动机叶片进行模态计算。
进一步的，步骤(5)中，利用Goodman直线公式将前缘部位应力转换为应力比R＝-1时的等效应力σ-1，即
进一步的，步骤(6)中，使用LS-DYNA动力学仿真软件，对前缘危险部位以该位置对应线速度作为冲击速度，根据外场统计的外物尺寸，进行有限元数值仿真，模拟外物损伤过程，得到危险部位硬物冲击后的损伤尺寸。
进一步的，步骤(6)中，如该部位在模拟外物损伤后在发动机服役对应的载荷下进行规定循环的疲劳试验，仍未发生疲劳破坏，则说明该叶片的设计满足该型发动机抗外物损伤性能的需求。
有益效果：与现有技术相比，本发明能够检测航空发动机叶片前缘各个部位工作载荷下的应力水平，可以准确得到航空发动机叶片前缘危险部位及服役时可能的损伤尺寸，可以为航空发动机叶片抗外物损伤设计以及外场外物损伤修理原则制定提供数据支撑。
附图说明
图1是本发明检测方法的流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的详细说明。
结合图1所示，本发明提供了一种航空发动机叶片前缘危险部位的判定方法，具体实现步骤如下：
(1)叶片动应力测试：航空发动机叶片的动应力测试在发动机试车台进行，航空发动机按照工作转速试车过程中，通过测点应变片监测的应变数据，并对所采集的动态应变信号进行分析处理，得到叶片工作状态下的实际固有频率fd，以及共振状态下的测点的真实振动应力σt[1]；该步骤中通过测点应变片监测应变数据的方法可参考公开号为103954384A的中国发明专利申请具体实施方式部分。
(2)叶片有限元模态计算：按照叶片工作状态下的边界条件及转速，建立叶片的有限元模型，使用ANSYS有限元分析软件，按照叶片工作状态下的边界条件及转速，建立叶片的有限元模型，并得到叶片的各阶动频，若计算的某阶动频结果与步骤(1)动应力测试实际测得的固有频率fd误差<5％，则可以确定该阶模态对应动应力测试得到的实际固有频率及动应力结果；并在ANSYS有限元分析软件后处理中提取叶片在该阶模态下的振动应力计算结果。
(3)叶片前缘动应力计算：在(2)的模态应力结果中提取(1)动应力测试测点的振动应力σt0以及前缘部位各节点的振动应力σa0，通过公式1计算得到前缘部位各节点的真实振动应力σa。
(4)叶片前缘静应力计算：按照与(2)中模态计算相同的边界条件和转速，对叶片的有限元模型进行静应力计算，得到叶片在工作状态下的静应力结果，并提取叶片前缘部位各节点的静应力σm。
(5)计算应力比为R下的等效应力：为了比较前缘部位各节点的应力水平，综合考虑静应力与振动应力的影响，需要将步骤(4)中的静应力σm与步骤(3)中的振动应力σa转换为应力比为R的等效应力。由于实验室中在对叶片前缘危险部位进行高周疲劳试验时一般采用振动台来对其施加恒幅对称循环的疲劳载荷，可利用Goodman直线公式(2)将前缘部位应力转换为应力比R＝-1时的等效应力σ-1。
(6)对比叶片前缘各个部位节点的等效应力，等效应力最大的节点部位即为该叶片前缘的危险部位，如该部位被发动机吸入的外来硬物打伤，则在服役过程中极易出现疲劳破坏。故在对叶片进行抗外物损伤性能的考核试验时，该部位则应被作为考核对象。由于叶片前缘不同位置的特征尺寸(如前缘半径)不同，且不同位置的叶高的差别会导致外物冲击前缘时的速度不同。使用LS-DYNA等动力学仿真软件，对前缘危险部位以该位置对应线速度作为冲击速度，根据外场统计的外物尺寸，进行有限元数值仿真，模拟外物损伤过程，得到危险部位硬物冲击后的损伤尺寸。如该部位在模拟外物损伤后在发动机服役对应的载荷下进行规定循环的疲劳试验，仍未发生疲劳破坏，则说明该叶片的设计满足该型发动机抗外物损伤性能的需求。在制定外物损伤叶片的维修原则时，应严格控制危险部位损伤的维修尺寸，较严重的损伤则应予以报废处理。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下还可以作出若干改进，这些改进也应视为本发明的保护范围。

Claims (5)

1.一种航空发动机叶片前缘危险部位的检测方法，其特征在于，包括以下步骤：
(1)、在航空发动机叶片的测点通过测点应变片监测的应变数据，获取叶片工作状态下的实际固有频率fd，以及共振状态下的测点的真实振动应力σt；
(2)、使用ANSYS有限元分析软件，按照叶片工作状态下的边界条件及转速，建立叶片的有限元模型，并得到叶片的各阶动频，若计算的某阶动频结果与步骤(1)动应力测试实际测得的固有频率fd误差<5％，则可以确定该阶模态对应动应力测试得到的实际固有频率及动应力结果；并在ANSYS有限元分析软件后处理中提取叶片在该阶模态下的振动应力计算结果；
(3)、叶片前缘动应力计算：在(2)的模态应力结果中提取(1)动应力测试测点的振动应力σt0以及前缘部位各节点的振动应力σa0，通过公式
计算得到前缘部位各节点的真实振动应力σa；
(4)、按照与公式(2)中模态计算相同的边界条件和转速，使用ANSYS有限元分析软件对叶片的有限元模型进行静应力计算，在ANSYS有限元分析软件后处理中得到叶片在工作状态下的静应力结果，并提取叶片前缘部位各节点的静应力σm；
(5)、将步骤(4)中的静应力σm与步骤(3)中的振动应力σa转换为应力比为R＝-1的等效应力；
(6)、对比叶片前缘各个部位节点的等效应力，等效应力最大的节点部位即为该叶片前缘的危险部位。
2.根据权利要求1所述的航空发动机叶片前缘危险部位的检测方法，其特征在于：步骤(2)中，使用ANSYS有限元分析软件对航空发动机叶片进行模态计算。
3.根据权利要求1所述的航空发动机叶片前缘危险部位的检测方法，其特征在于：步骤(5)中，利用Goodman直线公式将前缘部位应力转换为应力比R＝-1时的等效应力σ-1，即
4.根据权利要求1所述的航空发动机叶片前缘危险部位的检测方法，其特征在于：步骤(6)中，使用LS-DYNA动力学仿真软件，对前缘危险部位以该位置对应线速度作为冲击速度，根据外场统计的外物尺寸，进行有限元数值仿真，模拟外物损伤过程，得到危险部位硬物冲击后的损伤尺寸。
5.根据权利要求5所述的航空发动机叶片前缘危险部位的检测方法，其特征在于：步骤(6)中，如该部位在模拟外物损伤后在发动机服役对应的载荷下进行规定循环的疲劳试验，仍未发生疲劳破坏，则说明该叶片的设计满足该型发动机抗外物损伤性能的需求。
CN201910264098.4A 2019-04-03 2019-04-03 航空发动机叶片前缘危险部位的检测方法 Active CN110134990B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910264098.4A CN110134990B (zh) 2019-04-03 2019-04-03 航空发动机叶片前缘危险部位的检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910264098.4A CN110134990B (zh) 2019-04-03 2019-04-03 航空发动机叶片前缘危险部位的检测方法

Publications (2)

Publication Number Publication Date
CN110134990A true CN110134990A (zh) 2019-08-16
CN110134990B CN110134990B (zh) 2020-09-25

Family

ID=67569200

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910264098.4A Active CN110134990B (zh) 2019-04-03 2019-04-03 航空发动机叶片前缘危险部位的检测方法

Country Status (1)

Country Link
CN (1) CN110134990B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110987389A (zh) * 2019-11-27 2020-04-10 南京航空航天大学 一种应力集中系数不小于3的外物损伤缺口标定方法
CN111723510A (zh) * 2020-06-28 2020-09-29 中国航发湖南动力机械研究所 叶片危险模态的识别方法
CN111850442A (zh) * 2020-07-06 2020-10-30 中国人民解放军空军工程大学 一种钛合金整体叶盘叶片预防高阶振型诱发叶尖掉块的强化方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060067830A1 (en) * 2004-09-29 2006-03-30 Wen Guo Method to restore an airfoil leading edge
CN107991080A (zh) * 2017-12-05 2018-05-04 中国人民解放军总参谋部第六十研究所 一种基于非接触测振与仿真计算的高频叶片模态分析方法
CN109211698A (zh) * 2018-03-21 2019-01-15 中国航空制造技术研究院 一种三点弯曲试样以及测试金属材料局部性能的方法
CN109374449A (zh) * 2018-09-25 2019-02-22 南京航空航天大学 一种考虑高、低周疲劳的叶片前后缘裂纹型硬物损伤可用极限确定方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060067830A1 (en) * 2004-09-29 2006-03-30 Wen Guo Method to restore an airfoil leading edge
CN107991080A (zh) * 2017-12-05 2018-05-04 中国人民解放军总参谋部第六十研究所 一种基于非接触测振与仿真计算的高频叶片模态分析方法
CN109211698A (zh) * 2018-03-21 2019-01-15 中国航空制造技术研究院 一种三点弯曲试样以及测试金属材料局部性能的方法
CN109374449A (zh) * 2018-09-25 2019-02-22 南京航空航天大学 一种考虑高、低周疲劳的叶片前后缘裂纹型硬物损伤可用极限确定方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HERBERT J SUTHERLAND等: "Effect of mean stress on the damage of wind turbine blades", 《JOURNAL OF SOLAR ENERGY ENGINEERING-TRANSACTIONS OF THE ASME》 *
孙护国等: "前缘半径对钛合金叶片抗外物损伤能力影响的数值分析", 《航空发动机》 *
岳鹏: "发动机涡轮叶片高低周复合疲劳寿命预测与可靠性分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
葛宁: "发动机叶片抗外物损伤能力评估技术研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110987389A (zh) * 2019-11-27 2020-04-10 南京航空航天大学 一种应力集中系数不小于3的外物损伤缺口标定方法
CN111723510A (zh) * 2020-06-28 2020-09-29 中国航发湖南动力机械研究所 叶片危险模态的识别方法
CN111723510B (zh) * 2020-06-28 2023-02-17 中国航发湖南动力机械研究所 叶片危险模态的识别方法
CN111850442A (zh) * 2020-07-06 2020-10-30 中国人民解放军空军工程大学 一种钛合金整体叶盘叶片预防高阶振型诱发叶尖掉块的强化方法
CN111850442B (zh) * 2020-07-06 2022-06-10 中国人民解放军空军工程大学 一种钛合金整体叶盘叶片预防高阶振型诱发叶尖掉块的强化方法

Also Published As

Publication number Publication date
CN110134990B (zh) 2020-09-25

Similar Documents

Publication Publication Date Title
CN109815521B (zh) 一种航空发动机叶片抗fod能力的评估方法
Zeng et al. Vibration response analysis of a cracked rotating compressor blade during run-up process
CN110134990A (zh) 航空发动机叶片前缘危险部位的检测方法
CN107084844B (zh) 一种叶片结构模拟件的试验方法
JP4202400B1 (ja) き裂進展予測方法及びプログラム
CN105352586B (zh) 通过检测共振频率实现叶片疲劳裂纹检测的方法
Yang et al. Nonlinear dynamic behavior of rotating blade with breathing crack
CN102087139B (zh) 汽轮发电机组低频振动频率成份实时分析方法
Xie et al. Comprehensive fatigue estimation and fault diagnosis based on Refined Generalized Multi-Scale Entropy method of centrifugal fan blades
CN115795744B (zh) 航空涡扇发动机零组件级低循环疲劳寿命载荷谱编制方法
CN111046490B (zh) 一种外物损伤缺口分析中的网格尺寸反演方法
Hanschke et al. The effect of foreign object damage on compressor blade high cycle fatigue strength
CN107976482A (zh) 在寿命考核试验中的盘类零件的裂纹在位监测系统和方法
CN112733398B (zh) 一种凹坑型硬物冲击损伤免修极限确定方法
CN114580094A (zh) 一种金属叶片旋转状态鸟撞损伤工程预估方法
Kim et al. Modal characteristics and fatigue strength of compressor blades
CN110083886A (zh) 一种基于有限元的直升机振动载荷仿真分析方法
Szczepanik et al. Crack initiation of rotor blades in the first stage of SO-3 compressor
CN110987388B (zh) 一种基于缺口疲劳强度的等效机械加工缺口的方法
CN112926698A (zh) 一种大型旋转装备振动预测与装配评价方法
Beck et al. Probabilistic mistuning assessment using nominal and geometry based mistuning methods
Kielb et al. Experimental and numerical study of forced response in a full-scale rotating turbine
CN114139276A (zh) 一种盘轴一体式整体叶盘结构疲劳寿命分析方法
CN110516270B (zh) 风力发电机中叶片连接螺栓的疲劳强度评估方法及系统
CN115391929B (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