一种确定可靠性增长试验计划增长曲线的方法与流程

本发明属于装备可靠性工程技术领域，尤其涉及一种确定可靠性增长试验计划增长曲线的方法。

背景技术：

针对装备可靠性增长试验需求，现有方法以duane模型为基础制定可靠性增长试验计划增长曲线，需要确定的参数为平均故障间隔时间(mtbf)的起始水平mi、平均故障间隔时间(mtbf)增长目标mo、可靠性增长起始试验时间ti、可靠性增长总试验时间to以及可靠性增长率m。

在工程实际中，根据产品可靠性指标、研制进度、研制经费等约束，可以较容易地确定增长计划中的增长目标mo和总试验时间to。但可靠性计划增长曲线的起始水平mi往往较难确定。目前常用的方法为在可靠性增长试验起始阶段开展短时间摸底试验以确定被试产品的可靠性初始水平mi。

根据现有文献、标准规范，一般将可靠性摸底试验的时间纳入可靠性增长时间统计，称为第一阶段试验；后续试验称为第二阶段试验。

第一阶段试验中可能出现的情况包括：

情况1：仅发生了a类故障(定义引自gjb/z77-95《可靠性增长管理手册》，工程实践中，通常是指那些由于经费、时间、技术条件限制或其它原因，确定不进行纠正的系统性故障及所有的残余性故障)；

情况2：同时发生了a类故障和b类故障(定义引自gjb/z77-95《可靠性增长管理手册》，工程实践中，通常是指被确定为需要进行纠正的那些系统性故障)；

情况3：未发生故障；

情况4：仅发生了b类故障。

现有方法仅考虑了上述四种情况的“情况1”，即假设在第一阶段试验期间产品发生的全部为b类故障，产品的可靠性得到增长。在这种假设成立的前提下，按现有方法确定mtbf起始水平mi、mtbf增长目标mo、起始试验时间ti、总试验时间to，进而得到增长率m，进而得到计划增长曲线的函数表达式和计划增长曲线。参见图1，该图为现有技术获取的可靠性增长试验计划增长曲线示意图。

因此，对于上述四种情况中“情况2”、“情况3”、“情况4”，现有方法无法根据第一阶段试验结果直接确定mtbf起始水平mi，进而无法得到计划增长曲线的函数表达式和计划增长曲线。

技术实现要素：

为解决上述问题，本发明提供一种确定可靠性增长试验计划增长曲线的方法，能够获得更精确的计划增长曲线。

一种确定可靠性增长试验计划增长曲线的方法，包括以下步骤：

s1：获取计划增长曲线的名义增长起始点(ti,mi,名义)与实际增长起始点(ti,mi,实际)的函数关系γ：

其中，mi,实际为计划增长曲线的实际增长起始点纵坐标，mi,名义为计划增长曲线的名义增长起始点纵坐标，m为计划增长曲线的增长率，lambertw(·)为朗泊w函数，to为预设的可靠性增长试验的总试验时间，ti为可靠性增长试验中第一阶段试验的累计试验时间，mo为预设的平均故障间隔时间的增长目标；

s2：根据可靠性增长试验的第一阶段试验中发生的责任故障次数r、r次责任故障中a类故障次数ra以及b类故障次数rb，得到实际增长起始点纵坐标mi,实际或名义增长起始点纵坐标mi,名义；

s3：若步骤s2中得到的是实际增长起始点纵坐标mi,实际，则根据函数关系γ确定名义增长起始点纵坐标mi,名义；若步骤s2中得到的是名义增长起始点纵坐标mi,名义，根据duane模型函数关系确定实际增长起始点纵坐标mi,实际；

s4：若步骤s2中得到的是实际增长起始点纵坐标mi,实际，则将步骤s3得到的名义增长起始点纵坐标mi,名义代入函数关系γ，得到计划增长曲线的增长率m；若步骤s2中得到的是名义增长起始点纵坐标mi,名义，根据duane模型函数关系确定计划增长曲线的增长率m；

s5：根据名义增长起始点(ti,mi,名义)、实际增长起始点(ti,mi,实际)以及增长率m，绘制可靠性增长试验计划增长曲线。

进一步地，步骤s1中所述获取计划增长曲线的名义增长起始点(ti,mi,名义)与实际增长起始点(ti,mi,实际)的函数关系γ，具体为：

s101：获取计划增长曲线的增长率m的表达式：

s102：根据duane模型函数关系，获取计划增长曲线的名义增长起始点(ti,mi,名义)与实际增长起始点(ti,mi,实际)的函数关系γ*：

s103：将步骤s101获取的增长率m代入步骤s102中的函数关系γ*，得到计划增长曲线的名义增长起始点(ti,mi,名义)与实际增长起始点(ti,mi,实际)的函数关系γ。

进一步地，步骤s101所述的获取计划增长曲线的增长率m的表达式，具体为：

s101a：根据duane模型函数关系，获取如下关系式：

s102b：将上述关系式中m作为未知数求解，得到计划增长曲线的增长率m的表达式。

进一步地，步骤s2中所述得到实际增长起始点纵坐标mi,实际或名义增长起始点纵坐标mi,名义，具体为：

若第一阶段试验中发生责任故障数r>0，且r次责任故障全部为a类故障，则步骤s2得到实际增长起始点纵坐标mi,实际，具体为：

mi,实际＝ti/r

若第一阶段试验中发生责任故障数r>0，且r次责任故障中同时存在a类故障与b类故障，则步骤s2得到实际增长起始点纵坐标mi,实际，具体为：

mi,实际＝ti/ra

若第一阶段试验中发生责任故障数r＝0，则骤s2得到实际增长起始点纵坐标mi,实际，具体为：

mi,实际＝ti

若第一阶段试验中发生责任故障数r>0，且r次责任故障全部为b类故障，则骤s2得到名义增长起始点纵坐标mi,名义，具体为：

mi,名义＝ti/r。

有益效果：

本发明提供一种确定可靠性增长试验计划增长曲线的方法，能够解决在可靠性增长试验第一阶段期间出现“仅发生a类故障”、“同时发生a类故障和b类故障”、“未发生故障”情况时，无法确定计划增长曲线名义增长起始点的问题，以及可靠性增长试验第二阶段计划增长曲线实际增长起始点的问题，从而能够获得更精确的计划增长曲线。

附图说明

图1为现有技术获取的可靠性增长试验计划增长曲线示意图；

图2为本发明提供的一种确定可靠性增长试验计划增长曲线的方法的流程图；

图3为本发明提供的一种名义增长起始点与实际增长起始点的示意图；

图4为本发明提供的一种可靠性增长试验计划增长曲线示意图；

图5(a)本发明提供的采用本发明方法获取的名义起始点纵坐标和实际起始点纵坐标函数关系示意图；

图5(b)本发明提供的采用传统方法获取的名义起始点纵坐标和实际起始点纵坐标函数关系示意图；

图6为本发明提供的又一种可靠性增长试验计划增长曲线示意图。

具体实施方式

为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述。

实施例一

参见图2，该图为本实施例提供的一种确定可靠性增长试验计划增长曲线的方法的流程图。一种确定可靠性增长试验计划增长曲线的方法，具体包括以下步骤：

s1：获取计划增长曲线的名义增长起始点(ti,mi,名义)与实际增长起始点(ti,mi,实际)的函数关系γ：

其中，mi,实际为计划增长曲线的实际增长起始点纵坐标，mi,名义为计划增长曲线的名义增长起始点纵坐标，m为计划增长曲线的增长率，lambertw(·)为朗泊w函数，to为预设的可靠性增长试验的总试验时间，ti为可靠性增长试验中第一阶段试验的累计试验时间，mo为预设的平均故障间隔时间的增长目标。

需要说明的是，可靠性增长试验中，第一阶段试验计划增长曲线结束点对应的纵坐标定义为名义增长起始点(ti,mi,名义)，第二阶段试验计划增长曲线开始点对应的纵坐标定义为实际增长起始点(ti,mi,实际)，其中，名义增长起始点(ti,mi,名义)与实际增长起始点(ti,mi,实际)的横坐标(时间坐标)相同为第一阶段试验时间ti，而纵坐标不同。参见图3，该图为本实施例提供的名义增长起始点与实际增长起始点的示意图。

需要说明的是，平均故障间隔时间的增长目标mo可以根据根据装备研制要求规定的可靠性指标m规定设定，具体的：

mo＝αm规定(2)

其中，1≤α≤1.5。

需要说明的是，可靠性增长试验的总试验时间to可以根据根据装备研制进度和经费约束设定，可选的，to取值为mo的5～25倍。

下面介绍函数关系γ的一种获取方法，包括以下步骤：

s101：获取计划增长曲线的增长率m的表达式：

s102：根据duane模型函数关系，获取计划增长曲线的名义增长起始点(ti,mi,名义)与实际增长起始点(ti,mi,实际)的函数关系γ*：

s103：将步骤s101获取的增长率m代入步骤s102中的函数关系γ*，得到计划增长曲线的名义增长起始点(ti,mi,名义)与实际增长起始点(ti,mi,实际)的函数关系γ。

需要说明的是，将步骤s101获取的增长率m代入步骤s102中的函数关系γ*后，函数关系γ，即公式(1)中不含变量m，则当增长目标mo、试验时间to、第一阶段试验结束时间ti确定后，公式(1)为实际增长起始点纵坐标mi,实际关于名义增长起始点纵坐标mi,名义的函数。

进一步地，计划增长曲线的增长率m的的表达式获取方法具体为：

s101a：根据duane模型函数关系，获取如下关系式：

s102b：将上述关系式中m作为未知数求解，得到计划增长曲线的增长率m的表达式，即公式(3)。

s2：根据可靠性增长试验的第一阶段试验中发生的责任故障次数r、r次责任故障中a类故障次数ra以及b类故障次数rb，得到实际增长起始点纵坐标mi,实际或名义增长起始点纵坐标mi,名义。

需要说明的是，由于第一阶段试验可能出现的情况包括四种，因此，实际增长起始点纵坐标mi,实际或名义增长起始点纵坐标mi,名义的确定也分为四种情况，具体为：

情况1：若第一阶段试验中发生责任故障数r>0，且r次责任故障全部为a类故障，即ra＝r，则步骤s2得到实际增长起始点纵坐标mi,实际，具体为：

mi,实际＝ti/r(6)

情况2：若第一阶段试验中发生责任故障数r>0，且r次责任故障中同时存在a类故障与b类故障，即ra>0、rb>0，则步骤s2得到实际增长起始点纵坐标mi,实际，具体为：

mi,实际＝ti/ra(7)

情况3：若第一阶段试验中发生责任故障数r＝0，则骤s2得到实际增长起始点纵坐标mi,实际，具体为：

mi,实际＝ti(8)

情况4：若第一阶段试验中发生责任故障数r>0，且r次责任故障全部为b类故障，即rb＝r，则骤s2得到名义增长起始点纵坐标mi,名义，具体为：

mi,名义＝ti/r(9)

s3：若步骤s2中确定的是实际增长起始点纵坐标mi,实际，则根据函数关系γ确定名义增长起始点纵坐标mi,名义；若步骤s2中得到的是名义增长起始点纵坐标mi,名义，根据duane模型函数关系确定实际增长起始点纵坐标mi,实际。

需要说明的是，针对步骤s2中的“情况1”、“情况2”、“情况3”三种情况，根据可靠性增长试验第一阶段试验情况，首先确定的是实际增长起始点(ti,mi,实际)。即公式(1)中mi,实际已知。为绘制可靠性计划增长曲线，需要进一步确定名义增长起始点(ti,mi,名义)的纵坐标mi,名义。可以通过数值遍历法求解公式(1)中的mi,名义，进而得到名义增长起始点(ti,mi,名义)的值。

其中，可靠性增长试验第一阶段试验情况包括：

①累计试验时间ti；

②发生责任故障次数r；

③统计r次责任故障中a类故障(条件不允许改进提高)次数为ra，b类故障次数为rb，其中ra+rb＝r。

需要说明的是，针对步骤s2中的“情况4”，名义增长起始点(ti,mi,名义)已知，可以按传统方法，即公式(1)的前半部分，公式(10)来计算实际增长起始点(ti,mi,实际)。

s4：若步骤s2中确定的是实际增长起始点纵坐标mi,实际，则将步骤s3得到名义增长起始点纵坐标mi,名义代入函数关系γ，得到计划增长曲线的增长率m；若步骤s2中得到的是名义增长起始点纵坐标mi,名义，根据duane模型函数关系确定计划增长曲线的增长率m。

需要说明的是，根据函数关系γ可以得到：

则若步骤s2中确定的是实际增长起始点纵坐标mi,实际，将名义增长起始点纵坐标mi,名义代入公式(11)，即可得到增长率m的精确解。

反之，若步骤s2中确定的是名义增长起始点坐标mi,名义，则根据duane模型函数关系确定计划增长曲线的增长率m；也就是说，增长率m可以按传统方法，即公式(12)计算：

s5：根据名义增长起始点(ti,mi,名义)、实际增长起始点(ti,mi,实际)以及增长率m，绘制可靠性增长试验计划增长曲线。

下面介绍可靠性增长试验计划增长曲线的一种绘制方式，具体包括以下步骤：

①准备双对数坐标纸，要求横坐标最大值大于to，纵坐标最大值大于mo。

②将实际增长起始点(ti,mi,实际)、增长目标点(to,mo)、名义增长起始点(ti,mi,名义)依次绘制在双对数坐标纸上，如图4所示。

③用直线连接实际增长起始点(ti,mi,实际)和增长目标点(to,mo)；用虚线连接名义增长起始点(ti,mi,名义)和实际增长起始点(ti,mi,实际)；以名义增长起始点(ti,mi,名义)为终点绘制与横坐标平行的直线，交于纵坐标，可靠性计划增长曲线绘制完毕，如图4所示。

④写出计划增长曲线m(t)的函数表达式：

自此，就完成了/实现了可靠性增长试验计划增长曲线起始点的确定及计划增长曲线的绘制。

实施例二

根据可靠性研制需求，制定某舰船动力系统电动燃油泵可靠性增长试验计划增长曲线，具体实施步骤如下：

步骤一、分析统计可靠性增长第一阶段试验数据：

①电动燃油泵第一阶段可靠性增长试验累计试验时间t＝300h；

②发生责任故障r＝2次；

③其中a类故障(条件不允许改进提高)次数ra＝2，b类故障次数rb＝0。

步骤二、计算可靠性计划增长曲线实际增长起始点(ti,mi,实际)：

①将可靠性计划增长曲线的名义起始点表示为(ti,mi,名义)、实际增长起始点表示为(ti,mi,实际)；

②名义增长起始点(ti,mi,名义)与实际增长起始点(ti,mi,实际)的横坐标(时间坐标)相同，ti＝300h；

③根据可靠性增长第一阶段试验故障情况确定实际增长起始点的纵坐标：故障数r>0，且全部为a类故障，ra＝2，将实际增长起始点纵坐标设置为，

mi,实际＝ti/ra＝150h

步骤三、确定可靠性增长目标mo和总试验时间to：

①根据燃油泵研制合同：要求其mtbf规定值，m规定≥800h，确定mtbf增长目标mo，

mo＝αm规定＝1.25×800h＝1000h

其中α取值为1.25。

②根据装备研制进度和经费约束确定可靠性增长试验时间to＝5000h。

步骤四、计算名义增长起始点(ti,mi,名义)和增长率m：

①ti＝300h、mi,实际＝150h、to＝5000h、mo＝1000h代入下式，

得到mi,实际与mi,名义的函数关系如图5(a)所示，与mi,实际＝150h对应的名义增长起始点纵坐标mi,名义＝48.82h。也就是说，mi,实际关于mi,名义的函数不存在极小值点。

如果采用gjb/z77《可靠性增长管理手册》等文献中的现有方法，将得到mi,实际与mi,名义的函数关系为，

参见图5(b)，该图为本实施例提供的采用传统方法，即公式(15)计算名义起始点纵坐标和实际起始点纵坐标时，名义起始点纵坐标和实际起始点纵坐标函数关系示意图。由图5(b)可知，mi,实际关于mi,名义的函数存在极小值点(mi,名义＝37.5h、mi,实际＝171.3h)，也就是说，实际起始点纵坐标最小也只能取到171.3h，而现实情况是mi,实际＝150h。由此可见，现有方法在该情况下不满足工程实际需求。

将名义增长起始点纵坐标mi,名义＝48.82h代入下式，

得到，m＝0.6744。

步骤五、绘制可靠性计划增长曲线：

①通过步骤一至步骤四，得到电动燃油泵计划增长曲线参数如下：

●实际增长起始点(ti,mi,实际)＝(300h,150h)；

●增长目标点(to,mo)＝(5000h,1000h)；

●名义增长起始点(ti,mi,名义)＝(300h,48.82h)；

●增长率m＝0.6744。

②在双对数坐标纸中绘制计划增长曲线，如图6所示。

②划增长曲线函数m(t)的表达式为，

由此可见，本实施例首先统计可靠性增长试验第一阶段试验数据：含累计试验时间、累计故障次数、故障类型。然后，根据第一阶段试验数据，评估试验对象可靠性增长初始水平；并根据计划增长曲线特点定义实际增长起始点和名义增长起始点。其次，根据产品可靠性指标、研制进度、经费约束等确定可靠性增长目标、累计试验时间。最后，将计划增长曲线实际起始点设置为第一阶段试验后mtbf评估结果，并根据公式解算名义增长起始点和增长率，得到计划增长曲线函数并绘制计划增长曲线。

当然，本发明还可有其他多种实施例，在不背离本发明精神及其实质的情况下，熟悉本领域的技术人员当然可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。