CN110705103B - 一种伽玛型通用件的备件保障概率计算方法及装置 - Google Patents

一种伽玛型通用件的备件保障概率计算方法及装置 Download PDF

Info

Publication number
CN110705103B
CN110705103B CN201910944873.0A CN201910944873A CN110705103B CN 110705103 B CN110705103 B CN 110705103B CN 201910944873 A CN201910944873 A CN 201910944873A CN 110705103 B CN110705103 B CN 110705103B
Authority
CN
China
Prior art keywords
equal
probability
equipment
parts
gamma
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
Application number
CN201910944873.0A
Other languages
English (en)
Other versions
CN110705103A (zh
Inventor
阮旻智
王俊龙
肖剑波
李华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Naval University of Engineering PLA
Original Assignee
Naval University of Engineering PLA
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Naval University of Engineering PLA filed Critical Naval University of Engineering PLA
Priority to CN201910944873.0A priority Critical patent/CN110705103B/zh
Publication of CN110705103A publication Critical patent/CN110705103A/zh
Application granted granted Critical
Publication of CN110705103B publication Critical patent/CN110705103B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Complex Calculations (AREA)

Abstract

本发明涉及一种伽玛型通用件的备件保障概率计算方法及装置,属于装备备件保障技术领域。本发明首先利用各装备通用件的故障概率建立概率矩阵,然后依次选取列装通用件的装备作为当前装备,利用卷积运算依次计算其他装备中通用件故障之和的概率分布情况,最后根据所述的其他装备中通用件故障之和的概率分布情况计算当前装备的备件保障概率。本发明充分考虑到了通用件在不同装备中的故障概率,利用其他装备的故障概率计算当前装备的保障概率,统筹规划通用件的整体需求量,在满足备件资源的充分性的同时,减少资源浪费。

Description

一种伽玛型通用件的备件保障概率计算方法及装置
技术领域
本发明涉及装备备件保障技术领域,具体涉及一种伽玛型通用件的备件保障概率计算方法及装置。
背景技术
伽玛分布常用来描述类似“冲击”引起的故障,假若单元能经受若干次外界冲击,但当单元受冲击次数累积到一定次数时就产生故障。例如电网中存在着电涌现象,一些电子器件当承受的电涌冲击次数超过一定数量时会发生故障。伽玛型通用件是在不同类型或同类型不同规格的装备中具有互换性的寿命分布服从伽玛分布的零部件单元,伽玛型通用件的寿命X服从伽玛分布记作X~Ga(a,b),其中a>0为形状参数,b>0为尺度参数,X的密度函数为
Figure BDA0002223862840000011
其中Γ(a)为伽玛函数,且
Figure BDA0002223862840000012
通用件由于具有较高的通用性,因此在综合考虑所有装备(内含通用件)的总体需求后,对其备件资源常常采取集中贮存、统一供应的保障方式。在常规的备件保障方式中,是分别计算各装备的备件需求量,然后再进行汇总。如果在通用件的备件保障中采用常规方法,则因没有考虑到通用件的互换性而造成对备件的过度采购和储备,形成浪费。在对通用件的备件进行统筹规划的过程中,站在装备使用者的立场,需要回答一个关键问题:“统一供应保障时,如何评估各装备的备件保障效果?”目前,还没有较好的方法回答该问题。
发明内容
针对伽玛型通用件,当统一供应的备件数量为S时,本发明提供一种伽玛型通用件的备件保障概率计算方法及装置,用于计算内含该类型通用件各装备的备件保障概率,从而量化对各装备的保障效果,解决背景技术中所提到的问题。
本发明解决上述技术问题的技术方案如下:
第一方面,本发明提供一种伽玛型通用件的备件保障概率计算方法,包括以下步骤:
S1,建立并计算概率矩阵dyPs,矩阵中元素dyPs(i,j)表示第i号装备的通用件发生j-1次故障的概率,1≤i≤N,1≤j≤S+2,N为列装伽玛型通用件的装备台数,S为备件仓库中伽玛型通用件的备件数量;矩阵元素dyPs(i,j)取值如下:
Figure BDA0002223862840000021
式中Ti为伽玛型通用件在第i号装备一个保障周期内的计划工作时间;
S2,令当前装备编号i=1;
S3,其余装备编号保存在数组i23中,初始化i23m=m+1,1≤m≤N-1;若i>1,则修正i23i-1=1;
S4,利用卷积运算依次计算其他装备中通用件故障之和的概率分布情况,并根据所述的其他装备中通用件故障之和的概率分布情况计算当前装备的备件保障概率dyPi
S5,令i=i+1,若i≤N则跳转至步骤S3,否则终止计算dyPi即为各装备中该通用件的备件保障概率,1≤i≤N。
进一步的,所述的步骤S4包括:
S41,令j=1;
S42,令dp1=dyPs(i,j);
S43,若j=1,则令dyPi=dp1并跳转至步骤S46,否则跳转至步骤S44;
S44,计算其他装备通用件故障之和的概率分布情况,令m=1;
S441,令i2=i23m,令数组p1等于矩阵dyPs中的第i2行向量,即p1=dyPs(i2,:);
S442,若m=1,则令p2=p1;否则计算p1和p2的卷积,并将结果赋给p2,即:令p2=p1*p2,此处*为卷积计算符号,p2记录了其他装备通用件故障之和的概率分布情况;
S443,更新m,令m=m+1,若m≤N-1则跳转至步骤S441,否则跳转至步骤S45;
S45,令
Figure BDA0002223862840000031
pr反映了其他装备通用件故障之和不大于S-j+1的概率,更新
Figure BDA0002223862840000032
式中
Figure BDA0002223862840000033
近似描述了第i号装备通用件需要j-1个备件并且得到j-1个备件的概率;
S46,更新j,令j=j+1,若j≤S+1则跳转至步骤S42,否则跳转至步骤S5。
本发明的有益效果是:由于通用型备件具有较高的通用性,因此通用性备件资源常常采取集中贮存、统一供应的保障方式。常规的备件保障方式由于没有考虑到通用件的互换性,会造成对备件的过度采购和储备,形成浪费。本发明通过卷积运算计算除当前装备外的其他装备通用件故障之和的概率分布情况,然后利用其他装备通用件故障之和的概率分布情况计算当前装备的保障概率,上述方法充分考虑到了通用件在不同装备中的故障概率,利用其他装备的故障概率计算当前装备的保障概率,统筹规划通用件的整体需求量,在满足备件资源的充分性的同时,减少资源浪费。
第二方面,本发明提供一种伽玛型通用件的备件保障概率计算装置,包括:
概率矩阵建立模块,用于建立并计算概率矩阵dyPs,矩阵中元素dyPs(i,j)表示第i号装备的通用件发生j-1次故障的概率,1≤i≤N,1≤j≤S+2,N为列装伽玛型通用件的装备台数,S为备件仓库中伽玛型通用件的备件数量;矩阵元素dyPs(i,j)取值如下:
Figure BDA0002223862840000041
式中,Ti为伽玛型通用件在第i号装备一个保障周期内的计划工作时间;
备件保障概率计算模块,用于利用卷积运算依次计算出当前装备外的其他装备中通用件故障之和的概率分布情况,并根据所述的其他装备中通用件故障之和的概率分布情况计算当前装备的备件保障概率。
第三方面,本发明提供一种伽玛型通用件的备件保障概率计算装置,包括存储器和处理器;
所述存储器,用于存储计算机软件程序;
所述处理器,用于读取并执行所述存储器中存储的计算机软件程序,实现本发明第一方面所述的一种伽玛型通用件的备件保障概率计算方法。
第四方面,本发明还提供一种非暂态计算机可读存储介质,所述存储介质中存储有用于实现本发明第一方面所述的一种伽玛型通用件的备件保障概率计算方法的计算机软件程序。
附图说明
图1为本发明实施例一提供的一种伽玛型通用件的备件保障概率计算方法流程图;
图2~图4是上述算例备件数量从0~9时,采用本发明方法和仿真方法,各装备的通用件备件保障概率结果情况。
图5为本发明实施例二提供的一种伽玛型通用件的备件保障概率计算装置结构图;
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
在对本发明做进一步解释说明之前,需要约定:某通用件分别在N台装备中使用(列装数为1),这些装备可能是不同类型或同类型不同规格。已知该通用件在各装备一个保障周期内的计划工作时间Ti,1≤i≤N;当该通用件发生故障时,采用换件维修的方式排除故障(换件修理时间忽略不计)。备件仓库中该通用件的备件数量记为S,统一对这N台装备进行保障。
实施例一
如图1所示,本发明实施例提供一种伽玛型通用件的备件保障概率计算方法,包括以下步骤:
S1,利用各装备通用件的故障概率建立概率矩阵dyPs,矩阵中元素dyPs(i,j)表示第i号装备的通用件发生j-1次故障的概率,1≤i≤N,1≤j≤S+2,N为列装伽玛型通用件的装备台数,S为备件仓库中伽玛型通用件的备件数量;矩阵元素dyPs(i,j)取值如下:
Figure BDA0002223862840000051
式中,Ti为伽玛型通用件在第i号装备一个保障周期内的计划工作时间;
S2,令当前装备编号i=1;
S3,其余装备编号保存在数组i23中,初始化i23m=m+1,1≤m≤N-1;若i>1,则修正i23i-1=1;
S4,计算当前装备的备件保障概率dyPi
S41,令j=1;
S42,令dp1=dyPs(i,j);
S43,若j=1,则令dyPi=dp1并跳转至步骤S46,否则跳转至步骤S44;
S44,计算其他装备通用件故障之和的概率分布情况,令m=1;
S441,令i2=i23m,令数组p1等于矩阵dyPs中的第i2行向量,即p1=dyPs(i2,:);
S442,若m=1,则令p2=p1;否则计算p1和p2的卷积,并将结果赋给p2,即:令p2=p1*p2,此处*为卷积计算符号,p2记录了其他装备通用件故障之和的概率分布情况;
S443,更新m,令m=m+1,若m≤N-1则跳转至步骤S441,否则跳转至步骤S45;
S45,令
Figure BDA0002223862840000061
pr反映了其他装备通用件故障之和不大于S-j+1的概率,更新
Figure BDA0002223862840000062
式中
Figure BDA0002223862840000063
近似描述了第i号装备通用件需要j-1个备件并且得到j-1个备件的概率;
S46,更新j,令j=j+1,若j≤S+1则跳转至步骤S42,否则跳转至步骤S5;
S5,令i=i+1,若i≤N则跳转至步骤S3,否则终止计算,dyPi即为各装备中该通用件的备件保障概率,1≤i≤N。
算例:
3台不同类型的装备中都有一件某伽玛型通用件,该通用件的寿命服从伽玛分布Ga(2.2,200),该通用件在3台装备中每年的工作时间预计为600、900、1200小时。该通用件的备件数量为4,试计算各装备该通用件的备件保障概率。
解:1)计算概率矩阵dyPs,结果如表1。
表1
Figure BDA0002223862840000071
2)令当前装备编号i=1。
3)其余装备编号保存在数组i23中,初始化i23m=m+1,1≤m≤N-1。若i>1,则修正i23i-1=1。
4)计算当前装备的备件保障概率dyPi
5)令i=i+1,若i≤N则转3),否则转6)。步骤2)~5)的dyPi计算结果如表2。
表2各装备的通用件备件保障概率
装备编号 1 2 3
备件保障概率 0.605 0.531 0.500
6)终止计算,输出结果为:当备件数量为4时,三台装备的该通用件备件保障概率分别为0.605、0.531和0.500。
图2~图4是上述算例备件数量从0~9时,采用本发明方法和仿真方法,1#、2#、3#装备的通用件备件保障概率结果情况。大量仿真结果表明:本发明的方法能对伽玛型通用件的备件保障效果进行准确地定性评估,当定性评估结论为中等以上时,本发明的备件保障概率误差在工程应用范围以内。
实施例二
如图2所示,本发明实施例提供一种伽玛型通用件的备件保障概率计算装置,包括:
概率矩阵建立模块,用于建立并计算概率矩阵dyPs,矩阵中元素dyPs(i,j)表示第i号装备的通用件发生j-1次故障的概率,1≤i≤N,1≤j≤S+2,N为列装伽玛型通用件的装备台数,S为备件仓库中伽玛型通用件的备件数量;矩阵元素dyPs(i,j)取值如下:
Figure BDA0002223862840000081
式中,Ti为伽玛型通用件在第i号装备一个保障周期内的计划工作时间;
备件保障概率计算模块,用于利用卷积运算依次计算出当前装备外的其他装备中通用件故障之和的概率分布情况,并根据所述的其他装备中通用件故障之和的概率分布情况计算当前装备的备件保障概率。
实施例三
本发明实施例提供一种伽玛型通用件的备件保障概率计算装置,包括存储器和处理器;
所述存储器,用于存储计算机软件程序;
所述处理器,用于读取并执行所述存储器中存储的计算机软件程序,实现如本发明实施例一所述的一种伽玛型通用件的备件保障概率计算方法。
实施例四
本发明实施例提供一种非暂态计算机可读存储介质,所述存储介质中存储有用于实现如本发明实施例一所述的一种伽玛型通用件的备件保障概率计算方法的计算机软件程序。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.一种伽玛型通用件的备件保障概率计算方法,其特征在于,包括以下步骤:
S1,建立并计算概率矩阵dyPs,矩阵中元素dyPs(i,j)表示第i号装备的通用件发生j-1次故障的概率,1≤i≤N,1≤j≤S+2,N为列装伽玛型通用件的装备台数,S为备件仓库中伽玛型通用件的备件数量;矩阵元素dyPs(i,j)取值如下:
Figure FDA0003558809830000011
式中,Ti为伽玛型通用件在第i号装备一个保障周期内的计划工作时间;
S2,令当前装备编号i=1;
S3,其余装备编号保存在数组i23中,初始化i23m=m+1,1≤m≤N-1;若i>1,则修正i23i-1=1;
S4,利用卷积运算依次计算其他装备中通用件故障之和的概率分布情况,并根据所述的其他装备中通用件故障之和的概率分布情况计算当前装备的备件保障概率dyPi
S5,令i=i+1,若i≤N则跳转至步骤S3,否则终止计算dyPi即为各装备中该通用件的备件保障概率,1≤i≤N;
所述的步骤S4包括:
S41,令j=1;
S42,令dp1=dyPs(i,j);
S43,若j=1,则令dyPi=dp1并跳转至步骤S46,否则跳转至步骤S44;
S44,令m=1;
S441,令i2=i23m,令数组p1等于矩阵dyPs中的第i2行向量,即p1=dyPs(i2,:);
S442,若m=1,则令p2=p1;否则计算p1和p2的卷积,并将结果赋给p2,即:令p2=p1*p2,此处*为卷积计算符号,p2记录了其他装备通用件故障之和的概率分布情况;
S443,更新m,令m=m+1,若m≤N-1则跳转至步骤S441,否则跳转至步骤S45;
S45,令
Figure FDA0003558809830000021
更新
Figure FDA0003558809830000022
S46,更新j,令j=j+1,若j≤S+1则跳转至步骤S42,否则跳转至步骤S5。
2.一种伽玛型通用件的备件保障概率计算装置,其特征在于,包括存储器和处理器;
所述存储器,用于存储计算机软件程序;
所述处理器,用于读取并执行所述存储器中存储的计算机软件程序,实现权利要求1所述的一种伽玛型通用件的备件保障概率计算方法。
3.一种非暂态计算机可读存储介质,其特征在于,所述存储介质中存储有用于实现权利要求1所述的一种伽玛型通用件的备件保障概率计算方法的计算机软件程序。
CN201910944873.0A 2019-09-30 2019-09-30 一种伽玛型通用件的备件保障概率计算方法及装置 Active CN110705103B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910944873.0A CN110705103B (zh) 2019-09-30 2019-09-30 一种伽玛型通用件的备件保障概率计算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910944873.0A CN110705103B (zh) 2019-09-30 2019-09-30 一种伽玛型通用件的备件保障概率计算方法及装置

Publications (2)

Publication Number Publication Date
CN110705103A CN110705103A (zh) 2020-01-17
CN110705103B true CN110705103B (zh) 2022-06-14

Family

ID=69196602

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910944873.0A Active CN110705103B (zh) 2019-09-30 2019-09-30 一种伽玛型通用件的备件保障概率计算方法及装置

Country Status (1)

Country Link
CN (1) CN110705103B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116579494B (zh) * 2023-05-23 2024-03-19 中国人民解放军海军工程大学 基于维修耗时下机电设备的备件库存预测方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105095552A (zh) * 2014-05-23 2015-11-25 中国石油化工股份有限公司 一种通用概率分布参数估计流程与分析方法
CN107072624A (zh) * 2014-06-18 2017-08-18 Elekta股份有限公司 用于自动治疗计划的系统和方法
CN108509390A (zh) * 2018-03-25 2018-09-07 中国人民解放军海军工程大学 一种伽玛型串联部件的备件需求量计算方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8832497B2 (en) * 2012-02-07 2014-09-09 A.L.D. Advanced Logistics Development Ltd Methods, apparatus and systems for performing dynamic fault tree analysis
CN108564270B (zh) * 2018-04-09 2021-11-02 中国人民解放军海军工程大学 一种贮存失效风险下的伽玛型单元备件需求量计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105095552A (zh) * 2014-05-23 2015-11-25 中国石油化工股份有限公司 一种通用概率分布参数估计流程与分析方法
CN107072624A (zh) * 2014-06-18 2017-08-18 Elekta股份有限公司 用于自动治疗计划的系统和方法
CN108509390A (zh) * 2018-03-25 2018-09-07 中国人民解放军海军工程大学 一种伽玛型串联部件的备件需求量计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于可靠性增长的武器系统可靠性鉴定试验方案研究;阮旻智等;《兵工学报》;20171231;全文 *

Also Published As

Publication number Publication date
CN110705103A (zh) 2020-01-17

Similar Documents

Publication Publication Date Title
CN106844953B (zh) 一种威布尔型有寿件备件的保障概率计算方法
US11586490B2 (en) System and method of identifying self-healing actions for computing systems using reinforcement learning
CN110717266B (zh) 一种对数正态型通用件的备件保障概率计算方法及装置
CN113821364A (zh) 内存故障的处理方法、装置、设备及存储介质
US20090271235A1 (en) Apparatus and method for generating survival curve used to calculate failure probability
CN110705103B (zh) 一种伽玛型通用件的备件保障概率计算方法及装置
Saraf et al. Generalized multi‐release modelling of software reliability growth models from the perspective of two types of imperfect debugging and change point
CN110738007B (zh) 一种电子类通用件的备件保障概率计算方法及装置
Bauer et al. Practical system reliability
CN111860645B (zh) 挥发性有机物观测数据中缺省值的修复方法及装置
CN110738008B (zh) 一种机电类通用件的备件保障概率计算方法及装置
CN110705102B (zh) 一种机械类通用件的备件保障概率计算方法及装置
CN110674044B (zh) 功能自动化测试的覆盖率获取方法、系统、设备及介质
JP2008027284A (ja) 障害処理システム、障害処理方法、障害処理装置およびプログラム
CN115686909A (zh) 内存故障的预测方法和装置、存储介质及电子装置
CN110688759A (zh) 表决部件的备件量计算方法、仿真方法、终端及存储介质
US11681576B2 (en) Anomaly coping support apparatus, method, and program
US20140250138A1 (en) Method and apparatus for optimizing log file filtering
Cha et al. An information‐based burn‐in procedure for minimally repaired items from mixed population
CN113988789A (zh) 基于流程自动筛选数据的方法、系统、设备及存储介质
CN113886116A (zh) 数据质量校验方法、装置、计算机设备和存储介质
CN109543276B (zh) 大型货轮长期保障任务多Gamma单元备件需求量确定方法
CN110717267B (zh) 一种伽玛型通用件备件需求量计算方法及装置
CN110704884B (zh) 一种对数正态型通用件备件需求量计算方法及装置
CN110489208A (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