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

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

【技术实现步骤摘要】
一种基于Gamma分布的任务可靠性指标MTBCF验证方法


[0001]本专利技术属于可靠性数据分析
，尤其涉及一种基于Gamma分布的任务可靠性指标MTBCF验证方法。

技术介绍

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

技术实现思路

[0003]为了解决上述问题，本专利技术提出了一种基于Gamma分布的任务可靠性指标MTBCF验证方法，包括以下步骤：
[0004]采用Gamma分布计算产品的接收概率，所述产品为电子设备，电子设备任务可靠性的失效率拟合分布为固定形状参数后的Gamma分布；
[0005]确定接收/拒收产品的验证方案；
[0006]根据风险要求确定验证方案的参数。
[0007]进一步地，所述采用Gamma分布计算产品的接收概率具体包括以下步骤：
[0008]步骤1、以Gamma分布拟合电子设备的失效率分布，该电子设备包含冗余，电子设备的失效率函数具体为：
[0009]ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（1）
[0010]其中，为概率密度函数，为电子设备可靠度；Gamma分布概率密度函数的具体公式为：
[0011]ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（2）
[0012]式中，为形状参数，为尺度参数， 为伽马函数，t表示
时间；
[0013]步骤2、由于Gamma分布形状参数m反映失效机理和尺度参数反映失效时间量级，因此，固定形状参数；
[0014]步骤3、通过最小化误差平方和确定形状参数和尺度参数的拟合参数，最小化误差平方和,的公式具体为：
[0015]ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（3）
[0016]其中，为Gamma分布失效率函数，为电子设备失效率，为根据可靠性预计得到的多个时间点失效率预计值；
[0017]步骤4、计算产品的接收概率，其具体公式为：
[0018]ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（4）
[0019]其中，为固定形状参数后的Gamma分布，即电子设备任务可靠性的拟合分布，为电子设备在试验期间发生的整机故障总数；
[0020]电子设备在试验期间发生的整机故障总数的概率分布为：
[0021]ꢀꢀꢀꢀ
（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分布可能具有这样的特点。考虑到威布尔失效率为幂函数
，为了确认威布尔是否合适，可以观察失效率的双对数曲线图。如果幂函数拟合效果好，则双对数曲线应为线性；否则代表本文档来自技高网...

【技术保护点】

【技术特征摘要】
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分布，即电子设备任务可靠性的拟合分布，为电子设备在试验期间发生的整机故障总数；电子设备在试验期间发...

【专利技术属性】
技术研发人员：王宏，王宇歆，
申请(专利权)人：中国电子科技集团公司第十四研究所，
类型：发明
国别省市：