一种按次使用武器装备的可靠度评估方法与流程

本发明属于武器装备可靠性工程
技术领域：
，尤其涉及一种按次使用武器装备的可靠度评估方法。
背景技术：
：针对按次使用的武器装备，现有的解决方案是忽略因装备使用过程的离散性导致其故障概率分布离散化这一现象，依然按照gjb899a-2009《可靠性鉴定和验收试验》等标准、文献中给出的基于指数分布假设的方法计算按次使用武器装备的可靠度。指数分布是一种连续概率分布。通常连续时间工作的电子产品(如舰用计算机、舰用通信设备等)及大型复杂可修装备(如舰船电力系统、动力系统等)的故障时间服从指数分布，其故障时间可以是装备累计工作时间轴上的任意时刻(如101.032小时、1000.43小时等)。但对于诸如舰载机弹射器、舰炮等按次使用的武器装备，其故障发生的时间刻度以次数计，为正整数(如100次、400次、10000次等)。因此，目前工程中利用基于指数分布的方法评估按次使用武器装备的可靠性存在以下三个方面的不足：一是理论上不完备，指数分布为连续变量分布函数不能准确描述按次使用武器装备非连续的故障发生次数变量；二是物理意义不匹配，指数分布所蕴含的概率意义与按次使用武器装备故障发生的概率事件不匹配；三是计算结果不准确，由于存在理论上的不完备和物理意义的不匹配，基于指数分布假设的方法给出的评估结果为近似结果。技术实现要素：为解决上述问题，本发明提供一种按次使用武器装备的可靠度评估方法，能够建立按次使用武器装备完备的可靠度评估模型，提高评估结果的准确度。一种按次使用武器装备的可靠度评估方法，包括以下步骤：获取可靠度的点估计值其中，为所述武器装备使用到n次时发生故障的概率，为所述武器装备单次使用不发生故障的概率v的极大似然估计值；获取可靠度的单侧置信下限值r单,l(n)：其中，v单,l为所述武器装备单次使用不发生故障的概率v的单侧置信下限；获取可靠度的双侧置信区间的下限值r双,l(n)与上限值r双,u(n)：其中，v双,l为所述武器装备单次使用不发生故障的概率v的双侧置信区间的下限，v双,u为所述武器装备单次使用不发生故障的概率v的双侧置信区间的上限；根据可靠度的点估计值可靠度的单侧置信下限值r单,l(n)以及可靠度的双侧置信区间的下限值r双,l(n)与上限值r双,u(n)评估可靠度是否满足产品可靠度研制要求。进一步地，所述武器装备单次使用不发生故障的概率v的极大似然估计值的获取方式为：构造极大似然函数l(v)：l(v)＝v(n-z)wz其中，n为所述武器装备进行的有效可靠性试验的总次数，z为所述武器装备在n次可靠性试验中累计发生的故障次数，v为所述武器装备单次使用不发生故障的概率，w为所述武器装备单次使用发生故障的概率；获取极大似然函数l(v)的导数l′(v)：l′(v)＝v(n-z-1)(1-v)z-1[(n-z)(1-v)-zv]令l′(v)＝0，则有：(n-z)(1-v)-z×v＝0求解上式，得到所述武器装备单次使用不发生故障的概率v的极大似然估计值为：进一步地，所述武器装备单次使用不发生故障的概率v的单侧置信下限v单,l的获取方式为：假设置信度为γ，基于二项分布置信下限计算方法得到单侧置信下限v单,l满足如下关系式：其中，表示在n次有效可靠性试验中随机出现k次故障的组合数；采用数值遍历法求解上式，得到置信度为γ时，所述武器装备单次使用不发生故障的概率v的单侧置信下限v单,l。进一步地，所述武器装备单次使用不发生故障的概率v的双侧置信区间的下限v双,l与上限v双,u的获取方式为：假设置信度为γ，基于二项分布置信区间计算方法得到双侧置信区间的下限v双,l与上限v双,u分别满足如下关系式：其中，表示在n次有效可靠性试验中随机出现k次故障的组合数；采用数值遍历法求解上式，得到置信度为γ时，所述武器装备单次使用不发生故障的概率v的双侧置信区间的下限v双,l与上限v双,u。有益效果：本发明提供一种按次使用武器装备的可靠度评估方法，给出可靠度的点估计值、可靠度的单侧置信下限值、可靠度的双侧置信区间的下限值以及上限值的计算方法，得到的点估计值、单侧置信下限值、双侧置信区间的下限值以及上限值，能够全面评估按次使用武器装备的可靠度，从而解决诸如舰载机弹射器等按次使用武器装备的可靠度无法准确建模、物理意义不匹配等问题，能够建立按次使用武器装备完备的可靠度评估模型，提高评估结果的准确度。附图说明图1为本发明提供的一种按次使用武器装备的可靠度评估方法的流程图；图2为本发明提供的武器装备使用次数与可靠度的函数关系示意图。具体实施方式为了使本
技术领域：
的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述。实施例一参见图1，该图为本实施例提供的一种按次使用武器装备的可靠度评估方法的流程图。一种按次使用武器装备的可靠度评估方法，包括以下步骤：s1：获取可靠度的点估计值其中，为所述武器装备使用到n次时发生故障的概率，为所述武器装备单次使用不发生故障的概率v的极大似然估计值。进一步地，所述武器装备单次使用不发生故障的概率v的极大似然估计值的获取方式为：s101：构造极大似然函数l(v)：l(v)＝v(n-z)wz(2)其中，n为所述武器装备进行的有效可靠性试验的总次数，z为所述武器装备在n次可靠性试验中累计发生的故障次数，v为所述武器装备单次使用不发生故障的概率，w为所述武器装备单次使用发生故障的概率，其中w＝1-v。下面对第z次故障发生所述武器装备已经进行的有效可靠性试验次数fz进行详细说明。假设某按次使用的武器装备共开展n次有效可靠性试验，期间发生z(z≤n)次故障。参见表1，该表为每次故障发生时，武器装备已经进行的有效可靠性试验次数。表1故障第1次故障第2次故障…第z次故障试验次数f1f2…fz其中，{f1,f2,…,fz}小于等于n。根据表1，将第1至第z次故障发生事件依次计为概率事件a1,a2,…,az，则能够得到每次故障发生前无故障连续试验次数，如表2所示。表2基于几何分布，建立各概率事件a1,a2,…,az的概率模型。第1次故障发生对应概率事件a1，发生概率为：其中，其中p(·)代表求概率事件发生概率的函数，v为所述武器装备单次使用不发生故障的概率，w为所述武器装备单次使用发生故障的概率，显然w＝1-v。同理，对于第i(1≤i≤z)次故障事件ai，其发生概率为：如果有效可靠性试验在发生第z次故障后截止，则fz＝n。如果可靠性试验在发生第z次故障后继续试验，即fz<n，则装备在发生第z次故障后又连续无故障试验n-fz次，将该事件表示为ao，其发生概率为：至此，得到各概率事件a1,a2,…,az的概率模型。需要说明的是，对于有效可靠性试验在发生第z次故障后全部结束的情况(fz＝n)，极大似然函数为：其中，p(ai)为所述武器装备发生第i次故障事件对应的概率，si为第i次故障发生前无故障连续试验次数。对于有效可靠性试验在发生第z次故障后继续试验的情况(fz<n)，极大似然函数为：其中，p(ao)为所述武器装备在发生第z次故障后又连续无故障试验n-fz次的概率，fz为第z次故障发生所述武器装备已经进行的有效可靠性试验次数，si为第i次故障发生前无故障连续试验次数。由此可见，无论可靠性试验在发生第z次故障后截止，还是继续试验，极大似然函数的形式均为：l(v)＝v(n-z)wz(8)s102：获取极大似然函数l(v)的导数l′(v)：l′(v)＝v(n-z-1)(1-v)z-1[(n-z)(1-v)-zv](9)令l′(v)＝0，则有：(n-z)(1-v)-z×v＝0(10)需要说明的是，对极大似然函数l(v)求导的具体过程为：上式中v(n-z-1)和(1-v)z-1显然不为0，因此l′(v)＝0等价于：(n-z)(1-v)-z×v＝0(12)s103：求解上式，得到所述武器装备单次使用不发生故障的概率v的极大似然估计值为：s2：获取可靠度的单侧置信下限值r单,l(n)：其中，v单,l为所述武器装备单次使用不发生故障的概率v的单侧置信下限。进一步地，所述武器装备单次使用不发生故障的概率v的单侧置信下限v单,l的获取方式为：s201：假设置信度为γ，基于二项分布置信下限计算方法得到单侧置信下限v单,l满足如下关系式：其中，表示在n次有效可靠性试验中随机出现k次故障的组合数；s202：采用数值遍历法求解上式，得到置信度为γ时，所述武器装备单次使用不发生故障的概率v的单侧置信下限v单,l。s3：获取可靠度的双侧置信区间的下限值r双,l(n)与上限值r双,u(n)：其中，v双,l为所述武器装备单次使用不发生故障的概率v的双侧置信区间的下限，v双,u为所述武器装备单次使用不发生故障的概率v的双侧置信区间的上限。进一步地，所述武器装备单次使用不发生故障的概率v的双侧置信区间的下限v双,l与上限v双,u的获取方式为：s301：假设置信度为γ，基于二项分布置信区间计算方法得到双侧置信区间的下限v双,l与上限v双,u分别满足如下关系式：其中，表示在n次有效可靠性试验中随机出现k次故障的组合数；s302：采用数值遍历法求解上式，得到置信度为γ时，所述武器装备单次使用不发生故障的概率v的双侧置信区间的下限v双,l与上限v双,u。s4：根据可靠度的点估计值可靠度的单侧置信下限值r单,l(n)以及可靠度的双侧置信区间的下限值r双,l(n)与上限值r双,u(n)评估可靠度是否满足产品可靠度研制要求。实施例二基于以上实施例，本实施例对某航母舰载机弹射器可靠度进行评估(示例需要，非真实数据)。步骤一、分析并处理可靠性试验数据①假设某航母舰载机弹射器共开展可靠性试验n＝2000次，期间发生故障z＝5次(示例需要，非真实数据)。②第1次至第5次故障发生对应的总试验次数如表3所示：表3故障发生次序12345总试验次数f1＝438f2＝981f3＝1415f4＝1792f5＝1996根据上表，将第1至第5次故障发生事件依次计为概率事件a1,a2,…,a5，并计算每次故障发生前无故障连续试验次数，如表4所示：表4步骤二、建立故障发生事件的概率模型①假设装备单次试验任务不发生故障的概率为v，发生故障的概率为w，显然w＝1-v。②则，第1次故障发生对应概率事件a1，发生概率为：第2次故障发生对应概率事件a2，发生概率为：第3次故障发生对应概率事件a3，发生概率为：第4次故障发生对应概率事件a4，发生概率为：第5次故障发生对应概率事件a5，发生概率为：第5次故障后继续进行4次无故障试验对应概率事件ao，发生概率为：步骤三、构建极大似然函数步骤四、求解似然函数①根据l′(v)＝0得到等价公式(n-z)(1-v)-z×v＝0，求解弹射器单次使用不发生故障的概率v的极大似然估计为：②根据公式求解弹射器单次使用不发生故障的概率的单侧置信下限(置信度取80％)。将试验数据代入公式得：利用数值遍历法解上述方程，得到置信度为0.8时的单侧置信下限v单,l估计为：v单,l＝0.99664180(27)③根据公式求解弹射器单次使用不发生故障的概率的双侧置信区间(置信度取80％)。将试验数据代入公式得：利用遍历法解上述方程，得到置信度80％时的双侧置信区间[vl,vu]估计为：步骤五、可靠性参数估计①可靠度点估计弹射器在作战任务中连续使用n次(n≥1)的可靠度点估计为(置信度取80％)：在n取值为1,2,…,3000范围内，可靠度点估计随连续使用次数n的变化趋势如图2所示。部分典型连续使用次数n对应的可靠度如表5所示：表5②可靠度单侧置信下限估计弹射器在作战任务中连续使用n次(n≥1)的可靠度单侧置信下限估计为(置信度取80％)：在n取值为1,2,…,3000范围内，可靠度单侧置信下限估计r单,l(n)随连续使用次数n的变化趋势如图2所示。部分典型连续使用次数n对应的可靠度见表6：表6③可靠度双侧置信区间估计弹射器在作战任务中连续使用n次(n≥1)的可靠度双侧置信区间估计为(置信度取80％)：在n取值为1,2,…,3000范围内，可靠度双侧置信区间估计[r双,l(n),r双,u(n)]随连续使用次数n的变化趋势如图2所示，图2中可靠度函数为使用次数n的离散函数，n取正整数。部分典型连续使用次数n对应的可靠度如表7所示：表7由此可见，本实施例首先分析按次使用武器装备的可靠性试验数据，确定数据的有效性，掌握武器装备累计试验次数、总故障次数、故障分别发生在第几次试验等数据。然后，基于几何分布，即一种无记忆的离散随机变量概率分布，其在离散随机变量空间与指数分布在连续随机变量空间具有基本相同的特征，建立按次使用武器装备故障发生事件的概率模型。其次，以故障发生的概率事件依据，基于几何分布模型建立极大似然函数。最后，通过求解极大似然函数，得到按次使用武器装备的可靠度的评估结果。当然，本发明还可有其他多种实施例，在不背离本发明精神及其实质的情况下，熟悉本领域的技术人员当然可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。当前第1页12