[发明专利]基于备件保障数据的正态型单元寿命分布参数估计方法

说明书

本发明涉及一种基于备件保障数据的正态型单元寿命分布参数估计方法，该方法包括：首先根据保障任务的相关数据生成n组候选分布参数，然后初始化似然值；根据k次保障任务的相关数据计算保障任务概率，更新似然值；最后，在更新后的似然值L_j(1≤j≤n)中找到最大似然值，则最大似然值对应的均值和根方差即为正态分布均值参数和根方差参数的估计结果。本发明提出的寿命分布参数估计方法，能在因缺少专门的可靠性试验导致缺少寿命值数据时，有效利用实际工作产生的备件保障数据，较为准确地估计寿命分布参数。

技术领域

本发明涉及产品可靠性预测技术领域，具体涉及一种基于备件保障数据的正态型单元寿命分布参数估计方法。

背景技术

产品由各种单元组成。正态分布常用于描述机械件的寿命分布情况，例如汇流环、齿轮箱、减速器等。正态型单元指寿命服从正态分布的单元，寿命X服从正态分布记作X～N(μ,σ²)，其中μ为均值参数，σ为根方差参数，X的密度函数为

现有技术中，通常以产品的寿命分布函数来定量描述产品的可靠性。寿命分布函数由函数类别和分布参数决定，本质上是一种描述统计规律的函数。因此，只有当寿命样本值达到一定数量规模时，该函数才能较为准确的描述产品寿命的分布规律。在实际工作中，工作环境往往对产品可靠性有较大影响，因此需要关注工作环境下产品的寿命分布规律。在实际工作中，产品工况数据的记录质量与专门的可靠性试验相比往往较差，常常由于记录不及时或记录不全，导致故障发生时间、产品累积工作时间等直接反映产品寿命值的相关数据丢失，只留有保障任务开始前的备件数量、任务成功与否、任务计划时间等数据。这使得基于寿命值数据的分布参数常规估计方法难以从这种备件保障数据中分析产品的可靠性。

发明内容

本发明针对现有技术中存在的技术问题，提供一种基于备件保障数据的正态型单元寿命分布参数估计方法,该方法约定：共有k次保障任务的相关数据，第i次保障任务时该单元的计划工作时间记为Tw_i(以下简称为任务时间)，任务开始前的备件数量记为Nbj_i，保障任务结果记为F_i，1≤i≤k。在保障任务期间，以备件更换故障件的方式来排除故障，若发生的故障次数不大于备件数量，则本次保障任务成功(F_i记为1)，否则视为失败(F_i记为0)。即：以形如[Tw_i Nbj_i F_i]的备件保障数据来描述一次保障任务。在积累了k组备件保障数据后，估计寿命分布参数。

本发明解决上述技术问题的技术方案如下：一种基于备件保障数据的正态型单元寿命分布参数估计方法，包括：

步骤1，根据保障任务的相关数据生成n组候选分布参数(μ_j,σ_j)，1≤j≤n，其中，μ_j表示正态分布的均值参数，σ_j表示正态分布的根方差参数，n为正整数；

步骤2，初始化似然值L_j，令L_j＝0，1≤j≤n；

步骤3，根据k次保障任务的相关数据计算保障任务概率P1，更新似然值L_j；

步骤4，在更新后的似然值L_j(1≤j≤n)中找到最大似然值，记为L_M，则似然值L_M对应的μ_M、σ_M分别为正态分布均值参数和根方差参数的估计结果。

进一步的，所述步骤1包括：

步骤1.1，确定正态分布的均值参数μ1_j1＝μ_min+(j1-1)d1,1≤j1≤n1，其中，μ_max表示正态分布的均值参数上限，μ_min表示正态分布的均值参数下限，n1为正整数，且n1≥2；