CN108520145B - 一种贮存失效风险下的威布尔型单元备件需求量计算方法 - Google Patents
一种贮存失效风险下的威布尔型单元备件需求量计算方法 Download PDFInfo
- Publication number
- CN108520145B CN108520145B CN201810313157.8A CN201810313157A CN108520145B CN 108520145 B CN108520145 B CN 108520145B CN 201810313157 A CN201810313157 A CN 201810313157A CN 108520145 B CN108520145 B CN 108520145B
- Authority
- CN
- China
- Prior art keywords
- spare parts
- weibull
- spare part
- spare
- probability
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出了一种贮存失效风险下的威布尔型单元备件需求量计算方法,包括初始化步骤、计算保障概率和判断步骤三个步骤,计算准确、误差小,计算快速方便,与传统采用的一次保障任务的执行情况仿真模拟结果较为接近,满足工程应用要求。
Description
技术领域
本发明涉及备件需求量计算领域,尤其涉及一种贮存失效风险下的威布尔型单元备件需求量计算方法。
背景技术
准确计算备件需求量,能从备件采购费用这种经济成本的角度,定量描述装备的保障性优劣程度。
“备件在投入使用前始终是“良品”,不会贮存失效”是当前所有备件需求量计算方法的一个常见假定。如果备件贮存在环境良好的专业仓库中,则该假定是合理的,与现实中备件的实际表现也是极为相符的。但对于那些平时并非在专业仓库中贮存的备件而言,例如随装备配置的随机备件,有时受工作环境空间有限等各种条件限制,不能提供满足贮存标准的备件存放环境;尤其是对于那些长期在外场工作的装备,例如舰船上的相关装备常年处于海上,往往处于高湿、高盐度以及机电设备工作或船体摇晃带来的振动状态等较为恶劣的工作环境,如果此时的工作环境就是随机备件的贮存环境的话,则备件存在着在贮存期间失效的风险。此时,如果无视贮存失效风险,仍然采用常规方法计算备件需求量,则会因部分备件在贮存期间失效而数量不够,进而导致保障任务失败的结果。
发明内容
有鉴于此,本发明提出了一种计算准确、误差小的贮存失效风险下的威布尔型单元备件需求量计算方法。
在本文中假定:威布尔型单元备件的贮存寿命服从威布尔分布W(α1,b1),工作寿命服从威布尔分布W(α2,b2);保障任务开始时刻同时也是备件开始贮存的时刻,且贮存寿命和工作寿命二者相互独立;当保障任务时间为Tw时,要求配置一定数量的备件,使得备件保障概率Pok不得低于保障概率指标P0。
本发明的技术方案是这样实现的:本发明提供了一种贮存失效风险下的威布尔型单元备件需求量计算方法,包括以下步骤,
如果Pok≥P0,则备件需求量为0,计算终止,其中P0为保障概率指标;否则,令备件数量S=1,转步骤S1.2;
S1.2,令i=1,Snow=S;
S2.1,计算故障概率gPi,
S2.3,计算贮存效果Pzs
S2.4,令Pok=Pok+gPi×Pzs,令Snow=(Snow-1)×(1-Pz0);
S2.5,令i=i+1,若i≤S,则转步骤S2.1,否则转步骤S3;
本发明的贮存失效风险下的威布尔型单元备件需求量计算方法相对于现有技术具有以下有益效果:
(1)计算准确、误差小,计算快速方便,与传统采用的一次保障任务的执行情况仿真模拟结果较为接近,满足工程应用要求。
具体实施方式
下面将结合本发明实施方式,对本发明实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式仅仅是本发明一部分实施方式,而不是全部的实施方式。基于本发明中的实施方式,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式,都属于本发明保护的范围。
实施例1
某威布尔型单元的贮存寿命服从威布尔分布W(4000,1.4),工作寿命服从威布尔分布W(3000,1.6),保障任务时间Tw=5000h时,保障概率指标P0为0.85,试计算备件需求量步骤如下:
S1,初始化步骤
因为Pok<P0,因此令备件数量S=1,转步骤1.2。
S1.2,令i=1,Snow=S;
S2,计算保障概率Pok
S2.1,计算故障概率gP1,
令
S2.3,计算贮存效果Pzs
S2.4,令Pok=Pok+gP1×Pzs=0.3407,令Snow=(Snow-1)×(1-Pz0)=0;
S2.5,令i=i+1=2,由于i>S,因此转步骤3。
S3,判断步骤
重复以上过程,最终当S=6时,Pok=0.8913>P0满足要求,因此,备件需求量为6。
实施例2
某威布尔型单元的贮存寿命服从威布尔分布W(4000,1.4),工作寿命服从威布尔分布W(3000,1.6),保障任务时间Tw=5000h时,保障概率指标P0为0.85,试计算备件需求量步骤如下:
1)初始化单元的已工作时间simTw=0;
2)产生1个随机数t0,用于模拟装备中该单元的工作寿命,t0服从伽玛分布Ga(α2,b2);令simTw=t0;
3)产生k个随机数t1m(1≤m≤k),用于模拟备件的贮存寿命,t1m服从伽玛分布Ga(α1,b1);
4)产生k个随机数t2m(1≤m≤k),用于模拟备件的工作寿命,t2m服从伽玛分布Ga(α2,b2);
5)比较simTw与Tw的大小情况。
若simTw>Tw,则本次任务保障成功,记simFlag=1;
若simTw<Tw,则发生一次故障,在所有可用备件中寻找贮存寿命t1m大于simTw的备件,那些贮存寿命t1m不大于simTw的备件为贮存失效备件,将其从备件库中剔除。
如果不能找到贮存寿命大于simTw的备件,则本次任务保障失败,记simFlag=0;如果能找到一个贮存寿命大于simTw的备件(记其工作寿命为t2'),则令simTw=simTw+t2',然后把该备件从备件库中剔除,转5)。
按照上述流程,可多次模拟保障任务,对得到的所有仿真结果simFlag进行统计,其均值即为保障任务成功率,也是备件保障概率。
以下表1列出了实施例1中的计算过程中,S=1~6时,其对应的备件保障概率计算结果,以及实施例2中模拟结果。
表1备件保障概率的模拟结果和实施例1计算结果的对比情况
表1的结果表明,尽管步骤S2.3中的取整操作会导致计算误差,但是本发明结果仍然与模拟结果较为接近,满足工程应用要求。
以上所述仅为本发明的较佳实施方式而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种贮存失效风险下的威布尔型单元备件需求量计算方法,其特征在于:包括以下步骤,
如果Pok≥P0,则备件需求量为0,计算终止,其中P0为保障概率指标;否则,令备件数量S=1,转步骤S1.2;
S1.2,令i=1,Snow=S;
S2.1,计算故障概率gPi,
S2.3,计算贮存效果Pzs
S2.4,令Pok=Pok+gPi×Pzs,令Snow=(Snow-1)×(1-Pz0);
S2.5,令i=i+1,若i≤S,则转步骤S2.1,否则转步骤S3;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810313157.8A CN108520145B (zh) | 2018-04-09 | 2018-04-09 | 一种贮存失效风险下的威布尔型单元备件需求量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810313157.8A CN108520145B (zh) | 2018-04-09 | 2018-04-09 | 一种贮存失效风险下的威布尔型单元备件需求量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108520145A CN108520145A (zh) | 2018-09-11 |
CN108520145B true CN108520145B (zh) | 2022-04-05 |
Family
ID=63432288
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810313157.8A Active CN108520145B (zh) | 2018-04-09 | 2018-04-09 | 一种贮存失效风险下的威布尔型单元备件需求量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108520145B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109614583B (zh) * | 2018-10-24 | 2022-10-28 | 中国人民解放军海军工程大学 | 一种威布尔型单元备件损耗数量的计算方法 |
CN109583017B (zh) * | 2018-10-24 | 2022-10-25 | 中国人民解放军海军工程大学 | 一种正态型单元备件损耗数量的计算方法 |
CN109492914A (zh) * | 2018-11-15 | 2019-03-19 | 青岛航讯网络技术服务有限公司 | 大型货轮短期任务备件需求量确定方法 |
CN110288120B (zh) * | 2019-05-16 | 2023-07-14 | 中国人民解放军海军工程大学 | 整体贮存模式下多批次部件的备件方案优化方法及装置 |
CN110688760A (zh) * | 2019-09-30 | 2020-01-14 | 青岛航讯网络技术服务有限公司 | 表决部件的备件量计算方法、仿真方法、终端及存储介质 |
CN110727902B (zh) * | 2019-09-30 | 2023-07-18 | 中国人民解放军海军工程大学 | 一种威布尔型通用件备件需求量计算方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2557530A1 (en) * | 2011-08-11 | 2013-02-13 | General Electric Company | System and method for dynamic spare part management |
CN105844050A (zh) * | 2016-04-12 | 2016-08-10 | 吉林大学 | 基于时间相关的数控机床组件更换时间方法 |
CN105975739A (zh) * | 2016-03-02 | 2016-09-28 | 航天科工防御技术研究试验中心 | 整机产品贮存寿命建模方法 |
CN106874238A (zh) * | 2017-01-20 | 2017-06-20 | 中国人民解放军海军工程大学 | 一种威布尔型单元备件需求量的计算方法 |
CN107220216A (zh) * | 2017-05-16 | 2017-09-29 | 中国人民解放军海军工程大学 | 一种利用特征数的威布尔型备件需求量的近似计算方法 |
CN107437116A (zh) * | 2017-08-02 | 2017-12-05 | 中国航空工业集团公司西安飞机设计研究所 | 一种面向使用和维修任务的装备保障性分析方法 |
-
2018
- 2018-04-09 CN CN201810313157.8A patent/CN108520145B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2557530A1 (en) * | 2011-08-11 | 2013-02-13 | General Electric Company | System and method for dynamic spare part management |
CN105975739A (zh) * | 2016-03-02 | 2016-09-28 | 航天科工防御技术研究试验中心 | 整机产品贮存寿命建模方法 |
CN105844050A (zh) * | 2016-04-12 | 2016-08-10 | 吉林大学 | 基于时间相关的数控机床组件更换时间方法 |
CN106874238A (zh) * | 2017-01-20 | 2017-06-20 | 中国人民解放军海军工程大学 | 一种威布尔型单元备件需求量的计算方法 |
CN107220216A (zh) * | 2017-05-16 | 2017-09-29 | 中国人民解放军海军工程大学 | 一种利用特征数的威布尔型备件需求量的近似计算方法 |
CN107437116A (zh) * | 2017-08-02 | 2017-12-05 | 中国航空工业集团公司西安飞机设计研究所 | 一种面向使用和维修任务的装备保障性分析方法 |
Non-Patent Citations (2)
Title |
---|
Optimizing the damage detection cycle of the spare parts based on the failure data;Jianjun Qi ET AL.;《AIP Conference Proceedings 1864, 020201(2017)》;20170803;第020201-1至020201-6页 * |
Weibull型多部件备件需求量的改进算法;刘天华,等;《海军工程大学学报》;20120831;第41-45页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108520145A (zh) | 2018-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108520145B (zh) | 一种贮存失效风险下的威布尔型单元备件需求量计算方法 | |
CN108564270B (zh) | 一种贮存失效风险下的伽玛型单元备件需求量计算方法 | |
CN108536960B (zh) | 一种贮存失效风险下的指数型单元备件需求量计算方法 | |
CN108334720B (zh) | 一种贮存失效风险下的正态型单元备件需求量计算方法 | |
CN106874238B (zh) | 一种威布尔型单元备件需求量的计算方法 | |
CN107220216B (zh) | 一种利用特征数的威布尔型备件需求量的近似计算方法 | |
CN108304971B (zh) | 一种多台套设备的备件需求量计算方法 | |
CN108509389B (zh) | 一种威布尔型串联部件的备件需求量计算方法 | |
CN106844953B (zh) | 一种威布尔型有寿件备件的保障概率计算方法 | |
CN108304665B (zh) | 一种多台套设备的备件需求量计算方法 | |
CN108470250B (zh) | 一种正态型串联部件的备件需求量计算方法 | |
CA2501273A1 (en) | Process for determining competing cause event probability and/or system availability during the simultaneous occurrence of multiple events | |
CN109492974A (zh) | 大型货轮整体换件维修多Weibull单元部件备件需求量确定方法 | |
CN114415982B (zh) | 一种数据存储方法、装置、设备及可读存储介质 | |
CN108932358B (zh) | 一种确定威布尔型单元后续备件需求量的方法 | |
CN101571820B (zh) | 一种冗余保护方法、装置及系统 | |
CN109583017B (zh) | 一种正态型单元备件损耗数量的计算方法 | |
CN109614583B (zh) | 一种威布尔型单元备件损耗数量的计算方法 | |
CN114416437B (zh) | 一种自适应冗余级别的数据存储方法 | |
CN110399637B (zh) | 一种指数型备件寿命分布参数估计方法及系统 | |
CN109636086B (zh) | 一种伽玛型单元备件损耗数量的计算方法 | |
CN109376334B (zh) | 一种正态型单元的寿命分布参数估计方法 | |
CN110288119B (zh) | 整体贮存模式下补充贮存方案的优化方法及装置 | |
CN109145502B (zh) | 一种威布尔型单元寿命分布参数估计方法 | |
CN109388860B (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 |