一种正态型单元的寿命分布参数估计方法与流程

本发明涉及产品质量检测
技术领域：
，具体涉及一种正态型单元的寿命分布参数估计方法。
背景技术：
：产品可靠性是一种描述产品质量的核心属性，常用产品寿命的分布类型和参数来表达产品的可靠性。准确获知产品的可靠性，是开展产品的可靠性增长、维修性/保障性设计等工作的前提。在专门的可靠性试验中，一般能实时、在线监测产品的完好状态：一旦产品发生故障能立刻被发现，因此能获得产品寿命x的准确数值。在获得足够数量的寿命数据后即可分析出产品的寿命分布类型及参数。但在工作场景下，并不见得能针对产品配备在线监测设备，因而不能实时监测产品的完好状态。工作场景下更常见的做法是定期或不定期对产品做完好性检查。假定产品投入使用时刻为零时刻，如果在检查时刻tc，产品状态为完好，这意味着该产品的寿命x大于tc；如果在检查时刻tc，产品状态为故障，这意味着该产品的寿命x小于tc。与有寿命准确数值的x相比，[检查时刻tc状态(完好或故障)]是删失了部分信息的寿命数据。目前，在理论上还没有利用这种删失型数据准确估计寿命分布参数的方法。相对于标准的可靠性试验场景，当工作环境、使用方式等发生变化时，产品的实际寿命往往可能随之发生变化，因此，即便掌握了产品在可靠性试验场景下的寿命分布规律，也仍有必要去了解产品在工作场景下的实际寿命分布规律。技术实现要素：本发明针对现有技术中存在的技术问题，提出了一种利用删失型数据估计正态型单元的寿命分布参数的近似方法，其估计精度能满足工程要求。本发明解决上述技术问题的技术方案如下：一种正态型单元的寿命分布参数估计方法。产品由各种单元组成。正态分布常用于描述机械件的寿命分布情况，正态型单元指寿命服从正态分布的单元，寿命x服从正态分布记作x～n(μ,σ2)，其中μ为均值参数，σ为根方差参数，x的密度函数f(x)为：假定：单元投入使用时刻为零时刻，同批次的单元同时投入使用，且各批次单元的工作场景相似。在第i次检查时，检查时刻记为tci，该批次产品中，完好单元的数量记为nri，故障单元的数量记为nfi。共完成了k次检查。本方法包含以下步骤：步骤1，根据正态型单元使用寿命数据生成n组候选分布参数(μ2j,σ2j)，1≤j≤n，其中，μ2j表示正态分布的均值参数，σ2j表示正态分布的根方差参数，n为正整数；步骤2，初始化似然度pj，令步骤3，根据k次检查结果获取的完好单元的数量nri、故障单元的数量nfi以及检查时刻tci，依次更新似然度pj；步骤4，在更新后的似然度pj(1≤j≤n)中找到最大似然度，记为pm，则似然度pm对应的μ2m、σ2m分别为正态分布均值参数和根方差参数的估计结果。进一步，所述步骤1具体包括：步骤1.1，确定正态分布的均值参数μj1＝μmin+(j1-1)d1，1≤j1≤n1，其中，μmax表示正态型单元寿命分布的均值参数上限，μmin表示正态型单元寿命分布的均值参数下限，n1为正整数，且n1≥2；步骤1.2，确定正态分布的根方差参数σj2＝σmin+(j2-1)d2，1≤j2≤n2，其中，σmax表示正态型单元寿命分布的根方差参数上限，σmin表示正态型单元寿命分布的根方差参数下限，n2为正整数，且n2≥2；步骤1.3，取n＝n1×n2，由μj1和σj2进行遍历组合获得n组候选的分布参数(μ2j,σ2j)，1≤j≤n；进一步，步骤1.3中所述的遍历采用以下方式实现：令j＝1；第一层循环中遍历j1＝1:n1，第二层循环中遍历j2＝1:n2，令μ2j＝μj1；σ2j＝σj2；j＝j+1；其中，μmax≥μj1≥μmin，σmax≥σj2≥σmin。进一步，所述步骤3具体包括。步骤3.1，令i＝1，i表示检查次数；步骤3.2，遍历计算wj，1≤j≤n，令其中tci表示第i次检查时的检查时刻，nri为第i次检查时完好单元的数量，nfi为第i次检查时故障单元的数量；步骤3.3，遍历更新似然度pj，令步骤3.4，令i＝i+1，若i≤k则转步骤3.2，否则转步骤4，k为总检查次数。本发明的有益效果是：本方法的参数估计结果在总体上能“跟随”理论成熟方法的参数估计结果，其估计精度能满足工程要求。附图说明图1为本发明方法流程图；图2为采用理论成熟方法和本发明方法的仿真结果示意图。具体实施方式以下结合实施例对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。产品由各种单元组成。正态分布常用于描述机械件的寿命分布情况，例如汇流环、齿轮箱、减速器等。正态型单元指寿命服从正态分布的单元，寿命x服从正态分布记作x～n(μ,σ2)，其中μ为均值参数，σ为根方差参数，x的密度函数为假定：单元投入使用时刻为零时刻，同批次的单元同时投入使用，且各批次单元的工作场景相似。在第i次检查时，检查时刻记为tci，该批次产品中，完好单元的数量记为nri，故障单元的数量记为nfi。共完成了k次检查。实施例1本方法，如图1所示，步骤如下：1、确定候选的寿命分布参数生成n组候选分布参数(μ2j,σ2j),1≤j≤n，其中，μ2j表示正态分布的均值参数，σ2j表示正态分布的根方差参数，n为正整数；生成候选分布参数的具体实现方法为：1)确定正态分布的均值参数μj1＝μmin+(j1-1)d1,1≤j1≤n1，其中，μmax表示正态型单元寿命分布的均值参数上限，μmin表示正态型单元寿命分布的均值参数下限，n1为正整数，且n1≥2；2)确定正态分布的根方差参数σj2＝σmin+(j2-1)d2,1≤j2≤n2，其中，σmax表示正态型单元寿命分布的根方差参数上限，σmin表示正态型单元寿命分布的根方差参数下限，n2为正整数，且n2≥2；3)取n＝n1×n2，由μj1和σj2进行遍历组合获得n组候选的分布参数(μ2j,σ2j),1≤j≤n。其中，遍历的方式可以通过以下方式实现：令j＝1；第一层循环中遍历j1＝1:n1，第二层循环中遍历j2＝1:n2，令μ2j＝μj1；σ2j＝σj2；j＝j+1。其中，μmax≥μj1≥μmin，σmax≥σj2≥σmin。2、初始化似然度初始化似然度pj,1≤j≤n，令3、遍历调整似然度3.1令i＝13.2遍历计算wj,1≤j≤n，令其中3.3遍历更新似然度pj,1≤j≤n，令3.4更新i，令i＝i+1，若i≤k则转步骤3.2，否则转步骤4。4、输出寿命分布参数估计结果在所有似然度pj(1≤j≤n)中找到最大似然度，记为pm，则μ2m、σ2m分别为正态分布均值参数和根方差参数的估计结果。实施例2某正态型单元的10次状态检查结果如下表，试估计其寿命分布的均值参数和根方差参数。检查序号检查时刻h故障单元的数量完好单元的数量1970232550053140054830055280056420147111041812504191380501069014计算过程如下：1、确定候选的寿命分布参数按照以往经验，估计该单元的均值参数在500h～3000h范围内，以500为步长；估计该单元的根方差参数在100～500范围内，以100为步长；共生成30个候选的分布参数(μ2j,σ2j),1≤j≤30。2、初始化似然度初始化似然度pj,1≤j≤30，令3、遍历调整似然度3.1令i＝13.2遍历计算wj,1≤j≤30，令其中3.3遍历更新似然度pj,1≤j≤30，令3.4更新i，令i＝i+1，若i≤10则转3.2，否则转4。表1列出了第i次更新后的似然度。4、输出寿命分布参数估计结果在所有似然度pj(1≤j≤30)中最大似然度为p8，则μ28＝1000,σ28＝300分别为寿命分布均值参数和根方差参数的估计结果。表1第i次更新后的似然度实施例3可建立以下仿真模型来模拟对单元的检查过程。假定单元的实际寿命服从正态分布n(μ,σ2)，共进行k次检查，记第i次检查时刻为tci，同批次的单元同时投入使用，第i批次的单元数量为ni。1)令i＝12)随机产生ni个随机数simtij,1≤j≤ni，这些随机数服从正态分布n(μ,σ2)。3)在simtij(1≤j≤ni)中，找到大于tci的随机数，其数量为完好单元的数量记为nri，故障单元数量nfi为ni-nri。4)更新i，令i＝i+1。若i≤k则转2)，否则本次模拟k次检查结束。对于以上仿真模型得到的tci、nri、nfi，可采用本发明方法用于估计单元寿命的分布参数。对于以上仿真模型得到的simtij，可采用理论上成熟的方法用于估计单元寿命的分布参数。以单元的实际寿命服从正态分布n(1000,3002)为例，共进行10次检查，同批次的单元同时投入使用，第i批次的单元数量为5，利用上述仿真模型开展多次仿真，得到大量模拟检查结果并进行分布参数估计，并对多次估计结果进行统计分析。针对模拟的删失型检查结果数据tci、nri、nfi，采用本文方法估计的正态分布均值参数μ的均值为1005.0、根方差为50.0，根方差参数σ的均值为296.0、根方差为74.4。针对模拟的寿命数据simtij，采用理论方法估计的正态分布均值参数μ的均值为1003.6、根方差为42.2，根方差参数σ的均值为296.3、根方差为24.1。图2显示了利用10次仿真结果，分别针对寿命型数据和对应的删失型数据，各自采用理论成熟方法和本文方法得到的寿命分布参数估计结果。从图2来看，本文方法的参数估计结果在总体上能“跟随”理论成熟方法的参数估计结果。大量仿真验证结果表明，本文方法有较好的估计精度，满足工程应用要求。以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。当前第1页12