CN115685921A - 一种多agv和多机器联合调度的混合整数线性规划模型 - Google Patents

一种多agv和多机器联合调度的混合整数线性规划模型 Download PDF

Info

Publication number
CN115685921A
CN115685921A CN202211319802.XA CN202211319802A CN115685921A CN 115685921 A CN115685921 A CN 115685921A CN 202211319802 A CN202211319802 A CN 202211319802A CN 115685921 A CN115685921 A CN 115685921A
Authority
CN
China
Prior art keywords
foil
roll
agv
machine
time
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.)
Pending
Application number
CN202211319802.XA
Other languages
English (en)
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.)
Liaocheng University
Original Assignee
Liaocheng University
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 Liaocheng University filed Critical Liaocheng University
Priority to CN202211319802.XA priority Critical patent/CN115685921A/zh
Publication of CN115685921A publication Critical patent/CN115685921A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Metal Rolling (AREA)

Abstract

一种多AGV和多机器联合调度的混合整数线性规划模型,涉及车间调度技术领域,特别是属于一种在自动化电极箔生产车间的多AGV与多机器联合调度的混合整数线性规划模型。包括步骤1,定义混合整数线性规划模型的模型参数和决策变量,将求测电极箔生产所需的最大完工时间的最小值Minimise Cmax,作为混合整数线性规划模型输出的目标值;步骤2,在确定优化目标的基础上,对混合整数线性规划模型在时间和AGV容量上进行线性约束;步骤3,对建立的自动化电极箔生产车间的多AGV与多机器调度的混合整数线性规划模型,运用GUROBI求解器进行验证并根据实例得到甘特图。本发明能够实现AGV最大运载量的充分利用,具有提高电极箔生产效率的积极效果。

Description

一种多AGV和多机器联合调度的混合整数线性规划模型
技术领域
本发明涉及车间调度技术领域,特别是属于一种在自动化电极箔生产车间的多AGV与多机器联合调度的混合整数线性规划模型。
背景技术
在即将到来的工业4.0时代,自动导引车AGV作为一种灵活高效的运输设备被广泛应用于制造业。在自动化电极箔生产车间中,AGV被用来运输箔卷。整个生产线可描述如下:待加工的箔卷全部存放在物料库,需要由AGV运输至机器进行加工,加工完成后,AGV将成品运输至成品仓库。为了提高自动化电极箔生产车间的生产效率,就必须保证车间内的箔卷分配合理,能够在最快的时间内完成加工,这就需要制定高效的调度解决方案。在电极箔生产过程中,涉及两个调度问题:一个是将箔卷分配给哪个AGV进行运输,另一个是将箔卷分配到哪个机器进行加工处理,即AGV调度问题和并行机调度问题。综合考虑上述两个问题是真正提高自动电极箔生产厂效率的唯一途径,即考虑AGV和并行机的联合调度。
在车间调度研究领域,通常会以最小化最大完工时间作为研究目标,即一批工件被分配到不同的机器上进行加工,每个机器上所加工的最后一个工件加工完成的时间即该机器的完工时间,那么所有机器中最大的完工时间即为这批工件的完工时间,以最小化最大完工时间为目标是为了制定高效的机器调度方案来提高整个生产车间的生产效率。在自动化电极箔生产车间调度研究中,箔卷被AGV运输至加工机器的时间越短,机器开始加工的时间越早,进而箔卷从仓库运输直至加工完成的生产周期越短,由此将AGV的运输时间考虑到智能生产线的调度中更加贴合实际生产需要。但目前的自动化电极箔生产车间调度方案中,大多数研究只考虑机器的调度问题,即只考虑如何将箔卷更合理的分配给机器来得到最小化最大完工时间,忽略了箔卷如何搬运到机器的运输时间。由此,建立带有运输时间的最大完工时间最小为优化目标的问题模型,是当前自动化电极箔生产车间多AGV与多机器联合调度研究中亟待解决的问题。
发明内容
本发明的目的即在于提供一种多AGV和多机器联合调度的混合整数线性规划模型,以解决如何贴合实际生产需要满足电极箔生产的最大完工时间最小的优化目标的问题。
本发明所提供的一种多AGV和多机器联合调度的混合整数线性规划模型,其特征在于,包括以下步骤,
步骤1,定义混合整数线性规划模型的模型参数和决策变量,将求测电极箔生产所需的最大完工时间的最小值Minimise Cmax,作为混合整数线性规划模型输出的目标值;
步骤2,在确定优化目标的基础上,对混合整数线性规划模型在时间和AGV容量上进行线性约束;
步骤3,对建立的自动化电极箔生产车间的多AGV与多机器调度的混合整数线性规划模型,运用GUROBI求解器进行验证并根据实例得到甘特图。
进一步的,模型参数包括以下参数,
J-箔卷的数量;
j、j′-箔卷的编号,其中,j、j′∈[0,J];
M-机器的数量;
m、m′-机器的编号,其中,m、m′∈[1,M];
K-AGV的数量;
k-AGV的编号,其中,k∈[1,K];
xm、ym-机器m的位置坐标;
tu-卸货时间,tu=15s;
pj-箔卷j的加工处理时间;
Q-AGV的最大载重量;
v-AGV的速度;
P-AGV最多运输的趟数,其中,P=J;
p-AGV运输的第p趟;
N-AGV某一趟运输的箔卷总数;
n-AGV在某一趟运输到第n个位置;
h-一个无限大的数;
Tm,m'-AGV从机器m到机器m'的运输时间;
决策变量包括以下变量,
tj,m-有0和1两种取值,若箔卷j被分配在编号为m的机器上进行加工,tj,m=1;若箔卷j没有被分配给编号为m的机器加工,即分给机器m以外的机器加工,tj,m=0;
uj,j',m-有0和1两种取值,若编号为j和j'的箔卷都被安排到编号为m的机器上加工,并且箔卷j'是j的后继,即箔卷j'在j后面进行加工,则uj,j',m=1;若编号为j和j'的箔卷,其中有一个没有分配给编号为m的机器,或者两个箔卷都分配给编号为m的机器,但箔卷j'在j前面加工,则uj,j′,m=0;
wj,k,p,n-有0和1两种取值,若编号为j的箔卷安排到编号为k的AGV小车上第p趟运输的第n个位置,则wj,k,p,n=1,若满足下述三种情况:(1)箔卷j不被编号为k的AGV运输;(2)箔卷j不被编号为k的AGV在第p趟运输;(3)箔卷j不在编号为k的AGV在第p趟运输的第n个位置,则wj,k,p,n=0;
zk,p,n-有0和1两种取值;如果编号为k的AGV小车上第p趟运输的第n个位置为最后放置箔卷的位置,则zk,p,n=1;如果第n+1位置上还有箔卷,则zk,p,n=0。
进一步的,在步骤1中,电极箔加工所需的最大完工时间通过以下公式(1)-公式(5)计算,
Figure BDA0003909842080000041
Figure BDA0003909842080000042
Figure BDA0003909842080000043
Figure BDA0003909842080000044
Cmax≥cj′,j′=1,2,…,J (5)
其中,公式(1)表示AGV从机器m到机器m′的运输时间Tm,m′;公式(2)表示箔卷j和j′由同一台AGV运输时,且箔卷j的顺序
在箔卷j′前一个,箔卷j′到达加工机器m的到达时间Aj',m;公式(3)表示箔卷j'到达在机器m上开始加工的时间Sj′,m;公式(4)表示箔卷j'的完工时间cj';公式(5)表示电极箔加工所需的最大完工时间Cmax
进一步的,对未加工的箔卷进行机器分配时,需要保证一个箔卷的整体性,即一个箔卷同一时刻只能在一台机器上进行加工处理,而且一个箔卷一旦开始加工必须加工完成,不可停止,约束条件如公式(6)所示;两个未加工的箔卷在同一台机器上进行加工处理时按先后顺序进行,约束条件如公式(7)-公式(8)所示,
Figure BDA0003909842080000051
Figure BDA0003909842080000052
Figure BDA0003909842080000053
进一步的,某一台机器上需要加工处理的箔卷序列的第一个箔卷之前和最后一个箔卷之后都设置一个虚拟箔卷,虚拟箔卷记为箔卷0,那么每一个虚拟箔卷在机器m上有唯一的前驱和和唯一的后继,约束条件如公式(9)-公式(10)所示,
Figure BDA0003909842080000054
Figure BDA0003909842080000055
进一步的,一个箔卷同一时刻只能在一台AGV的某一趟的一个位
置上进行运输配送,从箔卷仓库运输至加工机器,中途不可卸下,约束条件如公式(11)所示;两个未加工的箔卷被同一台AGV运输一定有先后顺序的,即被安排的某一个小车的某一趟运输配送的某一个位置,最多只能放置1个箔卷,约束条件如公式(12)所示;如果在编号为k的AGV小车第p趟运输的第n个位置没有箔卷,则本趟运输的下一个位置也不能有箔卷约束条件如公式(13)所示;如果在编号为k的AGV小车AGV小车第p趟没有运输箔卷,则下一趟也不能运输箔卷(也就是说不存在下一趟),约束条件如公式(14)所示;判断编号为k的AGV小车上第p趟运输的最后一个位置(第N个位置)是否为最后放置箔卷的位置以及最后一趟运输是否只有虚拟箔卷,约束条件如公式(15)-(16)所示,
Figure BDA0003909842080000061
Figure BDA0003909842080000062
Figure BDA0003909842080000063
Figure BDA0003909842080000064
Figure BDA0003909842080000065
Figure BDA0003909842080000066
进一步的,箔卷在机器上进行加工处理时,箔卷的完工时间是根据其加工顺序决定的,如果箔卷j'是机器上加工的第一个箔卷,约束条件如公式(17)所示;如果箔卷j'是j的加工后继,约束条件如公式(18)所示,箔卷j的完工时间与箔卷的到达时间、卸货时间以及箔卷的加工时间的约束条件,如公式(19)所示,
Figure BDA0003909842080000071
Figure BDA0003909842080000072
cj≥Aj,m+tu+pj,j=1,2,…,J (19)。
进一步的,箔卷到达加工机器的时间包括AGV的运输时间和卸货时间,约束条件如公式(20)-公式(21)所示,
Figure BDA0003909842080000073
Aj′,m′≥Aj,m+tu+Tm,m′+(tj,m+tj′,m′+wj,k,p,n+wj′,k,p,n+1-4)·h (21)
其中,j′=1,2,…,J,j=1,2,…,J,j′≠j,m=1,2,…M,m′=1,2,…,M,k=1,2,…,K,p=1,2,…,P,n=1,2,…,N-1。
进一步的,目标值最大完工时间大于等于任何一个箔卷的完工时间,约束条件如公式(22)所示,
Cmax≥cj,j=1,2,…,J (22)。
进一步的,运用GUROBI求解器进行验证,具体包括以下过程,
获取各个箔卷任务的编号、重量和加工时间,以及AGV的最大负载量;
将上述数据输入GUROBI求解器中,得到最优的生产调度顺序。
本发明所提供的一种多AGV和多机器联合调度的混合整数线性规划模型,针对自动化电极箔生产线建立带有运输时间的最大完工时间最小为优化目标的问题模型,并通过GUROBI求解器进行验证。在自动化电极箔生产车间的AGV与机器的联合调度问题中通过AGV利用率的考虑,实现运输过程中AGV的最大运载量的充分利用,具有提高电极箔的生产效率的积极效果。
附图说明
图1为本发明构建的矩阵车间多AGV与多机器联合调度的场景图;
图2为本发明的原理图;
图3为本发明的小规模算例中6个箔卷的参数设置表;
图4为本发明的GUROBI求解器求解结果甘特图。
具体实施方式
如图1-2所示,本发明所提供的一种多AGV和多机器联合调度的混合整数线性规划模型,适用于矩阵生产车间的自动化电极箔生产线。执行操作设备,如电解和箔卷干燥等机器,均按照图1所示的矩阵分布进行排列,未加工的铝箔卷全部放置在物料库,需要由AGV运输到机器进行加工,AGV将成品电极箔运输至成品库。本发明的目标值为最小化最大完工时间,即计算所有机器上最后一个箔卷加工完成的时
间,不考虑将成品运回成品库这一阶段。
如图2所示,本发明的具体实现步骤如下:
步骤1,定义混合整数线性规划模型的模型参数和决策变量,将求测电极箔加工所需的最大完工时间的最小值Minimise Cmax,作为混合整数线性规划模型输出的目标值。具体地,模型参数包括以下参数设置,
J-箔卷的数量;
j、j′-箔卷的编号,其中,j、j′∈[0,J];
M-机器的数量;
m、m′-机器的编号,其中,m、m′∈[1,M];
K-AGV的数量;
k-AGV的编号,其中,k∈[1,K];
xm、ym-机器m的位置坐标;
tu-卸货时间,tu=15s;
pj-箔卷j的加工处理时间;
Q-AGV的最大载重量;
v-AGV的速度;
P-AGV最多运输的趟数,其中,P=J;
p-AGV运输的第p趟;
N-AGV某一趟运输的箔卷总数;
n-AGV在某一趟运输到第n个位置;
h-一个无限大的数;
Tm,m′-AGV从机器m到机器m′的运输时间;
决策变量包括以下变量设置,
tj,m-有0和1两种取值,若箔卷j被分配在编号为m的器上进行加工,tj,m=1;若箔卷j没有被分配给编号为m的机器加工,即分给机器m以外的机器加工,tj,m=0;
uj,j',m-有0和1两种取值,若编号为j和j'的箔卷都被安排到编号为m的机器上加工,并且箔卷j′是j的后继,即箔卷j'在j后面进行加工,则uj,j',m=1;若编号为j和j'的箔卷,其中有一个没有分配给编号为m的机器,或者两个箔卷都分配给编号为m的机器,但箔卷j'在j前面加工,则uj,j',m=0;
wj,k,p,n-有0和1两种取值,若编号为j的箔卷安排到编号为k的AGV小车上第p趟运输的第n个位置,则wj,k,p,n=1,若满足下述三种情况:(1)箔卷j不被编号为k的AGV运输;(2)箔卷j不被编号为k的AGV在第p趟运输;(3)箔卷j不在编号为k的AGV在第p趟运输的第n个位置,则wj,k,p,n=0;
zk,p,n-有0和1两种取值;如果编号为k的AGV小车上第p趟运输的第n个位置为最后放置箔卷的位置,则zk,p,n=1;如果第n+1位置上还有箔卷,则zk,p,n=0。
另外,在步骤1中,电极箔加工所需的最大完工时间通过以下公式(1)-公式(5)计算,
Figure BDA0003909842080000101
Figure BDA0003909842080000102
Figure BDA0003909842080000111
Figure BDA0003909842080000112
Cmax≥cj',j'=1,2,…,J (5)
其中,公式(1)表示AGV从机器m到机器m′的运输时间Tm,m′;公式(2)表示箔卷j和j′由同一台AGV运输时,且箔卷j的顺序在箔卷j′前一个,箔卷j′到达加工机器m的到达时间Aj′,m;公式(3)表示箔卷j′到达在机器m上开始加工的时间Sj′,m;公式(4)表示箔卷j′的完工时间cj′;公式(5)表示电极箔加工所需的最大完工时间Cmax
步骤2,在确定优化目标的基础上,对混合整数线性规划模型在时间和AGV容量上进行线性约束。具体地线性约束条件如下所述,
(1)箔卷j只能被安排的一个机器上加工,即,
Figure BDA0003909842080000113
(2)如果箔卷j′在机器m上加工,则箔卷j′在机器m上有唯一的前驱(不能是自己),也就是说在机器m上,与箔卷j′相邻且位置在箔卷j′前边的箔卷j是唯一的。如果箔卷j′是该机器上的第一个加工的箔卷,那么它的前驱是虚拟箔卷0,即,
Figure BDA0003909842080000114
(3)如果箔卷j在机器m上加工,则箔卷j在机器m上有唯一的后继(不能是自己),也就是说在机器m上,与箔卷j相邻且位置在箔卷j后边的箔卷j′是唯一的。如果箔卷j是该机器上的最后一个加工的箔卷,那么它的后继是虚拟箔卷0,即,
Figure BDA0003909842080000121
(4)对虚拟箔卷(箔卷0)的前驱和后继的约束。箔卷0在机器m上有唯一的前驱(可以是自己,表示机器m上没有加工箔卷,即u0,0,m=1),即,
Figure BDA0003909842080000122
(5)箔卷0在机器m上有唯一的后继(可以是自己,表示机器m上没有加工箔卷,即u0,0,m=1),即,
Figure BDA0003909842080000123
(6)箔卷j只能被安排某一个小车的某一趟运输的某一个位置,即,
Figure BDA0003909842080000124
(7)被安排某一个小车的某一趟运输的某一个位置最多只能放置1个箔卷,即,
Figure BDA0003909842080000125
(8)因为每台AGV的最大负载量Q=200,要保证AGV的单趟运输不超载,即编号为k的AGV第p趟运输的所有位置的箔卷的重量之和不超过Q,即,
Figure BDA0003909842080000131
(9)如果在编号为k的AGV小车第p趟运输的第n个位置没有箔卷,则本趟运输的下一个位置也不能有箔卷。例如,该AGV这一趟只运输3个箔卷,第四个位置是虚拟箔卷0,那么第五个位置也不可能有箔卷,排序的时候就是1,2,3的顺序,不会是1,2,4,即,
Figure BDA0003909842080000132
(10)如果在编号为k的AGV小车AGV小车第p趟运输没有运输箔卷,则下一趟运输也不能运输箔卷,即,
Figure BDA0003909842080000133
(11)如果箔卷j′是机器上加工的第一个箔卷,则箔卷j′的完工时间大于等于箔卷j′的加工时间pj,即,
Figure BDA0003909842080000134
当箔卷j′为机器上加工的第一个箔卷,那么
Figure BDA0003909842080000135
此时cj′≥pj′
(12)如果箔卷j′是j的加工后继,则工件j′的完工时间大于等于箔卷j的完工时间加上箔卷j′的加工时间pj,即,
Figure BDA0003909842080000136
Figure BDA0003909842080000137
当箔卷j′是j的加工后继,
Figure BDA0003909842080000138
那么cj′≥cj+pj′
(13)箔卷j的完工时间大于等于箔卷j的到达时间加上卸货时间加上箔卷j的加工时间,即,
cj≥Aj,m+tu+pj,j=1,2,…,J,
第一个箔卷到达机器后就开始加工处理,此时cj=Aj,m+tu+pj;其余的箔卷分为两种情况:一种是箔卷到达,但上一个箔卷并未加工完,需要等待,此时cj>Aj,m+tu+pj;另一种是上一个箔卷加工完成,机器空闲,那么cj=Aj,m+tu+pj
(14)如果箔卷j是AGV小车k第1趟运输的第1个位置上的箔卷,则箔卷j的到达时间大于等于AGV小车k从原点到达加工工件j的机器的时间加上卸货时间,即,
Figure BDA0003909842080000141
因为如果箔卷j是AGV小车k第1趟运输的第1个位置上的箔卷,那么
Figure BDA0003909842080000142
tj,m=1,即Aj,m≥Tm,m′+tu
(15)如果箔卷j在AGV小车k的第p趟运输的第n个位置,j′在AGV小车k的第p趟运输的第n+1个位置,则箔卷j′的到达时间大于等于箔卷j的到达时间加上卸货时间加上AGV小车从加工箔卷j的机器到达加工箔卷j′的机器的时间,即,
Aj′,m′≥Aj,m+tu+Tm,m′+(tj,m+tj′,m′+wj,k,p,n+wj′,k,p,n+1-4)·h
其中,j′=1,2,…,J,j=1,2,…,J,j′≠j,m=1,2,…M,m′=1,2,…,M,k=1,2,…,K,p=1,2,…,P,n=1,2,…,N-1。
(16)判断编号为k的AGV小车上第p趟运输的最后一个位置(第N个位置)是否为最后放置箔卷的位置,即,
Figure BDA0003909842080000143
判断编号为k的AGV小车上第p趟运输除最后一个位置外的其他位置是否为最后放置箔卷的位置,即,
Figure BDA0003909842080000151
Figure BDA0003909842080000152
如果位置n是最后一个,那么
Figure BDA0003909842080000153
Figure BDA0003909842080000154
所以zk,p,n=1。
(18)如果箔卷j排序在AGV小车k的第p趟运输的最后1个位置,j′排序在AGV小车k的第p+1趟运输的第1个位置,则工件j′的到达时间大于等于工件j的到达时间+卸货时间+AGV小车从加工箔卷j的机器回到原点的时间+AGV小车从原点到达加工箔卷j'的机器的时间,即,
Figure BDA0003909842080000155
Figure BDA0003909842080000156
其中,j'=1,2,…,J,j=1,2,…,J,j'≠j,k=1,2,…,K,p=1,2,…,P-1,n=1,2,…,N。
(19)目标值Cmax是指所有机器上最后一个箔卷加工完成的时间,所以它应该大于等于任意一个箔卷的完工时间,即,
Cmax≥cj,j=1,2,…,J。
步骤3,对建立的自动化电极箔生产车间的多AGV与多机器调度的混合整数线性规划模型,运用GUROBI求解器进行验证并根据实例得到甘特图。
为了验证该MILP模型的正确性,本发明通过具体实例进行验证,考虑到时间限制以及解空间的搜索,在本发明的具体实施例中,采用小规模算例-6个箔卷,并运用GUROBI求解器进行验证,验证过程包括:获取各个箔卷任务的编号、重量和加工时间,以及AGV的最大负载量;将上述数据输入GUROBI求解器中,得到最优的生产调度顺序。通过仿真内容对本发明做进一步的详细说明。本实验的箔卷任务规模
为6的小规模算例的调度模型,采用AGV的最大载重量为200KG,假设矩阵车间内有4台可用机器,两台可用AGV。6个箔卷任务的参数值的具体设置如图3所示,而后通过GUROBI求解器可得:箔卷1和箔卷4在机器1上加工,箔卷3在箔卷2在机器2上加工;箔卷5在机器3上加工,箔卷2和箔卷6在机器4上加工。AGV1运输箔卷(1,3,4),配送顺序是(3,4,1),AGV2运输箔卷(2,5,6),配送顺序是(6,5,2)。由图4所示的GUROBI求解器求解结果甘特图,可以看出甘特图中机器的等待时间为0,说明GUROBI求解器找到了最优解,由此证明了本发明所提出的一种在自动化电极箔生产车间的多AGV与多机器调度的混合整数线性规划模型的正确性,符合实际生产需要。

Claims (10)

1.一种多AGV和多机器联合调度的混合整数线性规划模型,其特征在于,包括以下步骤,
步骤1,定义混合整数线性规划模型的模型参数和决策变量,将求测电极箔生产所需的最大完工时间的最小值Minimise Cmax,作为混合整数线性规划模型输出的目标值;
步骤2,在确定优化目标的基础上,对混合整数线性规划模型在时间和AGV容量上进行线性约束;
步骤3,对建立的自动化电极箔生产车间的多AGV与多机器调度的混合整数线性规划模型,运用GUROBI求解器进行验证并根据实例得到甘特图。
2.根据权利要求1所述的一种多AGV和多机器联合调度的混合整数线性规划模型,其特征还在于,模型参数包括以下参数,
J-箔卷的数量;
j、j'-箔卷的编号,其中,j、j'∈[0,J];
M-机器的数量;
m、m'-机器的编号,其中,m、m'∈[1,M];
K-AGV的数量;
k-AGV的编号,其中,k∈[1,K];
xm、ym-机器m的位置坐标;
tu-卸货时间,tu=15s;
pj-箔卷j的加工处理时间;
Q-AGV的最大载重量;
v-AGV的速度;
P-AGV最多运输的趟数,其中,P=J;
p-AGV运输的第p趟;
N-AGV某一趟运输的箔卷总数;
n-AGV在某一趟运输到第n个位置;
h-一个无限大的数;
Tm,m'-AGV从机器m到机器m'的运输时间;
决策变量包括以下变量,
tj,m-有0和1两种取值,若箔卷j被分配在编号为m的机器上进行加工,tj,m=1;若箔卷j没有被分配给编号为m的机器加工,即分给机器m以外的机器加工,tj,m=0;
uj,j',m-有0和1两种取值,若编号为j和j'的箔卷都被安排到编号为m的机器上加工,并且箔卷j'是j的后继,即箔卷j'在j后面进行加工,则uj,j',m=1;若编号为j和j'的箔卷,其中有一个没有分配给编号为m的机器,或者两个箔卷都分配给编号为m的机器,但箔卷j'在j前面加工,则uj,j',m=0;
wj,k,p,n-有0和1两种取值,若编号为j的箔卷安排到编号为k的AGV小车上第p趟运输的第n个位置,则wj,k,p,n=1,若满足下述三种情况:(1)箔卷j不被编号为k的AGV运输;(2)箔卷j不被编号为k的AGV在第p趟运输;(3)箔卷j不在编号为k的AGV在第p趟运输的第n个位置,则wj,k,p,n=0;
zk,p,n-有0和1两种取值;如果编号为k的AGV小车上第p趟
运输的第n个位置为最后放置箔卷的位置,则zk,p,n=1;如果第n+1位置上还有箔卷,则zk,p,n=0。
3.根据权利要求2所述的一种多AGV和多机器联合调度的混合整数线性规划模型,其特征还在于,电极箔加工所需的最大完工时间通过以下公式(1)-公式(5)计算,
Figure FDA0003909842070000031
Figure FDA0003909842070000032
Figure FDA0003909842070000033
Figure FDA0003909842070000034
Cmax≥cj',j'=1,2,…,J (5)
其中,公式(1)表示AGV从机器m到机器m'的运输时间Tm,m';公式(2)表示箔卷j和j'由同一台AGV运输时,且箔卷j的顺序在箔卷j'前一个,箔卷j'到达加工机器m的到达时间Aj',m;公式(3)表示箔卷j'到达在机器m上开始加工的时间Sj′,m;公式(4)表示箔卷j'的完工时间cj';公式(5)表示电极箔加工所需的最大完工时间Cmax
4.根据权利要求3所述的一种多AGV和多机器联合调度的混合整数线性规划模型,其特征还在于,对未加工的箔卷进行机器分配时,需要保证一个箔卷的整体性,即一个箔卷同一时刻只能在一台机器上进行加工处理,而且一个箔卷一旦开始加工必须加工完成,不可停止,约束条件如公式(6)所示;两个未加工的箔卷在同一台机器上进行加工处理时按先后顺序进行,约束条件如公式(7)-公式(8)所示,
Figure FDA0003909842070000041
Figure FDA0003909842070000042
Figure FDA0003909842070000043
5.根据权利要求3所述的一种多AGV和多机器联合调度的混合整数线性规划模型,其特征还在于,某一台机器上需要加工处理的箔卷序列的第一个箔卷之前和最后一个箔卷之后都设置一个虚拟箔卷,虚拟箔卷记为箔卷0,那么每一个虚拟箔卷在机器m上有唯一的前驱和和唯一的后继,约束条件如公式(9)-公式(10)所示,
Figure FDA0003909842070000044
Figure FDA0003909842070000045
6.根据权利要求3所述的一种多AGV和多机器联合调度的混合整数线性规划模型,其特征还在于,一个箔卷同一时刻只能在一台AGV的某一趟的一个位置上进行运输配送,从箔卷仓库运输至加工机器,中途不可卸下,约束条件如公式(11)所示;两个未加工的箔卷被同一台AGV运输一定有先后顺序的,即被安排的某一个小车的某一趟运输配送的某一个位置,最多只能放置1个箔卷,约束条件如公式(12)所示;如果在编号为k的AGV小车第p趟运输的第n个位置没有箔卷,则本趟运输的下一个位置也不能有箔卷约束条件如公式(13)所示;如果在编号为k的AGV小车AGV小车第p趟没有运
输箔卷,则下一趟也不能运输箔卷(也就是说不存在下一趟),约束条件如公式(14)所示;判断编号为k的AGV小车上第p趟运输的最后一个位置(第N个位置)是否为最后放置箔卷的位置以及最后一趟运输是否只有虚拟箔卷,约束条件如公式(15)-(16)所示,
Figure FDA0003909842070000051
Figure FDA0003909842070000052
Figure FDA0003909842070000053
Figure FDA0003909842070000054
Figure FDA0003909842070000055
Figure FDA0003909842070000056
7.根据权利要求3所述的一种多AGV和多机器联合调度的混合整数线性规划模型,其特征还在于,箔卷在机器上进行加工处理时,箔卷的完工时间是根据其加工顺序决定的,如果箔卷j'是机器上加工的第一个箔卷,约束条件如公式(17)所示;如果箔卷j′是j的加工后继,约束条件如公式(18)所示,箔卷j的完工时间与箔卷的到达时间、卸货时间以及箔卷的加工时间的约束条件,如公式(19)所示,
Figure FDA0003909842070000061
Figure FDA0003909842070000062
cj≥Aj,m+tu+pj,j=1,2,…,J (19)。
8.根据权利要求3所述的一种多AGV和多机器联合调度的混合整数线性规划模型,其特征还在于,箔卷到达加工机器的时间包括AGV的运输时间和卸货时间,约束条件如公式(20)-公式(21)所示,
Figure FDA0003909842070000063
Aj',m'≥Aj,m+tu+Tm,m'+(tj,m+tj',m'+wj,k,p,n+wj′,k,p,n+1-4)·h (21)
其中,j′=1,2,…,J,j=1,2,…,J,j′≠j,m=1,2,…M,m′=1,2,…,M,k=1,2,…,K,p=1,2,…,P,n=1,2,…,N-1。
9.根据权利要求3所述的一种多AGV和多机器联合调度的混合整数线性规划模型,其特征还在于,目标值最大完工时间大于等于任何一个箔卷的完工时间,约束条件如公式(22)所示,
Cmax≥cj,j=1,2,…,J (22)。
10.根据权利要求4-9中任意一项所述的一种多AGV和多机器联合调度的混合整数线性规划模型,其特征还在于,运用GUROBI求解器进行验证,具体包括以下过程,
获取各个箔卷任务的编号、重量和加工时间,以及AGV的最大负载量;
将上述数据输入GUROBI求解器中,得到最优的生产调度顺序。
CN202211319802.XA 2022-10-26 2022-10-26 一种多agv和多机器联合调度的混合整数线性规划模型 Pending CN115685921A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211319802.XA CN115685921A (zh) 2022-10-26 2022-10-26 一种多agv和多机器联合调度的混合整数线性规划模型

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211319802.XA CN115685921A (zh) 2022-10-26 2022-10-26 一种多agv和多机器联合调度的混合整数线性规划模型

Publications (1)

Publication Number Publication Date
CN115685921A true CN115685921A (zh) 2023-02-03

Family

ID=85099235

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211319802.XA Pending CN115685921A (zh) 2022-10-26 2022-10-26 一种多agv和多机器联合调度的混合整数线性规划模型

Country Status (1)

Country Link
CN (1) CN115685921A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116736820A (zh) * 2023-08-16 2023-09-12 聊城大学 一种考虑加工、运输和装配集成调度的建模方法
CN116736819A (zh) * 2023-08-14 2023-09-12 聊城大学 一种基于混合整数规划模型的fjsp-agv的建模方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116736819A (zh) * 2023-08-14 2023-09-12 聊城大学 一种基于混合整数规划模型的fjsp-agv的建模方法
CN116736819B (zh) * 2023-08-14 2023-10-20 聊城大学 一种基于混合整数规划模型的fjsp-agv的建模方法
CN116736820A (zh) * 2023-08-16 2023-09-12 聊城大学 一种考虑加工、运输和装配集成调度的建模方法
CN116736820B (zh) * 2023-08-16 2023-10-17 聊城大学 一种考虑加工、运输和装配集成调度的建模方法

Similar Documents

Publication Publication Date Title
CN115685921A (zh) 一种多agv和多机器联合调度的混合整数线性规划模型
CN112561194B (zh) 一种混合流水车间生产与物流集成调度方法及系统
CN110598941A (zh) 一种基于仿生策略的粒子群优化制造系统双目标排产方法
CN105956689A (zh) 一种基于改进粒子群优化的运输和生产协同调度方法
CN109383969B (zh) 一种汽车零部件货到人的方法
Jiang et al. Uniform parallel batch machines scheduling considering transportation using a hybrid DPSO-GA algorithm
CN115034444B (zh) 基于学习效应的多目标双柔性作业车间调度方法及系统
CN109823757A (zh) 一种板件出库方法、系统及存储介质
CN115657616B (zh) 一种基于agv调度系统的任务分配方法
CN114022028B (zh) 一种自动化混合流水线调度布局集成优化方法
CN110750079A (zh) 一种允许工序跳跃的混合流水车间调度优化方法
CN116258308A (zh) 一种基于混合遗传算法的动态柔性作业车间调度方法
CN111340268A (zh) 一种多工序加工单元移动机器人的调度优化方法
CN113205277B (zh) 一种基于时空规则的车间生产与天车协同调度方法及装置
CN113177781B (zh) 基于变邻域和遗传算子的生产组装协同调度方法及系统
CN113110331B (zh) 一种订单式机器人冲压单元排版选配与调度集成优化方法
CN102478841B (zh) 一种平衡机台产能的方法和装置
CN111007821B (zh) 批加工时间受总重量和最大尺寸双重约束的车间调度方法
Zhang et al. Layout design of a mixed-flow production line based on processing energy consumption and buffer configuration
Yan et al. Flexible Job-Shop Scheduling with Setups and Variable Sublots
JP5444801B2 (ja) 板状鉄鋼製品の熱処理炉への投入計画作成方法及び投入計画作成装置
CN104915820A (zh) 运输工具虚拟为供应点或需求点的物资调度方法及系统
JPS62173155A (ja) 投入制御最適化方式
Ge et al. Research on online scheduling method for flexible assembly workshop of multi-agv system based on assembly island mode
CN116477263B (zh) 一种板式定制家具自动分拣方法

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