[go: up one dir, main page]

CN115130744A - A production scheduling method for a job shop with parallel machines under shift constraints - Google Patents

A production scheduling method for a job shop with parallel machines under shift constraints Download PDF

Info

Publication number
CN115130744A
CN115130744A CN202210728425.9A CN202210728425A CN115130744A CN 115130744 A CN115130744 A CN 115130744A CN 202210728425 A CN202210728425 A CN 202210728425A CN 115130744 A CN115130744 A CN 115130744A
Authority
CN
China
Prior art keywords
shift
equipment
time
solution
job shop
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
CN202210728425.9A
Other languages
Chinese (zh)
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.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN202210728425.9A priority Critical patent/CN115130744A/en
Publication of CN115130744A publication Critical patent/CN115130744A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/04Manufacturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • General Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Manufacturing & Machinery (AREA)
  • Primary Health Care (AREA)
  • Educational Administration (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention belongs to the technical field of workshop scheduling, and particularly discloses a scheduling method of a workshop with parallel machines under constraint of shift, which comprises the following steps: s1, constructing a job shop scheduling model under constraint conditions by taking the minimum lag as a target, wherein the constraint conditions comprise a shift constraint and a workpiece priority constraint; the shift constraint is as follows: one working procedure of the workpiece can only finish processing within one shift on the selected equipment, and the processing finishing time does not exceed the finishing time of the shift; the workpiece priority constraints are: all workpieces are divided into emergency groups and common groups, the emergency groups can not be delivered in an overdue mode, and the common groups can be delivered in an overdue mode; and S2, solving the job shop scheduling model according to the production information of the job shop with the parallel machine to obtain the optimal scheduling scheme of the shop. The invention can generate a better scheduling scheme according to the scheduling target, and can reasonably arrange overtime to ensure timely delivery of the emergency.

Description

一种班次约束下带并行机作业车间的排产方法A production scheduling method for a job shop with parallel machines under shift constraints

技术领域technical field

本发明属于车间调度技术领域,更具体地,涉及一种班次约束下带并行机作业车间的排产方法。The invention belongs to the technical field of workshop scheduling, and more particularly relates to a production scheduling method for a workshop with parallel machines under shift constraints.

背景技术Background technique

智能制造是全球制造业发展的趋势,其中的车间调度技术直接影响着制造的效率与成本。随着数控机床、加工中心等多功能设备的广泛应用和个性化、定制化生产模式的兴起,柔性作业车间调度问题也成为了研究热点。Smart manufacturing is the development trend of the global manufacturing industry, and the workshop scheduling technology directly affects the efficiency and cost of manufacturing. With the wide application of multi-functional equipment such as CNC machine tools and machining centers and the rise of individualized and customized production models, flexible job shop scheduling has also become a research hotspot.

带并行机的作业车间作为一种特殊的柔性作业车间,广泛存在于实际生产中。由于在带并行机的作业车间调度问题中,一个设备组内的每台设备加工同一道工序的时长可能是相等的,如此在设备选择子问题上,求解经典柔性作业车间调度问题的算法可能不再适用,故需要在此基础上探索新的求解方法。As a special kind of flexible workshop, the workshop with parallel machine widely exists in actual production. Since in the job shop scheduling problem with parallel machines, the processing time of each equipment in an equipment group may be equal, so on the equipment selection sub-problem, the algorithm for solving the classical flexible job shop scheduling problem may not be It is applicable again, so it is necessary to explore a new solution method on this basis.

此外,现实生产车间中,加工设备往往需要工人来操作,因此需要设计工人时间表,以规定工人上下班时间。在流水车间里,通常采用工人三班倒的方式保证生产线的连续性。但在柔性作业车间中,因为没有生产线的概念,不需要保持24小时不间断加工,因而通常会给工人安排休息时间,采用工人时间表将设备生产时间按班次进行划分,如此也对调度过程产生了极大影响。因此,班次约束是现实车间调度问题中必须要考虑的因素。In addition, in the real production workshop, the processing equipment often needs workers to operate, so it is necessary to design the worker schedule to stipulate the workers' commute time. In the assembly line, workers usually work in three shifts to ensure the continuity of the production line. However, in the flexible workshop, because there is no concept of a production line, there is no need to maintain 24-hour uninterrupted processing. Therefore, workers are usually arranged for rest time, and the production time of equipment is divided by shifts using the worker schedule, which also affects the scheduling process. great influence. Therefore, the shift constraint is a factor that must be considered in the real shop scheduling problem.

发明内容SUMMARY OF THE INVENTION

针对现有技术的以上缺陷或改进需求,本发明提供了一种班次约束下带并行机作业车间的排产方法,其目的在于,解决班次约束下带并行机的作业车间调度问题中,工件拖期严重的问题,提高车间生产效率。In view of the above defects or improvement needs of the prior art, the present invention provides a production scheduling method for a job shop with parallel machines under shift constraints, the purpose of which is to solve the problem of job shop scheduling with parallel machines under shift constraints. Serious problems in the long term and improve workshop production efficiency.

为实现上述目的,本发明提出了一种班次约束下带并行机作业车间的排产方法,包括如下步骤:In order to achieve the above purpose, the present invention proposes a production scheduling method with a parallel machine workshop under shift constraints, comprising the following steps:

S1、以最小化拖期为目标构建约束条件下的作业车间调度模型,所述约束条件包括班次约束和工件优先级约束;S1. A job shop scheduling model under constraints is constructed with the goal of minimizing delays, and the constraints include shift constraints and workpiece priority constraints;

所述班次约束为:工件的一道工序只能在选定设备上的一个班次内完成加工,且加工完成时刻不超过所在班次的结束时刻;The shift constraint is: a process of the workpiece can only be processed in one shift on the selected equipment, and the processing completion time does not exceed the end time of the shift;

所述工件优先级约束为:所有工件分为急件组和普通组,急件组不可逾期交付,普通组可逾期交付;The priority constraints of the workpieces are: all workpieces are divided into urgent groups and ordinary groups, the urgent group cannot be overdue for delivery, and the ordinary group can be overdue;

S2、根据带并行机作业车间的生产信息,对作业车间调度模型进行求解,得到车间最佳排产方案。S2. According to the production information of the job shop with parallel machines, the job shop scheduling model is solved to obtain the best production scheduling scheme of the shop.

作为进一步优选的,在作业车间调度模型中,目标函数f为:As a further preference, in the job shop scheduling model, the objective function f is:

Figure BDA0003711687940000021
Figure BDA0003711687940000021

所述班次约束如下:The shift constraints are as follows:

Sij≥Ci(j-1) S ij ≥C i(j-1)

Bkd≥Ek(d-1) B kd ≥E k(d-1)

所述工件优先级约束如下:The artifact priority constraints are as follows:

Figure BDA0003711687940000022
Figure BDA0003711687940000022

其中,n表示工件总数,Ci表示工件i的完工时间,Di表示工件i的交货期;Sij表示工件i的第j道工序Oij的开始加工时间,Ci(j-1)表示工序Oij上一道工序的完工时间;Bkd表示设备k第d个班次的开始时间,Ekd表示设备k第d个班次的结束时间。Among them, n represents the total number of workpieces, C i represents the completion time of workpiece i, D i represents the delivery date of workpiece i; S ij represents the start processing time of the jth operation O ij of workpiece i, and C i(j-1) represents the completion time of the previous process of process O ij ; B kd represents the start time of the d-th shift of equipment k, and E kd represents the end time of the d-th shift of equipment k.

作为进一步优选的,通过蚁群算法对作业车间调度模型进行求解。As a further preference, the job shop scheduling model is solved by an ant colony algorithm.

作为进一步优选的,通过蚁群算法对作业车间调度模型进行求解时,先对急件组进行求解,若急件组解为不可行解,即急件存在拖期,则进行班次替换;然后在急件组解的基础上,对普通组进行求解。As a further preference, when solving the job shop scheduling model through the ant colony algorithm, first solve the urgent group, if the urgent group is an infeasible solution, that is, there is a delay in the urgent, then the shift is replaced; and then the urgent group is solved. On the basis of , solve the common group.

作为进一步优选的,进行班次替换时,采用班次替换期望指数表示班次被期待延长时间的指数高低,据此进行班次替换;班次替换期望指数的计算式为:As a further preference, when the shift replacement is performed, the shift replacement expectation index is used to indicate the level of the index of the expected extended time of the shift, and the shift replacement is performed accordingly; the calculation formula of the shift replacement expectation index is:

Figure BDA0003711687940000031
Figure BDA0003711687940000031

其中,ESkd表示设备k班次d的班次替换期望指数;Etimekdij表示工序Oij不可安排在班次SFTkd上时,班次SFTkd延长多长时间可满足工序Oij的需求;Pstaykdij表示若延长时间满足Oij,该工序留在该班次加工的概率指数。Among them, ES kd represents the shift replacement expectation index of equipment k shift d; Etime kdij represents how long the shift SFT kd can be extended to meet the needs of the process O ij when the process O ij cannot be arranged on the shift SFT kd ; Pstay kdij represents if the process is extended If the time satisfies O ij , the probability index of the operation remaining in the shift processing.

作为进一步优选的,进行班次替换的具体步骤如下:As further preferred, the concrete steps of carrying out shift replacement are as follows:

(1)将每个班次的初始班次替换期望指数设置为0,按照不可行解中的工序顺序依次遍历每个工序;(1) Set the initial shift replacement expectation index of each shift to 0, and traverse each process in turn according to the sequence of processes in the infeasible solution;

(2)对工序Oij,从其可开始加工时刻到当前开工时刻的区间DG内,遍历其当前选定加工设备上的班次,形成待选班次集合GTij(2) for operation O ij , in the interval DG from the moment it can start processing to the current start moment, traverse the shifts on its currently selected processing equipment to form a set of shifts to be selected GT ij ;

(3)对于GTij中的每个班次,计算若Oij在其上加工则需要延长的时间Etimekdij;若延长时间超过下一个班次开始时间,则Etimekdijj设置为0;否则,保留Etimekdij的值,并计算Pstaykdij;进而得到GTij中每个班次的ESkd(3) For each shift in GT ij , calculate the time Etime kdij that needs to be extended if O ij is processed on it; if the extended time exceeds the start time of the next shift, then Etime kdijj is set to 0; otherwise, keep Etime kdij , and calculate Pstay kdij ; then obtain ES kd of each shift in GT ij ;

(4)重复步骤(2)和(3),直至所有工序遍历完成;(4) Repeat steps (2) and (3) until all processes are traversed;

(5)遍历所有ESkd大于0的班次,根据ESkd大小,进行班次替换;ESkd越大,选择的替换班次时长越长;若某设备所有工作日上的班次均已被替换为最长时间班次,仍无法满足交货期,则在需要替换班次所在周的周末安排休息日加班一日;若本周周末排班已满,则往后顺延;(5) Traverse all shifts with ES kd greater than 0, and perform shift replacement according to the size of ES kd ; the larger the ES kd , the longer the selected replacement shift time; if the shifts on all working days of a device have been replaced with the longest shift If the time shift is still unable to meet the delivery time, it will be arranged to work overtime on a rest day on the weekend of the week where the shift needs to be replaced; if the weekend schedule is full, it will be postponed later;

(6)班次替换完成后,若当前解中仍存在拖期,则转回步骤(1)进行下一轮班次替换,直至产生可行解,或超过一定的迭代次数放弃当前解,结束。(6) After the shift replacement is completed, if there is still a delay in the current solution, go back to step (1) for the next round of shift replacement, until a feasible solution is generated, or the current solution is abandoned after a certain number of iterations, and the end.

作为进一步优选的,概率指数Pstaykdij的计算式为:As a further preference, the calculation formula of the probability index Pstay kdij is:

Figure BDA0003711687940000041
Figure BDA0003711687940000041

其中,R为一个大于所有STDij的常数,STDij为从工序Oij可开始加工时刻起的天数。Among them, R is a constant greater than all STD ij , and STD ij is the number of days from the moment when the process O ij can start processing.

作为进一步优选的,获取急件组解或普通组解后,对当前解Vk进行局部搜索,得到解Vk′;分别计算相应的目标函数值,保留Vk′和Vk中较优的一个作为解;As a further preference, after obtaining the emergency solution or the common solution, perform a local search on the current solution V k to obtain the solution V k '; calculate the corresponding objective function values respectively, and keep the better one of V k ' and V k as a solution;

局部搜索的策略为:在每个设备组中,随机选择一台设备,然后随机选择该设备上的一道工序;对该工序,从其所属工件的上一道工序结束时间开始,到该工序完成时间为止,遍历设备组中的其他设备的空闲时间段;若组内其他设备有可容纳该工序且开工时间早于当前设备的时间段,则将该工序调整至最早可开工的设备上;若所有选择出的工序均未被调整加工设备,则随机选择这些工序所在设备组其他设备上的任意一道工序,交换两道工序的加工设备。The strategy of local search is: in each equipment group, randomly select a piece of equipment, and then randomly select a process on the equipment; for this process, start from the end time of the previous process of the workpiece to which it belongs, to the completion time of the process traverse the idle time period of other equipment in the equipment group; if other equipment in the group has a time period that can accommodate the process and the start time is earlier than the current equipment, the process will be adjusted to the earliest equipment that can start; If none of the selected processes has been adjusted to the processing equipment, then randomly select any process on other equipment in the equipment group where these processes are located, and exchange the processing equipment of the two processes.

作为进一步优选的,通过嵌套蚁群算法对作业车间调度模型进行求解,将调度过程分为工序排序和设备选择两个阶段,工序排序阶段的蚁群为大蚁群,设备选择阶段的蚁群为小蚁群;大蚁群中的每只大蚁对工序进行遍历,获得工序排序方案,小蚁群根据每只大蚁获得的方案,对设备进行遍历,获得设备选择方案,从而形成完整解。As a further preference, the job shop scheduling model is solved by the nested ant colony algorithm, and the scheduling process is divided into two stages: process sorting and equipment selection. The ant colony in the process sorting stage is the large ant colony, and the ant colony in the equipment selection stage It is a small ant colony; each large ant in the large ant colony traverses the process to obtain a process sequencing plan, and the small ant colony traverses the equipment according to the plan obtained by each large ant to obtain the equipment selection plan, thus forming a complete solution. .

作为进一步优选的,在工序排序阶段,优先安排加工时长与班次时长之比更大的工序在设备上加工;在设备选择阶段,优先选择同一设备组内开工时间早的设备。As a further preference, in the process sorting stage, the process with a larger ratio of processing time to shift time is preferentially arranged to be processed on the equipment; in the equipment selection stage, the equipment with the earliest start time in the same equipment group is preferentially selected.

总体而言,通过本发明所构思的以上技术方案与现有技术相比,主要具备以下的技术优点:In general, compared with the prior art, the above technical solutions conceived by the present invention mainly have the following technical advantages:

1.本发明在传统带并行机的作业车间调度问题上,考虑班次约束,以最小化拖期为目标建立班次约束下带并行机的作业车间调度问题数学模型;考虑工件优先级问题,按工件优先级将工件分为急件组和普通组,避免急件拖期。本发明可解决工件拖期严重的问题,提高车间生产效率,并且能够合理安排加班,保证急件按时交付。1. On the traditional job shop scheduling problem with parallel machines, the present invention considers shift constraints, and establishes a mathematical model of job shop scheduling problems with parallel machines under shift constraints with the goal of minimizing delays; Priority divides workpieces into urgent groups and common groups to avoid delays in urgent items. The invention can solve the problem of serious delay of workpieces, improve workshop production efficiency, and can reasonably arrange overtime to ensure timely delivery of urgent items.

2.本发明结合工件优先级设计的班次替换规则,可有效进行班次替换,达到尽可能合理安排加班以解决急件拖期问题的效果,大大降低了因急件无法按时交付而产生的违约成本。2. The present invention combines the shift replacement rules designed by the priority of the workpieces, which can effectively carry out shift replacement, achieve the effect of arranging overtime as reasonably as possible to solve the problem of delayed delivery of urgent items, and greatly reduce the cost of breach of contract caused by the failure to deliver urgent items on time.

3.本发明将调度问题分为工序排序和设备选择两个阶段,通过嵌套蚁群算法进行分阶段求解;同时,根据问题特点,设计启发式规则和局部搜索策略,优化求解结果,并提高迭代效率。3. The present invention divides the scheduling problem into two stages: process sequencing and equipment selection, and is solved in stages by nested ant colony algorithm; at the same time, according to the characteristics of the problem, heuristic rules and local search strategies are designed to optimize the solution results and improve the results. Iterative efficiency.

附图说明Description of drawings

图1为本发明实施例局部搜索时第一种工序调整方式示意图;1 is a schematic diagram of a first process adjustment mode during local search according to an embodiment of the present invention;

图2为本发明实施例局部搜索时第二种工序调整方式示意图;2 is a schematic diagram of a second process adjustment mode during local search according to an embodiment of the present invention;

图3为本发明实施例求解方案流程图;3 is a flowchart of a solution solution according to an embodiment of the present invention;

图4是本发明实施例基于蚁群算法求解初始调度方案的流程图;4 is a flowchart of solving an initial scheduling scheme based on an ant colony algorithm according to an embodiment of the present invention;

图5是本发明实施例工序排序阶段蚂蚁游历路径图;5 is a diagram of an ant travel path in a process sorting stage according to an embodiment of the present invention;

图6是本发明实施例设备选择阶段蚂蚁游历路径图;6 is a diagram of an ant travel path in a device selection stage according to an embodiment of the present invention;

图7是本发明实施例班次替换规则流程图。FIG. 7 is a flowchart of a shift replacement rule according to an embodiment of the present invention.

具体实施方式Detailed ways

为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。In order to make the objectives, technical solutions and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention, but not to limit the present invention. In addition, the technical features involved in the various embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

本发明实施例提供的一种班次约束下带并行机作业车间的排产方法,该车间的基本情况为:设备按类型分组,每组设备中有多台设备,同一设备组内,每台设备加工同种工件同一工序的时长相同;同一设备每日的班次开始、结束时间不一定相同。此外,班次是可变的,一台机器有多种班次,可根据实际排产情况,进行班次替换,包括将班次替换为时长更长的,以安排加班;将时长过长、空闲时间过多的班次替换为时长较短的,以减少工人不必要的工作时间。The embodiment of the present invention provides a production scheduling method for a workshop with parallel machines under shift constraints. The basic situation of the workshop is as follows: equipment is grouped by type, there are multiple pieces of equipment in each group of equipment, and within the same equipment group, each piece of equipment The duration of the same process for processing the same workpiece is the same; the start and end times of the daily shift of the same equipment are not necessarily the same. In addition, the shifts are variable. A machine has multiple shifts, which can be replaced according to the actual production schedule, including replacing the shifts with longer ones to arrange overtime; Replaced shifts with shorter durations to reduce unnecessary working hours for workers.

排产方法包括如下步骤:The scheduling method includes the following steps:

S1、获取带并行机作业车间的生产信息,包括带调度工件的加工时间、工件优先级、加工设备和班次时间等数据;根据工件优先级将工件分为急件组和普通组。S1. Obtain the production information of the workshop with parallel machines, including the processing time of the workpiece with scheduling, the priority of the workpiece, the processing equipment and the shift time and other data; according to the priority of the workpiece, the workpieces are divided into urgent groups and ordinary groups.

S2、考虑工件优先级,以最小化拖期为目标构建班次约束下带并行机的作业车间调度模型。该作业车间调度模型具体如下:S2. Considering the priority of the workpiece, a job shop scheduling model with parallel machines under shift constraints is constructed with the goal of minimizing the delay. The job shop scheduling model is as follows:

参数符号:Parameter notation:

n:工件总数n: total number of workpieces

m:设备组总数m: total number of device groups

h:设备总数h: total number of devices

σ:排产周期最大天数σ: the maximum number of days in the production scheduling cycle

i:工件序号,i={1,2,…,n}i: workpiece serial number, i={1,2,…,n}

j:工件的工序序号,j={1,2,…,Ji}j: the process number of the workpiece, j={1,2,…,J i }

w:设备组编号,w={1,2,…,m}w: device group number, w={1,2,…,m}

k:设备编号,k={1,2,…,h}k: device number, k={1,2,…,h}

d:班次序号,d={1,2,…,σ}d: class number, d={1,2,…,σ}

Oij:工件i的第j道工序O ij : the jth operation of workpiece i

Di:工件i的交货期D i : delivery date of workpiece i

Mij:可加工Oij的设备组M ij : equipment group that can process O ij

tij:工件i的工序j的加工时间t ij : machining time of process j of workpiece i

SFTkd:设备k第d个班次SFT kd : the d-th shift of equipment k

Bkd:设备k第d个班次的开始时间B kd : the start time of the d-th shift of equipment k

Ekd:设备k第d个班次的结束时间E kd : the end time of the d-th shift of equipment k

BKPw:设备组w的备选班次集合BKP w : candidate shift set for equipment group w

Tw:设备组w的班次时长Tw : shift duration of equipment group w

Sij:Oij的开始加工时间S ij : start processing time of O ij

Cij:Oij的完工时间C ij : Completion time of O ij

Ci:工件i的完工时间C i : Completion time of workpiece i

L:一个非常大的正数

Figure BDA0003711687940000071
L: a very large positive number
Figure BDA0003711687940000071

决策变量:Decision variables:

Figure BDA0003711687940000072
Figure BDA0003711687940000072

Figure BDA0003711687940000073
Figure BDA0003711687940000073

Figure BDA0003711687940000074
Figure BDA0003711687940000074

目标函数:Objective function:

Figure BDA0003711687940000075
Figure BDA0003711687940000075

约束条件:Restrictions:

Figure BDA0003711687940000076
Figure BDA0003711687940000076

Figure BDA0003711687940000077
Figure BDA0003711687940000077

Figure BDA0003711687940000078
Figure BDA0003711687940000078

Figure BDA0003711687940000081
Figure BDA0003711687940000081

Figure BDA0003711687940000082
Figure BDA0003711687940000082

Figure BDA0003711687940000083
Figure BDA0003711687940000083

Figure BDA0003711687940000084
Figure BDA0003711687940000084

Figure BDA0003711687940000085
Figure BDA0003711687940000085

Figure BDA0003711687940000086
Figure BDA0003711687940000086

Figure BDA0003711687940000087
Figure BDA0003711687940000087

Figure BDA0003711687940000088
Figure BDA0003711687940000088

其中,式(1)表示目标函数,为了避免加班过度问题,对于急件组的评估指标为提前完工时长,普通组评估指标为拖期;式(2)保证了每道工序只能分配给一个设备,并只能在该设备的一个班次上加工;式(3)给出了完工时间的定义;式(4)保证了工序Oij在未选择的设备上,开工时间和完工时间为0;式(5)规定了一个工件的上一道工序加工完成后才可加工下一道;式(6)表示上一个班次结束,下一个班次才可以开始,班次时段不可重叠;式(7)、式(8)将工序Oij的加工时段限制在选定的加工班次内,即开工时间不能在班次开始时间之前,完工时间不能在班次结束时间之后;式(9)、式(10)保证了两道工序不可同时在同一设备上加工;式(11)、式(12)限制了决策变量的取值范围。Among them, Equation (1) represents the objective function. In order to avoid the problem of excessive overtime, the evaluation index for the urgent group is the early completion time, and the evaluation index for the ordinary group is the delay; Equation (2) ensures that each process can only be assigned to one equipment , and can only be processed on one shift of the equipment; Equation (3) gives the definition of the completion time; Equation (4) ensures that the operation O ij is on the unselected equipment, and the start time and completion time are 0; (5) It is stipulated that the next process can be processed only after the previous process of a workpiece is completed; Equation (6) indicates that the previous shift can be completed before the next shift can begin, and the shift periods cannot overlap; Equation (7), Equation (8) ) limits the processing period of the operation O ij to the selected processing shift, that is, the start time cannot be before the start time of the shift, and the completion time cannot be after the end time of the shift; Equations (9) and (10) guarantee two procedures It cannot be processed on the same equipment at the same time; Equations (11) and (12) limit the value range of decision variables.

S3、根据带并行机作业车间的生产信息,对作业车间调度模型进行求解,调度过程中,先对急件组排产,若存在急件拖期则进行班次替换,再对普通组进行排产,如图3所示。S3. Solve the job shop scheduling model according to the production information of the job shop with parallel machines. During the scheduling process, first schedule the emergency group, if there is an urgent delay, replace the shift, and then schedule the ordinary group, such as shown in Figure 3.

本实施例优选采用蚁群算法进行求解,在每次迭代中,先对急件组求解,在急件组解的基础上再对普通组求解,并根据此次解对信息素进行更新,完成一次迭代;然后进行下一次迭代,直至达到终止条件,输出此时解为最终解,即最终排产方案。具体包括如下步骤:In this embodiment, the ant colony algorithm is preferably used to solve the problem. In each iteration, the emergency group is solved first, and then the common group is solved on the basis of the emergency group solution, and the pheromone is updated according to the solution to complete one iteration. ; and then proceed to the next iteration until the termination condition is reached, and the output at this time is the final solution, that is, the final scheduling plan. Specifically include the following steps:

S31、利用嵌套蚁群算法对急件组求解,如图4所示;S31, using the nested ant colony algorithm to solve the emergency group, as shown in Figure 4;

(1)具体步骤为:(1) The specific steps are:

Step 1工序排序阶段的初始化。输入调度信息,根据各工件的工序信息建立大蚁群遍历路径网,初始化其信息素表,使得每条路径上信息素浓度相等;将M只大蚁分散在任意一可选节点。Step 1 Initialization of the sequence sequence stage. Input the scheduling information, establish a large ant colony traversing path network according to the process information of each workpiece, initialize its pheromone table, so that the pheromone concentration on each path is equal; M large ants are scattered in any optional node.

Step 2工序排序阶段的路径选择。利用本文设计的启发式规则和路径上的信息素引导大蚁遍历每道工序。Step 2: Path selection in the process sorting stage. The heuristic rules designed in this paper and the pheromone on the path are used to guide the ants to traverse each process.

Step 3设备选择阶段的初始化。当某只大蚁ANT遍历结束后,根据大蚁ANT获得的工序排序,建立小蚁群遍历路径网,初始化其信息素表,使得每条路径上信息素浓度相等,并将m只小蚁分散在任意一个可选的设备节点。Step 3 Initialization of the device selection phase. When the traversal of a large ant ANT is completed, according to the process sequence obtained by the large ant ANT, a small ant colony traversal path network is established, and its pheromone table is initialized so that the pheromone concentration on each path is equal, and m small ants are dispersed. on any optional device node.

Step 4设备选择阶段的路径选择。利用本文设计的启发式规则和路径上的信息素,引导每只小蚁ant依次遍历每道工序对应设备组中的设备。Step 4. Path selection in the device selection stage. Using the heuristic rules designed in this paper and the pheromone on the path, each ant is guided to traverse the equipment in the equipment group corresponding to each process in turn.

Step 5设备选择阶段的信息素释放量计算。当某只小蚁ant遍历结束后,则可得到一个完整解,对该解进行评估,计算拖期值。再对该解进行局部搜索,对比之前获得解的拖期值,保留较优的一个,并据此计算信息素释放量。Step 5: Calculate the amount of pheromone released in the equipment selection stage. When the traversal of a small ant is completed, a complete solution can be obtained, the solution can be evaluated, and the delay value can be calculated. Then perform a local search on the solution, compare the delay value of the previously obtained solution, keep the better one, and calculate the pheromone release amount accordingly.

Step 6设备选择阶段的信息素更新。若所有小蚁到达终点,则对小蚁群遍历路径网上的每条路径进行信息素的蒸发和释放;否则,转Step 4。Step 6: Pheromone update in device selection stage. If all the small ants reach the end point, evaporate and release pheromone for each path on the small ant colony's traversed path network; otherwise, go to Step 4.

Step 7工序排序阶段的信息素释放量计算。若小蚁群搜寻达到终止条件,根据小蚁群获得的最优解,决定大蚁ANT释放的信息素量;否则,转Step 4。Step 7: Calculate the amount of pheromone released in the process sequencing stage. If the small ant colony search reaches the termination condition, according to the optimal solution obtained by the small ant colony, determine the amount of pheromone released by the large ant ANT; otherwise, go to Step 4.

Step 8工序排序阶段的信息素更新。若所有大蚁到达终点,则对大蚁群遍历路径上的每条路径进行信息素的蒸发和释放;否则,转Step 2。Step 8: The pheromone update in the process sequencing stage. If all the giant ants reach the end point, evaporate and release pheromone for each path on the traversed path of the giant ant colony; otherwise, go to Step 2.

Step 9判断大蚁群的搜索是否到达终止条件,如果到达则输出最优解;否则,转Step 2。Step 9: Determine whether the search of the large ant colony has reached the termination condition, and if so, output the optimal solution; otherwise, go to Step 2.

(2)算法设计(2) Algorithm design

根据问题的特点,将调度过程分为工序排序和设备选择两个阶段,两个阶段分别用蚁群算法求解,组成嵌套蚁群算法。将工序排序阶段的蚁群称为大蚁群,将设备选择阶段的蚁群称为小蚁群;大蚁群中的每只大蚁对工序进行遍历,获得工序排序方案;小蚁群根据每只大蚁获得的方案,对设备进行遍历,获得最优的设备选择方案,形成完整解并以此更新大蚁遍历路径上的信息素,由此使得两个阶段的算法形成嵌套关系。According to the characteristics of the problem, the scheduling process is divided into two stages: process sequencing and equipment selection. The two stages are solved by ant colony algorithm respectively to form a nested ant colony algorithm. The ant colony in the process sorting stage is called the large ant colony, and the ant colony in the equipment selection stage is called the small ant colony; each large ant in the large ant colony traverses the process to obtain the process sorting plan; Only the solution obtained by the big ants, traverse the equipment, obtain the optimal equipment selection scheme, form a complete solution and update the pheromone on the traversal path of the big ants, thus making the two-stage algorithm form a nested relationship.

1)工序排序阶段1) Process sequencing stage

本阶段仅确定工件的加工先后顺序,不考虑设备的选择,因此可看做一个JSP问题。算法过程可描述为,令每一只蚂蚁在同一工件的工序先后顺序约束下,遍历所有工序。因此可把一道工序看做蚂蚁游历图中的一个节点,以3×2的JSP问题为例,如图所示,节点0和节点7为虚节点,表示蚂蚁遍历时的开始和结束节点,所有蚂蚁从节点0出发,遍历所有工序节点后到节点7结束。首先,将所有工件的下一道待加工工序添加至下一可遍历节点集合中;当遍历过一个节点后,将其从可选节点集合中删除;若该节点代表工序所属工件并未加工完成,则添加其下一道工序的节点至可选节点集合,直至所有节点遍历完成。当遍历结束,蚂蚁走过的节点序列即本次循环得到的一个工序排序。At this stage, only the processing sequence of the workpiece is determined, and the selection of equipment is not considered, so it can be regarded as a JSP problem. The algorithm process can be described as, let each ant traverse all the processes under the constraints of the process sequence of the same workpiece. Therefore, a process can be regarded as a node in the ant traversal graph. Take the 3×2 JSP problem as an example. As shown in the figure, node 0 and node 7 are virtual nodes, indicating the start and end nodes of the ant traversal. All The ant starts from node 0, traverses all process nodes, and ends at node 7. First, add the next process to be processed of all workpieces to the next traversable node set; after traversing a node, delete it from the optional node set; if the node represents the workpiece to which the process belongs and has not been processed, Then add the node of its next process to the optional node set until all nodes are traversed. When the traversal ends, the sequence of nodes traversed by the ants is a sequence of processes obtained in this cycle.

蚂蚁根据如下所示的伪随机概率选择下一个遍历节点:The ant chooses the next traversed node according to a pseudo-random probability as shown below:

Figure BDA0003711687940000101
Figure BDA0003711687940000101

式中,

Figure BDA0003711687940000102
为第k只蚂蚁在节点i时选择节点j的概率,τij(t)为节点i到j的路径上信息素残余量,ηij(t)为启发式规则,α为残余信息素的重要程度,β为启发式规则的重要程度。In the formula,
Figure BDA0003711687940000102
is the probability that the kth ant selects node j at node i, τ ij (t) is the residual amount of pheromone on the path from node i to j, η ij (t) is the heuristic rule, and α is the importance of the residual pheromone degree, β is the importance degree of the heuristic rule.

选择工序Oij时的启发式规则为:The heuristic rules for selecting the operation O ij are:

Figure BDA0003711687940000111
Figure BDA0003711687940000111

式中,

Figure BDA0003711687940000112
——工序Oij加工设备组的班次时长。In the formula,
Figure BDA0003711687940000112
——The shift duration of the processing equipment group in the process O ij .

该启发式规则决定了,工序的加工时间越接近班次时长、交货期越靠前时,该工序被优先选择的概率越大。优先选择加工时间与班次时长之比更接近的工序,是由于在同一班次中,短工序若在长工序前加工,有很大的概率使得,长工序在可开始加工时刻之后的很长一段时间内,找不到大于长工序加工时长的时间段,使得拖期愈加严重。The heuristic rule determines that the closer the processing time of an operation is to the shift duration and the earlier the delivery date, the higher the probability of this operation being selected first. The process with a closer ratio of processing time to shift duration is preferred because in the same shift, if the short process is processed before the long process, there is a high probability that the long process will be processed for a long period of time after the start of processing. In this case, a time period longer than the processing time of a long process cannot be found, which makes the delay more serious.

2)设备选择阶段2) Equipment selection stage

根据上一阶段获取的工件序列,按顺序确定每个工序的加工设备。作为示例,如图所示,图中的节点表示设备,同一列的节点表示同一个设备组的几台设备。According to the workpiece sequence obtained in the previous stage, the processing equipment of each process is determined in sequence. As an example, as shown in the figure, the nodes in the figure represent devices, and the nodes in the same column represent several devices in the same device group.

设备按照设备组顺序进行编号,如表1所示:Devices are numbered in the order of device groups, as shown in Table 1:

表1设备编号Table 1 Device No.

Figure BDA0003711687940000113
Figure BDA0003711687940000113

该阶段算法基本思想是,每只蚂蚁按照确定工序顺序,依次确定每道工序Oij的加工设备。每道工序有一个对应的设备组Mij,蚂蚁对应设备组中的每个设备进行遍历,寻找出每台设备可最早开工的时间段;然后根据伪随机概率公式计算出来的概率,随机选择设备组中的一台设备,将工序Oij安排在该设备上,即在图中每列节点中选择一个;接着,遍历下一道工序的设备节点;如此循环往复,直至为所有工序都确定加工设备。The basic idea of the algorithm in this stage is that each ant determines the processing equipment of each process O ij in turn according to the determined process sequence. Each process has a corresponding equipment group M ij , and the ants traverse each equipment in the equipment group to find the earliest time period when each equipment can be started; and then randomly select the equipment according to the probability calculated by the pseudo-random probability formula One device in the group, arrange the process O ij on this device, that is, select one node in each column of the graph; then, traverse the device nodes of the next process; and so on, until the processing equipment is determined for all processes .

在本阶段,加工设备选择的启发式规则为:At this stage, the heuristic rules for processing equipment selection are:

Figure BDA0003711687940000121
Figure BDA0003711687940000121

该启发式规则决定了优先选择可开工早的设备。This heuristic rule determines the preference for devices that can start early.

S32、在蚂蚁遍历结束后,对急件组获取的当前解Vk进行局部搜索,得到解Vk′。局部搜索得到解Vk′后,分别计算Vk′、Vk对应的目标函数值,保存Vk′和Vk中较优的一个,并更新信息素。S32. After the ant traversal ends, perform a local search on the current solution V k obtained by the emergency group to obtain the solution V k ′. After the solution V k ′ is obtained by local search, the objective function values corresponding to V k ′ and V k are calculated respectively, the better one of V k ′ and V k is saved, and the pheromone is updated.

局部搜索策略为:在每个设备组中,随机选择一台设备,然后随机选择该设备上的一道工序。对该工序,从其所属工件的上一道工序结束时间开始,到该工序完成时间为止,遍历设备组中的其他设备的空闲时间段。The local search strategy is: in each equipment group, randomly select a piece of equipment, and then randomly select a process on this equipment. For the process, from the end time of the previous process of the workpiece to which it belongs, to the completion time of the process, traverse the idle time period of other equipment in the equipment group.

若组内其他设备有可容纳该工序且开工时间早于当前设备的时间段,则将该工序调整至最早可开工的设备上。作为示例,如图1所示,对于设备组2,随机选择设备M22上的工序O22,遍历设备M21和设备M23上从第一个班次时刻2到时刻8上的空闲时间段;发现设备M21上从时刻6到时刻8可安排工序O22,则将工序O22调整至设备M21的这一时间段上加工。If other equipment in the group has a time period that can accommodate the process and the start time is earlier than the current equipment, the process will be adjusted to the earliest equipment that can start work. As an example, as shown in FIG. 1, for equipment group 2, process O 22 on equipment M 22 is randomly selected, and the idle time period from time 2 to time 8 of the first shift on equipment M 21 and M 23 is traversed; It is found that the process O 22 can be arranged on the equipment M 21 from time 6 to time 8, then the process O 22 is adjusted to be processed in this time period of the equipment M 21 .

若所有选择出的工序均未被调整加工设备,则随机选择这些工序所在设备组其他设备上的任意一道工序,交换两道工序的加工设备,如图2所示。If all the selected processes have not adjusted the processing equipment, randomly select any process on other equipment in the equipment group where these processes are located, and exchange the processing equipment of the two processes, as shown in Figure 2.

S33、若当前的急件组解为不可行解,则还需进行班次替换,通过安排加班将完工时间提前,从而将不可行解转化为可行解,以解决急件拖期问题。如图7所示,班次替换规则如下:S33. If the current urgent group solution is an infeasible solution, it is necessary to perform shift replacement, and advance the completion time by arranging overtime, so as to convert the infeasible solution into a feasible solution to solve the problem of urgent delay. As shown in Figure 7, the shift replacement rules are as follows:

Step 1初始化。用班次替换期望指数ES的概念表示班次被期待延长时间的指数高低,将每个班次的ES初始化设置为0,按照不可行解中的工序顺序依次遍历每个工序。Step 1Initialization. The concept of replacing the expected index ES with the shift represents the index of the expected extended time of the shift. The ES of each shift is initialized to 0, and each process is traversed in turn according to the sequence of the processes in the infeasible solution.

Step 2对工序Oij,从其可开始加工时刻到当前开工时刻的区间DG内,遍历其当前选定加工设备上的班次,形成待选班次集合GTijStep 2: For the process O ij , traverse the shifts on the currently selected processing equipment within the interval DG from the time when it can start processing to the current start time, and form a set of shifts to be selected GT ij .

Step 3对于GTij中的每个班次,计算若Oij在其上加工则需要延长的时间Etimekdij。若延长时间超过下一个班次开始时间,则Etimekdijj设置为0;否则,计算Pstaykdij,公式为:Step 3: For each shift in GT ij , calculate the extended time Etime kdij if O ij is processed on it. If the extended time exceeds the start time of the next shift, Etime kdijj is set to 0; otherwise, Pstay kdij is calculated with the formula:

Figure BDA0003711687940000131
Figure BDA0003711687940000131

式中,R——常数,STDij——工序Oij从可开始加工时刻开始的第STDij天。该式表示,当班次越靠前,班次替换后工序Oij选择的概率越高。In the formula, R——constant, STD ij ——the STD ij day of the process O ij from the time when the processing can be started. This formula indicates that the higher the shift is, the higher the probability of selecting the process O ij after the shift is replaced.

Step 4计算GTij中每个班次的班次替换期望指数ESkd。若所有工序遍历完成,则转Step 5;否则,转Step 2遍历下一道工序。ESkd的计算式为:Step 4 Calculate the shift replacement expectation index ES kd for each shift in GT ij . If all processes are traversed, go to Step 5; otherwise, go to Step 2 to traverse the next process. The formula for calculating ES kd is:

Figure BDA0003711687940000132
Figure BDA0003711687940000132

其中,ESkd表示设备k班次d的替换期望指数;Etimekdij表示工序Oij不可安排在班次SFTkd上时,班次SFTkd延长多长时间可满足工序Oij的需求;Pstaykdij表示若延长时间满足Oij,则该工序留在该班次加工的概率指数。Among them, ES kd represents the replacement expectation index of equipment k shift d; Etime kdij represents how long the shift SFT kd can be extended to meet the needs of the process O ij when the process O ij cannot be arranged on the shift SFT kd ; Pstay kdij represents if the time is extended Satisfying O ij , the probability index of the operation remaining in the shift processing.

Step 5遍历所有期望指数ESkd大于0的班次,根据ESkd大小,进行班次替换。ESkd越大,选择的替换班次时长越长。若某设备所有工作日上的班次均已被替换为最长时间班次,仍无法满足交货期,则在需要替换班次所在周的周末安排休息日加班一日;若本周周末排班已满,则往后顺延。Step 5 Traverse all the shifts with the expected index ES kd greater than 0, and perform shift replacement according to the size of ES kd . The larger the ES kd , the longer the replacement shift duration is selected. If the shifts on all working days of a certain equipment have been replaced with the longest shifts, and the delivery time is still not met, then a rest day will be arranged on the weekend of the week where the shift needs to be replaced; , it will be postponed later.

Step 6班次替换完成后,将当前解重新进行主动式解码,若存在拖期,则转Step 1进行下一轮班次替换,直至产生可行解,或超过一定的迭代次数放弃当前解,结束。Step 6 After the shift replacement is completed, the current solution is re-actively decoded. If there is a delay, go to Step 1 to perform the next round of shift replacement until a feasible solution is generated, or the current solution is abandoned after a certain number of iterations, and the end.

S34、在急件组解的基础上,进行普通组求解,求解方法与急件组相同;具体的,对普通组求解时,利用嵌套蚁群算法求解并进行局部搜索,可不进行班次替换。S34. On the basis of the urgent group solution, perform the ordinary group solution, and the solution method is the same as that of the urgent group; specifically, when solving the ordinary group, use the nested ant colony algorithm to solve and perform a local search, and shift replacement may not be performed.

S35、剔除不必要的班次替换。对于急件组解和普通组解构成的整体解,若发现替换前班次可满足当前班次生产任务的,则消除此类不必要的替换。S35. Eliminate unnecessary shift replacements. For the overall solution composed of urgent group solution and common group solution, if it is found that the pre-replacement shift can meet the production task of the current shift, such unnecessary replacement will be eliminated.

具体的,对急件组求解时进行了班次替换,其中会产生一些不必要的班次替换(剔除后不影响急件组),但其中一部分班次替换后可以满足普通组的生产,即剔除后会对普通组的生产产生影响,所以这里只剔除对急件组和普通组都不产生影响的班次替换。Specifically, when the urgent group is solved, shift replacement is carried out, and some unnecessary shift replacements will occur (the urgent group will not be affected after removal). The production of the group has an impact, so only shift replacements that have no impact on the urgent group and the ordinary group are excluded here.

S36、评估当前解。此时完成一次迭代,根据当前解计算急件组工件提前的总时长和普通组拖期总时长,以此更新蚂蚁遍历网上的信息素。S36. Evaluate the current solution. At this point, one iteration is completed, and according to the current solution, the total advance time of the emergency group workpiece and the total delay time of the common group are calculated, so as to update the pheromone on the ant traversal network.

S37、若未达到终止条件(即预设的迭代次数),则回到步骤S31继续下一轮迭代;否则,以当前解为最优解,得到车间最佳排产方案。S37. If the termination condition (ie the preset number of iterations) is not reached, go back to step S31 to continue the next round of iterations; otherwise, take the current solution as the optimal solution to obtain the optimal production scheduling plan for the workshop.

必须说明的是,上述实施例中,方法并不必然按照序号顺序依次执行,只要从执行逻辑中不能推定必然按某一顺序执行,则意味着可以以其他任何可能的顺序执行。It must be noted that, in the above-mentioned embodiments, the methods are not necessarily executed sequentially according to the sequence number, as long as it cannot be inferred from the execution logic that the methods must be executed in a certain sequence, it means that the methods can be executed in any other possible sequence.

综上,在问题层面上,本发明在传统带并行机的作业车间调度问题上,考虑班次约束,以最小化拖期为目标建立班次约束下带并行机的作业车间调度问题数学模型;考虑工件优先级问题,按工件优先级将工件分为急件组和普通组,通过进行合理的班次替换解决急件拖期问题。在算法层面上,将问题分为工序排序和设备选择两个阶段,设计嵌套蚁群进行分阶段求解;根据问题特点,设计启发式规则和局部搜索策略;设计班次替换规则以解决急件拖期问题。本发明可根据调度目标生成较优的调度方案,并且能够合理安排加班,保证急件按时交付。To sum up, at the problem level, the present invention considers the shift constraints in the traditional job shop scheduling problem with parallel machines, and aims to minimize the delay time to establish a mathematical model of the job shop scheduling problem with parallel machines under shift constraints; considering the workpiece Priority problem, divide the workpieces into urgent groups and ordinary groups according to the priority of the workpieces, and solve the problem of delaying urgent items by replacing them with reasonable shifts. At the algorithm level, the problem is divided into two stages: process sequencing and equipment selection, and nested ant colonies are designed to solve them in stages; according to the characteristics of the problem, heuristic rules and local search strategies are designed; shift replacement rules are designed to solve urgent delays question. The present invention can generate an optimal scheduling scheme according to the scheduling target, and can reasonably arrange overtime to ensure timely delivery of urgent items.

本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。Those skilled in the art can easily understand that the above are only preferred embodiments of the present invention, and are not intended to limit the present invention. Any modifications, equivalent replacements and improvements made within the spirit and principles of the present invention, etc., All should be included within the protection scope of the present invention.

Claims (10)

1.一种班次约束下带并行机作业车间的排产方法,其特征在于,包括如下步骤:1. a production scheduling method with parallel machine workshop under a shift constraint, is characterized in that, comprises the steps: S1、以最小化拖期为目标构建约束条件下的作业车间调度模型,所述约束条件包括班次约束和工件优先级约束;S1. A job shop scheduling model under constraints is constructed with the goal of minimizing delays, and the constraints include shift constraints and workpiece priority constraints; 所述班次约束为:工件的一道工序只能在选定设备上的一个班次内完成加工,且加工完成时刻不超过所在班次的结束时刻;The shift constraint is: a process of the workpiece can only be processed in one shift on the selected equipment, and the processing completion time does not exceed the end time of the shift; 所述工件优先级约束为:所有工件分为急件组和普通组,急件组不可逾期交付,普通组可逾期交付;The priority constraints of the workpieces are: all workpieces are divided into urgent groups and ordinary groups, the urgent group cannot be overdue for delivery, and the ordinary group can be overdue; S2、根据带并行机作业车间的生产信息,对作业车间调度模型进行求解,得到车间最佳排产方案。S2. According to the production information of the job shop with parallel machines, the job shop scheduling model is solved to obtain the best production scheduling scheme of the shop. 2.如权利要求1所述的班次约束下带并行机作业车间的排产方法,其特征在于,在作业车间调度模型中,目标函数f为:2. the production scheduling method with parallel machine job shop under shift constraint as claimed in claim 1, it is characterized in that, in job shop scheduling model, objective function f is:
Figure FDA0003711687930000011
Figure FDA0003711687930000011
所述班次约束如下:The shift constraints are as follows: Sij≥Ci(j-1) S ij ≥C i(j-1) Bkd≥Ek(d-1)B kd ≥E k ( d-1 ) 所述工件优先级约束如下:The artifact priority constraints are as follows:
Figure FDA0003711687930000012
Figure FDA0003711687930000012
其中,n表示工件总数,Ci表示工件i的完工时间,Di表示工件i的交货期;Sij表示工件i的第j道工序Oij的开始加工时间,Ci(j-1)表示工序Oij上一道工序的完工时间;Bkd表示设备k第d个班次的开始时间,Ekd表示设备k第d个班次的结束时间。Among them, n represents the total number of workpieces, C i represents the completion time of workpiece i, D i represents the delivery date of workpiece i; S ij represents the start processing time of the jth operation O ij of workpiece i, and C i(j-1) represents the completion time of the previous process of process O ij ; B kd represents the start time of the d-th shift of equipment k, and E kd represents the end time of the d-th shift of equipment k.
3.如权利要求2所述的班次约束下带并行机作业车间的排产方法,其特征在于,通过蚁群算法对作业车间调度模型进行求解。3 . The production scheduling method for a job shop with parallel machines under shift constraints according to claim 2 , wherein the job shop scheduling model is solved by an ant colony algorithm. 4 . 4.如权利要求3所述的班次约束下带并行机作业车间的排产方法,其特征在于,通过蚁群算法对作业车间调度模型进行求解时,先对急件组进行求解,若急件组解为不可行解,即急件存在拖期,则进行班次替换;然后在急件组解的基础上,对普通组进行求解。4. The production scheduling method with parallel machine job shop under shift constraint as claimed in claim 3, it is characterized in that, when the job shop scheduling model is solved by ant colony algorithm, the urgent group is solved first, if the urgent group is solved If the solution is infeasible, that is, there is a delay in the emergency, then the shift is replaced; then, on the basis of the emergency group solution, the common group is solved. 5.如权利要求4所述的班次约束下带并行机作业车间的排产方法,其特征在于,进行班次替换时,采用班次替换期望指数表示班次被期待延长时间的指数高低,据此进行班次替换;班次替换期望指数的计算式为:5. The production scheduling method with parallel machine workshops under shift constraints as claimed in claim 4, characterized in that, when performing shift replacement, a shift replacement expectation index is used to represent the index level of the expected extended time of the shift, and the shift is performed accordingly. Replacement; the formula for the shift replacement expectation index is:
Figure FDA0003711687930000021
Figure FDA0003711687930000021
其中,ESkd表示设备k班次d的班次替换期望指数;Etimekdij表示工序Oij不可安排在班次SFTkd上时,班次SFTkd延长多长时间可满足工序Oij的需求;Pstaykdij表示若延长时间满足Oij,该工序留在该班次加工的概率指数。Among them, ES kd represents the shift replacement expectation index of equipment k shift d; Etime kdij represents how long the shift SFT kd can be extended to meet the needs of the process O ij when the process O ij cannot be arranged on the shift SFT kd ; Pstay kdij represents if the process is extended If the time satisfies O ij , the probability index of the operation remaining in the shift processing.
6.如权利要求5所述的班次约束下带并行机作业车间的排产方法,其特征在于,进行班次替换的具体步骤如下:6. the production scheduling method with parallel machine workshop under shift constraint as claimed in claim 5, is characterized in that, the concrete steps of carrying out shift replacement are as follows: (1)将每个班次的初始班次替换期望指数设置为0,按照不可行解中的工序顺序依次遍历每个工序;(1) Set the initial shift replacement expectation index of each shift to 0, and traverse each process in turn according to the sequence of processes in the infeasible solution; (2)对工序Oij,从其可开始加工时刻到当前开工时刻的区间DG内,遍历其当前选定加工设备上的班次,形成待选班次集合GTij(2) for operation O ij , in the interval DG from the moment it can start processing to the current start moment, traverse the shifts on its currently selected processing equipment to form a set of shifts to be selected GT ij ; (3)对于GTij中的每个班次,计算若Oij在其上加工则需要延长的时间Etimekdij;若延长时间超过下一个班次开始时间,则Etimekdijj设置为0;否则,保留Etimekdij的值,并计算Pstaykdij;进而得到GTij中每个班次的ESkd(3) For each shift in GT ij , calculate the time Etime kdij that needs to be extended if O ij is processed on it; if the extended time exceeds the start time of the next shift, then Etime kdijj is set to 0; otherwise, keep Etime kdij , and calculate Pstay kdij ; then obtain ES kd of each shift in GT ij ; (4)重复步骤(2)和(3),直至所有工序遍历完成;(4) Repeat steps (2) and (3) until all processes are traversed; (5)遍历所有ESkd大于0的班次,根据ESkd大小,进行班次替换;ESkd越大,选择的替换班次时长越长;若某设备所有工作日上的班次均已被替换为最长时间班次,仍无法满足交货期,则在需要替换班次所在周的周末安排休息日加班一日;若本周周末排班已满,则往后顺延;(5) Traverse all shifts with ES kd greater than 0, and perform shift replacement according to the size of ES kd ; the larger the ES kd , the longer the selected replacement shift time; if the shifts on all working days of a device have been replaced with the longest shift If the time shift is still unable to meet the delivery time, it will be arranged to work overtime on a rest day on the weekend of the week where the shift needs to be replaced; if the weekend schedule is full, it will be postponed later; (6)班次替换完成后,若当前解中仍存在拖期,则转回步骤(1)进行下一轮班次替换,直至产生可行解,或超过一定的迭代次数放弃当前解,结束。(6) After the shift replacement is completed, if there is still a delay in the current solution, go back to step (1) for the next round of shift replacement, until a feasible solution is generated, or the current solution is abandoned after a certain number of iterations, and the end. 7.如权利要求6所述的班次约束下带并行机作业车间的排产方法,其特征在于,概率指数Pstaykdij的计算式为:7. the production scheduling method with parallel machine workshop under the shift constraint as claimed in claim 6, is characterized in that, the calculation formula of probability index Pstay kdij is:
Figure FDA0003711687930000031
Figure FDA0003711687930000031
其中,R为一个大于所有STDij的常数,STDij为从工序Oij可开始加工时刻起的天数。Among them, R is a constant greater than all STD ij , and STD ij is the number of days from the moment when the process O ij can start processing.
8.如权利要求4所述的班次约束下带并行机作业车间的排产方法,其特征在于,获取急件组解或普通组解后,对当前解Vk进行局部搜索,得到解Vk′;分别计算相应的目标函数值,保留Vk′和Vk中较优的一个作为解;8. The production scheduling method with a parallel machine job shop under shift constraints as claimed in claim 4, characterized in that, after obtaining the emergency solution or the common solution, a local search is performed on the current solution V k to obtain the solution V k ′ ; Calculate the corresponding objective function values respectively, and keep the better one of V k ′ and V k as the solution; 局部搜索的策略为:在每个设备组中,随机选择一台设备,然后随机选择该设备上的一道工序;对该工序,从其所属工件的上一道工序结束时间开始,到该工序完成时间为止,遍历设备组中的其他设备的空闲时间段;若组内其他设备有可容纳该工序且开工时间早于当前设备的时间段,则将该工序调整至最早可开工的设备上;若所有选择出的工序均未被调整加工设备,则随机选择这些工序所在设备组其他设备上的任意一道工序,交换两道工序的加工设备。The strategy of local search is: in each equipment group, randomly select a piece of equipment, and then randomly select a process on the equipment; for this process, start from the end time of the previous process of the workpiece to which it belongs, to the completion time of the process traverse the idle time period of other equipment in the equipment group; if other equipment in the group has a time period that can accommodate the process and the start time is earlier than the current equipment, the process will be adjusted to the earliest equipment that can start; If none of the selected processes has been adjusted to the processing equipment, then randomly select any process on other equipment in the equipment group where these processes are located, and exchange the processing equipment of the two processes. 9.如权利要求3所述的班次约束下带并行机作业车间的排产方法,其特征在于,通过嵌套蚁群算法对作业车间调度模型进行求解,将调度过程分为工序排序和设备选择两个阶段,工序排序阶段的蚁群为大蚁群,设备选择阶段的蚁群为小蚁群;大蚁群中的每只大蚁对工序进行遍历,获得工序排序方案,小蚁群根据每只大蚁获得的方案,对设备进行遍历,获得设备选择方案,从而形成完整解。9. The production scheduling method with parallel machine job shop under shift constraints as claimed in claim 3, wherein the job shop scheduling model is solved by nested ant colony algorithm, and the scheduling process is divided into process sequence and equipment selection In two stages, the ant colony in the process sorting stage is the large ant colony, and the ant colony in the equipment selection stage is the small ant colony; each large ant in the large ant colony traverses the process to obtain the process sorting plan, and the small ant colony is based on each ant colony. Only the solution obtained by the big ants is to traverse the equipment to obtain the equipment selection scheme, thereby forming a complete solution. 10.如权利要求9所述的班次约束下带并行机作业车间的排产方法,其特征在于,在工序排序阶段,优先安排加工时长与班次时长之比更大的工序在设备上加工;在设备选择阶段,优先选择同一设备组内开工时间早的设备。10. The production scheduling method with parallel machine workshops under shift constraints as claimed in claim 9, characterized in that, in the process sequencing stage, processes with a larger ratio of processing duration to shift duration are preferentially arranged to be processed on the equipment; In the equipment selection stage, the equipment with the earliest start time in the same equipment group is preferentially selected.
CN202210728425.9A 2022-06-24 2022-06-24 A production scheduling method for a job shop with parallel machines under shift constraints Pending CN115130744A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210728425.9A CN115130744A (en) 2022-06-24 2022-06-24 A production scheduling method for a job shop with parallel machines under shift constraints

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210728425.9A CN115130744A (en) 2022-06-24 2022-06-24 A production scheduling method for a job shop with parallel machines under shift constraints

Publications (1)

Publication Number Publication Date
CN115130744A true CN115130744A (en) 2022-09-30

Family

ID=83380170

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210728425.9A Pending CN115130744A (en) 2022-06-24 2022-06-24 A production scheduling method for a job shop with parallel machines under shift constraints

Country Status (1)

Country Link
CN (1) CN115130744A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118551873A (en) * 2024-01-31 2024-08-27 广州云谷数字科技有限公司 An intelligent production scheduling method for power switchgear based on PJS model

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118551873A (en) * 2024-01-31 2024-08-27 广州云谷数字科技有限公司 An intelligent production scheduling method for power switchgear based on PJS model

Similar Documents

Publication Publication Date Title
Shen et al. Mathematical modeling and multi-objective evolutionary algorithms applied to dynamic flexible job shop scheduling problems
CN104268722B (en) Dynamic flexible job-shop scheduling method based on multi-objective Evolutionary Algorithm
CN102385364B (en) Cross-operation-unit control method under flexible path
CN109872091B (en) A method and device for workpiece scheduling based on ant colony algorithm
CN117391423B (en) A multi-constraint automated scheduling method for chip high multi-layer ceramic packaging substrate production lines
CN111144710A (en) Construction and dynamic scheduling method of a sustainable mixed flow workshop
CN111401616B (en) Double-layer scheduling method for precast concrete component in supply chain environment
CN112907057A (en) Production scheduling optimization method and system based on improved MOPSO algorithm
CN110458326A (en) A Mixed Swarm Intelligent Optimization Method for Distributed Blocking Pipeline Scheduling
CN117891220A (en) A distributed hybrid flow shop scheduling method based on multi-agent deep reinforcement learning
CN115130744A (en) A production scheduling method for a job shop with parallel machines under shift constraints
CN118917567A (en) Multi-target distributed flexible job shop scheduling method based on hierarchical selection type deep reinforcement learning
Wang et al. Real-time decision support with reinforcement learning for dynamic flowshop scheduling
CN107437121A (en) Handle the production process control method of either simplex part simultaneously suitable for more machines
CN117872980A (en) A method, device and storage medium for scheduling a blocking mixed group flow shop
CN104636610B (en) A kind of manufacture system being applied under dynamic environment sends work Information revision method
CN118011978A (en) Flexible flow shop dynamic scheduling and maintenance optimization method considering random interference
CN114185312B (en) A multi-objective intelligent optimization algorithm for solving dynamic production scheduling problems in distributed flow shops
Jain et al. An approach for optimisation of flexible flow shop group scheduling with sequence dependent set-up time and preventive maintenance
Shalaby et al. New routing rules for dynamic flexible job shop scheduling with sequence-dependent setup times
Hung et al. Sensitivity search for the rescheduling of semiconductor photolithography operations
Li et al. AMulti-objective lot-streaming optimization scheduling model considering the blocking effect
Sha et al. A modified particle swarm optimization for multi-objective open shop scheduling
CN110991907B (en) Order service quality evaluation method based on finite state automaton
CN118333313B (en) Multi-objective flexible workshop dynamic scheduling method, device and electronic equipment under dual resource constraints

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