CN109460538B - 一种舰载机弹射器的平均故障间隔次数评估方法 - Google Patents
一种舰载机弹射器的平均故障间隔次数评估方法 Download PDFInfo
- Publication number
- CN109460538B CN109460538B CN201811155878.7A CN201811155878A CN109460538B CN 109460538 B CN109460538 B CN 109460538B CN 201811155878 A CN201811155878 A CN 201811155878A CN 109460538 B CN109460538 B CN 109460538B
- Authority
- CN
- China
- Prior art keywords
- carrier
- double
- probability
- lower limit
- based aircraft
- 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.)
- Active
Links
- 238000011156 evaluation Methods 0.000 title claims description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 238000012360 testing method Methods 0.000 claims description 43
- 238000000034 method Methods 0.000 claims description 26
- 238000007476 Maximum Likelihood Methods 0.000 claims description 21
- 238000011161 development Methods 0.000 claims description 3
- 238000013210 evaluation model Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Complex Calculations (AREA)
Abstract
本发明提供一种按次使用武器装备的平均故障间隔次数评估方法,给出平均故障间隔次数的点估计值、平均故障间隔次数的单侧置信下限值、平均故障间隔次数的双侧置信区间的下限值以及上限值的计算方法,得到的点估计值、单侧置信下限值、双侧置信区间的下限值以及上限值能够全面评估按次使用武器装备的平均故障间隔次数,从而解决诸如舰载机弹射器等按次使用武器装备的平均故障间隔次数无法准确建模、物理意义不匹配等问题,能够建立按次使用武器装备完备的平均故障间隔次数评估模型,提高评估结果的准确度。
Description
技术领域
本发明属于武器装备可靠性工程技术领域,尤其涉及一种舰载机弹射器的平均故障间隔次数评估方法。
背景技术
针对按次使用的武器装备,现有的解决方案是忽略因装备使用过程的离散性导致其故障概率分布离散化这一现象,依然按照GJB 899A-2009《可靠性鉴定和验收试验》等标准、文献中给出的基于指数分布假设的方法计算按次使用武器装备的MCBF(平均故障间隔次数)。
指数分布是一种连续概率分布。通常连续时间工作的电子产品(如舰用计算机、舰用通信设备等)及大型复杂可修装备(如舰船电力系统、动力系统等)的故障时间服从指数分布,其故障时间可以是装备累计工作时间轴上的任意时刻(如101.032小时、1000.43小时等)。但对于诸如舰载机弹射器、舰炮等按次使用的武器装备,其故障发生的时间刻度以次数计,为正整数(如100次、400次、10000次等)。因此,目前工程中利用基于指数分布的方法评估按次使用武器装备的可靠性存在以下三个方面的不足:一是理论上不完备,指数分布为连续变量分布函数不能准确描述按次使用武器装备非连续的故障发生次数变量;二是物理意义不匹配,指数分布所蕴含的概率意义与按次使用武器装备故障发生的概率事件不匹配;三是计算结果不准确,由于存在理论上的不完备和物理意义的不匹配,基于指数分布假设的方法给出的评估结果为近似结果。
发明内容
为解决上述问题,本发明提供一种舰载机弹射器的平均故障间隔次数评估方法,能够建立舰载机弹射器完备的平均故障间隔次数评估模型,提高评估结果的准确度。
一种舰载机弹射器的平均故障间隔次数评估方法,包括以下步骤:
获取平均故障间隔次数的单侧置信下限值MCBF单,L:
其中,V单,L为所述舰载机弹射器单次使用不发生故障的概率V的单侧置信下限;
获取平均故障间隔次数的双侧置信区间的下限值MCBF双,L与上限值MCBF双,U:
其中,V双,L为所述舰载机弹射器单次使用不发生故障的概率V的双侧置信区间的下限,V双,U为所述舰载机弹射器单次使用不发生故障的概率V的双侧置信区间的上限;
根据所述平均故障间隔次数的点估计值平均故障间隔次数的单侧置信下限值MCBF单,L以及平均故障间隔次数的双侧置信区间的下限值MCBF双,L与上限值MCBF双,U评估平均故障间隔次数是否满足产品MCBF研制要求。
构造极大似然函数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单,L的获取方式为:
假设置信度为γ,基于二项分布置信下限计算方法得到单侧置信下限V单,L满足如下关系式:
采用数值遍历法求解上式,得到置信度为γ时,所述舰载机弹射器单次使用不发生故障的概率V的单侧置信下限V单,L。
进一步地,所述舰载机弹射器单次使用不发生故障的概率V的双侧置信区间的下限V双,L与上限V双,U的获取方式为:
假设置信度为γ,基于二项分布置信区间计算方法得到双侧置信区间的下限V双,L与上限V双,U分别满足如下关系式:
采用数值遍历法求解上式,得到置信度为γ时,所述舰载机弹射器单次使用不发生故障的概率V的双侧置信区间的下限V双,L与上限V双,U。
有益效果:
本发明提供一种舰载机弹射器的平均故障间隔次数评估方法,给出平均故障间隔次数的点估计值、平均故障间隔次数的单侧置信下限值、平均故障间隔次数的双侧置信区间的下限值以及上限值的计算方法,得到的点估计值、单侧置信下限值、双侧置信区间的下限值以及上限值,能够全面评估舰载机弹射器的平均故障间隔次数,从而解决诸如舰载机弹射器等舰载机弹射器的平均故障间隔次数无法准确建模、物理意义不匹配等问题,能够建立舰载机弹射器完备的平均故障间隔次数评估模型,提高评估结果的准确度。
附图说明
图1为本发明提供的一种舰载机弹射器的平均故障间隔次数评估方法的流程图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述。
实施例一
参见图1,该图为本实施例提供的一种舰载机弹射器的平均故障间隔次数评估方法的流程图。一种舰载机弹射器的平均故障间隔次数评估方法,包括以下步骤:
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
其中,{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)
S2:获取平均故障间隔次数的单侧置信区间下限值MCBF单,L:
其中,V单,L为所述舰载机弹射器单次使用不发生故障的概率V的单侧置信下限。
进一步地,所述舰载机弹射器单次使用不发生故障的概率V的单侧置信下限V单,L的获取方式为:
S201:假设置信度为γ,基于二项分布置信下限计算方法得到单侧置信下限V单,L满足如下关系式:
S202:采用数值遍历法求解上式,得到置信度为γ时,所述舰载机弹射器单次使用不发生故障的概率V的单侧置信下限V单,L。
S3:获取平均故障间隔次数的双侧置信区间的下限值MCBF双,L与上限值MCBF双,U:
其中,V双,L为所述舰载机弹射器单次使用不发生故障的概率V的双侧置信区间的下限,V双,U为所述舰载机弹射器单次使用不发生故障的概率V的双侧置信区间的上限。
进一步地,所述舰载机弹射器单次使用不发生故障的概率V的双侧置信区间的下限V双,L与上限V双,U的获取方式为:
S301:假设置信度为γ,基于二项分布置信区间计算方法得到双侧置信区间的下限V双,L与上限V双,U分别满足如下关系式:
S302:采用数值遍历法求解上式,得到置信度为γ时,所述舰载机弹射器单次使用不发生故障的概率V的双侧置信区间的下限V双,L与上限V双,U。
S4:根据所述平均故障间隔次数的点估计值平均故障间隔次数的单侧置信下限值MCBF单,L以及平均故障间隔次数的双侧置信区间的下限值MCBF双,L与上限值MCBF双,U评估平均故障间隔次数是否满足产品MCBF研制要求。
实施例二
基于以上实施例,本实施例对某航母舰载机弹射器的平均故障间隔次数进行评估(示例需要,非真实数据)。
步骤一、分析并处理可靠性试验数据
①假设某航母舰载机弹射器共开展可靠性试验N=2000次,期间发生故障Z=5次(示例需要,非真实数据)。
②第1次至第5次故障发生对应的总试验次数如表3所示:
表3
根据上表,将第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的极大似然估计为:
利用数值遍历法解上述方程,得到置信度为0.8时的单侧置信下限V单,L估计为:
V单,L=0.99664180 (27)
利用遍历法解上述方程,得到置信度80%时的双侧置信区间[VL,VU]估计为:
步骤五、可靠性参数估计
①MCBF点估计
基于可靠性试验结果,弹射器MCBF点估计为:
②MCBF单侧置信下限估计
基于可靠性试验结果,弹射器MCBF单侧置信下限估计值为(置信度取80%):
③MCBF双侧置信区间估计
基于可靠性试验结果,弹射器MCBF双侧置信区间估计值为(置信度取80%):
由此可见,本实施例首先分析舰载机弹射器的可靠性试验数据,确定数据的有效性,掌握舰载机弹射器累计试验次数、总故障次数、故障分别发生在第几次试验等数据。然后,基于几何分布,即一种无记忆的离散随机变量概率分布,其在离散随机变量空间与指数分布在连续随机变量空间具有基本相同的特征,建立舰载机弹射器故障发生事件的概率模型。其次,以故障发生的概率事件依据,基于几何分布模型建立极大似然函数。最后,通过求解极大似然函数,得到舰载机弹射器的平均故障间隔次数的评估结果。
当然,本发明还可有其他多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员当然可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (4)
1.一种舰载机弹射器的平均故障间隔次数评估方法,其特征在于,包括以下步骤:
获取平均故障间隔次数的单侧置信下限值MCBF单,L:
其中,V单,L为所述舰载机弹射器单次使用不发生故障的概率V的单侧置信下限;
获取平均故障间隔次数的双侧置信区间的下限值MCBF双,L与上限值MCBF双,U:
其中,V双,L为所述舰载机弹射器单次使用不发生故障的概率V的双侧置信区间的下限,V双,U为所述舰载机弹射器单次使用不发生故障的概率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
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811155878.7A CN109460538B (zh) | 2018-09-30 | 2018-09-30 | 一种舰载机弹射器的平均故障间隔次数评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811155878.7A CN109460538B (zh) | 2018-09-30 | 2018-09-30 | 一种舰载机弹射器的平均故障间隔次数评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109460538A CN109460538A (zh) | 2019-03-12 |
CN109460538B true CN109460538B (zh) | 2023-06-16 |
Family
ID=65607230
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811155878.7A Active CN109460538B (zh) | 2018-09-30 | 2018-09-30 | 一种舰载机弹射器的平均故障间隔次数评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109460538B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7050935B1 (en) * | 1999-03-08 | 2006-05-23 | Bombardier Transportation Gmbh | Method for assessing the reliability of technical systems |
CN103838619A (zh) * | 2014-03-17 | 2014-06-04 | 宋佰超 | 一种可修系统故障次数的确定方法 |
CN104614195A (zh) * | 2015-02-05 | 2015-05-13 | 中国人民解放军海军工程大学 | 一种纵向抽样型可靠性鉴定试验方法 |
CN105787247A (zh) * | 2014-12-25 | 2016-07-20 | 北京电子工程总体研究所 | 一种指挥控制系统的可靠性评估方法 |
CN106202938A (zh) * | 2016-07-14 | 2016-12-07 | 熊猫电子集团有限公司 | 一种通信装备可靠性评估方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7315805B2 (en) * | 2004-02-05 | 2008-01-01 | Raytheon Company | Operations and support discrete event stimulation system and method |
-
2018
- 2018-09-30 CN CN201811155878.7A patent/CN109460538B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7050935B1 (en) * | 1999-03-08 | 2006-05-23 | Bombardier Transportation Gmbh | Method for assessing the reliability of technical systems |
CN103838619A (zh) * | 2014-03-17 | 2014-06-04 | 宋佰超 | 一种可修系统故障次数的确定方法 |
CN105787247A (zh) * | 2014-12-25 | 2016-07-20 | 北京电子工程总体研究所 | 一种指挥控制系统的可靠性评估方法 |
CN104614195A (zh) * | 2015-02-05 | 2015-05-13 | 中国人民解放军海军工程大学 | 一种纵向抽样型可靠性鉴定试验方法 |
CN106202938A (zh) * | 2016-07-14 | 2016-12-07 | 熊猫电子集团有限公司 | 一种通信装备可靠性评估方法 |
Non-Patent Citations (2)
Title |
---|
杨立健 ; 赵敏 ; 甘祖旺 ; .一种装载可靠度评估方法.水下无人系统学报.2017,(第06期),92-95. * |
王宇勇 ; 刘立明 ; 石宪华 ; .某复杂武器系统的L-M系统可靠性综合评估及应用.水雷战与舰船防护.2013,(第03期),79-82. * |
Also Published As
Publication number | Publication date |
---|---|
CN109460538A (zh) | 2019-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109376407B (zh) | 一种舰载机弹射器的可靠度评估方法 | |
CN103471591B (zh) | 基于逻辑法、全局最近邻和目标航向信息的机动多目标数据互联方法 | |
CN107391335B (zh) | 一种用于检查集群健康状态的方法和设备 | |
CN105445624A (zh) | 结合小波变换及曲线拟合的电缆故障定位方法 | |
CN102207998B (zh) | 基于无失效数据的评估产品寿命及其可靠性的方法 | |
CN111444035A (zh) | 基于故障传播网络的复杂电子系统故障诊断方法 | |
CN104156612B (zh) | 基于粒子滤波的正向与逆向预测误差的故障预报方法 | |
CN104392086B (zh) | 一种基于皮尔逊秩次变量相关系数的信号检测电路及方法 | |
CN104502922A (zh) | 一种神经网络辅助粒子滤波的gps接收机自主完好性监测方法 | |
CN113762525A (zh) | 一种具有差分隐私保护的联邦学习模型训练方法 | |
CN109460538B (zh) | 一种舰载机弹射器的平均故障间隔次数评估方法 | |
WO2019019429A1 (zh) | 一种虚拟机异常检测方法、装置、设备及存储介质 | |
Sun et al. | Remaining useful life prediction for bivariate deteriorating systems under dynamic operational conditions | |
CN106410791A (zh) | 一种不良数据的前推回代追踪方法 | |
Ansell et al. | Practical reliability data analysis | |
CN117590172A (zh) | 应用于变压器的局部放电声电联合定位方法、装置和设备 | |
CN112884353B (zh) | 一种任务可靠度评价系统 | |
WO2014024640A1 (ja) | ロジック図面誤り箇所推定装置及びその方法 | |
Shi et al. | Statistical inference for a hybrid system model with incomplete observed data under adaptive progressive hybrid censoring | |
CN107562695B (zh) | 一种配电变压器的负载数据处理方法及装置 | |
CN111965609A (zh) | 雷达可靠度评估方法、装置、电子设备及可读存储介质 | |
Mosin | Quality improvement of analog circuits fault diagnosis based on ANN using clusterization as preprocessing | |
CN103258101A (zh) | 模拟电路焊点机械蠕变失效机理的表征参数确定方法 | |
CN113205215A (zh) | 一种基于知识的战场态势预测方法 | |
CN109459629B (zh) | 一种基于恢复率的恢复能力评估方法 |
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 |