【装备制造业发展研究】
用户需求的不确定性,使得其对于装备制造企业的产品需求也不相同。这种不确定性包括产品品种不确定、数量不确定、质量不确定、交货期不确定等。由于客户订单的不确定性以及各种情况的发生,虽然客户需求的产品品种不同,需求的时间、数量、质量以及交货期等都不相同,但是这些订单需求存在一定的相似性,如果能找到它们之间的相似性,就可以把相似性大的订单进行并单处理或者组合处理。这样就能有效地利用有限的资源满足各种复杂的订单需求,同时也可节省大量时间成本和生产成本。因此,研究不确定需求下订单相似性十分必要,可以对出现相似性的订单需求进行相似性智能处理,从而方便企业制定生产计划或者修改、取消计划,并制定出相应的订单排产策略。
针对制造企业订单相似性,很多学者和专家都进行了相应的工作和研究。韩文民等(2018)研究了紧急订单下序位相似性问题[1]。张于贤等(2018)研究了多品种小批量的产品交货期和加工工艺相似性订单成组问题[2]。吴天行等(2017)基于巷道相似性和通道相似性研究了订单分批问题[3-4]。陈俊等(2017)建立了产品质量属性相似度模型和优化方法[5]。程贤福等(2017)研究了产品相似性等方面的适应性指标[6]。徐新胜等(2016)提出了解决产品结构相似性的判断方法[7]。胡小林等(2019)研究了需求不均衡下拣选截止时间和通道的相似性[8]。禹春霞等(2019)为了满足用户需求,研究了用户与候选服务商的相似性[9]。李若岩等(2019)采用模糊数学聚类法分析青海石脑油相似性[10]。施伟锋等(2019)采用模糊C均值聚类算法分析属性空间相似性[11]。王永泉等(2019)采用模糊聚类多层次分析法研究了零件间加工工艺相似性[12]。张青山等(2018)研究了双渠道采购模式下装备制造业选址生产问题[13]。韩世莲(2016)研究了可以满足多种客户需求属性的模糊系统聚类方法[14]。岳文军(2015)采用RSCF分析法分析了订单的相似性等特征对物流体系的影响[15]。
虽然很多学者对于产品或订单相似性进行研究,但是没有从需求的角度分析,研究在多重不确定需求的情况下订单相似性问题。虽然有研究订单的交货期、工艺路线和数量,但是没有将各种情况进行综合研究。本文针对装备制造企业多品种生产中基于多重不确定性需求的产品订单相似性进行分析。
假设装备制造企业共有N种产品,某段时间内有H个产品订单需求,UncertainDemandh表示某一时刻的需求,包括需求到达的时间、需求的品种、需求的数量、需求的质量和需求的交货日期等,具体可表示为UncertainDemandh=(DNh,DVh,DQh,DWh,DTh,PSh)。式中:h=1,2,…,H;DNh表示第h订单需求到达的时间,其取值为0,1,2,…的非负整数,假设为随机离散变量;DVh表示为第h订单需求的产品品种,其取值为1,2,…,N,N+1的正整数,假设为随机离散变量;DQh表示第h订单需求的需求数量,其取值为正整数或者为批量的倍数,假设为随机离散变量;DWh表示第h订单需求的质量需求,假设为模糊区间需求;DTh表示订单交货期需求,假设交货期需求为模糊区间需求;PSh表示订单需求的情景,取值为:PSh=1表示正常订单,PSh=2表示紧急订单,PSh=3表示修改订单。不确定需求的构成如表1所示。
表1 不确定需求的构成
不确定需求DNhDVhDQhDWhDThPShUncertainDemand1DN1DV1DQ1DW1DT1PS1UncertainDemand2DN2DV2DQ2DW2DT2PS2UncertainDemand3DN3DV3DQ3DW3DT3PS3…………………UncertainDemandHDNHDVHDQHDWHDTHPSH
对于产品的相似性.主要是考虑产品加工工艺的相似性。对于装备制造企业来说,可以利用加工工艺的工序步骤来区分不同品种产品的相似性。对产品的细化,可以根据组件或者部件相似性来定义。产品相似性主要是由部件组成,包括组成部件的种类、部件的相似性、组件个数和产品组装工艺。而产品组件的相似性由部件或组件的加工工艺、加工生产车间、加工时间等来定义。
设装备制造企业的工件或组件有M个,Componentm表示第m个工件或组件。假设装备制造企业共有Z个加工工序,则Processz表示第z个加工工序,并设ProcessNoz=z,z=1,2,…,Z,具体如表2所示。
表2 工序内容信息
序号12…Z工件或组件Process1Process2…Processz
设第m个工件或组件需要加工OPNumberm个工序,用OPm(g)表示第m个组件的第g个加工工序,则
OPm(g)=ProcessNoz (z=1,2,…,Z)
(1)
另外,设第m个工件或组件Componentm的加工工艺路线为Routingm,设其包括RoutingNumberm个工序,则
Routingm=OPm(1)→OPm(2)→ … →
OPm(OPNumberm)
(2)
则企业的工件或组件的内容信息构成如表3所示。
表3 工件或组件的内容信息
序号工件或组件工序数量工件或组件的加工工艺路线1Component1RoutingNumber1Routing1=OP1(1)→OP1(2)→…→OP1(OPNumber1)2Component2RoutingNumber2Routing2=OP2(1)→OP2(2)→…→OP2(OPNumber2)…………MComponentMRoutingNumberMRoutingM=OPM(1)→OPM(2)→…→OPM(OPNumberM)
根据文献[16]的工序相似性算法计算组件的相似性距离,记为Distance(m1,m2),假设m1≠m2,表示工件或组件m2到工件或组件m1的加工工艺相似性距离。
计算工件或组件加工的步骤如下:
(1) 以m1组件为基准,对于其每一个工序i=1,2,…,RoutingNumberm1计算其在m2组件工序中的序号,用PosNumi(m1,m2)表示m1组件第i个工序OPm1(i)在m2组件工序Routingm2中的序号。如果m1组件的OPm1(i)在m2组件工序中则得到其序号,否则记为0。
(2) 对于步骤(1)中序号PosNumi(OPm1(i),Routingm)之间的距离则有
(3)
其中:i=1,2,…,RoutingNumberm1-1
(3) 组件的相似性距离
RoutingNumberm1+1
(4)
则构建组件间相似性距离M×M方形矩阵,组件间相似性距离如表4所示。
由表4可知,以组件中距离最远的组件与该行组件最不相似,距离最近的组件与该行组件最为相似。
表4 组件间相似性距离
组件Component1Component2…ComponentMComponent10Distance(2,1)…Distance(M,1)Component2Distance(1,2)0…Distance(M,2)……………ComponentMDistance(1,M)Distance(2,M)…0
假设similiar(i,j)表示组件组件i到j的相似性,则
(k=1,2,…,M)
(5)
式中,dt是产品相似性偏移量。
本文采用取平均值最小的行所在的组件作为基准组件,以该组件为对比组件计算组件的相似性。假设取的平均值最小的行号为minLine,则组件相似性为similiar(minLine,j),其中j=1,2,…,M。
首先,由于表1中多重不确定性需求的各项数据单位并不一致,所以需要进行归一化处理,可以采用如下处理:
DNh取最低值为1,其他值为与该最低值的比值;DVh取最低值为1,其他值为与该最低值的比值;DQh取最低值为1,其他值为与该最低值的比值;DWh采用三角模糊进行处理,然后取模糊数最低值为1,其他值为与该最低值的比值;DTh采用梯形模糊进行处理,然后取模糊数最低值为1,其他值为与该最低值的比值;PSh以正常订单为1,其他值为与该最低值的比值。
对UncertainDemandh进行归一化处理。采用组件相似性为similiar(minLine,j)对产品种类DVh进行归一化处理。到达时间DNh和交货期DTh取最小的时间为1,其他值为与该最低值的比值。设归一化后的矩阵UncertainDemand′={UncertainDemand′1,UncertainDemand′2,…,Uncertain-Demand′h},其中UncertainDemand′h={DN′h,DV′h,DQ′h,DW′h,PS′h}表示第h=1,2,…,H个需求。
对于UncertainDemand′具有H个订单需求,假设将这H个订单需求分成type_num(2≤type_num≤H)个类别,则UncertainDemand′被分成{SK1,SK2,…,SKtype_num}的type_num个类别,U为相似分类值,并设各类别的聚类中心为{center1,center2,…,centertype_num},设μhk为订单需求Uncertain-Demand′h对于SKk的隶属度,则可以将目标函数Jb表示为
(6)
dhk=d(UncertainDemand′h-centerk)=
(7)
式中:dhk表示第h个订单需求UncertainDemand′h与SKk的欧几里得距离;m为订单需求的特征数;b为加权参数,其取值为1≤b≤∞。
为寻找一个最佳的分类,模糊C-均值聚类方法(FCM)要求该订单需求对于类别的隶属度值和为1,从而产生最小的函数值Jb,即满足
(8)
利用式(5)、(6)来分别计算订单需求Uncertain-Demand′h对于SKk的隶属度μhk和type_num个聚类中心{centerk},则有
(9)
设Hk={h|2≤type_num<H;dhk=0},对于所有的h类,h∈Ik,μik=0。则有
(10)
利用式(7)、(8)重复修改{center1,center2,…,centertype_num}、μhk和分类。当算法收敛时模糊聚类结束,从而获得各类别的{center1,center2,…,centertype_num}以及μhk。
假设某装备制造企业工件或组件有10个,Componentm表示的是第m个工件或组件。假设装备制造企业共有Z个加工工序,则Processz表示第z个加工工序,并设ProcessNoz=z,z=1,2,…,Z。设第m个工件或组件需要加工OPNumberm个工序,用OPm(g)表示第m个组件的第g个加工工序。工件或组件加工工艺如表5所示。
表5 工件或组件加工工艺
序号工件或组件工序数量工件或组件的加工工艺路线1组件16模锻→热处理→车→铣→磨→终检2组件27模锻→热处理→车→铣→磨→滚齿→终检3组件36模锻→热处理→车→插→滚齿→终检4组件49铸造→热处理→划线→铣→钻→铣面钻孔→上下箱合箱→铣端面镗中心孔→终检5组件510铸造→热处理→划线→钻→铣→铣面钻孔→滚齿→上下箱合箱→铣端面镗中心孔→终检6组件68铸造→热处理→铣→插→划线→钻→磨→终检7组件78模锻→热处理→车→磨→插→钻→滚齿→终检8组件89下料→锻→热处理→车→铣→钻→精车→划线→终检9组件99热处理→粗车→磨→划线→钻→铣→滚齿→钳工→终检10组件109铸造→清砂→热处理→粗车→精车→检验→滚齿→钳工→终检
计算PosNumi(m1,m2)来表示m1组件第i个工序OPm1(i)在m2组件工序Routingm2中的序号,并以每个组件为基准计算每个组件的Pos-Numi(m1,m2)。表6为以m1=组件1为基准的每个组件的位置号,其他限于篇幅不一一列举。
表6 以m1=组件1为基准的每个组件的位置号PosNumi(1,m2)
工件/组件工序1工序2工序3工序4工序5工序6工序7工序8工序9工序10组件11.01.01.01.01.00 组件21.01.01.01.02.03.50组件31.01.03.06.03.00组件44.54.54.54.54.54.5990组件55.05.05.05.05.05.01010100组件64.04.04.04.01.04.080组件71.01.04.04.04.04.080组件84.51.01.04.54.54.5990组件94.54.54.53.06.04.5990组件104.54.59.09.04.54.5990
计算组件i到组件j的相似性距离,结果如表7所示。
计算组件i到j的相似性,其中产品相似性偏移量dt=1,结果如表8所示。
表7 组件间相似性距离
组件组件1组件2组件3组件4组件5组件6组件7组件8组件9组件10组件103.5937.051221930.037.046.0组件220937.044211130.028.039.0组件397.5046.054251238.039.039.0组件41322.01907223332.529.548.0组件51318.5199.00212536.024.043.5组件61520.51319.53101923.024.043.5组件7812.0637.04718029.023.039.0组件8917.01430.5392426029.541.5组件91317.02231.030192131.0038.0组件101925.51946.046373341.533.50
表8 组件间工艺相似性距离
组件组件1组件2组件3组件4组件5组件6组件7组件8组件9组件10组件11.00000.93270.82690.28850.01920.57690.63460.42310.28850.1154组件20.95561.00000.80000.17780.02220.53330.75560.33330.37780.1333组件30.83640.86361.00000.16360.01820.54550.78180.30910.29090.2909组件40.73470.55100.61221.00000.85710.55100.32650.33670.39800.0204组件50.70790.58430.57300.79781.00000.52810.43820.19100.46070.0225组件60.66290.53930.70790.56180.30341.00000.57300.48310.46070.0225组件70.83330.75000.87500.22920.02080.62501.00000.39580.52080.1875组件80.78820.60000.67060.28240.08240.43530.38821.00000.30590.0235组件90.66670.56410.43590.20510.23080.51280.46150.20511.00000.0256组件100.59570.45740.59570.02130.02130.21280.29790.11700.28721.0000
将平均值最小的行所在的组件作为基准组件,即组件6,以其为对比组件计算组件的相似性。假设取平均值最小的行号为minLine=6,则组件相似性如表9所示。
假设某装备制造企业的客户需求如表10所示。进行归一化处理,结果如表11所示。
表9 组件的工艺相似性
组件1组件2组件3组件4组件5组件6组件7组件8组件9组件100.66290.53930.70790.56180.30341.00000.57300.48310.46070.0225
表10 某装备制造企业的客户需求
序号品种提交时间模糊隶属数量质量模糊隶属11566989102604792242567588302483973861101224045746324461742034854263861746036964384973847023476566786101102378972861061192047119966781031162014710856679210560456111065779710830256125364872835036813267791041157013614438527387502341525872961081034516962759110340357172647693105104711181395475908025719430446983103462074053738540467
表11 某装备制造企业的客户需求归一化处理
序号品种提交时间模糊隶属数量质量模糊隶属10.66290.04120.04340.04620.04686.00000.02940.03170.037020.53930.05490.05350.05490.05423.00000.05880.05540.0416
表11(续)
序号品种提交时间模糊隶属数量质量模糊隶属30.46070.03160.03480.03740.03914.00000.02940.04430.047641.00000.07210.06810.06750.06452.00000.03920.05540.041650.56180.08870.07880.06750.06456.00000.03920.03690.037060.56180.06070.06110.05640.05687.00000.05880.07390.083371.00000.04120.04470.04790.04721.00000.05880.07390.047680.46070.03200.03480.03880.04012.00000.02940.03170.030390.46070.03490.03840.04000.04112.00000.11760.05540.0476100.48310.04120.04470.04470.04546.00000.02940.04430.0555110.02250.03550.03890.04240.04423.00000.05880.04430.0555120.30340.06410.06240.05720.05755.00000.03920.03690.0416130.53930.03440.03790.03960.04157.00000.11760.07390.0555140.56180.06070.05760.05640.05485.00000.05880.07390.0833150.53930.03980.04160.04290.04421.00000.03920.05540.0666160.46070.03720.03990.04520.04634.00000.03920.04430.0476170.53930.03600.03940.04430.04541.00000.02940.03170.0303180.66290.05910.05550.05490.05308.00000.05880.04430.0476190.56180.07690.06810.05960.05751.00000.03920.05540.0555200.57300.05770.05650.05640.05614.00000.02940.03690.0476
采用模糊C-均值聚类算法(FMC)对需求进行聚类。可以将20组需求按序号如下的方法进行聚类:第一类为2、8、9、11;第二类为1、5、6、10、13、18;第三类为3、12、14、16、20;第四类为4、7、15、17、19。
运行之后得到的目标函数值Jb=2.064 8,而且每次运行结果都不一样。这跟初始聚类中的中心点有关系,其中某次聚类结果如图1所示,三角形为各类的聚类中心点。
图1 多产品订单需求模糊C-均值聚类
本文根据产品组件的加工工艺过程进行相似性分析,找到相似性基准组件,并计算各组件之间的相似性。同时,针对订单需求的多重不确定性,包括品种、数量、交货期、质量等,采用模糊C-均值聚类算法对订单需求进行聚类。采用算例对产品订单需求的相似性进行分析,能够解决现实中的订单分组或成组生产问题。未来对相似订单排产问题可以进行进一步研究。
[1]韩文民,陈婷,高龙龙,等.急件订单干扰下虚拟单元重调度 [J].运筹与管理,2018,27(2):68-78.
[2]张于贤,丁修坤.面向多品种小批量的订单成组模型研究 [J].系统科学学报,2018,26(1):92-95.
[3]吴天行,郭键.基于图论的聚类算法在订单分批问题中的应用 [J].物流技术,2017,36(8):112-116.
[4]吴天行,郭键.基于“反学习”理论的人工蜂群算法在订单分批问题中的应用 [J].物流技术,2017,36(12):90-94.
[5]陈俊,梅文辉.基于属性相似性的产品质量模型的建立及优化方法研究 [J].信息技术与信息化,2017(Z1):129-136.
[6]程贤福,李文杰,王浩伦.基于相似性、重用性和定制柔性的产品设计适应性评价方法 [J].现代制造工程,2017(6):156-161.
[7]徐新胜,王诚,肖颖.基于非负矩阵分解的产品结构相似性判断及其应用 [J].中国机械工程,2016,27(8):1072-1077.
[8]胡小林,徐菱.动态并行分区拣选系统下订单分批问题研究 [J].物流技术,2019,38(5):79-87.
[9]禹春霞,刘玉茹,高旭.面向用户个性化需求的云服务商推荐 [J].工业工程,2019,22(2):1-9.
[10] 李若岩,景丽,张友鹏.基于模糊聚类分析石脑油性质相似性研究 [J].兰州交通大学学报,2019,38(5):131-134.
[11] 施伟锋,卓金宝,兰莹.一种基于属性空间相似性的模糊聚类算法 [J].电子与信息学报,2019,41(11):2722-2728.
[12] 王永泉,张广鹏,张倩.基于模糊聚类的可重构制造工艺规划方法 [J].西安理工大学学报,2019,35(2):127-137,124.
[13] 刘丽娜,张青山,徐伟.双渠道采销模式下装备制造企业选址影响因素分析 [J].沈阳工业大学学报(社会科学版),2018,11(4):320-327.
[14] 韩世莲.基于客户动态需求属性的物流配送线路聚类优化 [J].系统管理学报,2016,25(6):1146-1153.
[15] 岳文军.RSCF分析法的特征研究 [J].商场现代化,2015(14):46-54.
[16] 赵改平,刘丽兰,王森,等.运用聚类算法的订单归并模型研究 [J].现代制造工程,2011(12):52-56,68.