基于退化数据的动量轮可靠性分析与剩余寿命预测方法技术

本发明专利技术公开了一种基于退化数据的动量轮可靠性分析与剩余寿命预测方法。分析动量轮的性能退化机理，确定关键性能退化参数，并对退化数据进行预处理；分析退化参数的特点，建立合适的性能退化数据模型，利用极大似然法估计出模型中的未知参数，再根据模型的性质推出关于动量轮寿命的概率密度函数、可靠度函数以及剩余寿命概率密度函数；将估计出的参数值代入各个函数表达式中，根据函数分析动量轮各个测量时刻的可靠度以及剩余寿命。本发明专利技术方法可以实现对动量轮的可靠性估计与剩余寿命预测，解决了现有动量轮寿命预测中存在退化数据选取不合理以及没有考虑实际测量误差的问题。

【技术实现步骤摘要】

本专利技术属于预测与健康管理(PrognosticsandHealthManagement，PHM)研究
，具体涉及一种基于退化数据的动量轮可靠性分析与剩余寿命预测方法。
技术介绍
传统的剩余寿命预测方法是以产品的历史寿命数据为基础，通过对历史寿命数据的分析从而确定产品的剩余寿命分布。随着科技的发展,电子工业、航天工业、机械工业等众多领域出现了大量的长寿命、高可靠性产品。对于这类产品，传统的寿命实验以及加速寿命实验都很难在较短的时间内获得产品的失效数据或寿命数据，这就使得传统的剩余寿命预测很难得到比较满意的结果。同时，该类产品具有成本高、数量少、失效模式以及失效机理都非常复杂的特点，对该类产品的剩余寿命预测已成为国内外关注和研究的技术难题。由于动量轮这类航天产品的可靠性极高，寿命很长，失效机理复杂、价格昂贵，难以进行大样本加速寿命试验，因此目前并没有成熟的动量轮失效机理模型。研究表明高可靠长寿命产品的失效往往与其性能退化存在着必然的联系，因此可以通过识别其性能退化过程中的关键特征参数，进而利用监测的性能退化数据进行可靠性建模，且利用设备退化数据建立其性能退化模型是目前被广泛使用的寿命预测方法。鉴于此，采用基于性能退化数据的分析方法，研究动量轮可靠性建模与剩余寿命预测问题具有重要的理论和实际应用价值。动量轮是一种复杂的设备，由驱动电机、轴承组件、轮体组件、壳体组件四部分组成，通过对量轮退化机理分析可知，其失效与轴承温度、润滑剂系统失效、微振动、电流、转速等因素有关。现有关于动量轮的可靠性分析与寿命预测方法中，大多都是以动量轮的轴温或电流为关键性能参数，通过对其分别进行建模分析，估计出动量轮的可靠性与剩余寿命。卫星运转与太阳的距离存在以年为周期的周期性变化。在空间环境中，太阳的辐射强度与距离密切相关，距离的周期性变化导致的太阳辐射强度的改变最终会影响到动量轮各个部位温度的测量数据，这样使得温度测量值会存在较大的波动。由于动量轮电路中的电流较小，当存在较大的测量噪声时，电流的真实值会被测量噪声所淹没难以得出较为准确的电流值。相比较而言，润滑剂的剩余重量受周围环境影响因素小，工程上易于测量，并且润滑剂过量或缺乏都会影响动量轮轴温和电流值。因此，将润滑剂的剩余重量视为关键性能参数较为合理。
技术实现思路
本专利技术的目的是提供一种基于退化数据的动量轮可靠性分析与剩余寿命预测方法，解决了现有动量轮寿命预测中存在退化数据选取不合理以及没有考虑实际测量误差的问题。本专利技术所采用的技术方案是，基于退化数据的动量轮可靠性分析与剩余寿命预测方法，具体按照以下步骤实施：步骤1、分析动量轮的性能退化机理，确定关键性能退化参数，并对退化数据进行预处理；步骤2、分析步骤1中得出的退化参数的特点，建立合适的性能退化数据模型；步骤3、根据步骤2得到的退化模型，利用极大似然法估计出模型中的未知参数；步骤4、分析步骤2退化模型的特点以及性质，再根据退化模型的性质推出关于动量轮寿命的概率密度函数、可靠度函数以及剩余寿命概率密度函数；步骤5、将步骤3估计出的参数值代入步骤4中的各个函数表达式中，得出函数的具体表达式，从而根据函数分析动量轮各个测量时刻的可靠度以及剩余寿命。本专利技术的有益效果是，由于润滑剂剩余重量容易测量，受周围环境影响小，将动量轮轴系中的润滑剂剩余重量作为关键性能退化参数，使得退化数据比以往将动量轮温度以及电流作为关键性能退化参数更容易获取。实际测量过程中误差是无法避免的，将测量误差考虑其中，使得建立的模型以及估计结果会更加接近真实值，更有利于系统的维护。因此，通过本专利技术可以实现对动量轮的可靠性估计与剩余寿命预测，从而避免因寿命的不可预测而带来的巨大利益损失。附图说明图1是本专利技术基于退化数据的动量轮可靠性分析与剩余寿命预测方法的总体流程图；图2是本专利技术基于退化数据的动量轮可靠性分析与剩余寿命预测方法中不考虑测量误差和考虑测量误差两种情形下动量轮寿命的概率密度函数图；图3是本专利技术基于退化数据的动量轮可靠性分析与剩余寿命预测方法中不考虑测量误差和考虑测量误差两种情形下动量轮寿命的可靠度函数图；图4是本专利技术基于退化数据的动量轮可靠性分析与剩余寿命预测方法中不考虑测量误差和考虑测量误差两种情形下动量轮寿命的剩余寿命概率密度函数图。具体实施方式下面结合附图和具体实施方式对本专利技术进行详细说明。本专利技术基于退化数据的动量轮可靠性分析与剩余寿命预测方法，总体流程图如图1所示，具体按照以下步骤实施：步骤1、分析动量轮的性能退化机理，确定关键性能退化参数。通过失效机理分析，本专利技术将动量轮轴系润滑剂的剩余重量视为关键性能退化参数。由于本专利技术所用数据并非遥测数据，因此并未对数据进行预处理。步骤2、通过分析步骤1中得出的退化参数的特点，建立合适的性能退化数据模型。步骤2中对动量轮的退化建模具体步骤如下：2.1)将动量轮中的轴系润滑剂作为关键性能参数，通过观测退化数据发现，润滑剂剩余重量具有减小的趋势，且具有波动。2.2)根据维纳过程的性质可知，维纳过程不具有严格单调性，适合对具有波动单调的性能退化过程进行建模且该过程已被广泛的用于退化建模中，而动量轮退化数据正符合维纳过程的性质。2.3)根据步骤2.2)的分析，建立如下带有测量误差的退化模型：式(1)中：a为漂移参数，满足a～N(μa,δa2)，X(t)为不含误差的退化量，Y(t)为带有误差的退化量，B(t)用于表征退化过程随机性的标准布朗运动过程，δB为扩散参数，用于表示同类产品的共性，ε为随机测量误差，满足ε～N(0,δε2)，进一步假设a,B(t)和ε是相互独立的。步骤3、根据步骤2得到的退化模型，利用极大似然法估计模型中的未知参数。步骤3中对模型参数估计的具体过程如下：假设Yi(tj)表示第i个动量轮在第j时刻测量得到的性能退化量，则Yi(tj)＝aitj+δBB(tj)+δεεij，1＜i＜n,1＜j＜m(2)令t＝(t1,t2...tm)',Yi＝(Yi(t1),Yi(t2)...Yi(tm))',Y＝(Y1',...Yn')'，其中，Yi服从维正态分布，其均值和方差为μ0＝μat，∑＝δa2tt'+Ω(3)其中Ω＝δB2Q+δε2P根据极大似然原理写出对数似然函数为：其中对式(4)关于μa和δa分别求一阶偏导，并分别令其等于零，得：将式(5)和(6)中的代入式(4)中，得到：通过极大化式(7)，得到δB和δε的极大似然估计值，再将δB和δε的估计值代入式(5)和(6)中得到μa和δa的极大似然估计值，最后估计出模型中的所有未知数。步骤4、分析步骤2退化模型的特点以及性质，再根据模型的性质推出关于动量轮寿命的概率密度函数、可靠度函数以及剩余寿命概率密度函数。步骤4中关于动量轮寿命的概率密度函数、可靠度函数以及剩余寿命概率密度函数各自的具体表达式求解过程如下：4.1)根据首达时间(FirstHittingTime,FHT)的概念来定义寿命，即如果退化过程{X(t),t≥0本文档来自技高网...

【技术保护点】
基于退化数据的动量轮可靠性分析与剩余寿命预测方法，具体按照以下步骤实施：步骤1、分析动量轮的性能退化机理，确定关键性能退化参数，并对退化数据进行预处理；步骤2、分析步骤1中得出的退化参数的特点，建立合适的性能退化数据模型；步骤3、根据步骤2得到的退化模型，利用极大似然法估计出模型中的未知参数；步骤4、分析步骤2退化模型的特点以及性质，再根据模型的性质推出关于动量轮寿命的概率密度函数、可靠度函数以及剩余寿命概率密度函数；步骤5、将步骤3估计出的参数值代入步骤4中的各个函数表达式中，得出函数的具体表达式，从而根据函数分析动量轮各个测量时刻的可靠度以及剩余寿命。

【技术特征摘要】
1.基于退化数据的动量轮可靠性分析与剩余寿命预测方法，具体按照以下步骤实施：步骤1、分析动量轮的性能退化机理，确定关键性能退化参数，并对退化数据进行预处理；步骤2、分析步骤1中得出的退化参数的特点，建立合适的性能退化数据模型；步骤3、根据步骤2得到的退化模型，利用极大似然法估计出模型中的未知参数；步骤4、分析步骤2退化模型的特点以及性质，再根据模型的性质推出关于动量轮寿命的概率密度函数、可靠度函数以及剩余寿命概率密度函数；步骤5、将步骤3估计出的参数值代入步骤4中的各个函数表达式中，得出函数的具体表达式，从而根据函数分析动量轮各个测量时刻的可靠度以及剩余寿命。2.根据权利要求1所述的基于退化数据的动量轮可靠性分析与剩余寿命预测方法，其特征在于，步骤2中对动量轮的退化建模具体步骤如下：2.1)将动量轮中的轴系润滑剂作为关键性能参数，通过观测退化数据发现，润滑剂剩余重量具有单调减小的趋势，且具有波动；2.2)根据维纳过程的性质可知，维纳过程不具有严格单调性，适合对具有单调波动的性能退化过程进行建模且该过程已被广泛的用于退化建模中，而动量轮退化数据正符合维纳过程的性质；2.3)根据步骤2.2)的分析，建立如下带有测量误差的退化数据模型：X(t)=at+δBB(t)Y(t)=X(t)+ϵ---(1)]]>式(1)中：a为漂移参数，满足X(t)为不含误差的退化量，Y(t)为带有误差的退化量，B(t)用于表征退化过程随机性的标准布朗运动过程，δB为扩散参数，用于表示同类产品的共性，ε为随机测量误差，满足进一步假设a,B(t)和ε是相互独立的。3.根据权利要求1所述的基于退化数据的动量轮可靠性分析与剩余寿命预...

【专利技术属性】
技术研发人员：钱富才，彭才华，杜许龙，
申请(专利权)人：西安理工大学，
类型：发明
国别省市：陕西;61