CN108536960A - 一种贮存失效风险下的指数型单元备件需求量计算方法 - Google Patents
一种贮存失效风险下的指数型单元备件需求量计算方法 Download PDFInfo
- Publication number
- CN108536960A CN108536960A CN201810313141.7A CN201810313141A CN108536960A CN 108536960 A CN108536960 A CN 108536960A CN 201810313141 A CN201810313141 A CN 201810313141A CN 108536960 A CN108536960 A CN 108536960A
- Authority
- CN
- China
- Prior art keywords
- spare parts
- calculates
- demand amount
- enables
- spare
- 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
Abstract
本发明提出了一种贮存失效风险下的指数型单元备件需求量计算方法,包括初始化步骤、计算保障概率和判断步骤三个步骤,计算准确、误差小,计算快速方便,与传统采用的一次保障任务的执行情况仿真模拟结果较为接近,满足工程应用要求。
Description
技术领域
本发明涉及备件需求量计算领域,尤其涉及一种贮存失效风险下的指数型单元备件需求量计算方法。
背景技术
准确计算备件需求量,能从备件采购费用这种经济成本的角度,定量描述装备的保障性优劣程度。
“备件在投入使用前始终是“良品”,不会贮存失效”是当前所有备件需求量计算方法的一个常见假定。如果备件贮存在环境良好的专业仓库中,则该假定是合理的,与现实中备件的实际表现也是极为相符的。但对于那些平时并非在专业仓库中贮存的备件而言,例如随装备配置的随机备件,有时受工作环境空间有限等各种条件限制,不能提供满足贮存标准的备件存放环境;尤其是对于那些长期在外场工作的装备,例如舰船上的相关装备常年处于海上,往往处于高湿、高盐度以及机电设备工作或船体摇晃带来的振动状态等较为恶劣的工作环境,如果此时的工作环境就是随机备件的贮存环境的话,则备件存在着在贮存期间失效的风险。此时,如果无视贮存失效风险,仍然采用常规方法计算备件需求量,则会因部分备件在贮存期间失效而数量不够,进而导致保障任务失败的结果。
一般来说,电子零部件等单元的寿命服从指数分布,如:印制电路板插件、电子部件、电阻、电容、集成电路等。指数型单元指寿命服从指数分布的单元,寿命X服从指数分布记作X~Exp(λ),X的密度函数为f(x)=λe-λx。在本文中把寿命服从指数分布的单元称之为指数型单元。
针对寿命服从指数分布的产品单元,式(1)常被用于计算备件需求量、评估备件保障效果。式中S是备件数量,Tw为保障任务时间,λ为指数分布Exp(λ)的失效率参数,用于描述该单元的寿命分布。Pok为备件保障概率,用于描述该备件保障方案的保障效果。在确定备件数量后、执行任务前,Pok是我们对某保障任务在该批备件支持下能得以成功的预期概率(任务期间的单元故障次数不大于备件数量记为成功),也被称为保障任务成功率。式(1)中没有关于贮存的相关参数,显然不能用于计算存在贮存失效风险的备件需求量。
发明内容
有鉴于此,本发明提出了一种计算准确、误差小的贮存失效风险下的指数型单元备件需求量计算方法。
在本文中假定:指数型单元的贮存寿命服从指数分布Exp(λ1),工作寿命服从指数分布Exp(λ2);保障任务开始时刻同时也是备件开始贮存的时刻,且贮存寿命和工作寿命二者相互独立;当保障任务时间为Tw时,要求配置一定数量的备件,使得备件保障概率Pok不得低于保障概率指标P0。
本发明的技术方案是这样实现的:本发明提供了一种贮存失效风险下的指数型单元备件需求量计算方法,包括以下步骤,
S1.1,令备件保障概率其中,λ2为备件工作寿命指数分布的失效率参数,Tw为保障任务时间;
如果Pok≥P0,则备件需求量为0,计算终止;否则,令备件数量S=1,转步骤S1.2;
S1.2,令i=1,Snow=S;
S2.1,计算故障概率gPi,
S2.2,计算故障时刻gTi,
S2.3,计算贮存时间Tzc,若i>1,则令Tzc=gTi-gTi-1;否则令Tzc=gT1
S2.4,计算贮存效果Pzs
式中,λ1为备件贮存寿命指数分布的失效率参数,S1为对Snow进行取整操作后的整数,取整操作为四舍五入、向上取整或向下取整中的任意一种;
S2.5,令Pok=Pok+gPi×Pzs,令Snow=(Snow-1)×(1-Pz0);
S2.6,令i=i+1,若i≤S,则转步骤S2.1,否则转步骤S3;
S3,如果Pok≥P0,则备件需求量为S,计算终止;否则,令备件数量S=S+1,备件保障概率转步骤S1.2。
本发明的贮存失效风险下的指数型单元备件需求量计算方法相对于现有技术具有以下有益效果:
(1)计算准确、误差小,计算快速方便,与传统采用的一次保障任务的执行情况仿真模拟结果较为接近,满足工程应用要求。
具体实施方式
下面将结合本发明实施方式,对本发明实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式仅仅是本发明一部分实施方式,而不是全部的实施方式。基于本发明中的实施方式,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式,都属于本发明保护的范围。
实施例1
假设某指数型单元的贮存寿命服从指数分布Exp(0.0002),工作寿命服从指数分布Exp(0.001),保障任务时间Tw=5000h时,保障概率指标P0为0.85,试计算备件需求量步骤如下:
S1,初始化步骤
S1.1,令备件保障概率Pok=e-0.001×5000=0.0067;
因为Pok<P0,因此令备件数量S=1,转步骤1.2。
S1.2,令i=1,Snow=S;
S2,计算保障概率Pok
S2.1,计算故障概率gP1,
令
S2.2,计算故障时刻gT1,令
S2.3,计算贮存时间Tzc,因为不满足i>1,因此令Tzc=gT1=993.3
S2.4,计算贮存效果Pzs,由于此时S1=1,因此Pzs=1-Pz0=e-0.0002×993.3=0.8198;
2.5,令Pok=Pok+gP1×Pzs=0.0344,令Snow=(Snow-1)×(1-Pz0)=0;
2.6,令i=i+1=2,由于i>S,因此转步骤3;
S3,判断步骤
因为此时Pok<P0,因此令备件数量S=S+1=2,备件保障概率转步骤1.2。
重复以上过程,最终当S=13时,Pok=0.8608>P0满足要求,因此,备件需求量为13。
实施例2
假设某指数型单元的贮存寿命服从指数分布Exp(0.0002),工作寿命服从指数分布Exp(0.001),保障任务时间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~13时,其对应的备件保障概率计算结果,以及实施例2中模拟结果。
表1备件保障概率的模拟结果和实施例1计算结果的对比情况
表1的结果表明,尽管步骤S2.3中的取整操作会导致计算误差,但是本发明结果仍然与模拟结果较为接近,满足工程应用要求。
以上所述仅为本发明的较佳实施方式而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种贮存失效风险下的指数型单元备件需求量计算方法,其特征在于:包括以下步骤,
S1.1,令备件保障概率其中,λ2为备件工作寿命指数分布的失效率参数,Tw为保障任务时间;
如果Pok≥P0,则备件需求量为0,计算终止;否则,令备件数量S=1,转步骤S1.2;
S1.2,令i=1,Snow=S;
S2.1,计算故障概率gPi,
S2.2,计算故障时刻gTi,
S2.3,计算贮存时间Tzc,若i>1,则令Tzc=gTi-gTi-1;否则令Tzc=gT1
S2.4,计算贮存效果Pzs
式中,λ1为备件贮存寿命指数分布的失效率参数,S1为对Snow进行取整操作后的整数,取整操作为四舍五入、向上取整或向下取整中的任意一种;
S2.5,令Pok=Pok+gPi×Pzs,令Snow=(Snow-1)×(1-Pz0);
S2.6,令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 |
---|---|---|---|
CN201810313141.7A CN108536960B (zh) | 2018-04-09 | 2018-04-09 | 一种贮存失效风险下的指数型单元备件需求量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810313141.7A CN108536960B (zh) | 2018-04-09 | 2018-04-09 | 一种贮存失效风险下的指数型单元备件需求量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108536960A true CN108536960A (zh) | 2018-09-14 |
CN108536960B CN108536960B (zh) | 2022-04-05 |
Family
ID=63479564
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810313141.7A Active CN108536960B (zh) | 2018-04-09 | 2018-04-09 | 一种贮存失效风险下的指数型单元备件需求量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108536960B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109583017A (zh) * | 2018-10-24 | 2019-04-05 | 中国人民解放军海军工程大学 | 一种正态型单元备件损耗数量的计算方法 |
CN109598464A (zh) * | 2018-10-24 | 2019-04-09 | 中国人民解放军海军工程大学 | 一种指数型单元备件损耗数量的计算方法 |
CN109614583A (zh) * | 2018-10-24 | 2019-04-12 | 中国人民解放军海军工程大学 | 一种威布尔型单元备件损耗数量的计算方法 |
CN110287523A (zh) * | 2019-05-16 | 2019-09-27 | 中国人民解放军海军工程大学 | 模块化贮存模式下多批次部件的备件方案优化方法及装置 |
CN110288120A (zh) * | 2019-05-16 | 2019-09-27 | 中国人民解放军海军工程大学 | 整体贮存模式下多批次部件的备件方案优化方法及装置 |
CN110738007A (zh) * | 2019-09-30 | 2020-01-31 | 中国人民解放军海军工程大学 | 一种电子类通用件的备件保障概率计算方法及装置 |
CN116843119A (zh) * | 2023-05-23 | 2023-10-03 | 中国人民解放军海军工程大学 | 一种考虑维修耗时的电子单元备件设计方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103426073A (zh) * | 2013-08-21 | 2013-12-04 | 易程科技股份有限公司 | 一种基于设备寿命预测的客服系统备件库存管理方法 |
CN104462755A (zh) * | 2014-10-30 | 2015-03-25 | 中国船舶重工集团公司第七二六研究所 | 基于可靠性模型的电子设备备件配置计算方法 |
WO2015175570A1 (en) * | 2014-05-12 | 2015-11-19 | Virginia Tech Intellectual Properties, Inc. | Selective modulation of intracellular effects of cells using pulsed electric fields |
CN106845109A (zh) * | 2017-01-20 | 2017-06-13 | 中国人民解放军海军工程大学 | 一种指数型有寿件备件需求量的计算方法 |
CN106971048A (zh) * | 2017-04-14 | 2017-07-21 | 中国人民解放军海军工程大学 | 一种伽马型有寿件备件需求量计算方法 |
CN107885926A (zh) * | 2017-11-06 | 2018-04-06 | 中国人民解放军海军工程大学 | 一种无维护状态下舰艇装备的备件配置优化方法 |
-
2018
- 2018-04-09 CN CN201810313141.7A patent/CN108536960B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103426073A (zh) * | 2013-08-21 | 2013-12-04 | 易程科技股份有限公司 | 一种基于设备寿命预测的客服系统备件库存管理方法 |
WO2015175570A1 (en) * | 2014-05-12 | 2015-11-19 | Virginia Tech Intellectual Properties, Inc. | Selective modulation of intracellular effects of cells using pulsed electric fields |
CN104462755A (zh) * | 2014-10-30 | 2015-03-25 | 中国船舶重工集团公司第七二六研究所 | 基于可靠性模型的电子设备备件配置计算方法 |
CN106845109A (zh) * | 2017-01-20 | 2017-06-13 | 中国人民解放军海军工程大学 | 一种指数型有寿件备件需求量的计算方法 |
CN106971048A (zh) * | 2017-04-14 | 2017-07-21 | 中国人民解放军海军工程大学 | 一种伽马型有寿件备件需求量计算方法 |
CN107885926A (zh) * | 2017-11-06 | 2018-04-06 | 中国人民解放军海军工程大学 | 一种无维护状态下舰艇装备的备件配置优化方法 |
Non-Patent Citations (2)
Title |
---|
EWELINA KOSICKA ET AL.: "Intelligent Systems of Forecasting the Failure of Machinery Park and Supporting Fulfilment of Orders of Spare Parts", 《INTELLIGENT SYSTEMS IN PRODUCTION ENGINEERING AND MAINTENANCE – ISPEM 2017》 * |
刘任洋,等: "串件拼修策略下指数型不完全修复件的可用度评估", 《系统工程理论与实践》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109583017A (zh) * | 2018-10-24 | 2019-04-05 | 中国人民解放军海军工程大学 | 一种正态型单元备件损耗数量的计算方法 |
CN109598464A (zh) * | 2018-10-24 | 2019-04-09 | 中国人民解放军海军工程大学 | 一种指数型单元备件损耗数量的计算方法 |
CN109614583A (zh) * | 2018-10-24 | 2019-04-12 | 中国人民解放军海军工程大学 | 一种威布尔型单元备件损耗数量的计算方法 |
CN110287523A (zh) * | 2019-05-16 | 2019-09-27 | 中国人民解放军海军工程大学 | 模块化贮存模式下多批次部件的备件方案优化方法及装置 |
CN110288120A (zh) * | 2019-05-16 | 2019-09-27 | 中国人民解放军海军工程大学 | 整体贮存模式下多批次部件的备件方案优化方法及装置 |
CN110288120B (zh) * | 2019-05-16 | 2023-07-14 | 中国人民解放军海军工程大学 | 整体贮存模式下多批次部件的备件方案优化方法及装置 |
CN110738007A (zh) * | 2019-09-30 | 2020-01-31 | 中国人民解放军海军工程大学 | 一种电子类通用件的备件保障概率计算方法及装置 |
CN110738007B (zh) * | 2019-09-30 | 2022-10-28 | 中国人民解放军海军工程大学 | 一种电子类通用件的备件保障概率计算方法及装置 |
CN116843119A (zh) * | 2023-05-23 | 2023-10-03 | 中国人民解放军海军工程大学 | 一种考虑维修耗时的电子单元备件设计方法及系统 |
CN116843119B (zh) * | 2023-05-23 | 2024-02-20 | 中国人民解放军海军工程大学 | 一种考虑维修耗时的电子单元备件设计方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108536960B (zh) | 2022-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108536960A (zh) | 一种贮存失效风险下的指数型单元备件需求量计算方法 | |
CN108564270A (zh) | 一种贮存失效风险下的伽玛型单元备件需求量计算方法 | |
CN108509389B (zh) | 一种威布尔型串联部件的备件需求量计算方法 | |
Jin et al. | Efficient board-level functional fault diagnosis with missing syndromes | |
CN108304665B (zh) | 一种多台套设备的备件需求量计算方法 | |
CN105868543B (zh) | 基于逆高斯寿命分布的贮存寿命试验加速因子评估方法 | |
CN108304971A (zh) | 一种多台套设备的备件需求量计算方法 | |
CN108470250A (zh) | 一种正态型串联部件的备件需求量计算方法 | |
CN104468691B (zh) | 用于飞行器的控制系统和飞行器 | |
CN108572793A (zh) | 数据写入和数据恢复方法、装置、电子设备及存储介质 | |
CN105512457B (zh) | 用于提供关于维护活动的信息的计算系统和方法 | |
CN108984381A (zh) | 数据异常分析的方法及系统 | |
Garland et al. | Modern electronic maintenance principles | |
Youssef et al. | Availability assessment of multi-state manufacturing systems using universal generating function | |
US20190333294A1 (en) | Detecting fault states of an aircraft | |
CN115907476A (zh) | 一种制氢装置hazop半定量分析方法和装置 | |
Vintr | Reliability assessment for components of complex mechanisms and machines | |
CN109598464A (zh) | 一种指数型单元备件损耗数量的计算方法 | |
JP6302493B2 (ja) | 回路基板の修正作業を支援する装置 | |
US7873883B2 (en) | Method for scrubbing storage in a computer memory | |
Ablanedo-Rosas et al. | Benchmarking the operational efficiency of major container ports in five developing regions | |
CN108038317A (zh) | 精密仪器性能参数保持期预测方法和系统 | |
Yaghoubi et al. | A new method for computation the success probability of coverage for switch unit in the switching systems | |
CN109636086B (zh) | 一种伽玛型单元备件损耗数量的计算方法 | |
CN109188331A (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 |