[go: up one dir, main page]

CN102566560B - A kind of production line scheduling method based on structure type heuritic approach - Google Patents

A kind of production line scheduling method based on structure type heuritic approach Download PDF

Info

Publication number
CN102566560B
CN102566560B CN201210061685.1A CN201210061685A CN102566560B CN 102566560 B CN102566560 B CN 102566560B CN 201210061685 A CN201210061685 A CN 201210061685A CN 102566560 B CN102566560 B CN 102566560B
Authority
CN
China
Prior art keywords
sum
processing sequence
matrix
workpiece
workpieces
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.)
Expired - Fee Related
Application number
CN201210061685.1A
Other languages
Chinese (zh)
Other versions
CN102566560A (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.)
Chengdu Information Technology Co Ltd of CAS
Original Assignee
Chengdu Information Technology Co Ltd of CAS
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 Chengdu Information Technology Co Ltd of CAS filed Critical Chengdu Information Technology Co Ltd of CAS
Priority to CN201210061685.1A priority Critical patent/CN102566560B/en
Publication of CN102566560A publication Critical patent/CN102566560A/en
Application granted granted Critical
Publication of CN102566560B publication Critical patent/CN102566560B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • General Factory Administration (AREA)

Abstract

本发明提供一种基于构造型启发式算法的生产线调度方法,该方法包括如下步骤:S1、若n个工件在m台机器上加工,设pi,j为第j个工件在第i台机器上的执行时间,构成矩阵P,其中i=1,2,Λ,m;j=1,2,Λ,n;S2、任意选择矩阵P中的两列元素,即工件a,b分别在m台机器上的执行时间Pa和Pb,其中1≤a,b≤n,a≠b;S3、确定工件a,b的加工顺序;S4、判断矩阵P中的n列元素是否均已两两比较,若是,则结束判断,按确定的加工顺序对工件进行调整,并依次在m台机器上进行加工,否则,返回步骤S2。本发明以总完工时间最小为目标来实现流水车间生产线的调度,通过对工件加工顺序的调整来减小每个工件在加工前的等待时间。相对于现有技术,本发明计算复杂度低,计算时间短,且具有较好的调度性能。

The present invention provides a production line scheduling method based on a constructional heuristic algorithm, which includes the following steps: S1. If n workpieces are processed on m machines, set p i, j as the jth workpiece on the i machine Execution time above constitutes a matrix P, where i=1, 2, Λ, m; j=1, 2, Λ, n; S2, select two columns of elements in the matrix P arbitrarily, that is, workpieces a and b are respectively in m Execution time P a and P b on two machines, where 1≤a, b≤n, a≠b; S3. Determine the processing sequence of workpieces a and b; S4. Determine whether the elements in the n columns of the matrix P Two comparisons, if yes, end the judgment, adjust the workpiece according to the determined processing sequence, and process it on m machines in turn, otherwise, return to step S2. The invention realizes the scheduling of the flow shop production line with the goal of minimizing the total completion time, and reduces the waiting time of each workpiece before processing by adjusting the processing sequence of the workpieces. Compared with the prior art, the present invention has low calculation complexity, short calculation time and better scheduling performance.

Description

一种基于构造型启发式算法的生产线调度方法A Production Line Scheduling Method Based on Constructive Heuristic Algorithm

技术领域 technical field

本发明涉及自动控制与信息技术领域,尤其涉及一种基于构造型启发式算法的生产线调度方法。The invention relates to the field of automatic control and information technology, in particular to a production line scheduling method based on a constructional heuristic algorithm.

背景技术 Background technique

生产线调度是制造企业生产过程中非常重要的问题,良好的调度策略将极大的提高生产效率。总完工时间(makespan)是调度过程中一个非常重要的性能指标,总完工时间最小可使得资源更加有效利用、任务更迅速传递及在制品库存最小。目前的生产线调度方法分为两种类型,一种为穷举法,如动态规划法、分支定界法等,虽然可以对工件的加工顺序得到最优的调度,但是这些方法的搜索空间会随着工件数量的增加呈指数式急剧增长,计算复杂度高,对机器硬件的要求较高,难以应用在大规模的生产线调度上;另一种方法为启发式算法,包括元启发式算法与构造型启发式算法,目前提出的构造型启发式算法如(1)Palmer算法——Palmer D S.Sequencing Jobs through a Multi-Stage Process in theMinimum Total Time-a Quick Method of Obtaining a near Optimum[J].Operational Research Quarterly,1965,16:101-107;(2)Gupta算法——Gupta J.A Functional Heuristic Algorithm for the Flowshop SchedulingProblem[J].Operational Research Quarterly,1971,22:39-47;(3)CDS算法——Campbell H G,Dudek R A,Smith M L.A Heuristic Algorithm for then-Job,m-Machine Scheduling Problem.[J].Management Science,1970,16:630-637;(4)RA算法——Dannenbring D G.An Evaluation of Flow ShopSequencing Heuristics[J].Management Science,1977,23(11):1174-1182;(5)NEH算法——Nawaz M,Enscore E,Ham I.A Heuristic Algorithm forthe m Machine,n Job Flow Shop[J].OMEGA:The International Journal ofManagement Sciences,1983,11(1):91-95等,以上几种构造型启发式算法中以NEH算法的性能最佳。然而,NEH算法在实现过程中需要通过多次计算拟定工件序列的总完工时间并进行比较,因此,计算复杂度将远大于其他构造型启发式算法。Production line scheduling is a very important issue in the production process of manufacturing enterprises. A good scheduling strategy will greatly improve production efficiency. The total completion time (makespan) is a very important performance index in the scheduling process. The minimum total completion time can lead to more efficient use of resources, faster delivery of tasks and minimum WIP inventory. The current production line scheduling methods are divided into two types, one is exhaustive method, such as dynamic programming method, branch and bound method, etc., although the processing sequence of workpieces can be optimally scheduled, but the search space of these methods will vary with time. As the number of workpieces increases exponentially, the calculation complexity is high, and the requirements for machine hardware are high, so it is difficult to apply to large-scale production line scheduling; another method is heuristic algorithm, including meta-heuristic algorithm and construction Constructive heuristic algorithms currently proposed such as (1) Palmer algorithm——Palmer D S. Sequencing Jobs through a Multi-Stage Process in the Minimum Total Time-a Quick Method of Obtaining a near Optimum[J]. Operational Research Quarterly, 1965, 16: 101-107; (2) Gupta algorithm - Gupta J.A Functional Heuristic Algorithm for the Flowshop Scheduling Problem [J]. Operational Research Quarterly, 1971, 22: 39-47; (3) CDS algorithm - —Campbell H G, Dudek R A, Smith M L.A Heuristic Algorithm for then-Job, m-Machine Scheduling Problem.[J].Management Science, 1970, 16:630-637; (4) RA Algorithm——Dannenbring D G .An Evaluation of Flow ShopSequencing Heuristics[J].Management Science, 1977, 23(11): 1174-1182; (5) NEH algorithm——Nawaz M, Enscore E, Ham I.A Heuristic Algorithm for the m Machine, n Job Flow Shop [J]. OMEGA: The International Journal of Management Sciences, 1983, 11(1): 91-95, etc. Among the above several constructive heuristic algorithms, the NEH algorithm has the best performance. However, the NEH algorithm needs to calculate and compare the total completion time of the proposed workpiece sequence multiple times during the implementation process, so the computational complexity will be much greater than other stereotyped heuristic algorithms.

发明内容 Contents of the invention

针对现有技术存在的问题,本发明的主要目的在于提供一种具有较高的调度性能,同时计算复杂度较低的基于构造型启发式算法的生产线调度方法。In view of the problems existing in the prior art, the main purpose of the present invention is to provide a production line scheduling method based on a constructional heuristic algorithm with high scheduling performance and low computational complexity.

为实现上述目的,本发明提供一种基于构造型启发式算法的生产线调度方法的实施例,该方法包括如下步骤:In order to achieve the above object, the present invention provides an embodiment of a production line scheduling method based on a constructive heuristic algorithm, the method comprising the following steps:

S1、若n个工件在m台机器上加工,设pi,j为第j个工件在第i台机器上的执行时间,构成矩阵P,其中i=1,2,Λ,m;j=1,2,Λ,n;S1. If n workpieces are processed on m machines, set p i, j as the execution time of the jth workpiece on the i machine, forming a matrix P, where i=1, 2, Λ, m; j= 1, 2, Λ, n;

S2、任意选择矩阵P中的两列元素,即工件a,b分别在m台机器上的执行时间Pa和Pb,其中1≤a,b≤n,a≠b;S2. Arbitrarily select two columns of elements in the matrix P, that is, the execution times P a and P b of the workpieces a and b on m machines respectively, where 1≤a, b≤n, a≠b;

S3、确定工件a,b的加工顺序;S3. Determine the processing sequence of workpieces a and b;

S4、判断矩阵P中的n列元素是否均已两两比较,若是,则结束判断,否则,返回步骤S2,S4, judge whether the n column elements in the matrix P have been compared in pairs, if so, then end the judgment, otherwise, return to step S2,

其中步骤S3包含如下步骤:Wherein step S3 comprises the following steps:

S31、将矩阵P中两列元素Pa,Pb的值分别代入计算得到Sa和SbS31. Substituting the values of elements P a and P b in the two columns of matrix P into Calculate S a and S b ;

S32、判断Sa*Sb<=0是否成立,若成立,则进入步骤S321,若不成立,则进入步骤S33;S32. Judging whether S a * S b <= 0 is true, if true, go to step S321, if not, go to step S33;

S321、判断Sa>0是否成立,若成立,则将元件a的加工顺序置于元件b之前,若不成立,则将元件b的加工顺序置于元件a之前;S321. Determine whether S a > 0 is true, if true, place the processing sequence of component a before component b, and if not, place the processing sequence of component b before component a;

S33、判断Sa,Sb>0是否成立,若成立,则进入步骤S34,若不成立,则进入步骤S35;S33. Judging whether S a , S b > 0 is true, if true, go to step S34, if not, go to step S35;

S34、将任意选择的矩阵P中的两列元素Pa,Pb的值分别代入计算得到sum_ca和sum_cbS34. Substituting the values of the two columns of elements P a and P b in the arbitrarily selected matrix P into Calculate sum_c a and sum_c b ;

S341、判断sum_ca<sum_cb是否成立,若成立,则将元件a的加工顺序置于元件b之前;若不成立,则进入步骤S342;S341. Determine whether sum_c a < sum_c b is true, if true, place the processing sequence of component a before component b; if not, proceed to step S342;

S342、判断sum_ca>sum_cb是否成立,若成立,则将元件b的加工顺序置于元件a之前,若不成立,则去掉Pa和Pb列的最末一个元素并返回步骤S31;S342, judging whether sum_c a > sum_c b is true, if true, place the processing sequence of component b before component a, if not true, remove the last element in the P a and P b columns and return to step S31;

S35、将任意选择的矩阵P中的两列元素Pa,Pb的值分别代入计算得到sum_fa和sum_fbS35. Substituting the values of the two columns of elements P a and P b in the arbitrarily selected matrix P into Calculate sum_f a and sum_f b ;

S351、判断sum_fa>sum_fb是否成立,若成立,则将元件a的加工顺序置于元件b之前,若不成立,则进入步骤S352;S351. Determine whether sum_f a > sum_f b is true, if true, place the processing sequence of component a before component b, if not, go to step S352;

S352、判断sum_fa<sum_fb是否成立,若成立,则将元件b的加工顺序置于元件a之前,若不成立,则去掉Pa和Pb列的首个元素并返回步骤S31。S352. Determine whether sum_f a < sum_f b is true. If true, place the processing sequence of component b before component a. If not, remove the first element in columns P a and P b and return to step S31.

进一步地,当判断矩阵P中的n列元素均已两两比较,则按确定的加工顺序对工件进行调整,并依次在m台机器上进行加工。Further, when the n columns of elements in the judgment matrix P have been compared pairwise, the workpieces are adjusted according to the determined processing sequence, and processed on m machines in sequence.

本发明以总完工时间最小为目标来实现流水车间生产线的调度,通过对工件加工顺序的调整来减小每个工件在加工前的等待时间。相对于现有技术,本发明计算复杂度低,计算时间短,且具有较好的调度性能。The invention realizes the scheduling of the flow shop production line with the goal of minimizing the total completion time, and reduces the waiting time of each workpiece before processing by adjusting the processing sequence of the workpieces. Compared with the prior art, the present invention has low calculation complexity, short calculation time and better scheduling performance.

附图说明 Description of drawings

图1是流水车间生产线调度示意图。Figure 1 is a schematic diagram of the production line scheduling in the flow shop.

图2是本发明基于构造型启发式算法的生产线调度方法流程图。Fig. 2 is a flow chart of the production line scheduling method based on the stereotype heuristic algorithm of the present invention.

图3是本发明基于构造型启发式算法的生产线调度方法中步骤S3的流程图。Fig. 3 is a flowchart of step S3 in the production line scheduling method based on the stereotype heuristic algorithm of the present invention.

具体实施方式 Detailed ways

下面结合附图,详细说明本发明的具体实施方式。The specific implementation manner of the present invention will be described in detail below in conjunction with the accompanying drawings.

如图1所示,为本发明的流水车间生产线调度示意图,表述为n个工件在m台机器上加工,即每个工件需要经过m道工序,每道工序要求不同的机器,但n个工件通过m台机器的顺序相同,即n个工件在每台机器上的加工顺序相同。即工件1、工件2......工件n按顺序分别经过机器1、机器2......机器m的加工,但这种加工顺序有可能并不是最佳的,为了使总加工时间最短,必须调整工件的加工顺序。As shown in Figure 1, it is a schematic diagram of the production line scheduling of the flow shop of the present invention, which is expressed as n workpieces are processed on m machines, that is, each workpiece needs to go through m processes, and each process requires different machines, but n workpieces The order of passing through m machines is the same, that is, the processing order of n workpieces on each machine is the same. That is, workpiece 1, workpiece 2...workpiece n are respectively processed by machine 1, machine 2...machine m in sequence, but this processing sequence may not be optimal. In order to make the total The processing time is the shortest, and the processing sequence of the workpieces must be adjusted.

定义Oi,j为第j个工件在第i台机器上操作,pi,j为Oi,j的执行时间,ci,j为Oi,j的完成时间,ci,j为Oi,j加等待时间。其中i=1,2,Λ,m;j=1,2,Λ,n。为了确定n个工件的加工任务在每台机器上的最优加工顺序,使所有工件的加工任务全部完工的时间最短。在进行工件加工时,每个加工任务在机器上的加工顺序相同,均为1,2,Λ,m;每台机器同时只能进行一个加工任务;一个加工任务不能同时在不同的机器上进行;各加工任务在加工完后立即送下一道工序;任务在机器上开始加工,必须一直进行到该工序完工,中途不允许停下来插入其它任务;所有任务在0时刻已准备就绪,机器调整时间包括在加工时间内;允许任务在工序之间等待;允许机器在任务未到达时闲置。最大完成时间makespan(Cmax)也即是最后一个工件加工的完成时间cm,n,找到一个工件的加工顺序,工件a,工件b......工件n,使得总完工时间(makespan)最小。Define O i, j as the j-th workpiece operating on the i machine, p i, j is the execution time of O i, j, c i, j is the completion time of O i, j, and c i, j is O i, j plus waiting time. where i=1, 2, Λ, m; j=1, 2, Λ, n. In order to determine the optimal processing sequence of the processing tasks of n workpieces on each machine, the time to complete the processing tasks of all workpieces is the shortest. When processing workpieces, the processing sequence of each processing task on the machine is the same, which is 1, 2, Λ, m; each machine can only perform one processing task at the same time; one processing task cannot be performed on different machines at the same time ; Each processing task is sent to the next process immediately after processing; the task starts to process on the machine, and must continue until the process is completed, and it is not allowed to stop and insert other tasks in the middle; all tasks are ready at time 0, and the machine adjusts the time Included in processing time; allowing tasks to wait between operations; allowing machines to sit idle while tasks do not arrive. The maximum completion time makespan(C max ) is also the completion time c m, n of the last workpiece processing. Find the processing sequence of a workpiece, workpiece a, workpiece b...job n, so that the total completion time (makespan ) minimum.

设工件数量为i时的总完工时间(makespan)为Ci,其中1≤i≤n,当i=n时,定义Ci=Cn=C;设第j个工件在第i台机器上加工前的等待时间为ti,j,其中ti,1=0,1≤i≤m,1≤j≤n;设第n个工件在第i台机器上加工前的等待时间为Di,其中D1=0,即所有工件在第一台机器上的等待时间均为0,1≤i≤m;设n个工件在m台机器上的加工时间矩阵为P,其中: P = p 1,1 p 1,2 &Lambda; p 1 , n p 2,1 p 2,2 &Lambda; p 2 , n M M O M p m , 1 p m , 2 &Lambda; p m , n , 在矩阵P中,pi,j表示第j个工件在第i台机器上的加工时间,1≤i≤m,1≤j≤n。Let the total completion time (makespan) when the number of workpieces is i be C i , where 1≤i≤n, when i=n, define C i =C n =C; suppose the jth workpiece is on the i machine The waiting time before processing is t i, j , where t i, 1 = 0, 1≤i≤m, 1≤j≤n; let the waiting time of the nth workpiece before being processed on the i machine be D i , where D 1 =0, that is, the waiting time of all workpieces on the first machine is 0, 1≤i≤m; let the processing time matrix of n workpieces on m machines be P, where: P = p 1,1 p 1,2 &Lambda; p 1 , no p 2,1 p 2,2 &Lambda; p 2 , no m m o m p m , 1 p m , 2 &Lambda; p m , no , In the matrix P, p i, j represent the processing time of the jth workpiece on the i machine, 1≤i≤m, 1≤j≤n.

当工件数量为1,即n=1时,工件在每个机器上加工前的等待时间为0,因此总完工时间(makespan)为该工件在所有机器上的加工时间之和: C 1 = p 1,1 + p 2,1 + &Lambda; + p m , 1 = &Sigma; i = 1 m p i , 1 ; When the number of workpieces is 1, that is, when n=1, the waiting time before the workpiece is processed on each machine is 0, so the total completion time (makespan) is the sum of the processing time of the workpiece on all machines: C 1 = p 1,1 + p 2,1 + &Lambda; + p m , 1 = &Sigma; i = 1 m p i , 1 ;

当工件数量为2,即n=2时,总完工时间(makespan)为第1个工件在所有机器上的加工时间之和、第二个工件在最后一台机器上的加工时间与第二个工件在最后一台机器上的加工前等待时间之和:C2=C1+pm,2+D2,而D2=t2,n,可按如下方式进行计算:t2,2=max(p2,1-p1,2,0),t2,3=max(t2,2+p2,2-p′1,3,0),Λ,t2,n=max(t2,n-1+p2,n-1-p′1,n,0),其中p′i,j=pi,j+max(pi-1,j-pj,i-1,0),2≤i≤m,1≤j≤n;When the number of workpieces is 2, that is, when n=2, the total completion time (makespan) is the sum of the processing time of the first workpiece on all machines, the processing time of the second workpiece on the last machine and the second The sum of the waiting time of the workpiece on the last machine before processing: C 2 =C 1 +p m,2 +D 2 , and D 2 =t 2,n can be calculated as follows: t 2,2 = max(p 2,1 -p 1,2,0 ), t 2,3 =max(t 2,2 +p 2,2 -p' 1,3,0 ), Λ,t 2,n =max( t 2,n-1 +p 2,n-1 -p′ 1,n ,0), where p′ i,j =p i,j +max(p i-1,j -p j,i-1 , 0), 2≤i≤m, 1≤j≤n;

以此类推,当工件数量为n时,Cn=Cn-1+pm,n+Dn,而Dn=tm,n,可按如下方式计算:tm,2=max(pm,1-p′m-1,2,0),tm,3=max(tm,2+pm,2-p′m-1,3,0),Λ,tm,n=max(tm,n-1+pm,n-1-p′m-1,n,0)。即:By analogy, when the number of workpieces is n, C n =C n-1 +p m,n +D n , and D n =t m,n can be calculated as follows: t m,2 =max(p m, 1 -p' m-1, 2 , 0), t m, 3 = max(t m, 2 +p m, 2 -p' m-1, 3 , 0), Λ, t m, n = max(t m, n-1 + p m, n-1 -p' m-1, n , 0). Right now:

makespanmake span == CC ==

CC nno == CC nno -- 11 ++ pp mm ,, nno ++ DD. nno

== CC nno -- 22 ++ pp mm -- 11 ,, nno ++ pp mm ,, nno ++ DD. nno -- 11 ++ DD. nno

== CC nno -- 33 ++ pp mm -- 22 ,, nno ++ pp mm -- 11 ,, nno ++ pp mm ,, nno ++ DD. nno -- 22 ++ DD. nno -- 11 ++ DD. nno

Mm

== pp 1,11,1 ++ pp 2,12,1 ++ &Lambda;&Lambda; ++ pp mm ,, 11 ++ pp mm ,, 22 ++ pp mm ,, 33 ++ &Lambda;&Lambda; ++ pp mm ,, nno ++ DD. 11 ++ DD. 22 ++ &Lambda;&Lambda; ++ DD. nno

== &Sigma;&Sigma; ii == 11 mm pp ii ,, 11 ++ &Sigma;&Sigma; jj == 22 nno pp mm ,, jj ++ &Sigma;&Sigma; kk == 11 nno DD. kk

若要使总完工时间(makespan)最小,则需找一个工件加工序列,使得第n个工件在第m台机器上的完成时间尽量接近或等于即要缩小第一个工件在所有机器的加工时间之和、除第一个工件之外的其他工件在最后一台机器的加工时间之和以及除第一个工件之外的其他工件在加工前的等待时间之和。To minimize the total completion time (makespan), it is necessary to find a workpiece processing sequence so that the completion time of the nth workpiece on the mth machine is as close to or equal to That is to reduce the sum of the processing time of the first workpiece on all machines, the sum of the processing time of other workpieces except the first workpiece on the last machine, and the processing time of other workpieces except the first workpiece before processing sum of waiting times.

如图2所示,为本发明基于构造型启发式算法的生产线调度方法流程图,该方法包括如下步骤:As shown in Figure 2, it is a flow chart of the production line scheduling method based on the constructional heuristic algorithm of the present invention, and the method includes the following steps:

S1、若n个工件在m台机器上加工,设pi,j为第j个工件在第i台机器上的执行时间,构成矩阵P, P = p 1,1 p 1,2 &Lambda; p 1 , n p 2,1 p 2,2 &Lambda; p 2 , n M M O M p m , 1 p m , 2 &Lambda; p m , n , 其中i=1,2,Λ,m;j=1,2,Λ,n;S1. If n workpieces are processed on m machines, let p i, j be the execution time of the jth workpiece on the i machine, forming a matrix P, P = p 1,1 p 1,2 &Lambda; p 1 , no p 2,1 p 2,2 &Lambda; p 2 , no m m o m p m , 1 p m , 2 &Lambda; p m , no , Wherein i=1,2,Λ,m; j=1,2,Λ,n;

S2、任意选择矩阵P中的两列元素,即工件a,b分别在m台机器上的执行时间Pa和Pb,其中1≤a,b≤n,a≠b;S2. Arbitrarily select two columns of elements in the matrix P, that is, the execution times P a and P b of the workpieces a and b on m machines respectively, where 1≤a, b≤n, a≠b;

S3、确定工件a,b的加工顺序;S3. Determine the processing sequence of workpieces a and b;

S4、判断矩阵P中的n列元素是否均已两两比较,若是,则结束判断,按确定的加工顺序对工件进行调整,并依次在m台机器上进行加工,否则,返回步骤S2。S4. Judging whether the elements in the n columns in the matrix P have been compared in pairs. If so, end the judgment, adjust the workpiece according to the determined processing sequence, and process it on m machines in turn, otherwise, return to step S2.

如图3所示,为本发明基于构造型启发式算法的生产线调度方法中步骤S3的流程图。步骤S3,即确定工件a,b的加工顺序,具体包含如下步骤:As shown in FIG. 3 , it is a flowchart of step S3 in the production line scheduling method based on the stereotype heuristic algorithm of the present invention. Step S3, that is to determine the processing sequence of workpieces a and b, specifically includes the following steps:

S31、将矩阵P中两列元素Pa,Pb的值分别代入(Sj为矩阵P中第j列的加权值),计算得到Sa和SbS31. Substituting the values of elements P a and P b in the two columns of matrix P into (S j is the weighted value of column j in the matrix P), calculate S a and S b ;

S32、判断Sa*Sb<=0是否成立,若成立,则进入步骤S321,若不成立,则进入步骤S33;S32. Judging whether S a * S b <= 0 is true, if true, go to step S321, if not, go to step S33;

S321、判断Sa>0是否成立,若成立,则进入步骤S322,即将元件a的加工顺序置于元件b之前,若不成立,则进入步骤S323,即将元件b的加工顺序置于元件a之前;S321. Judging whether S a > 0 is true, if true, proceed to step S322, that is, place the processing sequence of component a before component b, and if not, proceed to step S323, that is, place the processing sequence of component b before component a;

S33、判断Sa,Sb>0是否成立,若成立,则进入步骤S34,若不成立,则进入步骤S35;S33. Judging whether S a , S b > 0 is true, if true, go to step S34, if not, go to step S35;

S34、将任意选择的矩阵P中的两列元素Pa,Pb,代入(sum_cj为不断去掉矩阵P中第j列最后一个元素的和),计算得到sum_ca和sum_cbS34. Substituting the elements P a and P b in the two columns of the matrix P selected arbitrarily into (sum_c j is the sum of constantly removing the last element of column j in the matrix P), calculate sum_c a and sum_c b ;

S341、判断sum_ca<sum_cb是否成立,若成立,则进入步骤S343,即将元件a的加工顺序置于元件b之前;若不成立,则进入步骤S342;S341. Judging whether sum_c a < sum_c b is true, if true, proceed to step S343, that is, place the processing sequence of component a before component b; if not, proceed to step S342;

S342、判断sum_ca>sum_cb是否成立,若成立,则进入步骤S345,即将元件b的加工顺序置于元件a之前,若不成立,可知sum_ca=sum_cb,则进入步骤S344,即去掉Pa和Pb列的最末一个元素并返回步骤S31;S342. Determine whether sum_c a > sum_c b is true, if true, enter step S345, that is, place the processing sequence of component b before component a, if not, it can be seen that sum_c a = sum_c b , then enter step S344, that is, remove P a and the last element of the P b column and return to step S31;

S35、将任意选择的矩阵P中的两列元素Pa,Pb,代入(sum_fj为不断去掉矩阵P中第j列最前一个元素的和),计算得到sum_fa和sum_fbS35. Substituting the elements P a and P b in the two columns of the matrix P selected arbitrarily into (sum_f j is to constantly remove the sum of the first element in the jth column in the matrix P), calculate sum_f a and sum_f b ;

S351、判断sum_fa>sum_fb是否成立,若成立,则进入步骤S353,即将元件a的加工顺序置于元件b之前,若不成立,则进入步骤S352;S351. Determine whether sum_f a > sum_f b is true, if true, proceed to step S353, that is, place the processing sequence of component a before component b, and if not, proceed to step S352;

S352、判断sum_fa<sum_fb是否成立,若成立,则进入步骤S355,即将元件b的加工顺序置于元件a之前,若不成立,可知sum_fa=sum_fb,则进入步骤S354,即去掉Pa和Pb列的首个元素并返回步骤S31。S352. Determine whether sum_f a < sum_f b is true. If true, proceed to step S355, that is, place the processing sequence of component b before component a. If not, it can be known that sum_f a = sum_f b , then proceed to step S354, that is, remove P a and the first element of column P b and return to step S31.

通过实验结果说明本发明的基于构造型启发式算法的生产线调度方法的实际调度性能。其中Y为最优率,即在特定样本规模中,使用某种调度方法得出最优解占所有样本数量的比例,该数值越高,则表明使用该调度方法得出最优解的比例越高,性能更优。另外,最优偏差率E=(MK-MKp)·100/MKp,其中MK为使用某种调度方法排序后工件的实际总完工时间,MKp为该组工件在最优加工顺序时的总完工时间,最优偏差率E越小则表明使用某种调度方法对工件加工顺序的排列越好,而E=0则表明某次工件的加工顺序排列为最优。The actual scheduling performance of the production line scheduling method based on the stereotype heuristic algorithm of the present invention is illustrated by the experimental results. Among them, Y is the optimal rate, that is, in a specific sample size, the proportion of the optimal solution obtained by using a certain scheduling method accounts for the proportion of all samples. The higher the value, the higher the proportion of the optimal solution obtained by using the scheduling method. High, better performance. In addition, the optimal deviation rate E=(MK-MKp)·100/MKp, where MK is the actual total completion time of the workpieces sorted by a certain scheduling method, and MKp is the total completion time of the group of workpieces in the optimal processing order , the smaller the optimal deviation rate E is, the better the arrangement of the workpiece processing sequence is by using a certain scheduling method, and E=0 indicates that the processing sequence of a certain workpiece is optimal.

实验所用到的工件在各机器上的加工时间矩阵P中的各个元素Pi,j为0到9的随机整数(1≤i≤m,1≤j≤n),下表给出了多个以样本规模20为单位的测试数据,从下表可以看出,利用本发明的基于构造型启发式算法的生产线调度方法最优偏差率E>10%的非常少,平均最优偏差率基本都在5%以下,表明该方法的调度性能较优。另外,本发明的基于构造型启发式算法的生产线调度方法算法复杂度为nlog(n)+nm,计算复杂度低,计算时间短。Each element Pi in the processing time matrix P of the workpiece used in the experiment on each machine, j is a random integer from 0 to 9 (1≤i≤m, 1≤j≤n), the following table gives a number of The test data with a sample size of 20 as a unit can be seen from the following table, using the production line scheduling method based on the constructional heuristic algorithm of the present invention, the optimal deviation rate E>10% is very small, and the average optimal deviation rate is basically in the Below 5%, it shows that the scheduling performance of this method is better. In addition, the algorithm complexity of the production line scheduling method based on the constructional heuristic algorithm of the present invention is nlog(n)+nm, the calculation complexity is low, and the calculation time is short.

实施例一Embodiment one

假设5个工件在4台机器上的加工时间矩阵P如下,Pi表示P的第i列。Assuming that the processing time matrix P of 5 workpieces on 4 machines is as follows, P i represents the i-th column of P.

PP == 11 11 33 77 77 33 99 99 22 33 99 77 33 77 77 1010 66 88 77 22

任意选择矩阵P中的两列元素P1、P2Arbitrarily select two columns of elements P 1 and P 2 in the matrix P;

将P1、P2的值分别代入计算得到S1=33,S2=13;Substitute the values of P 1 and P 2 into Calculated to get S 1 =33, S 2 =13;

因为S1,S2>0,将P1、P2的值分别代入计算得到sum_C1=13,sum_C2=17;Since S 1 and S 2 >0, substitute the values of P 1 and P 2 into Calculate sum_C 1 =13, sum_C 2 =17;

因为sum_C1<sum_C2,因此将P1的加工顺序置于P2之前;Because sum_C 1 < sum_C 2 , the processing sequence of P 1 is placed before P 2 ;

任意选择矩阵P中的两列元素P1、P3Arbitrarily select two columns of elements P 1 and P 3 in the matrix P;

将P1、P3的值分别代入计算得到S1=33,S3=9;Substitute the values of P 1 and P 3 into Calculated to get S 1 =33, S 3 =9;

因为S1,S3>0,将P1、P3的值分别代入计算得到sum_C1=13,sum_C3=15;Since S 1 and S 3 >0, substitute the values of P 1 and P 3 into Calculate sum_C 1 =13, sum_C 3 =15;

因为sum_C1<sum_C3,因此将P1的加工顺序置于P3之前;Because sum_C 1 < sum_C 3 , the processing sequence of P 1 is placed before P3;

任意选择矩阵P中的两列元素P1、P4Arbitrarily select two columns of elements P 1 and P 4 in the matrix P;

将P1、P4的值分别代入计算得到S1=33,S4=4;Substitute the values of P 1 and P 4 into Calculated to get S 1 =33, S 4 =4;

因为S1,S4>0,将P1、P4的值分别代入计算得到sum_C1=13,sum_C4=16;Since S 1 and S 4 >0, substitute the values of P 1 and P 4 into Calculate sum_C 1 =13, sum_C 4 =16;

因为sum_C1<sum_C4,因此将P1的加工顺序置于P4之前;Because sum_C 1 < sum_C 4 , the processing sequence of P 1 is placed before P 4 ;

任意选择矩阵P中的两列元素P1、P5Arbitrarily select two columns of elements P 1 and P 5 in the matrix P;

将P1、P5的值分别代入计算得到S1=33,S5=-10;Substitute the values of P 1 and P 5 into Calculated to get S 1 =33, S 5 =-10;

因为S1·S5<=0,且S1>0,因此将P1的加工顺序置于P5之前;Since S 1 ·S 5 <=0, and S 1 >0, the processing sequence of P 1 is placed before P 5 ;

任意选择矩阵P中的两列元素P2、P3Arbitrarily select two columns of elements P 2 and P 3 in the matrix P;

将P2、P3的值分别代入计算得到S2=13,S3=9;Substitute the values of P 2 and P 3 into Calculated to get S 2 =13, S 3 =9;

因为S2,S3>0,将P2、P3的值分别代入计算得到sum_C2=17,sum_C3=15;Since S 2 and S 3 >0, substitute the values of P 2 and P 3 into Calculate sum_C 2 =17, sum_C 3 =15;

因为sum_C2>sum_C3,因此将P3的加工顺序置于P2之前;Because sum_C 2 > sum_C 3 , the processing sequence of P 3 is placed before P 2 ;

任意选择矩阵P中的两列元素P2、P4Arbitrarily select two columns of elements P 2 and P 4 in the matrix P;

将P2、P4的值分别代入计算得到S2=13,S4=4;Substitute the values of P 2 and P 4 into It is calculated that S 2 =13, S 4 =4;

因为S2,S4>0,将P2、P4的值分别代入计算得到sum_C2=17,sum_C4=16;Because S 2 , S 4 >0, substitute the values of P 2 and P 4 into Calculate sum_C 2 =17, sum_C 4 =16;

因为sum_C2>sum_C4,因此将P4的加工顺序置于P2之前;Because sum_C 2 > sum_C 4 , the processing sequence of P 4 is placed before P 2 ;

任意选择矩阵P中的两列元素P2、P5Arbitrarily select two columns of elements P 2 and P 5 in the matrix P;

将P2、P5的值分别代入计算得到S2=13,S5=-10;Substitute the values of P 2 and P 5 into Calculated to get S 2 =13, S 5 =-10;

因为S2·S5<=0,且S2>0,因此将P2的加工顺序置于P5之前;Since S 2 ·S 5 <=0, and S 2 >0, the processing sequence of P 2 is placed before P 5 ;

任意选择矩阵P中的两列元素P3、P4Arbitrarily select two columns of elements P 3 and P 4 in the matrix P;

将P3、P4的值分别代入计算得到S3=9,S4=4;Substitute the values of P 3 and P 4 into Calculated to get S 3 =9, S 4 =4;

因为S3,S4>0,将P3、P4的值分别代入计算得到sum_C3=15,sum_C4=16;Since S 3 and S 4 >0, substitute the values of P 3 and P 4 into Calculate sum_C 3 =15, sum_C 4 =16;

因为sum_C3<sum_C4,因此将P3的加工顺序置于P4之前;Because sum_C 3 < sum_C 4 , the processing sequence of P 3 is placed before P 4 ;

任意选择矩阵P中的两列元素P3、P5Arbitrarily select two columns of elements P 3 and P 5 in the matrix P;

将P3、P5的值分别代入计算得到S3=9,S5=-10;Substitute the values of P 3 and P 5 into Calculated to get S 3 =9, S 5 =-10;

因为S3·S5<=0,且S3>0,因此将P3的加工顺序置于P5之前;Since S 3 ·S 5 <=0, and S 3 >0, the processing sequence of P 3 is placed before P 5 ;

任意选择矩阵P中的两列元素P4、P5Arbitrarily select two columns of elements P 4 and P 5 in the matrix P;

将P4、P5的值分别代入计算得到S4=4,S5=-10;Substitute the values of P 4 and P 5 into Calculated to get S 4 =4, S 5 =-10;

因为S4·S5<=0,且S4>0,因此将P4的加工顺序置于P5之前;Since S 4 ·S 5 <=0, and S 4 >0, the processing sequence of P 4 is placed before P 5 ;

矩阵P中的5列元素均已两两比较,各工件的最终加工顺序为P1P3P4P2P5The elements in the 5 columns in the matrix P have been compared in pairs, and the final processing sequence of each workpiece is P 1 P 3 P 4 P 2 P 5 .

调整工件的加工顺序为P1P3P4P2P5,依次在m台机器上进行加工。Adjust the processing sequence of workpieces to P 1 P 3 P 4 P 2 P 5 , and process them on m machines in turn.

以上介绍了基于构造型启发式算法的生产线调度方法。本发明并不限定于以上实施例,任何未脱离本发明技术方案,即仅仅对其进行本领域普通技术人员所知悉的改进或变更,均属于本发明的保护范围之内。The production line scheduling method based on the constructional heuristic algorithm is introduced above. The present invention is not limited to the above embodiments, and any improvement or change that does not deviate from the technical solution of the present invention, that is, only improvements or changes known to those skilled in the art, falls within the protection scope of the present invention.

Claims (2)

1., based on a production line scheduling method for structure type heuritic approach, comprise following steps:
If a S1 n workpiece is processed on m platform machine, if p i,jfor the execution time of a jth workpiece on i-th machine, form matrix P, wherein i=1,2 ..., m; J=1,2 ..., n;
Two column elements in S2, arbitrarily selection matrix P, i.e. workpiece a, the b execution time P respectively on m platform machine aand P b, wherein 1≤a≤n, 1≤b≤n, a ≠ b;
S3, determine workpiece a, the processing sequence of b;
Whether the n column element in S4, judgment matrix P all compares between two, if so, then terminates to judge, otherwise, return step S2,
It is characterized in that, described step S3 comprises following steps:
S31, by two column element P in matrix P a, P bvalue substitute into respectively calculate S aand S b;
S32, judge S a* S bwhether <=0 sets up, if set up, then enters step S321, if be false, then enters step S33;
S321, judge S awhether >0 sets up, if set up, then before the processing sequence of element a being placed in element b, if be false, then before the processing sequence of element b being placed in element a;
S33, judge S a>0 and S bwhether >0 sets up, if set up, then enters step S34, if be false, then enters step S35;
S34, by two column element P in optional matrix P a, P bvalue substitute into respectively calculate sum_c aand sum_c b;
S341, judge sum_c a<sum_c bwhether set up, if set up, then before the processing sequence of element a being placed in element b; If be false, then enter step S342;
S342, judge sum_c a>sum_c bwhether set up, if set up, then, before the processing sequence of element b being placed in element a, if be false, then remove P aand P bthe last element arranged also returns step S31;
S35, by two column element P in optional matrix P a, P bvalue substitute into respectively calculate sum_f aand sum_f b;
S351, judge sum_f a>sum_f bwhether set up, if set up, then, before the processing sequence of element a being placed in element b, if be false, then enter step S352;
S352, judge sum_f a<sum_f bwhether set up, if set up, then, before the processing sequence of element b being placed in element a, if be false, then remove P aand P bthe first element arranged also returns step S31.
2. as claimed in claim 1 based on the production line scheduling method of structure type heuritic approach, it is characterized in that: when the n column element in judgment matrix P all compares between two, then by the processing sequence determined, workpiece is adjusted, and process on m platform machine successively.
CN201210061685.1A 2012-03-11 2012-03-11 A kind of production line scheduling method based on structure type heuritic approach Expired - Fee Related CN102566560B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210061685.1A CN102566560B (en) 2012-03-11 2012-03-11 A kind of production line scheduling method based on structure type heuritic approach

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210061685.1A CN102566560B (en) 2012-03-11 2012-03-11 A kind of production line scheduling method based on structure type heuritic approach

Publications (2)

Publication Number Publication Date
CN102566560A CN102566560A (en) 2012-07-11
CN102566560B true CN102566560B (en) 2015-07-29

Family

ID=46412204

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210061685.1A Expired - Fee Related CN102566560B (en) 2012-03-11 2012-03-11 A kind of production line scheduling method based on structure type heuritic approach

Country Status (1)

Country Link
CN (1) CN102566560B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103246240B (en) * 2012-10-30 2015-09-30 中国科学院沈阳自动化研究所 A kind of PREDICTIVE CONTROL dispatching method ratifying standby single process continuously
CN103676902B (en) * 2013-12-20 2016-08-17 东北大学 A kind of Flow Shop rescheduling method
CN103926910B (en) * 2014-05-06 2016-08-17 安徽工程大学 There is equipment unavailable time section constraint production line scheduling control method
CN104391488A (en) * 2014-11-18 2015-03-04 广东工业大学 Optimizing and dispatching method of energy consumption of flexible flow shop with associated adjustment time and sequence
CN105235271B (en) * 2015-11-20 2017-04-12 合肥合锻智能制造股份有限公司 Hydropress automatic production line robot dispatching method based on minimum waiting time
CN108107848B (en) * 2016-11-24 2020-05-22 江苏创源电子有限公司 A method of assembly line shop scheduling based on minimum idle time
CN108108829A (en) * 2016-11-24 2018-06-01 江苏创源电子有限公司 A kind of job-shop scheduling method based on improvement drosophila algorithm
CN111105062B (en) * 2018-10-26 2024-06-25 北京北方华创微电子装备有限公司 Material scheduling method and material scheduling system based on equipment state
CN118605429B (en) * 2024-08-08 2024-11-22 天翼物联科技有限公司 Production line scheduling method based on heuristic algorithm apparatus, device and medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1487381A (en) * 2003-07-17 2004-04-07 上海交通大学 Heuristic Scheduling Method of Mixed Flow Production Line Based on Parameter Space Search
CN1776554A (en) * 2005-10-20 2006-05-24 同济大学 Recombinative production line scheduling method based on genetic algorithm

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2739858B2 (en) * 1996-03-25 1998-04-15 日本電気株式会社 Production control system and device
US5706200A (en) * 1996-08-15 1998-01-06 The Board Of Trustees Of The University Of Il. Scheduling system and scheduling method for reentrant line processes
AU2001281194A1 (en) * 2000-08-09 2002-02-18 Myteam.Com, Inc. Electronic concession stand

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1487381A (en) * 2003-07-17 2004-04-07 上海交通大学 Heuristic Scheduling Method of Mixed Flow Production Line Based on Parameter Space Search
CN1776554A (en) * 2005-10-20 2006-05-24 同济大学 Recombinative production line scheduling method based on genetic algorithm

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于自适应构件的工作流流程动态变更模型;魏乐等;《计算机集成制造系统》;20101231;第16卷(第12期);第2603-2610页 *

Also Published As

Publication number Publication date
CN102566560A (en) 2012-07-11

Similar Documents

Publication Publication Date Title
CN102566560B (en) A kind of production line scheduling method based on structure type heuritic approach
Wang et al. A hybrid local-search algorithm for robust job-shop scheduling under scenarios
Ying Solving non-permutation flowshop scheduling problems by an effective iterated greedy heuristic
Laha et al. A Hungarian penalty-based construction algorithm to minimize makespan and total flow time in no-wait flow shops
CN106294288B (en) A kind of distribution non-negative matrix factorization method
Hidri et al. Bounding strategies for the hybrid flow shop scheduling problem
İşler et al. Scheduling in a two-machine flow-shop for earliness/tardiness under learning effect
CN115640898A (en) A Large-Scale Flexible Job Shop Scheduling Method Based on DDQN Algorithm
DE112022000723T5 (en) BRANCHING PROCESS FOR A CIRCUIT OF A NEURONAL PROCESSOR
Li et al. The distributed permutation flowshop scheduling problem with different transport timetables and loading capacities
CN114595633A (en) A multi-constraint-based multi-objective flexible job shop energy-saving scheduling method
CN110750079A (en) Hybrid flow shop scheduling optimization method allowing process jump
Riahi et al. Tailoring customer order scheduling search algorithms
Balaji et al. Artificial immune system algorithm and simulated annealing algorithm for scheduling batches of parts based on job availability model in a multi-cell flexible manufacturing system
Wang et al. A variable-representation discrete artificial bee colony algorithm for a constrained hybrid flow shop
Deng et al. Minimizing mean completion time in a batch processing system
CN117872980A (en) A method, device and storage medium for scheduling a blocking mixed group flow shop
Hadda et al. The two-stage assembly flow shop scheduling with an availability constraint: worst case analysis
CN108107848B (en) A method of assembly line shop scheduling based on minimum idle time
CN115438929A (en) Scheduling method for reentrant process of aeroengine assembly workshop
Piroozfard et al. A hybrid harmony search algorithm for the job shop scheduling problems
Laha et al. An efficient heuristic approach to total flowtime minimization in permutation flowshop scheduling
Futatsuishi et al. A study of the multi-stage flowshop scheduling problem with alternative operation assignments
Wang et al. A hybrid forward/backward approach for single batch scheduling problems with non-identical job sizes
CN104881359A (en) Method for automatically generating test data realizing path covering through linear fitting

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150729

Termination date: 20180311