CN108334720B - 一种贮存失效风险下的正态型单元备件需求量计算方法 - Google Patents

一种贮存失效风险下的正态型单元备件需求量计算方法 Download PDF

Info

Publication number
CN108334720B
CN108334720B CN201810313142.1A CN201810313142A CN108334720B CN 108334720 B CN108334720 B CN 108334720B CN 201810313142 A CN201810313142 A CN 201810313142A CN 108334720 B CN108334720 B CN 108334720B
Authority
CN
China
Prior art keywords
spare parts
spare
calculating
storage
demand
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
CN201810313142.1A
Other languages
English (en)
Other versions
CN108334720A (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 CN201810313142.1A priority Critical patent/CN108334720B/zh
Publication of CN108334720A publication Critical patent/CN108334720A/zh
Application granted granted Critical
Publication of CN108334720B publication Critical patent/CN108334720B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/08—Logistics, e.g. warehousing, loading or distribution; Inventory or stock management
    • G06Q10/087—Inventory or stock management, e.g. order filling, procurement or balancing against orders

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Finance (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Geometry (AREA)
  • Accounting & Taxation (AREA)
  • Evolutionary Computation (AREA)
  • Development Economics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提出了一种贮存失效风险下的正态型单元备件需求量计算方法，包括初始化步骤、计算保障概率和判断步骤三个步骤，计算准确、误差小，计算快速方便，与传统采用的一次保障任务的执行情况仿真模拟结果较为接近，满足工程应用要求。

Description

一种贮存失效风险下的正态型单元备件需求量计算方法
技术领域
本发明涉及备件需求量计算领域，尤其涉及一种贮存失效风险下的正态型单元备件需求量计算方法。
背景技术
准确计算备件需求量，能从备件采购费用这种经济成本的角度，定量描述装备的保障性优劣程度。
“备件在投入使用前始终是“良品”，不会贮存失效”是当前所有备件需求量计算方法的一个常见假定。如果备件贮存在环境良好的专业仓库中，则该假定是合理的，与现实中备件的实际表现也是极为相符的。但对于那些平时并非在专业仓库中贮存的备件而言，例如随装备配置的随机备件，有时受工作环境空间有限等各种条件限制，不能提供满足贮存标准的备件存放环境；尤其是对于那些长期在外场工作的装备，例如舰船上的相关装备常年处于海上，往往处于高湿、高盐度以及机电设备工作或船体摇晃带来的振动状态等较为恶劣的工作环境，如果此时的工作环境就是随机备件的贮存环境的话，则备件存在着在贮存期间失效的风险。此时，如果无视贮存失效风险，仍然采用常规方法计算备件需求量，则会因部分备件在贮存期间失效而数量不够，进而导致保障任务失败的结果。
机械件寿命一般服从正态分布，如汇流环、齿轮箱、减速器等。正态型单元指寿命服从正态分布的单元，寿命X分布记作X～N(μ,σ2)，其中μ为寿命的均值，σ2为寿命的方差，X的密度函数为
Figure BDA0001623040720000011
发明内容
有鉴于此，本发明提出了一种计算准确、误差小的贮存失效风险下的正态型单元备件需求量计算方法。
在本文中假定：正态型单元的贮存寿命服从正态分布N(μ1,σ1 2)，工作寿命服从正态分布N(μ2,σ2 2)；保障任务开始时刻同时也是备件开始贮存的时刻，且贮存寿命和工作寿命二者相互独立；当保障任务时间为Tw时，要求配置一定数量的备件，使得备件保障概率Pok不得低于保障概率指标P0。
本发明的技术方案是这样实现的：本发明提供了一种贮存失效风险下的正态型单元备件需求量计算方法，包括以下步骤，
S1.1，令备件保障概率
Figure BDA0001623040720000021
其中，σ2为备件工作寿命标准差，μ2为备件工作寿命的均值，Tw为保障任务时间；
如果Pok≥P0，则备件需求量为0，计算终止；否则，令备件数量S＝1，转步骤S1.2；
S1.2，令i＝1，Snow＝S；
S2.1，计算故障概率gPi，
Figure BDA0001623040720000022
S2.2，计算故障时刻gTi，
Figure BDA0001623040720000023
S2.3，计算贮存效果Pzs
Figure BDA0001623040720000024
式中，
Figure BDA0001623040720000025
σ1为备件贮存寿命标准差，μ1为备件贮存寿命的均值，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，备件保障概率
Figure BDA0001623040720000031
转步骤S1.2。
本发明的贮存失效风险下的正态型单元备件需求量计算方法相对于现有技术具有以下有益效果：
(1)计算准确、误差小，计算快速方便，与传统采用的一次保障任务的执行情况仿真模拟结果较为接近，满足工程应用要求。
具体实施方式
下面将结合本发明实施方式，对本发明实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式仅仅是本发明一部分实施方式，而不是全部的实施方式。基于本发明中的实施方式，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式，都属于本发明保护的范围。
实施例1
假设某正态型单元的贮存寿命服从正态分布N(4000,6502)，工作寿命服从正态分布N(2000,5502)，保障任务时间Tw＝5000h时，保障概率指标P0为0.85，，试计算备件需求量步骤如下：
S1，初始化步骤
S1.1，令备件保障概率
Figure BDA0001623040720000032
因为Pok＜P0，因此令备件数量S＝1，转步骤1.2。
S1.2，令i＝1，Snow＝S；
S2，计算保障概率Pok
S2.1，计算故障概率gP1，
令
Figure BDA0001623040720000041
S2.2，计算故障时刻gT1，令
Figure BDA0001623040720000042
S2.3，计算贮存效果Pzs
由于此时S1＝1，因此
Figure BDA0001623040720000043
2.4，令Pok＝Pok+gP1×Pzs＝0.0992，令Snow＝(Snow-1)×(1-Pz0)＝0；
2.5，令i＝i+1＝2，由于i＞S，因此转步骤3；
S3，判断步骤
因为此时Pok＜P0，因此令备件数量S＝S+1＝2，备件保障概率
Figure BDA0001623040720000044
转步骤1.2。
重复以上过程，最终当S＝6时，Pok＝0.8532＞P0满足要求，因此，备件需求量为6。
实施例2
假设某正态型单元的贮存寿命服从正态分布N(4000,6502)，工作寿命服从正态分布N(2000,5502)，保障任务时间Tw＝5000h时，保障概率指标P0为0.85，当备件数量为k时，一次保障任务的执行情况仿真流程步骤如下：
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计算结果的对比情况
Figure BDA0001623040720000051
表1的结果表明，尽管步骤S2.3中的取整操作会导致计算误差，但是本发明结果仍然与模拟结果较为接近，满足工程应用要求。
以上所述仅为本发明的较佳实施方式而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.一种贮存失效风险下的正态型单元备件需求量计算方法，其特征在于：包括以下步骤，
S1.1，令备件保障概率
Figure FDA0001623040710000011
其中，σ2为备件工作寿命标准差，μ2为备件工作寿命的均值，Tw为保障任务时间；
如果Pok≥P0，则备件需求量为0，计算终止；否则，令备件数量S＝1，转步骤S1.2；
S1.2，令i＝1，Snow＝S；
S2.1，计算故障概率gPi，
Figure FDA0001623040710000012
S2.2，计算故障时刻gTi，
Figure FDA0001623040710000013
S2.3，计算贮存效果Pzs
Figure FDA0001623040710000014
式中，
Figure FDA0001623040710000015
σ1为备件贮存寿命标准差，μ1为备件贮存寿命的均值，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，备件保障概率
Figure FDA0001623040710000021
转步骤S1.2。
CN201810313142.1A 2018-04-09 2018-04-09 一种贮存失效风险下的正态型单元备件需求量计算方法 Active CN108334720B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810313142.1A CN108334720B (zh) 2018-04-09 2018-04-09 一种贮存失效风险下的正态型单元备件需求量计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810313142.1A CN108334720B (zh) 2018-04-09 2018-04-09 一种贮存失效风险下的正态型单元备件需求量计算方法

Publications (2)

Publication Number Publication Date
CN108334720A CN108334720A (zh) 2018-07-27
CN108334720B true CN108334720B (zh) 2021-08-31

Family

ID=62933998

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810313142.1A Active CN108334720B (zh) 2018-04-09 2018-04-09 一种贮存失效风险下的正态型单元备件需求量计算方法

Country Status (1)

Country Link
CN (1) CN108334720B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109583017B (zh) * 2018-10-24 2022-10-25 中国人民解放军海军工程大学 一种正态型单元备件损耗数量的计算方法
CN109508792A (zh) * 2018-10-31 2019-03-22 中航航空服务保障(天津)有限公司 一种飞机定检耗材供应清单的确定方法
CN110287523B (zh) * 2019-05-16 2023-07-18 中国人民解放军海军工程大学 模块化贮存模式下多批次部件的备件方案优化方法及装置
CN116955910B (zh) * 2023-07-19 2024-04-02 中国人民解放军海军工程大学 一种机械串联部件的保障任务成功率计算方法和系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8342279B1 (en) * 2009-09-21 2013-01-01 The Boeing Company Modular vehicle and associated method of construction
CN103426073A (zh) * 2013-08-21 2013-12-04 易程科技股份有限公司 一种基于设备寿命预测的客服系统备件库存管理方法
CN106777819A (zh) * 2017-01-20 2017-05-31 中国人民解放军海军工程大学 一种正态型有寿件更换周期的计算方法
CN106844953A (zh) * 2017-01-20 2017-06-13 中国人民解放军海军工程大学 一种威布尔型有寿件备件的保障概率计算方法
CN107220216A (zh) * 2017-05-16 2017-09-29 中国人民解放军海军工程大学 一种利用特征数的威布尔型备件需求量的近似计算方法
CN107436983A (zh) * 2017-07-28 2017-12-05 南京理工大学 一种基于多元样本差异的o型橡胶密封圈寿命预测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8342279B1 (en) * 2009-09-21 2013-01-01 The Boeing Company Modular vehicle and associated method of construction
CN103426073A (zh) * 2013-08-21 2013-12-04 易程科技股份有限公司 一种基于设备寿命预测的客服系统备件库存管理方法
CN106777819A (zh) * 2017-01-20 2017-05-31 中国人民解放军海军工程大学 一种正态型有寿件更换周期的计算方法
CN106844953A (zh) * 2017-01-20 2017-06-13 中国人民解放军海军工程大学 一种威布尔型有寿件备件的保障概率计算方法
CN107220216A (zh) * 2017-05-16 2017-09-29 中国人民解放军海军工程大学 一种利用特征数的威布尔型备件需求量的近似计算方法
CN107436983A (zh) * 2017-07-28 2017-12-05 南京理工大学 一种基于多元样本差异的o型橡胶密封圈寿命预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Optimal Inventory Control of Spare Parts Based on Storage Availability;ZHOU Liang ET AL.;《2016 28th Chinese Control and Decision Conference (CCDC)》;20161231;第6409-6414页 *
正态型有寿件的备件方案确定方法;邵松世，等;《航空学报》;20180125;第221147-1至221147-6页 *

Also Published As

Publication number Publication date
CN108334720A (zh) 2018-07-27

Similar Documents

Publication Publication Date Title
CN108564270B (zh) 一种贮存失效风险下的伽玛型单元备件需求量计算方法
CN108334720B (zh) 一种贮存失效风险下的正态型单元备件需求量计算方法
CN108536960B (zh) 一种贮存失效风险下的指数型单元备件需求量计算方法
CN108509390B (zh) 一种伽玛型串联部件的备件需求量计算方法
CN107220216B (zh) 一种利用特征数的威布尔型备件需求量的近似计算方法
CN106844953B (zh) 一种威布尔型有寿件备件的保障概率计算方法
CN108304971B (zh) 一种多台套设备的备件需求量计算方法
CN108304665B (zh) 一种多台套设备的备件需求量计算方法
CN108470250B (zh) 一种正态型串联部件的备件需求量计算方法
CN106201338A (zh) 数据存储方法及装置
KR101572038B1 (ko) 메모리 어레이의 에러 수정 방법
CN101739306A (zh) 数据错误处理方法、数据错误检查和纠正装置及系统
CN111930538A (zh) 一种基于kafka集群的生产与消费的方法
CN116469448B (zh) 一种闪存颗粒的筛选方法和装置
CN101571820A (zh) 一种冗余保护方法、装置及系统
CN109583017B (zh) 一种正态型单元备件损耗数量的计算方法
CN114416437B (zh) 一种自适应冗余级别的数据存储方法
CN109376334B (zh) 一种正态型单元的寿命分布参数估计方法
CN109003643A (zh) 一种数据处理方法及装置
Allen Model specification in regression analysis
CN110399637B (zh) 一种指数型备件寿命分布参数估计方法及系统
CN109271275A (zh) 一种固态硬盘中的坏块的标记方法、装置、设备及介质
CN109636086B (zh) 一种伽玛型单元备件损耗数量的计算方法
CN112906951A (zh) 一种多海工模块组合装船方案的筛选方法
Thas et al. The embedding of (0, α)-geometries in PG (n, q): part II

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