本发明涉及一种电子整机贮存寿命评估与预测方法，属于系统贮存延寿研究技术领域。。
背景技术：
电子类整机产品是航天器的核心组成部分，正确预估其贮存寿命具有重要的理论、实际意义，军事、经济价值。因为贮存时间不足，贮存期评定一般采取加速试验的方法，即通过施加高应力，使产品在既定加速模型下发生失效，之后结合定时或定数截尾等统计分析方法完成贮存期评定。而现役航天器往往具有高可靠特性，贮存加速试验已经难以产生失效。
传统可靠性理论只关注失效时间、失效次数等突发失效特征，未考虑因产品性能退化产生的可用信息。实际上，性能退化是产品的一种自然属性，包含丰富的可靠性信息。近年来，基于敏感参数变化来表征性能退化的研究已成为新的方向，如何利用性能退化数据设计一种整机贮存寿命评估与预测方法是延寿领域亟需解决的问题。
技术实现要素：
本发明的目的在于克服现有技术的不足，提供一种电子整机贮存寿命评估与预测方法，解决电子整机贮存寿命加速因子和贮存时间参数确定不规范的问题，有针对性根据产品的特性开展加速贮存试验。
本发明目的通过如下技术方案予以实现：
提供一种电子整机贮存寿命评估与预测方法，包括如下步骤：
(1)计算电子整机的加速因子aft：
其中ni为第i种元器件的个数；λui为第i种元器件的通用失效率；afi为第i种元器件的加速因子；β为修正因子，初值为1，m为元器件种类数量；
(2)开展温度为te，时长为t/aft的加速试验，并获得每个整机寿命参数x(t)，t为目标贮存期；获取现有自然贮存数据；
(3)建立寿命参数优化灰度模型：
a1为发展系数，a2为影响因子，a3为灰作用量，t为贮存时间；
采用利用步骤(2)中获得的每个整机寿命参数x(t)替换x(0)(t)，获得加速试验年均退化量α2；
利用步骤(2)中获取的现有自然贮存数据中每个整机寿命参数替换x(0)(t)，获得自然贮存年均退化量α1；
(4)计算修正因子β＝α2/α1；
(5)利用步骤(4)计算的修正因子修正加速试验的整机的加速因子aft，得到加速试验对应的修正后的等效自然贮存时间为βt。
优选的，步骤(3)中获得加速试验年均退化量α2的具体方法为：采用利用步骤(2)中获得的每个整机寿命参数x(t)替换x(0)(t)，获取n台整机加速试验条件下的拟合与预测曲线，对n条拟合与预测曲线取均值，获得整机加速试验条件下的均值拟合与预测曲线，根据该曲线获得加速试验年均退化量α2。
优选的，在获得加速试验年均退化量α2的同时获得n条曲线对应的极大值的拟合与预测曲线和极小值的拟合与预测曲线。
优选的，获得自然贮存年均退化量α1的方法为：利用步骤(2)中获取的现有自然贮存数据中每个整机寿命参数替换x(0)(t)，分别获得n台整机自然贮存条件下的拟合与预测曲线，对n条拟合与预测曲线取均值，获得整机自然贮存条件下的均值拟合与预测曲线；根据该曲线获得自然贮存年均退化量α1。
优选的，当加速试验数条件下的均值拟合与预测曲线数据退化为下降趋势时，利用加速试验极小值的拟合与预测曲线进行寿命预测，当加速试验数条件下的均值拟合与预测曲线数据退化为上升趋势时，利用加速试验极大值的拟合与预测曲线进行寿命预测，修正用于寿命预测的极大值的拟合与预测曲线或极小值的拟合与预测曲线，预测剩余寿命或者预测给定贮存期内是否会失效。
优选的，第i种元器件的加速因子afi的计算方法如下：
其中eai为第i种元器件的激活能，tu为自然贮存环境温度，te为加速试验温度，k为速率常数。
优选的，步骤(1)中ni通过查看元器件清单得到；λui和eai通过《电子设备非工作状态可靠性预计手册》查阅。
优选的，步骤(5)中还包括获得自然贮存条件下n条曲线对应的极大值的拟合曲线与预测曲线和极小值的拟合与预测曲线，用于辅助监控自然贮存条件下是否出现退化失效。
优选的，辅助监控自然贮存条件下是否出现退化失效的具体方法为：如果自然贮存条件下n条曲线对应的极大值的拟合曲线与预测曲线和极小值的拟合与预测曲线与加速试验条件下的极大值的拟合与预测曲线和极小值的拟合与预测曲线趋势不同，则重新进行加速试验，并重新选取寿命参数x(t)。
优选的，寿命参数x(t)为电压或电流。
本发明与现有技术相比具有如下优点：
(1)本发明在分析典型电子整机失效机理的基础上，首先提出了基于arrhenius温度加速模型的整机加速因子计算方法，该方法只需获得3个常见参数即可计算整机加速因子，便于工程应用。
(2)本发明根据加速因子计算结果开展加速试验并获得数据，利用修正gm(1,1)灰色模型预测方法，对自然贮存和加速试验数据的退化规律分别开展研究，并对退化趋势的一致性进行分析，该方法计算简单，拟合精度高。
(3)利用一致性分析结果实现加速模型的闭环修正，进而利用修正后的模型实现整机贮存寿命预测。研究结果可知，基于gm(1,1)灰色模型预测导弹整机贮存寿命的方法简单且具有较高的预测精度,适用于工程应用。
附图说明
图1为本发明加速试验及寿命预测方法原理框图；
图2为本发明自然贮存数据的拟合与预测曲线；
图3为本发明加速试验数据的拟合与预测曲线。
具体实施方式
本发明提供一种弹载电子整机贮存寿命评估与预测方法，结合图1，具体包括如下步骤：
(1)基于元器件信息的整机加速因子计算方法，步骤如下：
温度应力加速试验模型中，最常用的是arrhenius模型：
然而该模型只适用于元器件，不能直接用于整机，基于此新提出一种利用元器件信息计算“整机加速因子”的方法：
式中aft为整机加速因子；λat为加速应力下设备总平均失效率；λut为实际使用时设备总平均失效率。
由于元器件失效率之和为“总失效率”:
由式(2)、(3)可得:
式中m为整机中元器件种数；ni为第i种元器件的数量；λui表示实际使用时，第i种元器件的“通用失效率”；λai为加速试验下，第i种元器件“平均失效率”。
由“加速因子“定义可得，每种元器件有:
λai＝λui·afi(5)
式中afi为第i种元器件的加速因子。
由此可得:
β为修正因子，初值设为1，可根据后续试验结果修正。
式(6)中：
式(7)中eai为第i种元器件的激活能。
上述式(6)即为基于arrhenius模型的“整机加速因子”计算方法，该方法只需获得ni,λui,eai3个参数即可。其中ni为某类元器件的个数，可以通过查阅“元器件清单”获得；λui与eai分别是某种元器件的“失效率”和“加速因子”，可通过查阅gjb/z108a-2006《电子设备非工作状态可靠性预计手册》获得，tu为自然贮存温度，te为加速试验温度，k为速率常数。该方法利用元器件的“数量”、“通用失效率”和“激活能”作为加权因子，实现了整机加速因子的计算，为整机级加速试验奠定了基础。
(2)开展加速温度为te，为目标贮存期为t，时长为t/aft的加速试验，并获得每台整机寿命参数x(t)；同时提取已有的自然贮存数据，用于数据分析。
(3)基于优化gm(1,1)灰色模型的数据拟合与预测：
gm(1,1)模型是灰色预测理论能够很好地处理小子样数据的拟合与预测问题，但常会出现精度较差的现象，已有的改进方法可以一定程度上改善拟合精度，但无法根除由拟合方法缺陷所造成的固有误差。
本发明提出了一种优化的gm(1,1)灰色模型，该模型是对服从非齐次指数增长规律的数据建模，克服了原有gm(1,1)模型指数规律的缺点。本发明首先给出gm(1,1)模型的推广微分形式，考虑初值的影响，建立了优化gm(1,1)模型，并通过精确的离散化求取参数，实例显示新模型的拟合与预测精度较原模型高。
gm(1,1)模型的建模基础是假定原始数据序列服从近似指数规律，即x(0)(k)≈be-ak，但实际情况是，许多数列并不服从上述规律，若仍用上述方法建模，会影响预测效果。对服从非齐次指数增长规律的数列，即x(0)(k)≈be-ak+c，给出推广gm(l,1)模型的白化方程可表成如下的初值问题：
其中a1为发展系数，a2为影响因子，a3为灰作用量。
一阶还原计算得到的预测序列为：
经优化gm(l,1)灰色模型的精确离散形式为：
整机数量为n，对于每台整机的使用实际的自然贮存数据作为x(0)(t)带入式10，t为时间；对x(0)(t)进行一阶累加，获得x(1)(t)数据；由最小二乘估计得到式10中a1、a2、a3的估计值，分别获得n台整机自然贮存条件下的拟合与预测曲线，对n条拟合与预测曲线取均值，获得整机自然贮存拟合与预测的均值曲线；同时获得n条曲线对应的极大值、极小值的拟合与预测均值曲线，参见图2，两条极值曲线用于监控是否会出现退化失效。
对每台整机使用步骤(2)获得的寿命参数数据作为x(0)(t)带入式10，t为时间；对x(0)(t)进行一阶累加，获得x(1)(t)数据；由最小二乘估计得到式10中a1、a2、a3的估计值，分别获得加速试验条件下n台整机的拟合与预测曲线，进而获得拟合与预测的均值曲线，根据该曲线可获得加速试验年均退化量α2；同时获得n条曲线对应的极大值、极小值的拟合与预测均值曲线，参见图3。
(4)对比自然贮存和加速试验数据的一致性，实现加速模型的闭环修正，进一步实现贮存寿命预测，具体方法如下：
根据自然贮存拟合与预测的均值曲线获得自然贮存年均退化量α1，根据加速试验拟合与预测的均值曲线获得加速试验数年均退化量α2。
更新修正因子：
β＝α2/α1
(5)修正整机加速试验的加速因子aft，得到修正后的等效自然贮存时间为βt。当数据退化趋势为下降时，可利用极小值的拟合与预测曲线进行寿命预测，当数据退化趋势为上升时，可利用极大值的拟合与预测曲线进行寿命预测，修正寿命预测的极大值或极小值的拟合与预测曲线(横轴坐标t修改为βt)，可预测剩余寿命或给定贮存期内整机是否会退化失效。
通过图2、3可对比知：自然贮存年均退化量0.0167v；恒定应力加速试验设定80℃、7天一循环为一年伪寿命，实测每循环循环平均退化量为0.033v；恒加试验每周期平均退化量是自然贮存的2.0倍。即80℃高温下加速因子为理论值的2.0倍，2.0即为修正因子，可用于修正加速试验结果。由修正后的恒加试验预测曲线可知：在30年贮存期内，该整机不会产生退化失效。
本发明利用优化gm(1,1)预测方法，对加速试验和自然贮存数据分别进行了拟合和预测，通过自然贮存和加速试验数据退化趋势的一致性比对，可实现加速模型的闭环修正，利用修正后的模型可以实现整机贮存寿命的预测。经过“残差检验法”和“后验差检验法”对拟合预测精度进行误差比较可知，优化gm(1,1)模型在5年内预测精度为一级，显著优于传统方法。
本发明适应性较强，可用于为各型武器系统电子整机延寿研究。
以上所述，仅为本发明最佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。