CN109583017B - 一种正态型单元备件损耗数量的计算方法 - Google Patents

一种正态型单元备件损耗数量的计算方法 Download PDF

Info

Publication number
CN109583017B
CN109583017B CN201811247378.6A CN201811247378A CN109583017B CN 109583017 B CN109583017 B CN 109583017B CN 201811247378 A CN201811247378 A CN 201811247378A CN 109583017 B CN109583017 B CN 109583017B
Authority
CN
China
Prior art keywords
spare parts
spare
spare part
now
calculating
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
CN201811247378.6A
Other languages
English (en)
Other versions
CN109583017A (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 CN201811247378.6A priority Critical patent/CN109583017B/zh
Publication of CN109583017A publication Critical patent/CN109583017A/zh
Application granted granted Critical
Publication of CN109583017B publication Critical patent/CN109583017B/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

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~N(μ,σ2),其中μ为寿命的均值,σ2为寿命的方差,X的密度函数为
Figure BDA0001840832720000021
发明内容
有鉴于此,本发明提出了一种计算准确、误差小的贮存失效风险下的正态型单元备件损耗数量的计算方法。
在本文中假定:正态型单元的贮存寿命服从正态分布X~N(μ11 2),工作寿命服从正态分布X~N(μ22 2);正态型单元在工作期间为连续工作模式,当该单元出现故障时,需要备件进行换件维修,此时对备件仓库进行清点,找出可用的备件、剔除贮存失效的备件;当保障任务结束时,会对备件仓库进行清点,剔除贮存失效的备件;保障任务开始时刻同时也是备件开始贮存的时刻,且贮存寿命和工作寿命二者相互独立。
本发明的技术方案是这样实现的:本发明提供了一种正态型单元备件损耗数量的计算方法,包括以下步骤,
S1,令当前备件数量Snow=Sn,Nf=0,i=1,其中,配置的备件数量为Sn,备件损耗数量为Nf,i为清点仓库的次数;
S2.1,计算第i次清点仓库的时刻cTi
Figure BDA0001840832720000022
Tw为保障任务时间,σ2为备件工作寿命的标准差,μ2为备件工作寿命的均值;
S2.2,若i>1,则令Tzc=cTi-cTi-1;否则令Tzc=cT1,Tzc为阶段性贮存时间;
S2.3,计算贮存失效概率Pz0
Figure BDA0001840832720000031
σ1为备件贮存寿命的标准差,μ1为备件贮存寿命的均值;
S2.4,更新Nf、Snow,令Nf=Nf+Snow×Pz0,令Snow为(Snow-1)×(1-Pz0)和0之间的较大值;
S2.5,令i=i+1,若i≤Sn,则转步骤S2.1,否则转步骤S3;
S3,输出Nf为最终的备件损耗数量。
本发明的正态型单元备件损耗数量的计算方法相对于现有技术具有以下有益效果:
(1)计算准确、误差小,计算快速方便,与传统采用的存在贮存失效风险的备件保障过程仿真模拟结果较为接近,满足工程应用要求。
具体实施方式
下面将结合本发明实施方式,对本发明实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式仅仅是本发明一部分实施方式,而不是全部的实施方式。基于本发明中的实施方式,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式,都属于本发明保护的范围。
实施例1
某正态型单元的贮存寿命服从正态分布N(4000,16002),工作寿命服从正态分布N(1800,4002),保障任务时间Tw=5000h时,备件数量Sn为6,试计算备件损耗数量Nf步骤如下:
S1,初始化步骤
令当前备件数量Snow=6,Nf=0,i=1;
S2,循环计算步骤
S2.1,计算第i次清点仓库的时刻cTi
Figure BDA0001840832720000041
S2.2,计算阶段性贮存时间Tzc
因为i=1,所以令Tzc=cT1=1800.0;
S2.3,计算贮存失效概率Pz0
Figure BDA0001840832720000042
2.4,更新Nf、Snow
令Nf=Nf+Snow×Pz0=0.51,令Snow=max([(Snow-1)×(1-Pz0)0])=4.58;
2.5,令i=i+1=2,由于i≤Sn,因此转步骤S2.1,重复以上步骤;
当i=i+1=7时,由于i>Sn,因此转步骤S3;
S3,输出结果步骤
输出最终结果Nf=3.38。
实施例2
某正态型单元的贮存寿命服从正态分布N(4000,16002),工作寿命服从正态分布N(1800,4002),保障任务时间Tw=5000h时,配置的备件数量k为6,存在贮存失效风险的备件保障过程仿真流程步骤如下:
1)初始化单元的已工作时间simTw=0;
2)产生1个随机数t0,用于模拟装备中该单元的工作寿命,t0服从正态分布N(μ22 2);令simTw=t0
3)产生k个随机数t1m(1≤m≤k),用于模拟备件的贮存寿命,t1m服从正态分布N(μ11 2);
4)产生k个随机数t2m(1≤m≤k),用于模拟备件的工作寿命,t2m服从正态分布N(μ22 2);
5)比较simTw与Tw的大小情况
若simTw>Tw,则本次模拟结束,转6);
若simTw<Tw,则产生一次备件需求,清点仓库,如果不能找到贮存寿命大于simTw的备件,则本次模拟结束,转6);如果能找到一个贮存寿命大于simTw的备件(记其工作寿命为t2'),则令simTw=simTw+t2',然后把该备件从备件库中剔除,转5);
6)清点仓库,查找所有贮存寿命小于Tw的备件,这些备件为贮存失效备件,其数量为损耗数量。
按照上述流程,可多次模拟保障任务,对得到的所有备件损耗数量仿真结果进行统计,其均值即为模拟的备件损耗数量。
以下表1列出了实施例1中的计算过程中,Sn=1~10时,其对应的备件损耗数量计算结果和实施例2模拟结果。
表1备件损耗数量的模拟结果和实施例1计算结果的对比情况
Figure BDA0001840832720000051
表1的结果表明,尽管步骤S2.4会导致计算误差,本发明结果仍然与模拟结果较为接近,满足工程应用要求。
以上所述仅为本发明的较佳实施方式而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种正态型单元备件损耗数量的计算方法,其特征在于:包括以下步骤,
S1,令当前备件数量Snow=Sn,Nf=0,i=1,其中,配置的备件数量为Sn,备件损耗数量为Nf,i为清点仓库的次数;
S2.1,计算第i次清点仓库的时刻cTi
Figure FDA0001840832710000011
Tw为保障任务时间,σ2为备件工作寿命的标准差,μ2为备件工作寿命的均值;
S2.2,若i>1,则令Tzc=cTi-cTi-1;否则令Tzc=cT1,Tzc为阶段性贮存时间;
S2.3,计算贮存失效概率Pz0
Figure FDA0001840832710000012
σ1为备件贮存寿命的标准差,μ1为备件贮存寿命的均值;
S2.4,更新Nf、Snow,令Nf=Nf+Snow×Pz0,令Snow为(Snow-1)×(1-Pz0)和0之间的较大值;
S2.5,令i=i+1,若i≤Sn,则转步骤S2.1,否则转步骤S3;
S3,输出Nf为最终的备件损耗数量。
CN201811247378.6A 2018-10-24 2018-10-24 一种正态型单元备件损耗数量的计算方法 Active CN109583017B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811247378.6A CN109583017B (zh) 2018-10-24 2018-10-24 一种正态型单元备件损耗数量的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811247378.6A CN109583017B (zh) 2018-10-24 2018-10-24 一种正态型单元备件损耗数量的计算方法

Publications (2)

Publication Number Publication Date
CN109583017A CN109583017A (zh) 2019-04-05
CN109583017B true CN109583017B (zh) 2022-10-25

Family

ID=65920454

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811247378.6A Active CN109583017B (zh) 2018-10-24 2018-10-24 一种正态型单元备件损耗数量的计算方法

Country Status (1)

Country Link
CN (1) CN109583017B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116843231B (zh) * 2023-07-20 2024-03-19 中国人民解放军海军工程大学 一种考虑维修耗时的机械设备使用可用度量化方法和系统

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05290058A (ja) * 1992-04-07 1993-11-05 Nec Corp 通信システムの見積もりシステム
CN108304971B (zh) * 2018-02-07 2021-07-16 中国人民解放军海军工程大学 一种多台套设备的备件需求量计算方法
CN108470250B (zh) * 2018-03-25 2021-11-02 中国人民解放军海军工程大学 一种正态型串联部件的备件需求量计算方法
CN108334720B (zh) * 2018-04-09 2021-08-31 中国人民解放军海军工程大学 一种贮存失效风险下的正态型单元备件需求量计算方法
CN108520145B (zh) * 2018-04-09 2022-04-05 中国人民解放军海军工程大学 一种贮存失效风险下的威布尔型单元备件需求量计算方法
CN108536960B (zh) * 2018-04-09 2022-04-05 中国人民解放军海军工程大学 一种贮存失效风险下的指数型单元备件需求量计算方法

Also Published As

Publication number Publication date
CN109583017A (zh) 2019-04-05

Similar Documents

Publication Publication Date Title
CN106874238B (zh) 一种威布尔型单元备件需求量的计算方法
CN111429048B (zh) 确定补货信息的方法、装置及设备
CN106971048B (zh) 一种伽马型有寿件备件需求量计算方法
Perron et al. Searching for additive outliers in nonstationary time series
CN106777819B (zh) 一种正态型有寿件更换周期的计算方法
US20200242628A1 (en) System and method for determining a grouping of segments within a market
CN108520145B (zh) 一种贮存失效风险下的威布尔型单元备件需求量计算方法
US20090271235A1 (en) Apparatus and method for generating survival curve used to calculate failure probability
CN109583017B (zh) 一种正态型单元备件损耗数量的计算方法
CN108334720B (zh) 一种贮存失效风险下的正态型单元备件需求量计算方法
CN110598363A (zh) 表决部件的备件量计算方法、仿真方法、终端及存储介质
US20140032605A1 (en) Selection of data paths
Eilon et al. Adaptive limits in inventory control
CN107330709B (zh) 确定目标对象的方法及装置
CN109636086B (zh) 一种伽玛型单元备件损耗数量的计算方法
CN109614583B (zh) 一种威布尔型单元备件损耗数量的计算方法
CN109598464B (zh) 一种指数型单元备件损耗数量的计算方法
Sabri-Laghaie et al. Combining a maintenance center M/M/c/m queue into the economic production quantity model with stochastic machine breakdowns and repair
WO2023123982A1 (zh) 一种提高光刻胶模型精度的方法
CN107783896B (zh) 一种数据处理模型的优化方法和装置
CN106325227A (zh) 制程管控线的生成方法和装置及制程管控系统
CN110309533B (zh) 基于备件保障数据的正态型单元寿命分布参数估计方法
CN109145502B (zh) 一种威布尔型单元寿命分布参数估计方法
CN109388860B (zh) 一种伽玛型单元寿命分布参数估计方法
CN105656858B (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