CN109948116A - 一种成败型和联系统可靠度置信下限计算方法 - Google Patents
一种成败型和联系统可靠度置信下限计算方法 Download PDFInfo
- Publication number
- CN109948116A CN109948116A CN201910189503.0A CN201910189503A CN109948116A CN 109948116 A CN109948116 A CN 109948116A CN 201910189503 A CN201910189503 A CN 201910189503A CN 109948116 A CN109948116 A CN 109948116A
- Authority
- CN
- China
- Prior art keywords
- success
- failure type
- connection system
- reliability
- lower limit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 19
- 238000012360 testing method Methods 0.000 claims abstract description 79
- 230000004927 fusion Effects 0.000 claims abstract description 6
- 238000004519 manufacturing process Methods 0.000 claims description 14
- 230000001186 cumulative effect Effects 0.000 claims description 10
- 238000000034 method Methods 0.000 description 6
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004513 sizing Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种成败型和联系统可靠度置信下限计算方法,通过引入单元产品可靠性试验数据转换为成败型和联系统等效可靠性试验数据的折合系数,将各单元产品融合等效为成败型和联系统后再计算可靠度置信下限,能够解决现有计算方法理论不严谨、计算结果存在误差等问题,能够提高成败型和联系统可靠度置信下限的计算准确度。
Description
技术领域
本发明属于军事装备可靠性工程技术领域,具体涉及一种成败型和联系统可靠度置信下限计算方法。
背景技术
装备鉴定定型时通常要求武器装备以规定的置信度达到规定的可靠度要求,即要求某置信度下可靠度置信下限达到研制要求规定的可靠度指标。成败型和联系统是武器装备特有的一种可靠性逻辑关系。一般情况下难以针对成败型和联系统直接开展可靠性试验,而是针对成败型和联系统的单元产品开展可靠性试验以间接验证成败型和联系统的可靠度。发射类武器装备作战单元(如鱼雷、导弹、制导炸弹、舰炮等)为可靠性成败型产品,在开展其可靠性试验时,不同单元产品的可靠性试验次数、故障次数均存在差异,如何利用这些单元产品的可靠性试验数据计算成败型和联系统的可靠度置信下限一直是可靠性领域较难解决的问题。
目前,常用的方法是先基于不同单元产品可靠性试验数据计算各自可靠度置信下限,然后根据成败型和联系统数学表达式计算成败型和联系统的可靠度置信下限。该方法存在的问题是:由不同单元产品可靠度置信下限加权求和计算成败型和联系统可靠度置信下限在数学理论上并不完全成立,是一种近似算法,计算结果存在误差。
发明内容
有鉴于此,本发明提供了成败型和联系统可靠度置信下限计算方法,通过引入单元产品可靠性试验数据转换为成败型和联系统等效可靠性试验数据的折合系数,将各单元产品融合等效为成败型和联系统后再计算可靠度置信下限,能够有效降低计算的误差。
本发明提供的一种成败型和联系统可靠度置信下限计算方法,针对构成成败型和联系统的每个单元产品,计算其融合等效为成败型和联系统的折合系数;
基于各单元产品的可靠性试验数据,利用所述折合系数,实现对成败型和联系统可靠度置信下限计算。
进一步地,所述折合系数为:单元产品不可靠度点估计与成败型和联系统不可靠度点估计之比。
进一步地,所述单元产品不可靠度点估计为:1≤i≤N,其中,N为组成成败型和联系统的单元产品的总数,N≥2;i为单元产品的编号;表示第i个单元产品的不可靠度点估计,ni>0为第i个单元产品可靠性试验样本累计试验次数,ri>0为第i个单元产品可靠性试验样本累计故障次数;
所述成败型和联系统不可靠度点估计为:其中,αi表示第i个单元产品的和联系数。
进一步地,所述对成败型和联系统可靠度置信下限计算为:根据各单元产品可靠性试验样本累计试验次数ni,计算成败型和联系统等效可靠性试验累计次数n和;根据各单元产品可靠性试验样本累计故障次数ri及所述折合系数,计算成败型和联系统等效故障次数r和;基于所述等效可靠性试验累计次数n和和等效故障次数r和,计算成败型和联系统可靠度置信下限。
进一步地,所述成败型和联系统等效可靠性试验累计次数n和为:
所述成败型和联系统等效故障次数为r和为:其中,βi为第i个单元产品可靠性试验数据转换为成败型和联系统等效可靠性试验数据的折合系数,
进一步地,所述计算成败型和联系统可靠度置信下限:
由于引入了折合系数,r和会存在为整数或非整数两种情况:
当r和为整数时,采用公式(1)计算成败型和联系统在置信度为γ时的可靠度置信下限R和,γ,
式中R和,γ为成败型和联系统在置信度为γ(0<γ<1)时的可靠度置信下限,表示在n和次有效可靠性试验中随机发生r和次故障的组合数;
当r和为非整数时,将r和向负无穷方向取最大整数表示为将r和向正无穷方向取最小整数表示为将n和和代入公式(2)中,计算
将n和和代入公式(3)中,计算
再将计算得到和按插值法计算成败型和联系统可靠度置信下限,
进一步地,所述R和,γ、及均采用数值遍历法进行求解。
有益效果:
本发明根据提出的单元产品可靠性试验数据转换为成败型和联系统等效可靠性试验数据的折合系数,将各单元产品融合等效为成败型和联系统,计算成败型和联系统等效可靠性试验累计次数、成败型和联系统等效故障次数,再计算等效的成败型和联系统的可靠度置信下限,能够解决现有计算方法理论不严谨、计算结果存在误差等问题;能够提高成败型和联系统可靠度置信下限的计算准确度。
附图说明
图1为本发明提供的成败型和联系统可靠度置信下限计算方法流程图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种成败型和联系统可靠度置信下限计算方法,其基本思想是:首先,确定成败型和联系统可靠性模型及置信度,统计各单元产品的可靠性试验数据,计算各单元产品不可靠度点估计及成败型和联系统不可靠度点估计;然后,计算单元产品可靠性试验数据转换为成败型和联系统等效可靠性试验数据的折合系数;最后,基于折合系数计算成败型和联系统等效可靠性试验数据,并计算成败型和联系统可靠度置信下限。通过采用本发明提出的方法,能够有效地降低计算误差。
定义单元产品可靠性试验数据转换为成败型和联系统等效可靠性试验数据的折合系数为:单元产品不可靠度点估计与成败型和联系统不可靠度点估计之比,其计算过程如下:
步骤1、确定成败型和联系统可靠性模型及置信度
假设某成败型和联系统由N(N≥2)个单元产品组成,各单元产品的和联系数αi,其中,i为单元产品的个数,1≤i≤N,0<αi<1,成败型和联系统可靠性理论数学模型为公式(1):
式中R和表示成败型和联系统可靠度的真值,F和表示成败型和联系统不可靠度真值,αi表示第i个单元产品的和联系数,Ri表示第i个单元产品的可靠度真值,N表示单元产品数量。但是,在工程实际中单元产品可靠度Ri与成败型和联系统可靠度R和的真值通常是不可知的,只能基于样本可靠性试验数据得到其估计值,式(1)对于Ri与R和的点估计成立,但对其置信下限不成立。
根据装备研制要求,确定成败型和联系统可靠度计算结果的置信度γ,0<γ<1。
步骤2、统计成败型单元产品可靠性试验数据
统计成败型和联系统中各单元产品的可靠性试验数据,表示为ni,ri,1≤i≤N,其中ni>0为第i个单元产品可靠性试验样本累计试验次数,ri>0为第i个单元产品可靠性试验样本累计故障次数。
步骤3、计算各单元产品不可靠度点估计及成败型和联系统不可靠度点估计
根据步骤2统计的各单元产品可靠性试验数据ni,ri,1≤i≤N,根据公式(2)计算各单元产品可靠度点估计:
式中,表示第i个单元产品可靠度的点估计。
根据公式(3)计算各单元产品不可靠度点估计为:
式中,表示第i个单元产品不可靠度点估计。
根据公式(1),计算成败型和联系统不可靠度的点估计为:
式中,表示成败型和联系统不可靠度的点估计。
步骤4、根据公式(5)计算单元产品可靠性试验数据转换为成败型和联系统等效可靠性试验数据的折合系数,与现有技术中先针对每个单元产品计算可靠度置信下限,再根据单元产品的可靠度置信下限计算成败型和联系统的置信度下限相比,本发明通过引入折合系数先将单元产品融合等效为成败型和联系统,再针对成败型和联系统计算置信度下限,提高了成败型和联系统置信度下限计算的准确性:
式中,βi为第i个单元产品可靠性试验数据转换为成败型和联系统等效可靠性试验数据的折合系数。
基于上述步骤,根据各单元产品可靠性试验样本累计试验次数ni,计算成败型和联系统等效可靠性试验累计次数n和;根据各单元产品可靠性试验样本累计故障次数ri及所述折合系数,计算成败型和联系统等效故障次数r和;基于所述等效可靠性试验累计次数n和和等效故障次数r和,计算成败型和联系统可靠度置信下限,步骤如下:
步骤5、将步骤2中获取的各单元产品的可靠性试验数据ni,ri,1≤i≤N,按公式(6)和公式(7)转换为成败型和联系统的等效可靠性试验数据。
其中,根据公式(6)计算出成败型和联系统等效可靠性试验累计次数n和:
根据公式(7)计算出成败型和联系统等效故障次数为r和:
步骤6、根据步骤5中计算得到的成败型和联系统等效可靠性试验数据,计算成败型和联系统可靠度置信下限,由于引入了折合系数,r和会存在为整数或非整数两种情况:
当r和为整数时,将n和和r和代入如下公式,计算成败型和联系统在置信度为γ时的可靠度置信下限R和,γ,
式中R和,γ为成败型和联系统在置信度为γ(0<γ<1)时的可靠度置信下限,n和由公式(6)计算得到,r和由公式(7)计算得到,表示在n和次有效可靠性试验中随机发生r和次故障的组合数。公式(8)可以采用数值遍历法进行求解。
当r和为非整数时:将向负无穷方向取最大整数表示为r和 -,将向正无穷方向取最小整数表示为r和 +。
将n和和r和 -代入公式(9)中,计算R- 和,γ,
公式(9)中R- 和,γ可以采用数值遍历法进行求解。
将n和和r和 +代入公式(10)中,计算R+ 和,γ,
公式(10)中R+ 和,γ可以采用数值遍历法进行求解。
由公式(9)和公式(10)计算得到R- 和,γ和R+ 和,γ后,根据公式(11)按插值法计算成败型和联系统可靠度置信下限。
步骤1.1、确定成败型和联系统可靠性模型及置信度
下面以某防空系统由3个不同的成败型作战单元组成可靠性成败型和联系统为例对本发明进行详细说明。其中,单元1、单元2、单元3的和联系数依次为α1=0.6、α2=0.3、α3=0.2,根据成败型和联系统可靠性理论数学模型公式(1),可得:
R和=1-F和=0.6R1+0.2R2+0.2R3 (12)
假设该成败型和联系统的置信度为γ=90%。
步骤1.2、统计成败型单元产品可靠性试验数据
成败型和联系统中单元1、单元2、单元3的可靠性试验数据依次为n1=17、r1=2;n2=15、r2=1;n3=21、r3=4。其中n1、n2、n3表示可靠性试验次数,r1、r2、r3表示试验期间累计故障次数。
步骤1.3、计算各单元产品不可靠度点估计及成败型和联系统不可靠度点估计
将步骤1.2中各单元可靠性试验数据代入公式(2),得到各单元的可靠度点估计为式(13):
式中分别表示单元1、单元2、单元3的可靠度点估计。
将步骤1.2中各单元可靠性试验数据代入公式(3),得到各单元的不可靠度点计值为式(14):
式中分别表示单元1、单元2、单元3的不可靠度点估计。
根据公式(4),得到成败型和联系统的不可靠度点估计值为,
步骤1.4、计算单元产品可靠性试验数据转换为成败型和联系统等效可靠性试验数据的折合系数
将式(14)、式(15)的计算结果代入式(5),得到单元1、单元2、单元3可靠性试验数据转换为成败型和联系统等效可靠性试验数据的折合系数为,
步骤1.5、计算成败型和联系统等效可靠性试验数据
将步骤1.2中各单元的可靠性试验数据n1=17、r1=2;n2=15、r2=1;n3=21、r3=4代入公式(6)和公式(7),得到:
成败型和联系统等效可靠性试验次数n和为:
成败型和联系统等效故障次数为r和为:
步骤1.6、计算成败型和联系统可靠度置信下限
根据步骤1.1,置信度取γ=90%,根据式(18)r和=6.4669,为非整数。r和=6.4669向负无穷方向取最大整数为向正无穷方向取最小整数为
将n和=53,代入公式(9),得到
利用数值遍历法求解上式得到
将n和=53,代入公式(10),得到
利用数值遍历法求解上式得到
将代入公式(11),得到该示例成败型和联系统可靠度置信下限,
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种成败型和联系统可靠度置信下限计算方法,其特征在于:
针对构成成败型和联系统的每个单元产品,计算其融合等效为成败型和联系统的折合系数;
基于各单元产品的可靠性试验数据,利用所述折合系数,实现对成败型和联系统可靠度置信下限计算。
2.根据权利要求1所述的一种成败型和联系统可靠度置信下限计算方法,其特征在于,所述折合系数为:单元产品不可靠度点估计与成败型和联系统不可靠度点估计之比。
3.根据权利要求2所述的一种成败型和联系统可靠度置信下限计算方法,其特征在于:
所述单元产品不可靠度点估计为:1≤i≤N,其中,N为组成成败型和联系统的单元产品的总数,N≥2;i为单元产品的编号;表示第i个单元产品的不可靠度点估计,ni>0为第i个单元产品可靠性试验样本累计试验次数,ri>0为第i个单元产品可靠性试验样本累计故障次数;
所述成败型和联系统不可靠度点估计为:其中,αi表示第i个单元产品的和联系数。
4.根据权利要求1或2或3所述的一种成败型和联系统可靠度置信下限计算方法,其特征在于:
所述对成败型和联系统可靠度置信下限计算为:根据各单元产品可靠性试验样本累计试验次数ni,计算成败型和联系统等效可靠性试验累计次数n和;根据各单元产品可靠性试验样本累计故障次数ri及所述折合系数,计算成败型和联系统等效故障次数r和;基于所述等效可靠性试验累计次数n和和等效故障次数r和,计算成败型和联系统可靠度置信下限。
5.根据权利要求4所述的一种成败型和联系统可靠度置信下限计算方法,其特征在于:
所述成败型和联系统等效可靠性试验累计次数n和为:
所述成败型和联系统等效故障次数为r和为:其中,βi为第i个单元产品可靠性试验数据转换为成败型和联系统等效可靠性试验数据的折合系数,
6.根据权利要求5所述的一种成败型和联系统可靠度置信下限计算方法,其特征在于:
所述计算成败型和联系统可靠度置信下限:
由于引入了折合系数,r和会存在为整数或非整数两种情况:
当r和为整数时,采用公式(1)计算成败型和联系统在置信度为γ时的可靠度置信下限R和,γ,
式中R和,γ为成败型和联系统在置信度为γ(0<γ<1)时的可靠度置信下限,表示在n和次有效可靠性试验中随机发生r和次故障的组合数;
当r和为非整数时,将r和向负无穷方向取最大整数表示为将r和向正无穷方向取最小整数表示为将n和和代入公式(2)中,计算
将n和和代入公式(3)中,计算
再将计算得到和按插值法计算成败型和联系统可靠度置信下限,
7.根据权利要求6所述的系统,其特征在于,所述R和,γ、及均采用数值遍历法进行求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910189503.0A CN109948116B (zh) | 2019-03-13 | 2019-03-13 | 一种成败型和联系统可靠度置信下限计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910189503.0A CN109948116B (zh) | 2019-03-13 | 2019-03-13 | 一种成败型和联系统可靠度置信下限计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109948116A true CN109948116A (zh) | 2019-06-28 |
CN109948116B CN109948116B (zh) | 2024-01-30 |
Family
ID=67008765
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910189503.0A Active CN109948116B (zh) | 2019-03-13 | 2019-03-13 | 一种成败型和联系统可靠度置信下限计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109948116B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012173096A (ja) * | 2011-02-21 | 2012-09-10 | Tokyo Electric Power Co Inc:The | 電力機器の保全システム、電力機器の保全プログラム、および電力機器の保全方法 |
CN107132834A (zh) * | 2017-04-12 | 2017-09-05 | 北京航空航天大学 | 基于统计量排序的和联系统可靠度精确置信下限估计方法 |
CN107526715A (zh) * | 2017-07-26 | 2017-12-29 | 武汉船舶通信研究所(中国船舶重工集团公司第七二二研究所) | 一种可靠性评估方法和装置 |
CN109376407A (zh) * | 2018-09-30 | 2019-02-22 | 中国人民解放军92942部队 | 一种按次使用武器装备的可靠度评估方法 |
-
2019
- 2019-03-13 CN CN201910189503.0A patent/CN109948116B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012173096A (ja) * | 2011-02-21 | 2012-09-10 | Tokyo Electric Power Co Inc:The | 電力機器の保全システム、電力機器の保全プログラム、および電力機器の保全方法 |
CN107132834A (zh) * | 2017-04-12 | 2017-09-05 | 北京航空航天大学 | 基于统计量排序的和联系统可靠度精确置信下限估计方法 |
CN107526715A (zh) * | 2017-07-26 | 2017-12-29 | 武汉船舶通信研究所(中国船舶重工集团公司第七二二研究所) | 一种可靠性评估方法和装置 |
CN109376407A (zh) * | 2018-09-30 | 2019-02-22 | 中国人民解放军92942部队 | 一种按次使用武器装备的可靠度评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109948116B (zh) | 2024-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mendershausen | On the significance of Professor Douglas' production function | |
Garrido et al. | A new look at Horn’s parallel analysis with ordinal variables. | |
Maddox et al. | Value-based assessment of dod acquisitions programs | |
Triantafyllou et al. | Signature and IFR preservation of 2-within-consecutive $ k $-out-of-$ n $: $ F $ systems | |
CN108985373B (zh) | 一种多传感器数据加权融合方法 | |
Ewens | The pseudo-transient distribution and its uses in genetics | |
CN108132914A (zh) | 有限数量产品性能失效率的高置信度统计推断和确认方法 | |
CN111428793A (zh) | 一种基于改进证据相异度的证据融合方法 | |
CN106056098A (zh) | 一种基于类别合并的脉冲信号聚类分选方法 | |
Holzman | Politics and guesswork: CIA and DIA estimates of Soviet military spending | |
CN109948116A (zh) | 一种成败型和联系统可靠度置信下限计算方法 | |
CN109376407A (zh) | 一种按次使用武器装备的可靠度评估方法 | |
Zwick et al. | Application of an empirical Bayes enhancement of Mantel-Haenszel differential item functioning analysis to a computerized adaptive test | |
Murphy et al. | A comparison of item selection techniques for testlets | |
Kvålseth | Cautionary note about the Herfindahl-Hirschman index of market (industry) concentration | |
CN109684713A (zh) | 基于贝叶斯的复杂系统可靠性分析方法 | |
CN109978349A (zh) | 一种指数型和联系统可靠度置信下限计算方法 | |
Rabaiah et al. | Generalizations of power function and type-i half logistic distributions using quadratic transmutation map | |
Almawajdeh et al. | Generalization of Rani and generalized Akash distributions | |
Özönder | Long range azimuthal correlations of multiple gluons in gluon saturation limit | |
Ismail et al. | Estimation of 2-and 3-parameter Burr Type XII distribution using EM algorithm | |
Pfisterer | Algorithmic Fairness | |
CN113076661B (zh) | 一种不确定数据驱动的雷达预警探测建模方法 | |
Lifshits | On representation of Levy’s fields by indicators | |
Anggraeni et al. | Pengaruh Inflasi Dan Nilai Tukar Rupiah Terhadap Dollar Dengan Kurs Transaksi Tengah Sebagai Gejala Resesi Terhadap Perubahan Ihsg |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |