CN116300763B - Mixed flow shop mathematical heuristic scheduling method and system considering machine configuration - Google Patents

Mixed flow shop mathematical heuristic scheduling method and system considering machine configuration Download PDF

Info

Publication number
CN116300763B
CN116300763B CN202310337907.6A CN202310337907A CN116300763B CN 116300763 B CN116300763 B CN 116300763B CN 202310337907 A CN202310337907 A CN 202310337907A CN 116300763 B CN116300763 B CN 116300763B
Authority
CN
China
Prior art keywords
machine
workpiece
stage
solution
sequence
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.)
Active
Application number
CN202310337907.6A
Other languages
Chinese (zh)
Other versions
CN116300763A (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.)
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 CN202310337907.6A priority Critical patent/CN116300763B/en
Publication of CN116300763A publication Critical patent/CN116300763A/en
Application granted granted Critical
Publication of CN116300763B publication Critical patent/CN116300763B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/418Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
    • G05B19/41865Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by job scheduling, process planning, material flow
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/32Operator till task planning
    • G05B2219/32252Scheduling production, machining, job shop

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Quality & Reliability (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • General Factory Administration (AREA)

Abstract

The invention belongs to the technical field of workshop scheduling, and discloses a mathematical heuristic scheduling method and system for a mixed flow workshop taking machine configuration into consideration, wherein the method comprises the following steps: (1) Constructing a mixed integer linear programming model for a mixed flow shop considering machine configuration; (2) Determining a coding structure of a mathematical heuristic algorithm based on decision variables of the model, and decoding to obtain an optimal scheduling solution of the current round; (3) Updating the solutions in each iteration elite population to replace the last suboptimal solution; (4) Determining a transmission rule of an elite population to a mixed integer linear programming model transmission variable based on an example scale, and determining a sub-model corresponding to the transmission rule; (5) Solving the sub-model, wherein each solved value is used as the fitness function value of the current transfer chromosome and enters the population to perform superior and inferior elimination, so as to obtain the final dispatching optimal solution. The invention improves the searching quality of the algorithm and realizes a scheduling scheme for efficiently acquiring the optimal target.

Description

Mixed flow shop mathematical heuristic scheduling method and system considering machine configuration
Technical Field
The invention belongs to the technical field related to workshop scheduling, and particularly relates to a mathematical heuristic scheduling method and system for a mixed flow workshop taking machine configuration into consideration.
Background
In modern manufacturing, order scheduling in a production management process is closely related to the operating conditions of the entire plant. Scheduling is a scheduling scheme meeting actual requirements under limited production resources and workshop resources. At present, due to the small-batch and customized production mode of modern factories, orders often have the characteristics of scattered and large fluctuation along with time, and enterprises generally preferentially meet the production requirement of the largest order in order to ensure the maximization of benefits, so that workshop resources are redundant for the order resources. Meanwhile, a Hybrid Flow Shop (HFS) type is taken as one of typical workshop types of a modern automatic factory, and is widely applied to industries such as food, electronic manufacturing and chemical industry, so that the scheduling problem (Hybrid Flow Shop Scheduling WITH MACHINE Configurations, HFSP-MC) of the Hybrid Flow Shop of machine configuration is considered, the energy is saved by optimizing the use quantity of machines, the purpose of solving the fluctuation of orders and the difference of the saturation of workshops is achieved, the number of machines per stage, the workpiece sequence and the machine selection are optimized in the Hybrid Flow Shop, and the purposes of optimizing the delivery period of the whole order and the number of machines are achieved.
In the conventional mixed flow shop scheduling problem and its variant, the conventional algorithm is often designed around the same coding and decoding mode, and in the decoding stage, whether mixed rule decoding or positive and negative mixed decoding is adopted, a part of solution space cannot be searched. The algorithm makes the solution space complex by adding a complex domain searching mechanism and preferential performance, and the stability of algorithm solution is challenged.
The current algorithm and the model are separately solved, and constraint properties of the modulation model are not analyzed and solved.
Disclosure of Invention
Aiming at the above defects or improvement demands of the prior art, the invention provides a mathematical heuristic scheduling method and a mathematical heuristic scheduling system for a Mixed flow shop taking into consideration machine configuration.
To achieve the above object, according to one aspect of the present invention, there is provided a mathematical heuristic scheduling method of a hybrid flow shop considering machine configuration, the method comprising the steps of:
(1) Constructing a mixed integer linear programming model for a mixed flow shop considering machine configuration, wherein the optimization goal of the mixed integer linear programming model is to determine the number of machines in each stage and the processing sequence of work procedures on the machines, so that the weighted sum f=loga∑∑Mj,k+logb∑(Wi E·Ei+Wi T·Ti) of the delivery period deviation time and the number of machines is minimum, M j,k is whether the kth machine in the jth stage participates in processing, 1 represents that the machine participates in processing, or not; a, b are logarithmic parameters, and are obtained through parameter adjustment experiments; w i E is a weight value for advancing the finishing time; e i is the advance completion time; w i T is a weight value for delaying the completion time; t i is the delay completion time;
(2) Determining a coding structure of a mathematical heuristic algorithm based on decision variables of a mixed integer linear programming model, adding a coding sequence into a population as a chromosome, and decoding to obtain an optimal scheduling solution of a current round;
(3) Adding the optimal scheduling solution in the chromosome population of each iteration to the elite population, namely updating the solution in the elite population of each iteration to replace the last suboptimal solution so as to ensure that the solution in the elite population is a historical optimal solution;
(4) Determining a transmission rule of an elite population to a mixed integer linear programming model transmission variable based on an example scale, and determining a sub-model corresponding to the transmission rule, namely parameterizing the mixed integer linear programming model selected by the machine number of each stage and the machine of each stage of the workpiece;
(5) And solving the obtained submodel, wherein each solved value is used as the fitness function value of the current transfer chromosome and enters the population to perform superior and inferior elimination search, so as to obtain the final dispatching optimal solution.
Further, the coding structure is a double-layer coding structure, the first layer is a workpiece sequence, and the first layer represents a processing machine which can be used in each stage of workpiece procedure; the second layer is a workpiece sequence representing the priority relationship of a workpiece process at a certain stage.
Further, the double-layer coding sequence is added into the population as a chromosome, and the sequence of the coding structure is as follows:
[M1 M2 M3 ... Mj-1 Mj]
[S1 S2 ... Si-2 Si-1 Si]。
Further, during decoding, the chromosome coding sequence is used as a workpiece sequence at the first stage of the workpiece, and the earliest idle machine is orderly arranged for processing according to the machine number of each stage represented by the machine number coding sequence, so that a scheduling solution G Positive direction after forward decoding is obtained; then, the chromosome coding sequence is used as a workpiece sequence at the tail end stage of the workpiece, and the processing is sequentially carried out according to the latest idle machine according to the machine number of each stage represented by the machine number coding sequence, so as to obtain a scheduling solution G Reverse-rotation after reverse decoding; comparing the merits of the forward decoded scheduling solution G Positive direction and the backward decoded scheduling solution G Reverse-rotation , and selecting a better solution as the decoded scheduling solution G 0; and processing the decoded scheduling solution by using the idle time block movement rule with the delivery date deviation of each workpiece reduced as a target to obtain a better scheduling solution after processing, wherein the better scheduling solution is used as a final scheduling solution G Terminal (A) of the round.
Further, the forward decoding is to process the workpiece sequence provided by the chromosome of the workpiece sequence as a first stage workpiece sequence according to the rule of idle machine first processing; the reverse decoding is to process the workpiece sequence of the chromosome from right to left by taking the workpiece sequence of the chromosome as a processing sequence from right to left in the last stage, and also process the workpiece sequence by adopting a rule of idle machine processing first, and select a better solution in the two decoding processes as a decoding result.
Further, for the continuous time block on each machine at the end stage of the decoded feasible solution, judging whether the continuous time block is an independent time block or not sequentially from right to left, judging whether the independent continuous time block meets the condition of left shift or right shift, carrying out shifting by one unit each time, if the integral objective function value is reduced after shifting, shifting effectively, otherwise, not shifting.
Further, in the iterative process of elite population generation, firstly, performing crossover operation on the chromosome according to the randomly generated probability p, namely selecting two chromosome sequences, randomly selecting two crossover positions for a machine number layer and a workpiece process layer respectively, exchanging genes between crossover points, removing genes exchanged from another chromosome in an original chromosome, then sequentially filling residual genes to vacant positions, then, performing mutation operation on the chromosome according to the randomly generated probability p, namely randomly selecting one gene position in the sequences for the machine number layer and the workpiece process layer, randomly generating a value to replace the value of the original gene position by a standard machine number, and keeping the integrity of the sequence.
Further, the transmission rule comprises a rule A, and the machine number of each stage and the machine selection of each stage workpiece are transmitted to the MILP model, namely, the sub-model to be solved corresponding to the rule A is the MILP model parameterized by the machine number of each stage and the machine selection of each stage workpiece.
Further, decision variables include number of machines per stage, work piece sequence, and machine selection; solving the sub-model, taking the value obtained by each solving as the fitness function value of the transfer chromosome, entering the population to perform superior and inferior elimination search, continuing to perform rule judgment, sub-model determination and model solving on the chromosomes in the elite population in the generation until the elite population in the generation is completely optimized, judging whether the iteration time threshold of the algorithm is met, outputting the optimal solution if the iteration time threshold is met, and if the iteration time threshold is not met, continuing to generate part of elite population of the genetic algorithm, and circulating.
The invention provides a mathematical heuristic scheduling system of a mixed flow shop taking machine configuration into consideration, which comprises a memory and a processor, wherein the memory stores a computer program, and the processor executes the mathematical heuristic scheduling method of the mixed flow shop taking machine configuration into consideration when executing the computer program.
In general, compared with the prior art, the mixed flow shop mathematical heuristic scheduling method and system taking the machine configuration into consideration mainly have the following beneficial effects:
1. Compared with the existing hybrid algorithm, the hybrid meta-heuristic algorithm provided by the invention acquires elite population information by adopting the meta-heuristic algorithm, designs a transmission rule for a subsequence thereof, determines a specified corresponding submodel, and accurately solves the submodel, so that an optimal solution in the current field range can be obtained, the searching quality of the algorithm is improved, and a scheduling scheme for efficiently acquiring an optimal target is realized.
2. According to the invention, through model refinement of the feasible solution subsequence of the elite population, the elite population can acquire the current optimal solution, and on the basis, the model can search the optimal solution in the field range of the feasible solution.
3. The invention ensures the optimality of the decoding process by mixing the decoding method and the idle time block movement rule, and can obtain a better solution of the corresponding scheduling result in the coding space by adopting the mode because the objective function comprises the minimum advance/delay time weighted sum, thereby improving the searching efficiency.
4. For the mixed flow shop scheduling problem considering the machine configuration, when the scale of an example is increased, the searching efficiency of using other mixed algorithms is reduced, and the local optimization is possibly involved, and when the solving method is used, more feasible domain solutions can be searched, so that the global optimality is ensured.
Drawings
FIG. 1 is a flow chart of a mathematical heuristic scheduling method for a hybrid flow shop that takes into account machine configuration provided by the present invention;
FIG. 2 (a), (b), (c) and (d) are diagrams of chromosome encoding and decoding, respectively;
fig. 3 (a) and (b) are schematic diagrams of idle time block movement, respectively;
FIG. 4 is a schematic representation of elite population subsequence variable transfer.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention. In addition, the technical features of the embodiments of the present invention described below may be combined with each other as long as they do not collide with each other.
The invention provides a mathematical heuristic scheduling method of a mixed flow shop considering machine configuration, which mixes a mathematical solving method and a (meta) heuristic solving method by the mathematical heuristic method to realize graded solving of decision variables and is used for solving the technical problem that the algorithm is difficult to solve due to multiple decision variables and large solution space in the scheduling problem. The method comprises the steps of firstly modeling a mixed flow shop scheduling considering machine configuration, establishing a Mixed Integer Linear Programming (MILP) model, and carrying out fusion solution on the established MILP model by adopting a mathematical heuristic algorithm based on a genetic algorithm so as to solve the problems of regular decoding and solution space explosion in the traditional scheduling problem and improve the scheduling solution quality in the limit solution time.
Referring to fig. 1, the scheduling method mainly includes the following steps:
Step one, constructing a mixed integer linear programming model for a mixed flow shop considering machine configuration, wherein the optimization goal of the mixed integer linear programming model is to determine the number of machines in each stage and the processing sequence of work procedures on the machines, so that the weighted sum f=loga∑∑Mj,k+logbΣ(Wi E·Ei+Wi T·Ti) of the delivery date deviation time and the number of machines to be used is minimum, M j,k is whether the kth machine in the jth stage participates in processing, 1 represents that the machine participates in processing, and otherwise the machine does not participate in processing; a, b are logarithmic parameters, and are obtained through parameter adjustment experiments; w i E is a weight value for advancing the finishing time; e i is the advance completion time; w i T is a weight value for delaying the completion time; t i is the delay completion time.
Linearizing a traditional mixed flow shop scheduling model, realizing linearization of the model by introducing 0-1 variable and a large M method, taking machine configuration of each stage into consideration, and introducing machine number decision of each stage as model constraint to obtain an MILP model taking the mixed flow shop scheduling problem of the machine configuration into consideration, wherein the optimization goal of the model is to determine the number of machines of each stage and the processing sequence of workpiece procedures on the machines, so that the weighted sum of delivery period deviation time and the number of used machines is achieved:
f=loga∑∑Mj,k+logbΣ(Wi E·Ei+Wi T·Ti) Minimum.
The constraint formula of the workpiece priority relation on each stage of machine is as follows:
Wherein the variable x i,j,j' represents the priority relationship between the workpiece i and the workpiece i' at the j-th stage; i is the work set, i= {1, 2., i..n }, n is the number of work pieces; j is a set of phases, j= {1, 2., j..m, m being the number of stages.
The corresponding constraints on the number of machines used per stage of the work-piece process are:
Wherein, the variable M j,k is equal to 1, which indicates that the kth machine in the jth stage can be used, and the total number of the kth machine is smaller than the standard machine number; Standard machine number per stage; j is a set of phases, j= {1, 2., j..m, m being the number of stages.
The constraint formula corresponding to the machining machine number selected by the workpiece at each stage is as follows:
Wherein the variable y i,j,k represents whether the workpiece i is processed by the kth machine selected in the jth stage, and the total number of the workpiece i is smaller than the number of machines capable of being used; k is the machine label; k is the set of machines, k= {1,2, once again, K, M j; n is the number of work pieces; i is the workpiece set, i= {1,2, i., n }; j is the stage index, J is the stage set, j= {1,2,..j,..m }.
The constraint formula of the workpiece processing time of the same processing machine without conflict is as follows:
In each stage of selected processing machines, wherein the variable y i,j,k represents whether the workpiece i is processed in the j-th stage of selected kth machine, the total number of which is less than the number of machines that can be used; c i,j is the finishing time of the workpiece i in the j-th stage; m is the number of stages; c i,0 is the finishing time of the workpiece i in the first stage; t i',0 is the finishing time of the workpiece i' in the first stage; y i,0,k is whether the workpiece i is machined on machine k in the first stage, if machined, 1, otherwise 0; y i',0,k is whether the workpiece i' is machined on machine k in the first stage, 1 if machined, or 0 otherwise; x i,i',0 is the priority relationship between the first stage workpiece i and workpiece i'; a is a maximum positive number; t i',j is the processing time of the workpiece i' at stage j; y i',j,k is whether the workpiece i' is machined on machine k in the j-th stage, if so, 1, otherwise 0; x i,i',j is the priority relationship between the j-th stage workpiece i and workpiece i'.
The formula for determining the lead time variable and the lead delay time variable is:
in the formula, the variables And/>Respectively representing the lower limit and the upper limit of the delivery period of the workpiece i; e i is the advanced finishing time of the workpiece i; c i,m is the finishing processing of the workpiece i in the final stage; t i is the delayed finishing time of the workpiece i;
and secondly, encoding and decoding of a mathematical heuristic algorithm based on a mixed integer linear programming model.
And determining a coding structure of a mathematical heuristic algorithm based on decision variables of the mixed integer linear programming model, adding the coding sequence into the population as a chromosome, and decoding to obtain an optimal scheduling solution of the current round.
Encoding: designing a coding structure based on a genetic algorithm according to decision variables of a mixed integer linear programming model (a main model), wherein the main model comprises three decision variables of the number of machines in each stage, a workpiece sequence and machine selection, and the coding structure is designed into a double-layer coding structure, wherein a first layer is a workpiece sequence and represents a processing machine which can be used by a workpiece procedure in each stage; the second layer is a workpiece sequence and represents the priority relation of a workpiece process at a certain stage; adding the double-layer coding sequence as chromosome into population, and the coding structure sequence is as follows:
[M1 M2 M3 ... Mj-1 Mj]
[S1 S2 ... Si-2 Si-1 Si]
Decoding: firstly, a chromosome coding sequence is used as a workpiece sequence at the first stage of a workpiece, and the earliest idle machine is orderly arranged for processing according to the machine number of each stage represented by the machine number coding sequence, so as to obtain a scheduling solution G Positive direction after forward decoding; then, the chromosome coding sequence is used as a workpiece sequence at the tail end stage of the workpiece, and the processing is sequentially carried out according to the latest idle machine according to the machine number of each stage represented by the machine number coding sequence, so as to obtain a scheduling solution G Reverse-rotation after reverse decoding; comparing the merits of the forward decoded scheduling solution G Positive direction and the backward decoded scheduling solution G Reverse-rotation , and selecting a better solution as the decoded scheduling solution G 0; and processing the decoded scheduling solution by using the idle time block movement rule with the delivery date deviation of each workpiece reduced as a target to obtain a better scheduling solution after processing, wherein the better scheduling solution is used as a final scheduling solution G Terminal (A) of the round.
The decoding process of the 6×3 example shown in fig. 2, forward decoding, is to process according to the rule of idle machine first processing by taking the workpiece sequence provided by the workpiece sequence chromosome as the first stage workpiece sequence; the reverse decoding is to process the workpiece sequence of the chromosome from right to left by taking the workpiece sequence of the chromosome as a processing sequence from right to left in the last stage, and also process the workpiece sequence by adopting a rule of idle machine processing first, and select a better solution in the two decoding processes as a decoding result.
Judging whether the continuous time blocks on each machine at the end stage of the decoded feasible solution are independent time blocks from right to left in sequence, judging whether the independent continuous time blocks meet the condition of left shift or right shift, carrying out shifting by one unit each time, if the integral objective function value is reduced after shifting, shifting effectively, otherwise, not shifting;
Assuming, as in the 6 x 3 example shown in fig. 3, that the current workpiece time block profile on the last stage machine is as shown, then in M 3,1 the workpiece 3 does not meet the move condition, i.e. the start machining time cannot exceed the finish time of the stage on the workpiece and so cannot move; in M 3,2, the time block consisting of workpiece 6 and workpiece 1, after shifting one unit to the right, can be optimized to minimize the total weighted advance/retard time by one unit, so the shift is effective.
And thirdly, adding the optimal scheduling solution in the chromosome population of each iteration round into the elite population, namely updating the solution in the elite population of each iteration round to replace the last suboptimal solution so as to ensure that the solution in the elite population is the historical optimal solution.
Specifically, adding the optimal scheduling solution in the chromosome population of each iteration round into the elite population, and keeping the elite population size to be 10, namely updating the solution in the elite population of each iteration round, and exchanging the last suboptimal solution to ensure that the solution in the elite population is the historical optimal solution.
Iterative process of elite population generation: firstly, performing crossover operation on chromosomes according to randomly generated probability p, namely selecting two chromosome sequences, randomly selecting two crossover positions for a machine number layer and a workpiece process layer respectively, exchanging genes between crossover points, removing genes exchanged from another chromosome in an original chromosome, sequentially filling residual genes to vacant positions, then performing mutation operation on the chromosomes according to randomly generated probability p, namely randomly selecting one gene position in the sequence for the machine number layer and the workpiece process layer, randomly generating a value to replace the value of the original gene position by a standard machine number, and maintaining the integrity of the sequence.
And step four, determining a transmission rule of an elite population to a mixed integer linear programming model based on an example scale, and determining a sub-model corresponding to the transmission rule, namely, parameterizing the mixed integer linear programming model selected by the machine number of each stage and the machine of each stage of workpiece.
Determining the transmission rule of elite population to main model transmission variable according to the scale, wherein the larger the scale, the algorithm tends to select transmission of more variable, and the smaller the determined decision variable type of sub-model optimization. According to the number of machines at each stage in the problem characteristics and models, the priority relation among work procedures of the workpieces and the machine selection three decision variables of the workpieces at each stage, A, B, C three transmission rules and corresponding submodel representations to be solved are designed.
Wherein, rule A: and transmitting the number of machines per stage and the machine selection of the workpieces per stage to the MILP model, namely, parameterizing the MILP model of the machine selection of the machines per stage and the workpieces per stage according to the submodel to be solved corresponding to the rule A, wherein the constraint of the submodel is as follows:
constraining workpiece priority relationships on each machine per stage, where variables The priority relationship between the processes XY j,k,i and XY j,k,i' on the kth machine in the jth stage is shown.
The workpiece processing times on the same machine are constrained from conflicting, where the variable XY j,k,i represents the workpiece number that the elite population selected machine k to process at stage j.
In the method, in the process of the invention,The processing completion time of the workpiece XY j,k,i in the j stage; /(I)Processing time of the workpiece XY j,k,i' in the j stage; /(I)A priority relationship between the workpiece XY j,k,i and the workpiece XY j,k,i' on the kth machine of the jth stage; /(I)The processing completion time of the workpiece XY j,k,i' at stage j; i 1 is a workpiece set after machine distribution is completed; i' is the workpiece label.
Rule B: and transmitting the number of machines per stage, the machine selection of the workpieces per stage and the workpiece priority relation on each machine of the non-end stage to the MILP model, namely, parameterizing the sub-model to be solved corresponding to the rule B, namely, the MILP model with the number of machines per stage, the machine selection of the workpieces per stage and the workpiece priority relation on each machine of the non-end stage, wherein the constraint of the sub-model refers to the model after the parameterization of A.
Rule C: rule C, the number of machines per stage, the machine selection of the workpieces per stage and the workpiece processing sequence on the non-end stage and non-end stage machine are transferred to the MILP model, namely the sub-model to be solved corresponding to the rule C is the MILP model parameterized by the number of machines per stage, the machine selection of the workpieces per stage and the workpiece processing sequence on the non-end stage and non-end stage machine, and the constraint of the sub-model refers to the model parameterized by A.
After the decoding rule and the idle time block movement rule are determined, a better elite population can be ensured, but the solution represented by the elite population at the moment is not necessarily the optimal solution, and when the solution space is increased, the quality of the solution is more difficult to ensure, so that the implementation mode continuously optimizes the elite population, specifically:
Elite population, the transmitted elite population subsequence rule is determined according to the size of the solution space corresponding to the scale of the example and the dimension of the variable, and when the solution space is too large, the model receives more transmitted variables.
As shown in FIG. 4, in a 10X 3 embodiment, among the feasible solutions in elite populations of the current genetic algorithm, a 3X 1 size machine quantity sequence [ 23 2], a workpiece sequence [ 61 94 810 725 3] and a machine selection sequence [ 122 122 121 1] are included by rational evaluation of the selection rules.
After the rule of the transmitted subsequence is obtained, on the basis of the main model, the parameterization of the variables is set, the obtained decision variables are fixed as model parameters, the variable range can be reduced, the variable dimension is reduced, and the MILP submodel with smaller corresponding solution space is obtained through a relaxation model.
And fifthly, solving the obtained submodel, wherein each solved value is used as an fitness function value of the current transfer chromosome and enters a population to perform superior and inferior elimination search, so that a final dispatching optimal solution is obtained.
In this embodiment, the fifth step includes the following sub-steps:
1. For the machine number variables and machine selection variables of the embodiment of fig. 4, they are used to parameterize the main model, convert the sequence variables into binary 0-1 variables in the model, determine the variable y i,j,k, and convert the parameterized decision variables represented by the MILP submodel into the priority relationships of the workpieces on each machine, the variable dimension is reduced, and the specific constraints of the MILP submodel are as follows:
(1) The priority relation constraint of the workpieces on each machine at each stage is used for constraining the workpiece priority on each machine:
(2) The work pieces on each machine occupy time constraint, and only one work piece can be processed on each machine at the same time:
2. Model solving, taking a function f=log bΣ(Wi E·Ei+Wi T·Ti) as an objective function, solving the submodel by calling a linear solving module in the IBM CPLEX, taking the value solved each time as an fitness function value of the transfer chromosome, and entering a population to perform quality elimination search.
3. And continuing to perform rule judgment, sub-model determination and model solving on the subsequences in the elite population in the generation until the elite population in the generation is completely optimized, judging whether the iteration time threshold of the algorithm is met, outputting an optimal solution if the iteration time threshold is met, and if the iteration time threshold is not met, continuing to generate part of elite population of the genetic algorithm, and circulating.
Based on the model and the solving algorithm, in order to verify the practical application effect of the invention, the invention carries out repeated tests on a self-generated random computing example and compares the repeated tests with other common algorithms for solving the mixed flow shop scheduling, wherein GMHA is the result of the invention.
As the MILP result shown in table 1 is the result of directly solving the problem by the main model, it can be seen that the mixed integer linear programming submodel established by the invention has good solving effect, and the solving example has larger scale and smaller time; the GMHA algorithm results shown in table 1 are mathematical heuristic algorithm results of the present invention, and the solution quality is significantly superior compared to pure GA and hybrid algorithms employing full encoding; table 2 shows the results of solving the problem by several common scheduling algorithms, and the average solution of the optimal solution for multiple repeated experiments is not as good as the mathematical heuristic algorithm proposed by the invention. Therefore, the method provided by the invention has obvious advantages in solving the scheduling problem of the mixed flow shop considering the machine configuration, has high calculation result and has guiding effect in actual production.
Table 1 smaller scale example test results
TABLE 2 larger Scale example test results
As shown in tables 1 and 2, the numerical values in the tables represent the optimal solution for solving the corresponding algorithm and the average solution for ten repeated experiments, and it can be seen that the method provided by the invention can solve the optimal solution better than the MILP, is better than other meta-heuristic algorithms and mixing algorithms, achieves advantages in the scheduling problem of the mixed flow shop considering the machine configuration, and can be used for guiding actual production.
The invention also provides a mathematical heuristic scheduling system of the mixed flow shop considering the machine configuration, which comprises a memory and a processor, wherein the memory stores a computer program, and the processor executes the mathematical heuristic scheduling method of the mixed flow shop considering the machine configuration when executing the computer program.
It will be readily appreciated by those skilled in the art that the foregoing description is merely a preferred embodiment of the invention and is not intended to limit the invention, but any modifications, equivalents, improvements or alternatives falling within the spirit and principles of the invention are intended to be included within the scope of the invention.

Claims (10)

1. A mathematical heuristic scheduling method for a hybrid flow shop taking machine configuration into account, the method comprising the steps of:
(1) Constructing a mixed integer linear programming model for a mixed flow shop considering machine configuration, wherein the optimization goal of the mixed integer linear programming model is to determine the number of machines in each stage and the processing sequence of work procedures on the machines, so that the weighted sum f=loga∑∑Mj,k+logb∑(Wi E·Ei+Wi T·Ti) of the delivery period deviation time and the number of machines is minimum, M j,k is whether a kth machine in a jth stage participates in processing, 1 represents participation in processing, or not; a, b are logarithmic parameters, and are obtained through parameter adjustment experiments; w i E is a weight value for advancing the finishing time; e i is the advance completion time; w i T is a weight value for delaying the completion time; t i is the delay completion time;
(2) Determining a coding structure of a mathematical heuristic algorithm based on decision variables of a mixed integer linear programming model, adding a coding sequence into a population as a chromosome, and decoding to obtain an optimal scheduling solution of a current round;
(3) Adding the optimal scheduling solution in the chromosome population of each iteration to the elite population, namely updating the solution in the elite population of each iteration to replace the last suboptimal solution so as to ensure that the solution in the elite population is a historical optimal solution;
(4) Determining a transmission rule of an elite population to a mixed integer linear programming model transmission variable based on an example scale, and determining a sub-model corresponding to the transmission rule, namely parameterizing the mixed integer linear programming model selected by the machine number of each stage and the machine of each stage of the workpiece;
(5) Solving the obtained submodel, wherein each solved value is used as an fitness function value of the current transfer chromosome and enters a population to perform superior and inferior elimination search, so as to obtain a final dispatching optimal solution;
The constraint formula of the workpiece priority relation on each stage of machine is as follows:
Wherein the variable x i,j,j' represents the priority relationship between the workpiece i and the workpiece i' at the j-th stage; i is the work set, i= {1, 2., i..n }, n is the number of work pieces; j is the set of stages, j= {1,2,..j,..once, m }, m is the number of stages;
the corresponding constraints on the number of machines used per stage of the work-piece process are:
Wherein, the variable M j,k is equal to 1, which indicates that the kth machine in the jth stage can be used, and the total number of the kth machine is smaller than the standard machine number; Standard machine number per stage;
the constraint formula corresponding to the machining machine number selected by the workpiece at each stage is as follows:
Wherein the variable y i,j,k represents whether the workpiece i is processed by the kth machine selected in the jth stage, and the total number of the workpiece i is smaller than the number of machines capable of being used; k is the machine label; k is the set of machines, k= {1,2, once again, K, M j; i is the workpiece set, i= {1,2, i., n }; j is a stage index;
the determination formulas of the advance completion time E i and the retard completion time T i are:
in the formula, the variables And/>Respectively representing the lower limit and the upper limit of the delivery period of the workpiece i; c i,m is the finishing of the workpiece i in the final stage.
2. The hybrid flow shop mathematical heuristic scheduling method considering machine configuration as claimed in claim 1, wherein: the coding structure is a double-layer coding structure, the first layer is a workpiece sequence, and represents a processing machine which can be used in each stage of workpiece procedure; the second layer is a workpiece sequence representing the priority relationship of a workpiece process at a certain stage.
3. The hybrid flow shop mathematical heuristic scheduling method taking into account machine configuration as set forth in claim 2, wherein: adding the double-layer coding sequence as chromosome into population, and the coding structure sequence is as follows:
[M1 M2 M3 ... Mj-1 Mj]
[S1 S2 ... Si-2 Si-1 Si]。
4. The hybrid flow shop mathematical heuristic scheduling method considering machine configuration as claimed in claim 1, wherein: when decoding, the chromosome coding sequence is used as a workpiece sequence at the first stage of the workpiece, and the earliest idle machine is orderly arranged for processing according to the machine number of each stage represented by the machine number coding sequence, so as to obtain a scheduling solution G Positive direction after forward decoding; then, the chromosome coding sequence is used as a workpiece sequence at the tail end stage of the workpiece, and the processing is sequentially carried out according to the latest idle machine according to the machine number of each stage represented by the machine number coding sequence, so as to obtain a scheduling solution G Reverse-rotation after reverse decoding; comparing the merits of the forward decoded scheduling solution G Positive direction and the backward decoded scheduling solution G Reverse-rotation , and selecting a better solution as the decoded scheduling solution G 0; and processing the decoded scheduling solution by using the idle time block movement rule with the delivery date deviation of each workpiece reduced as a target to obtain a better scheduling solution after processing, wherein the better scheduling solution is used as a final scheduling solution G Terminal (A) of the round.
5. The hybrid flow shop mathematical heuristic scheduling method considering machine configuration as claimed in claim 1, wherein: the forward decoding is to process the workpiece sequence provided by the chromosome of the workpiece sequence as a first stage workpiece sequence according to the rule of idle machine first processing; the reverse decoding is to process the workpiece sequence of the chromosome from right to left by taking the workpiece sequence of the chromosome as a processing sequence from right to left in the last stage, and also process the workpiece sequence by adopting a rule of idle machine processing first, and select a better solution in the two decoding processes as a decoding result.
6. A machine configuration-aware hybrid flow shop mathematical heuristic scheduling method according to any of claims 1-5, wherein: and judging whether the continuous time blocks on each machine at the end stage of the decoded feasible solution are independent time blocks from right to left, judging whether the independent continuous time blocks meet the condition of left shift or right shift, carrying out shifting by one unit each time, and if the integral objective function value is reduced after shifting, enabling shifting, otherwise, not shifting.
7. A machine configuration-aware hybrid flow shop mathematical heuristic scheduling method according to any of claims 1-5, wherein: in the iteration process of elite population generation, firstly, carrying out crossover operation on chromosomes according to randomly generated probability p, namely selecting two chromosome sequences, randomly selecting two crossover positions for a machine number layer and a workpiece process layer respectively, exchanging genes between crossover points, removing genes exchanged from another chromosome in an original chromosome, then sequentially filling residual genes to vacant positions, then carrying out mutation operation on the chromosomes according to randomly generated probability p, namely randomly selecting one gene position in the sequences for the machine number layer and the workpiece process layer, randomly generating a value for replacing the value of the original gene position by a standard machine number, and maintaining the integrity of the sequences.
8. A machine configuration-aware hybrid flow shop mathematical heuristic scheduling method according to any of claims 1-5, wherein: the transmission rule comprises a rule A, and the number of machines in each stage and the machine selection of the workpieces in each stage are transmitted to the MILP model, namely, the sub-model to be solved corresponding to the rule A is the MILP model parameterized by the number of machines in each stage and the machine selection of the workpieces in each stage.
9. A machine configuration-aware hybrid flow shop mathematical heuristic scheduling method according to any of claims 1-5, wherein: decision variables include number of machines per stage, work piece sequence, and machine selection; solving the sub-model, taking the value obtained by each solving as the fitness function value of the transfer chromosome, entering the population to perform superior and inferior elimination search, continuing to perform rule judgment, sub-model determination and model solving on the chromosomes in the elite population in the generation until the elite population in the generation is completely optimized, judging whether the iteration time threshold of the algorithm is met, outputting the optimal solution if the iteration time threshold is met, and if the iteration time threshold is not met, continuing to generate part of elite population of the genetic algorithm, and circulating.
10. A mathematical heuristic scheduling system for a hybrid flow shop taking machine configuration into account, characterized by: the system comprising a memory storing a computer program and a processor executing the computer program to perform the hybrid flow shop mathematical heuristic scheduling method taking into account machine configuration of any of claims 1-9.
CN202310337907.6A 2023-03-31 2023-03-31 Mixed flow shop mathematical heuristic scheduling method and system considering machine configuration Active CN116300763B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310337907.6A CN116300763B (en) 2023-03-31 2023-03-31 Mixed flow shop mathematical heuristic scheduling method and system considering machine configuration

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310337907.6A CN116300763B (en) 2023-03-31 2023-03-31 Mixed flow shop mathematical heuristic scheduling method and system considering machine configuration

Publications (2)

Publication Number Publication Date
CN116300763A CN116300763A (en) 2023-06-23
CN116300763B true CN116300763B (en) 2024-05-24

Family

ID=86795890

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310337907.6A Active CN116300763B (en) 2023-03-31 2023-03-31 Mixed flow shop mathematical heuristic scheduling method and system considering machine configuration

Country Status (1)

Country Link
CN (1) CN116300763B (en)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018018270A (en) * 2016-07-27 2018-02-01 株式会社東芝 Production schedule preparation system and production schedule preparation method
CN107678411A (en) * 2017-10-16 2018-02-09 华中科技大学 A kind of modeling method of uncorrelated parallel machine hybrid flow shop scheduling
CN110221585A (en) * 2019-06-14 2019-09-10 同济大学 A kind of energy-saving distribution control method considering plant maintenance for hybrid flowshop
CN110782085A (en) * 2019-10-23 2020-02-11 武汉晨曦芸峰科技有限公司 Casting production scheduling method and system
CN111353604A (en) * 2018-12-24 2020-06-30 南京理工大学 Flexible assembly multi-objective dynamic optimization method
CN111966049A (en) * 2020-06-22 2020-11-20 同济大学 Scheduling control method for production equipment of mixed flow shop
CN111967654A (en) * 2020-07-27 2020-11-20 西安工程大学 Method for solving flexible job shop scheduling based on hybrid genetic algorithm
CN112561194A (en) * 2020-12-22 2021-03-26 华中科技大学 Production and logistics integrated scheduling method and system for hybrid flow shop
US10970682B1 (en) * 2015-06-04 2021-04-06 Incontact, Inc. System and method for agent scheduling using mixed integer linear programming
CN113361860A (en) * 2021-05-12 2021-09-07 同济大学 Flexible flow shop scheduling control method, medium and equipment considering fatigue degree
CN114493337A (en) * 2022-02-15 2022-05-13 河南中烟工业有限责任公司 Flexible job shop scheduling method based on improved particle swarm genetic hybrid algorithm

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108320057B (en) * 2018-02-07 2021-06-18 江南大学 Flexible job shop scheduling method based on limited stable pairing strategy

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10970682B1 (en) * 2015-06-04 2021-04-06 Incontact, Inc. System and method for agent scheduling using mixed integer linear programming
JP2018018270A (en) * 2016-07-27 2018-02-01 株式会社東芝 Production schedule preparation system and production schedule preparation method
CN107678411A (en) * 2017-10-16 2018-02-09 华中科技大学 A kind of modeling method of uncorrelated parallel machine hybrid flow shop scheduling
CN111353604A (en) * 2018-12-24 2020-06-30 南京理工大学 Flexible assembly multi-objective dynamic optimization method
CN110221585A (en) * 2019-06-14 2019-09-10 同济大学 A kind of energy-saving distribution control method considering plant maintenance for hybrid flowshop
CN110782085A (en) * 2019-10-23 2020-02-11 武汉晨曦芸峰科技有限公司 Casting production scheduling method and system
CN111966049A (en) * 2020-06-22 2020-11-20 同济大学 Scheduling control method for production equipment of mixed flow shop
CN111967654A (en) * 2020-07-27 2020-11-20 西安工程大学 Method for solving flexible job shop scheduling based on hybrid genetic algorithm
CN112561194A (en) * 2020-12-22 2021-03-26 华中科技大学 Production and logistics integrated scheduling method and system for hybrid flow shop
CN113361860A (en) * 2021-05-12 2021-09-07 同济大学 Flexible flow shop scheduling control method, medium and equipment considering fatigue degree
CN114493337A (en) * 2022-02-15 2022-05-13 河南中烟工业有限责任公司 Flexible job shop scheduling method based on improved particle swarm genetic hybrid algorithm

Also Published As

Publication number Publication date
CN116300763A (en) 2023-06-23

Similar Documents

Publication Publication Date Title
CN110796355B (en) Flexible job shop scheduling method based on dynamic decoding mechanism
CN112561194B (en) Integrated scheduling method and system for production and logistics of mixed flow shop
CN110598941A (en) Bionic strategy-based dual-target scheduling method for particle swarm optimization manufacturing system
CN101901425A (en) Flexible job shop scheduling method based on multi-species coevolution
CN110288185B (en) Distributed flexible pipeline scheduling method
CN112053037B (en) Flexible PCB workshop scheduling optimization method and system
CN113705978B (en) Static and dynamic integrated decision-making method and system for multi-machine task cutter
CN111291969B (en) Automobile reordering method based on genetic algorithm
CN110378583B (en) Method for interchanging adjacent procedures of pseudo-critical path and equipment
Cheng et al. New benchmark algorithms for No-wait flowshop group scheduling problem with sequence-dependent setup times
CN110929930A (en) Scheduling and scheduling optimization method for marine crankshaft production line
CN114595633A (en) Multi-constraint-based multi-target flexible job shop energy-saving scheduling method
CN117132181B (en) Distributed flexible production and transportation cooperative scheduling method
CN116300763B (en) Mixed flow shop mathematical heuristic scheduling method and system considering machine configuration
CN113762811A (en) Method and system for solving non-stalled Job Shop scheduling problem considering overtime
CN117555305A (en) NSGAII-based multi-target variable sub-batch flexible workshop job scheduling method
CN110705844A (en) Robust optimization method of job shop scheduling scheme based on non-forced idle time
CN116468137A (en) Distributed process planning and workshop scheduling integrated optimization method
CN113177667A (en) Electrical equipment manufacturing resource optimal configuration method based on improved genetic algorithm
Xiuli et al. Greedy simulated annealing algorithm for solving hybrid flow shop scheduling problem with re-entrant batch processing machine
CN113011767A (en) Hybrid genetic method for scheduling multi-target flexible job shop
CN115700647A (en) Workshop flexible operation scheduling method based on tabu search genetic algorithm
CN114330940A (en) Multi-target mixed spider monkey optimization method for PCB electroplating task sequencing problem
CN114676987B (en) Intelligent flexible job shop active scheduling method based on hyper-heuristic algorithm
CN117495041A (en) Multi-target flexible job shop scheduling method based on mixed wolf algorithm

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
GR01 Patent grant
GR01 Patent grant