CN108520145A - 一种贮存失效风险下的威布尔型单元备件需求量计算方法 - Google Patents
一种贮存失效风险下的威布尔型单元备件需求量计算方法 Download PDFInfo
- Publication number
- CN108520145A CN108520145A CN201810313157.8A CN201810313157A CN108520145A CN 108520145 A CN108520145 A CN 108520145A CN 201810313157 A CN201810313157 A CN 201810313157A CN 108520145 A CN108520145 A CN 108520145A
- Authority
- CN
- China
- Prior art keywords
- spare parts
- demand amount
- spare
- calculates
- enables
- 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
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
技术领域
本发明涉及备件需求量计算领域,尤其涉及一种贮存失效风险下的威布尔型单元备件需求量计算方法。
背景技术
准确计算备件需求量,能从备件采购费用这种经济成本的角度,定量描述装备的保障性优劣程度。
“备件在投入使用前始终是“良品”,不会贮存失效”是当前所有备件需求量计算方法的一个常见假定。如果备件贮存在环境良好的专业仓库中,则该假定是合理的,与现实中备件的实际表现也是极为相符的。但对于那些平时并非在专业仓库中贮存的备件而言,例如随装备配置的随机备件,有时受工作环境空间有限等各种条件限制,不能提供满足贮存标准的备件存放环境;尤其是对于那些长期在外场工作的装备,例如舰船上的相关装备常年处于海上,往往处于高湿、高盐度以及机电设备工作或船体摇晃带来的振动状态等较为恶劣的工作环境,如果此时的工作环境就是随机备件的贮存环境的话,则备件存在着在贮存期间失效的风险。此时,如果无视贮存失效风险,仍然采用常规方法计算备件需求量,则会因部分备件在贮存期间失效而数量不够,进而导致保障任务失败的结果。
机电件寿命一般服从威布尔型分布,如:滚珠轴承、继电器、蓄电池、液压泵、齿轮、材料疲劳件等,该分布适于描述老化导致的故障。威布尔型单元指寿命服从威布尔分布的单元,寿命X分布记作记作X~W(α,b),其中尺度参数α>0,在工程上形状参数b≥1,X的密度函数为
发明内容
有鉴于此,本发明提出了一种计算准确、误差小的贮存失效风险下的威布尔型单元备件需求量计算方法。
在本文中假定:威布尔型单元备件的贮存寿命服从威布尔分布W(α1,b1),工作寿命服从威布尔分布W(α2,b2);保障任务开始时刻同时也是备件开始贮存的时刻,且贮存寿命和工作寿命二者相互独立;当保障任务时间为Tw时,要求配置一定数量的备件,使得备件保障概率Pok不得低于保障概率指标P0。
本发明的技术方案是这样实现的:本发明提供了一种贮存失效风险下的威布尔型单元备件需求量计算方法,包括以下步骤,
S1.1,令备件保障概率其中,α2为备件工作寿命威布尔分布形状参数,b2为备件工作寿命威布尔分布尺度参数,Tw为保障任务时间;
如果Pok≥P0,则备件需求量为0,计算终止,其中P0为保障概率指标;否则,令备件数量S=1,转步骤S1.2;
S1.2,令i=1,Snow=S;
S2.1,计算故障概率gPi,
式中,
S2.2,计算故障时刻gTi,
S2.3,计算贮存效果Pzs
式中,α1为备件贮存寿命威布尔分布形状参数,b1为备件贮存寿命威布尔分布尺度参数,S1为对Snow进行取整操作后的整数,取整操作为四舍五入、向上取整或向下取整中的任意一种;
S2.4,令Pok=Pok+gPi×Pzs,令Snow=(Snow-1)×(1-Pz0);
S2.5,令i=i+1,若i≤S,则转步骤S2.1,否则转步骤S3;
S3,如果Pok≥P0,则备件需求量为S,计算终止;否则,令备件数量S=S+1,备件保障概率转步骤S1.2。
本发明的贮存失效风险下的威布尔型单元备件需求量计算方法相对于现有技术具有以下有益效果:
(1)计算准确、误差小,计算快速方便,与传统采用的一次保障任务的执行情况仿真模拟结果较为接近,满足工程应用要求。
具体实施方式
下面将结合本发明实施方式,对本发明实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式仅仅是本发明一部分实施方式,而不是全部的实施方式。基于本发明中的实施方式,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式,都属于本发明保护的范围。
实施例1
某威布尔型单元的贮存寿命服从威布尔分布W(4000,1.4),工作寿命服从威布尔分布W(3000,1.6),保障任务时间Tw=5000h时,保障概率指标P0为0.85,试计算备件需求量步骤如下:
S1,初始化步骤
S1.1,令备件保障概率
因为Pok<P0,因此令备件数量S=1,转步骤1.2。
S1.2,令i=1,Snow=S;
S2,计算保障概率Pok
S2.1,计算故障概率gP1,
令
S2.2,计算故障时刻gT1,令
S2.3,计算贮存效果Pzs
由于此时S1=1,因此
S2.4,令Pok=Pok+gP1×Pzs=0.3407,令Snow=(Snow-1)×(1-Pz0)=0;
S2.5,令i=i+1=2,由于i>S,因此转步骤3。
S3,判断步骤
因为此时Pok<P0,因此令备件数量S=S+1=2,备件保障概率转步骤1.2。
重复以上过程,最终当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.一种贮存失效风险下的威布尔型单元备件需求量计算方法,其特征在于:包括以下步骤,
S1.1,令备件保障概率其中,α2为备件工作寿命威布尔分布形状参数,b2为备件工作寿命威布尔分布尺度参数,Tw为保障任务时间;
如果Pok≥P0,则备件需求量为0,计算终止,其中P0为保障概率指标;否则,令备件数量S=1,转步骤S1.2;
S1.2,令i=1,Snow=S;
S2.1,计算故障概率gPi,
式中,
S2.2,计算故障时刻gTi,
S2.3,计算贮存效果Pzs
式中,α1为备件贮存寿命威布尔分布形状参数,b1为备件贮存寿命威布尔分布尺度参数,S1为对Snow进行取整操作后的整数,取整操作为四舍五入、向上取整或向下取整中的任意一种;
S2.4,令Pok=Pok+gPi×Pzs,令Snow=(Snow-1)×(1-Pz0);
S2.5,令i=i+1,若i≤S,则转步骤S2.1,否则转步骤S3;
S3,如果Pok≥P0,则备件需求量为S,计算终止;否则,令备件数量S=S+1,备件保障概率转步骤S1.2。
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 true CN108520145A (zh) | 2018-09-11 |
CN108520145B 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) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109492914A (zh) * | 2018-11-15 | 2019-03-19 | 青岛航讯网络技术服务有限公司 | 大型货轮短期任务备件需求量确定方法 |
CN109583017A (zh) * | 2018-10-24 | 2019-04-05 | 中国人民解放军海军工程大学 | 一种正态型单元备件损耗数量的计算方法 |
CN109614583A (zh) * | 2018-10-24 | 2019-04-12 | 中国人民解放军海军工程大学 | 一种威布尔型单元备件损耗数量的计算方法 |
CN110288120A (zh) * | 2019-05-16 | 2019-09-27 | 中国人民解放军海军工程大学 | 整体贮存模式下多批次部件的备件方案优化方法及装置 |
CN110688760A (zh) * | 2019-09-30 | 2020-01-14 | 青岛航讯网络技术服务有限公司 | 表决部件的备件量计算方法、仿真方法、终端及存储介质 |
CN110727902A (zh) * | 2019-09-30 | 2020-01-24 | 中国人民解放军海军工程大学 | 一种威布尔型通用件备件需求量计算方法及装置 |
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 |
---|
JIANJUN QI ET AL.: "Optimizing the damage detection cycle of the spare parts based on the failure data", 《AIP CONFERENCE PROCEEDINGS 1864, 020201(2017)》 * |
刘天华,等: "Weibull型多部件备件需求量的改进算法", 《海军工程大学学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109583017A (zh) * | 2018-10-24 | 2019-04-05 | 中国人民解放军海军工程大学 | 一种正态型单元备件损耗数量的计算方法 |
CN109614583A (zh) * | 2018-10-24 | 2019-04-12 | 中国人民解放军海军工程大学 | 一种威布尔型单元备件损耗数量的计算方法 |
CN109492914A (zh) * | 2018-11-15 | 2019-03-19 | 青岛航讯网络技术服务有限公司 | 大型货轮短期任务备件需求量确定方法 |
CN110288120A (zh) * | 2019-05-16 | 2019-09-27 | 中国人民解放军海军工程大学 | 整体贮存模式下多批次部件的备件方案优化方法及装置 |
CN110288120B (zh) * | 2019-05-16 | 2023-07-14 | 中国人民解放军海军工程大学 | 整体贮存模式下多批次部件的备件方案优化方法及装置 |
CN110688760A (zh) * | 2019-09-30 | 2020-01-14 | 青岛航讯网络技术服务有限公司 | 表决部件的备件量计算方法、仿真方法、终端及存储介质 |
CN110727902A (zh) * | 2019-09-30 | 2020-01-24 | 中国人民解放军海军工程大学 | 一种威布尔型通用件备件需求量计算方法及装置 |
CN110727902B (zh) * | 2019-09-30 | 2023-07-18 | 中国人民解放军海军工程大学 | 一种威布尔型通用件备件需求量计算方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108520145B (zh) | 2022-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108520145A (zh) | 一种贮存失效风险下的威布尔型单元备件需求量计算方法 | |
CN108564270A (zh) | 一种贮存失效风险下的伽玛型单元备件需求量计算方法 | |
CN108536960A (zh) | 一种贮存失效风险下的指数型单元备件需求量计算方法 | |
CN108334720A (zh) | 一种贮存失效风险下的正态型单元备件需求量计算方法 | |
CN106874238B (zh) | 一种威布尔型单元备件需求量的计算方法 | |
CN108509389A (zh) | 一种威布尔型串联部件的备件需求量计算方法 | |
CN106844953B (zh) | 一种威布尔型有寿件备件的保障概率计算方法 | |
CN108304665B (zh) | 一种多台套设备的备件需求量计算方法 | |
CN110502365A (zh) | 数据存储及恢复的方法、装置及计算机设备 | |
CN108304971B (zh) | 一种多台套设备的备件需求量计算方法 | |
CN108572793A (zh) | 数据写入和数据恢复方法、装置、电子设备及存储介质 | |
CN108470250A (zh) | 一种正态型串联部件的备件需求量计算方法 | |
CN109492974A (zh) | 大型货轮整体换件维修多Weibull单元部件备件需求量确定方法 | |
CN109523160A (zh) | 大型货轮长期保障任务多正态单元部件备件需求量确定方法 | |
CN109492914A (zh) | 大型货轮短期任务备件需求量确定方法 | |
CN108932358B (zh) | 一种确定威布尔型单元后续备件需求量的方法 | |
CN104316457B (zh) | 飞机结构防护体系日历寿命可靠度的确定方法 | |
CN106158041A (zh) | 一种列再次修复方法和装置 | |
Bowo et al. | The development of marine accidents human reliability assessment approach: heart methodology and MOP model | |
CN117132144A (zh) | 一种考虑维修耗时的机械单元使用可用度估计方法和系统 | |
CN105975739B (zh) | 整机产品贮存寿命建模方法 | |
CN109583017B (zh) | 一种正态型单元备件损耗数量的计算方法 | |
CN107703912A (zh) | 一种基于航电设备的故障诊断方法 | |
Blokus et al. | Reliability analysis of the crude oil transfer system in the oil port terminal | |
CN109614583B (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 |