一种基于Gamma分布的任务可靠性指标MTBCF验证方法与流程

一种基于gamma分布的任务可靠性指标mtbcf验证方法
技术领域
1.本发明属于可靠性数据分析技术领域，尤其涉及一种基于gamma分布的任务可靠性指标mtbcf验证方法。


背景技术：

2.gjb 899a《可靠性鉴定和验收试验》中的可靠性验证方案仅适用于服从指数分布的序贯/定时/全数试验，核心是计算接收概率，提出接收/拒收的检验策略，最后根据生产方和使用方风险确定验证试验策略参数。然而将该方法应用于平均严重故障间隔时间mtbcf（任务可靠性指标）的指标验证并不严谨，因为gjb 899a中的验证方案是基于电子设备失效率服从指数分布这一前提，在一定试验总时间下，比较失效样本个数和接收/拒收阈值大小，最后根据生产方和使用方两类风险确定策略参数。对于全串联可靠性模型的平均故障间隔时间 mtbf（基本可靠性指标），由于全串联可靠性模型的电子设备失效率服从指数分布，因而gjb 899a完全适用于mtbf指标验证。然而，电子设备任务可靠性模型包含了冗余等系统，并非全串联可靠性模型，因此，非串联可靠性模型的电子设备失效率并不服从指数分布，采用gjb 899a进行mtbcf指标验证不满足失效率服从指数分布的前提，gjb 899a并不适用于mtbcf指标验证。


技术实现要素：

3.为了解决上述问题，本发明提出了一种基于gamma分布的任务可靠性指标mtbcf验证方法，包括以下步骤：
4.采用gamma分布计算产品的接收概率，所述产品为电子设备，电子设备任务可靠性的失效率拟合分布为固定形状参数后的gamma分布；
5.确定接收/拒收产品的验证方案；
6.根据风险要求确定验证方案的参数。
7.进一步地，所述采用gamma分布计算产品的接收概率具体包括以下步骤：
8.步骤1、以gamma分布拟合电子设备的失效率分布，该电子设备包含冗余，电子设备的失效率函数具体为：
9.ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（1）
10.其中，为概率密度函数，为电子设备可靠度；gamma分布概率密度函数的具体公式为：
11.ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（2）
12.式中，为形状参数，为尺度参数， 为伽马函数，t表示
时间；
13.步骤2、由于gamma分布形状参数m反映失效机理和尺度参数反映失效时间量级，因此，固定形状参数；
14.步骤3、通过最小化误差平方和确定形状参数和尺度参数的拟合参数，最小化误差平方和,的公式具体为：
15.ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（3）
16.其中，为gamma分布失效率函数，为电子设备失效率，为根据可靠性预计得到的多个时间点失效率预计值；
17.步骤4、计算产品的接收概率，其具体公式为：
18.ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（4）
19.其中，为固定形状参数后的gamma分布，即电子设备任务可靠性的拟合分布，为电子设备在试验期间发生的整机故障总数；
20.电子设备在试验期间发生的整机故障总数的概率分布为：
21.ꢀꢀꢀꢀ
（5）。
[0022][0023]
进一步地，所述验证方案具体为，在试验停止前样机故障次数为，则拒收产品；否则接收产品；记验证方案为，其中，为试验总时间，r为故障次数。
[0024]
进一步地，所述所述风险要求包括生产方风险要求和使用方风险要求，具体为：当mtbcf指标达到目标值时，拒收产品的概率不超过，当低于最低可接受值时，接收产品的概率不超过。
[0025]
进一步地，所述根据风险要求确定验证方案的参数具体包括以下步骤：
[0026]
将mtbcf指标的目标值和最低可接受值分别转换为尺度参数值：
[0027]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（6）
[0028]
根据如下方程确定验证方案参数：
[0029]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（7）
[0030]
由于故障次数阈值是正整数，通过依次选取整数，在不超过生产方和使用方风险的条件下，寻找最小的试验总时间，即转化为以下优化问题：
[0031]
ꢀꢀꢀꢀꢀꢀꢀꢀ
（8）。
[0032][0033]
进一步地，寻找最小的试验总时间的求解方法如下：
[0034]
a)令；
[0035]
b)计算和；
[0036]
如果，则令，返回b)步；否则取，，计算终止，其中，为最小的试验总时间，为最小试验总时间对应的r值。
[0037]
本发明与现有技术相比，具备的优点在于：
[0038]
本发明通过采用gamma分布计算接收概率，提出接收/拒收的验证方案，最后根据生产方和使用方风险确定验证试验策略参数，可解决mtbcf指标验证目前尚缺严谨有效方法的瓶颈，同时本发明做出mtbcf指标验证判断的试验时间较短，可大幅节约试验的人力物力等资源。
附图说明
[0039]
图1 为gamma分布失效率曲线图。
[0040]
图2 为电子设备整机预计失效率曲线图。
[0041]
图3 为多倍部件失效率预计下电子设备整机失效率曲线图。
[0042]
图4 为多倍部件失效率预计下电子设备整机失效率双对数曲线图。
具体实施方式
[0043]
为了便于理解本技术，下面将参照相关附图对本技术进行更全面的描述。附图中给出了本技术的实施例。但是，本技术可以以许多不同的形式来实现，并不限于本文所描述的实施例。相反地，提供这些实施例的目的是使本技术的公开内容更加透彻全面。
[0044]
除非另有定义，本文所使用的所有的技术和科学术语与属于本技术的技术领域的技术人员通常理解的含义相同。本文中在本技术的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在限制本技术。
[0045]
在进行电子设备mtbcf分布拟合时，由于不同随机分布的可靠度曲线形状非常相似，并不敏感，但是它们的失效率曲线通常差异较大，可以反映各自分布的特点。因此，通过选择失效率曲线拟合效果最好的分布作为拟合分布。
[0046]
由图2可见，预计得到的电子设备整机失效率为增函数，且失效率在后期增长逐渐放缓，因此选择gamma分布作为拟合分布较为合适。考虑到预计值与实际值可能存在偏差，为增强分布选择的稳健性，可选取1倍、2倍、4倍、6倍和10倍的部件失效率预计值分别计算电子设备整机失效率，结果如图3所示，多倍部件失效率下的电子设备整机失效率曲线形状和单倍情况（如图2）没有显著差别，均为开始加速上升后增速放缓趋于平稳，常用随机分布中的威布尔分布、gamma分布可能具有这样的特点。考虑到威布尔失效率为幂函数
，为了确认威布尔是否合适，可以观察失效率的双对数曲线图。如果幂函数拟合效果好，则双对数曲线应为线性；否则代表幂函数拟合效果不好。电子设备整机失效率双对数曲线如图4所示，呈现非线性，所以，电子设备整机失效率不适用威布尔分布拟合，因此，选用gamma分布。gamma分布形状参数大于1时，失效率斜率先增后减，最后趋于平稳，这与电子设备整机（含冗余等）失效率趋势相吻合。
[0047]
因此，本实施例提出了一种基于gamma分布的任务可靠性指标mtbcf验证方法，通过采用gamma分布计算产品的接收概率（所述产品为电子设备），提出接收/拒收的验证方案，最后根据生产方和使用方风险确定验证试验策略参数，可解决mtbcf指标验证目前尚缺严谨有效方法的瓶颈，同时本发明做出mtbcf指标验证判断的试验时间较短，可大幅节约试验的人力物力等资源，具体技术方案如下：
[0048]
1）以gamma分布拟合电子设备（含冗余）的失效率分布
[0049]
gamma分布概率密度函数为：
[0050]
ꢀꢀꢀꢀꢀꢀꢀ
（1）
[0051]
式中，为形状参数，为尺度参数， 为伽马函数，t表示时间。
[0052]
gamma分布期望和方差分别为：
[0053]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（2）
[0054]
电子设备的失效率为：
[0055]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（3）
[0056]
为概率密度函数，为电子设备可靠度；
[0057]
gamma失效率有如下性质（如图1所示）：
[0058]
a)；
[0059]
b)当时，为减函数；
[0060]
c)当时，为常数，此时gamma分布退化为指数分布；
[0061]
d)当时，为增函数。
[0062]
gamma分布的形状参数m的大小决定了失效率曲线的增减性，无论取值，失效率均趋于恒定，为尺度参数的倒数，剩余寿命分布趋近于指数分布。
[0063]
2）一般地，gamma分布形状参数m反映失效机理，尺度参数反映失效时间量级，因此，可以固定形状参数。
[0064]
3）设gamma分布失效率函数为，令电子设备失效率为，根据可靠性预计得到多个时间点失效率预计值，通过最小化误差平方和确定形状参
数和尺度参数的拟合参数，其中，最小化误差平方和,的公式为：
[0065]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（4）
[0066]
多倍部件失效率下整机失效率gamma分布参数拟合结果如表1所示。
[0067]
表1
[0068][0069]
4）试验停止前样机故障次数为，则拒收产品；否则接收产品。记验证方案为，其中，为试验总时间，r为故障次数。
[0070]
电子设备任务可靠性的拟合分布为固定形状参数后的gamma分布，该分布函数记为，和分别为形状参数和尺度参数。因此，电子设备在试验期间发生的整机故障总数的概率分布为：
[0071]
ꢀꢀꢀꢀ
（5）
[0072]
当给定时，产品的接收概率为：
[0073]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（6）
[0074]
将gamma分布形状参数m作为固定参数，其由整机失效率预计值拟合结果确定。由表1，本例中取整，gamma分布待定参数只有。
[0075]
5）进一步，将mtbcf指标的目标值和最低可接受值分别转换为尺度参数值：
[0076]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（7）
[0077]
当时，公式（7）具体为。
[0078]
6）为满足生产方和使用方风险要求，当mtbcf指标达到目标值时，拒收产品的概率不超过，当低于最低可接受值时，接收产品的概率不超过。因此，可列如下方程确定验证方案参数：
[0079]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（8）
[0080]
7）由于故障次数阈值是整数，通过依次选取整数，在不超过生产方和使用方风险的条件下，寻找最小的试验总时间，即转化为以下优化问题：
[0081]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（9）
[0082]
8）最优解求解算法如下：
[0083]
a)令；
[0084]
b)计算和；（其中，m=14）
[0085]
如果，则，返回b)步；否则取，，计算终止，其中，为最小的试验总时间，为最小试验总时间对应的r值。
[0086]
本实施例的mtbcf验证方案算例如表2所示。
[0087]
表2
[0088][0089]
最后，基于gamma分布拟合的验证方案与指数分布方案比较如表3所示。
[0090]
表3
[0091][0092]
由表3可见，基于gamma分布拟合的验证方案试验时间小于指数分布方案；基于gamma分布拟合的验证方案拒收阈值也小于指数分布方案，因此，做出mtbcf检验判断的验证试验时间较短，可节约人力物力等试验资源。
[0093]
本发明的前提为：gjb 899a中的验证方案是基于电子设备失效率服从指数分布这一前提，在一定试验总时间下，比较失效样本个数和接收/拒收阈值大小，最后根据生产方和使用方两类风险确定验证方案参数。因此，对于全串联可靠性模型的平均故障间隔时间mtbf（基本可靠性指标），由于全串联可靠性模型的电子设备失效率服从指数分布，因而gjb 899a适用于mtbf指标验证。然而，电子设备任务可靠性模型包含了冗余等系统，因此，非串联可靠性模型的电子设备失效率并不服从指数分布，使用gjb 899a进行mtbcf指标验证不满足失效率服从指数分布的前提。
[0094]
gamma分布（如图1所示，其中从上而下的线条分别表示形状参数m为0.5、0.8、1、3和6时的失效率的线条）无论形状参数m取值，失效率均趋于恒定，剩余寿命分布趋近于指数分布。由此可见，gamma分布与含冗余等系统的电子设备失效规律相符合，开始时由于冗余部分的作用，电子设备失效率较低（m＞1），随着冗余部分的失效，电子设备可靠性模型逐渐向串联模型转化，即失效率分布趋近于指数分布。
[0095]
因此，使用gamma拟合分布计算接收概率，提出接收/拒收的验证方案较为合理，可解决mtbcf指标验证目前尚缺严谨有效方法的瓶颈，同时本发明做出mtbcf指标验证判断的试验时间较短，可大幅节约试验的人力物力等资源。
[0096]
以上所述仅为本发明的较佳实施例、并不用以限制本发明、凡在本发明的精神和原则之内、所作的任何修改、等同替换、改进等、均应包含在本发明的保护范围之内。

技术特征：
1.一种基于gamma分布的任务可靠性指标mtbcf验证方法，其特征在于，包括以下步骤：采用gamma分布计算产品的接收概率，所述产品为电子设备，电子设备任务可靠性的失效率拟合分布为固定形状参数后的gamma分布；确定接收/拒收产品的验证方案；根据风险要求确定验证方案的参数。2.根据权利要求1所述的基于gamma分布的任务可靠性指标mtbcf验证方法，其特征在于，所述采用gamma分布计算产品的接收概率具体包括以下步骤：步骤1、以gamma分布拟合电子设备的失效率分布，该电子设备包含冗余，电子设备的失效率函数具体为：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（1）其中，为概率密度函数，为电子设备可靠度；gamma分布概率密度函数的具体公式为：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（2）式中，为形状参数，为尺度参数， 为伽马函数，t表示时间；步骤2、由于gamma分布形状参数m反映失效机理和尺度参数反映失效时间量级，因此，固定形状参数；步骤3、通过最小化误差平方和确定形状参数和尺度参数的拟合参数，最小化误差平方和,的公式具体为：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（3）其中，为gamma分布失效率函数，为电子设备失效率，为根据可靠性预计得到的多个时间点失效率预计值；步骤4、计算产品的接收概率，其具体公式为：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（4）其中，为固定形状参数后的gamma分布，即电子设备任务可靠性的拟合分布，为电子设备在试验期间发生的整机故障总数；电子设备在试验期间发生的整机故障总数的概率分布为：
ꢀꢀꢀꢀ
（5）。3.根据权利要求2所述的基于gamma分布的任务可靠性指标mtbcf验证方法，其特征在于，所述验证方案具体为，在试验停止前样机故障次数为，则拒收产品；否则接收产品；记验证方案为，其中，为试验总时间，r为故障次数。4.根据权利要求3所述的基于gamma分布的任务可靠性指标mtbcf验证方法，其特征在
于，所述风险要求包括生产方风险要求和使用方风险要求，具体为：当mtbcf指标达到目标值时，拒收产品的概率不超过，当低于最低可接受值时，接收产品的概率不超过。5.根据权利要求4所述的基于gamma分布的任务可靠性指标mtbcf验证方法，其特征在于，所述根据风险要求确定验证方案的参数具体包括以下步骤：将mtbcf指标的目标值和最低可接受值分别转换为尺度参数值：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（6）根据如下方程确定验证方案参数：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（7）由于故障次数阈值是正整数，通过依次选取整数，在不超过生产方和使用方风险的条件下，寻找最小的试验总时间，即转化为以下优化问题：
ꢀꢀꢀꢀꢀꢀꢀꢀ
（8）。6.根据权利要求5所述的基于gamma分布的任务可靠性指标mtbcf验证方法，其特征在于，寻找最小的试验总时间的求解方法如下：a)令；b)计算和；如果，则令，返回b)步；否则取，，计算终止，其中，为最小的试验总时间，为最小试验总时间对应的r值。

技术总结
本发明公开了一种基于Gamma分布的任务可靠性指标MTBCF验证方法，包括：采用Gamma分布计算产品的接收概率，所述产品为电子设备，电子设备任务可靠性的失效率拟合分布为固定形状参数后的Gamma分布；确定接收/拒收产品的验证方案；根据风险要求确定验证方案的参数。本发明做出MTBCF指标验证判断的试验时间较短，可大幅节约试验的人力物力等资源。可大幅节约试验的人力物力等资源。可大幅节约试验的人力物力等资源。


技术研发人员：王宏 王宇歆
受保护的技术使用者：中国电子科技集团公司第十四研究所
技术研发日：2023.08.18
技术公布日：2023/9/16