CN109492914A - 大型货轮短期任务备件需求量确定方法 - Google Patents

Abstract

用于大型货轮短期任务备件需求量确定方法，属于机械设备保障性的设计方法，解决机械设备保障性设计过程中备件需求量计算困难、不准确的问题。本发明以可用度为保障要求，包括设置初值步骤、计算失效率步骤、计算保障概率步骤和判断步骤。本发明能针对Weibull单元的可靠性特性、立足多单元k/N表决部件的真实可靠性结构，进行快速计算备件数量，解决现有机械设备保障性设计中备件需求量计算不准确的问题。在机械设备的保障性设计阶段和列装运行阶段，能准确确定备件需求量进而能较为准确地估计全寿命周期机械设备相关费用，为评估机械设备费效比状态、实现机械设备“好保障、保障好”提供技术支撑。

Description

大型货轮短期任务备件需求量确定方法

技术领域

本发明属于机械设备保障性的设计方法，特别涉及换件维修模式下，可用为保障要求的大型货轮短期任务备件需求量确定方法。

背景技术

备件需求量是机械设备保障性中的一种具体属性，备件所对应的费用是全寿命期间机械设备费用的重要组成部分。准确计算备件需求量不论是对机械设备的保障性设计还是机械设备有效运行都具有重要的意义。短期次保障任务在本文中指机械设备工作的时间与单元平均寿命的比例在0～3之间，本文中该比例不大于3。

在机械设备的层次结构中：机械设备由多种部件组成，部件由N个单元构成(N≥1)；单元与单元之间以串联或某种可靠性结构形式进行连接，构成部件，该部件称之为单元装机数为N的部件，简称多单元部件。Weibull单元指该单元的寿命服从Weibull分布，单元的寿命记为W(a1,b1),a1、b1分别为W(a1,b1)的尺度参数和形状参数，其失效密度函数对单元或部件进行可靠性寿命试验，利用寿命试验数据进行Weibull分布拟合计算，得到a1、b1。Weibull寿命件主要适用于机电件，如：滚珠轴承、断路器、电位计、发动机、蓄电池、液压泵、材料疲劳件等。多Weibull单元部件指由多个Weibull单元组成的部件。

k/N表决部件是由N个单元组成的多单元部件，其中至少有k个单元正常工作，部件才正常工作，1≤k≤N，记为k/N(G)，它的可靠性框图如图1所示；N/N(G)为串联部件，1/N(G)为并联部件。与串联、并联结构相比，表决结构是一种更为一般、普遍的可靠性结构形式。表决结构常见于两种情况：一种在机械设备的可靠性设计中，针对单元可靠度达不到要求的情况，采用表决结构作为一种冗余设计来提高部件的整体可靠度；另一种情况，则是在部件出现劣化或部分故障导致性能下降后，通过定义k来决定所能接受的部件最低性能，小于k则认为部件整体性能已不能满足工作要求。

备件需求量与维修模式紧密相关。换件维修主要针对不可修单元，或由于条件限制只能在现场进行换件维修的情况----大型货轮在海上运输期间，由于受维修条件、维修能力等限制，对出现故障的单元，船员只能采取换件维修的方法排除故障、恢复机械设备性能。换件维修模式下，单元装机数为N的表决部件，现有确定备件需求量的方法有：

(1)卷积法。在数学理论上换件维修等同于卷积计算。备件需求量为m的计算过程意味着一个(m+1)重卷积计算。对多重卷积而言，此时卷积表达式形式极其复杂、难以推导；并且，利用计算机进行多重卷积的数值求解计算，也随着备件需求量m每增加1个,其计算量以10倍以上方式递增，不仅耗时长且难以保证计算精度。

(2)模拟方法。该方法以多次模拟表决部件在保障期间工作运行的方式，来统计计算保障概率要求下所需的备件需求数量。与计算备件需求量相比，该方法更适合用于备件数量确定后的保障效果评估。这是由于随机因素的影响，使得备件数量对应的保障概率同时具有数学期望值和标准差两种属性，模拟方法得到的保障概率往往在保障概率期望值附近波动，这种波动在某些情况下(例如较长任务时间)会导致备件数量出现一个较大的波动，而导致模拟方法往往不能得到稳定的备件需求量结果。

(3)近似法。由一个Weibull单元构成的部件，其备件需求量S的近似计算公式：

式中：

E---平均寿命，a1为尺度参数，b1为形状参数。

k---变异系数，

u_p---正态分布分位数，可从GB/T4086.1从查出。

T为保障任务时间。该公式不能用于由多单元组成部件的情况。

发明内容

本发明提供的一种大型货轮短期任务备件需求量确定方法，以可用度为保障要求，解决机械设备保障性设计过程中备件需求量计算困难、不准确的问题。

本发明所提供的一种大型货轮短期任务备件需求量确定方法，以可用度为保障要求，包括设置初值步骤、计算失效率步骤、计算保障概率步骤和判断步骤，其特征在于：

(1)设置初值步骤：

设置保障概率目标值P0＝0.5～0.99，置备件需求变量j＝0；

(2)计算失效率步骤：

将Weibull单元近似为指数单元，计算指数单元的失效率λ：

式中，a0、b0分别为Weibull单元寿命W(a0,b0)的尺度参数和形状参数，由Weibull单元可靠性寿命试验数据进行Weibull分布拟合计算得到；T为保障任务时间，Gamma函数

(3)计算保障概率步骤：

计算保障概率P：

其中，R(x)为等效Gamma部件在x时刻的可靠度，a＝N-k+1，N为单元装机数，k为保证部件正常工作的最少单元数量，均由部件特性确定，1≤k≤N；T为保障任务时间，i为部件正常工作允许的失效单元数量，最大为N－k；j为备件需求变量，x、y为时间变量；

(4)判断步骤：

判断是否P≥P0，是则得到备件需求量j×(N-k+1)；否则j＝j+1，转步骤(3)。

所述的大型货轮短期任务备件需求量确定方法，其特征在于：

所述计算失效率步骤中，利用Matlab中的gamma函数计算指数单元失效率λ：

所述的大型货轮短期任务备件需求量确定方法，其特征在于：

所述计算保障概率步骤中，先利用Matlab中的gamcdf函数计算R(x)，R(x)＝1-gamcdf(x,a(j+1),b)；再利用Matlab中的quad积分函数来计算保障概率P，P＝quad(@R(x),0,T)/T。

本发明将多Weibull单元表决部件等效成Gamma部件，其寿命记为G(a,b)，服从Gamma分布，G(a,b)的失效密度函数以可用度为保障概率指标，再根据失效密度函数计算保障概率P。

按照冲击失效理论：寿命分布服从G(a,b)的产品，能经受若干次外界冲击，但当产品受冲击次数累积到a次冲击作用在时就会产生失效。根据k/N表决部件定义：当失效单元数量累积到(N-k+1)时，该部件会失效。因此，G(a,b)产品和k/N表决部件二者在失效条件上具有等效的物理意义。本发明步骤(2)将多Gamma单元k/N表决部件近似为多指数单元k/N表决部件，本发明步骤(3)进一步将近似的多指数单元表决部件等效为Gamma部件，该部件寿命G(a,b)服从Gamma分布，其失效密度函数其中，Gamma函数保障概率指标可用度P即为等效Gamma部件的平均可靠度，为等效Gamma部件失效密度函数的二重积分：从而将以往部件备件数量的多重卷积计算过程转换为等效Gamma部件失效密度函数的二重积分计算，避免了多重卷积数值计算过程中迭代计算导致的误差大、耗时长的问题。

本发明能针对针对大型货轮短期任务Weibull单元的可靠性特性、立足多单元k/N表决部件的真实可靠性结构，进行快速计算备件数量，解决现有机械设备保障性设计中备件需求量计算不准确的问题。在机械设备的保障性设计阶段和列装运行阶段，能确定备件需求量进而能较为准确地估计单次保障任务机械设备运行相关费用，为评估机械设备费效比状态、实现机械设备“好保障、保障好”提供技术支撑。

附图说明

图1是表决结构示意图；

图2是本发明的流程框图；

图3是模拟方法的流程框图。

具体实施方式

以下结合附图和实施例对本发明进一步说明。

实施例1，某12/15部件由15个蓄电池单元组成，当有效单元数量小于12时，认定该部件出现故障，对失效单元进行换件维修。蓄电池单元服从Weibull分布W(1000,2)，单元平均寿命为886h。在可用度P不低于0.9的要求下，计算保障时间500h所需的备件数量：

本实施例包括设置初值步骤、计算失效率步骤、计算保障概率步骤和判断步骤：

(1)设置初值步骤：

设置保障概率目标值P0＝0.9，置备件需求变量j＝0；

(2)计算失效率步骤：

将Gamma单元近似为指数单元，计算指数单元失效率λ：

式中，a0＝1000、b0＝2，保障任务时间T＝500，Gamma函数

利用Matlab中gamma函数计算指数单元失效率λ：

λ＝5.0000×10^-4；

(3)计算保障概率步骤：

计算保障概率P：

其中，a＝4，N＝15为单元装机数，k＝12为保证部件正常工作的最少单元数量，均由部件特性确定；i为部件正常工作允许的失效单元数量，最大为3；j为备件需求变量，x、y为时间变量；

先利用Matlab中的gamcdf函数计算R(x)，R(x)＝1-gamcdf(x,4(j+1),149.2)；再利用Matlab中的quad积分函数来计算保障概率P，P＝quad(@R(x),0,500)/500。

计算结果如下表：

j	0	1
			保障概率P	0.864	0.997

(4)判断步骤：

判断是否P≥P0，是则得到备件需求量j×4；否则j＝j+1，转步骤(3)。

当j＝1时，满足保障要求，则备件需求量为4。

当保障任务时间为1500h，蓄电池备件数量为0～20，分别以本发明和模拟方法计算其对应的可用度保障概率，结果表明本发明的方法计算结果准确。本发明的计算耗时不到模拟方法的1/10。计算结果见下表。

备件数量	0	4	8	12	16	20
							本发明：可用度	0.224	0.449	0.666	0.848	0.953	0.991
模拟方法：可用度	0.357	0.541	0.708	0.877	0.974	0.998

实施例2，某3/5部件由5个电位计单元组成，当有效单元数量小于3时，认定该部件出现故障，对失效单元进行换件维修。电位计单元服从Weibull分布W(1500,2.5)，单元平均寿命为1331h。在可用度P不低于0.9的要求下，计算保障时间1500h所需的备件数量：

本实施例包括设置初值步骤、计算失效率步骤、计算保障概率步骤和判断步骤：

(1)设置初值步骤：

设置保障概率目标值P0＝0.9，置备件需求变量j＝0；

(2)计算失效率步骤：

将Gamma单元近似为指数单元，计算指数单元失效率λ：

式中，a0＝1500、b0＝2.5，保障任务时间T＝1500，Gamma函数

利用Matlab中gamma函数计算指数单元失效率λ：

λ＝5.5717×10^-4；

(3)计算保障概率步骤：

计算保障概率P：

其中，a＝3，N＝5为单元装机数，k＝3为保证部件正常工作的最少单元数量，均由部件特性确定；i为部件正常工作允许的失效单元数量，最大为2；j为备件需求变量，x、y为时间变量；

先利用Matlab中的gamcdf函数计算R(x)，R(x)＝1-gamcdf(x,3(j+1),468.6)；再利用Matlab中的quad积分函数来计算保障概率P，P＝quad(@R(x),0,1500)/1500。

计算结果如下表：

j	0	1
			保障概率P	0.752	0.978

(4)判断步骤：

判断是否P≥P0，是则得到备件需求量j×3；否则j＝j+1，转步骤(3)。

当j＝1时，满足保障要求，则备件需求量为3。

保障任务时间为2500h，电位计备件数量为0～9，分别以本发明和模拟方法计算其对应的可用度保障概率，结果表明本发明的方法计算结果准确。本发明的计算耗时不到模拟方法的1/10。计算结果见下表。

备件数量	0	3	6	9
					本发明：可用度	0.537	0.882	0.985	0.999
模拟方法：可用度	0.524	0.872	0.996	1.000

上述实施例所涉及的模拟方法，系采用离散事件系统仿真技术来模拟多单元部件在工作过程中，随着仿真时间推进，某单元随机出现故障引发部件失效，进而实施换件维修、部件恢复功能继续工作，直至任务执行完毕或因备件消耗完毕无法排除故障造成部件停机，达到仿真终止条件。如该次仿真以任务执行完毕终止，则认为本次保障成功，部件工作时间为任务时间；如仿真因部件故障停机而终止，则认为本次保障失败，部件工作时间为故障停机时间(一般假定换件维修耗时为0)。对于多单元k/N表决部件，采用所述模拟方法进行仿真的流程如图3所示：

(1)仿真初始化：设置保障任务时间T、备件数量，置仿真时间gT＝0；

(2)根据单元的寿命分布特性随机生成部件中各单元寿命T_i；

(3)仿真时间向前推进一个步长；

(4)判断是否到达任务时间终点，是则转(8)，否则转(5)；

(5)更新各单元寿命：令各单元寿命值减1；

(6)判断部件中，寿命值大于0的单元数量是否小于k，是则部件出现故障，转(7)，否则部件正常工作，转(3)；

(7)判断该部件的单元是否还有备件，是则进行换件维修，对更换后的单元随机产生单元寿命，并更新备件数量(如只更换失效单元，则库存单元备件数量减(N-k+1)；如是部件整体更换维修，库存单元备件数量减N)，转(3)；

(8)仿真终止。

成功率和可用度是机械设备使用方最为关心的两种常用保障概率形式。在多次仿真后，保障成功次数与仿真次数之间的比例即为保障成功率，部件工作时间平均值与任务时间之间的比例即为可用度。

Claims (3)

1.一种大型货轮短期任务备件需求量确定方法，以可用度为保障要求，包括设置初值步骤、计算失效率步骤、计算保障概率步骤和判断步骤，其特征在于：

(1)设置初值步骤：

设置保障概率目标值P0＝0.5～0.99，置备件需求变量j＝0；

(2)计算失效率步骤：

将Weibull单元近似为指数单元，计算指数单元的失效率λ：

式中，a0、b0分别为Weibull单元寿命W(a0,b0)的尺度参数和形状参数，由Weibull单元可靠性寿命试验数据进行Weibull分布拟合计算得到；T为保障任务时间，Gamma函数

(3)计算保障概率步骤：

计算保障概率P：

其中，R(x)为等效Gamma部件在x时刻的可靠度，a＝N-k+1，N为单元装机数，k为保证部件正常工作的最少单元数量，均由部件特性确定，1≤k≤N；T为保障任务时间，i为部件正常工作允许的失效单元数量，最大为N－k；j为备件需求变量，x、y为时间变量；

(4)判断步骤：

判断是否P≥P0，是则得到备件需求量j×(N-k+1)；否则j＝j+1，转步骤(3)。

2.如权利要求1所述的大型货轮短期任务备件需求量确定方法，其特征在于：

所述计算失效率步骤中，利用Matlab中的gamma函数计算指数单元失效率λ：

3.如权利要求1所述的大型货轮短期任务备件需求量确定方法，其特征在于：

所述计算保障概率步骤中，先利用Matlab中的gamcdf函数计算R(x)，R(x)＝1-gamcdf(x,a(j+1),b)；再利用Matlab中的quad积分函数来计算保障概率P，P＝quad(@R(x),0,T)/T。