一种确定伽玛型单元后续备件需求量的方法技术

一种确定伽玛型单元后续备件需求量的方法，属于设备保障性能设计方法，解决现有伽玛型单元后续备件需求量方法所存在的不能准确计算“(旧品)单元+(新品)备件”情况下的后续备件需求量问题。本发明专利技术包括设置初始值步骤、计算剩余寿命失效度步骤、计算保障概率步骤和判断步骤。本发明专利技术考虑了已工作一段时间的“旧品单元”对后续备件需求量的影响，能准确计算后续备件需求量及其备件数量对应的保障概率，为实现装备的精确保障奠定基础。

【技术实现步骤摘要】
一种确定伽玛型单元后续备件需求量的方法
本专利技术属于设备保障性能设计方法，特别涉及一种确定伽玛型单元后续备件需求量的方法。
技术介绍
“单元和备件都是新品”是绝大多数备件需求量计算方法的前提假定，此时计算所得的备件需求量结果称之为初始备件方案，所谓“初始”即单元处于初始零时刻、刚开始使用的意思。如果把已工作一段时间的单元视为新品，用于计算后续备件需求量的话，则会导致备件数量偏少、保障不足的风险；如果把已工作一段时间的单元索性视为故障件，用于计算后续备件需求量的话，则会导致备件数量偏多、过度保障的情况。实际上，在大多数情况下，某单元已经工作了一段时间且还未发生故障，仍能在下一任务期间继续使用，可称为“旧品单元”。因此，如何计算“旧品单元+新品备件”情况下的后续备件需求量，相比初始备件的计算问题，是一个更为普遍的问题。理论上，只有寿命服从指数分布的单元，其后续备件可以采用初始备件需求量计算方法。伽玛分布是一种常见的分布类型，适用于描述工程实际中性能逐步连续退化的过程，例如刀具的磨损是一个典型的连续时间、连续状态的性能退化过程，其寿命可用伽玛分布来表示。对于寿命服从伽玛分布的单元，本专利技术简称伽玛型单元。现有确定伽玛型单元后续备件需求量的方法中，一种常见的思路是将已工作一段时间的单元仍然视为新品，在“新品单元+新品备件”假定下，用计算初始备件需求量的方法来确定后续备件需求量，这会导致备件数量偏少、保障不足的风险；另一种常见的思路是索性把已工作一段时间的单元视为故障件，在“故障单元+新品备件”假定下计算后续备件需求量，这会导致备件数量偏多、过度保障的情况。记伽玛型单元的寿命z服从伽玛分布Ga(α,λ)，其中α＞0为形状参数，λ＞0为尺度参数，针对具体的伽玛型单元，这两个参数是确定的，伽玛分布的密度函数f(z)：式中Γ(α)为伽玛函数，且假定单元已经正常工作了时间t1且未发生故障，剩余寿命记为T1，即单元寿命为t1+T1，剩余寿命指单元已经工作了一段时间且未发生故障条件下，还能再次正常工作的时间。
技术实现思路
本专利技术提供一种确定伽玛型单元后续备件需求量的方法，解决现有伽玛型单元后续备件需求量方法所存在的不能准确计算“(旧品)单元+(新品)备件”情况下的后续备件需求量问题。本专利技术所提供的一种确定伽玛型单元后续备件需求量的方法，包括设置初始值步骤、计算剩余寿命失效度步骤、计算保障概率步骤和判断步骤，其特征在于：(1)设置初始值步骤：设置保障概率目标值P0，0<P0<1，置备件需求变量j＝0；(2)计算剩余寿命失效度步骤：计算剩余寿命失效度F(x|t1)：式中，x为时间变量，且x>0，t1为单元已经正常工作时间，λ＞0为尺度参数，α＞0为形状参数，Γ(α)为伽玛函数，且y为自变量；(3)计算保障概率步骤：计算备件数量为j时的保障概率P：上式中，F(x|t1)为步骤(2)中的剩余寿命失效度，Tw为单元下一任务的计划工作时间，Tw＞0；(4)判断步骤：判断是否P≥P0，是则得到的j值即为备件需求量，否则将j+1的值赋予j，转步骤(3)。步骤(3)中，剩余寿命失效度F(x|t1)体现了“旧品单元”已工作时间t1对后续备件需求量的影响，且步骤(3)把备件需求量计算中原本复杂的j+1重积分问题简化成一个二重积分问题。后续备件的数量较多，都采用保障过度的计算会造成浪费，本专利技术能够排除保障过度和保障不足的风险，利用旧品的剩余寿命和新品备件的寿命之和的计算实际是执行的卷积运算。本专利技术考虑了已工作一段时间的“旧品单元”对后续备件需求量的影响，能准确计算后续备件需求量及其备件数量对应的保障概率，为实现装备的精确保障奠定基础。具体实施方式以下结合实施例对本专利技术进一步说明。实施例：某单元的寿命服从伽玛分布Ga(2.1,0.002)，该单元已累计正常工作时间t1为650h，预计下一任务该单元的计划工作时间Tw为1500h，要求备件保障概率≥0.85，计算为该任务准备的后续备件需求量，包括以下步骤：(1)设置初始值步骤：设置保障概率目标值P0＝0.85，置备件需求变量j＝0；(2)计算剩余寿命失效度步骤：计算剩余寿命失效度F(x|650)：(3)计算保障概率步骤：计算备件数量为j时的保障概率P：上式中，Tw为单元下一任务的计划工作时间；(4)判断步骤：判断是否P≥P0，是则得到的j值即为备件需求量，否则将j+1的值赋予j，转步骤(3)。步骤(3)计算结果如表1所示：表1计算结果表1中的模拟结果是采用仿真方法模拟得到的。从表1可知，备件需求量为2时，保障概率能满足不低于0.85的保障概率要求，且由步骤(3)计算的保障概率与模拟结果极为吻合。当备件数量为j时，模拟一次保障过程的仿真方法如下：(1)产生1个随机数simT1，用于模拟已正常工作t1的单元寿命，simT1服从伽玛分布Ga(α,λ)，且simT1＞t1；(2)判断是否j＞0，是则进行步骤(3)，否则转步骤(4)；(3)产生j个随机数simTk(1≤k≤j)，用于模拟j个备件的寿命，simTk服从伽玛分布Ga(α,λ)，令转步骤(5)；(4)如果j＝0，则令simT＝simT1-t1；转步骤(5)；(5)判断是否simT≥Tw，是则令flag＝1，本次备件保障任务成功，否则令flag＝0，本次备件保障任务失败。模拟多次上述保障过程，得到大量的flag值，多次flag的统计均值为模拟的备件数量为j时对应的保障概率。本文档来自技高网...

【技术保护点】
1.一种确定伽玛型单元后续备件需求量的方法，包括设置初始值步骤、计算剩余寿命失效度步骤、计算保障概率步骤和判断步骤，其特征在于：(1)设置初始值步骤：设置保障概率目标值P0，0

【技术特征摘要】
1.一种确定伽玛型单元后续备件需求量的方法，包括设置初始值步骤、计算剩余寿命失效度步骤、计算保障概率步骤和判断步骤，其特征在于：(1)设置初始值步骤：设置保障概率目标值P0，0<P0<1，置备件需求变量j＝0；(2)计算剩余寿命失效度步骤：计算剩余寿命失效度F(x|t1)：式中，x为时间变量，且x>0，t1为单元已经正...

【专利技术属性】
技术研发人员：邵松世，谷高全，刘凯，李华，
申请(专利权)人：中国人民解放军海军工程大学，
类型：发明
国别省市：湖北,42