CN113833641A - 机载燃油泵退化试验方案设计与寿命预测方法 - Google Patents
机载燃油泵退化试验方案设计与寿命预测方法 Download PDFInfo
- Publication number
- CN113833641A CN113833641A CN202111065694.3A CN202111065694A CN113833641A CN 113833641 A CN113833641 A CN 113833641A CN 202111065694 A CN202111065694 A CN 202111065694A CN 113833641 A CN113833641 A CN 113833641A
- Authority
- CN
- China
- Prior art keywords
- fuel pump
- wear
- test
- stress
- degradation
- 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
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04B—POSITIVE-DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS
- F04B51/00—Testing machines, pumps, or pumping installations
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
- B64F5/60—Testing or inspecting aircraft components or systems
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种机载燃油泵退化试验方案设计与寿命预测方法,包括以下步骤:通过对机载燃油泵的失效模式对应的失效机理的分析,获得造成燃油泵失效机理的主要应力,并根据相关标准确定应力水平;基于所选择的应力类型及水平,设计并搭建复杂应力条件下机载燃油泵退化试验系统;基于正交试验思想设计机载燃油泵复杂应力退化试验方案,并采集试验过程中燃油泵出口压力数据;失效寿命在线预测,采用失效物理与数据驱动融合的方法建立燃油泵退化模型,并采用无迹卡尔曼滤波器对模型参数进行更新,使用更新后的模型对失效寿命进行预测。本发明基于正交试验思想制定相应的试验方案,基于失效物理与数据驱动融合的方法对失效寿命进行预测。
Description
技术领域
本发明涉及飞机燃油泵健康管理领域,具体涉及一种机载燃油泵复杂应力退化试验方案设计与寿命预测方法。
背景技术
机载燃油泵用于向散热分系统输油和供油箱输油,给发动机和燃油系统提供一定的流量和压力。作为飞机燃油系统的重要组成部件,机载燃油泵的性能对飞机的飞行安全和作战任务的顺利完成起着至关重要的作用。
目前国内外对机载燃油泵研究较少,仅有的一些关于机载燃油泵的试验研究也仅限于电应力这一单一应力水平。然而,机载燃油泵实际工作环境复杂,工作过程中会受到振动、过载、高低温转换、油液侵蚀等环境应力的影响,单一的电应力条件下的退化试验不能正确地反映其实际退化状态,因此有必要对复杂应力条件下机载燃油泵的退化规律进行研究。
发明内容
针对现有技术存在的问题,本发明提出一种机载燃油泵退化试验方案设计与寿命预测方法,具体包括以下步骤:
步骤1:试验应力与水平选取,通过对机载燃油泵的失效模式对应的失效机理进行分析,获得造成燃油泵失效机理的主要应力,并根据相关标准确定应力水平;
磨损失效是燃油泵最主要的失效机理;选取电应力和机械振动作为主要应力,开展电应力和机械振动耦合应力作用下的机载燃油泵性能退化试验;
根据可靠性鉴定与验收试验GJB899A-2009中有关直升机给燃油泵提供的电应力规定,确定标称电压为115V,高电压为122V,低电压为104V,将其作为试验施加的电应力水平;
按照GJB150.16-1986和GJB150.16A-2009中规定进行振动试验设计;将施加的振动载荷按频率分为低10~100Hz、中200~300Hz、高500~600Hz,三个等级,按功率谱密度PSD分为低1g2/Hz、中10g2/Hz、高20g2/Hz,三个等级,将其作为试验施加的振动应力水平;
步骤2:试验平台搭建,基于上述选择的应力类型及水平,设计并搭建复杂应力条件下机载燃油泵退化试验系统;
机载燃油泵退化试验系统包括机载燃油泵退化试验平台和振动试验台两部分;机载燃油泵退化试验平台用于给燃油泵施加电压大小为104V、115V及122V 的电应力,振动试验台用于给燃油泵施加机械振动应力,具体而言包括9个等级的应力,各等级水平如表2所示;
步骤3:试验方案设计,基于正交试验思想设计机载燃油泵复杂应力退化试验方案,并采集试验过程中燃油泵出口压力数据;
机载燃油泵实际工作中会受到电应力和振动的共同作用,为描述这种作用影响,根据步骤1所选应力类型及其水平,将燃油泵的工作环境描述为不同的电应力梯度和不同的振动量级的耦合作用,耦合作用的实施方法就是正交试验设计的方法;选择能够表征电应力载荷与振动载荷等级程度的电压U、功率谱密度PSD与频率F作为因素设计正交试验,各电压与振动量级的因素水平如表3 所示;
表3正交试验因素水平表
表3详细列出了电应力与机械振动耦合作用的应力水平,为三因素三水平,为此选用L9(34)混合水平正交表设计机载燃油泵复杂应力退化试验,如表2 所示,只需进行9组试验;
表2正交试验方案表
根据表2所列的正交试验方案,在搭建的机载燃油泵退化试验系统上开展各应力水平组合条件下的退化试验并采集燃油泵出口压力退化数据;根据所选试验载荷量级确定两种应力试验时长近似为a:b,因此设定每个试验循环时长为 (a+b)或(a+b)的倍数,其中前a或a倍数时长在燃油泵退化实验台上进行电应力试验,后b或b倍数时长在振动台上进行随机振动试验;
步骤4:失效寿命在线预测,采用失效物理与数据驱动融合的方法建立燃油泵退化模型,并采用无迹卡尔曼滤波器对模型参数进行更新,使用更新后的模型对失效寿命进行预测;
基于Archard磨损模型建立燃油泵出口压力退化模型;轴承磨损中存在黏着磨损和磨粒磨损两种磨损方式,基于Archard模型的黏着磨损公式为:
其中,V表示磨损体积,ks表示磨损系数,x表示相对滑动距离,FN表示接触面的法向载荷,δs表示接触对中较软材料的受压屈服极限;
基于Archard模型的磨粒磨损公式为:
其中,H表示两个接触对象中较软材料的硬度;由以上磨损模型可以看出,磨损体积V与接触面的法向载荷FN、相对滑动距离x以及接触面的特性有关;磨损系数ks与接触面的接触状况有关,根据磨损类型将其分为黏着磨损系数和磨粒磨损系数;
由于在轴承磨损过程中两种磨损方式同时存在,而它们的Archard公式又具有相似的表达形式,因此将二者进行融合,使用强度δs作为衡量两个接触对象中较软材料的磨损特性,将轴承磨损公式统一定义如下:
由于轴承位于电机内部,实际工作过程中无法直接采集其接触面的接触压力,因此通过仿真模拟的方法获得其大小,强度δs的值与材料的特性有关;
轴承接触面的接触压力不会维持恒定,而是随磨损量的变化而变化;而在一个微小的相对滑动距离dx内,其接触压力FN可以认为是恒定的,因此轴承磨损公式进一步写为:
其中,dV为相对滑动距离dx内的磨损体积;
在燃油泵工作过程中,其电机内部的石墨轴承高速旋转,假设转速为n,轴承的半径为R,因此,在dt时间内石墨轴承与轴相对滑动距离dx表示为:
dx=2πR·n·dt
ks、FN与磨损体积V有关,即表示为ks(V)、FN(V),n与t有关,即表示为n(t),于是轴承磨损公式进一步表示为:
对上式进一步变形,得到单位时间内的磨损体积变化,即轴承的磨损率ω:
由上式知,磨损率ω随接触状态的变化而变化;其中磨损体积V为时间的函数,即V(t),因此ks(V)、FN(V)也为时间的函数,即能够表示为ks(V(t))、FN(V(t)),则磨损率ω也能够表示为时间的函数ω(t);于是上式能够表示为:
这一部分中,ks、FN、n表示基础公式里的符号,所以没有下标,进一步地,ks、FN与磨损体积V有关,n与t有关,因此后续写成ks(V)、FN(V)、n(t),再进一步地,由于磨损体积V为时间的函数,即V(t),因此ks(V)、FN(V)也为时间的函数,因此写作ks(V(t))、FN(V(t)),它们之间是层层递进关系;
因此,在时间t0至tT内,轴承磨损体积表示为:
因为该型燃油泵的电机为恒速电机,因此其转速一定,也即n(t)与t无关,于是上式能够进一步表示为:
磨损系数ks表示磨擦副的本质特性,它与材料的性质以及磨擦副的接触特性有关,会随着磨损体积V的变化而变化;在同一个磨损阶段内磨损系数ks能够认为保持恒定,因此上式进一步表示为:
由上式可知,只需知道轴承转动过程中各个时刻磨擦副的接触应力,便能够计算出轴承的磨损状态,然后根据燃油泵出口压力失效阈值所对应的轴承磨损状态便能够计算出燃油泵的寿命;
根据石墨轴承应力仿真分析的结果知,接触压力FN与相对滑动距离x之间近似满足双指数模型,而x=n·t,即接触压力FN与时间t也满足双指数模型,假设其关系式为FN(t)=aebt+cedt,则轴承磨损公式能够写成:
其中,a、b、c、d、、m、p1、p2、p3、p4、p5均为待拟合参数;
燃油泵性能退化规律符合轴承磨损退化规律,因此建立出口压力的退化模型为:
其中,q1,q2,q3,q4,q5为待拟合参数,y为燃油泵出口压力;
使用快速退化阶段的燃油泵出口压力退化数据建立退化模型并对燃油泵失效寿命进行预测,具体为:首先选取快速退化阶段的前W次测量数据对建立的出口压力退化模型参数进行初始化,具体而言,就是使用参数拟合方法对进行计算,得到模型参数q1、q2、q3、q4、q5的初始值;然后通过新获取的燃油泵出口压力退化数据,对模型参数q1、q2、q3、q4、q5不断更新,并预测燃油泵失效寿命,具体而言,就是每获取一次燃油泵出口压力退化数据,就采用模型参数更新方法对模型参数q1、q2、q3、q4、q5进行一次更新,并使用模型参数q1、q2、q3、q4、q5更新后的退化模型对燃油泵的失效寿命进行预测,直到新获取的燃油泵出口压力退化数据达到失效阈值M时,停止对模型参数q1、q2、q3、q4、q5的更新,同时停止对燃油泵失效寿命的预测。
在本发明的一个具体实施例中,a=9,b=1。
在本发明的另一个具体实施例中,W=100。
在本发明的又一个具体实施例中,q1=-0.453,q2=7.4e-3,q3=86.471, q4=2.262e-4,q5=-26.867。
在本发明的再一个具体实施例中,数据拟合方法采用最小二乘法。
在本发明的还一个具体实施例中,模型参数更新方法采用无迹卡尔曼滤波器。
在本发明的其他一个具体实施例中,M=62.3KPa。
本发明的方法能够解决以往关于机载燃油泵试验应力因素单一的缺点,并且能够实现燃油泵失效寿命预测的实时更新。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是本发明实施例的实现流程;
图2是燃油泵的详细结构;
图3是燃油泵退化试验平台结构;
图4是振动试验台结构;
图5是试验剖面;
图6是出口压力退化曲线;
图7是接触压力与磨损体积的关系;
图8是寿命预测结果与误差百分比随数据更新数量变化曲线。
具体实施方式
下面详细描述本发明的实施例,实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
参见图1,本发明提供了一种基于正交试验的机载燃油泵复杂应力条件下的退化试验方案设计和寿命预测的方法,包括以下步骤:
步骤1:试验应力与水平选取,通过对机载燃油泵的失效模式对应的失效机理进行分析,获得造成燃油泵失效机理的主要应力,并根据相关标准确定应力水平。
其中,机载燃油泵为机载燃油系统的核心部件之一,用于向散热分系统输油和供油箱输油。整个工作过程,燃油泵需要保持额定的流量和压力要求。燃油泵的详细结构参见图2(该型燃油泵由南京机电液压工程研究中心提供,文献焦晓璇,景博,黄以锋,羌晓清,刘晓东.基于小波包BP_AdaBoost算法的机载燃油泵故障诊断研究[J].仪器仪表学报,2016,37(09):1978-1988.),主要由电机和离心叶轮两部分组成。电机采用鼠笼式三相异步电动机,由它驱动安装在电机轴上的离心叶轮旋转,将油泵入口处的燃油抽吸到离心叶轮进行增压,经过增压的燃油经蜗室汇集后进入燃油出口管,最后输送至输油管。虽然以该型燃油泵作为示例,但本发明广泛适用于多种型号的机载燃油泵。
机载燃油泵的失效机理主要包括疲劳、磨损和裂纹等。表1给出燃油泵的失效机理及其所影响的材料和相关应力,磨损失效是燃油泵最主要的失效机理。通过对故障的燃油泵拆解分析,发现磨损部位一般发生在轴承、壳体组件以及转子组件等部位,其中电动机的轴承磨损最严重,造成燃油泵出口压力降低、输送燃油能力下降。
表1燃油泵失效机理分析
由表1可知,影响轴承磨损的主要应力包括转速、机械振动、油液污染物、电应力和温度等。由于该型燃油泵的电动机为交流恒速电动机,同时燃油系统中只有部分油路循环,燃油清洁度好,以及燃油具有较好的冷却作用,因此可忽略转速、油液污染度以及温度对轴承磨损的影响。在飞机飞行过程中,燃油泵会接收到飞机供电系统输送来的不同幅值的中频电压,同时会受到由于气流扰动、发动机振动、武器发射等产生的振动的影响,因此本发明实例选取电应力和机械振动作为主要应力,开展电应力和机械振动耦合应力作用下的机载燃油泵性能退化试验。
各应力水平的选取需保证机载燃油泵的失效机理具有一致性。电应力水平优先选用燃油泵在实际使用中的实测电应力数据来制定,根据可靠性鉴定与验收试验GJB899A-2009中有关直升机给燃油泵提供的电应力规定确定标称电压为 115V,高电压为122V,低电压为104V,将其作为试验施加的电应力水平。按照 GJB150.16-1986和GJB150.16A-2009中规定进行振动试验设计。振动载荷的主要参数有振动频率和振动量值,对于机身、发动机舱等区域的振动环境均为随机振动,通常用功率谱密度或加速度均方根值来表示随机振动的振动量值。由于机载燃油泵实际工作过程中,特别是起飞、降落或是导弹发射时,处于强振动环境中,因此结合燃油泵所处工作环境的现场实测振动功率谱,并适当提高其载荷量级,将施加的振动载荷按频率分为低(10~100Hz)、中(200~300Hz)、高 (500~600Hz)三个等级,按功率谱密度(Power Spectral Density,PSD)分为低 (1g2/Hz)、中(10g2/Hz)、高(20g2/Hz)三个等级,将其作为试验施加的振动应力水平。
步骤2:试验平台搭建,基于上述选择的应力类型及水平,设计并搭建复杂应力条件下机载燃油泵退化试验系统。
其中,机载燃油泵退化试验系统包括机载燃油泵退化试验平台和振动试验台两部分。机载燃油泵退化试验平台用于给燃油泵施加电压大小为104V、115V 及122V的电应力,振动试验台用于给燃油泵施加机械振动应力,具体而言包括 9个等级的应力,各等级水平如表2所示。
表2正交试验方案表
机载燃油泵退化试验平台由多个功能模块组成,根据燃油泵的外形、安装方式、技术性能等要求设计了专用测试模块,另外还包括油源模块、测试台架和相关功能附件,其结构参见图3。(中国专利申请“一种飞机燃油泵加速退化试验装置”,申请号202010410066.3)
振动试验台主要由振动台台体、振动控制仪、功率放大器以及冷却系统等部分构成,其实物参见图4(中国专利“振动、温湿度耦合环境中微焊点失效试验系统”,专利号ZL201618001651.2)。通过振动控制仪完成试验类型选择,控制参数、目标谱以及进度表的设置,其控制信号通过功率放大器放大后通过励磁效应作用于台体。振动台通过冷却系统进行散热,以保证系统能够长期稳定的工作。
步骤3:试验方案设计,基于正交试验思想设计机载燃油泵复杂应力退化试验方案,并采集试验过程中燃油泵出口压力数据。
机载燃油泵实际工作中会受到电应力和振动的共同作用,为了描述这种作用影响,根据步骤1所选应力类型及其水平,可将燃油泵的工作环境描述为不同的电应力梯度和不同的振动量级的耦合作用,耦合作用的实施方法就是本发明中提到的正交试验设计的方法。例如,高电压与低频、低PSD值耦合的载荷,标称电压与中频、低PSD值耦合的载荷等。选择能够表征电应力载荷与振动载荷等级程度的电压(U)、功率谱密度(PSD)与频率(F)作为因素设计正交试验,各电压与振动量级的因素水平如表3所示。
表3正交试验因素水平表
表3详细列出了电应力与机械振动耦合作用的应力水平,为三因素三水平,为此选用L9(34)混合水平正交表设计机载燃油泵复杂应力退化试验(董如何,肖必华,方永水.正交试验设计的理论分析方法及应用[J].安徽建筑工业学院学报 (自然科学版),2004(06):103-106.),如表2所示,只需进行9组试验。在本发明的一个具体实施例中,考虑到样本个体的差异性,每组载荷条件选择3台燃油泵进行试验,共需要27台燃油泵。
根据表2所列的正交试验方案,在搭建的机载燃油泵退化试验系统上开展各应力水平组合条件下的退化试验并采集燃油泵出口压力退化数据。根据所选试验载荷量级确定两种应力试验时长近似为9/1,因此设定每个试验循环为10 小时,其中前9个小时在燃油泵退化实验台上进行电应力试验,后1个小时在振动台上进行随机振动试验,共进行87个试验循环。由于振动试验时燃油泵无法运转,因此只在进行电应力试验时采集燃油泵出口流量为12000L/h时的出口压力值,设定每小时采集1次,整个试验共采集783次,其试验剖面参见图5。以上所选每个循环试验时长以及试验循环次数是根据本实施例制定的,根据不同试验对象可以制定不同的循环试验时长以及试验循环次数。图6给出了4号试验的燃油泵出口压力退化曲线。
步骤4:失效寿命在线预测,采用失效物理与数据驱动融合的方法建立燃油泵退化模型,并采用无迹卡尔曼滤波器对模型参数进行更新,使用更新后的模型对失效寿命进行预测。
基于Archard磨损模型建立燃油泵出口压力退化模型。轴承磨损中主要存在黏着磨损和磨粒磨损两种磨损方式,基于Archard模型的黏着磨损公式为:
其中,V表示磨损体积,ks表示磨损系数,x表示相对滑动距离,FN表示接触面的法向载荷,δs表示接触对中较软材料的受压屈服极限。
基于Archard模型的磨粒磨损公式为:
其中,H表示两个接触对象中较软材料的硬度。由以上磨损模型可以看出,磨损体积V与接触面的法向载荷FN、相对滑动距离x以及接触面的特性有关。磨损系数ks与接触面的接触状况有关,根据磨损类型可以将其分为黏着磨损系数和磨粒磨损系数。
由于在轴承磨损过程中两种磨损方式同时存在,而它们的Archard公式又具有相似的表达形式,因此本文将二者进行融合,使用强度δs作为衡量两个接触对象中较软材料的磨损特性,将轴承磨损公式统一定义如下:
由于轴承位于电机内部,实际工作过程中无法直接采集其接触面的接触压力,因此可通过有限元仿真模拟的方法获得其大小,强度δs的值与材料的特性有关。
通过有限元仿真模拟的方法可以获得接触压力与磨损体积之间的关系如图 7所示。从图中可以看出,轴承接触面的接触压力不会维持恒定,而是随磨损量的变化而变化。而在一个微小的相对滑动距离dx内,其接触压力FN可以认为是恒定的,因此轴承磨损公式可进一步写为:
其中,dV为相对滑动距离dx内的磨损体积。
在燃油泵工作过程中,其电机内部的石墨轴承高速旋转,假设转速为n,轴承的半径为R,因此,在dt时间内石墨轴承与轴相对滑动距离dx可表示为:
dx=2πR·n·dt
ks、FN与磨损体积V有关,即表示为ks(V)、FN(V),n与t有关,即表示为n(t),于是轴承磨损公式可进一步表示为:
对上式进一步变形,可得单位时间内的磨损体积变化,即轴承的磨损率ω:
由上式可知,磨损率ω会随着接触状态的变化而变化。其中磨损体积V为时间的函数,即V(t),因此ks(V)、FN(V)也为时间的函数,即可表示为ks(V(t))、 FN(V(t)),则磨损率ω也可表示为时间的函数ω(t)。于是上式可表示为:
这一部分中,ks、FN、n表示基础公式里的符号,所以没有下标,进一步的,在本发明中,ks、FN与磨损体积V有关,n与t有关,因此后续写成 ks(V)、FN(V)、n(t),再进一步,由于磨损体积V为时间的函数,即V(t),因此 ks(V)、FN(V)也为时间的函数,因此写作ks(V(t))、FN(V(t)),它们之间是层层递进关系。
因此,在时间t0至tT内,轴承磨损体积可表示为:
因为该型燃油泵的电机为恒速电机,因此其转速一定,也即n(t)与t无关,于是上式可进一步表示为:
磨损系数ks表示磨擦副的本质特性,它与材料的性质以及磨擦副的接触特性有关,会随着磨损体积V的变化而变化。有研究表明,在同一个磨损阶段内磨损系数ks可认为保持恒定,因此上式可进一步表示为:
由上式可知,只需知道轴承转动过程中各个时刻磨擦副的接触应力,便可以计算出轴承的磨损状态,然后根据燃油泵出口压力失效阈值所对应的轴承磨损状态便可以计算出燃油泵的寿命。
根据石墨轴承应力仿真分析的结果可知,接触压力FN与相对滑动距离x之间近似满足双指数模型,而x=n·t,即接触压力FN与时间t也满足双指数模型,假设其关系式为FN(t)=aebt+cedt,则轴承磨损公式可以写成:
其中,a、b、c、d、、m、p1、p2、p3、p4、p5均为待拟合参数。
燃油泵性能退化规律符合轴承磨损退化规律,因此可建立出口压力的退化模型为:
其中,q1,q2,q3,q4,q5为待拟合参数,y为燃油泵出口压力。
使用快速退化阶段的燃油泵出口压力退化数据建立退化模型并对燃油泵失效寿命进行预测,具体为:首先选取快速退化阶段的前100次测量数据对建立的出口压力退化模型参数进行初始化,具体而言,就是使用最小二乘法对进行计算,得到模型参数的初始值为:q1=-0.453,q2=7.4e-3, q3=86.471,q4=2.262e-4,q5=-26.867;然后通过新获取的燃油泵出口压力退化数据,对模型参数q1、q2、q3、q4、q5不断更新,并预测燃油泵失效寿命,具体而言,就是每获取一次燃油泵出口压力退化数据,就采用无迹卡尔曼滤波器(阙子俊,金晓航,孙毅.基于UKF的轴承剩余寿命预测方法研究[J].仪器仪表学报,2016,37(9):2036-2043.)对模型参数q1、q2、q3、q4、q5进行一次更新,并使用模型参数q1、q2、q3、q4、q5更新后的退化模型对燃油泵的失效寿命进行预测,直到新获取的燃油泵出口压力退化数据达到失效阈值 62.3KPa时,停止对模型参数q1、q2、q3、q4、q5的更新,同时停止对燃油泵失效寿命的预测。图8给出了寿命预测结果及百分比误差随数据更新数量变化的曲线。可以看出,本发明所提方法随着新的退化数据的不断获取与模型参数的更新,对燃油泵失效寿命的预测也越来越准确。
本发明与现有技术方法相比,解决了以往关于机载燃油泵试验应力因素单一的缺点,通过选取应力类型与水平,搭建复杂应力试验系统,并基于正交试验思想设计试验方案,通过构建失效物理与数据驱动融合的退化模型,并采用无迹卡尔曼滤波器对模型参数不断更新,实时更新对失效寿命的预测。经过对本发明实验验证,实验结果良好。
本发明所提机载燃油泵退化试验方案设计与寿命预测方法不只局限于本实施例,也可用于燃油泵疲劳、裂纹、腐蚀等失效机理的分析,根据具体的应力种类确定应力水平大小,然后确定施加应力的试验系统,并基于正交试验思想制定相应的试验方案,最后基于失效物理与数据驱动融合的方法对失效寿命进行预测。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (7)
1.一种机载燃油泵退化试验方案设计与应力影响分析方法,其特征在于,具体包括以下步骤:
步骤1:试验应力与水平选取,通过对机载燃油泵的失效模式对应的失效机理进行分析,获得造成燃油泵失效机理的主要应力,并根据相关标准确定应力水平;
磨损失效是燃油泵最主要的失效机理;选取电应力和机械振动作为主要应力,开展电应力和机械振动耦合应力作用下的机载燃油泵性能退化试验;
根据可靠性鉴定与验收试验GJB899A-2009中有关直升机给燃油泵提供的电应力规定,确定标称电压为115V,高电压为122V,低电压为104V,将其作为试验施加的电应力水平;
按照GJB150.16-1986和GJB150.16A-2009中规定进行振动试验设计;将施加的振动载荷按频率分为低10~100Hz、中200~300Hz、高500~600Hz,三个等级,按功率谱密度PSD分为低1g2/Hz、中10g2/Hz、高20g2/Hz,三个等级,将其作为试验施加的振动应力水平;
步骤2:试验平台搭建,基于上述选择的应力类型及水平,设计并搭建复杂应力条件下机载燃油泵退化试验系统;
机载燃油泵退化试验系统包括机载燃油泵退化试验平台和振动试验台两部分;机载燃油泵退化试验平台用于给燃油泵施加电压大小为104V、115V及122V的电应力,振动试验台用于给燃油泵施加机械振动应力,具体而言包括9个等级的应力,各等级水平如表2所示;
表2 正交试验方案表
步骤3:试验方案设计,基于正交试验思想设计机载燃油泵复杂应力退化试验方案,并采集试验过程中燃油泵出口压力数据;
机载燃油泵实际工作中会受到电应力和振动的共同作用,为描述这种作用影响,根据步骤1所选应力类型及其水平,将燃油泵的工作环境描述为不同的电应力梯度和不同的振动量级的耦合作用,耦合作用的实施方法就是正交试验设计的方法;选择能够表征电应力载荷与振动载荷等级程度的电压U、功率谱密度PSD与频率F作为因素设计正交试验,各电压与振动量级的因素水平如表3所示;
表3 正交试验因素水平表
表3详细列出了电应力与机械振动耦合作用的应力水平,为三因素三水平,为此选用L9(34)混合水平正交表设计机载燃油泵复杂应力退化试验,如表2所示,只需进行9组试验;
根据表2所列的正交试验方案,在搭建的机载燃油泵退化试验系统上开展各应力水平组合条件下的退化试验并采集燃油泵出口压力退化数据;根据所选试验载荷量级确定两种应力试验时长近似为a∶b,因此设定每个试验循环时长为(a+b)或(a+b)的倍数,其中前a或a倍数时长在燃油泵退化实验台上进行电应力试验,后b或b倍数时长在振动台上进行随机振动试验;
步骤4:失效寿命在线预测,采用失效物理与数据驱动融合的方法建立燃油泵退化模型,并采用无迹卡尔曼滤波器对模型参数进行更新,使用更新后的模型对失效寿命进行预测;
基于Archard磨损模型建立燃油泵出口压力退化模型;轴承磨损中存在黏着磨损和磨粒磨损两种磨损方式,基于Archard模型的黏着磨损公式为:
其中,V表示磨损体积,ks表示磨损系数,x表示相对滑动距离,FN表示接触面的法向载荷,δs表示接触对中较软材料的受压屈服极限;
基于Archard模型的磨粒磨损公式为:
其中,H表示两个接触对象中较软材料的硬度;由以上磨损模型可以看出,磨损体积V与接触面的法向载荷FN、相对滑动距离x以及接触面的特性有关;磨损系数ks与接触面的接触状况有关,根据磨损类型将其分为黏着磨损系数和磨粒磨损系数;
由于在轴承磨损过程中两种磨损方式同时存在,而它们的Archard公式又具有相似的表达形式,因此将二者进行融合,使用强度δs作为衡量两个接触对象中较软材料的磨损特性,将轴承磨损公式统一定义如下:
由于轴承位于电机内部,实际工作过程中无法直接采集其接触面的接触压力,因此通过仿真模拟的方法获得其大小,强度δs的值与材料的特性有关;
轴承接触面的接触压力不会维持恒定,而是随磨损量的变化而变化;而在一个微小的相对滑动距离dx内,其接触压力FN可以认为是恒定的,因此轴承磨损公式进一步写为:
其中,dV为相对滑动距离dx内的磨损体积;
在燃油泵工作过程中,其电机内部的石墨轴承高速旋转,假设转速为n,轴承的半径为R,因此,在dt时间内石墨轴承与轴相对滑动距离dx表示为:
dx=2πR·n·dt
ks、FN与磨损体积V有关,即表示为ks(V)、FN(V),n与t有关,即表示为n(t),于是轴承磨损公式进一步表示为:
对上式进一步变形,得到单位时间内的磨损体积变化,即轴承的磨损率ω:
由上式知,磨损率ω随接触状态的变化而变化;其中磨损体积V为时间的函数,即V(t),因此ks(V)、FN(V)也为时间的函数,即能够表示为ks(V(t))、FN(V(t)),则磨损率ω也能够表示为时间的函数ω(t);于是上式能够表示为:
这一部分中,ks、FN、n表示基础公式里的符号,所以没有下标,进一步地,ks、FN与磨损体积V有关,n与t有关,因此后续写成ks(V)、FN(V)、n(t),再进一步地,由于磨损体积V为时间的函数,即V(t),因此ks(V)、FN(V)也为时间的函数,因此写作ks(V(t))、FN(V(t)),它们之间是层层递进关系;
因此,在时间t0至tT内,轴承磨损体积表示为:
因为该型燃油泵的电机为恒速电机,因此其转速一定,也即n(t)与t无关,于是上式能够进一步表示为:
磨损系数ks表示磨擦副的本质特性,它与材料的性质以及磨擦副的接触特性有关,会随着磨损体积V的变化而变化;在同一个磨损阶段内磨损系数ks能够认为保持恒定,因此上式进一步表示为:
由上式可知,只需知道轴承转动过程中各个时刻磨擦副的接触应力,便能够计算出轴承的磨损状态,然后根据燃油泵出口压力失效阈值所对应的轴承磨损状态便能够计算出燃油泵的寿命;
根据石墨轴承应力仿真分析的结果知,接触压力FN与相对滑动距离x之间近似满足双指数模型,而x=n·t,即接触压力FN与时间t也满足双指数模型,假设其关系式为FN(t)=aebt+cedt,则轴承磨损公式能够写成:
其中,a、b、c、d、、m、p1、p2、p3、p4、p5均为待拟合参数;
燃油泵性能退化规律符合轴承磨损退化规律,因此建立出口压力的退化模型为:
其中,q1,q2,q3,q4,q5为待拟合参数,y为燃油泵出口压力;
使用快速退化阶段的燃油泵出口压力退化数据建立退化模型并对燃油泵失效寿命进行预测,具体为:首先选取快速退化阶段的前W次测量数据对建立的出口压力退化模型参数进行初始化,具体而言,就是使用参数拟合方法对进行计算,得到模型参数q1、q2、q3、q4、q5的初始值;然后通过新获取的燃油泵出口压力退化数据,对模型参数q1、q2、q3、q4、q5不断更新,并预测燃油泵失效寿命,具体而言,就是每获取一次燃油泵出口压力退化数据,就采用模型参数更新方法对模型参数q1、q2、q3、q4、q5进行一次更新,并使用模型参数q1、q2、q3、q4、q5更新后的退化模型对燃油泵的失效寿命进行预测,直到新获取的燃油泵出口压力退化数据达到失效阈值M时,停止对模型参数q1、q2、q3、q4、q5的更新,同时停止对燃油泵失效寿命的预测。
2.如权利要求1所述的机载燃油泵退化试验方案设计与寿命预测方法,其特征在于,a=9,b=1。
3.如权利要求1所述的机载燃油泵退化试验方案设计与寿命预测方法,其特征在于,W=100。
4.如权利要求1所述的机载燃油泵退化试验方案设计与寿命预测方法,其特征在于,q1=-0.453,q2=7.4e-3,q3=86.471,q4=2.262e-4,q5=-26.867。
5.如权利要求1所述的机载燃油泵退化试验方案设计与寿命预测方法,其特征在于,数据拟合方法采用最小二乘法。
6.如权利要求1所述的机载燃油泵退化试验方案设计与寿命预测方法,其特征在于,模型参数更新方法采用无迹卡尔曼滤波器。
7.如权利要求1所述的机载燃油泵退化试验方案设计与寿命预测方法,其特征在于,M=62.3KPa。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111065694.3A CN113833641B (zh) | 2021-09-10 | 2021-09-10 | 机载燃油泵退化试验方案设计与寿命预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111065694.3A CN113833641B (zh) | 2021-09-10 | 2021-09-10 | 机载燃油泵退化试验方案设计与寿命预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113833641A true CN113833641A (zh) | 2021-12-24 |
CN113833641B CN113833641B (zh) | 2023-06-30 |
Family
ID=78959048
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111065694.3A Active CN113833641B (zh) | 2021-09-10 | 2021-09-10 | 机载燃油泵退化试验方案设计与寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113833641B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114966297A (zh) * | 2022-08-01 | 2022-08-30 | 成都国营锦江机器厂 | 一种电子调节器综合试验台的应用方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013086822A1 (zh) * | 2011-12-13 | 2013-06-20 | 盛瑞传动股份有限公司 | 自动变速箱控制油泵试验台 |
US20130179104A1 (en) * | 2012-01-09 | 2013-07-11 | GM Global Technology Operations LLC | Dc-motor and fuel pump faults and brush-wear prognosis |
CN105205276A (zh) * | 2015-10-09 | 2015-12-30 | 中国人民解放军国防科学技术大学 | 关节轴承磨损失效物理建模与分析方法 |
CN105352588A (zh) * | 2015-09-08 | 2016-02-24 | 北京航空航天大学 | 无刷直流电机振动检测系统设计 |
JP5990729B1 (ja) * | 2015-04-03 | 2016-09-14 | トライボテックス株式会社 | 汎用劣化曲線作成方法及び機械寿命予測方法、並びに汎用劣化曲線作成プログラム及び機械寿命予測プログラム |
CN106021719A (zh) * | 2016-05-19 | 2016-10-12 | 浙江工业大学 | 一种基于无迹卡尔曼滤波算法的轴承剩余寿命预测方法 |
CN108869263A (zh) * | 2018-05-09 | 2018-11-23 | 中国人民解放军空军工程大学 | 飞机燃油泵综合性能退化试验台 |
CN109870305A (zh) * | 2019-02-28 | 2019-06-11 | 昆明理工大学 | 一种伺服电机试验台 |
CN112067293A (zh) * | 2020-09-29 | 2020-12-11 | 中国航空综合技术研究所 | 自润滑关节轴承磨损寿命预测模型修正方法 |
US20210027001A1 (en) * | 2019-07-23 | 2021-01-28 | Beihang University | Method for Predicting Sealing Reliability of Soft Packing Lithium Ion Battery |
CN113125888A (zh) * | 2021-04-20 | 2021-07-16 | 中国航空综合技术研究所 | 基于故障行为的航空机电产品加速寿命试验方法 |
CN113138337A (zh) * | 2021-04-25 | 2021-07-20 | 西安交通大学 | 基于状态划分与频段同步校正的三相电机故障诊断方法 |
-
2021
- 2021-09-10 CN CN202111065694.3A patent/CN113833641B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013086822A1 (zh) * | 2011-12-13 | 2013-06-20 | 盛瑞传动股份有限公司 | 自动变速箱控制油泵试验台 |
US20130179104A1 (en) * | 2012-01-09 | 2013-07-11 | GM Global Technology Operations LLC | Dc-motor and fuel pump faults and brush-wear prognosis |
JP5990729B1 (ja) * | 2015-04-03 | 2016-09-14 | トライボテックス株式会社 | 汎用劣化曲線作成方法及び機械寿命予測方法、並びに汎用劣化曲線作成プログラム及び機械寿命予測プログラム |
CN105352588A (zh) * | 2015-09-08 | 2016-02-24 | 北京航空航天大学 | 无刷直流电机振动检测系统设计 |
CN105205276A (zh) * | 2015-10-09 | 2015-12-30 | 中国人民解放军国防科学技术大学 | 关节轴承磨损失效物理建模与分析方法 |
CN106021719A (zh) * | 2016-05-19 | 2016-10-12 | 浙江工业大学 | 一种基于无迹卡尔曼滤波算法的轴承剩余寿命预测方法 |
CN108869263A (zh) * | 2018-05-09 | 2018-11-23 | 中国人民解放军空军工程大学 | 飞机燃油泵综合性能退化试验台 |
CN109870305A (zh) * | 2019-02-28 | 2019-06-11 | 昆明理工大学 | 一种伺服电机试验台 |
US20210027001A1 (en) * | 2019-07-23 | 2021-01-28 | Beihang University | Method for Predicting Sealing Reliability of Soft Packing Lithium Ion Battery |
CN112067293A (zh) * | 2020-09-29 | 2020-12-11 | 中国航空综合技术研究所 | 自润滑关节轴承磨损寿命预测模型修正方法 |
CN113125888A (zh) * | 2021-04-20 | 2021-07-16 | 中国航空综合技术研究所 | 基于故障行为的航空机电产品加速寿命试验方法 |
CN113138337A (zh) * | 2021-04-25 | 2021-07-20 | 西安交通大学 | 基于状态划分与频段同步校正的三相电机故障诊断方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114966297A (zh) * | 2022-08-01 | 2022-08-30 | 成都国营锦江机器厂 | 一种电子调节器综合试验台的应用方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113833641B (zh) | 2023-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108733891B (zh) | 用于预测使用寿命周期和低周期疲劳消耗的系统和方法 | |
Khan et al. | Toward an accurate physics-based UAV thruster model | |
EP3184780A1 (en) | Electric windmill pump for improved gearbox durability of a geared fan engine | |
RU2389998C1 (ru) | Способ оценки технического состояния авиационного газотурбинного двигателя | |
CN103528776A (zh) | 一种高阶旋转叶片动力学相似测试实验台及测试方法 | |
CN107271194A (zh) | 燃气轮机压气机振动模拟试验台 | |
CN113833641A (zh) | 机载燃油泵退化试验方案设计与寿命预测方法 | |
CN113574259A (zh) | 用于监控燃气涡轮机的减速器的状态的方法和系统 | |
CA2834022C (en) | Method and system for integrating gas turbine trim balancing system into electronic engine controls | |
Scanavino et al. | An experimental analysis on propeller performance in a climate-controlled facility | |
RU2499240C1 (ru) | Способ вибродиагностики газотурбинного двигателя | |
CN110160946A (zh) | 一种深度调峰汽轮机叶片侵蚀试验系统及其试验方法 | |
CN207516036U (zh) | 一种燃气轮机压气机振动模拟试验台 | |
CN209542085U (zh) | 一种涡轮冷却器等效寿命试验设备 | |
KR101757986B1 (ko) | 가스터빈 시험장치 및 이를 이용한 가스터빈 시험방법 | |
CN208109410U (zh) | 动车组齿轮箱组成高低温试验系统 | |
RU2395723C1 (ru) | Способ эксплуатации насосного агрегата в процессе закачки жидкости в пласт | |
CN109987251A (zh) | 基于应力等效方法的针对薄弱环节的等效寿命试验方法及设备 | |
CN112985721A (zh) | 涡轮转子叶片振动特性的检测装置和检测方法 | |
Vodička et al. | Wear behaviour of vanes for a rotary vane expander with various graphite materials under dry sliding conditions | |
Komshin et al. | Modern Diagnostics of Aircraft Gas Turbine Engines | |
CN113702219A (zh) | 一种金属材料在超低温环境下振动疲劳性能的检测装置及方法 | |
Štorch et al. | Experimental setup for measurement of contra-rotating propellers | |
CN115356122B (zh) | 一种压气机的加温加压试验方法 | |
RU2426902C2 (ru) | Способ определения осевой нагрузки, действующей на упорный подшипник ротора авиационного газотурбинного двигателя |
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 |