中国指挥与控制学会会刊 
近年来,随着部队编制体制改革的不断深入,军队的规模结构和力量编成发生了革命性重塑,形成了军委管总、战区主战、军种主建的新格局,一体化联合作战机制体制日益凸显。部队合成化步伐不断加快,综合作战能力不断提升。与此同时,部队装备技术复杂化、种类多样化的特点也随之日益突出,对装备保障工作提出了更新更高的要求。器材保障作为实施装备保障工作的重要物质基础,逐渐成为战时装备保障研究的热点问题。
器材是指装备维修中所需元器件、零部件、附品、工具、军械油料、装护具、擦拭材料和原材料等的统称[1]。按照管理种类可分为配套器材、正常周转器材和战备储备器材[2]三种。其中正常周转器材一般以年度为基准,进行集中统一供应。对于器材需求预测问题,国内外学者做了大量研究,也建立了一些行之有效的数学模型,其大致可以分为四类:
一是基于数据统计方法,在历史统计数据的基础上,根据历史数据的变化趋势预测将来一段时间内器材的消耗量。其前提是在预测的这段时间内,影响器材消耗的主要因素不发生大的变化,且每类器材自装备服役以来的年实际消耗量和该类器材的年供应量等为已知。该方法使用范围广,主要有时间序列模型、回归分析模型、指数平滑模型、小样本模型、灰色预计模型、Bayesian模型、Croston模型、Bootstrap模型等,上述模型由于数据样本量的不同,标准计算的精度和结果有所不同[3],刘慎洋[4]基于军械装备使用及消耗分析,给出了军械器材消耗的常用模型及其消耗分布的获取方法,依据分段消耗规律给出了基于先验分布的军械器材消耗预测模型。崔亦斌[5]等针对二炮备件消耗标准制定问题,对影响备件消耗的因素进行分析,给出了单项备件的年度消耗标准以及战时备件的消耗标准模型。
二是基于寿命分布方法,该类模型建立在单元可靠性的基础上,从单元的损坏机理出发,依据单元在不同工作时间长度内发生故障的概率来预测备件在将来一段时间内的器材消耗量,备件保障度模型是其典型代表。Larry S.Mickel[6]提出寿命服从指数分布的可修备件需求量预计模型。夏长俊[7]提出寿命服从威布尔分布的可修备件需求量预计模型。
三是基于仿真预测方法,主要是在充分考虑部件消耗影响因素的基础上,做出一系列假设,运用系统仿真技术,建立装备和装备维修保障系统的仿真模型,模拟部件的使用、故障和维修过程,从而能够在最贴近实际的情况下得出部件消耗预测值。当前常用的建模方法主要有实体流图法、DEF3、Petri 网和ARIS法等。郭会军[8]等采用ARIS方法,对战时装备维修过程进行了分析,建立了战时备件需求的仿真模型。
以上方法均是从不同的角度采用不同的方法对器材进行了预测,呈现出多元化、多样化的特点,每种模型都存在一定的限定条件和使用范围,因此其预测模型的准确度和适用性也存在较大的差别。
1 问题描述
当前由于器材需求的随机性与供应保障的周期性的矛盾,导致部队有的器材积压严重,造成人力、物力和财力的极大浪费,并且器材的“保鲜”问题也随之突出。而有的器材供不应求,后续补充周期较长,造成了工作效率的大幅下降。目前已有的预测模型主要存在模型复杂、操作难度大等特点,缺乏一定的实用性。因此,寻求一种科学可靠、操作性强的方法,准确制定器材供应保障标准数量模型是当前部队维修保障亟须解决的一大难题。
2 假设与符号
2.1 基本假设
为了便于模型建立和求解,对装备器材作出如下假设。
1)装备器材适用于部队装备维修使用;
2)各单位装备器材不能调剂使用;
3)装备器材总体需求数量(分布)已知;
4)部队装备器材每年按装备器材标准补充一次。
2.2 数学符号
Xj:第j种装备器材年度需求数量,分布密度函数为fj(x)(已知);
Bj:第j种装备器材的标准数量(待求);
Mj:第j种装备器材满足数量(未知);
Qj:第j种装备器材短缺数量(未知);
Uj:第j种装备器材冗余数量(未知);
Dj:第j种装备器材误差数量(未知)。
3 模型建立
由于装备器材需求的随机性和装备器材标准的稳定性,在装备器材需求数量(实数)与装备器材标准数量(整数)之间总是存在不一致性,具体可用满足数量、短缺数量、冗余数量和误差数量来描述。
3.1 满足数量
装备器材满足数量是指按照装备器材标准进行装备器材保障,能够满足装备器材需求的数量。
第j种装备器材的满足数量为装备器材需求数量与装备器材标准数量之间的最小值。
Mj=min(Xj, Bj)
满足数量Mj的积分平均值为
E(Mj)=E[min(Xj, Bj)]= xfj(x)dx+ Bjfj(x)dx
3.2 短缺数量
装备器材短缺数量是指按照装备器材标准进行装备器材保障,不能满足装备器材需求的数量。
第j种装备器材的短缺数量为装备器材需求数量与装备器材标准数量之差。
Qj=Xj-Bj(非负)
短缺数量Qj的积分平均值为
E(Qj)=E(Xj-Bj)= (x-Bj)fj(x)dx
3.3 冗余数量
装备器材冗余数量是指按照装备器材标准进行装备器材保障,超出装备器材需求的数量。
第j种装备器材的冗余数量为装备器材标准数量与装备器材需求数量之差。
Uj=Bj-Xj(非负)
冗余数量Uj的积分平均值为
E(Uj)=E(Bj-Xj)= (Bj-x)fj(x)dx
3.4 误差数量
装备器材误差数量是指按照装备器材标准进行装备器材保障,保障的装备器材数量与需求数量之间差额的绝对值。
第j种装备器材的误差数量为装备器材需求数量与装备器材标准数量之差的绝对值。
Dj=|Bj-Xj|
误差数量的积分平均值为
E(Dj)=E(|Xj-Bj|)=E(Qj)+E(Uj)= (Bj-x)fj(x)dx+ (x-Bj)fj(x)dx
4 优化方法
4.1 满足数量优化
拟制的装备器材标准应尽量满足装备器材需求,达到装备器材满足数量最大。即
max Mj=min(Xj, Bj)
通过调整装备器材标准Bj值,使目标函数取最大值。
4.2 短缺数量优化
拟制的装备器材标准应尽量减少装备器材短缺,达到装备器材短缺数量最小。即
min Qj=Xj-Bj
通过调整装备器材标准Bj值,使目标函数取最小值。
4.3 冗余数量优化
拟制的装备器材标准应尽量减少装备器材冗余,达到装备器材冗余数量最小。即
min Uj=Bj-Xj
通过调整装备器材标准Bj值,使目标函数取最小值。
4.4 误差数量优化
拟制的装备器材标准应尽量减少装备器材保障误差,达到装备器材误差数量最小。即
min Dj=|Bj-Xj|
通过调整装备器材标准Bj值,使目标函数取最小值即可得出最优供应标准数量。因解析求解过程过于复杂,所以利用simpson公式进行数值积分,并应用Matlab编程进行求解可得,短缺数量积分均值是关于标准数量Bj的一个单调不增函数,冗余数量积分均值是关于标准数量Bj的一个单调不减函数,因此无法求出最优解,而误差数量积分均值是关于标准数量Bj的单调递减而后单调递增的一个函数,因此存在最小值,可以求出最优解。
5 算例分析
对于已知需求规律的器材,其寿命分布主要为指数分布、正态分布和威布尔分布三种类型,求解不同器材的最优供应数量都可以采用此方法,只需将分布规律进行替换即可。下面,本文以服从威布尔分布的某型器材为例,运用此模型进行分析求解。
某合成旅修理连维修某型装备所需一种器材j,上级单位每年为合成旅供应一次此类器材,且此器材各部队之间不能调剂。已知该种器材的寿命分布服从k=10、λ=1 000的威布尔分布,上级供应数量上、下限为(500,1 000)。求解上级每年供应该种器材j最优数量。
由以上数据可知,器材j服从威布尔分布:
f(x)=
a=500,c=1 000,并且Xj服从威布尔分布,其中k=10,λ=1 000,因此根据simpson公式,分别计算满足数量均值Mj、短缺数量均值Qj、冗余数量均值Uj及误差数量均值Dj可得:
E(Mj)=E[min(Xj, Bj)]= xfj(x)dx+ Bjfj(x)dx= xfj(x)dx+ Bjfj(x)dx
由simpson数值积分,并用Matlab进行求解可得图1 。
由simpson数值积分,并用Matlab进行求解可得图2 。
由simpson数值积分,并用Matlab进行求解可得图3 。
由simpson数值积分,并用Matlab进行求解可得图4 。
由图4 可知,当Bj=908时,取得最小值,误差均值最小。
由此可以得出,当上级供应器材j的数量为908时,能够最大限度满足该部队对器材j使用需求的同时,也避免了因过量储存而造成库存积压和运力浪费等问题,为部队器材供应保障筹划工作提供决策依据。
6 结束语
本文通过simpson积分公式,对器材需求与供应保障之间的误差函数进行积分和建模,并运用Matlab编程对积分均值进行求解,从而得出最优解,即供应保障最大满足标准数量,为决策者科学制定供应保障标准数量,解决部分器材积压或短缺所带来的矛盾问题,提供了解决方法,并且可以针对不同的需求规律的器材分别进行求解,最大程度保证了部队器材供应与需求的精准对接。