CN107122531B - 一种基于加速寿命试验的高速叶轮寿命评估方法 - Google Patents

一种基于加速寿命试验的高速叶轮寿命评估方法 Download PDF

Info

Publication number
CN107122531B
CN107122531B CN201710240350.9A CN201710240350A CN107122531B CN 107122531 B CN107122531 B CN 107122531B CN 201710240350 A CN201710240350 A CN 201710240350A CN 107122531 B CN107122531 B CN 107122531B
Authority
CN
China
Prior art keywords
speed impeller
speed
life
service life
impeller
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.)
Expired - Fee Related
Application number
CN201710240350.9A
Other languages
English (en)
Other versions
CN107122531A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201710240350.9A priority Critical patent/CN107122531B/zh
Publication of CN107122531A publication Critical patent/CN107122531A/zh
Application granted granted Critical
Publication of CN107122531B publication Critical patent/CN107122531B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04Ageing analysis or optimisation against ageing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)

Abstract

本发明公开了一种基于加速寿命试验的高速叶轮寿命评估方法,以涡轮机高速叶轮在加速寿命试验下的寿命预测问题作为切入点,分析小样本数据下试验个数是否满足观察要求,构建了高速叶轮在机械应力作用下寿命特征与应力的关系框架,推算出其寿命评估模型,并详细的说明了各个关键步骤逐步实施的策略,着重展开说明了如何利用威布尔分布描述其寿命分布并得到模型关键参数的估计值的策略。通过在所获得模型中关键参数估计值的基础上结合寿命分布模型得到高速叶轮在给定可靠度和置信水平下的工作寿命置信下限模型,最后结合高速叶轮的实际工作频率等现场数据,推算出高速叶轮的使用寿命年限。

Description

一种基于加速寿命试验的高速叶轮寿命评估方法
技术领域
本发明属于可靠性分析寿命评估技术领域,具体涉及一种基于加速寿命试验的高速叶轮寿命评估方法的设计。
背景技术
高速叶轮是离心压缩机,空气涡轮机中的最关键旋转部件,在船舶工程、电气工程、机械工程中,都得到了广泛应用,在军用、民用、家用机械或电器的各个方面,也都发挥着极其重要的作用。作为主要承力部件、叶轮的主要功能是传输功率,工作时承受离心力、气动力、激振力、外物冲击等循环交变载荷与动载荷作用,其主要失效模式是循环疲劳断裂、塑性变形超限、振动开裂以及这些模式之间相互作用的结果。在外部载荷的作用下,叶轮的材料性能或结构有可能会产生不可逆的变化。在十分严酷恶劣的工作环境下,叶片除受到高速旋转带来的离心载荷的作用外,还受到复杂的气动力和热应力的影响,以及振动、腐蚀和氧化等作用,因此叶片故障时有发生,其性能的优劣及寿命的长短直接影响到发动机的工作效率和运行可靠性,且一旦发生故障,就有很大可能造成毁灭性的后果。众所周知,塑性变形在工程结构中是广泛存在的现象,由塑性变形造成的构件破坏,会造成巨大的经济损失甚至危及生命。
作为空气涡轮机的核心部件,涡轮机对气体做功、提高气体能量的唯一部件,高速叶轮性能的好坏直接关系到涡轮机的性能。高速叶轮的变形损坏问题一直是影响涡轮机可靠性的主要因素之一。在涡轮机中,叶轮通常承受很大的应力,工作时受力极其复杂,主要包括旋转产生的惯性力、气动载荷和热载荷。由于气流分布不均勾、不稳定及转子的不平衡等均会使叶轮发生振动,导致叶轮在工作中可靠性差,出现塑性变形,甚至产生疲劳断裂。随着近年来科学技术的发展,叶轮的设计取得了很大的进步,按照传统的寿命试验技术进行评估很难发现失效。涡轮机高速叶轮系统的失效模式多样,在所有失效模式动态参数条件下对其进行可靠性及寿命评估存在较大困难,现有的故障树分析方法很难表述复杂的高速叶轮系统。
发明内容
本发明的目的是为了解决现有的可靠性寿命评估方法在分析复杂、动态、失效模式多样的系统时表达复杂、计算效率低、求解困难的问题,提出了一种基于加速寿命试验的高速叶轮寿命评估方法。
本发明的技术方案为:一种基于加速寿命试验的高速叶轮寿命评估方法,包括以下步骤:
S1、对高速叶轮进行加速寿命试验,得到每次试验下高速叶轮的相关试验数据;
S2、根据《材料疲劳试验统计分析方法》,验证中值疲劳寿命测定中的观察样本个数是否符合要求;
S3、根据涡轮机高速叶轮系统试验的历史统计数据,并结合高速叶轮加速寿命试验数据,通过威布尔分布参数的矩估计求解出高速叶轮的形状参数;
S4、用逆幂率模型描述高速叶轮在机械应力作用下寿命特征与应力的关系,推算出高速叶轮在加速寿命试验中的寿命评估模型;
S5、将步骤S3求解出的形状参数以及步骤S1中的相关试验数据代入寿命评估模型中进行计算,得到高速叶轮的使用寿命下限;
S6、根据高速叶轮的实际工作现场数据,推算出高速叶轮的使用寿命年限。
本发明的有益效果是:
(1)本发明通过将高速叶轮的历史数据与失效模式相结合,运用威布尔分布参数的矩估计,通过对试验对象的可靠性数据分析,建立了高速叶轮基于加速实验的寿命评估模型,可根据系统的使用频率,工作状况直观得出系统寿命下限。本发明克服了传统可靠性寿命评估方法在分析复杂、动态、失效模式多样的系统时表达复杂、计算效率低、求解困难等问题。
(2)本发明将加速寿命试验引入高速叶轮系统的寿命评估方法中,通过对叶轮的威布尔分布形状参数的估计,结合叶轮的现场试验寿命数据建立寿命模型进行评估,极大缩短传统可靠性试验的持续时间,只需记录加速寿命试验的现场失效数据,就能得出高速叶轮系统的寿命年限,表述方式更加方便简洁,无需分析所有失效模式动态参数,提高了计算效率,因此能显著地提高空气涡轮机高速叶轮系统的可靠性分析和寿命评估效率。
附图说明
图1为本发明提供的一种基于加速寿命试验的高速叶轮寿命评估方法流程图。
图2为本发明实施例的高速叶轮加速寿命试验原理示意图。
具体实施方式
下面结合附图对本发明的实施例作进一步的说明。
本发明提供了一种基于加速寿命试验的高速叶轮寿命评估方法,如图1所示,包括以下步骤:
S1、对高速叶轮进行加速寿命试验,起始叶轮以正常转速工作,每工作一定时间停机一段时间,然后提高转速进行下一次试验,直到高速叶轮在极限转速下发生形变失效时停止试验。记下加速寿命试验中高速叶轮的转速作为观察样本,其中正常工作转速记为S0,正常工作持续时间记为t0,最高无失效转速记为S1,最高无失效转速持续时间记为t1,极限转速记为Sm,极限转速持续时间记为tm
高速叶轮的加速寿命试验原理如图2所示,其中,横轴为使用寿命(对数化),竖轴为应力水平。t0'为高速叶轮在正常应力水平(工作状态)S0下的真实寿命,t1为叶轮样本在加速应力水平S1下的试验寿命。本发明实施例中将试验所得数据t1换算到正常应力水平S0下的使用寿命t0。通过加速寿命试验使t0尽可能向真实寿命t0'逼近,从而得到更准确的高速叶轮寿命评估结果。
S2、根据《材料疲劳试验统计分析方法》,验证中值疲劳寿命测定中的观察样本个数是否符合要求。该步骤具体包括以下分步骤:
S21、由高速叶轮系统历史统计数据样本,结合使用频率和设计寿命,取ni(i≥3)个样本的寿命循环次数,取其对数,得到高速叶轮寿命循环次数对数,如表1所示。
表1
Figure BDA0001269211550000031
根据公式:
Figure BDA0001269211550000032
Figure BDA0001269211550000033
计算子样平均值
Figure BDA0001269211550000036
和子样标准差s,得出样本的变异系数
Figure BDA0001269211550000034
S22:根据HBZ 112-1986《材料疲劳试验统计分析方法》中值疲劳寿命的测定,验证观察样本个数是否满足下式:
Figure BDA0001269211550000035
其中δmax为误差限度,本发明实施例中取5%,γ为置信水平,根据HBZ 112-1986《材料疲劳试验统计分析方法》中tγ数值表与观察样本个数N,进行三次多项式拟合,得到小样本修正下的tγ数值表,如表2所示。
表2
tr数值表(γ=95%)
Figure BDA0001269211550000041
tr数值表(γ=90%)
Figure BDA0001269211550000042
根据公式(3)检验观察样本个数是否满足要求。
S3、根据涡轮机高速叶轮系统试验的历史统计数据,并结合高速叶轮加速寿命试验数据,通过威布尔分布参数的矩估计求解出高速叶轮的形状参数。
把采集到的n个正常工作及加速寿命试验下高速叶轮转速的样本,记为X1,X2,...,Xn,由矩估计
Figure BDA0001269211550000043
得形状参数β的估计为:
Figure BDA0001269211550000044
式中
Figure BDA0001269211550000045
表示通过威布尔分布参数的矩估计求出的形状参数β的估计值,
Figure BDA0001269211550000046
λ为任意常数,本发明实施例中λ取为1。
S4、用逆幂率模型描述高速叶轮在机械应力作用下寿命特征与应力的关系,推算出高速叶轮在加速寿命试验中的寿命评估模型。
根据高速叶轮在机械应力作用下寿命特征与应力的关系以及加速寿命试验理论,由对数线性化的逆幂率模型:
lnξ=a+b ln S (5)
推导出正常应力水平下的等效寿命t:
Figure BDA0001269211550000051
结合威布尔分布,给出高速叶轮在正常应力水平下,可靠度为R,置信水平为γ的寿命置信下限tR,low的评估模型:
Figure BDA0001269211550000052
S5、由威布尔分布描述的寿命评估模型,将步骤S3求解出的形状参数β以及步骤S1中的相关试验数据代入公式(7)进行计算,得到高速叶轮在给定可靠度R,置信水平γ下的使用寿命下限T=tR,low
S6、根据高速叶轮的实际工作现场数据,推算出高速叶轮的使用寿命年限。
根据高速叶轮实际平均年化工作频率f′,实际正常工作转速S′以及实际正常工作持续时间t′,给出高速叶轮基于加速寿命试验的年化寿命公式:
T=S′×t′×f′×N′ (8)
由公式(8)可以计算得到高速叶轮在加速寿命试验下的使用寿命年限
Figure BDA0001269211550000053
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。

Claims (7)

1.一种基于加速寿命试验的高速叶轮寿命评估方法,其特征在于,包括以下步骤:
S1、对高速叶轮进行加速寿命试验,得到每次试验下高速叶轮的相关试验数据;
S2、验证中值疲劳寿命测定中的观察样本个数是否符合要求,具体包括以下分步骤:
S21、由高速叶轮系统历史统计数据样本,结合使用频率和设计寿命,取ni(i≥3)个样本的寿命循环次数,取其对数,得到高速叶轮寿命循环次数对数;
S22:根据HBZ 112-1986《材料疲劳试验统计分析方法》中值疲劳寿命的测定,确定观察样本个数是否符合要求;
S3、根据涡轮机高速叶轮系统试验的历史统计数据,并结合高速叶轮加速寿命试验数据,通过威布尔分布参数的矩估计求解出高速叶轮的形状参数;
S4、用逆幂率模型描述高速叶轮在机械应力作用下寿命特征与应力的关系,推算出高速叶轮在加速寿命试验中的寿命评估模型;
S5、将步骤S3求解出的形状参数以及步骤S1中的相关试验数据代入寿命评估模型中进行计算,得到高速叶轮的使用寿命下限;
S6、根据高速叶轮的实际工作现场数据,推算出高速叶轮的使用寿命年限。
2.根据权利要求1所述的高速叶轮寿命评估方法,其特征在于,所述步骤S1中的相关试验数据包括正常工作转速S0,正常工作持续时间t0,最高无失效转速S1,最高无失效转速持续时间t1,极限转速Sm,极限转速持续时间tm
3.根据权利要求2所述的高速叶轮寿命评估方法,其特征在于,所述步骤S3具体为:
把采集到的n个正常工作及加速寿命试验下高速叶轮转速的样本,记为X1,X2,...,Xn,由矩估计
Figure FDA0002461288600000011
得形状参数β的估计为:
Figure FDA0002461288600000012
式中
Figure FDA0002461288600000013
表示通过威布尔分布参数的矩估计求出的形状参数β的估计值,
Figure FDA0002461288600000014
λ为任意常数。
4.根据权利要求3所述的高速叶轮寿命评估方法,其特征在于,所述步骤S4中的寿命评估模型具体为:
Figure FDA0002461288600000021
式中R为可靠度,γ为置信水平。
5.根据权利要求4所述的高速叶轮寿命评估方法,其特征在于,所述步骤S5中的使用寿命下限T=tR,low
6.根据权利要求5所述的高速叶轮寿命评估方法,其特征在于,所述步骤S6中高速叶轮的实际工作现场数据包括实际平均年化工作频率f′,实际正常工作转速S′以及实际正常工作持续时间t′。
7.根据权利要求6所述的高速叶轮寿命评估方法,其特征在于,所述步骤S6中高速叶轮的使用寿命年限
Figure FDA0002461288600000022
CN201710240350.9A 2017-04-13 2017-04-13 一种基于加速寿命试验的高速叶轮寿命评估方法 Expired - Fee Related CN107122531B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710240350.9A CN107122531B (zh) 2017-04-13 2017-04-13 一种基于加速寿命试验的高速叶轮寿命评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710240350.9A CN107122531B (zh) 2017-04-13 2017-04-13 一种基于加速寿命试验的高速叶轮寿命评估方法

Publications (2)

Publication Number Publication Date
CN107122531A CN107122531A (zh) 2017-09-01
CN107122531B true CN107122531B (zh) 2020-06-09

Family

ID=59724678

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710240350.9A Expired - Fee Related CN107122531B (zh) 2017-04-13 2017-04-13 一种基于加速寿命试验的高速叶轮寿命评估方法

Country Status (1)

Country Link
CN (1) CN107122531B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109682580B (zh) * 2017-10-18 2020-07-28 宁波方太厨具有限公司 油烟机叶轮保养期限判定测试方法
CN108776730A (zh) * 2018-05-30 2018-11-09 中国航发动力股份有限公司 一种燃气轮机叶片断裂故障快速定位方法
CN109987251B (zh) * 2019-04-08 2020-08-28 中国航空综合技术研究所 基于应力等效方法的针对薄弱环节的等效寿命试验方法及设备
CN110472275B (zh) * 2019-07-03 2023-05-30 上海卫星工程研究所 飞轮产品不同工况下数据处理与寿命评估方法
CN110260907B (zh) * 2019-07-11 2021-07-16 中国电子科技集团公司第四十九研究所 一种用于传感器的温度应力无失效加速寿命试验方法
CN110375971B (zh) * 2019-07-11 2020-08-14 北京理工大学 径流式涡轮叶轮和压气机叶轮的加速寿命试验装置及方法
CN112444725B (zh) * 2019-09-05 2024-01-26 中车株洲电力机车研究所有限公司 一种通孔焊点对比加速寿命试验方法
CN112199776A (zh) * 2020-10-10 2021-01-08 中车大同电力机车有限公司 机车全寿命周期的评估方法、评估系统及优化方法
CN117129332B (zh) * 2023-08-29 2024-05-10 中国航发湖南动力机械研究所 发动机整机腐蚀试验的强度寿命输入数据获取方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103617364A (zh) * 2013-12-04 2014-03-05 南京工业大学 基于小样本的大型回转支承剩余使用寿命预测方法
CN103970965A (zh) * 2014-05-26 2014-08-06 中国航空动力机械研究所 燃气涡轮发动机加速寿命试验试车方法
CN105404756A (zh) * 2015-12-17 2016-03-16 北京航空航天大学 一种涡轮叶盘结构寿命可靠性设计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103617364A (zh) * 2013-12-04 2014-03-05 南京工业大学 基于小样本的大型回转支承剩余使用寿命预测方法
CN103970965A (zh) * 2014-05-26 2014-08-06 中国航空动力机械研究所 燃气涡轮发动机加速寿命试验试车方法
CN105404756A (zh) * 2015-12-17 2016-03-16 北京航空航天大学 一种涡轮叶盘结构寿命可靠性设计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于威布尔分布的疲劳剩余寿命可靠性预测方法;凌丹,何俐萍等;《机械设计》;20110731;第28卷(第7期);全文 *

Also Published As

Publication number Publication date
CN107122531A (zh) 2017-09-01

Similar Documents

Publication Publication Date Title
CN107122531B (zh) 一种基于加速寿命试验的高速叶轮寿命评估方法
US8332164B2 (en) Method for determining fatigue damage in a power train of a wind turbine
CN111597673B (zh) 一种随机振动疲劳的加速试验方法及系统
CN112685939A (zh) 一种基于实测的海上风力机基础疲劳损伤分析方法
Shukla et al. An experimental and FEM modal analysis of cracked and normal steam turbine blade
CN115640666B (zh) 一种基于损伤等效的航空发动机加速任务试车谱编制方法
JP2017096922A (ja) 回転機械の静止部品の欠陥を検出するシステムおよび方法
CN114239175B (zh) 基于非线性蠕变-疲劳交互作用的高温结构件寿命预测方法
CN110595709A (zh) 一种透平发动机叶片许用振幅确定方法
CN110332080B (zh) 一种基于共振响应的风机叶片健康实时监测方法
CN108204342A (zh) 风力发电机的叶片结冰识别方法和装置
CN116542095A (zh) 涡轮叶片寿命评估模型的获取方法及寿命评估方法
KR102253855B1 (ko) 모터 손실 진단방법 및 그 방법을 이용한 모터 손실 진단시스템
CN110646205A (zh) 一种长寿命航空发动机的加速等效试验方法
Ahmad et al. Structural integrity analysis and life estimation of a gas turbine bladed-disc
JP2012167550A (ja) 異常検出装置
CN110516270B (zh) 风力发电机中叶片连接螺栓的疲劳强度评估方法及系统
KR101651414B1 (ko) 풍력터빈의 성능 향상 시스템 및 그 방법
Kurstak et al. Mistuning identification for rotating bladed disks using stationary measurements and reduced order models
CN114542402A (zh) 一种基于多参量分析的风电叶片故障类型在线诊断方法及系统
Basirico et al. Testing of full speed no load operating conditions in a subscale steam turbine test vehicle
CN111188742A (zh) 基于光纤加速度传感器的风力发电机组叶片结冰检测方法
Peng et al. Non-linear vibration behaviour of an axial compressor rotor blade
CN113449431B (zh) 一种基于卸载能的低周疲劳寿命预测方法
Zhang et al. Aerodynamic performance and fatigue damage assessment of wind turbine composite blades using corrected BEM method

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200609

CF01 Termination of patent right due to non-payment of annual fee