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 PDFInfo
- 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
Links
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total 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
技术领域 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和Sb;S31. 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_cb;S34. 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_fb;S35. 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,其中:
当工件数量为1,即n=1时,工件在每个机器上加工前的等待时间为0,因此总完工时间(makespan)为该工件在所有机器上的加工时间之和:
当工件数量为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:
若要使总完工时间(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,
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和Sb;S31. 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_cb;S34. 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_fb;S35. 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.
任意选择矩阵P中的两列元素P1、P2;Arbitrarily 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、P3;Arbitrarily 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、P4;Arbitrarily 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、P5;Arbitrarily 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、P3;Arbitrarily 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、P4;Arbitrarily 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、P5;Arbitrarily 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、P4;Arbitrarily 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、P5;Arbitrarily 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、P5;Arbitrarily 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列元素均已两两比较,各工件的最终加工顺序为P1P3P4P2P5。The 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)
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)
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)
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)
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 |
-
2012
- 2012-03-11 CN CN201210061685.1A patent/CN102566560B/en not_active Expired - Fee Related
Patent Citations (2)
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)
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 |