CN108932358A - 一种确定威布尔型单元后续备件需求量的方法 - Google Patents

一种确定威布尔型单元后续备件需求量的方法 Download PDF

Info

Publication number
CN108932358A
CN108932358A CN201810339513.3A CN201810339513A CN108932358A CN 108932358 A CN108932358 A CN 108932358A CN 201810339513 A CN201810339513 A CN 201810339513A CN 108932358 A CN108932358 A CN 108932358A
Authority
CN
China
Prior art keywords
demand
unit
spares
replenished
spare part
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
Application number
CN201810339513.3A
Other languages
English (en)
Other versions
CN108932358B (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 CN201810339513.3A priority Critical patent/CN108932358B/zh
Publication of CN108932358A publication Critical patent/CN108932358A/zh
Application granted granted Critical
Publication of CN108932358B publication Critical patent/CN108932358B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04Ageing analysis or optimisation against ageing

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)
  • Debugging And Monitoring (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种确定威布尔型单元后续备件需求量的方法,本发明属于设备保障性能设计方法,解决现有威布尔型单元后续备件需求量方法所存在的不能准确计算“(旧品)单元+(新品)备件”情况下的后续备件需求量问题。后续备件的数量较多,都采用过度保障的计算会造成浪费,本发明能够排除保障不足或者过度保障的风险,利用旧品的剩余寿命和新品备件的寿命之和的计算实际是执行的卷积运算,本发明考虑了已工作一段时间的“旧品单元”对后续备件需求量的影响,能准确计算后续备件需求量及其备件数量对应的保障概率,为实现装备的精确保障奠定基础。

Description

一种确定威布尔型单元后续备件需求量的方法
技术领域
本发明属于设备保障性能设计方法,特别涉及一种确定威布尔型单元后续备件需求量的方法。
背景技术
“单元和备件都是新品”是绝大多数备件需求量计算方法的前提假定,此时计算所得的备件需求量结果称之为初始备件方案,所谓“初始”即单元处于初始零时刻、刚开始工作的意思。实际上,在大多数情况下,该单元已经工作了一段时间且还未发生故障,仍能在下一任务期间继续使用,可称为“旧品单元”。因此,如何计算“旧品单元+新品备件”情况下的后续备件需求量,相比初始备件的计算问题,是一个更为普遍的问题。理论上,只有寿命服从指数分布的单元,其后续备件可以采用初始备件需求量计算方法。现有确定威布尔型单元后续备件需求量的方法中,一种常见的思路是将已工作一段时间的单元仍然视为新品,在“新品单元+新品备件”假定下,用计算初始备件需求量的方法来确定后续备件需求量,这会导致备件数量偏少、保障不足的风险;另一种常见的思路是索性把已工作一段时间的单元视为故障件,在“故障单元+新品备件”假定下计算后续备件需求量,这会导致备件数量偏多、过度保障的情况。
威布尔分布用来描述那些失效率随时间变化的产品单元,解释因老化、磨损而导致的故障统计规律,主要适用于机电类,是一种应用范围广泛的分布类型。对于寿命服从威布尔分布的单元,例如:滚珠轴承、继电器、开关、断路器、某些电容器、电子管、磁控管、电位计、陀螺、电动机、航空发电机、蓄电池、液压泵、空气涡轮发动机、齿轮、活门、材料疲劳件等,本文简称威布尔型单元,本发明针对威布尔型单元给出了一种后续备件需求量计算方法。
记威布尔型单元的寿命Z服从威布尔分布W(α,b),其中α>0为尺度参数,b>0为形状参数,针对具体的威布尔型单元,这两个参数是确定的,威布尔分布的概率密度函数f(Z):
假定单元已经正常工作了时间t1且未发生故障,剩余寿命记为T1,即单元寿命为t1+T1,剩余寿命指单元已经工作了一段时间且未发生故障条件下,还能再次正常工作的时间。
发明内容
本发明提供一种确定威布尔型单元后续备件需求量的方法,解决现有威布尔型单元后续备件需求量方法所存在的不能准确计算“(旧品)单元+(新品)备件”情况下的后续备件需求量问题。
本发明所提供的一种确定威布尔型单元后续备件需求量的方法,包括以下步骤,其特征在于:
(1)设置初始值步骤:
设置保障概率目标值P0,0<P0<1,置备件需求变量j=0;
(2)计算剩余寿命失效度步骤:
计算剩余寿命失效度F(x|t1):
式中,x为时间变量,且x>0,t1为单元已经正常工作时间,α>0为尺度参数,b>0为形状参数;
(3)计算保障概率步骤:
计算备件数量为j时的保障概率P:
式中,Γ(·)为伽玛函数,且x为时间变量;F(x|t1)为步骤(2)中的剩余寿命失效度,Tw为单元下一任务的计划工作时间,Tw>0;
(4)判断步骤:
判断是否P≥P0,是则得到的j值即为备件需求量,否则将j+1的值赋予j,转步骤(3)。
步骤(3)中,剩余寿命失效度F(x|t1)体现了“旧品单元”已工作时间t1对后续备件需求量的影响,且步骤(3)把备件需求量计算中原本复杂的j+1重积分问题简化成一个二重积分问题。
后续备件的数量较多,都采用过度保障的计算会造成浪费,本发明能够排除保障不足或者过度保障的风险,利用旧品的剩余寿命和新品备件的寿命之和的计算实际是执行的卷积运算,本发明考虑了已工作一段时间的“旧品单元”对后续备件需求量的影响,能准确计算后续备件需求量及其备件数量对应的保障概率,为实现装备的精确保障奠定基础。
具体实施方式
以下结合实施例对本发明进一步说明。
实施例1:某单元的寿命服从威布尔分布W(1000,1.9),该单元已累计正常工作时间t1为650h,预计下一任务该单元的计划工作时间Tw为2000h,要求备件保障概率≥0.85,计算为该任务准备的后续备件需求量,包括以下步骤:
(1)设置初始值步骤:
设置保障概率目标值P0=0.85,置备件需求变量j=0;
(2)计算剩余寿命失效度步骤:
计算剩余寿命失效度F(x|t1):
(3)计算保障概率步骤:
计算备件数量为j时的保障概率P:
式中,Tw=2000为单元下一任务的计划工作时间;
(4)判断步骤:
判断是否P≥P0,是则得到的j值即为备件需求量,否则将j+1的值赋予j,转步骤(3)。
步骤(3)计算结果如表1所示:
表1实施例1计算结果
表1中的模拟结果是采用仿真方法模拟得到的。
从表1可知,备件需求量为3时,保障概率能满足不低于0.85的保障概率要求,且由步骤(3)计算的保障概率与模拟结果极为吻合。
算例2:某单元的寿命服从威布尔分布W(1000,3.9),该单元已累计正常工作时间t1为650h,预计下一任务该单元的计划工作时间Tw为2000h,要求备件保障概率不≥0.85,计算为该任务准备的后续备件需求量,包括以下步骤:
(1)设置初始值步骤:
设置保障概率目标值P0=0.85,置备件需求变量j=0;
(2)计算剩余寿命失效度步骤:
计算剩余寿命失效度F(x|t1):
(3)计算保障概率步骤:
计算备件数量为j时的保障概率P:
式中,Tw=2000为单元下一任务的计划工作时间;
(4)判断步骤:
判断是否P≥P0,是则得到的j值即为备件需求量,否则将j+1的值赋予j,转步骤(3)。
步骤(3)计算结果如表2所示:
表2实施例2计算结果
表2中的模拟结果是采用仿真方法模拟得到的。
从表2可知,备件需求量为3时,保障概率能满足不低于0.85的保障概率要求,且由步骤(3)计算的保障概率与模拟结果极为吻合。
当备件数量为j时,模拟一次保障过程的仿真方法如下:
(1)产生1个随机数simT1,用于模拟已正常工作t1的单元寿命,simT1服从威布尔分布W(α,b),且simT1>t1
(2)判断是否j>0,是则进行步骤(3),否则转步骤(4);
(3)产生j个随机数simTk(1≤k≤j),用于模拟j个备件的寿命,simTk服从威布尔分布W(α,b),令转步骤(5);
(4)如果j=0,则令simT=simT1-t1;转步骤(5);
(5)判断是否simT≥Tw,是则令flag=1,本次备件保障任务成功,否则令flag=0,本次备件保障任务失败。

Claims (1)

1.一种确定威布尔型单元后续备件需求量的方法,包括以下步骤,其特征在于:
(1)设置初始值步骤:
设置保障概率目标值P0,0<P0<1,置备件需求变量j=0;
(2)计算剩余寿命失效度步骤:
计算剩余寿命失效度F(x|t1):
式中,x为时间变量,且x>0,t1为单元已经正常工作时间,α>0为尺度参数,b>0为形状参数;
(3)计算保障概率步骤:
计算备件数量为j时的保障概率P:
式中,Γ(·)为伽玛函数,且x为时间变量;F(x|t1)为步骤(2)中的剩余寿命失效度,Tw为单元下一任务的计划工作时间,Tw>0;
(4)判断步骤:
判断是否P≥P0,是则得到的j值即为备件需求量,否则将j+1的值赋予j,转步骤(3)。
CN201810339513.3A 2018-04-16 2018-04-16 一种确定威布尔型单元后续备件需求量的方法 Active CN108932358B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810339513.3A CN108932358B (zh) 2018-04-16 2018-04-16 一种确定威布尔型单元后续备件需求量的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810339513.3A CN108932358B (zh) 2018-04-16 2018-04-16 一种确定威布尔型单元后续备件需求量的方法

Publications (2)

Publication Number Publication Date
CN108932358A true CN108932358A (zh) 2018-12-04
CN108932358B CN108932358B (zh) 2021-08-13

Family

ID=64448088

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810339513.3A Active CN108932358B (zh) 2018-04-16 2018-04-16 一种确定威布尔型单元后续备件需求量的方法

Country Status (1)

Country Link
CN (1) CN108932358B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110717267A (zh) * 2019-09-30 2020-01-21 中国人民解放军海军工程大学 一种伽玛型通用件备件需求量计算方法及装置
CN110727902A (zh) * 2019-09-30 2020-01-24 中国人民解放军海军工程大学 一种威布尔型通用件备件需求量计算方法及装置
CN110738008A (zh) * 2019-09-30 2020-01-31 中国人民解放军海军工程大学 一种机电类通用件的备件保障概率计算方法及装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101169801A (zh) * 2007-08-18 2008-04-30 燕山大学 备件寿命分布在事后维修中对系统可靠性影响的仿真方法
US20080154458A1 (en) * 2006-12-21 2008-06-26 Brandstetter Jason C System, method and program product for predicting fleet reliability and maintaining a fleet of vehicles
CN105302098A (zh) * 2015-11-11 2016-02-03 同济大学 一种基于ietm的地铁车辆互操作维修保障平台及其搭建方法
CN105607628A (zh) * 2016-03-18 2016-05-25 中国直升机设计研究所 一种直升机使用维修保障一体化系统
CN106844953A (zh) * 2017-01-20 2017-06-13 中国人民解放军海军工程大学 一种威布尔型有寿件备件的保障概率计算方法
CN106874238A (zh) * 2017-01-20 2017-06-20 中国人民解放军海军工程大学 一种威布尔型单元备件需求量的计算方法
CN107220216A (zh) * 2017-05-16 2017-09-29 中国人民解放军海军工程大学 一种利用特征数的威布尔型备件需求量的近似计算方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080154458A1 (en) * 2006-12-21 2008-06-26 Brandstetter Jason C System, method and program product for predicting fleet reliability and maintaining a fleet of vehicles
CN101169801A (zh) * 2007-08-18 2008-04-30 燕山大学 备件寿命分布在事后维修中对系统可靠性影响的仿真方法
CN105302098A (zh) * 2015-11-11 2016-02-03 同济大学 一种基于ietm的地铁车辆互操作维修保障平台及其搭建方法
CN105607628A (zh) * 2016-03-18 2016-05-25 中国直升机设计研究所 一种直升机使用维修保障一体化系统
CN106844953A (zh) * 2017-01-20 2017-06-13 中国人民解放军海军工程大学 一种威布尔型有寿件备件的保障概率计算方法
CN106874238A (zh) * 2017-01-20 2017-06-20 中国人民解放军海军工程大学 一种威布尔型单元备件需求量的计算方法
CN107220216A (zh) * 2017-05-16 2017-09-29 中国人民解放军海军工程大学 一种利用特征数的威布尔型备件需求量的近似计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
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)》 *
刘任洋,等: "任意寿命分布单元表决系统备件需求量的解析算法", 《系统工程与电子技术》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110717267A (zh) * 2019-09-30 2020-01-21 中国人民解放军海军工程大学 一种伽玛型通用件备件需求量计算方法及装置
CN110727902A (zh) * 2019-09-30 2020-01-24 中国人民解放军海军工程大学 一种威布尔型通用件备件需求量计算方法及装置
CN110738008A (zh) * 2019-09-30 2020-01-31 中国人民解放军海军工程大学 一种机电类通用件的备件保障概率计算方法及装置
CN110717267B (zh) * 2019-09-30 2022-08-05 中国人民解放军海军工程大学 一种伽玛型通用件备件需求量计算方法及装置
CN110727902B (zh) * 2019-09-30 2023-07-18 中国人民解放军海军工程大学 一种威布尔型通用件备件需求量计算方法及装置

Also Published As

Publication number Publication date
CN108932358B (zh) 2021-08-13

Similar Documents

Publication Publication Date Title
CN108509389A (zh) 一种威布尔型串联部件的备件需求量计算方法
CN108932358A (zh) 一种确定威布尔型单元后续备件需求量的方法
EP3163706B1 (en) Voltage stability monitoring device and method
CN108733883A (zh) 一种确定正态型单元后续备件需求量的方法
CN108509390A (zh) 一种伽玛型串联部件的备件需求量计算方法
CN109787243B (zh) 一种人工紧急减负荷在线优化决策方法、系统及存储介质
CN108470250A (zh) 一种正态型串联部件的备件需求量计算方法
CN104779608A (zh) 基于离线策略的电力系统安全稳定在线附加紧急控制方法
CN109492974A (zh) 大型货轮整体换件维修多Weibull单元部件备件需求量确定方法
US11038352B1 (en) Control system with load selection strategy that minimizes load overshed
US20170338653A1 (en) Demand control apparatus, group generation apparatus, and group generation method
CN103310307A (zh) 基于概率风险评估的电力系统规划方案灵活性评估方法
CN109523160A (zh) 大型货轮长期保障任务多正态单元部件备件需求量确定方法
JP2007189840A (ja) 電力系統安定化装置
CN111736076B (zh) 电池系统状态判断方法、装置、可读存储介质和电子设备
CN109672170B (zh) 一种区域电网有功备用评估方法和系统
CN117132144A (zh) 一种考虑维修耗时的机械单元使用可用度估计方法和系统
CN109583779B (zh) 一种电-气互联系统的薄弱元件确定方法
CN105044518B (zh) 应用于电压监测中的自动适应门限确定方法
CN104795816B (zh) 一种基于负荷重要性的切负荷方法
CN109376380A (zh) 一种确定伽玛型单元后续备件需求量的方法
CN110119995A (zh) 考虑调峰和断面双重约束的新能源消纳方法及系统
CN110197320B (zh) 一种含多端柔性多状态开关的配电网过载风险评估方法
CN110970899B (zh) 一种多区域紧急减负荷协同决策方法、系统及存储介质
CN114330865A (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